"Gap hunting" to characterize clustered probe signals in Illumina methylation array data

Methylation patterns associated with C-reactive protein in racially and ethnically diverse populations

Systemic low-grade inflammation is a feature of chronic disease. C-reactive protein (CRP) is a common biomarker of inflammation and used as an indicator of disease risk; however, the role of inflammation in disease is not completely understood. Methylation is an epigenetic modification in the DNA which plays a pivotal role in gene expression. In this study we evaluated differential DNA methylation patterns associated with blood CRP level to elucidate biological pathways and genetic regulatory mechanisms to improve the understanding of chronic inflammation. The racially and ethnically diverse participants in this study were included as 50% White, 41% Black or African American, 7% Hispanic or Latino/a, and 2% Native Hawaiian, Asian American, American Indian, or Alaska Native (total n = 13,433) individuals. We replicated 113 CpG sites from 87 unique loci, of which five were novel (CADM3, NALCN, NLRC5, ZNF792, and cg03282312), across a discovery set of 1,150 CpG sites associated with CRP level (p < 1.2E-7). The downstream pathways affected by DNA methylation included the identification of IFI16 and IRF7 CpG-gene transcript pairs which contributed to the innate immune response gene enrichment pathway along with NLRC5, NOD2, and AIM2. Gene enrichment analysis also identified the nuclear factor-kappaB transcription pathway. Using two-sample Mendelian randomization (MR) we inferred methylation at three CpG sites as causal for CRP levels using both White and Black or African American MR instrument variables. Overall, we identified novel CpG sites and gene transcripts that could be valuable in understanding the specific cellular processes and pathogenic mechanisms involved in inflammation.

The Pace of Biological Aging Predicts Nonspecific Chronic Low Back Pain Severity

This study aimed to determine if and how the pace of biological aging was associated with nonspecific chronic low back pain (cLBP) and compare what measure of epigenetic age acceleration most strongly predicts cLBP outcomes. We used the Dunedin Pace of Aging from the Epigenome (DunedinPACE), Horvath's, Hannum's, and PhenoAge clocks to determine the pace of biological aging in 69 cLBP, and 49 pain-free controls (PFCs) adults, ages 18 to 85 years. On average, participants with cLBP had higher DunedinPACE (P < .001) but lower Horvath (P = .04) and Hannum (P = .02) accelerated epigenetic age than PFCs. There was no significant difference in PhenoAge acceleration between the cLBP and PFC groups (P = .97). DunedinPACE had the largest effect size (Cohen's d = .78) on group differences. In univariate regressions, a unit increase in DunedinPACE score was associated with 265.98 times higher odds of cLBP than the PFC group (P < .001). After controlling for sex, race, and body mass index (BMI), the odds ratio of cLBP to PFC group was 149.62 (P < .001). Furthermore, among participants with cLBP, DunedinPACE scores positively correlated with pain severity (rs = .385, P = .001) and interference (rs = .338, P = .005). Epigenetic age acceleration from Horvath, Hannum, and PhenoAge clocks were not significant predictors of cLBP. The odds of a faster pace of biological aging are higher among adults with cLBP, and this was associated with greater pain severity and disability. Future interventions to slow the pace of biological aging may improve cLBP outcomes. PERSPECTIVE: Accelerated epigenetic aging is common among adults with nonspecific cLBP. Higher DunedinPACE scores positively correlate with pain severity and interference, and better predict cLBP than other DNA methylation clocks. Interventions to slow the pace of biological aging may be viable targets for improving pain outcomes.

Faecal Microbiota transplantation affects liver DNA methylation in Non-alcoholic fatty liver disease: a multi-omics approach

Individuals with nonalcoholic fatty liver disease (NAFLD) have an altered gut microbiota composition. Moreover, hepatic DNA methylation may be altered in the state of NAFLD. Using a fecal microbiota transplantation (FMT) intervention, we aimed to investigate whether a change in gut microbiota composition relates to altered liver DNA methylation in NAFLD. Moreover, we assessed whether plasma metabolite profiles altered by FMT relate to changes in liver DNA methylation. Twenty-one individuals with NAFLD underwent three 8-weekly vegan allogenic donor (n = 10) or autologous (n = 11) FMTs. We obtained hepatic DNA methylation profiles from paired liver biopsies of study participants before and after FMTs. We applied a multi-omics machine learning approach to identify changes in the gut microbiome, peripheral blood metabolome and liver DNA methylome, and analyzed cross-omics correlations. Vegan allogenic donor FMT compared to autologous FMT induced distinct differential changes in I) gut microbiota profiles, including increased abundance of Eubacterium siraeum and potential probiotic Blautia wexlerae; II) plasma metabolites, including altered levels of phenylacetylcarnitine (PAC) and phenylacetylglutamine (PAG) both from gut-derived phenylacetic acid, and of several choline-derived long-chain acylcholines; and III) hepatic DNA methylation profiles, most importantly in Threonyl-TRNA Synthetase 1 (TARS) and Zinc finger protein 57 (ZFP57). Multi-omics analysis showed that Gemmiger formicillis and Firmicutes bacterium_CAG_170 positively correlated with both PAC and PAG. E siraeum negatively correlated with DNA methylation of cg16885113 in ZFP57. Alterations in gut microbiota composition by FMT caused widespread changes in plasma metabolites (e.g. PAC, PAG, and choline-derived metabolites) and liver DNA methylation profiles in individuals with NAFLD. These results indicate that FMTs might induce metaorganismal pathway changes, from the gut bacteria to the liver.

Novel DNA methylome biomarkers associated with adalimumab response in rheumatoid arthritis patients

Background and aims

Rheumatoid arthritis (RA) patients are currently treated with biological agents mostly aimed at cytokine blockade, such as tumor necrosis factor-alpha (TNFα). Currently, there are no biomarkers to predict therapy response to these agents. Here, we aimed to predict response to adalimumab (ADA) treatment in RA patients using DNA methylation in peripheral blood (PBL).

Methods

DNA methylation profiling on whole peripheral blood from 92 RA patients before the start of ADA treatment was determined using Illumina HumanMethylationEPIC BeadChip array. After 6 months, treatment response was assessed according to the European Alliance of Associations for Rheumatology (EULAR) criteria for disease activity. Patients were classified as responders (Disease Activity Score in 28 Joints (DAS28) < 3.2 or decrease of 1.2 points) or as non-responders (DAS28 > 5.1 or decrease of less than 0.6 points). Machine learning models were built through stability-selected gradient boosting to predict response prior to ADA treatment with predictor DNA methylation markers.

Results

Of the 94 RA patients, we classified 49 and 43 patients as responders and non-responders, respectively. We were capable of differentiating responders from non-responders with a high performance (area under the curve (AUC) 0.76) using a panel of 27 CpGs. These classifier CpGs are annotated to genes involved in immunological and pathophysiological pathways related to RA such as T-cell signaling, B-cell pathology, and angiogenesis.

Conclusion

Our findings indicate that the DNA methylome of PBL provides discriminative capabilities in discerning responders and non-responders to ADA treatment and may therefore serve as a tool for therapy prediction.

Gaps in current methods to detect polymorphic CpGs from Illumina Infinium human methylation microarrays and exploring their potential impact in multi-EWAS analyses

DNA methylation (DNAm) epigenome-wide association studies (EWAS) have been performed on diverse ethnicities to discover novel biomarkers associated with various diseases, such as cancers, autoimmune diseases, and neurological disorders. However, genetic polymorphisms can influence DNAm levels resulting in methylation quantitative trait loci (meQTL). These can be either direct effects, by altering the sequence of the methylation (CpG) site itself, or, in the case of array-based measures, indirectly altering the detection probe-binding site interaction. Given that genetic variant frequencies associated with meQTL can differ between population groups, these have the potential to confound EWAS observations, particularly in multi-ethnic populations. In this study, we analysed publicly available DNA methylation profiles (450K array), consisting of 1342 individuals from 6 distinct ancestral groups. We investigate two distinct tools (GapHunter and MethylToSNP) specifically designed to identify CpG sites that may be influenced by genetic variation. Results from this aggregated trans-ancestral epigenome-wide dataset suggest that both tools fail to consistently identify not only rarer (MAF < 0.05) genetic variant effects but also more than half of sites predicted to be associated with variants with much higher allele frequencies (MAF >0.2). In addition, there is a relatively low concordance in the detection of polymorphic CpGs between GapHunter and MethylToSNP. Screening of CpG site associations from EWAS using either of these tools is unlikely to be a robust or comprehensive means of identifying all genetic variant confounding effects.

Multi-ancestry epigenome-wide analyses identify methylated sites associated with aortic augmentation index in TOPMed MESA

Despite the prognostic value of arterial stiffness (AS) and pulsatile hemodynamics (PH) for cardiovascular morbidity and mortality, epigenetic modifications that contribute to AS/PH remain unknown. To gain a better understanding of the link between epigenetics (DNA methylation) and AS/PH, we examined the relationship of eight measures of AS/PH with CpG sites and co-methylated regions using multi-ancestry participants from Trans-Omics for Precision Medicine (TOPMed) Multi-Ethnic Study of Atherosclerosis (MESA) with sample sizes ranging from 438 to 874. Epigenome-wide association analysis identified one genome-wide significant CpG (cg20711926-CYP1B1) associated with aortic augmentation index (AIx). Follow-up analyses, including gene set enrichment analysis, expression quantitative trait methylation analysis, and functional enrichment analysis on differentially methylated positions and regions, further prioritized three CpGs and their annotated genes (cg23800023-ETS1, cg08426368-TGFB3, and cg17350632-HLA-DPB1) for AIx. Among these, ETS1 and TGFB3 have been previously prioritized as candidate genes. Furthermore, both ETS1 and HLA-DPB1 have significant tissue correlations between Whole Blood and Aorta in GTEx, which suggests ETS1 and HLA-DPB1 could be potential biomarkers in understanding pathophysiology of AS/PH. Overall, our findings support the possible role of epigenetic regulation via DNA methylation of specific genes associated with AIx as well as identifying potential targets for regulation of AS/PH.

Associations between epigenome-wide DNA methylation and height-related traits among Sub-Saharan Africans: the RODAM study

