scRNA-seq Viewer
Info Studies were incorporated in a barcode, matrix, feature format. This format can be found on the 10x Genomics website for processing of single cell studies to obtain the gene expression matrices. For each study, the metadata incorporated in GEO were manually curated into profiles and the samples were separated based on applicable groups and conditions. The expression data from the cells for the samples within the same profile and condition were aggregated into an expression matrix with the cell barcodes having the sample name appended to it to ensure unique cell names.

Description (from GEO)

Submission Date: Apr 11, 2022

Summary: Type 1 diabetes (T1D) is an autoimmune disorder defined by CD8 T cell-mediated destruction of pancreatic β cells. Our previous work has shown that diabetogenic CD8 T cells in the islets of the non-obese diabetic (NOD) mouse model of T1D are phenotypically heterogenous, but CD8 T cell clonal heterogeneity in this model remains relatively unexplored. Here, we use paired single-cell RNA sequencing (scRNA-seq) and single-cell T cell receptor sequencing (scTCR-seq) to characterize autoreactive CD8 T cells from the islets and spleens of NOD mice. Clonal analysis of scTCR-seq data demonstrates that CD8 T cells targeting the immunodominant β cell epitope IGRP206-214 exhibit highly restricted TCR gene usage, with over 70% of IGRP206-214-reactive cells utilizing the same TCR V alpha, J alpha, and V beta genes. Despite this, we observe only 5% overlap of IGRP206-214-reactive CD8 T cell clones between two groups of 10 NOD mice, demonstrating the immense TCR heterogeneity generated by junctional diversity during V(D)J recombination. scRNA-seq identifies two new clusters of autoreactive CD8 T cells in the islets and six clusters of diabetogenic CD8 T cells in the spleen, including multiple memory-like clusters and a population of exhausted cells. Strong clonal overlap between IGRP206-214-reactive CD8 T cells in the islets and spleen suggests that these cells may exit the islets and enter the peripheral circulation. Finally, we identify correlations between TCR J beta gene usage, which is less restricted than that of other TCR genes, and CD8 T cell clonal expansion as well as effector fate. Collectively, our work demonstrates that IGRP206-214-specific CD8 T cells are phenotypically heterogeneous but clonally similar, raising the possibility of selectively targeting either conserved or divergent TCR structures of these cells for therapeutic benefit.

GEO Accession ID: GSE200608

PMID: 35667687

Select Condition to compute Visualization and Analysis
Info Preprocessing and downstream analysis were computed using the scanpy Python library and the steps of processing followed the Seurat vignette. Cells and genes with no expression or very low expression were removed from the dataset based on a predefined threshold. The data was then normalized across the expression within the cells and log normalized. The top 2000 highly variable genes were extracted to be used for downstream analysis. For each of these aggregated data matrices, the clusters were computed using the leiden algorithm. Scanpy was then used to compute the PCA, t-SNE, and UMAPs. The points in the plots are labelled by their corresponding cell type labels. The cell type labels were computed using the wilcoxon method as the differential gene expression method. The top 250 genes were then used for enrichment analysis against the CellMarker library in order to determine the most appropriate cell type label with the lowest p-value.
PCA, tSNE, and UMAP Visualizations Labeled by Cell Types:
Gene Expression Data Explorer
Enter gene symbol:

Select conditions below to toggle them from the plot:

Cell Types Cell Samples
Differential Gene Expression Analysis
Info Differential gene expression can be computed for a single cell type labeled group of cells vs the rest. These include wilcoxon, DESeq2, or characteristic direction.
Select Cluster for differential gene expression analysis:
Select differential expression analysis method: