Maternal dysglycaemia, changes in the infant's epigenome modified with a diet and physical activity intervention in pregnancy: Secondary analysis of a randomised control trial

BACKGROUND

Higher maternal plasma glucose (PG) concentrations, even below gestational diabetes mellitus (GDM) thresholds, are associated with adverse offspring outcomes, with DNA methylation proposed as a mediating mechanism. Here, we examined the relationships between maternal dysglycaemia at 24 to 28 weeks' gestation and DNA methylation in neonates and whether a dietary and physical activity intervention in pregnant women with obesity modified the methylation signatures associated with maternal dysglycaemia.

METHODS AND FINDINGS

We investigated 557 women, recruited between 2009 and 2014 from the UK Pregnancies Better Eating and Activity Trial (UPBEAT), a randomised controlled trial (RCT), of a lifestyle intervention (low glycaemic index (GI) diet plus physical activity) in pregnant women with obesity (294 contol, 263 intervention). Between 27 and 28 weeks of pregnancy, participants had an oral glucose (75 g) tolerance test (OGTT), and GDM diagnosis was based on diagnostic criteria recommended by the International Association of Diabetes and Pregnancy Study Groups (IADPSG), with 159 women having a diagnosis of GDM. Cord blood DNA samples from the infants were interrogated for genome-wide DNA methylation levels using the Infinium Human MethylationEPIC BeadChip array. Robust regression was carried out, adjusting for maternal age, smoking, parity, ethnicity, neonate sex, and predicted cell-type composition. Maternal GDM, fasting glucose, 1-h, and 2-h glucose concentrations following an OGTT were associated with 242, 1, 592, and 17 differentially methylated cytosine-phosphate-guanine (dmCpG) sites (false discovery rate (FDR) ≤ 0.05), respectively, in the infant's cord blood DNA. The most significantly GDM-associated CpG was cg03566881 located within the leucine-rich repeat-containing G-protein coupled receptor 6 (LGR6) (FDR = 0.0002). Moreover, we show that the GDM and 1-h glucose-associated methylation signatures in the cord blood of the infant appeared to be attenuated by the dietary and physical activity intervention during pregnancy; in the intervention arm, there were no GDM and two 1-h glucose-associated dmCpGs, whereas in the standard care arm, there were 41 GDM and 160 1-h glucose-associated dmCpGs. A total of 87% of the GDM and 77% of the 1-h glucose-associated dmCpGs had smaller effect sizes in the intervention compared to the standard care arm; the adjusted r2 for the association of LGR6 cg03566881 with GDM was 0.317 (95% confidence interval (CI) 0.012, 0.022) in the standard care and 0.240 (95% CI 0.001, 0.015) in the intervention arm. Limitations included measurement of DNA methylation in cord blood, where the functional significance of such changes are unclear, and because of the strong collinearity between treatment modality and severity of hyperglycaemia, we cannot exclude that treatment-related differences are potential confounders.

CONCLUSIONS

Maternal dysglycaemia was associated with significant changes in the epigenome of the infants. Moreover, we found that the epigenetic impact of a dysglycaemic prenatal maternal environment appeared to be modified by a lifestyle intervention in pregnancy. Further research will be needed to investigate possible medical implications of the findings.

TRIAL REGISTRATION

ISRCTN89971375.

TWIST1 DNA methylation is a cell marker of airway and parenchymal lung fibroblasts that are differentially methylated in asthma

Background

Mesenchymal fibroblasts are ubiquitous cells that maintain the extracellular matrix of organs. Within the lung, airway and parenchymal fibroblasts are crucial for lung development and are altered with disease, but it has been difficult to understand their roles due to the lack of distinct molecular markers. We studied genome-wide DNA methylation and gene expression in airway and parenchymal lung fibroblasts from healthy and asthmatic donors, to identify a robust cell marker and to determine if these cells are molecularly distinct in asthma.

Results

Airway (N = 8) and parenchymal (N = 15) lung fibroblasts from healthy individuals differed in the expression of 158 genes, and DNA methylation of 3936 CpGs (Bonferroni adjusted p value < 0.05). Differential DNA methylation between cell types was associated with differential expression of 42 genes, but no single DNA methylation CpG feature (location, effect size, number) defined the interaction. Replication of gene expression and DNA methylation in a second cohort identified TWIST1 gene expression, DNA methylation and protein expression as a cell marker of airway and parenchymal lung fibroblasts, with DNA methylation having 100% predictive discriminatory power. DNA methylation was differentially altered in parenchymal (112 regions) and airway fibroblasts (17 regions) with asthmatic status, with no overlap between regions.

Conclusions

Differential methylation of TWIST1 is a robust cell marker of airway and parenchymal lung fibroblasts. Airway and parenchymal fibroblast DNA methylation are differentially altered in individuals with asthma, and the role of both cell types should be considered in the pathogenesis of asthma.

The PedBE clock accurately estimates DNA methylation age in pediatric buccal cells