Human height and related traits are highly complex, and extensively research has shown that these traits are determined by both genetic and environmental factors. Such factors may partially affect these traits through epigenetic programing. Epigenetic programing is dynamic and plays an important role in controlling gene expression and cell differentiation during (early) development. DNA methylation (DNAm) is the most commonly studied epigenetic feature. In this study we conducted an epigenome-wide DNAm association analysis on height-related traits in a Sub-Saharan African population, in order to detect DNAm biomarkers across four height-related traits. DNAm profiles were acquired in whole blood samples of 704 Ghanaians, sourced from the Research on Obesity and Diabetes among African Migrants study, using the Illumina Infinium HumanMethylation450 BeadChip. Linear models were fitted to detect differentially methylated positions (DMPs) and regions (DMRs) associated with height, leg-to-height ratio (LHR), leg length, and sitting height. No epigenome-wide significant DMPs were recorded. However we did observe among our top DMPs five informative probes associated with the height-related traits: cg26905768 (leg length), cg13268132 (leg length), cg19776793 (height), cg23072383 (LHR), and cg24625894 (sitting height). All five DMPs are annotated to genes whose functions were linked to bone cell regulation and development. DMR analysis identified overlapping DMRs within the gene body of HLA-DPB1 gene, and the HOXA gene cluster. In this first epigenome-wide association studies of these traits, our findings suggest DNAm associations with height-related heights, and might influence development and maintenance of these traits. Further studies are needed to replicate our findings, and to elucidate the molecular mechanism underlying human height-related traits.

Prenatal exposures to per- and polyfluoroalkyl substances and epigenetic aging in umbilical cord blood: The Healthy Start study

BACKGROUND

Per- and polyfluoroalkyl substances (PFAS) are ubiquitous, environmentally persistent chemicals, and prenatal exposures have been associated with adverse child health outcomes. Prenatal PFAS exposure may lead to epigenetic age acceleration (EAA), defined as the discrepancy between an individual's chronologic and epigenetic or biological age.

OBJECTIVES

We estimated associations of maternal serum PFAS concentrations with EAA in umbilical cord blood DNA methylation using linear regression, and a multivariable exposure-response function of the PFAS mixture using Bayesian kernel machine regression.

METHODS

Five PFAS were quantified in maternal serum (median: 27 weeks of gestation) among 577 mother-infant dyads from a prospective cohort. Cord blood DNA methylation data were assessed with the Illumina HumanMethylation450 array. EAA was calculated as the residuals from regressing gestational age on epigenetic age, calculated using a cord-blood specific epigenetic clock. Linear regression tested for associations between each maternal PFAS concentration with EAA. Bayesian kernel machine regression with hierarchical selection estimated an exposure-response function for the PFAS mixture.

RESULTS

In single pollutant models we observed an inverse relationship between perfluorodecanoate (PFDA) and EAA (-0.148 weeks per log-unit increase, 95% CI: -0.283, -0.013). Mixture analysis with hierarchical selection between perfluoroalkyl carboxylates and sulfonates indicated the carboxylates had the highest group posterior inclusion probability (PIP), or relative importance. Within this group, PFDA had the highest conditional PIP. Univariate predictor-response functions indicated PFDA and perfluorononanoate were inversely associated with EAA, while perfluorohexane sulfonate had a positive association with EAA.

CONCLUSIONS

Maternal mid-pregnancy serum concentrations of PFDA were negatively associated with EAA in cord blood, suggesting a pathway by which prenatal PFAS exposures may affect infant development. No significant associations were observed with other PFAS. Mixture models suggested opposite directions of association between perfluoroalkyl sulfonates and carboxylates. Future studies are needed to determine the importance of neonatal EAA for later child health outcomes.

Multi-ancestry epigenome-wide analyses identify methylated sites associated with aortic augmentation index in TOPMed MESA

Despite the prognostic value of arterial stiffness (AS) and pulsatile hemodynamics (PH) for cardiovascular morbidity and mortality, epigenetic modifications that contribute to AS/PH remain unknown. To gain a better understanding of the link between epigenetics (DNA methylation) and AS/PH, we examined the relationship of eight measures of AS/PH with CpG sites and co-methylated regions using multi-ancestry participants from Trans-Omics for Precision Medicine (TOPMed) Multi-Ethnic Study of Atherosclerosis (MESA) with sample sizes ranging from 438 to 874. Epigenome-wide association analysis identified one genome-wide significant CpG (cg20711926-CYP1B1) associated with aortic augmentation index (AIx). Follow-up analyses, including gene set enrichment analysis, expression quantitative trait methylation analysis, and functional enrichment analysis on differentially methylated positions and regions, further prioritized three CpGs and their annotated genes (cg23800023-ETS1, cg08426368-TGFB3, and cg17350632-HLA-DPB1) for AIx. Among these, ETS1 and TGFB3 have been previously prioritized as candidate genes. Furthermore, both ETS1 and HLA-DPB1 have significant tissue correlations between Whole Blood and Aorta in GTEx, which suggests ETS1 and HLA-DPB1 could be potential biomarkers in understanding pathophysiology of AS/PH. Overall, our findings support the possible role of epigenetic regulation via DNA methylation of specific genes associated with AIx as well as identifying potential targets for regulation of AS/PH.

Integrative genomic analyses in adipocytes implicate DNA methylation in human obesity and diabetes

DNA methylation variations are prevalent in human obesity but evidence of a causative role in disease pathogenesis is limited. Here, we combine epigenome-wide association and integrative genomics to investigate the impact of adipocyte DNA methylation variations in human obesity. We discover extensive DNA methylation changes that are robustly associated with obesity (N = 190 samples, 691 loci in subcutaneous and 173 loci in visceral adipocytes, P < 1 × 10-7). We connect obesity-associated methylation variations to transcriptomic changes at >500 target genes, and identify putative methylation-transcription factor interactions. Through Mendelian Randomisation, we infer causal effects of methylation on obesity and obesity-induced metabolic disturbances at 59 independent loci. Targeted methylation sequencing, CRISPR-activation and gene silencing in adipocytes, further identifies regional methylation variations, underlying regulatory elements and novel cellular metabolic effects. Our results indicate DNA methylation is an important determinant of human obesity and its metabolic complications, and reveal mechanisms through which altered methylation may impact adipocyte functions.

The X-factor in ART: does the use of assisted reproductive technologies influence DNA methylation on the X chromosome?

BACKGROUND

Assisted reproductive technologies (ART) may perturb DNA methylation (DNAm) in early embryonic development. Although a handful of epigenome-wide association studies of ART have been published, none have investigated CpGs on the X chromosome. To bridge this knowledge gap, we leveraged one of the largest collections of mother-father-newborn trios of ART and non-ART (natural) conceptions to date to investigate sex-specific DNAm differences on the X chromosome. The discovery cohort consisted of 982 ART and 963 non-ART trios from the Norwegian Mother, Father, and Child Cohort Study (MoBa). To verify our results from the MoBa cohort, we used an external cohort of 149 ART and 58 non-ART neonates from the Australian 'Clinical review of the Health of adults conceived following Assisted Reproductive Technologies' (CHART) study. The Illumina EPIC array was used to measure DNAm in both datasets. In the MoBa cohort, we performed a set of X-chromosome-wide association studies ('XWASs' hereafter) to search for sex-specific DNAm differences between ART and non-ART newborns. We tested several models to investigate the influence of various confounders, including parental DNAm. We also searched for differentially methylated regions (DMRs) and regions of co-methylation flanking the most significant CpGs. Additionally, we ran an analogous model to our main model on the external CHART dataset.

RESULTS

In the MoBa cohort, we found more differentially methylated CpGs and DMRs in girls than boys. Most of the associations persisted after controlling for parental DNAm and other confounders. Many of the significant CpGs and DMRs were in gene-promoter regions, and several of the genes linked to these CpGs are expressed in tissues relevant for both ART and sex (testis, placenta, and fallopian tube). We found no support for parental DNAm-dependent features as an explanation for the observed associations in the newborns. The most significant CpG in the boys-only analysis was in UBE2DNL, which is expressed in testes but with unknown function. The most significant CpGs in the girls-only analysis were in EIF2S3 and AMOT. These three loci also displayed differential DNAm in the CHART cohort.

CONCLUSIONS

Genes that co-localized with the significant CpGs and DMRs associated with ART are implicated in several key biological processes (e.g., neurodevelopment) and disorders (e.g., intellectual disability and autism). These connections are particularly compelling in light of previous findings indicating that neurodevelopmental outcomes differ in ART-conceived children compared to those naturally conceived.

Examining the immunoepigenetic-gut microbiome axis in the context of self-esteem among Native Hawaiians and other Pacific Islanders

Introduction: Native Hawaiian and other Pacific Islander (NHPI) populations experience higher rates of immunometabolic diseases compared to other racial-ethnic groups in Hawaii. As annual NHPI mortality rates for suicide and type 2 diabetes mellitus (T2DM) exceed those of the state as a whole, understanding the social and biological mechanisms underlying these disparities are urgently needed to enable preventive strategies. Methods: A community-based approach was used to investigate the immunoepigenetic-gut microbiome axis in an NHPI-enriched cohort of Oahu residents (N = 68). Self-esteem (SE) data was collected using a modified Rosenberg self-esteem (SE) assessment as a proxy measure for mental wellbeing in consideration for cultural competency. T2DM status was evaluated using point-of-care A1c (%) tests. Stool samples were collected for 16s-based metagenomic sequencing analyses. Plasma from blood samples were isolated by density-gradient centrifugation. Peripheral blood mononuclear cells (PBMCs) were collected from the same samples and enriched for monocytes using negative selection techniques. Flow-cytometry was used for immunoprofiling assays. Monocyte DNA was extracted for Illumina EPIC array-based methylation analysis. Results: Compared to individuals with normal SE (NSE), those with low SE (LSE) exhibited significantly higher plasma concentrations (pg/ml) of proinflammatory cytokines IL-8 (p = 0.051) and TNF-α (p = 0.011). Metagenomic analysis revealed that the relative abundance (%) of specific gut bacteria significantly differed between SE groups - some of which directly correlated with SE scores. Gene ontology analysis revealed that 104 significantly differentially methylated loci (DML) between SE groups were preferentially located at genes involved in immunometabolic processes. Horvath clock analyses indicated epigenetic age (Epi-Age) deceleration in individuals with LSE and acceleration in individuals with NSE (p = 0.042), yet was not reproduced by other clocks. Discussion: These data reveal novel differences in the immunoepigenetic-gut microbiome axis with respect to SE, warranting further investigation into its relationship to brain activity and mental health in NHPI. Unexpected results from Epi-Age analyses warrant further investigation into the relationship between biological age and disparate health outcomes among the NHPI population. The modifiable component of epigenetic processes and the gut microbiome makes this axis an attractive target for potential therapeutics, biomarker discovery, and novel prevention strategies.

