Dataset information
This report has been verified by Polly as per framework v1.0 Learn More
| Dataset information | Value |
|---|---|
| Dataset ID | E-MTAB-10948_raw_custom_processed |
| Abstract | Regulatory T cells (Tregs) play an important role in controlling inflammation and limiting autoimmunity, but their phenotypes at inflammatory sites in human disease are poorly understood. We here analyze the single-cell transcriptome of >16,000 Tregs obtained from peripheral blood and synovial fluid of two patients with HLA-B27+ ankylosing spondylitis and three patients with psoriatic arthritis, closely related forms of inflammatory spondyloarthritis. We identify multiple Treg clusters with distinct transcriptomic profiles, including, among others, a regulatory CD8+ subset expressing cytotoxic markers/genes, and a Th17-like RORC+ Treg subset characterized by IL-10 and LAG-3 expression. Synovial Tregs show upregulation of interferon signature and TNF receptor superfamily genes, and marked clonal expansion, consistent with tissue adaptation and antigen contact respectively. Individual synovial Treg clones map to different clusters indicating cell fate divergence. Finally, we demonstrate that LAG-3 directly inhibits IL-12/23 and TNF secretion by patient-derived monocytes, a mechanism with translational potential in SpA. Our detailed characterization of Tregs at an important inflammatory site illustrates the marked specialization of Treg subpopulations. |
| Description | Single cell analysis of spondyloarthritis regulatory T cells identifies distinct synovial gene expression patterns and clonal fates |
| Number of cells | 13954 |
| Number of genes | 15874 |
| Number of samples | 4 |
| Organism | Homo Sapiens |
| Tissue | Blood, Synovia |
| Disease | Spondylitis, Ankylosing |
| Cell Lines | None |
| Cell Type | |
| Drug | None |
| Marker genes for cell type are available | True |
| Doublet detection method | scrublet |
| Normalization method | log1p: true; target_sum: none; scaling_applied: true; max_value: none; zero_center: false |
| Remove gene groups | none |
| Batch correction method and key | batch_removal_method: harmony; batch_key: sample_id |
| Regress covariates | none |
1. Distribution of Key Quality Control Metrics
Figure 1: These violin plots display the distribution of quality control metrics for each cell. Metrics include the number of genes detected, total transcript counts, and the percentage of mitochondrial transcripts.
| A good-quality dataset would typically have a reasonable number of genes detected per cell and a moderate total transcript count. High mitochondrial transcript percentages can indicate low-quality, dying cells. Please Note: certain datasets do not have mitochondrial genes (MT-), thus figure for percentage of mitochondrial transcripts may be empty. |
2. UMAP visualization of cells colored by sample
Figure 2: Sample level distribution of clustering pattern of cells with the help of UMAP embeddings.
| If cells from the same sample cluster together distinctly from cells of other samples, it may indicate the presence of batch effects. Ideally, cells should be mixed and group based on their biological characteristics rather than their originating sample, indicating that the data is free of significant batch effects and the samples are comparable. |
3. Stacked barplot of cell types distributed across samples
Figure 3: The bar plot showcases the distribution and abundance of different cell types within each sample. Each color in a bar represents a different cell type with the height of the color segment indicating the count of that cell type in the sample.
| A uniform distribution of cell types across samples, may suggest that the sample preparation and preprocessing methods used were effective and there was minimal bias or variation in the processing steps. In some cases, if the experiment design ensures enrichment of a cell-type in a sample, then a non-uniform distribution is also valid. |
4. Stacked barplot of clusters distributed across samples
Figure 4: The bar plot showcases the distribution and abundance of different clusters within each sample. Each color in a bar represents a different cluster with the height of the color segment indicating the count of that cluster in the sample.
| Generally, a uniform distribution of clusters across samples, suggests there was minimal bias or variation in the processing steps. |
5. Stacked barplot of cell-types distributed across clusters
Figure 5: The bar plot showcases the distribution and abundance of different cell types within each cluster. Each color in a bar represents a different cell-type with the height of the color segment indicating the count of that cell-type in the cluster.
| Generally, each cluster should have only one cell-type to indicate accurate cell-type annotation. A corner-cases are observed when the authors have only provided cell ID to cell-type mapping and no marker genes. These need to manually rectified. |
6. Distribution of (a) Cell Counts (b) Median Gene Counts (c) Median Mitochondrial Genes, across Samples
Figure 6a: The bar plot visualizes the total count of cells detected in each sample. Each bar corresponds to a different sample, with its height representing the number of cells.
| This plot provides an understanding of the sample distribution in terms of cellularity. A wide variance in cell numbers across samples might indicate inconsistencies in cell isolation, sample preparation, or sequencing depth. Consistent cell counts across samples, however, would suggest a more uniform sampling process. |
Figure 6b: The bar plot illustrates the median number of genes detected in each sample. Each bar represents a different sample, and its height corresponds to the median gene counts.
| Consistently low gene counts might indicate low sequencing depth or poor-quality samples. On the other hand, large variances between samples or cell types might point to technical biases or true biological differences. |
Figure 6c: The bar plot showcases the median percentage of mitochondrial gene transcripts across samples.
| Consistently high mitochondrial gene percentages across samples might indicate a widespread issue with cell viability, while sporadic high values could suggest sample-specific issues which can be removed before downstream analysis |
7. Gene Counts Distribution
Figure 7: The plot provides a smoothed representation of the distribution of detected genes across cells.
| This plot gives an idea about the average gene richness in cells. High variability might indicate a mix of high and low-quality cells. |
8. UMI Count Distribution
Figure 8: The plot provides a smoothed representation of the distribution of UMIs across cells.
| This plot offers insight into the typical transcriptomic depth of the dataset. A broad distribution might indicate variability in sequencing depth across cells. |
9. UMI vs Gene counts distribution scatter plots colored by density
Figure 9: The scatter plot provides a visual representation of the relationship between the number of unique molecular identifiers (UMIs) and the number of genes detected in single cells. The color intensity indicates the density of data points in a particular region of the plot, allowing for the identification of trends and patterns.
| Ideally, one would expect to see a positive correlation between UMIs and genes, indicating that cells with more transcripts also express more unique genes. Areas with higher density may represent the most typical cells in the dataset, while outliers could indicate low-quality cells or potential doublets. |
10. Batch Mixing Metrics
| NMI | ARI | PCR_batch | Graph_iLISI | kBET_accept_rate | batch_correction_score | |
|---|---|---|---|---|---|---|
| uncorrected | 0.5026850334362329 | 0.7189496302097514 | 0.8589766833708714 | 0.002700964609781901 | 0.002364913286512828 | 0.0 |
| corrected | 0.7562660955615742 | 0.8831431021745508 | 0.999641646701263 | 0.05307745933532715 | 0.04944818689981367 | 1.0 |
Table 1: Table displaying batch mixing metrics
| These metrics are adopted from a recent benchmarking study of single-cell integration methods (Lueken et al. 2022). Values closer to 1 indicate better mixing of cells from the different batches. |
11. Cell Type Annotation Metrics
| sc_cluster | prediction | sctype_score | sctype_confidence | diff_exp_cell_markers |
|---|---|---|---|---|
| 0 | Cytotoxic | 3.3581258736158235 | 2.520638196222706 | CCL4,CCL5,CD8A,CD8B,NKG7 |
| 1 | canonical Treg | 1.4430955098410294 | 1.9535556084922565 | CD27,FOXP3,TIGIT |
| 2 | canonical Treg | 0.7823790731528928 | 1.0885935426314584 | CD27,FOXP3,TIGIT |
| 3 | IFN Signature | 0.4133995863316864 | 0.8949390113944568 | IFI6,ISG15,OAS1 |
| 4 | CCR7+ | 1.010106859493755 | 1.9959946748262258 | CCR7,JUNB,TCF7 |
| 5 | CCR7+ | 1.4313989369235969 | 2.786438453818763 | CCR7,JUNB,TCF7 |
| 6 | CCR7+ | 0.28746454956110856 | 0.6401464428637319 | CCR7,TCF7 |
| 7 | Cytotoxic | 2.798557552009433 | 2.1186955573400614 | CCL4,CCL5,CD8A,CD8B,NKG7 |
| 8 | KLRB1+ | 0.8992559328340234 | 1.646537150845831 | CCR6,GPR25,KLRB1,RORC |
| 9 | KLRB1+ | 1.675548512689046 | 3.124456692207059 | CCR6,GPR25,KLRB1,MGAT4A,RORC |
| 10 | canonical Treg | 1.0988892980096865 | 1.502945786281512 | FOXP3,TIGIT |
| 11 | Cycling | 4.2316616825704525 | 3.4008272165033113 | CD74,HIST1H1E,MKI67,STMN1,TUBB,TYMS |
| 12 | canonical Treg | 1.0765041107763713 | 1.4736407217805236 | CD27,FOXP3,TIGIT |
| 13 | Cytotoxic | 2.9011051644643007 | 2.192356371166863 | CCL4,CCL5,CD8A,CD8B,NKG7 |
| 14 | ANXA1 high | 0.1029989098250452 | 0.2532713094300443 | ANXA1 |
| 15 | canonical Treg | 0.8563717750471573 | 1.1854594244761598 | CD27,FOXP3 |
| 16 | CCR7+ | 0.3724674746650708 | 0.7996320751464437 | JUNB |
| 17 | MALAT1 high | 0.024159981259278835 | -0.1604195233537828 | |
| 18 | CCR4/ Helios + | 2.0017923580032497 | 3.9045410189409626 | CCR4,IKZF2,RORA |
| 19 | Cycling | 2.527174771959541 | 2.008804656144139 | CD74,FOS,HIST1H1E,HLA-DPA1,HLA-DRA,HLA-DRB5,TUBB |
Table 2: Table displaying sctype score and differential expressed genes per cell annotation
| Cell type predictions are made using the author-reported cell types. Next to the predictions, the marker genes of the assigned cell type which are differentially expressed in the corresponding cluster are also highlighted (where found). Differentially expressed genes were identified by running the Scanpy rank_genes_groups function with the following settings:Log-fold change cutoff: 1.0, Statistical test: t-test Adjusted p-value cutoff (Benjamini-Hochberg): 0.05 By default, "normalized_counts" layer is used for DE testing. DE genes per cluster are identified separately within each batch, and the results from all batches are summarized at the cluster level. |
1. Metadata information
| Metadata information | Value |
|---|---|
| Polly curated metadata fields are present at dataset level ℹ | Pass |
| Polly curated metadata fields are present at sample level ℹ | Pass |
| Polly curated metadata fields are present in output file ℹ | Pass |
| Custom fields are present in output file ℹ | Pass |
| Publication Link is provided ℹ | Pass |
| Publication Link is valid ℹ | Pass |
| Dataset-Level vs. Sample-Level Metadata: concordance check ℹ | Fail | Accuracy of raw counts availability tag ℹ | Pass |
2. Data Matrix
| Data Matrix | Value |
|---|---|
| Unique Cell Barcodes ℹ | Pass |
| Unique Gene Identifiers ℹ | Pass |
| Embeddings are available ℹ | Fail |
| Gene Identifier Format ℹ | Pass |
| Raw counts are available in output file ℹ | Pass |
| Raw vs Processed Counts are different ℹ | Pass |
| Valid Raw Counts ℹ | Pass |
| Concordance of number of cells in raw and processed counts matrices in output file ℹ | Pass |
| Valid Columns ℹ | Pass |
| Highly Variable Genes is available ℹ | Pass |
| Valid Processed Counts ℹ | Pass |
| UMAP/tSNE Projections are available ℹ | Both present |
| QC Metrics are available ℹ | Pass |
| Reproducibility of Gene Counts ℹ | Pass |
| Reproducibility of UMI Counts ℹ | Pass |
| Cluster information is available ℹ | Pass |
| Number of Clusters ℹ | 20 |
| Minimum genes per cell threshold ℹ | 500 |
| Minimum cells per gene threshold ℹ | 2 |
3. Cell Clusters in umap Embeddings Colored by Samples: Re-Processed and Polly Datasets
Figure 1a: Sample level distribution of clustering pattern of cells with the
help of
umap embeddings on the existing on polly data.
Figure 1b: Sample level distribution of clustering pattern of cells with the
help of umap embeddings on the re - processed data to validate reproducibility of results.
Figure 1a: Sample level distribution of clustering pattern of cells with the help of umap embeddings on the existing on polly data.
Figure 1b: Sample level distribution of clustering pattern of cells with the help of umap embeddings on the re - processed data to validate reproducibility of results.
| The plot visualizes the distribution of samples across various clusters. For both Polly and
reprocessed dataset, these should appear very similar.
Additionally the plot for Polly datasets can be used to understand if there is any batch-effect.
Sample Clustering: If samples are grouped in a diverse manner, where cells from the same sample are not closely clustered together, this suggests no batch effects on samples. Batch Effect Evidence: If the opposite is true, with cells from the same sample clustering together, there might be evidence of batch effects on samples. Biological Variation Check: It's essential to ensure that any batch effects observed are not due to inherent biological differences between samples. Distribution Visualization: The plot also illustrates how samples are spread across different clusters, providing insights into their distribution. Limitation of Reprocessed dataset: Note that using the UMAP/tSNE plot for reprocessed dataset may not be a valid approach to assess batch effects on samples, particularly when dealing with re-processed data primarily focused on reproducibility checks. |
4. Cell Clusters in umap Embeddings Colored by 'Author Cell Types': Comparison Between Polly and Re-Processed Datasets
Figure 2a: Author cell type level distribution of clustering pattern of cells
with the help of umap embeddings on the existing on polly data.
Figure 2b: Author cell type level distribution of clustering pattern of cells
with the help of umap embeddings on the re - processed data to validate reproducibility of
results.
Figure 2a: Author cell type level distribution of clustering pattern of cells with the help of umap embeddings on the existing on polly data.
Figure 2b: Author cell type level distribution of clustering pattern of cells with the help of umap embeddings on the re - processed data to validate reproducibility of results.
|
Cell Type Distribution (author-defined): The plot visualizes the distribution of
author-defined cell types across various clusters. As a quality check, for both Polly and reprocessed
dataset, these should appear very similar.
Cell Type Similarity: UMAP plot also reveals the degree of similarity between different cell types. If cell types A and B are closely clustered, their gene expression patterns are similar, indicating biological similarities between these cell types. |
5. Cell Clusters in umap Embeddings Colored by 'Curated Cell Types': Comparison Between Polly Dataset and Re-Processed Data
Figure 5a: Curated cell type level distribution of clustering pattern of cells
with the help of umap embeddings on the existing on polly data.
Figure 5b: Curated cell type level distribution of clustering pattern of cells
with the help of umap embeddings on the re - processed data to validate reproducibility of
results.
Figure 5a: Curated cell type level distribution of clustering pattern of cells with the help of umap embeddings on the existing on polly data.
Figure 5b: Curated cell type level distribution of clustering pattern of cells with the help of umap embeddings on the re - processed data to validate reproducibility of results.
|
Cell Type Distribution by Elucidata (Curation Experts): The plot visualizes how
curated cell types are distributed among different clusters. As a quality check, For both Polly and
reprocessed dataset, these should appear very similar.
Cell Type Relationships: It shows the proximity of different cell types within the clusters. If cell types A and B cluster closely, it suggests similar gene expression patterns between them, indicating biological similarities between these cell types. |
6. Violin plot visualization for doublet
Figure 5: Sanity check of detected doublets
| To assess the validity of doublet predictions, we plot the distribution of detected genes in predicted doublets v/s singlets per sample (number of genes per count are expected to be typically higher in heterotypic doublets). If doublets are removed the plot only shows the distribution of genes per count in singlets. |
7. Cell Type Frequency Distribution
| Cell type (reported in publication) | Cell type (Polly curated) | Number of cells | |
|---|---|---|---|
| 0 | ["ANXA1 high"] | ["T cell"] | 351 |
| 1 | ["CCR4/ Helios +"] | ["T cell"] | 143 |
| 2 | ["CCR7+"] | ["T cell"] | 3537 |
| 3 | ["Cycling"] | ["stem cell"] | 636 |
| 4 | ["Cytotoxic"] | ["cytotoxic T cell"] | 2765 |
| 5 | ["IFN Signature"] | ["T cell"] | 1140 |
| 6 | ["KLRB1+"] | ["T cell"] | 1217 |
| 7 | ["MALAT1 high"] | ["T cell"] | 177 |
| 8 | ["canonical Treg"] | ["regulatory T cell"] | 3988 |
Table 2: Table displaying author cell types, curated cell types and the number of cells for each cell-type
| Authors frequently supply cell types that may not adhere to ontological standards or utilize abbreviations and marker gene names. These are substituted with ontological terms. The table offers insight into the degree of alignment between the ontological terms and the terms provided by the authors. |
1. Expression of Marker Genes Across Cell Types
Figure 1: The dot plot showcases the expression levels (often represented by dot size) and prevalence (often represented by dot color intensity) of specific marker genes across different cell types.
| Marker genes that are predominantly expressed in specific cell types validate the identified cell populations and help in characterizing and annotating them. |
2. Expression of Marker Genes Across Clusters
Figure 2: The dot plot showcases the expression levels (often represented by dot size) and prevalence (often represented by dot color intensity) of specific marker genes across different clusters.
| This visualization aids in understanding the heterogeneity within the dataset and can hint at different cellular states or subtypes within a cell type. |