The development of biological markers of aging has primarily focused on adult samples. Epigenetic clocks are a promising tool for measuring biological age that show impressive accuracy across most tissues and age ranges. In adults, deviations from the DNA methylation (DNAm) age prediction are correlated with several age-related phenotypes, such as mortality and frailty. In children, however, fewer such associations have been made, possibly because DNAm changes are more dynamic in pediatric populations as compared to adults. To address this gap, we aimed to develop a highly accurate, noninvasive, biological measure of age specific to pediatric samples using buccal epithelial cell DNAm. We gathered 1,721 genome-wide DNAm profiles from 11 different cohorts of typically developing individuals aged 0 to 20 y old. Elastic net penalized regression was used to select 94 CpG sites from a training dataset (n = 1,032), with performance assessed in a separate test dataset (n = 689). DNAm at these 94 CpG sites was highly predictive of age in the test cohort (median absolute error = 0.35 y). The Pediatric-Buccal-Epigenetic (PedBE) clock was characterized in additional cohorts, showcasing the accuracy in longitudinal data, the performance in nonbuccal tissues and adult age ranges, and the association with obstetric outcomes. The PedBE tool for measuring biological age in children might help in understanding the environmental and contextual factors that shape the DNA methylome during child development, and how it, in turn, might relate to child health and disease.

Methylation of cysteinyl leukotriene receptor 1 genes associates with lung function in asthmatics exposed to traffic-related air pollution

Air pollution is associated with early declines in lung function and increased levels of asthma-related cysteinyl leukotrienes (CysLT) but a biological pathway linking this rapid response has not been delineated. In this randomized controlled diesel exhaust (DE) challenge study of 16 adult asthmatics, increased exposure-attributable urinary leukotriene E4 (uLTE4, a biomarker of cysteinyl leukotriene production) was correlated (p = 0.04) with declines in forced expiratory volume in 1-second (FEV₁) within 6 hours of exposure. Exposure-attributable uLTE4 increases were correlated (p = 0.02) with increased CysLT receptor 1 (CysLTR1) methylation in peripheral blood mononuclear cells which, in turn, was marginally correlated (p = 0.06) with decreased CysLTR1 expression. Decreased CysLTR1 expression was, in turn, correlated (p = 0.0007) with FEV₁ declines. During the same time period, increased methylation of GPR17 (a negative regulator of CysLTR1) was observed after DE exposure (p = 0.02); this methylation increase was correlated (p = 0.001) with decreased CysLTR1 methylation which, in turn, was marginally correlated (p = 0.06) with increased CysLTR1 expression; increased CysLTR1 expression was correlated (p = 0.0007) with FEV₁ increases. Collectively, these data delineate a potential mechanistic pathway linking increased DE exposure-attributable CysLT levels to lung function declines through changes in CysLTR1-related methylation and gene expression.

5-Hydroxymethylcytosine in cord blood and associations of DNA methylation with sex in newborns

DNA methylation has been widely studied for associations with exposures and health outcomes. Both 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) are epigenetic marks that may function differently to impact gene expression; however, the most commonly used technology to assess methylation for population studies in blood use are the Illumina 450K and EPIC BeadChips, for which the traditional bisulfite conversion does not differentiate 5mC and 5hmC marks. We used a modified protocol originally developed by Stewart et al. to analyse oxidative bisulfite-converted and conventional bisulfite-converted DNA for the same subject in parallel by the EPIC chip, allowing us to isolate the two measures. We measured 5mC and 5hmC in cord blood of 41 newborn participants of the Center for Health Assessment of Mothers and Children of Salinas (CHAMACOS) birth cohort and investigated differential methylation of 5mC + 5hmC, isolated 5mC and isolated 5hmC with sex at birth as an example of a biological variable previously associated with DNA methylation. Results showed low levels of 5hmC throughout the epigenome in the cord blood samples in comparison to 5mC. The concordance of autosomal hits between 5mC + 5hmC and exclusive 5mC analyses were low (25%); however, overlap was larger with increased effect size difference. There were 43 autosomal cytosine nucleotide followed by a guanine nucleotide (CpG) sites where 5hmC was associated with sex, 21 of which were unique to 5hmC after adjustment for cell composition. 5hmC only accounts for a small portion of overall methylation in cord blood; however, it has the potential to impact interpretation of combined 5hmC + 5mC studies in cord blood, especially given that effect sizes of differential methylation analyses are often small. Several significant CpG sites were unique to 5hmC, suggesting some functions distinct from 5mC. More studies of genome-wide 5hmC in children are warranted.

Amygdala 5-HTT Gene Network Moderates the Effects of Postnatal Adversity on Attention Problems: Anatomo-Functional Correlation and Epigenetic Changes

Variations in serotoninergic signaling have been related to behavioral outcomes. Alterations in the genome, such as DNA methylation and histone modifications, are affected by serotonin neurotransmission. The amygdala is an important brain region involved in emotional responses and impulsivity, which receives serotoninergic input. In addition, studies suggest that the serotonin transporter gene network may interact with the environment and influence the risk for psychiatric disorders. We propose to investigate whether/how interactions between the exposure to early life adversity and serotonin transporter gene network in the amygdala associate with behavioral disorders. We constructed a co-expression-based polygenic risk score (ePRS) reflecting variations in the function of the serotonin transporter gene network in the amygdala and investigated its interaction with postnatal adversity on attention problems in two independent cohorts from Canada and Singapore. We also described how interactions between ePRS-5-HTT and postnatal adversity exposure predict brain gray matter density and variation in DNA methylation across the genome. We observed that the expression-based polygenic risk score, reflecting the function of the amygdala 5-HTT gene network, interacts with postnatal adversity, to predict attention and hyperactivity problems across both cohorts. Also, both postnatal adversity score and amygdala ePRS-5-HTT score, as well as their interaction, were observed to be associated with variation in DNA methylation across the genome. Variations in gray matter density in brain regions linked to attentional processes were also correlated to our ePRS score. These results confirm that the amygdala 5-HTT gene network is strongly associated with ADHD-related behaviors, brain cortical density, and epigenetic changes in the context of adversity in young children.