Maternal obesity and gestational diabetes reprogram the methylome of offspring beyond birth by inducing epigenetic signatures in metabolic and developmental pathways

BACKGROUND

Obesity is a negative chronic metabolic health condition that represents an additional risk for the development of multiple pathologies. Epidemiological studies have shown how maternal obesity or gestational diabetes mellitus during pregnancy constitute serious risk factors in relation to the appearance of cardiometabolic diseases in the offspring. Furthermore, epigenetic remodelling may help explain the molecular mechanisms that underlie these epidemiological findings. Thus, in this study we explored the DNA methylation landscape of children born to mothers with obesity and gestational diabetes during their first year of life.

METHODS

We used Illumina Infinium MethylationEPIC BeadChip arrays to profile more than 770,000 genome-wide CpG sites in blood samples from a paediatric longitudinal cohort consisting of 26 children born to mothers who suffered from obesity or obesity with gestational diabetes mellitus during pregnancy and 13 healthy controls (measurements taken at 0, 6 and 12 month; total N = 90). We carried out cross-sectional and longitudinal analyses to derive DNA methylation alterations associated with developmental and pathology-related epigenomics.

RESULTS

We identified abundant DNA methylation changes during child development from birth to 6 months and, to a lesser extent, up to 12 months of age. Using cross-sectional analyses, we discovered DNA methylation biomarkers maintained across the first year of life that could discriminate children born to mothers who suffered from obesity or obesity with gestational diabetes. Importantly, enrichment analyses suggested that these alterations constitute epigenetic signatures that affect genes and pathways involved in the metabolism of fatty acids, postnatal developmental processes and mitochondrial bioenergetics, such as CPT1B, SLC38A4, SLC35F3 and FN3K. Finally, we observed evidence of an interaction between developmental DNA methylation changes and maternal metabolic condition alterations.

CONCLUSIONS

Our observations highlight the first six months of development as being the most crucial for epigenetic remodelling. Furthermore, our results support the existence of systemic intrauterine foetal programming linked to obesity and gestational diabetes that affects the childhood methylome beyond birth, which involves alterations related to metabolic pathways, and which may interact with ordinary postnatal development programmes.

CeDAR: incorporating cell type hierarchy improves cell type-specific differential analyses in bulk omics data

Bulk high-throughput omics data contain signals from a mixture of cell types. Recent developments of deconvolution methods facilitate cell type-specific inferences from bulk data. Our real data exploration suggests that differential expression or methylation status is often correlated among cell types. Based on this observation, we develop a novel statistical method named CeDAR to incorporate the cell type hierarchy in cell type-specific differential analyses of bulk data. Extensive simulation and real data analyses demonstrate that this approach significantly improves the accuracy and power in detecting cell type-specific differential signals compared with existing methods, especially in low-abundance cell types.

Measuring the long arm of childhood in real-time: Epigenetic predictors of BMI and social determinants of health across childhood and adolescence

Children who are socioeconomically disadvantaged are at increased risk for high body mass index (BMI) and multiple diseases in adulthood. The developmental origins of health and disease hypothesis proposes that early life conditions affect later-life health in a manner that is only partially modifiable by later-life experiences. Epigenetic mechanisms may regulate the influence of early life conditions on later life health. Recent epigenetic studies of adult blood samples have identified DNA-methylation sites associated with higher BMI and worse health (epigenetic-BMI). Here, we used longitudinal and twin study designs to examine whether epigenetic predictors of BMI developed in adults are valid biomarkers of child BMI and are sensitive to early life social determinants of health. Salivary epigenetic-BMI was calculated from two samples: (1) N=1,183 8-to-19-year-olds (609 female, mean age=13.4) from the Texas Twin Project (TTP), and (2) N=2,020 children (1,011 female) measured at 9 and 15 years from the Future of Families and Child Well-Being Study (FFCWS). We found that salivary epigenetic-BMI is robustly associated with children's BMI (r=0.36 to r=0.50). Longitudinal analysis suggested that epigenetic-BMI is highly stable across adolescence, but remains both a leading and lagging indicator of BMI change. Twin analyses showed that epigenetic-BMI captures differences in BMI between monozygotic twins. Moreover, children from more disadvantaged socioeconomic status (SES) and marginalized race/ethnic groups had higher epigenetic-BMI, even when controlling for concurrent BMI, pubertal development, and tobacco exposure. SES at birth relative to concurrent SES best predicted epigenetic-BMI in childhood and adolescence. We show for the first time that epigenetic predictors of BMI calculated from pediatric saliva samples are valid biomarkers of childhood BMI that are sensitive to social inequalities. Our findings are in line with the hypothesis that early life conditions are especially important factors in epigenetic regulation of later life health. Research showing that health later in life is linked to early life conditions have important implications for the development of early-life interventions that could significantly extend healthy life span.

Long-term Temporal Stability of Peripheral Blood DNA Methylation Profiles in Patients With Inflammatory Bowel Disease

BACKGROUND & AIMS

There is great current interest in the potential application of DNA methylation alterations in peripheral blood leukocytes (PBLs) as biomarkers of susceptibility, progression, and treatment response in inflammatory bowel disease (IBD). However, the intra-individual stability of PBL methylation in IBD has not been characterized. Here, we studied the long-term stability of all probes located on the Illumina HumanMethylation EPIC BeadChip array.

METHODS

We followed a cohort of 46 adult patients with IBD (36 Crohn's disease [CD], 10 ulcerative colitis [UC]; median age, 44 years; interquartile range [IQR] 27-56 years; 50% female) that received standard care follow-up at the Amsterdam University Medical Centers. Paired PBL samples were collected at 2 time points with a median of 7 years (range, 2-9 years) in between. Differential methylation and intra-class correlation (ICC) analyses were used to identify time-associated differences and temporally stable CpGs, respectively.

RESULTS

Around 60% of all EPIC array loci presented poor intra-individual stability (ICC <0.50); 78.114 (≈9%) showed good (ICC, 0.75-0.89), and 41.274 (≈5%) showed excellent (ICC ≥0.90) stability, between both measured time points. Focusing on previously identified consistently differentially methylated positions indicated that 22 CD-, 11 UC-, and 24 IBD-associated loci demonstrated high stability (ICC ≥0.75) over time; of these, we observed a marked stability of CpG loci associated to the HLA genes.

CONCLUSIONS

Our data provide insight into the long-term stability of the PBL DNA methylome within an IBD context, facilitating the selection of biologically relevant and robust IBD-associated epigenetic biomarkers with increased potential for independent validation. These data also have potential implications in understanding disease pathogenesis.

Using epigenomics to understand cellular responses to environmental influences in diseases

It is a generally accepted model that environmental influences can exert their effects, at least in part, by changing the molecular regulators of transcription that are described as epigenetic. As there is biochemical evidence that some epigenetic regulators of transcription can maintain their states long term and through cell division, an epigenetic model encompasses the idea of maintenance of the effect of an exposure long after it is no longer present. The evidence supporting this model is mostly from the observation of alterations of molecular regulators of transcription following exposures. With the understanding that the interpretation of these associations is more complex than originally recognised, this model may be oversimplistic; therefore, adopting novel perspectives and experimental approaches when examining how environmental exposures are linked to phenotypes may prove worthwhile. In this review, we have chosen to use the example of nonalcoholic fatty liver disease (NAFLD), a common, complex human disease with strong environmental and genetic influences. We describe how epigenomic approaches combined with emerging functional genetic and single-cell genomic techniques are poised to generate new insights into the pathogenesis of environmentally influenced human disease phenotypes exemplified by NAFLD.

Gaseous air pollutants and DNA methylation in a methylome-wide association study of an ethnically and environmentally diverse population of U.S. adults

Epigenetic mechanisms may underlie air pollution-health outcome associations. We estimated gaseous air pollutant-DNA methylation (DNAm) associations using twelve subpopulations within Women's Health Initiative (WHI) and Atherosclerosis Risk in Communities (ARIC) cohorts (n = 8397; mean age 61.3 years; 83% female; 46% African-American, 46% European-American, 8% Hispanic/Latino). We used geocoded participant address-specific mean ambient carbon monoxide (CO), nitrogen oxides (NO2; NOx), ozone (O3), and sulfur dioxide (SO2) concentrations estimated over the 2-, 7-, 28-, and 365-day periods before collection of blood samples used to generate Illumina 450 k array leukocyte DNAm measurements. We estimated methylome-wide, subpopulation- and race/ethnicity-stratified pollutant-DNAm associations in multi-level, linear mixed-effects models adjusted for sociodemographic, behavioral, meteorological, and technical covariates. We combined stratum-specific estimates in inverse variance-weighted meta-analyses and characterized significant associations (false discovery rate; FDR<0.05) at Cytosine-phosphate-Guanine (CpG) sites without among-strata heterogeneity (PCochran's Q > 0.05). We attempted replication in the Cooperative Health Research in Region of Augsburg (KORA) study and Normative Aging Study (NAS). We observed a -0.3 (95% CI: -0.4, -0.2) unit decrease in percent DNAm per interquartile range (IQR, 7.3 ppb) increase in 28-day mean NO2 concentration at cg01885635 (chromosome 3; regulatory region 290 bp upstream from ZNF621; FDR = 0.03). At intragenic sites cg21849932 (chromosome 20; LIME1; intron 3) and cg05353869 (chromosome 11; KLHL35; exon 2), we observed a -0.3 (95% CI: -0.4, -0.2) unit decrease (FDR = 0.04) and a 1.2 (95% CI: 0.7, 1.7) unit increase (FDR = 0.04), respectively, in percent DNAm per IQR (17.6 ppb) increase in 7-day mean ozone concentration. Results were not fully replicated in KORA and NAS. We identified three CpG sites potentially susceptible to gaseous air pollution-induced DNAm changes near genes relevant for cardiovascular and lung disease. Further harmonized investigations with a range of gaseous pollutants and averaging durations are needed to determine the effect of gaseous air pollutants on DNA methylation and ultimately gene expression.

An epigenome-wide association study of insulin resistance in African Americans

