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.

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.

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.

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.


3. umap plots for categorical metadata

Figure 3: The umap visualization represents cells in a reduced dimensional space, with colors indicating various categorical attributes.



4. umap plots for Polly curated metadata

Figure 4: This umap visualization represents cells in a reduced dimensional space, with colors indicating the Polly curated fields.


5. Sunburst plots for metadata fields

Figure 5: A Sunburst plot illustrating the distribution of data. It reflects user-defined custom fields if specified; otherwise, it represents standard fields.