Constantly decreasing costs of high-throughput profiling on many molecular levels generate vast amounts of multi-omics data. Studying one biomedical question on two or more omic levels provides ...deeper insights into underlying molecular processes or disease pathophysiology. For the majority of multi-omics data projects, the data analysis is performed level-wise, followed by a combined interpretation of results. Hence the full potential of integrated data analysis is not leveraged yet, presumably due to the complexity of the data and the lacking toolsets. We propose a versatile approach, to perform a multi-level fully integrated analysis: The Knowledge guIded Multi-Omics Network inference approach, KiMONo ( https://github.com/cellmapslab/kimono ). KiMONo performs network inference by using statistical models for combining omics measurements coupled to a powerful knowledge-guided strategy exploiting prior information from existing biological sources. Within the resulting multimodal network, nodes represent features of all input types e.g. variants and genes while edges refer to knowledge-supported and statistically derived associations. In a comprehensive evaluation, we show that our method is robust to noise and exemplify the general applicability to the full spectrum of multi-omics data, demonstrating that KiMONo is a powerful approach towards leveraging the full potential of data sets for detecting biomarker candidates.
N6-methyladenosine (m6A) and N6,2′-O-dimethyladenosine (m6Am) are abundant mRNA modifications that regulate transcript processing and translation. The role of both, here termed m6A/m, in the stress ...response in the adult brain in vivo is currently unknown. Here, we provide a detailed analysis of the stress epitranscriptome using m6A/m-seq, global and gene-specific m6A/m measurements. We show that stress exposure and glucocorticoids region and time specifically alter m6A/m and its regulatory network. We demonstrate that deletion of the methyltransferase Mettl3 or the demethylase Fto in adult neurons alters the m6A/m epitranscriptome, increases fear memory, and changes the transcriptome response to fear and synaptic plasticity. Moreover, we report that regulation of m6A/m is impaired in major depressive disorder patients following glucocorticoid stimulation. Our findings indicate that brain m6A/m represents a novel layer of complexity in gene expression regulation after stress and that dysregulation of the m6A/m response may contribute to the pathophysiology of stress-related psychiatric disorders.
Display omitted
•m6A/m mRNA methylation in the adult mouse brain is regulated by stress•m6A/m mRNA regulation is brain region, time, and gene specific•Mettl3 and Fto cKO alter m6A/m, fear memory, expression, and synaptic plasticity•The m6A/m glucocorticoid response is impaired in major depressive disorder patients
Engel et al. demonstrate a region- and time-dependent role of brain m6A/m methylation in stress-response regulation. Manipulating m6A/m alters fear memory, transcriptome response, and synaptic plasticity. Altered m6A/m dynamics in depressed patients suggest importance of m6A/m modifications for stress-related psychiatric disorders.
Microarray technologies are established approaches for high throughput gene expression, methylation and genotyping analysis. An accurate mapping of the array probes is essential to generate reliable ...biological findings. However, manufacturers of the microarray platforms typically provide incomplete and outdated annotation tables, which often rely on older genome and transcriptome versions that differ substantially from up-to-date sequence databases. Here, we present the Re-Annotator, a re-annotation pipeline for microarray probe sequences. It is primarily designed for gene expression microarrays but can also be adapted to other types of microarrays. The Re-Annotator uses a custom-built mRNA reference database to identify the positions of gene expression array probe sequences. We applied Re-Annotator to the Illumina Human-HT12 v4 microarray platform and found that about one quarter (25%) of the probes differed from the manufacturer's annotation. In further computational experiments on experimental gene expression data, we compared Re-Annotator to another probe re-annotation tool, ReMOAT, and found that Re-Annotator provided an improved re-annotation of microarray probes. A thorough re-annotation of probe information is crucial to any microarray analysis. The Re-Annotator pipeline is freely available at http://sourceforge.net/projects/reannotator along with re-annotated files for Illumina microarrays HumanHT-12 v3/v4 and MouseRef-8 v2.
Although gene expression profiles in peripheral blood in major depression are not likely to identify genes directly involved in the pathomechanism of affective disorders, they may serve as biomarkers ...for this disorder. As previous studies using baseline gene expression profiles have provided mixed results, our approach was to use an in vivo dexamethasone challenge test and to compare glucocorticoid receptor (GR)-mediated changes in gene expression between depressed patients and healthy controls. Whole genome gene expression data (baseline and following GR-stimulation with 1.5 mg dexamethasone p.o.) from two independent cohorts were analyzed to identify gene expression pattern that would predict case and control status using a training (N=18 cases/18 controls) and a test cohort (N=11/13). Dexamethasone led to reproducible regulation of 2670 genes in controls and 1151 transcripts in cases. Several genes, including FKBP5 and DUSP1, previously associated with the pathophysiology of major depression, were found to be reliable markers of GR-activation. Using random forest analyses for classification, GR-stimulated gene expression outperformed baseline gene expression as a classifier for case and control status with a correct classification of 79.1 vs 41.6% in the test cohort. GR-stimulated gene expression performed best in dexamethasone non-suppressor patients (88.7% correctly classified with 100% sensitivity), but also correctly classified 77.3% of the suppressor patients (76.7% sensitivity), when using a refined set of 19 genes. Our study suggests that in vivo stimulated gene expression in peripheral blood cells could be a promising molecular marker of altered GR-functioning, an important component of the underlying pathology, in patients suffering from depressive episodes.
It is unknown whether indoles, metabolites of tryptophan that are derived entirely from bacterial metabolism in the gut, are associated with symptoms of depression and anxiety. Serum samples ...(baseline, 12 weeks) were drawn from participants (n = 196) randomized to treatment with cognitive behavioral therapy (CBT), escitalopram, or duloxetine for major depressive disorder. Baseline indoxyl sulfate abundance was positively correlated with severity of psychic anxiety and total anxiety and with resting state functional connectivity to a network that processes aversive stimuli (which includes the subcallosal cingulate cortex (SCC-FC), bilateral anterior insula, right anterior midcingulate cortex, and the right premotor areas). The relation between indoxyl sulfate and psychic anxiety was mediated only through the metabolite's effect on the SCC-FC with the premotor area. Baseline indole abundances were unrelated to post-treatment outcome measures, and changes in symptoms were not correlated with changes in indole concentrations. These results suggest that CBT and antidepressant medications relieve anxiety via mechanisms unrelated to modulation of indoles derived from gut microbiota; it remains possible that treatment-related improvement stems from their impact on other aspects of the gut microbiome. A peripheral gut microbiome-derived metabolite was associated with altered neural processing and with psychiatric symptom (anxiety) in humans, which provides further evidence that gut microbiome disruption can contribute to neuropsychiatric disorders that may require different therapeutic approaches. Given the exploratory nature of this study, findings should be replicated in confirmatory studies.Clinical trial NCT00360399 "Predictors of Antidepressant Treatment Response: The Emory CIDAR" https://clinicaltrials.gov/ct2/show/NCT00360399 .
Genome-wide association studies (GWAS) identify genetic variants associated with traits or diseases. GWAS never directly link variants to regulatory mechanisms. Instead, the functional annotation of ...variants is typically inferred by post hoc analyses. A specific class of deep learning-based methods allows for the prediction of regulatory effects per variant on several cell type-specific chromatin features. We here describe "DeepWAS", a new approach that integrates these regulatory effect predictions of single variants into a multivariate GWAS setting. Thereby, single variants associated with a trait or disease are directly coupled to their impact on a chromatin feature in a cell type. Up to 61 regulatory SNPs, called dSNPs, were associated with multiple sclerosis (MS, 4,888 cases and 10,395 controls), major depressive disorder (MDD, 1,475 cases and 2,144 controls), and height (5,974 individuals). These variants were mainly non-coding and reached at least nominal significance in classical GWAS. The prediction accuracy was higher for DeepWAS than for classical GWAS models for 91% of the genome-wide significant, MS-specific dSNPs. DSNPs were enriched in public or cohort-matched expression and methylation quantitative trait loci and we demonstrated the potential of DeepWAS to generate testable functional hypotheses based on genotype data alone. DeepWAS is available at https://github.com/cellmapslab/DeepWAS.
Substantial sex differences have been reported in the physiological response to stress at multiple levels, including the release of the stress hormone, cortisol. Here, we explore the genomic variants ...in 93 females and 196 males regulating the initial transcriptional response to cortisol via glucocorticoid receptor (GR) activation. Gene expression levels in peripheral blood were obtained before and after GR-stimulation with the selective GR agonist dexamethasone to identify differential expression following GR-activation. Sex stratified analyses revealed that while the transcripts responsive to GR-stimulation were mostly overlapping between males and females, the quantitative trait loci (eQTLs) regulation differential transcription to GR-stimulation was distinct. Sex-stratified eQTL SNPs (eSNPs) were located in different functional genomic elements and sex-stratified transcripts were enriched within postmortem brain transcriptional profiles associated with Major Depressive Disorder (MDD) specifically in males and females in the cingulate cortex. Female eSNPs were enriched among SNPs linked to MDD in genome-wide association studies. Finally, transcriptional sensitive genetic profile scores derived from sex-stratified eSNPS regulating differential transcription to GR-stimulation were predictive of depression status and depressive symptoms in a sex-concordant manner in a child and adolescent cohort (n = 584). These results suggest the potential of eQTLs regulating differential transcription to GR-stimulation as biomarkers of sex-specific biological risk for stress-related psychiatric disorders.
Epigenetic mechanisms may play a major role in the biological embedding of early-life stress (ELS). One proposed mechanism is that glucocorticoid (GC) release following ELS exposure induces ...long-lasting alterations in DNA methylation (DNAm) of important regulatory genes of the stress response. Here, we investigate the dynamics of GC-dependent methylation changes in key regulatory regions of the FKBP5 locus in which ELS-associated DNAm changes have been reported.
We repeatedly measured DNAm in human peripheral blood samples from 2 independent cohorts exposed to the GC agonist dexamethasone (DEX) using a targeted bisulfite sequencing approach, complemented by data from Illumina 450K arrays. We detected differentially methylated CpGs in enhancers co-localizing with GC receptor binding sites after acute DEX treatment (1 h, 3 h, 6 h), which returned to baseline levels within 23 h. These changes withstood correction for immune cell count differences. While we observed main effects of sex, age, body mass index, smoking, and depression symptoms on FKBP5 methylation levels, only the functional FKBP5 SNP (rs1360780) moderated the dynamic changes following DEX. This genotype effect was observed in both cohorts and included sites previously shown to be associated with ELS.
Our study highlights that DNAm levels within regulatory regions of the FKBP5 locus show dynamic changes following a GC challenge and suggest that factors influencing the dynamics of this regulation may contribute to the previously reported alterations in DNAm associated with current and past ELS exposure.
Cigarette smoking has severe adverse health consequences in adults and in the offspring of mothers who smoke during pregnancy. One of the most widely reported effects of smoking during pregnancy is ...reduced birth weight which is in turn associated with chronic disease in adulthood. Epigenome-wide association studies have revealed that smokers show a characteristic "smoking methylation pattern", and recent authors have proposed that DNA methylation mediates the impact of maternal smoking on birth weight. The aims of the present study were to replicate previous reports that methylation mediates the effect of maternal smoking on birth weight, and for the first time to investigate whether the observed mediation effects are sex-specific in order to account for known sex-specific differences in methylation levels.
Methylation levels in the cord blood of 313 newborns were determined using the Illumina HumanMethylation450K Beadchip. A total of 5,527 CpG sites selected on the basis of evidence from the literature were tested. To determine whether the observed association between maternal smoking and birth weight was attributable to methylation, mediation analyses were performed for significant CpG sites. Separate analyses were then performed in males and females.
Following quality control, 282 newborns eventually remained in the analysis. A total of 25 mothers had smoked consistently throughout the pregnancy. The birthweigt of newborns whose mothers had smoked throughout pregnancy was reduced by >200g. After correction for multiple testing, 30 CpGs showed differential methylation in the maternal smoking subgroup including top "smoking methylation pattern" genes AHRR, MYO1G, GFI1, CYP1A1, and CNTNAP2. The effect of maternal smoking on birth weight was partly mediated by the methylation of cg25325512 (PIM1); cg25949550 (CNTNAP2); and cg08699196 (ITGB7). Sex-specific analyses revealed a mediating effect for cg25949550 (CNTNAP2) in male newborns.
The present data replicate previous findings that methylation can mediate the effect of maternal smoking on birth weight. The analysis of sex-dependent mediation effects suggests that the sex of the newborn may have an influence. Larger studies are warranted to investigate the role of both the identified differentially methylated loci and the sex of the newborn in mediating the association between maternal smoking during pregnancy and birth weight.
Early-life adversity is an important risk factor for major depressive disorder (MDD) and schizophrenia (SCZ) that interacts with genetic factors to confer disease risk through mechanisms that are ...still insufficiently understood. One downstream effect of early-life adversity is the activation of glucocorticoid receptor (GR)-dependent gene networks that drive acute and long-term adaptive behavioral and cellular responses to stress. We have previously shown that genetic variants that moderate GR-induced gene transcription (GR-response eSNPs) are significantly enriched among risk variants from genome-wide association studies (GWASs) for MDD and SCZ. Here, we show that the 63 transcripts regulated by these disease-associated functional genetic variants form a tight glucocorticoid-responsive co-expression network (termed GCN). We hypothesized that changes in the correlation structure of this GCN may contribute to early-life adversity-associated disease risk. Therefore, we analyzed the effects of different qualities of social support and stress throughout life on GCN formation across distinct brain regions using a translational mouse model. We observed that different qualities of social experience substantially affect GCN structure in a highly brain region-specific manner. GCN changes were predominantly found in two functionally interconnected regions, the ventral hippocampus and the hypothalamus, two brain regions previously shown to be of relevance for the stress response, as well as psychiatric disorders. Overall, our results support the hypothesis that a subset of genetic variants may contribute to risk for MDD and SCZ by altering circuit-level effects of early and adult social experiences on GCN formation and structure.