Background

African Americans have a high risk for type 2 diabetes (T2D) and insulin resistance. Studies among other population groups have identified DNA methylation loci associated with insulin resistance, but data in African Americans are lacking. Using DNA methylation profiles of blood samples obtained from the Illumina Infinium® HumanMethylation450 BeadChip, we performed an epigenome-wide association study to identify DNA methylation loci associated with insulin resistance among 136 non-diabetic, unrelated African American men (mean age 41.6 years) from the Howard University Family Study.

Results

We identified three differentially methylated positions (DMPs) for homeostatic model assessment of insulin resistance (HOMA-IR) at 5% FDR. One DMP (cg14013695, HOXA5) is a known locus among Mexican Americans, while the other two DMPs are novel—cg00456326 (OSR1; beta = 0.027) and cg20259981 (ST18; beta = 0.010). Although the cg00456326 DMP is novel, the OSR1 gene has previously been found associated with both insulin resistance and T2D in Europeans. The genes HOXA5 and ST18 have been implicated in biological processes relevant to insulin resistance. Differential methylation at the significant HOXA5 and OSR1 DMPs is associated with differences in gene expression in the iMETHYL database. Analysis of differentially methylated regions (DMRs) did not identify any epigenome-wide DMRs for HOMA-IR. We tested transferability of HOMA-IR associated DMPs from five previous EWAS in Mexican Americans, Indian Asians, Europeans, and European ancestry Americans. Out of the 730 previously reported HOMA-IR DMPs, 47 (6.4%) were associated with HOMA-IR in this cohort of African Americans.

Conclusions

The findings from our study suggest substantial differences in DNA methylation patterns associated with insulin resistance across populations. Two of the DMPs we identified in African Americans have not been reported in other populations, and we found low transferability of HOMA-IR DMPs reported in other populations in African Americans. More work in African-ancestry populations is needed to confirm our findings as well as functional analyses to understand how such DNA methylation alterations contribute to T2D pathology.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-022-01309-4.

EWAS of post-COVID-19 patients shows methylation differences in the immune-response associated gene, IFI44L, three months after COVID-19 infection

Although substantial progress has been made in managing COVID-19, it is still difficult to predict a patient’s prognosis. We explored the epigenetic signatures of COVID-19 in peripheral blood using data from an ongoing prospective observational study of COVID-19 called the Norwegian Corona Cohort Study. A series of EWASs were performed to compare the DNA methylation profiles between COVID-19 cases and controls three months post-infection. We also investigated differences associated with severity and long-COVID. Three CpGs—cg22399236, cg03607951, and cg09829636—were significantly hypomethylated (FDR < 0.05) in COVID-19 positive individuals. cg03607951 is located in IFI44L which is involved in innate response to viral infection and several systemic autoimmune diseases. cg09829636 is located in ANKRD9, a gene implicated in a wide variety of cellular processes, including the degradation of IMPDH2. The link between ANKRD9 and IMPDH2 is striking given that IMPDHs are considered therapeutic targets for COVID-19. Furthermore, gene ontology analyses revealed pathways involved in response to viruses. The lack of significant differences associated with severity and long-COVID may be real or reflect limitations in sample size. Our findings support the involvement of interferon responsive genes in the pathophysiology of COVID-19 and indicate a possible link to systemic autoimmune diseases.

Further Introduction of DNA Methylation (DNAm) Arrays in Regular Diagnostics

Methylation tests have been used for decades in regular DNA diagnostics focusing primarily on Imprinting disorders or specific loci annotated to specific disease associated gene promotors. With the introduction of DNA methylation (DNAm) arrays such as the Illumina Infinium HumanMethylation450 Beadchip array or the Illumina Infinium Methylation EPIC Beadchip array (850 k), it has become feasible to study the epigenome in a timely and cost-effective way. This has led to new insights regarding the complexity of well-studied imprinting disorders such as the Beckwith Wiedemann syndrome, but it has also led to the introduction of tests such as EpiSign, implemented as a diagnostic test in which a single array experiment can be compared to databases with known episignatures of multiple genetic disorders, especially neurodevelopmental disorders. The successful use of such DNAm tests is rapidly expanding. More and more disorders are found to be associated with discrete episignatures which enables fast and definite diagnoses, as we have shown. The first examples of environmentally induced clinical disorders characterized by discrete aberrant DNAm are discussed underlining the broad application of DNAm testing in regular diagnostics. Here we discuss exemplary findings in our laboratory covering this broad range of applications and we discuss further use of DNAm tests in the near future.

Update of gene expression/methylation and MiRNA profiling in colorectal cancer; application in diagnosis, prognosis, and targeted therapy

Background

Colorectal cancer is one of the most deadliest malignancies worldwide. Due to the dearth of appropriate biomarkers, the diagnosis of this mortal disease is usually deferred, in its turn, culminating in the failure of prevention. By the same token, proper biomarkers are at play in determining the quality of prognosis. In other words, the survival rate is contingent upon the regulation of such biomarkers.

Materials and methods

The information regarding expression (GSE41258, and GSE31905), methylation (GSE101764), and miRNA (dbDEMC) were downloaded. MEXPRESS and GEPIA confirmed the validated differentially expressed/methylated genes using TCGA data. Taking advantage of the correlation plots and receiver-operating-characteristic (ROC) curves, expression and methylation profiles were compared. The interactions between validated differentially expressed genes and differentially expressed miRNA were recognized and visualized by miRTarBase and Cytoscape, respectively. Then, the protein-protein interaction (PPI) network and hub genes were established via STRING and Cytohubba plugin. Utilizing R packages (DOSE, Enrichplot, and clusterProfiler) and DAVID database, the Functional Enrichment analysis and the detection of KEGG pathways were performed. Ultimately, in order to recognize the prognostic value of found biomarkers, they were evaluated through drawing survival plots for CRC patients.

Results

In this research, we found an expression profile (with 13 novel genes), a methylation profile (with two novel genes), and a miRNA profile with diagnostic value. Concerning diagnosis, the expression profile was evaluated more powerful in comparison with the methylation profile. Furthermore, a prognosis-related expression profile was detected.

Conclusion

In addition to diagnostic- and prognostic-applicability, the discerned profiles can assist in targeted therapy and current therapeutic strategies.

Characterization of the Illumina EPIC Array for Optimal Applications in Epigenetic Research Targeting Diverse Human Populations

The Illumina EPIC array is widely used for high-throughput profiling of DNA cytosine modifications in human samples, covering more than 850,000 modification sites across various genomic features. The application of this platform is expected to provide novel insights into the epigenetic contribution to human complex traits and diseases. Considering the diverse inter-population genetic and epigenetic variation, it will benefit the research community with a comprehensive characterization of this platform for its applicability to major global populations. Specifically, we mapped 866,836 CpG probes from the EPIC array to the human genome reference. We detected 91,034 CpG probes that did not align reliably to the human genome reference. In addition, 21,256 CpG probes were found to ambiguously map to multiple loci in the human genome, and 448 probes showing inaccurate genomic information from the original Illumina annotations. We further characterized those uniquely mapped CpG probes in terms of whether they contained common genetic variants, i.e., single nucleotide polymorphisms (SNPs), in major global populations, by utilizing the 1000 Genomes Project data. A list of optimal CpG probes on the EPIC array was generated for major global populations, with the aim of providing a resource to facilitate future studies of diverse human populations. In conclusion, our analysis indicated that studies of diverse human populations using the EPIC array would be benefited by taking into account of the technical features of this platform.

Methylomic and transcriptomic predictors of one-month exposure to cortisol in healthy individuals

Allostatic load (AL) refers to the cumulative "wear and tear" on an organism throughout its lifetime. One of the primary contributing factors to AL is prolonged exposure to stress or its primary catabolic agent cortisol. Chronic exposure to stress or cortisol is associated with numerous diseases, including cardiovascular disease, metabolic disorders, and psychiatric disorders. Therefore, a molecular marker capable of integrating a past history of cortisol exposure would be of great utility for assessing disease risk. To this end, we recruited 87 healthy males and females of European ancestry between 18 and 60 years old, extracted genomic DNA and RNA from leukocytes, and implemented a gene-centric DNA enrichment method coupled with bisulfite sequencing and RNA-Seq of total RNA for the determination of genome-wide methylation and gene transcription, respectively. Sequencing data were analyzed against awakening and bedtime cortisol data to identify differentially methylated regions (DMRs) and CpGs (DMCs) and differentially expressed genes (DEGs). Six candidate DMCs (p_unadjusted < 0.005) and nine DEGs (p_unadjusted < 0.0005) were used to construct a prediction model that could capture past 30+ days of both bedtime and awakening cortisol levels. Utilizing a cross-validation approach, we obtained a regression coefficient of R² = 0.308 for predicting continuous awakening cortisol and an area under the curve (AUC) = 0.753 for dichotomous (high vs. low tertile) awakening cortisol, and R² = 0.224 and AUC = 0.723 for continuous and dichotomous bedtime cortisol levels, respectively. To our knowledge, the current study represents the first attempt to identify genome-wide predictors of cortisol exposure that utilizes both methylation and transcription targets. The utility of our approach needs to be replicated in an independent cohort of samples for which similar cortisol metrics are available.

Revisiting genetic artifacts on DNA methylation microarrays exposes novel biological implications

BACKGROUND

Illumina DNA methylation microarrays enable epigenome-wide analysis vastly used for the discovery of novel DNA methylation variation in health and disease. However, the microarrays' probe design cannot fully consider the vast human genetic diversity, leading to genetic artifacts. Distinguishing genuine from artifactual genetic influence is of particular relevance in the study of DNA methylation heritability and methylation quantitative trait loci. But despite its importance, current strategies to account for genetic artifacts are lagging due to a limited mechanistic understanding on how such artifacts operate.

RESULTS

To address this, we develop and benchmark UMtools, an R-package containing novel methods for the quantification and qualification of genetic artifacts based on fluorescence intensity signals. With our approach, we model and validate known SNPs/indels on a genetically controlled dataset of monozygotic twins, and we estimate minor allele frequency from DNA methylation data and empirically detect variants not included in dbSNP. Moreover, we identify examples where genetic artifacts interact with each other or with imprinting, X-inactivation, or tissue-specific regulation. Finally, we propose a novel strategy based on co-methylation that can discern between genetic artifacts and genuine genomic influence.