Genetic susceptibility to asthma increases the vulnerability to indoor air pollution

INTRODUCTION

Indoor air pollution and maternal smoking during pregnancy are associated with respiratory symptoms in infants, but little is known about the direct association with lung function or interactions with genetic risk factors. We examined associations of exposure to indoor particulate matter with a 50% cut-off aerodynamic diameter of 10 µm (PM₁₀) and maternal smoking with infant lung function and the role of gene-environment interactions.

METHODS

Data from the Drakenstein Child Health Study, a South African birth cohort, were analysed (n=270). Lung function was measured at 6 weeks and 1 year of age, and lower respiratory tract infection episodes were documented. We measured pre- and postnatal PM₁₀ exposures using devices placed in homes, and prenatal tobacco smoke exposure using maternal urine cotinine levels. Genetic risk scores determined from associations with childhood-onset asthma in the UK Biobank were used to investigate effect modifications.

RESULTS

Pre- and postnatal exposure to PM₁₀ as well as maternal smoking during pregnancy were associated with reduced lung function at 6 weeks and 1 year as well as with lower respiratory tract infection in the first year. Due to a significant interaction between the genetic risk score and prenatal exposure to PM₁₀, infants carrying more asthma-related risk alleles were more susceptible to PM₁₀-associated reduced lung function (p_interaction=0.007). This interaction was stronger in infants with Black African ancestry (p_interaction=0.001) and nonexistent in children with mixed ancestry (p_interaction=0.876).

CONCLUSIONS

PM₁₀ and maternal smoking exposures were associated with reduced lung function, with a higher susceptibility for infants with an adverse genetic predisposition for asthma that also depended on the infant's ancestry.

Neonatal Alcohol Exposure in Mice Induces Select Differentiation- and Apoptosis-Related Chromatin Changes Both Independent of and Dependent on Sex

Prenatal alcohol exposure (PAE) affects many aspects of physiology and behavior, including brain development. Specifically, ethanol can influence expression of genes important for brain growth, including chromatin modifiers. Ethanol can also increase apoptotic cell death in the brain and alter epigenetic profiles such as modifications to histones and DNA methylation. Although differential sex outcomes and disruptions to the function of multiple brain regions have been reported in fetal alcohol spectrum disorder (FASD), the majority of our knowledge on molecular epigenetic and apoptotic dysregulation in PAE is based on data from males and is sometimes limited to assessments of the whole brain or one brain region. Here, we examined histone modifications, DNA methylation, and expression of genes involved in differentiation and proliferation related-chromatin modifications and apoptosis in the cerebral cortex and cerebellum of C57BL/6J mice exposed to an acute alcohol challenge on postnatal day 7, with a focus on differential outcomes between sexes and brain regions. We found that neonatal alcohol exposure altered histone modifications, and impacted expression of a select few chromatin modifier and apoptotic genes in both the cortex and cerebellum. The results were observed primarily in a sex-independent manner, although some additional trends toward sexual dimorphisms were observed. Alcohol exposure induced trends toward increased bulk H3K4me3 levels, increased Kmt2e expression, and elevated levels of Casp6 mRNA and bulk γH2A.X. Additional trends indicated that ethanol may impact Kdm4a promoter DNA methylation levels and bulk levels of the histone variant H2A.Z, although further studies are needed. We comprehensively examined effects of ethanol exposure across different sexes and brain regions, and our results suggest that major impacts of ethanol on bulk chromatin modifications underlying differentiation and apoptosis may be broadly applicable across the rodent cortex and cerebellum in both sexes.

CoMeBack: DNA methylation array data analysis for co-methylated regions

Motivation

High-dimensional DNA methylation (DNAm) array coverage, while sparse in the context of the entire DNA methylome, still constitutes a very large number of CpG probes. The ensuing multiple-test corrections affect the statistical power to detect associations, likely contributing to prevalent limited reproducibility. Array probes measuring proximal CpG sites often have correlated levels of DNAm that may not only be biologically meaningful but also imply statistical dependence and redundancy. New methods that account for such correlations between adjacent probes may enable improved specificity, discovery and interpretation of statistical associations in DNAm array data.

Results

We developed a method named Co-Methylation with genomic CpG Background (CoMeBack) that estimates DNA co-methylation, defined as proximal CpG probes with correlated DNAm across individuals. CoMeBack outputs co-methylated regions (CMRs), spanning sets of array probes constructed based on all genomic CpG sites, including those not measured on the array, and without any phenotypic variable inputs. This approach can reduce the multiple-test correction burden, while enhancing the discovery and specificity of statistical associations. We constructed and validated CMRs in whole blood, using publicly available Illumina Infinium 450 K array data from over 5000 individuals. These CMRs were enriched for enhancer chromatin states, and binding site motifs for several transcription factors involved in blood physiology. We illustrated how CMR-based epigenome-wide association studies can improve discovery and reduce false positives for associations with chronological age.

