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.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-014-0088-0) contains supplementary material, which is available to authorized users.

Lactose transport system of Streptococcus thermophilus. The role of histidine residues

The lactose transport protein (LacS) of Streptococcus thermophilus is a chimeric protein consisting of an amino-terminal carrier domain and a carboxyl-terminal phosphoenolpyruvate:sugar phosphotransferase system (PTS) IIA protein domain. The histidine residues of LacS were changed individually into glutamine or arginine residues. Of the 11 histidine residues present in LacS, only the His-376 substitution in the carrier domain significantly affected sugar transport. The region around His-376 was found to exhibit sequence similarity to the region around His-322 of the lactose transport protein (LacY) of Escherichia coli, which has been implicated in sugar binding and in coupling of sugar and H+ transport. The H376Q mutation resulted in a reduced rate of uptake and altered affinity for lactose (beta-galactoside), melibiose (alpha-galactoside), and the lactose analog methyl-beta-D-thiogalactopyranoside. Similarly, the extent of accumulation of the galactosides by cells expressing LacS(H376Q) was highly reduced in comparison to cells bearing the wild-type protein. Nonequilibrium exchange of lactose and methyl-beta-D-thiogalactopyranoside by the H376Q mutant was approximately 2-fold reduced in comparison to the activity of the wild-type transport protein. The data indicate that His-376 is involved in sugar recognition and is important, but not essential, for the cotransport of protons and galactosides. The carboxyl-terminal domain of LacS contains 2 histidine residues (His-537 and His-552) that are conserved in seven homologous IIA protein(s) (domains) of PTSs. P-enolpyruvate-dependent phosphorylation of wild-type LacS, but not of the mutant H552Q, was demonstrated using purified Enzyme I and HPr, the general energy coupling proteins of the PTS, and inside-out membrane vesicles isolated from E. coli in which the lactose transport gene was expressed. The His-537 and His-552 mutations did not affect transport activity when the corresponding genes were expressed in E. coli.