CONCLUSIONS

We provide an atlas to navigate through the huge diversity of genetic artifacts encountered on DNA methylation microarrays. Overall, our study sets the ground for a paradigm shift in the study of the genetic component of epigenetic variation in DNA methylation microarrays.

Characterizing the properties of bisulfite sequencing data: maximizing power and sensitivity to identify between-group differences in DNA methylation

BACKGROUND

The combination of sodium bisulfite treatment with highly-parallel sequencing is a common method for quantifying DNA methylation across the genome. The power to detect between-group differences in DNA methylation using bisulfite-sequencing approaches is influenced by both experimental (e.g. read depth, missing data and sample size) and biological (e.g. mean level of DNA methylation and difference between groups) parameters. There is, however, no consensus about the optimal thresholds for filtering bisulfite sequencing data with implications for the reproducibility of findings in epigenetic epidemiology.

RESULTS

We used a large reduced representation bisulfite sequencing (RRBS) dataset to assess the distribution of read depth across DNA methylation sites and the extent of missing data. To investigate how various study variables influence power to identify DNA methylation differences between groups, we developed a framework for simulating bisulfite sequencing data. As expected, sequencing read depth, group size, and the magnitude of DNA methylation difference between groups all impacted upon statistical power. The influence on power was not dependent on one specific parameter, but reflected the combination of study-specific variables. As a resource to the community, we have developed a tool, POWEREDBiSeq, which utilizes our simulation framework to predict study-specific power for the identification of DNAm differences between groups, taking into account user-defined read depth filtering parameters and the minimum sample size per group.

CONCLUSIONS

Our data-driven approach highlights the importance of filtering bisulfite-sequencing data by minimum read depth and illustrates how the choice of threshold is influenced by the specific study design and the expected differences between groups being compared. The POWEREDBiSeq tool, which can be applied to different types of bisulfite sequencing data (e.g. RRBS, whole genome bisulfite sequencing (WGBS), targeted bisulfite sequencing and amplicon-based bisulfite sequencing), can help users identify the level of data filtering needed to optimize power and aims to improve the reproducibility of bisulfite sequencing studies.

Genome-wide DNA methylation analysis on C-reactive protein among Ghanaians suggests molecular links to the emerging risk of cardiovascular diseases

Molecular mechanisms at the intersection of inflammation and cardiovascular diseases (CVD) among Africans are still unknown. We performed an epigenome-wide association study to identify loci associated with serum C-reactive protein (marker of inflammation) among Ghanaians and further assessed whether differentially methylated positions (DMPs) were linked to CVD in previous reports, or to estimated CVD risk in the same population. We used the Illumina Infinium® HumanMethylation450 BeadChip to obtain DNAm profiles of blood samples in 589 Ghanaians from the RODAM study (without acute infections, not taking anti-inflammatory medications, CRP levels < 40 mg/L). We then used linear models to identify DMPs associated with CRP concentrations. Post-hoc, we evaluated associations of identified DMPs with elevated CVD risk estimated via ASCVD risk score. We also performed subset analyses at CRP levels ≤10 mg/L and replication analyses on candidate probes. Finally, we assessed for biological relevance of our findings in public databases. We subsequently identified 14 novel DMPs associated with CRP. In post-hoc evaluations, we found that DMPs in PC, BTG4 and PADI1 showed trends of associations with estimated CVD risk, we identified a separate DMP in MORC2 that was associated with CRP levels ≤10 mg/L, and we successfully replicated 65 (24%) of previously reported DMPs. All DMPs with gene annotations (13) were biologically linked to inflammation or CVD traits. We have identified epigenetic loci that may play a role in the intersection between inflammation and CVD among Ghanaians. Further studies among other Africans are needed to confirm our findings.

DNA methylation as the link between migration and the major noncommunicable diseases: the RODAM study

Aim: We assessed epigenome-wide DNA methylation (DNAm) differences between migrant and non-migrant Ghanaians. Materials & methods: We used the Illumina Infinium® HumanMethylation450 BeadChip to profile DNAm of 712 Ghanaians in whole blood. We used linear models to detect differentially methylated positions (DMPs) associated with migration. We performed multiple post hoc analyses to validate our findings. Results: We identified 13 DMPs associated with migration (delta-beta values: 0.2-4.5%). Seven DMPs in CPLX2, EIF4E3, MEF2D, TLX3, ST8SIA1, ANG and CHRM3 were independent of extrinsic genomic influences in public databases. Two DMPs in NLRC5 were associated with duration of stay in Europe among migrants. All DMPs were biologically linked to migration-related factors. Conclusion: Our findings provide the first insights into DNAm differences between migrants and non-migrants.

Combined epigenetic/genetic study identified an ALS age of onset modifier

Age at onset of amyotrophic lateral sclerosis (ALS) is highly variable (eg, 27-74 years in carriers of the G4C2-expansion in C9orf72). It might be influenced by environmental and genetic factors via the modulation of DNA methylation (DNAm) at CpG-sites. Hence, we combined an epigenetic and genetic approach to test the hypothesis that some common single nucleotide polymorphisms (SNPs) at CpG-sites (CpG-SNPs) could modify ALS age of onset. Our genome-wide DNAm analysis suggested three CpG-SNPs whose DNAm levels are significantly associated with age of onset in 249 ALS patients (q < 0.05). Next, genetic analysis validated the association of rs4970944 with age of onset in the discovery (n = 469; P = 0.025) and replication (n = 4160; P = 0.007) ALS cohorts. A meta-analysis of the cohorts combined showed that the median onset in AA-carriers is two years later than in GG-carriers (n = 4629; P = 0.0012). A similar association was observed with its tagging SNPs, implicating a 16 Kb region at the 1q21.3 locus as a modifier of ALS age of onset. Notably, rs4970944 genotypes are also associated with age of onset in C9orf72-carriers (n = 333; P = 0.025), suggesting that each A-allele delays onset by 1.6 years. Analysis of Genotype-Tissue Expression data revealed that the protective A-allele is linked with the reduced expression of CTSS in cerebellum (P = 0.00018), which is a critical brain region in the distributed neural circuits subserving motor control. CTSS encodes cathepsin S protein playing a key role in antigen presentation. In conclusion, we identified a 16 Kb locus tagged by rs4970944 as a modifier of ALS age of onset. Our findings support the role of antigen presenting processes in modulating age of onset of ALS and suggest potential drug targets (eg, CTSS). Future replication studies are encouraged to validate the link between the locus tagged by rs4970944 and age of onset in independent ALS cohorts, including different ethnic groups.

Cross-reactive probes on Illumina DNA methylation arrays: a large study on ALS shows that a cautionary approach is warranted in interpreting epigenome-wide association studies

Illumina DNA methylation arrays are a widely used tool for performing genome-wide DNA methylation analyses. However, measurements obtained from these arrays may be affected by technical artefacts that result in spurious associations if left unchecked. Cross-reactivity represents one of the major challenges, meaning that probes may map to multiple regions in the genome. Although several studies have reported on this issue, few studies have empirically examined the impact of cross-reactivity in an epigenome-wide association study (EWAS). In this paper, we report on cross-reactivity issues that we discovered in a large EWAS on the presence of the C9orf72 repeat expansion in ALS patients. Specifically, we found that that the majority of the significant probes inadvertently cross-hybridized to the C9orf72 locus. Importantly, these probes were not flagged as cross-reactive in previous studies, leading to novel insights into the extent to which cross-reactivity can impact EWAS. Our findings are particularly relevant for epigenetic studies into diseases associated with repeat expansions and other types of structural variation. More generally however, considering that most spurious associations were not excluded based on pre-defined sets of cross-reactive probes, we believe that the presented data-driven flag and consider approach is relevant for any type of EWAS.

Occurrence of Accelerated Epigenetic Aging and Methylation Disruptions in Human Immunodeficiency Virus Infection Before Antiretroviral Therapy

BACKGROUND

Whether accelerated aging develops over the course of chronic human immunodeficiency virus (HIV) infection or can be observed before significant immunosuppression on is unknown. We studied DNA methylation in blood to estimate cellular aging in persons living with HIV (PLWH) before the initiation of antiretroviral therapy (ART).

METHODS

A total of 378 ART-naive PLWH who had CD4 T-cell counts >500/µL and were enrolled in the Strategic Timing of Antiretroviral Therapy trial (Pulmonary Substudy) were compared with 34 HIV-negative controls. DNA methylation was performed using the Illumina MethylationEPIC BeadChip. Differentially methylated positions (DMPs) and differentially methylated regions (DMRs) in PLWH compared with controls were identified using a robust linear model. Methylation age was calculated using a previously described epigenetic clock.

RESULTS

There were a total of 56 639 DMPs and 6103 DMRs at a false discovery rate of <0.1. The top 5 DMPs corresponded to genes NLRC5, VRK2, B2M, and GPR6 and were highly enriched for cancer-related pathways. PLWH had significantly higher methylation age than HIV-negative controls (P = .001), with black race, low CD4 and high CD8 T-cell counts, and duration of HIV being risk factors for age acceleration.

CONCLUSIONS

PLWH before the initiation of ART and with preserved immune status show evidence of advanced methylation aging.

The DNA methylome of human sperm is distinct from blood with little evidence for tissue-consistent obesity associations

Epidemiological research suggests that paternal obesity may increase the risk of fathering small for gestational age offspring. Studies in non-human mammals indicate that such associations could be mediated by DNA methylation changes in spermatozoa that influence offspring development in utero. Human obesity is associated with differential DNA methylation in peripheral blood. It is unclear, however, whether this differential DNA methylation is reflected in spermatozoa. We profiled genome-wide DNA methylation using the Illumina MethylationEPIC array in a cross-sectional study of matched human blood and sperm from lean (discovery n = 47; replication n = 21) and obese (n = 22) males to analyse tissue covariation of DNA methylation, and identify obesity-associated methylomic signatures. We found that DNA methylation signatures of human blood and spermatozoa are highly discordant, and methylation levels are correlated at only a minority of CpG sites (~1%). At the majority of these sites, DNA methylation appears to be influenced by genetic variation. Obesity-associated DNA methylation in blood was not generally reflected in spermatozoa, and obesity was not associated with altered covariation patterns or accelerated epigenetic ageing in the two tissues. However, one cross-tissue obesity-specific hypermethylated site (cg19357369; chr4:2429884; P = 8.95 × 10-8; 2% DNA methylation difference) was identified, warranting replication and further investigation. When compared to a wide range of human somatic tissue samples (n = 5,917), spermatozoa displayed differential DNA methylation across pathways enriched in transcriptional regulation. Overall, human sperm displays a unique DNA methylation profile that is highly discordant to, and practically uncorrelated with, that of matched peripheral blood. We observed that obesity was only nominally associated with differential DNA methylation in sperm, and therefore suggest that spermatozoal DNA methylation is an unlikely mediator of intergenerational effects of metabolic traits.

