Single-Cell RNA Sequencing of Peripheral Blood Mononuclear Cells From Pediatric Coeliac Disease Patients Suggests Potential Pre-Seroconversion Markers

Celiac Disease (CeD) is a complex immune disorder involving villous atrophy in the small intestine that is triggered by gluten intake. Current CeD diagnosis is based on late-stage pathophysiological parameters such as detection of specific antibodies in blood and histochemical detection of villus atrophy and lymphocyte infiltration in intestinal biopsies. To date, no early onset biomarkers are available that would help prevent widespread villous atrophy and severe symptoms and co-morbidities. To search for novel CeD biomarkers, we used single-cell RNA sequencing (scRNAseq) to investigate PBMC samples from 11 children before and after seroconversion for CeD and 10 control individuals matched for age, sex and HLA-genotype. We generated scRNAseq profiles of 9559 cells and identified the expected major cellular lineages. Cell proportions remained stable across the different timepoints and health conditions, but we observed differences in gene expression profiles in specific cell types when comparing patient samples before and after disease development and comparing patients with controls. Based on the time when transcripts were differentially expressed, we could classify the deregulated genes as biomarkers for active CeD or as potential pre-diagnostic markers. Pathway analysis showed that active CeD biomarkers display a transcriptional profile associated with antigen activation in CD4+ T cells, whereas NK cells express a subset of biomarker genes even before CeD diagnosis. Intersection of biomarker genes with CeD-associated genetic risk loci pinpointed genetic factors that might play a role in CeD onset. Investigation of potential cellular interaction pathways of PBMC cell subpopulations highlighted the importance of TNF pathways in CeD. Altogether, our results pinpoint genes and pathways that are altered prior to and during CeD onset, thereby identifying novel potential biomarkers for CeD diagnosis in blood

Mass cytometric analysis unveils a disease-specific immune cell network in the bone marrow in acquired aplastic anemia

Idiopathic acquired aplastic anemia (AA) is considered an immune-mediated syndrome of bone marrow failure since approximately 70% of patients respond to immunosuppressive therapy (IST) consisting of a course of anti-thymocyte globulin (ATG) followed by long-term use of ciclosporin. However, the immune response that underlies the pathogenesis of AA remains poorly understood. In this study, we applied high-dimensional mass cytometry on bone marrow aspirates of AA patients pre-ATG, AA patients post-ATG and healthy donors to decipher which immune cells may be implicated in the pathogenesis of AA. We show that the bone marrow of AA patients features an immune cell composition distinct from healthy donors, with significant differences in the myeloid, B-cell, CD4+ and CD8+ T-cells lineages. Specifically, we discovered that AA pre-ATG is characterized by a disease-specific immune cell network with high frequencies of CD16+ myeloid cells, CCR6++ B-cells, Th17-like CCR6+ memory CD4+ T-cells, CD45RA+CCR7+CD38+ CD8+ T-cells and KLRG1+ terminally differentiated effector memory (EMRA) CD8+ T-cells, compatible with a state of chronic inflammation. Successful treatment with IST strongly reduced the levels of CD16+ myeloid cells and showed a trend toward normalization of the frequencies of CCR6++ B-cells, CCR6+ memory CD4+ T-cells and KLRG1+EMRA CD8+ T-cells. Altogether, our study provides a unique overview of the immune landscape in bone marrow in AA at a single-cell level and proposes CCR6 as a potential new therapeutic target in AA

Systematic Prioritization of Candidate Genes in Disease Loci Identifies TRAFD1 as a Master Regulator of IFN gamma Signaling in Celiac Disease

Celiac disease (CeD) is a complex T cell-mediated enteropathy induced by gluten. Although genome-wide association studies have identified numerous genomic regions associated with CeD, it is difficult to accurately pinpoint which genes in these loci are most likely to cause CeD. We used four different in silico approaches-Mendelian randomization inverse variance weighting, COLOC, LD overlap, and DEPICT-to integrate information gathered from a large transcriptomics dataset. This identified 118 prioritized genes across 50 CeD-associated regions. Co-expression and pathway analysis of these genes indicated an association with adaptive and innate cytokine signaling and T cell activation pathways. Fifty-one of these genes are targets of known drug compounds or likely druggable genes, suggesting that our methods can be used to pinpoint potential therapeutic targets. In addition, we detected 172 gene combinations that were affected by our CeD-prioritized genes in trans. Notably, 41 of these trans-mediated genes appear to be under control of one master regulator, TRAF-type zinc finger domain containing 1 (TRAFD1), and were found to be involved in interferon (IFN)gamma signaling and MHC I antigen processing/presentation. Finally, we performed in vitro experiments in a human monocytic cell line that validated the role of TRAFD1 as an immune regulator acting in trans. Our strategy confirmed the role of adaptive immunity in CeD and revealed a genetic link between CeD and IFN gamma signaling as well as with MHC I antigen processing, both major players of immune activation and CeD pathogenesis

Potential impact of celiac disease genetic risk factors on T cell receptor signaling in gluten-specific CD4+ T cells

Celiac disease is an auto-immune disease in which an immune response to dietary gluten leads to inflammation and subsequent atrophy of small intestinal villi, causing severe bowel discomfort and malabsorption of nutrients. The major instigating factor for the immune response in celiac disease is the activation of gluten-specific CD4+ T cells expressing T cell receptors that recognize gluten peptides presented in the context of HLA-DQ2 and DQ8. Here we provide an in-depth characterization of 28 gluten-specific T cell clones. We assess their transcriptional and epigenetic response to T cell receptor stimulation and link this to genetic factors associated with celiac disease. Gluten-specific T cells have a distinct transcriptional profile that mostly resembles that of Th1 cells but also express cytokines characteristic of other types of T-helper cells. This transcriptional response appears not to be regulated by changes in chromatin state, but rather by early upregulation of transcription factors and non-coding RNAs that likely orchestrate the subsequent activation of genes that play a role in immune pathways. Finally, integration of chromatin and transcription factor binding profiles suggest that genes activated by T cell receptor stimulation of gluten‑specific T cells may be impacted by genetic variation at several genetic loci associated with celiac disease.</p

