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.

Circulating miRNAs as Potential Biomarkers for Celiac Disease Development.

Background & Aims: Celiac disease (CeD), an immune-mediated disease with enteropathy triggered by gluten, affects ~1% of the general European population. Currently, there are no biomarkers to predict CeD development. MicroRNAs (miRNAs) are short RNAs involved in post-transcriptional gene regulation, and certain disease- and stage-specific miRNA profiles have been found previously. We aimed to investigate whether circulating miRNAs can predict the development of CeD. Methods: Using next-generation miRNA-sequencing, we determined miRNAs in >200 serum samples from 53 participants of the PreventCD study, of whom 33 developed CeD during follow-up. Following study inclusion at 3 months of age, samples were drawn at predefined ages, diagnosis (first anti-transglutaminase antibody (TGA) positivity or diagnostic biopsy) and after the start of a gluten-free diet (GFD). This allowed identification of circulating miRNAs that are deregulated before TGA positivity. For validation of the biomarkers for CeD and GFD response, two additional cohorts were included in subsequent meta-analyses. Additionally, miRNAs were measured in duodenal biopsies in a case-control cohort. Results: 53 circulating miRNAs were increased (27) or decreased (26) in CeD versus controls. We assessed specific trends in these individual miRNAs in the PreventCD cohort by grouping the pre-diagnostic samples of the CeD patients (all had negative TGA) by how close to seroconversion (first sample positive TGA) the samples were taken. 8/53 miRNAs differed significantly between controls and samples taken <1 year before TGA positivity: miR-21-3p, miR-374a-5p, 144-3p, miR-500a-3p, miR-486-3p let-7d-3p, let-7e-5p and miR-3605-3p. 6/26 downregulated miRNAs reconstituted upon GFD, including miR-150-5p/-3p, whereas no upregulated miRNAs were downregulated upon GFD. 15/53 biomarker candidates also differed between CeD biopsies and controls, with a concordant direction, indicating that these circulating miRNAs might originate from the intestine. Conclusions: We identified 53 circulating miRNAs that are potential early biomarkers for CeD, of which several can be detected more than a year before TGA positivity and some start to normalize upon GFD.

A Combined mRNA- and miRNA-Sequencing Approach Reveals miRNAs as Potential Regulators of the Small Intestinal Transcriptome in Celiac Disease.

Celiac disease (CeD) is triggered by gluten and results in inflammation and villous atrophy of the small intestine. We aimed to explore the role of miRNA-mediated deregulation of the transcriptome in CeD. Duodenal biopsies of CeD patients (n = 33) and control subjects (n = 10) were available for miRNA-sequencing, with RNA-sequencing also available for controls (n = 5) and CeD (n = 6). Differential expression analysis was performed to select CeD-associated miRNAs and genes. MiRNA‒target transcript pairs selected from public databases that also displayed a strong negative expression correlation in the current dataset (R < -0.7) were used to construct a CeD miRNA‒target transcript interaction network. The network includes 2030 miRNA‒target transcript interactions, including 423 experimentally validated pairs. Pathway analysis found that interactions are involved in immune-related pathways (e.g., interferon signaling) and metabolic pathways (e.g., lipid metabolism). The network includes 13 genes previously prioritized to be causally deregulated by CeD-associated genomic variants, including STAT1. CeD-associated miRNAs might play a role in promoting inflammation and decreasing lipid metabolism in the small intestine, thereby contributing unbalanced cell turnover in the intestinal crypt. Some CeD-associated miRNAs deregulate genes that are also affected by genomic CeD-risk variants, adding an additional layer of complexity to the deregulated transcriptome in CeD.

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.

Small RNA sequencing reveals a comprehensive miRNA signature of BRCA1-associated high-grade serous ovarian cancer.

AIMS: BRCA1 mutation carriers are at increased risk of developing high-grade serous ovarian cancer (HGSOC), a malignancy that originates from fallopian tube epithelium. We aimed to identify differentially expressed known and novel miRNAs in BRCA1-associated HGSOC. METHODS: Small RNA sequencing was performed on eight normal tubal and five HGSOC samples of BRCA1 carriers. Differential expression of a subset of known and novel miRNAs was validated by qRT-PCR on the samples used for small RNA sequencing and a second sample cohort comprising normal and HGSOC tissue of matched BRCA1 and non-BRCA carriers. Data from The Cancer Genome Atlas were used to determine the clinical relevance of the validated differentially expressed miRNAs. RESULTS: 59 known and 20 novel miRNAs showed a significant >fourfold expression difference between normal tubal tissue and HGSOC. qRT-PCR validation confirmed a significant difference in expression levels for 10 out of 11 known miRNAs. Upregulation of two novel miRNAs could not be confirmed. Interestingly, for seven miRNAs a significant increase in expression was observed when comparing normal tubal tissue of postmenopausal women with premenopausal women. Expression levels of miR-145-5p significantly increased with International Federation of Gynecology and Obstetrics stage, while the expression levels of the other nine validated miRNAs were not associated with clinical characteristics. CONCLUSIONS: We report a comprehensive expression signature including both known and novel miRNAs of BRCA1-associated HGSOC. Comparison with previous profiling studies showed a good overlap and a large number of miRNAs not reported to be differentially expressed in HGSOC before underscoring the importance of this study.