Placental Epigenome-Wide Association Study Identified Loci Associated with Childhood Adiposity at 3 Years of Age

The aim of this study was to identify placental DNA methylation (DNAm) variations associated with adiposity at 3 years of age. We quantified placental DNAm using the Infinium MethylationEPIC BeadChips. We assessed associations between DNAm at single-CpGs and skinfold thickness using robust linear regression models adjusted for gestational age, child's sex, age at follow-up and cellular heterogeneity. We sought replication of DNAm association with child adiposity in an independent cohort. We quantified placental mRNA levels for annotated gene using qRT-PCR and tested for correlation with DNAm. Lower DNAm at cg22593959 and cg22436429 was associated with higher adiposity (β = -1.18, q = 0.002 and β = -0.82, q = 0.04). The cg22593959 is located in an intergenic region (chr7q31.3), whereas cg22436429 is within the TFAP2E gene (1p34.3). DNAm at cg22593959 and cg22436429 was correlated with mRNA levels at FAM3C (r_s = -0.279, p = 0.005) and TFAP2E (r_s = 0.216, p = 0.03). In an independent cohort, the association between placental DNAm at cg22593959 and childhood adiposity was of similar strength and direction (β = -3.8 ± 4.1, p = 0.36), yet non-significant. Four genomic regions were also associated with skinfold thickness within FMN1, MAGI2, SKAP2 and BMPR1B genes. We identified placental epigenetic variations associated with adiposity at 3 years of age suggesting that childhood fat accretion patterns might be established during fetal life.

"Epigenome-wide methylation profile of chronic kidney disease-derived arterial DNA uncovers novel pathways in disease-associated cardiovascular pathology."

Chronic kidney disease (CKD) related cardiovascular disease (CVD) is characterized by vascular remodelling with well-established structural and functional changes in the vascular wall such as arterial stiffness, matrix deposition, and calcification. These phenotypic changes resemble pathology seen in ageing, and are likely to be mediated by sustained alterations in gene expression, which may be caused by epigenetic changes such as tissue-specific DNA methylation. We aimed to investigate tissue specific changes in DNA methylation that occur in CKD-related CVD. Genome-wide DNA methylation changes were examined in bisulphite converted genomic DNA isolated from the vascular media of CKD and healthy arteries. Methylation-specific PCR was used to validate the array data, and the association between DNA methylation and gene and protein expression was examined. The DNA methylation age was compared to the chronological age in both cases and controls. Three hundred and nineteen differentially methylated regions (DMR) were identified spread across the genome. Pathway analysis revealed that DMRs associated with genes were involved in embryonic and vascular development, and signalling pathways such as TGFβ and FGF. Expression of top differentially methylated gene HOXA5 showed a significant negative correlation with DNA methylation. Interestingly, DNA methylation age and chronological age were highly correlated, but there was no evidence of accelerated age-related DNA methylation in the arteries of CKD patients. In conclusion, we demonstrated that differential DNA methylation in the arterial tissue of CKD patients represents a potential mediator of arterial pathology and may be used to uncover novel pathways in the genesis of CKD-associated complications.

DNA methylation study of Huntington's disease and motor progression in patients and in animal models

Although Huntington's disease (HD) is a well studied Mendelian genetic disorder, less is known about its associated epigenetic changes. Here, we characterize DNA methylation levels in six different tissues from 3 species: a mouse huntingtin (Htt) gene knock-in model, a transgenic HTT sheep model, and humans. Our epigenome-wide association study (EWAS) of human blood reveals that HD mutation status is significantly (p < 10^-7) associated with 33 CpG sites, including the HTT gene (p = 6.5 × 10^-26). These Htt/HTT associations were replicated in the Q175 Htt knock-in mouse model (p = 6.0 × 10^-8) and in the transgenic sheep model (p = 2.4 × 10^-88). We define a measure of HD motor score progression among manifest HD cases based on multiple clinical assessments. EWAS of motor progression in manifest HD cases exhibits significant (p < 10^-7) associations with methylation levels at three loci: near PEX14 (p = 9.3 × 10^-9), GRIK4 (p = 3.0 × 10^-8), and COX4I2 (p = 6.5 × 10^-8). We conclude that HD is accompanied by profound changes of DNA methylation levels in three mammalian species.

Epigenome-wide association study of DNA methylation and adult asthma in the Agricultural Lung Health Study

Epigenome-wide studies of methylation in children support a role for epigenetic mechanisms in asthma; however, studies in adults are rare and few have examined non-atopic asthma. We conducted the largest epigenome-wide association study (EWAS) of blood DNA methylation in adults in relation to non-atopic and atopic asthma.

We measured DNA methylation in blood using the Illumina MethylationEPIC array among 2286 participants in a case-control study of current adult asthma nested within a United States agricultural cohort. Atopy was defined by serum specific immunoglobulin E (IgE). Participants were categorised as atopy without asthma (n=185), non-atopic asthma (n=673), atopic asthma (n=271), or a reference group of neither atopy nor asthma (n=1157). Analyses were conducted using logistic regression.

No associations were observed with atopy without asthma. Numerous cytosine–phosphate–guanine (CpG) sites were differentially methylated in non-atopic asthma (eight at family-wise error rate (FWER) p<9×10⁻⁸, 524 at false discovery rate (FDR) less than 0.05) and implicated 382 novel genes. More CpG sites were identified in atopic asthma (181 at FWER, 1086 at FDR) and implicated 569 novel genes. 104 FDR CpG sites overlapped. 35% of CpG sites in non-atopic asthma and 91% in atopic asthma replicated in studies of whole blood, eosinophils, airway epithelium, or nasal epithelium. Implicated genes were enriched in pathways related to the nervous system or inflammation.

We identified numerous, distinct differentially methylated CpG sites in non-atopic and atopic asthma. Many CpG sites from blood replicated in asthma-relevant tissues. These circulating biomarkers reflect risk and sequelae of disease, as well as implicate novel genes associated with non-atopic and atopic asthma.

Distinct methylation signals are found in non-atopic and atopic asthma. Most are related to gene expression and are replicated in asthma-relevant tissues, confirming the value of blood DNA methylation for identifying novel genes linked in asthma pathogenesis.https://bit.ly/2VnbJg3

PRDM8 reveals aberrant DNA methylation in aging syndromes and is relevant for hematopoietic and neuronal differentiation

BACKGROUND

Dyskeratosis congenita (DKC) and idiopathic aplastic anemia (AA) are bone marrow failure syndromes that share characteristics of premature aging with severe telomere attrition. Aging is also reflected by DNA methylation changes, which can be utilized to predict donor age. There is evidence that such epigenetic age predictions are accelerated in premature aging syndromes, but it is yet unclear how this is related to telomere length. DNA methylation analysis may support diagnosis of DKC and AA, which still remains a challenge for these rare diseases.

RESULTS

In this study, we analyzed blood samples of 70 AA and 18 DKC patients to demonstrate that their epigenetic age predictions are overall increased, albeit not directly correlated with telomere length. Aberrant DNA methylation was observed in the gene PRDM8 in DKC and AA as well as in other diseases with premature aging phenotype, such as Down syndrome and Hutchinson-Gilford-Progeria syndrome. Aberrant DNA methylation patterns were particularly found within subsets of cell populations in DKC and AA samples as measured with barcoded bisulfite amplicon sequencing (BBA-seq). To gain insight into the functional relevance of PRDM8, we used CRISPR/Cas9 technology to generate induced pluripotent stem cells (iPSCs) with heterozygous and homozygous knockout. Loss of PRDM8 impaired hematopoietic and neuronal differentiation of iPSCs, even in the heterozygous knockout clone, but it did not impact on epigenetic age.

CONCLUSION

Taken together, our results demonstrate that epigenetic aging is accelerated in DKC and AA, independent from telomere attrition. Furthermore, aberrant DNA methylation in PRDM8 provides another biomarker for bone marrow failure syndromes and modulation of this gene in cellular subsets may be related to the hematopoietic and neuronal phenotypes observed in premature aging syndromes.

Genome-wide DNA methylation differences and polychlorinated biphenyl (PCB) exposure in a US population

Exposure to polychlorinated biphenyls (PCBs), an endocrine-disrupting compound, is ubiquitous despite decades-old bans on the manufacture and use of PCBs. Increased exposure to PCBs is associated with adverse health consequences throughout life, including type 2 diabetes and cancer. PCB exposure is also associated with alterations in epigenetic marks and gene transcription, which could lead to adverse health outcomes, but many of these are population-specific. To further investigate the association between PCB and epigenetic marks, DNA methylation was measured at 787,684 CpG sites in 641 peripheral blood samples from the Michigan Polybrominated Biphenyl (PBB) Registry. 1345 CpGs were associated with increased total PCB level after controlling for age, sex, and 24 surrogate variables (FDR < 0.05). These CpGs were enriched in active promoter and transcription associated regions (p < 0.05), and in regions around the binding sites for transcription factors involved in xenobiotic metabolism and immune function (FDR < 0.05). PCB exposure also associated with proportions of CD4T, NK, and granulocyte cell types, and with the neutrophil to lymphocyte ratio (NLR) (p < 0.05), and the estimated effect sizes of PCB on the epigenome were correlated with the effect sizes previously reported in an epigenome-wide study of C-reactive protein (r = 0.29; p = 2.22e-5), supporting previous studies on the association between PCB and immune dysfunction. These results indicate that PCB exposure is associated with differences in epigenetic marks in active regions of the genome, and future work should investigate whether these may mediate the association between PCB and health consequences.

Gene-methylation interactions: discovering region-wise DNA methylation levels that modify SNP-associated disease risk

BACKGROUND