Availability and implementation

https://bitbucket.org/flopflip/comeback.

Supplementary information

Supplementary data are available at Bioinformatics online.

A Role of Oxytocin Receptor Gene Brain Tissue Expression Quantitative Trait Locus rs237895 in the Intergenerational Transmission of the Effects of Maternal Childhood Maltreatment

OBJECTIVE

Women exposed to childhood maltreatment (CM) are more likely to exhibit insensitive parenting, which may have consequences for their offspring's development. Variation in the oxytocin-receptor gene (OXTR) moderates risk of CM-associated long-term sequelae associated with mother-child attachment, although functionality of previously investigated single nucleotide polymorphisms (SNPs) remained elusive. Here, we investigated the role of OXTR rs237895, a brain tissue expression quantitative trait locus (eQTL), as a moderator of the relationship between CM and maternal behavior (MB) and the association between MB and offspring attachment security.

METHOD

Of 110 women with information on rs237895 genotype (T-allele = 64, CC = 46), 107 had information on CM (CTQ) and 99 on standardized observer-based ratings of MB at 6 months postpartum (responsivity and detachment), which were used in principal component analysis to obtain a latent factor representing MB. Offspring (n = 86) attachment was evaluated at 12 months of age. Analyses predicting MB were adjusted for socioeconomic status, age, postpartum depression, and genotype-based ethnicity. Analyses predicting child attachment were adjusted for infant sex, socioeconomic status, and postpartum depression.

RESULTS

rs237895 significantly moderated the relationship between CM and MB (F1;66 = 7.99, p < .01), indicating that CM was associated with maternal insensitivity only in high-OXTR-expressing T-allele carriers but not in low-OXTR-expressing CC homozygotes. Moreover, maternal insensitivity predicted offspring insecure attachment (B = -0.551; p < .05).

CONCLUSION

Women with a high OXTR expressing genotype are more susceptible to CM-related impairments in MB that, in turn, predict attachment security in their children, supporting the role of the OT system in the intergenerational transmission of risk associated with maternal CM.

Reversal of epigenetic aging and immunosenescent trends in humans

Epigenetic “clocks” can now surpass chronological age in accuracy for estimating biological age. Here, we use four such age estimators to show that epigenetic aging can be reversed in humans. Using a protocol intended to regenerate the thymus, we observed protective immunological changes, improved risk indices for many age‐related diseases, and a mean epigenetic age approximately 1.5 years less than baseline after 1 year of treatment (−2.5‐year change compared to no treatment at the end of the study). The rate of epigenetic aging reversal relative to chronological age accelerated from −1.6 year/year from 0–9 month to −6.5 year/year from 9–12 month. The GrimAge predictor of human morbidity and mortality showed a 2‐year decrease in epigenetic vs. chronological age that persisted six months after discontinuing treatment. This is to our knowledge the first report of an increase, based on an epigenetic age estimator, in predicted human lifespan by means of a currently accessible aging intervention.

Epigenetic aging in children from a small-scale farming society in The Congo Basin: Associations with child growth and family conflict

Developmental environments influence individuals' long-term health trajectories, and there is increasing emphasis on understanding the biological pathways through which this occurs. Epigenetic aging evaluates DNA methylation at a suite of distinct CpG sites in the genome, and epigenetic age acceleration (EAA) is linked to heightened chronic morbidity and mortality risks in adults. Consequently, EAA provides insights on trajectories of biological aging, which early life experiences may help shape. However, few studies have measured correlates of children's epigenetic aging, especially outside of the U.S. and Europe. In particular, little is known about how children's growth and development relate to EAA in ecologies in which energetic and pathogenic stressors are commonplace. We studied EAA from dried blood spots among Bondongo children (n = 54) residing in a small-scale, fisher-farmer society in a remote region of the Republic of the Congo. Here, infectious disease burdens and their resultant energy demands are high. Children who were heavier for height or taller for age, respectively, exhibited greater EAA, including intrinsic EAA, which is considered to measure EAA internal to cells. Furthermore, we found that children in families with more conflict between parents had greater intrinsic EAA. These results suggest that in contexts in which limited energy must be allocated to competing demands, more investment in growth may coincide with greater EAA, which parallels findings in European children who do not face similar energetic constraints. Our findings also indicate that associations between adverse family environments and greater intrinsic EAA were nonetheless observable but only after adjustment for covariates relevant to the energetically and immunologically demanding nature of the local ecology.

Airway epithelial cell isolation techniques affect DNA methylation profiles with consequences for analysis of asthma related perturbations to DNA methylation