T cell receptor cross-reactivity between gliadin and bacterial peptides in celiac disease

The human leukocyte antigen (HLA) locus is strongly associated with T cell-mediated autoimmune disorders. HLA-DQ2.5-mediated celiac disease (CeD) is triggered by the ingestion of gluten, although the relative roles of genetic and environmental risk factors in CeD is unclear. Here we identify microbially derived mimics of gliadin epitopes and a parental bacterial protein that is naturally processed by antigen-presenting cells and activated gliadin reactive HLA-DQ2.5-restricted T cells derived from CeD patients. Crystal structures of T cell receptors in complex with HLA-DQ2.5 bound to two distinct bacterial peptides demonstrate that molecular mimicry underpins cross-reactivity toward the gliadin epitopes. Accordingly, gliadin reactive T cells involved in CeD pathogenesis cross-react with ubiquitous bacterial peptides, thereby suggesting microbial exposure as a potential environmental factor in CeD

Single-Cell Analysis of Refractory Celiac Disease Demonstrates Inter- and Intra-Patient Aberrant Cell Heterogeneity

Background & Aims: Refractory celiac disease type II (RCDII) is a rare indolent lymphoma in the small intestine characterized by a clonally expanded intraepithelial intracellular CD3+surfaceCD3-CD7+CD56- aberrant cell population. However, RCDII pathogenesis is ill-defined. Here, we aimed at single-cell characterization of the innate and adaptive immune system in RCDII. Methods: Paired small intestinal and blood samples from 12 RCDII patients and 6 healthy controls were assessed by single-cell mass cytometry with a 39–cell surface marker antibody panel, designed to capture heterogeneity of the innate and adaptive immune system. A second single-cell mass cytometry panel that included transcription factors and immune checkpoints was used for analysis of paired samples from 5 RCDII patients. Single-cell RNA sequencing analysis was performed on duodenal samples from 2 RCDII patients. Finally, we developed a 40-marker imaging mass cytometry antibody panel to evaluate cell–cell interactions in duodenal biopsy specimens of RCDII patients. Results: We provide evidence for intertumoral and intratumoral cell heterogeneity within the duodenal and peripheral aberrant cell population present in RCDII. Phenotypic discrepancy was observed between peripheral and duodenal aberrant cells. In addition, we observed that part of the aberrant cell population proliferated and observed co-localization of aberrant cells with CD163+ antigen-presenting cells (APCs) in situ. In addition, we observed phenotypic discrepancy between peripheral and duodenal aberrant cells. Conclusions: Novel high-dimensional single-cell technologies show substantial intertumoral and intratumoral heterogeneity in the aberrant cell population in RCDII. This may underlie variability in refractory disease status between patients and responsiveness to therapy, pointing to the need for personalized therapy in RCDII based on patient-specific immune profiles

Single-Cell Analysis of Refractory Celiac Disease Demonstrates Inter- and Intra-Patient Aberrant Cell Heterogeneity

Background &amp; Aims: Refractory celiac disease type II (RCDII) is a rare indolent lymphoma in the small intestine characterized by a clonally expanded intraepithelial intracellular CD3+surfaceCD3-CD7+CD56- aberrant cell population. However, RCDII pathogenesis is ill-defined. Here, we aimed at single-cell characterization of the innate and adaptive immune system in RCDII. Methods: Paired small intestinal and blood samples from 12 RCDII patients and 6 healthy controls were assessed by single-cell mass cytometry with a 39–cell surface marker antibody panel, designed to capture heterogeneity of the innate and adaptive immune system. A second single-cell mass cytometry panel that included transcription factors and immune checkpoints was used for analysis of paired samples from 5 RCDII patients. Single-cell RNA sequencing analysis was performed on duodenal samples from 2 RCDII patients. Finally, we developed a 40-marker imaging mass cytometry antibody panel to evaluate cell–cell interactions in duodenal biopsy specimens of RCDII patients. Results: We provide evidence for intertumoral and intratumoral cell heterogeneity within the duodenal and peripheral aberrant cell population present in RCDII. Phenotypic discrepancy was observed between peripheral and duodenal aberrant cells. In addition, we observed that part of the aberrant cell population proliferated and observed co-localization of aberrant cells with CD163+ antigen-presenting cells (APCs) in situ. In addition, we observed phenotypic discrepancy between peripheral and duodenal aberrant cells. Conclusions: Novel high-dimensional single-cell technologies show substantial intertumoral and intratumoral heterogeneity in the aberrant cell population in RCDII. This may underlie variability in refractory disease status between patients and responsiveness to therapy, pointing to the need for personalized therapy in RCDII based on patient-specific immune profiles.Pattern Recognition and Bioinformatic

Peptides generated from the 26-mer and 33-mer.

<p>Numbers and letters refer to peaks in <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0128065#pone.0128065.g002" target="_blank">Fig 2</a></p><p>Sequences that contain epitopes are underlined</p><p>*peptides shorter than eight amino acids are considered too short to contain an epitope</p><p>Peptides generated from the 26-mer and 33-mer.</p