Current technology allows rapid assessment of DNA sequences and methylation levels at a single-site resolution for hundreds of thousands of sites in the human genome, in thousands of individuals simultaneously. This has led to an increase in epigenome-wide association studies (EWAS) of complex traits, particularly those that are poorly explained by previous genome-wide association studies (GWAS). However, the genome and epigenome are intertwined, e.g., DNA methylation is known to affect gene expression through, for example, genomic imprinting. There is thus a need to go beyond single-omics data analyses and develop interaction models that allow a meaningful combination of information from EWAS and GWAS.

RESULTS

We present two new methods for genetic association analyses that treat offspring DNA methylation levels as environmental exposure. Our approach searches for statistical interactions between SNP alleles and DNA methylation (G ×Me) and between parent-of-origin effects and DNA methylation (PoO ×Me), using case-parent triads or dyads. We use summarized methylation levels over nearby genomic region to ease biological interpretation. The methods were tested on a dataset of parent-offspring dyads, with EWAS data on the offspring. Our results showed that methylation levels around a SNP can significantly alter the estimated relative risk. Moreover, we show how a control dataset can identify false positives.

CONCLUSIONS

The new methods, G ×Me and PoO ×Me, integrate DNA methylation in the assessment of genetic relative risks and thus enable a more comprehensive biological interpretation of genome-wide scans. Moreover, our strategy of condensing DNA methylation levels within regions helps overcome specific disadvantages of using sparse chip-based measurements. The methods are implemented in the freely available R package Haplin ( https://cran.r-project.org/package=Haplin ), enabling fast scans of multi-omics datasets.

MethylToSNP: identifying SNPs in Illumina DNA methylation array data

Background

Current array-based methods for the measurement of DNA methylation rely on the process of sodium bisulfite conversion to differentiate between methylated and unmethylated cytosine bases in DNA. In the absence of genotype data this process can lead to ambiguity in data interpretation when a sample has polymorphisms at a methylation probe site. A common way to minimize this problem is to exclude such potentially problematic sites, with some methods removing as much as 60% of array probes from consideration before data analysis.

Results

Here, we present an algorithm implemented in an R Bioconductor package, MethylToSNP, which detects a characteristic data pattern to infer sites likely to be confounded by polymorphisms. Additionally, the tool provides a stringent reliability score to allow thresholding on SNP predictions. We calibrated parameters and thresholds used by the algorithm on simulated and real methylation data sets. We illustrate findings using methylation data from YRI (Yoruba in Ibadan, Nigeria), CEPH (European descent) and KhoeSan (southern African) populations. Our polymorphism predictions made using MethylToSNP have been validated through SNP databases and bisulfite and genomic sequencing.

Conclusions

The benefits of this method are threefold. First, it prevents extensive data loss by considering only SNPs specific to the individuals in the study. Second, it offers the possibility to identify new polymorphisms in samples for which there is little known about the genetic landscape. Third, it identifies variants as they exist in functional regions of a genome, such as in CTCF (transcriptional repressor) sites and enhancers, that may be common alleles or personal mutations with potential to deleteriously affect genomic regulatory activities. We demonstrate that MethylToSNP is applicable to the Illumina 450K and Illumina 850K EPIC array data and is also backwards compatible to the 27K methylation arrays. Going forward, this kind of nuanced approach can increase the amount of information derived from precious data sets by considering samples of the project individually to enable more informed decisions about data cleaning.

Methylome-wide association study provides evidence of particulate matter air pollution-associated DNA methylation

BACKGROUND

DNA methylation (DNAm) may contribute to processes that underlie associations between air pollution and poor health. Therefore, our objective was to evaluate associations between DNAm and ambient concentrations of particulate matter (PM) ≤2.5, ≤10, and 2.5-10 μm in diameter (PM2.5; PM10; PM2.5-10).

METHODS

We conducted a methylome-wide association study among twelve cohort- and race/ethnicity-stratified subpopulations from the Women's Health Initiative and the Atherosclerosis Risk in Communities study (n = 8397; mean age: 61.5 years; 83% female; 45% African American; 9% Hispanic/Latino American). We averaged geocoded address-specific estimates of daily and monthly mean PM concentrations over 2, 7, 28, and 365 days and 1 and 12 months before exams at which we measured leukocyte DNAm in whole blood. We estimated subpopulation-specific, DNAm-PM associations at approximately 485,000 Cytosine-phosphate-Guanine (CpG) sites in multi-level, linear, mixed-effects models. We combined subpopulation- and site-specific estimates in fixed-effects, inverse variance-weighted meta-analyses, then for associations that exceeded methylome-wide significance and were not heterogeneous across subpopulations (P < 1.0 × 10-7; PCochran's Q > 0.10), we characterized associations using publicly accessible genomic databases and attempted replication in the Cooperative Health Research in the Region of Augsburg (KORA) study.

RESULTS

Analyses identified significant DNAm-PM associations at three CpG sites. Twenty-eight-day mean PM10 was positively associated with DNAm at cg19004594 (chromosome 20; MATN4; P = 3.33 × 10-8). One-month mean PM10 and PM2.5-10 were positively associated with DNAm at cg24102420 (chromosome 10; ARPP21; P = 5.84 × 10-8) and inversely associated with DNAm at cg12124767 (chromosome 7; CFTR; P = 9.86 × 10-8). The PM-sensitive CpG sites mapped to neurological, pulmonary, endocrine, and cardiovascular disease-related genes, but DNAm at those sites was not associated with gene expression in blood cells and did not replicate in KORA.

CONCLUSIONS

Ambient PM concentrations were associated with DNAm at genomic regions potentially related to poor health among racially, ethnically and environmentally diverse populations of U.S. women and men. Further investigation is warranted to uncover mechanisms through which PM-induced epigenomic changes may cause disease.

Differentially methylated regions in bipolar disorder and suicide

The addition of a methyl group to, typically, a cytosine-guanine dinucleotide (CpG) creates distinct DNA methylation patterns across the genome that can regulate gene expression. Aberrant DNA methylation of CpG sites has been associated with many psychiatric disorders including bipolar disorder (BD) and suicide. Using the SureSelect^XT system, Methyl-Seq, we investigated the DNA methylation status of CpG sites throughout the genome in 50 BD individuals (23 subjects who died by suicide and 27 subjects who died from other causes) and 31 nonpsychiatric controls. We identified differentially methylated regions (DMRs) from three analyses: (a) BD subjects compared to nonpsychiatric controls (BD-NC), (b) BD subjects who died by suicide compared to nonpsychiatric controls (BDS-NC), and (c) BDS subjects compared to BD subjects who died from other causes (BDS-BDNS). One DMR from the BDS-NC analysis, located in ARHGEF38, was significantly hypomethylated (23.4%) in BDS subjects. This finding remained significant after multiple testing (P_Bootstrapped = 9.0 × 10^-3 ), was validated using pyrosequencing, and was more significant in males. A secondary analysis utilized Ingenuity Pathway Analysis to identify enrichment in nominally significant DMRs. This identified an association with several pathways including axonal guidance signaling, calcium signaling, β-adrenergic signaling, and opioid signaling. Our comprehensive study provides further support that DNA methylation alterations influence the risk for BD and suicide. However, further investigation is required to confirm these associations and identify their functional consequences.

RF_Purify: a novel tool for comprehensive analysis of tumor-purity in methylation array data based on random forest regression

Background

With the advent of array-based techniques to measure methylation levels in primary tumor samples, systematic investigations of methylomes have widely been performed on a large number of tumor entities. Most of these approaches are not based on measuring individual cell methylation but rather the bulk tumor sample DNA, which contains a mixture of tumor cells, infiltrating immune cells and other stromal components. This raises questions about the purity of a certain tumor sample, given the varying degrees of stromal infiltration in different entities. Previous methods to infer tumor purity require or are based on the use of matching control samples which are rarely available. Here we present a novel, reference free method to quantify tumor purity, based on two Random Forest classifiers, which were trained on ABSOLUTE as well as ESTIMATE purity values from TCGA tumor samples. We subsequently apply this method to a previously published, large dataset of brain tumors, proving that these models perform well in datasets that have not been characterized with respect to tumor purity .

Results

Using two gold standard methods to infer purity – the ABSOLUTE score based on whole genome sequencing data and the ESTIMATE score based on gene expression data- we have optimized Random Forest classifiers to predict tumor purity in entities that were contained in the TCGA project. We validated these classifiers using an independent test data set and cross-compared it to other methods which have been applied to the TCGA datasets (such as ESTIMATE and LUMP).

Using Illumina methylation array data of brain tumor entities (as published in Capper et al. (Nature 555:469-474,2018)) we applied this model to estimate tumor purity and find that subgroups of brain tumors display substantial differences in tumor purity.

Conclusions

Random forest- based tumor purity prediction is a well suited tool to extrapolate gold standard measures of purity to novel methylation array datasets. In contrast to other available methylation based tumor purity estimation methods, our classifiers do not need a priori knowledge about the tumor entity or matching control tissue to predict tumor purity.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-3014-z) contains supplementary material, which is available to authorized users.

Exposure to polybrominated biphenyl and stochastic epigenetic mutations: application of a novel epigenetic approach to environmental exposure in the Michigan polybrominated biphenyl registry

Endocrine-disrupting compounds are associated with altered epigenetic regulation and adverse health outcomes, although inconsistent results suggest that people have varied responses to the same exposure. Interpersonal variation in response to environmental exposures is not identified using standard, population-based methods. However, methods that capture an individual's response, such as analyzing stochastic epigenetic mutations (SEMs), may capture currently missed effects of environmental exposure. To test whether polybrominated biphenyl (PBB) was associated with SEMs, DNA methylation was measured using Illumina's MethylationEPIC array in PBB-exposed individuals, and SEMs were identified. Association was tested using a linear regression with robust sandwich variance estimators, controlling for age, sex, lipids, and cell types. The number of SEMs was variable (range: 119-18,309), and positively associated with age (p = 1.23e-17), but not with sex (p = 0.97). PBBs and SEMs were only positively associated in people who were older when they were exposed (p = 0.02 vs. p = 0.91). Many subjects had SEMs enriched in biological pathways, particularly in pathways involved with xenobiotic metabolism and endocrine function. Higher number of SEMs was also associated with higher age acceleration (intrinsic: p = 1.70e-3; extrinsic: p = 3.59e-11), indicating that SEMs may be associated with age-related health problems. Finding an association between environmental contaminants and higher SEMs may provide insight into individual differences in response to environmental contaminants, as well as into the biological mechanism behind SEM formation. Furthermore, these results suggest that people may be particularly vulnerable to epigenetic dysregulation from environmental exposures as they age.