The airway epithelium forms the interface between the inhaled environment and the lung. The airway epithelium is dysfunctional in asthma and epigenetic mechanisms are considered a contributory factor. We hypothesised that the DNA methylation profiles of cultured primary airway epithelial cells (AECs) would differ between cells isolated from individuals with asthma (n = 17) versus those without asthma (n = 16). AECs were isolated from patients by two different isolation techniques; pronase digestion (9 non-asthmatic, 8 asthmatic) and bronchial brushings (7 non-asthmatic and 9 asthmatic). DNA methylation was assessed using an Illumina Infinium HumanMethylation450 BeadChip array. DNA methylation of AECs clustered by isolation technique and linear regression identified 111 CpG sites differentially methylated between isolation techniques in healthy individuals. As a consequence, the effect of asthmatic status on DNA methylation was assessed within AEC samples isolated using the same technique. In pronase isolated AECs, 15 DNA regions were differentially methylated between asthmatics and non-asthmatics. In bronchial brush isolated AECs, 849 differentially methylated DNA regions were identified with no overlap to pronase regions. In conclusion, regardless of cell isolation technique, differential DNA methylation was associated with asthmatic status in AECs, providing further evidence for aberrant DNA methylation as a signature of epithelial dysfunction in asthma.

DNA methylation of HPA-axis genes and the onset of major depressive disorder in adolescent girls: a prospective analysis

The stress response system is disrupted in individuals with major depressive disorder (MDD) as well as in those at elevated risk for developing MDD. We examined whether DNA methylation (DNAm) levels of CpG sites within HPA-axis genes predict the onset of MDD. Seventy-seven girls, approximately half (n = 37) of whom were at familial risk for MDD, were followed longitudinally. Saliva samples were taken in adolescence (M age = 13.06 years [SD = 1.52]) when participants had no current or past MDD diagnosis. Diagnostic interviews were administered approximately every 18 months until the first onset of MDD or early adulthood (M age of last follow-up = 19.23 years [SD = 2.69]). We quantified DNAm in saliva samples using the Illumina EPIC chip and examined CpG sites within six key HPA-axis genes (NR3C1, NR3C2, CRH, CRHR1, CRHR2, FKBP5) alongside 59 genotypes for tagging SNPs capturing cis genetic variability. DNAm levels within CpG sites in NR3C1, CRH, CRHR1, and CRHR2 were associated with risk for MDD across adolescence and young adulthood. To rule out the possibility that findings were merely due to the contribution of genetic variability, we re-analyzed the data controlling for cis genetic variation within these candidate genes. Importantly, methylation levels in these CpG sites continued to significantly predict the onset of MDD, suggesting that variation in the epigenome, independent of proximal genetic variants, prospectively predicts the onset of MDD. These findings suggest that variation in the HPA axis at the level of the methylome may predict the development of MDD.

MRE11-RAD50-NBS1 promotes Fanconi Anemia R-loop suppression at transcription-replication conflicts

Ectopic R-loop accumulation causes DNA replication stress and genome instability. To avoid these outcomes, cells possess a range of anti-R-loop mechanisms, including RNaseH that degrades the RNA moiety in R-loops. To comprehensively identify anti-R-loop mechanisms, we performed a genome-wide trigenic interaction screen in yeast lacking RNH1 and RNH201. We identified >100 genes critical for fitness in the absence of RNaseH, which were enriched for DNA replication fork maintenance factors including the MRE11-RAD50-NBS1 (MRN) complex. While MRN has been shown to promote R-loops at DNA double-strand breaks, we show that it suppresses R-loops and associated DNA damage at transcription-replication conflicts. This occurs through a non-nucleolytic function of MRE11 that is important for R-loop suppression by the Fanconi Anemia pathway. This work establishes a novel role for MRE11-RAD50-NBS1 in directing tolerance mechanisms at transcription-replication conflicts.

Regulation of Skn7-dependent, oxidative stress-induced genes by the RNA polymerase II-CTD phosphatase, Fcp1, and Mediator kinase subunit, Cdk8, in yeast

Fcp1 is a protein phosphatase that facilitates transcription elongation and termination by dephosphorylating the C-terminal domain of RNA polymerase II. High-throughput genetic screening and gene expression profiling of fcp1 mutants revealed a novel connection to Cdk8, the Mediator complex kinase subunit, and Skn7, a key transcription factor in the oxidative stress response pathway. Briefly, Skn7 was enriched as a regulator of genes whose mRNA levels were altered in fcp1 and cdk8Δ mutants and was required for the suppression of fcp1 mutant growth defects by loss of CDK8 under oxidative stress conditions. Targeted analysis revealed that mutating FCP1 decreased Skn7 mRNA and protein levels as well as its association with target gene promoters but paradoxically increased the mRNA levels of Skn7-dependent oxidative stress-induced genes (TRX2 and TSA1) under basal and induced conditions. The latter was in part recapitulated via chemical inhibition of transcription in WT cells, suggesting that a combination of transcriptional and posttranscriptional effects underscored the increased mRNA levels of TRX2 and TSA1 observed in the fcp1 mutant. Interestingly, loss of CDK8 robustly normalized the mRNA levels of Skn7-dependent genes in the fcp1 mutant background and also increased Skn7 protein levels by preventing its turnover. As such, our work suggested that loss of CDK8 could overcome transcriptional and/or posttranscriptional alterations in the fcp1 mutant through its regulatory effect on Skn7. Furthermore, our work also implicated FCP1 and CDK8 in the broader response to environmental stressors in yeast.

