STCAT_Herrera_Visualization

1. load libraries

2. Load Data into Seurat

ss_herrera <- readRDS("STCAT/All_Patients_Integrated_ss_Herrera_annotated.RDS")

3. Visualizations

# Manually enter your data as a named vector
prediction_counts <- c(
  "None T" = 9622,
  "CD4 Tcm" = 4935,
  "CD4 Tn" = 3927,
  "Tgd" = 1531,
  "CD8 Tn" = 1381,
  "CD8 Temra" = 1264,
  "CD4 Trm cell-death" = 1134,
  "CD4 Tex" = 906,
  "CD8 Tc" = 743,
  "CD4 Tisg" = 702,
  "CD8 Tem" = 669,
  "CD4 Trm" = 615,
  "CD8 Trm naive-like" = 476,
  "MAIT" = 385,
  "CD4 Tc" = 361,
  "CD4 Tem" = 272,
  "CD4 Treg" = 228,
  "CD4 Tn adhesion" = 163,
  "CD8 Trm" = 142,
  "CD4 Treg naive-like" = 142,
  "CD4 Th17" = 137,
  "CD8 proliferation" = 134,
  "CD4 proliferation" = 108,
  "CD4 Temra" = 62,
  "CD4 activated" = 40,
  "Double Negative" = 32,
  "CD4 Tfh" = 31,
  "CD8 Tstr" = 26,
  "CD8 Tcm" = 19,
  "other CD8 T" = 9,
  "CD8 senescence" = 4,
  "CD8 Tex" = 2
)

# Convert to data frame
df <- data.frame(
  Prediction = factor(names(prediction_counts), levels = names(sort(prediction_counts))),
  Count = prediction_counts
)

# Load ggplot2
library(ggplot2)

# Horizontal barplot (flipped coordinates)
ggplot(df, aes(x = Count, y = Prediction)) +
  geom_bar(stat = "identity", fill = "steelblue") +
  theme_minimal() +
  labs(
    title = "Cell Type Predictions from STCAT",
    x = "Number of Cells",
    y = "Predicted Cell Type"
  ) +
  theme(axis.text.y = element_text(size = 8))

NA
NA

Dimplot

DimPlot(ss_herrera, group.by = "tissue", label = TRUE, repel = T,label.box = T)

DimPlot(ss_herrera, group.by = "sample_id", label = TRUE, repel = T,label.box = T)

DimPlot(ss_herrera, group.by = "RNA_snn_res.0.5", label = TRUE, repel = T, label.box = T)

DimPlot(ss_herrera, group.by = "STCAT_Prediction", label = F, repel = F, label.box = T)

Distributions

table(ss_herrera$STCAT_Prediction)  # should show all predicted labels

      CD4 activated   CD4 proliferation              CD4 Tc             CD4 Tcm             CD4 Tem           CD4 Temra 
                 40                 108                 361                4935                 272                  62 
            CD4 Tex             CD4 Tfh            CD4 Th17            CD4 Tisg              CD4 Tn     CD4 Tn adhesion 
                906                  31                 137                 702                3927                 163 
           CD4 Treg CD4 Treg naive-like             CD4 Trm  CD4 Trm cell-death   CD8 proliferation      CD8 senescence 
                228                 142                 615                1134                 134                   4 
             CD8 Tc             CD8 Tcm             CD8 Tem           CD8 Temra             CD8 Tex              CD8 Tn 
                743                  19                 669                1264                   2                1381 
            CD8 Trm  CD8 Trm naive-like            CD8 Tstr     Double Negative                MAIT              None T 
                142                 476                  26                  32                 385                9622 
        other CD8 T                 Tgd 
                  9                1531 

table(ss_herrera$sample_id)         # should show L1–L7

Healthy_Blood  Healthy_Skin     MF1_Blood      MF1_Skin     SS1_Blood      SS1_Skin     SS2_Blood      SS2_Skin 
         4386            33          4880           293          3631           374          6078           244 
    SS3_Blood      SS3_Skin     SS4_Blood      SS4_Skin     SS5_Blood     SS6_Blood 
         2109           272          2072           320          2832          2821 

table(ss_herrera$tissue) 

Blood  Skin 
28809  1536 

# Contingency table 
table_annotation_tissue <- table(ss_herrera$STCAT_Prediction, ss_herrera$tissue)
table_annotation_sample_id <- table(ss_herrera$STCAT_Prediction, ss_herrera$sample_id)
table_annotation_cluster <- table(ss_herrera$STCAT_Prediction, ss_herrera$RNA_snn_res.0.5)

library(pheatmap)

# Proportion table (color scale)
prop_table <- prop.table(table_annotation_cluster, margin = 1)

# Heatmap with raw counts as labels
pheatmap(prop_table,
         display_numbers = table_annotation_cluster,
         cluster_rows = TRUE,
         cluster_cols = TRUE,
         color = colorRampPalette(c("lightyellow", "yellow", "lightgreen", "white"))(100),
         fontsize_row = 10,
         fontsize_col = 10,
         main = "STCAT Prediction vs Clusters")

library(pheatmap)

# Proportion table (color scale)
prop_table <- prop.table(table_annotation_sample_id, margin = 1)

# Heatmap with raw counts as labels
pheatmap(prop_table,
         display_numbers = table_annotation_sample_id,
         cluster_rows = TRUE,
         cluster_cols = TRUE,
         color = colorRampPalette(c("lightyellow", "yellow", "lightgreen", "white"))(100),
         fontsize_row = 10,
         fontsize_col = 10,
         main = "STCAT Prediction vs Patients")

# Proportion table (color scale)
prop_table <- prop.table(table_annotation_tissue, margin = 1)

# Heatmap with raw counts as labels
pheatmap(prop_table,
         display_numbers = table_annotation_tissue,
         cluster_rows = TRUE,
         cluster_cols = TRUE,
         color = colorRampPalette(c("lightyellow", "yellow", "lightgreen", "white"))(100),
         fontsize_row = 10,
         fontsize_col = 10,
         main = "STCAT Prediction vs Patients")

NA
NA

library(ggplot2)
library(ggridges)

ggplot(ss_herrera@meta.data, aes(x = STCAT_Uncertainty, y = STCAT_Prediction, fill = STCAT_Prediction)) +
  geom_density_ridges(alpha = 0.7, scale = 2, rel_min_height = 0.01) +
  theme_minimal() +
  labs(
    x = "Cluster Entropy Score",
    y = "Predicted T Cell Type (STCAT)",
    title = "Density Ridgeline of Cluster Entropy by T Cell Type"
  ) +
  theme(
    legend.position = "none",
    axis.text.y = element_text(size = 9)
  )

NA
NA