Kidney cytosine methylation changes improve renal function decline estimation in patients with diabetic kidney disease

Epigenetic changes might provide the biological explanation for the long-lasting impact of metabolic alterations of diabetic kidney disease development. Here we examined cytosine methylation of human kidney tubules using Illumina Infinium 450 K arrays from 91 subjects with and without diabetes and varying degrees of kidney disease using a cross-sectional design. We identify cytosine methylation changes associated with kidney structural damage and build a model for kidney function decline. We find that the methylation levels of 65 probes are associated with the degree of kidney fibrosis at genome wide significance. In total 471 probes improve the model for kidney function decline. Methylation probes associated with kidney damage and functional decline enrich on kidney regulatory regions and associate with gene expression changes, including epidermal growth factor (EGF). Altogether, our work shows that kidney methylation differences can be detected in patients with diabetic kidney disease and improve kidney function decline models indicating that they are potentially functionally important.

Patients with diabetes commonly develop diabetic kidney disease (DKD). Here Gluck et al. identify a set of probes differentially methylated in renal samples from patients with DKD, and find that inclusion of these methylation probes improves current prediction models of renal function decline.

Epigenome-wide Analysis Identifies Genes and Pathways Linked to Neurobehavioral Variation in Preterm Infants

Neonatal molecular biomarkers of neurobehavioral responses (measures of brain-behavior relationships), when combined with neurobehavioral performance measures, could lead to better predictions of long-term developmental outcomes. To this end, we examined whether variability in buccal cell DNA methylation (DNAm) associated with neurobehavioral profiles in a cohort of infants born less than 30 weeks postmenstrual age (PMA) and participating in the Neonatal Neurobehavior and Outcomes in Very Preterm Infants (NOVI) Study (N = 536). We tested whether epigenetic age, age acceleration, or DNAm levels at individual loci differed between infants based on their NICU Network Neurobehavioral Scale (NNNS) profile classifications. We adjusted for recruitment site, infant sex, PMA, and tissue heterogeneity. Infants with an optimally well-regulated NNNS profile had older epigenetic age compared to other NOVI infants (β₁ = 0.201, p-value = 0.026), but no significant difference in age acceleration. In contrast, infants with an atypical NNNS profile had differential methylation at 29 CpG sites (FDR < 10%). Some of the genes annotated to these CpGs included PLA2G4E, TRIM9, GRIK3, and MACROD2, which have previously been associated with neurological structure and function, or with neurobehavioral disorders. These findings contribute to the existing evidence that neonatal epigenetic variations may be informative for infant neurobehavior.

DNA methylation signature of smoking in lung cancer is enriched for exposure signatures in newborn and adult blood

Smoking impacts DNA methylation genome-wide in blood of newborns from maternal smoking during pregnancy and adults from personal smoking. We compared smoking-related DNA methylation in lung adenocarcinoma (61 never smokers, 91 current smokers, and 238 former smokers) quantified with the Illumina450k BeadArray in The Cancer Genome Atlas with published large consortium meta-analyses of newborn and adult blood. We assessed whether CpG sites related to smoking in blood from newborns and adults were enriched in the lung adenocarcinoma methylation signal. Testing CpGs differentially methylated by smoke exposure, we identified 296 in lung adenocarcinoma meeting a P < 10-4 cutoff, while previous meta-analyses identified 3,042 in newborn blood, and 8,898 in adult blood meeting the same P < 10-4 cutoff. Lung signals were highly enriched for those seen in newborn (24 overlapping CpGs, Penrichment = 1.2 × 10-18) and adult blood (66 overlapping CpGs, Penrichment = 1.2 × 10-48). The 105 genes annotated to CpGs differentially methylated in lung tumors, but not blood, were enriched for RNA processing ontologies. Some epigenetic alterations associated with cigarette smoke exposure are tissue specific, but others are common across tissues. These findings support the value of blood-based methylation biomarkers for assessing exposure effects in target tissues.

Whole-genome methylation profiling of the retinal pigment epithelium of individuals with age-related macular degeneration reveals differential methylation of the SKI, GTF2H4, and TNXB genes

Background

Age-related macular degeneration (AMD) is a degenerative disorder of the central retina and the foremost cause of blindness. The retinal pigment epithelium (RPE) is a primary site of disease pathogenesis. The genetic basis of AMD is relatively well understood; however, this knowledge is yet to yield a treatment for the most prevalent non-neovascular disease forms. Therefore, tissue-specific epigenetic mechanisms of gene regulation are of considerable interest in AMD. We aimed to identify differentially methylated genes associated with AMD in the RPE and differentiate local DNA methylation aberrations from global DNA methylation changes, as local DNA methylation changes may be more amenable to therapeutic manipulation.

Methods

Epigenome-wide association study and targeted gene expression profiling were carried out in RPE cells from eyes of human donors. We performed genome-wide DNA methylation profiling (Illumina 450k BeadChip array) on RPE cells from 44 human donor eyes (25 AMD and 19 normal controls). We validated the findings using bisulfite pyrosequencing in 55 RPE samples (30 AMD and 25 normal controls) including technical (n = 38) and independent replicate samples (n = 17). Long interspersed nucleotide element 1 (LINE-1) analysis was then applied to assess global DNA methylation changes in the RPE. RT-qPCR on independent donor RPE samples was performed to assess gene expression changes.

Results

Genome-wide DNA methylation profiling identified differential methylation of multiple loci including the SKI proto-oncogene (SKI) (p = 1.18 × 10⁻⁹), general transcription factor IIH subunit H4 (GTF2H4) (p = 7.03 × 10⁻⁷), and Tenascin X (TNXB) (p = 6.30 × 10⁻⁶) genes in AMD. Bisulfite pyrosequencing validated the differentially methylated locus cg18934822 in SKI, and cg22508626 within GTF2H4, and excluded global DNA methylation changes in the RPE in AMD. We further demonstrated the differential expression of SKI, GTF2H4, and TNXB in the RPE of independent AMD donors.

Conclusions

We report the largest genome-wide methylation analysis of RPE in AMD along with associated gene expression changes to date, for the first-time reaching genome-wide significance, and identified novel targets for functional and future therapeutic intervention studies. The novel differentially methylated genes SKI and GTF2H4 have not been previously associated with AMD, and regulate disease pathways implicated in AMD, including TGF beta signaling (SKI) and transcription-dependent DNA repair mechanisms (GTF2H4).

Electronic supplementary material

The online version of this article (10.1186/s13148-019-0608-2) contains supplementary material, which is available to authorized users.

Epigenome-wide meta-analysis of DNA methylation and childhood asthma

BACKGROUND

Epigenetic mechanisms, including methylation, can contribute to childhood asthma. Identifying DNA methylation profiles in asthmatic patients can inform disease pathogenesis.

OBJECTIVE

We sought to identify differential DNA methylation in newborns and children related to childhood asthma.

METHODS

Within the Pregnancy And Childhood Epigenetics consortium, we performed epigenome-wide meta-analyses of school-age asthma in relation to CpG methylation (Illumina450K) in blood measured either in newborns, in prospective analyses, or cross-sectionally in school-aged children. We also identified differentially methylated regions.

RESULTS

In newborns (8 cohorts, 668 cases), 9 CpGs (and 35 regions) were differentially methylated (epigenome-wide significance, false discovery rate < 0.05) in relation to asthma development. In a cross-sectional meta-analysis of asthma and methylation in children (9 cohorts, 631 cases), we identified 179 CpGs (false discovery rate < 0.05) and 36 differentially methylated regions. In replication studies of methylation in other tissues, most of the 179 CpGs discovered in blood replicated, despite smaller sample sizes, in studies of nasal respiratory epithelium or eosinophils. Pathway analyses highlighted enrichment for asthma-relevant immune processes and overlap in pathways enriched both in newborns and children. Gene expression correlated with methylation at most loci. Functional annotation supports a regulatory effect on gene expression at many asthma-associated CpGs. Several implicated genes are targets for approved or experimental drugs, including IL5RA and KCNH2.

CONCLUSION

Novel loci differentially methylated in newborns represent potential biomarkers of risk of asthma by school age. Cross-sectional associations in children can reflect both risk for and effects of disease. Asthma-related differential methylation in blood in children was substantially replicated in eosinophils and respiratory epithelium.

Genetic variants influence on the placenta regulatory landscape

From genomic association studies, quantitative trait loci analysis, and epigenomic mapping, it is evident that significant efforts are necessary to define genetic-epigenetic interactions and understand their role in disease susceptibility and progression. For this reason, an analysis of the effects of genetic variation on gene expression and DNA methylation in human placentas at high resolution and whole-genome coverage will have multiple mechanistic and practical implications. By producing and analyzing DNA sequence variation (n = 303), DNA methylation (n = 303) and mRNA expression data (n = 80) from placentas from healthy women, we investigate the regulatory landscape of the human placenta and offer analytical approaches to integrate different types of genomic data and address some potential limitations of current platforms. We distinguish two profiles of interaction between expression and DNA methylation, revealing linear or bimodal effects, reflecting differences in genomic context, transcription factor recruitment, and possibly cell subpopulations. These findings help to clarify the interactions of genetic, epigenetic, and transcriptional regulatory mechanisms in normal human placentas. They also provide strong evidence for genotype-driven modifications of transcription and DNA methylation in normal placentas. In addition to these mechanistic implications, the data and analytical methods presented here will improve the interpretability of genome-wide and epigenome-wide association studies for human traits and diseases that involve placental functions.

The placenta is a critical organ playing multiple roles including oxygen and metabolite transfer from mother to fetus, hormone production, and vascular perfusion. With this study, we aimed to deliver a placenta-specific regulatory map based on a combination of publicly available and newly generated data. To complete this reference, we obtained genotype information (n = 303), DNA methylation (n = 303) and expression data (n = 80) for placentas from healthy women. Our analysis of methylation and expression quantitative trait loci (QTLs) and correlations between methylation and expression data were designed to identify fundamental associations between genome, transcriptome, and epigenome in this key fetal organ. The results provide high-resolution genetic and epigenetic maps specific to the placenta based on a representative ethnically diverse cohort. As interest and efforts are growing to better understand the etiology of placental disease and the impact of the environment on placental function these data will provide a reference and enhance future investigations.