Systematic evaluation and validation of reference and library selection methods for deconvolution of cord blood DNA methylation data

BACKGROUND

Umbilical cord blood (UCB) is commonly used in epigenome-wide association studies of prenatal exposures. Accounting for cell type composition is critical in such studies as it reduces confounding due to the cell specificity of DNA methylation (DNAm). In the absence of cell sorting information, statistical methods can be applied to deconvolve heterogeneous cell mixtures. Among these methods, reference-based approaches leverage age-appropriate cell-specific DNAm profiles to estimate cellular composition. In UCB, four reference datasets comprising DNAm signatures profiled in purified cell populations have been published using the Illumina 450 K and EPIC arrays. These datasets are biologically and technically different, and currently, there is no consensus on how to best apply them. Here, we systematically evaluate and compare these datasets and provide recommendations for reference-based UCB deconvolution.

RESULTS

We first evaluated the four reference datasets to ascertain both the purity of the samples and the potential cell cross-contamination. We filtered samples and combined datasets to obtain a joint UCB reference. We selected deconvolution libraries using two different approaches: automatic selection using the top differentially methylated probes from the function pickCompProbes in minfi and a standardized library selected using the IDOL (Identifying Optimal Libraries) iterative algorithm. We compared the performance of each reference separately and in combination, using the two approaches for reference library selection, and validated the results in an independent cohort (Generation R Study, n = 191) with matched Fluorescence-Activated Cell Sorting measured cell counts. Strict filtering and combination of the references significantly improved the accuracy and efficiency of cell type estimates. Ultimately, the IDOL library outperformed the library from the automatic selection method implemented in pickCompProbes.

CONCLUSION

These results have important implications for epigenetic studies in UCB as implementing this method will optimally reduce confounding due to cellular heterogeneity. This work provides guidelines for future reference-based UCB deconvolution and establishes a framework for combining reference datasets in other tissues.

Gestational intermittent hypoxia induces endothelial dysfunction, reduces perivascular adiponectin and causes epigenetic changes in adult male offspring

KEY POINTS

Obstructive sleep apnoea (OSA) is characterized by intermittent hypoxia, which causes oxidative stress and inflammation and increases the risk of cardiovascular disease. OSA during pregnancy causes adverse maternal and fetal outcomes. The effects of pre-existing OSA in pregnant women on cardiometabolic outcomes in the offspring are unknown. We evaluated basic metabolic parameters, as well as aortic vascular and perivascular adipose tissue (PVAT) function in response to adiponectin, and examined DNA methylation of adiponectin gene promoter in PVAT in 16-week-old adult offspring exposed to gestational intermittent hypoxia (GIH). GIH decreased body weights at week 1 in both male and female offspring, and caused subsequent increases in body weight and food consumption in male offspring only. Adult female offspring had normal levels of lipids, glucose and insulin, with no endothelial dysfunction. Adult male offspring exhibited dyslipidaemia, insulin resistance and hyperleptinaemia. Decreased endothelial-dependent vasodilatation, loss of anti-contractile activity of PVAT and low circulating PVAT adiponectin levels, as well as increased pro-inflammatory gene expression and DNA methylation of adiponectin gene promoter, occurred in adult male offspring. Our results suggest that male offspring of women with OSA could be at risk of developing cardiometabolic disease during adulthood.

ABSTRACT

Perturbations during pregnancy can program the offspring to develop cardiometabolic diseases later in life. Obstructive sleep apnoea (OSA) is a chronic condition that frequently affects pregnancies and leads to adverse fetal outcomes. We assessed the offspring of female mice experiencing gestational intermittent hypoxia (GIH), a hallmark of OSA, for changes in metabolic profiles, aortic nitric oxide (NO)-dependent relaxations, perivascular adipose tissue (PVAT) anti-contractile activities and the responses to adiponectin, and DNA methylation of the adiponectin gene promoter in PVAT tissue. Pregnant mouse dams were exposed to intermittent hypoxic cycles ( F I O 2 21-12%) for 18 days. GIH resulted in lower body weights of pups at week 1, followed by significant weight gain by week 16 of age in male but not female offspring. Plasma lipids, leptin and insulin resistance were higher in GIH male adult offspring. Endothelium-dependent relaxation in response to ACh and the anti-contractile activity of PVAT in the abdominal aorta was reduced in GIH adult male offspring. Incubation of arteries from GIH adult male offspring with adiponectin restored the anti-contractile activity of PVAT. Both circulating and PVAT tissue homogenate levels of adiponectin, as well as gene expression of adiponectin in PVAT, were lower in GIH male offspring, along with an increased gene expression of inflammatory cytokines. Pyrosequencing of adiponectin gene promoter in PVAT showed increased DNA methylation in GIH male offspring. Our results indicate that GIH leads to vascular disease in adult male offspring through PVAT dysfunction, which was associated with low adiponectin levels and epigenetic modifications on the adiponectin gene promoter.

Inflammation-related genes are associated with epigenetic aging in HIV