Expression profiles of long non-coding RNAs located in autoimmune disease-associated regions reveal immune cell-type specificity.

BACKGROUND: Although genome-wide association studies (GWAS) have identified hundreds of variants associated with a risk for autoimmune and immune-related disorders (AID), our understanding of the disease mechanisms is still limited. In particular, more than 90% of the risk variants lie in non-coding regions, and almost 10% of these map to long non-coding RNA transcripts (lncRNAs). lncRNAs are known to show more cell-type specificity than protein-coding genes. METHODS: We aimed to characterize lncRNAs and protein-coding genes located in loci associated with nine AIDs which have been well-defined by Immunochip analysis and by transcriptome analysis across seven populations of peripheral blood leukocytes (granulocytes, monocytes, natural killer (NK) cells, B cells, memory T cells, naive CD4(+) and naive CD8(+) T cells) and four populations of cord blood-derived T-helper cells (precursor, primary, and polarized (Th1, Th2) T-helper cells). RESULTS: We show that lncRNAs mapping to loci shared between AID are significantly enriched in immune cell types compared to lncRNAs from the whole genome (α <0.005). We were not able to prioritize single cell types relevant for specific diseases, but we observed five different cell types enriched (α <0.005) in five AID (NK cells for inflammatory bowel disease, juvenile idiopathic arthritis, primary biliary cirrhosis, and psoriasis; memory T and CD8(+) T cells in juvenile idiopathic arthritis, primary biliary cirrhosis, psoriasis, and rheumatoid arthritis; Th0 and Th2 cells for inflammatory bowel disease, juvenile idiopathic arthritis, primary biliary cirrhosis, psoriasis, and rheumatoid arthritis). Furthermore, we show that co-expression analyses of lncRNAs and protein-coding genes can predict the signaling pathways in which these AID-associated lncRNAs are involved. CONCLUSIONS: The observed enrichment of lncRNA transcripts in AID loci implies lncRNAs play an important role in AID etiology and suggests that lncRNA genes should be studied in more detail to interpret GWAS findings correctly. The co-expression results strongly support a model in which the lncRNA and protein-coding genes function together in the same pathways.

Correlation of genetic risk and messenger RNA expression in a Th17/IL23 pathway analysis in inflammatory bowel disease.

BACKGROUND: The Th17/IL23 pathway has both genetically and biologically been implicated in the pathogenesis of the inflammatory bowel diseases (IBD), Crohn's disease, and ulcerative colitis. So far, it is unknown whether and how associated risk variants affect expression of the genes encoding for Th17/IL23 pathway proteins. METHODS: Ten IBD-associated SNPs residing near Th17/IL23 genes were used to construct a genetic risk model in 753 Dutch IBD cases and 1045 controls. In an independent cohort of 40 Crohn's disease, 40 ulcerative colitis, and 40 controls, the genetic risk load and presence of IBD were correlated to quantitative PCR-generated messenger RNA (mRNA) expression of 9 representative Th17/IL23 genes in both unstimulated and PMA/CaLo stimulated peripheral blood mononuclear cells. In 1240 individuals with various immunological diseases with whole genome genotype and mRNA-expression data, we also assessed correlation between genetic risk load and differential mRNA expression and sought for SNPs affecting expression of all currently known Th17/IL23 pathway genes (cis-expression quantitative trait locus). RESULTS: The presence of IBD, but not the genetic risk load, was correlated to differential mRNA expression for IL6 in unstimulated peripheral blood mononuclear cells and to IL23A and RORC in response to stimulation. The cis-expression quantitative trait locus analysis showed little evidence for correlation between genetic risk load and mRNA expression of Th17/IL23 genes, because we identified for only 2 of 22 Th17/IL23 genes a cis-expression quantitative trait locus single nucleotide polymorphism that is also associated to IBD (STAT3 and CCR6). CONCLUSIONS: Our results suggest that only the presence of IBD and not the genetic risk load alters mRNA expression levels of IBD-associated Th17/IL23 genes.