Chronic inflammation is characteristic of both HIV and aging ("inflammaging") and may contribute to the accelerated aging observed in people living with HIV (PLWH). We examined whether three inflammation-related single-nucleotide polymorphisms (SNPs) were risk factors for accelerated aging and HIV-associated, non-AIDS (HANA) conditions among PLWH. We examined 155 postmortem cases with HIV (mean age = 47.3, 81% male, 68% self-reported White) from the National NeuroAIDS Tissue Consortium who had pre-mortem neurobehavioral/medical/virologic data and epigenomic data from occipital cortex tissue. Accelerated aging was measured according to the Epigenetic Clock; an aging biomarker based on DNA methylation levels. Past or current age-associated HANA conditions including cerebrovascular, liver and kidney disease, chronic obstructive pulmonary disease, cancer, and diabetes were determined via self-report. Epigenetic Aging Z-scores and likelihood of past/current HANA conditions were compared between major allele homozygotes and minor allele carriers for each SNP (IL-6 - 174G>C, IL-10 - 592C>A, TNF-α - 308 G>A) separately. Analyses were adjusted for relevant demographic/clinical factors. Epigenetic aging (e.g., higher Z-scores) was significantly greater in IL-6 C allele carriers (p = .002) and IL-10 CC homozygotes (p = .02) compared to other genotype groups. The likelihood of any past/current HANA condition did not differ by IL-10 genotype but was 3.36 times greater in IL-6 C allele carriers versus others (OR = 3.36, 95%CI = 1.09-10.34, p = .03). TNF-α genotype was not associated with epigenetic aging or HANA conditions. IL-6 and IL-10 SNPs may help to identify PLWH who are at high risk for accelerated aging. These insights into pathophysiological pathways may inform interventional approaches to treat rapid aging among PLWH.

Neonatal amygdalae and hippocampi are influenced by genotype and prenatal environment, and reflected in the neonatal DNA methylome

The amygdala and hippocampus undergo rapid development in early life. The relative contribution of genetic and environmental factors to the establishment of their developmental trajectories has yet to be examined. We performed imaging on neonates and examined how the observed variation in volume and microstructure of the amygdala and hippocampus varied by genotype, and compared with prenatal maternal mental health and socioeconomic status. Gene × Environment models outcompeted models containing genotype or environment only to best explain the majority of measures but some, especially of the amygdaloid microstructure, were best explained by genotype only. Models including DNA methylation measured in the neonate umbilical cords outcompeted the Gene and Gene × Environment models for the majority of amygdaloid measures and minority of hippocampal measures. This study identified brain region-specific gene networks associated with individual differences in fetal brain development. In particular, genetic and epigenetic variation within CUX1 was highlighted.

Dynamic DNA methylation changes in the maternal oxytocin gene locus (OXT) during pregnancy predict postpartum maternal intrusiveness

Maternal behavior (MB) is observable across mammals and represents an important feature of environmental variation during early postnatal development. Oxytocin (OT) plays a crucial role in MB. Even prior to childbirth, pregnancy induces epigenetic and other downstream changes in the maternal OT-system, likely mediated by the actions of steroid hormones. However, little is known about the nature and consequences of epigenetic modifications in the maternal OT-encoding gene (OXT) during pregnancy. Our study aims to investigate temporal dynamics of OXT promoter DNA methylation (DNAm) throughout pregnancy in predicting MB in humans. In 107 mother-child dyads, maternal OXT DNAm was serially analyzed in whole blood in early, mid and late pregnancy. MB was coded based on standardized mother-child interactions at six months postpartum. After controlling for cellular heterogeneity, race/ethnicity, age, and socioeconomic status, OXT-promoter DNAm exhibited a dynamic profile during pregnancy (b = 0.026, t=-3.37, p < .001), with decreases in DNAm from early to mid-pregnancy and no further change until late pregnancy. Moreover, dynamic DNAm trajectories of the OXT-promoter region predicted MB (intrusiveness) at six months postpartum (b = 0.006, t = 2.0, p < 0.05), with 6% higher OXT DNAm in late pregnancy in intrusive compared to non-intrusive mothers. We here demonstrate that OXT promoter DNAm changes significantly throughout gestation in peripheral blood and that these changes are associated with variability in MB, providing a novel potential biomarker predicting postnatal MB.

Glutaminase Deficiency Caused by Short Tandem Repeat Expansion in GLS

We report an inborn error of metabolism caused by an expansion of a GCA-repeat tract in the 5' untranslated region of the gene encoding glutaminase (GLS) that was identified through detailed clinical and biochemical phenotyping, combined with whole-genome sequencing. The expansion was observed in three unrelated patients who presented with an early-onset delay in overall development, progressive ataxia, and elevated levels of glutamine. In addition to ataxia, one patient also showed cerebellar atrophy. The expansion was associated with a relative deficiency of GLS messenger RNA transcribed from the expanded allele, which probably resulted from repeat-mediated chromatin changes upstream of the GLS repeat. Our discovery underscores the importance of careful examination of regions of the genome that are typically excluded from or poorly captured by exome sequencing.

A biologically-informed polygenic score identifies endophenotypes and clinical conditions associated with the insulin receptor function on specific brain regions

Background

Activation of brain insulin receptors modulates reward sensitivity, inhibitory control and memory. Variations in the functioning of this mechanism likely associate with individual differences in the risk for related mental disorders (attention deficit hyperactivity disorder or ADHD, addiction, dementia), in agreement with the high co-morbidity between insulin resistance and psychopathology. These neurobiological mechanisms can be explored using genetic studies. We propose a novel, biologically informed genetic score reflecting the mesocorticolimbic and hippocampal insulin receptor-related gene networks, and investigate if it predicts endophenotypes (impulsivity, cognitive ability) in community samples of children, and psychopathology (addiction, dementia) in adults.

Methods

Lists of genes co-expressed with the insulin receptor in the mesocorticolimbic system or hippocampus were created. SNPs from these genes (post-clumping) were compiled in a polygenic score using the association betas described in a conventional GWAS (ADHD in the mesocorticolimbic score and Alzheimer in the hippocampal score). Across multiple samples (n = 4502), the biologically informed, mesocorticolimbic or hippocampal specific insulin receptor polygenic scores were calculated, and their ability to predict impulsivity, risk for addiction, cognitive performance and presence of Alzheimer's disease was investigated.

Findings

The biologically-informed ePRS-IR score showed better prediction of child impulsivity and cognitive performance, as well as risk for addiction and Alzheimer's disease in comparison to conventional polygenic scores for ADHD, addiction and dementia.

Interpretation

This novel, biologically-informed approach enables the use of genomic datasets to probe relevant biological processes involved in neural function and disorders.

Fund

Toxic Stress Research network of the JPB Foundation, Jacobs Foundation (Switzerland), Sackler Foundation.

Genome-wide analysis of DNA methylation in relation to socioeconomic status during development and early adulthood

OBJECTIVES

Socioeconomic status (SES) is a powerful determinant of health, but the underlying biological mechanisms are poorly understood. This study investigates whether levels of DNA methylation at CpG sites across the genome are associated with SES in a cohort of young adults in the Philippines.

METHODS

DNA methylation was assayed with the Illumina HumanMethylation450 Bead Chip, in leukocytes from 489 participants in the Cebu Longitudinal Health and Nutrition Survey (mean age = 20.9 years). SES was measured in infancy/childhood and adulthood, and was based on composite measures of income, assets, and education. Genome-wide analysis of variable probes identified CpG sites significantly associated with SES after adjustment for multiple comparisons. Functional enrichment analysis was used to identify biological pathways associated with these sites.

RESULTS

A total of 2,546 CpG sites, across 1,537 annotated genes, were differentially methylated in association with SES. In comparison with high SES, low SES was associated with increased methylation at 1,777 sites, and decreased methylation at 769 sites. Functional enrichment analysis identified over-representation of biological pathways related to immune function, skeletal development, and development of the nervous system.

CONCLUSIONS

Socioeconomic status predicts DNA methylation at a large number of CpG sites across the genome. The scope of these associations is commensurate with the wide range of biological systems and health outcomes that are shaped by SES, and these findings suggest that DNA methylation may play an important role.

Analysis of two birth tissues provides new insights into the epigenetic landscape of neonates born preterm

Background

Preterm birth (PTB), defined as child birth before completion of 37 weeks of gestation, is a major challenge in perinatal health care and can bear long-term medical and financial burden. Over a million children die each year due to PTB complications, and those who survive can face developmental delays. Unfortunately, our understanding of the molecular pathways associated with PTB remains limited. There is a growing body of evidence suggesting the role of DNA methylation (DNAm) in mediating the effects of PTB on future health outcomes. Thus, epigenome-wide association studies (EWAS), where DNAm sites are examined for associations with PTB, can help shed light on the biological mechanisms linking the two.

Results

In an Asian cohort of 1019 infants (68 preterm, 951 full term), we examined and compared the associations between PTB and genome-wide DNAm profiles using both cord tissue (n = 1019) and cord blood (n = 332) samples on Infinium HumanMethylation450 arrays. PTB was significantly associated (P < 5.8e−7) with DNAm at 296 CpGs (209 genes) in the cord blood. Over 95% of these CpGs were replicated in other PTB/gestational age EWAS conducted in (cord) blood. This replication was apparent even across populations of different ethnic origin (Asians, Caucasians, and African Americans). More than a third of these 296 CpGs were replicated in at least 4 independent studies, thereby identifying a robust set of PTB-linked epigenetic signatures in cord blood. Interrogation of cord tissue in addition to cord blood provided novel insights into the epigenetic status of the neonates born preterm. Overall, 994 CpGs (608 genes, P < 3.7e−7) associated with PTB in cord tissue, of which only 10 of these CpGs were identified in the analysis using cord blood. Genes from cord tissue showed enrichment of molecular pathways related to fetal growth and development, while those from cord blood showed enrichment of immune response pathways. A substantial number of PTB-associated CpGs from both the birth tissues were also associated with gestational age.

Conclusions

Our findings provide insights into the epigenetic landscape of neonates born preterm, and that its status is captured more comprehensively by interrogation of more than one neonatal tissue in tandem. Both these neonatal tissues are clinically relevant in their unique ways and require careful consideration in identification of biomarkers related to PTB and gestational age.

Trial registration

This birth cohort is a prospective observational study designed to study the developmental origins of health and disease, and was retrospectively registered on 1 July 2010 under the identifier NCT01174875.

Electronic supplementary material

The online version of this article (10.1186/s13148-018-0599-4) contains supplementary material, which is available to authorized users.