RELIC: a novel dye-bias correction method for Illumina Methylation BeadChip

Fetal programming by the parental microbiome of offspring behavior, and DNA methylation and gene expression within the hippocampus

Background

The microorganisms colonizing the gastrointestinal tract of animals, collectively referred to as the gut microbiome, affect numerous host behaviors dependent on the central nervous system (CNS). Studies comparing germ-free mice to normally colonized mice have demonstrated influences of the microbiome on anxiety-related behaviors, voluntary activity, and gene expression in the CNS. Additionally, there is epidemiologic evidence supporting an intergenerational influence of the maternal microbiome on neurodevelopment of offspring and behavior later in life. There is limited experimental evidence however directly linking the maternal microbiome to long-term neurodevelopmental outcomes, or knowledge regarding mechanisms responsible for such effects.

Results

Here we show that that the maternal microbiome has a dominant influence on several offspring phenotypes including anxiety-related behavior, voluntary activity, and body weight. Adverse outcomes in offspring were associated with features of the maternal microbiome including bile salt hydrolase activity gene expression (bsh), abundance of certain bile acids, and hepatic expression of Slc10a1. In cross-foster experiments, offspring resembled their birth dam phenotypically, despite faithful colonization in the postnatal period with the surrogate dam microbiome. Genome-wide methylation analysis of hippocampal DNA identified microbiome-associated differences in methylation of 196 loci in total, 176 of which show conserved profiles between mother and offspring. Further, single-cell transcriptional analysis revealed accompanying differences in expression of several differentially methylated genes within certain hippocampal cell clusters, and vascular expression of genes associated with bile acid transport. Inferred cell-to-cell communication in the hippocampus based on coordinated ligand-receptor expression revealed differences in expression of neuropeptides associated with satiety.

Conclusions

Collectively, these data provide proof-of-principle that the maternal gut microbiome has a dominant influence on the neurodevelopment underlying certain offspring behaviors and activities, and selectively affects genome methylation and gene expression in the offspring CNS in conjunction with that neurodevelopment.

Epigenome-wide association study of fluoride exposure during early adolescence and DNA methylation among U.S. children

Exposure to fluoride in early childhood has been associated with altered cognition, intelligence, attention, and neurobehavior. Fluoride-related neurodevelopment effects have been shown to vary by sex and very little is known about the mechanistic processes involved. There is limited research on how fluoride exposure impacts the epigenome, potentially leading to changes in DNA methylation of specific genes regulating key developmental processes. In the Cincinnati Childhood Allergy and Air Pollution Study (CCAAPS), urine samples were analyzed using a microdiffusion method to determine childhood urinary fluoride adjusted for specific gravity (CUFsg) concentrations. Whole blood DNA methylation was assessed using the Infinium MethylationEPIC BeadChip 850 k Array. In a cross-sectional analysis, we interrogated epigenome-wide DNA methylation at 775,141 CpG loci across the methylome in relation to CUFsg concentrations in 272 early adolescents at age 12 years. Among all participants, higher concentrations of CUF were associated with differential methylation of one CpG (p < 6 × 10-8) located in the gene body of GBF1 (cg25435255). Among females, higher concentrations of CUFsg were associated with differential methylation of 7 CpGs; only three CpGs were differentially methylated among males with no overlap of significant CpGs observed among females. Secondary analyses revealed several differentially methylated regions (DMRs) and CpG loci mapping to genes with key roles in psychiatric outcomes, social interaction, and cognition, as well as immunologic and metabolic phenotypes. While fluoride exposure may impact the epigenome during early adolescence, the functional consequences of these changes are unclear warranting further investigation.

MethylGenotyper: Accurate Estimation of SNP Genotypes and Genetic Relatedness from DNA Methylation Data

Epigenome-wide association studies (EWAS) are susceptible to widespread confounding caused by population structure and genetic relatedness. Nevertheless, kinship estimation is challenging in EWAS without genotyping data. Here, we proposed MethylGenotyper, a method that for the first time enables accurate genotyping at thousands of single nucleotide polymorphisms (SNPs) directly from commercial DNA methylation microarrays. We modeled the intensities of methylation probes near SNPs with a mixture of three beta distributions corresponding to different genotypes and estimated parameters with an expectation-maximization algorithm. We conducted extensive simulations to demonstrate the performance of the method. When applying MethylGenotyper to the Infinium EPIC array data of 4662 Chinese samples, we obtained genotypes at 4319 SNPs with a concordance rate of 98.26%, enabling the identification of 255 pairs of close relatedness. Furthermore, we showed that MethylGenotyper allows for the estimation of both population structure and cryptic relatedness among 702 Australians of diverse ancestry. We also implemented MethylGenotyper in a publicly available R package (https://github.com/Yi-Jiang/MethylGenotyper) to facilitate future large-scale EWAS.

Genetic and multi-omic resources for Alzheimer disease and related dementia from the Knight Alzheimer Disease Research Center

The Knight-Alzheimer Disease Research Center (Knight-ADRC) at Washington University in St. Louis has pioneered and led worldwide seminal studies that have expanded our clinical, social, pathological, and molecular understanding of Alzheimer Disease. Over more than 40 years, research volunteers have been recruited to participate in cognitive, neuropsychologic, imaging, fluid biomarkers, genomic and multi-omic studies. Tissue and longitudinal data collected to foster, facilitate, and support research on dementia and aging. The Genetics and high throughput -omics core (GHTO) have collected of more than 26,000 biological samples from 6,625 Knight-ADRC participants. Samples available include longitudinal DNA, RNA, non-fasted plasma, cerebrospinal fluid pellets, and peripheral blood mononuclear cells. The GHTO has performed deep molecular profiling (genomic, transcriptomic, epigenomic, proteomic, and metabolomic) from large number of brain (n = 2,117), CSF (n = 2,012) and blood/plasma (n = 8,265) samples with the goal of identifying novel risk and protective variants, identify novel molecular biomarkers and causal and druggable targets. Overall, the resources available at GHTO support the increase of our understanding of Alzheimer Disease.

Blood-Based Epigenetic Age Acceleration and Incident Colorectal Cancer Risk: Findings from a Population-Based Case-Control Study

This study investigates the association between epigenetic age acceleration (EAA) derived from DNA methylation and the risk of incident colorectal cancer (CRC). We utilized data from a random population sample of 9,360 individuals (men and women, aged 45-69) from the HAPIEE Study who had been followed up for 16 years. A nested case-control design yielded 35 incident CRC cases and 354 matched controls. Six baseline epigenetic age (EA) measures (Horvath, Hannum, PhenoAge, Skin and Blood (SB), BLUP, and Elastic Net (EN)) were calculated along with their respective EAAs. After adjustment, the odds ratios (ORs) for CRC risk per decile increase in EAA ranged from 1.20 (95% CI: 1.04-1.39) to 1.44 (95% CI: 1.21-1.76) for the Horvath, Hannum, PhenoAge, and BLUP measures. Conversely, the SB and EN EAA measures showed borderline inverse associations with ORs of 0.86-0.87 (95% CI: 0.76-0.99). Tertile analysis reinforced a positive association between CRC risk and four EAA measures (Horvath, Hannum, PhenoAge, and BLUP) and a modest inverse relationship with EN EAA. Our findings from a prospective population-based-case-control study indicate a direct association between incident CRC and four markers of accelerated baseline epigenetic age. In contrast, two markers showed a negative association or no association. These results warrant further exploration in larger cohorts and may have implications for CRC risk assessment and prevention.

Translational toxicoepigenetic Meta-Analyses identify homologous gene DNA methylation reprogramming following developmental phthalate and lead exposure in mouse and human offspring

Although toxicology uses animal models to represent real-world human health scenarios, a critical translational gap between laboratory-based studies and epidemiology remains. In this study, we aimed to understand the toxicoepigenetic effects on DNA methylation after developmental exposure to two common toxicants, the phthalate di(2-ethylhexyl) phthalate (DEHP) and the metal lead (Pb), using a translational paradigm that selected candidate genes from a mouse study and assessed them in four human birth cohorts. Data from mouse offspring developmentally exposed to DEHP, Pb, or control were used to identify genes with sex-specific sites with differential DNA methylation at postnatal day 21. Associations of human infant DNA methylation in homologous mouse genes with prenatal DEHP or Pb were examined with a meta-analysis. Differential methylation was observed on 6 cytosines (adjusted-p < 0.05) and 90 regions (adjusted-p < 0.001). This translational approach offers a unique method that can detect conserved epigenetic differences that are developmentally susceptible to environmental toxicants.

Epigenetic aging in older people living with HIV in Eswatini: a pilot study of HIV and lifestyle factors and epigenetic aging

BACKGROUND

People living with HIV (PLHIV) on effective antiretroviral therapy are living near-normal lives. Although they are less susceptible to AIDS-related complications, they remain highly vulnerable to non-communicable diseases. In this exploratory study of older PLHIV (OPLHIV) in Eswatini, we investigated whether epigenetic aging (i.e., the residual between regressing epigenetic age on chronological age) was associated with HIV-related parameters, and whether lifestyle factors modified these relationships. We calculated epigenetic aging focusing on the Horvath, Hannum, PhenoAge and GrimAge epigenetic clocks, and a pace of biological aging biomarker (DunedinPACE) among 44 OPLHIV in Eswatini.

RESULTS

Age at HIV diagnosis was associated with Hannum epigenetic age acceleration (EAA) (β-coefficient [95% Confidence Interval]; 0.53 [0.05, 1.00], p = 0.03) and longer duration since HIV diagnosis was associated with slower Hannum EAA (- 0.53 [- 1.00, - 0.05], p = 0.03). The average daily dietary intake of fruits and vegetables was associated with DunedinPACE (0.12 [0.03, 0.22], p = 0.01). The associations of Hannum EAA with the age at HIV diagnosis and duration of time since HIV diagnosis were attenuated when the average daily intake of fruits and vegetables or physical activity were included in our models. Diet and self-perceived quality of life measures modified the relationship between CD4+ T cell counts at participant enrollment and Hannum EAA.

CONCLUSIONS

Epigenetic age is more advanced in OPLHIV in Eswatini in those diagnosed with HIV at an older age and slowed in those who have lived for a longer time with diagnosed HIV. Lifestyle and quality of life factors may differentially affect epigenetic aging in OPLHIV. To our knowledge, this is the first study to assess epigenetic aging in OPLHIV in Eswatini and one of the few in sub-Saharan Africa.

Circulating Leukocyte Subsets Before and After a Breast Cancer Diagnosis and Therapy

Importance

Changes in leukocyte composition often precede chronic disease onset. Patients with a history of breast cancer (hereinafter referred to as breast cancer survivors) are at increased risk for subsequent chronic diseases, but the long-term changes in peripheral leukocyte composition following a breast cancer diagnosis and treatment remain unknown.

Objective

To examine longitudinal changes in peripheral leukocyte composition in women who did and did not develop breast cancer and identify whether differences in breast cancer survivors were associated with specific treatments.

Design, Setting, and Participants

In this prospective cohort study, paired blood samples were collected from 2315 women enrolled in The Sister Study, a US-nationwide prospective cohort study of 50 884 women, at baseline (July 2003 to March 2009) and follow-up (October 2013 to March 2015) home visits, with a mean (SD) follow-up interval of 7.6 (1.4) years. By design, approximately half of the included women had been diagnosed and treated for breast cancer after enrollment and before the second blood draw. A total of 410 women were included in the present study, including 185 breast cancer survivors and 225 who remained free of breast cancer over a comparable follow-up period. Data were analyzed from April 21 to September 9, 2022.

Exposures

Breast cancer status and, among breast cancer survivors, cancer treatment type (chemotherapy, radiotherapy, endocrine therapy, or surgery).

Main Outcomes and Measures

Blood DNA methylation data were generated in 2019 using a genome-wide methylation screening tool and deconvolved to estimate percentages of 12 circulating leukocyte subsets.

Results

Of the 410 women included in the analysis, the mean (SD) age at enrollment was 56 (9) years. Compared with breast cancer-free women, breast cancer survivors had decreased percentages of circulating eosinophils (-0.45% [95% CI, -0.87% to -0.03%]; P = .03), total CD4+ helper T cells (-1.50% [95% CI, -2.56% to -0.44%]; P = .01), and memory B cells (-0.22% [95% CI, -0.34% to -0.09%]; P = .001) and increased percentages of circulating naive B cells (0.46% [95% CI, 0.17%-0.75%]; P = .002). In breast cancer survivor-only analyses, radiotherapy was associated with decreases in total CD4+ T cell levels, whereas chemotherapy was associated with increases in naive B cell levels. Surgery and endocrine therapy were not meaningfully associated with leukocyte changes.

Conclusions and Relevance

In this cohort study of 410 women, breast cancer survivors experienced lasting changes in peripheral leukocyte composition compared with women who remained free of breast cancer. These changes may be related to treatment with chemotherapy or radiotherapy and could influence future chronic disease risk.

Associations between CD70 methylation of T cell DNA and age in adults with systemic lupus erythematosus and population controls: The Michigan Lupus Epidemiology & Surveillance (MILES) Program

BACKGROUND

Environmental factors can influence epigenetic regulation, including DNA methylation, potentially contributing to systemic lupus erythematosus (SLE) development and progression. We compared methylation of the B cell costimulatory CD70 gene, in persons with lupus and controls, and characterized associations with age.

RESULTS

In 297 adults with SLE and 92 controls from the Michigan Lupus Epidemiology and Surveillance (MILES) Cohort, average CD70 methylation of CD4+ T cell DNA across 10 CpG sites based on pyrosequencing of the promoter region was higher for persons with SLE compared to controls, accounting for covariates [β = 2.3, p = 0.011]. Using Infinium MethylationEPIC array data at 18 CD70-annoted loci (CD4+ and CD8+ T cell DNA), sites within the promoter region tended to be hypomethylated in SLE, while those within the gene region were hypermethylated. In SLE but not controls, age was significantly associated with pyrosequencing-based CD70 methylation: for every year increase in age, methylation increased by 0.14 percentage points in SLE, accounting for covariates. Also within SLE, CD70 methylation approached a significantly higher level in Black persons compared to White persons (β = 1.8, p = 0.051).

CONCLUSIONS

We describe altered CD70 methylation patterns in T lymphocyte subsets in adults with SLE relative to controls, and report associations particular to SLE between methylation of this immune-relevant gene and both age and race, possibly a consequence of "weathering" or accelerated aging which may have implications for SLE pathogenesis and potential intervention strategies.

DNA methylation profiles reveal sex-specific associations between gestational exposure to ambient air pollution and placenta cell-type composition in the PRISM cohort study

BACKGROUND

Gestational exposure to ambient air pollution has been associated with adverse health outcomes for mothers and newborns. The placenta is a central regulator of the in utero environment that orchestrates development and postnatal life via fetal programming. Ambient air pollution contaminants can reach the placenta and have been shown to alter bulk placental tissue DNA methylation patterns. Yet the effect of air pollution on placental cell-type composition has not been examined. We aimed to investigate whether the exposure to ambient air pollution during gestation is associated with placental cell types inferred from DNA methylation profiles.

METHODS

We leveraged data from 226 mother-infant pairs in the Programming of Intergenerational Stress Mechanisms (PRISM) longitudinal cohort in the Northeastern US. Daily concentrations of fine particulate matter (PM2.5) at 1 km spatial resolution were estimated from a spatiotemporal model developed with satellite data and linked to womens' addresses during pregnancy and infants' date of birth. The proportions of six cell types [syncytiotrophoblasts, trophoblasts, stromal, endothelial, Hofbauer and nucleated red blood cells (nRBCs)] were derived from placental tissue 450K DNA methylation array. We applied compositional regression to examine overall changes in placenta cell-type composition related to PM2.5 average by pregnancy trimester. We also investigated the association between PM2.5 and individual cell types using beta regression. All analyses were performed in the overall sample and stratified by infant sex adjusted for covariates.

RESULTS

In male infants, first trimester (T1) PM2.5 was associated with changes in placental cell composition (p = 0.03), driven by a decrease [per one PM2.5 interquartile range (IQR)] of 0.037 in the syncytiotrophoblasts proportion (95% confidence interval (CI) [- 0.066, - 0.012]), accompanied by an increase in trophoblasts of 0.033 (95% CI: [0.009, 0.064]). In females, second and third trimester PM2.5 were associated with overall changes in placental cell-type composition (T2: p = 0.040; T3: p = 0.049), with a decrease in the nRBC proportion. Individual cell-type analysis with beta regression showed similar results with an additional association found for third trimester PM2.5 and stromal cells in females (decrease of 0.054, p = 0.024).

CONCLUSION

Gestational exposure to air pollution was associated with placenta cell composition. Further research is needed to corroborate these findings and evaluate their role in PM2.5-related impact in the placenta and consequent fetal programming.

Associations Between Epigenetic Age Acceleration and microRNA Expression Among U.S. Firefighters

Epigenetic changes may be biomarkers of health. Epigenetic age acceleration (EAA), the discrepancy between epigenetic age measured via epigenetic clocks and chronological age, is associated with morbidity and mortality. However, the intersection of epigenetic clocks with microRNAs (miRNAs) and corresponding miRNA-based health implications have not been evaluated. We analyzed DNA methylation and miRNA profiles from blood sampled among 332 individuals enrolled across 2 U.S.-based firefighter occupational studies (2015-2018 and 2018-2020). We considered 7 measures of EAA in leukocytes (PhenoAge, GrimAge, Horvath, skin-blood, and Hannum epigenetic clocks, and extrinsic and intrinsic epigenetic age acceleration). We identified miRNAs associated with EAA using individual linear regression models, adjusted for sex, race/ethnicity, chronological age, and cell type estimates, and investigated downstream effects of associated miRNAs with miRNA enrichment analyses and genomic annotations. On average, participants were 38 years old, 88% male, and 75% non-Hispanic white. We identified 183 of 798 miRNAs associated with EAA (FDR q < 0.05); 126 with PhenoAge, 59 with GrimAge, 1 with Horvath, and 1 with the skin-blood clock. Among miRNAs associated with Horvath and GrimAge, there were 61 significantly enriched disease annotations including age-related metabolic and cardiovascular conditions and several cancers. Enriched pathways included those related to proteins and protein modification. We identified miRNAs associated with EAA of multiple epigenetic clocks. PhenoAge had more associations with individual miRNAs, but GrimAge and Horvath had greater implications for miRNA-associated pathways. Understanding the relationship between these epigenetic markers could contribute to our understanding of the molecular underpinnings of aging and aging-related diseases.

Changes in methylation-based aging in women who do and do not develop breast cancer

BACKGROUND

Breast cancer survivors have increased incidence of age-related diseases, suggesting that some survivors may experience faster biological aging.

METHODS

Among 417 women enrolled in the prospective Sister Study cohort, DNA methylation data were generated on paired blood samples collected an average of 7.7 years apart and used to calculate 3 epigenetic metrics of biological aging (PhenoAgeAccel, GrimAgeAccel, and Dunedin Pace of Aging Calculated from the Epigenome [DunedinPACE]). Approximately half (n = 190) the women sampled were diagnosed and treated for breast cancer between blood draws, whereas the other half (n = 227) remained breast cancer-free. Breast tumor characteristics and treatment information were abstracted from medical records.

RESULTS

Among women who developed breast cancer, diagnoses occurred an average of 3.5 years after the initial blood draw and 4 years before the second draw. After accounting for covariates and biological aging metrics measured at baseline, women diagnosed and treated for breast cancer had higher biological aging at the second blood draw than women who remained cancer-free as measured by PhenoAgeAccel (standardized mean difference [β] = 0.13, 95% confidence interval [CI) = 0.00 to 0.26), GrimAgeAccel (β = 0.14, 95% CI = 0.03 to 0.25), and DunedinPACE (β = 0.37, 95% CI = 0.24 to 0.50). In case-only analyses assessing associations with different breast cancer therapies, radiation had strong positive associations with biological aging (PhenoAgeAccel: β = 0.39, 95% CI = 0.19 to 0.59; GrimAgeAccel: β = 0.29, 95% CI = 0.10 to 0.47; DunedinPACE: β = 0.25, 95% CI = 0.02 to 0.48).

CONCLUSIONS

Biological aging is accelerated following a breast cancer diagnosis and treatment. Breast cancer treatment modalities appear to differentially contribute to biological aging.

Epigenetic aging in older people living with HIV in Eswatini: a pilot study of HIV and lifestyle factors and epigenetic aging

Background

People living with HIV (PLHIV) on effective antiretroviral therapy (ART) are living near-normal lives. Although they are less susceptible to AIDS-related complications, they remain highly vulnerable to non-communicable diseases (NCD). In this exploratory study of older PLHIV (OPLHIV) in Eswatini, we investigated whether biological aging (i.e., the difference between epigenetic age and chronological age, termed 'epigenetic age acceleration [EAA]') was associated with HIV-related parameters, and whether lifestyle factors modified these relationships. We calculated EAA focusing on the second-generation epigenetic clocks, PhenoAge and GrimAge, and a pace of aging biomarker (DunedinPACE) among 44 OPLHIV in Eswatini.

Results

Among participants, the PhenoAge clock showed older epigenetic age (68 years old [63, 77]) but a younger GrimAge epigenetic age (median=56 years old [interquartile range=50, 61]) compared to the chronological age (59 years old [54, 66]). Participants diagnosed with HIV at an older age showed slower DunedinPACE (β-coefficient [95% Confidence Interval]; -0.02 [-0.04, -0.01], p=0.002) and longer duration since HIV diagnosis was associated with faster DunedinPACE (0.02 [0.01, 0.04], p=0.002). The average daily dietary intake of fruits and vegetables was associated with faster DunedinPACE (0.12 [0.03, 0.22], p=0.01) and modified the relationship between HIV status variables (number of years living with HIV since diagnosis, age at HIV diagnosis, CD4+ T cell counts) and PhenoAge EAA, and DunedinPACE.

Conclusions

Biological age is accelerated in OPLHIV in Eswatini, with those living with HIV for a longer duration at risk for faster biological aging. Lifestyle factors, especially healthier diets, may attenuate biological aging in OPLHIV. To our knowledge, this is the first study to assess biological aging in Eswatini and one of the few in sub-Saharan Africa.

Epigenomic response to albuterol treatment in asthma-relevant airway epithelial cells

BACKGROUND

Albuterol is the first-line asthma medication used in diverse populations. Although DNA methylation (DNAm) is an epigenetic mechanism involved in asthma and bronchodilator drug response (BDR), no study has assessed whether albuterol could induce changes in the airway epithelial methylome. We aimed to characterize albuterol-induced DNAm changes in airway epithelial cells, and assess potential functional consequences and the influence of genetic variation and asthma-related clinical variables.

RESULTS

We followed a discovery and validation study design to characterize albuterol-induced DNAm changes in paired airway epithelial cultures stimulated in vitro with albuterol. In the discovery phase, an epigenome-wide association study using paired nasal epithelial cultures from Puerto Rican children (n = 97) identified 22 CpGs genome-wide associated with repeated-use albuterol treatment (p < 9 × 10-8). Albuterol predominantly induced a hypomethylation effect on CpGs captured by the EPIC array across the genome (probability of hypomethylation: 76%, p value = 3.3 × 10-5). DNAm changes on the CpGs cg23032799 (CREB3L1), cg00483640 (MYLK4-LINC01600), and cg05673431 (KSR1) were validated in nasal epithelia from 10 independent donors (false discovery rate [FDR] < 0.05). The effect on the CpG cg23032799 (CREB3L1) was cross-tissue validated in bronchial epithelial cells at nominal level (p = 0.030). DNAm changes in these three CpGs were shown to be influenced by three independent genetic variants (FDR < 0.05). In silico analyses showed these polymorphisms regulated gene expression of nearby genes in lungs and/or fibroblasts including KSR1 and LINC01600 (6.30 × 10-14 ≤ p ≤ 6.60 × 10-5). Additionally, hypomethylation at the CpGs cg10290200 (FLNC) and cg05673431 (KSR1) was associated with increased gene expression of the genes where they are located (FDR < 0.05). Furthermore, while the epigenetic effect of albuterol was independent of the asthma status, severity, and use of medication, BDR was nominally associated with the effect on the CpG cg23032799 (CREB3L1) (p = 0.004). Gene-set enrichment analyses revealed that epigenomic modifications of albuterol could participate in asthma-relevant processes (e.g., IL-2, TNF-α, and NF-κB signaling pathways). Finally, nine differentially methylated regions were associated with albuterol treatment, including CREB3L1, MYLK4, and KSR1 (adjusted p value < 0.05).

CONCLUSIONS

This study revealed evidence of epigenetic modifications induced by albuterol in the mucociliary airway epithelium. The epigenomic response induced by albuterol might have potential clinical implications by affecting biological pathways relevant to asthma.

Associations between blood leukocyte DNA methylation and sustained attention in mid-to-late childhood

Aims: To identify associations between DNA methylation (DNAm) across the epigenome and symptoms related to attention-deficit/hyperactivity disorder in a population of Hispanic children. Materials & methods: Among 517 participants in the ELEMENT study aged 9-18 years, we conducted an epigenome-wide association study examining associations between blood leukocyte DNAm and performance on the Conners' continuous performance test (CPT3). Results: DNAm at loci in or near ZNF814, ELF4 and OR6K6 and functional enrichment for gene pathways pertaining to ferroptosis, inflammation, immune response and neurotransmission were significantly related to CPT3 scores. Conclusion: DNAm was associated with CPT3 performance. Further analysis is warranted to understand how these genes and enriched pathways contribute to attention-deficit/hyperactivity disorder.

In utero exposure to diethylstilbestrol and blood DNA methylation in adult women: Results from a meta-analysis of two cohort studies

BACKGROUND

Prenatal exposure to diethylstilbestrol (DES) is associated with several adverse health outcomes. Animal studies have shown associations between prenatal DES exposure and DNA methylation.

OBJECTIVE

The aim of this study was to explore blood DNA methylation in women exposed and unexposed to DES in utero.

METHODS

Sixty women (40 exposed and 20 unexposed) in the National Cancer Institute's Combined DES Cohort Study and 199 women (99 exposed and 100 unexposed women) in the Sister Study Cohort were included in this analysis. Within each study, robust linear regression models were used to assess associations between DES exposure and blood DNA methylation. Study-specific associations were combined using fixed-effect meta-analysis with inverse variance weights. Our analysis focused on CpG sites located within nine candidate genes identified in animal models. We further explored whether in utero DES exposure was associated with age acceleration.

RESULTS

Blood DNA methylation levels at 10 CpG sites in six of the nine candidate genes were statistically significantly associated with prenatal DES exposure (P < 0.05) in this meta-analysis. Genes included EGF, EMB, EGFR, WNT11, FOS, and TGFB1, which are related to cell proliferation and differentiation. The most statistically significant CpG site was cg19830739 in gene EGF, and it was associated with lower methylation levels in women prenatally exposed to DES compared with those not exposed (P < 0.0001; false discovery rate<0.05). The association between prenatal DES exposure in utero and age acceleration was not statistically significant (P = 0.07 for meta-analyzed results).

CONCLUSIONS

There are few opportunities to investigate the effects of prenatal DES exposure. These findings suggest that in utero DES exposure may be associated with differential blood DNA methylation levels, which could mediate the increased risk of several adverse health outcomes observed in exposed women. Our findings need further evaluation using larger data sets.

Long-term exposure to ambient fine particulate components and leukocyte epigenome-wide DNA Methylation in older men: the Normative Aging Study

BACKGROUND

Epigenome-wide association studies of ambient fine particulate matter (PM2.5) have been reported. However, few have examined PM2.5 components (PMCs) and sources or included repeated measures. The lack of high-resolution exposure measurements is the key limitation. We hypothesized that significant changes in DNA methylation might vary by PMCs and the sources.

METHODS

We predicted the annual average of 14 PMCs using novel high-resolution exposure models across the contiguous U.S., between 2000-2018. The resolution was 50 m × 50 m in the Greater Boston Area. We also identified PM2.5 sources using positive matrix factorization. We repeatedly collected blood samples and measured leukocyte DNAm with the Illumina HumanMethylation450K BeadChip in the Normative Aging Study. We then used median regression with subject-specific intercepts to estimate the associations between long-term (one-year) exposure to PMCs / PM2.5 sources and DNA methylation at individual cytosine-phosphate-guanine CpG sites. Significant probes were identified by the number of independent degrees of freedom approach, using the number of principal components explaining > 95% of the variation of the DNA methylation data. We also performed regional and pathway analyses to identify significant regions and pathways.

RESULTS

We included 669 men with 1,178 visits between 2000-2013. The subjects had a mean age of 75 years. The identified probes, regions, and pathways varied by PMCs and their sources. For example, iron was associated with 6 probes and 6 regions, whereas nitrate was associated with 15 probes and 3 regions. The identified pathways from biomass burning, coal burning, and heavy fuel oil combustion sources were associated with cancer, inflammation, and cardiovascular diseases, whereas there were no pathways associated with all traffic.

CONCLUSIONS

Our findings showed that the effects of PM2.5 on DNAm varied by its PMCs and sources.

Genetic impacts on DNA methylation help elucidate regulatory genomic processes

BACKGROUND

Pinpointing genetic impacts on DNA methylation can improve our understanding of pathways that underlie gene regulation and disease risk.

RESULTS

We report heritability and methylation quantitative trait locus (meQTL) analysis at 724,499 CpGs profiled with the Illumina Infinium MethylationEPIC array in 2358 blood samples from three UK cohorts. Methylation levels at 34.2% of CpGs are affected by SNPs, and 98% of effects are cis-acting or within 1 Mbp of the tested CpG. Our results are consistent with meQTL analyses based on the former Illumina Infinium HumanMethylation450 array. Both SNPs and CpGs with meQTLs are overrepresented in enhancers, which have improved coverage on this platform compared to previous approaches. Co-localisation analyses across genetic effects on DNA methylation and 56 human traits identify 1520 co-localisations across 1325 unique CpGs and 34 phenotypes, including in disease-relevant genes, such as USP1 and DOCK7 (total cholesterol levels), and ICOSLG (inflammatory bowel disease). Enrichment analysis of meQTLs and integration with expression QTLs give insights into mechanisms underlying cis-meQTLs (e.g. through disruption of transcription factor binding sites for CTCF and SMC3) and trans-meQTLs (e.g. through regulating the expression of ACD and SENP7 which can modulate DNA methylation at distal sites).

CONCLUSIONS

Our findings improve the characterisation of the mechanisms underlying DNA methylation variability and are informative for prioritisation of GWAS variants for functional follow-ups. The MeQTL EPIC Database and viewer are available online at https://epicmeqtl.kcl.ac.uk .

No impact of prenatal paracetamol and folic acid exposure on cord blood DNA methylation in children with attention-deficit/hyperactivity disorder

Pharmacoepigenetic studies are important to understand the mechanisms through which medications influence the developing fetus. For instance, we and others have reported associations between prenatal paracetamol exposure and offspring DNA methylation (DNAm). Additionally, folic acid (FA) intake during pregnancy has been associated with DNAm in genes linked to developmental abnormalities. In this study, we aimed to: (i) expand on our previous findings showing differential DNAm associated with long-term prenatal paracetamol exposure in offspring with attention-deficit/hyperactivity disorder (ADHD), and (ii) examine if there is an interaction effect of FA and paracetamol on DNAm in children with ADHD. We used data from the Norwegian Mother, Father and Child Cohort Study (MoBa) and the Medical Birth Registry of Norway (MBRN). We did not identify any impact of paracetamol or any interaction effect of paracetamol and FA on cord blood DNAm in children with ADHD. Our results contribute to the growing literature on prenatal pharmacoepigenetics, but should be replicated in other cohorts. Replication of pharmacoepigenetic studies is essential to ensure robust findings and to increase the clinical relevance of such studies.

Next-Generation Sequencing and Triple-Negative Breast Cancer: Insights and Applications

The poor survival of triple-negative breast cancer (TNBC) is due to its aggressive behavior, large heterogeneity, and high risk of recurrence. A comprehensive molecular investigation of this type of breast cancer using high-throughput next-generation sequencing (NGS) methods may help to elucidate its potential progression and discover biomarkers related to patient survival. In this review, the NGS applications in TNBC research are described. Many NGS studies point to TP53 mutations, immunocheckpoint response genes, and aberrations in the PIK3CA and DNA repair pathways as recurrent pathogenic alterations in TNBC. Beyond their diagnostic and predictive/prognostic value, these findings suggest potential personalized treatments in PD -L1-positive TNBC or in TNBC with a homologous recombination deficit. Moreover, the comprehensive sequencing of large genomes with NGS has enabled the identification of novel markers with clinical value in TNBC, such as AURKA, MYC, and JARID2 mutations. In addition, NGS investigations to explore ethnicity-specific alterations have pointed to EZH2 overexpression, BRCA1 alterations, and a BRCA2-delaAAGA mutation as possible molecular signatures of African and African American TNBC. Finally, the development of long-read sequencing methods and their combination with optimized short-read techniques promise to improve the efficiency of NGS approaches for future massive clinical use.

Genome-wide DNA methylation association study of recent and cumulative marijuana use in middle aged adults

Marijuana is a widely used psychoactive substance in the US and medical and recreational legalization has risen over the past decade. Despite the growing number of individuals using marijuana, studies investigating the association between epigenetic factors and recent and cumulative marijuana use remain limited. We therefore investigated the association between recent and cumulative marijuana use and DNA methylation levels. Participants from the Coronary Artery Risk Development in Young Adults Study with whole blood collected at examination years (Y) 15 and Y20 were randomly selected to undergo DNA methylation profiling at both timepoints using the Illumina MethylationEPIC BeadChip. Recent use of marijuana was queried at each examination and used to estimate cumulative marijuana use from Y0 to Y15 and Y20. At Y15 (n = 1023), we observed 22 and 31 methylation markers associated (FDR P ≤ 0.05) with recent and cumulative marijuana use and 132 and 16 methylation markers at Y20 (n = 883), respectively. We replicated 8 previously reported methylation markers associated with marijuana use. We further identified 640 cis-meQTLs and 198 DMRs associated with recent and cumulative use at Y15 and Y20. Differentially methylated genes were statistically overrepresented in pathways relating to cellular proliferation, hormone signaling, and infections as well as schizophrenia, bipolar disorder, and substance-related disorders. We identified numerous methylation markers, pathways, and diseases associated with recent and cumulative marijuana use in middle-aged adults, providing additional insight into the association between marijuana use and the epigenome. These results provide novel insights into the role marijuana has on the epigenome and related health conditions.

Methylation-Based Biological Age and Hypertension Prevalence and Incidence

BACKGROUND

Hypertension is common in older individuals and is a major risk factor for cardiovascular disease. Blood DNA methylation profiles have been used to derive metrics of biological age that capture age-related physiological change, disease risk, and mortality. The relationships between hypertension and DNA methylation-based biological age metrics have yet to be carefully described.

METHODS

Among 4419 women enrolled in the prospective Sister Study cohort, DNA methylation data generated from whole blood samples collected at baseline were used to calculate 3 biological age metrics (PhenoAgeAccel, GrimAgeAccel, DunedinPACE). Women were classified as hypertensive at baseline if they had high blood pressure (systolic blood pressure ≥140 mm Hg or diastolic blood pressure ≥90 mm Hg) or reported current use of antihypertensive medication. New incident cases of hypertension during follow-up were identified via self-report on annual health questionnaires.

RESULTS

All 3 DNA methylation metrics of biological age were positively associated with prevalent hypertension at baseline (per 1-SD increase; PhenoAgeAccel, adjusted odds ratio, 1.16 [95% CI, 1.05-1.28]; GrimAgeAccel, adjusted odds ratio, 1.28 [95% CI, 1.14-1.45]; DunedinPACE, adjusted odds ratio, 1.16 [95% CI, 1.03-1.30]). Among 2610 women who were normotensive at baseline, women with higher biological age were more likely to be diagnosed with incident hypertension (per 1-SD increase; PhenoAgeAccel, adjusted hazard ratio, 1.09 [95% CI, 0.97-1.23]; GrimAgeAccel, adjusted hazard ratio, 1.16 [95% CI, 0.99-1.36]; DunedinPACE, adjusted hazard ratio, 1.16 [95% CI, 1.01-1.33]).

CONCLUSIONS

Methylation-based biological age metrics increase before a hypertension diagnosis and appear to remain elevated in the years after clinical diagnosis and treatment.

Time to relapse in chronic lymphocytic leukemia and DNA-methylation-based biological age

Chronic lymphocytic leukemia (CLL) is a mature B cell neoplasm with a predilection for older individuals. While previous studies have identified epigenetic signatures associated with CLL, whether age-related DNA methylation changes modulate CLL relapse remains elusive. In this study, we examined the association between epigenetic age acceleration and time to CLL relapse in a publicly available dataset. DNA methylation profiling of 35 CLL patients prior to initiating chemoimmunotherapy was performed using the Infinium HumanMethylation450 BeadChip. Four epigenetic age acceleration metrics (intrinsic epigenetic age acceleration [IEAA], extrinsic epigenetic age acceleration [EEAA], PhenoAge acceleration [PhenoAA], and GrimAge acceleration [GrimAA]) were estimated from blood DNA methylation levels. Linear, quantile, and logistic regression and receiver operating characteristic curve analyses were conducted to assess the association between each epigenetic age metric and time to CLL relapse. EEAA (p = 0.011) and PhenoAA (p = 0.046) were negatively and GrimAA (p = 0.040) was positively associated with time to CLL relapse. Simultaneous assessment of EEAA and GrimAA in male patients distinguished patients who relapsed early from patients who relapsed later (p = 0.039). No associations were observed with IEAA. These findings suggest epigenetic age acceleration prior to chemoimmunotherapy initiation is associated with time to CLL relapse. Our results provide novel insight into the association between age-related DNA methylation changes and CLL relapse and may serve has biomarkers for treatment relapse, and potentially, treatment selection.

Prenatal phthalate exposure and cord blood DNA methylation

Exposure to phthalates has been shown to impede the human endocrine system, resulting in deleterious effects on pregnant women and their children. Phthalates modify DNA methylation patterns in infant cord blood. We examined the association between prenatal phthalate exposure and DNA methylation patterns in cord blood in a Korean birth cohort. Phthalate levels were measured in 274 maternal urine samples obtained during late pregnancy and 102 neonatal urine samples obtained at birth, and DNA methylation levels were measured in cord blood samples. For each infant in the cohort, associations between CpG methylation and both maternal and neonate phthalate levels were analyzed using linear mixed models. The results were combined with those from a meta-analysis of the levels of phthalates in maternal and neonatal urine samples, which were also analyzed for MEOHP, MEHHP, MnBP, and DEHP. This meta-analysis revealed significant associations between the methylation levels of CpG sites near the CHN2 and CUL3 genes, which were also associated with MEOHP and MnBP in neonatal urine. When the data were stratified by the sex of the infant, MnBP concentration was found to be associated with one CpG site near the OR2A2 and MEGF11 genes in female infants. In contrast, the concentrations of the three maternal phthalates showed no significant association with CpG site methylation. Furthermore, the data identified distinct differentially methylated regions in maternal and neonatal urine samples following exposure to phthalates. The CpGs with methylation levels that were positively associated with phthalate levels (particularly MEOHP and MnBP) were found to be enriched genes and related pathways. These results indicate that prenatal phthalate exposure is significantly associated with DNA methylation at multiple CpG sites. These alterations in DNA methylation may serve as biomarkers of maternal exposure to phthalates in infants and are potential candidates for investigating the mechanisms by which phthalates impact maternal and neonatal health.

Accelerated epigenetic aging and myopenia in young adult cancer survivors

BACKGROUND

Young adult cancer survivors experience early aging-related morbidities and mortality. Biological aging biomarkers may identify at-risk survivors and increase our understanding of mechanisms underlying this accelerated aging.

METHODS

Using an observational study design, we cross-sectionally measured DNA methylation-based epigenetic age in young adult cancer survivors at a tertiary, academic state cancer hospital. Participants were a convenience sample of consecutively enrolled survivors of childhood, adolescent, and young adult cancers treated with either an anthracycline or alkylating agent, and who were at least 3 months post-treatment. Similarly aged healthy comparators were consecutively enrolled. Cancer treatment and treatment intensity were compared to DNA methylation-based epigenetic age and pace of aging.

RESULTS

Sixty survivors (58 completing assessments, mean age 20.5 years, range 18-29) and 27 comparators (mean age 20 years, range 17-29) underwent DNA methylation measurement. Survivors were predominantly female (62%) and white (60%) and averaged nearly 6 years post-treatment (range 0.2-25 years). Both epigenetic age (AgeAccelGrim: 1.5 vs. -2.4, p < 0.0001; AgeAccelPheno 2.3 vs. -3.8, p = 0.0013) and pace of aging (DunedinPACE 0.99 vs. 0.83, p < 0.0001) were greater in survivors versus comparators. In case-case adjusted analysis, compared to survivors with normal muscle mass, myopenic survivors had higher AgeAccelGrim (2.2 years, 95% CI 0.02-4.33, p = 0.02), AgeAccelPheno (6.2 years, 2.36-10.09, p < 0.001), and DunedinPACE (0.11, 0.05-0.17, p < 0.001).

CONCLUSIONS

Epigenetic age is older and pace of aging is faster in young adult cancer survivors compared to noncancer peers, which is evident in the early post-therapy period. Survivors with physiological impairment demonstrate greater epigenetic age advancement. Measures of epigenetic age may identify young adult survivors at higher risk for poor functional and health outcomes.

Mediation effects of DNA methylation and hydroxymethylation on birth outcomes after prenatal per- and polyfluoroalkyl substances (PFAS) exposure in the Michigan mother-infant Pairs cohort

BACKGROUND

Per- and polyfluoroalkyl substances (PFAS) are chemicals that are resistant to degradation and ubiquitous in our environments. PFAS may impact the developing epigenome, but current human evidence is limited to assessments of total DNA methylation. We assessed associations between first trimester PFAS exposures with newborn DNA methylation, including 5-methylcytosine (5-mC) and 5-hydroxymethylcytosine (5-hmC). DNA methylation mediation of associations between PFAS and birth outcomes were explored in the Michigan Mother Infant Pairs cohort. Nine PFAS were measured in maternal first trimester blood. Seven were highly detected and included for analysis: PFHxS, PFOA, PFOS, PFNA, PFDA, PFUnDA, and MeFOSAA. Bisulfite-converted cord blood DNA (n = 141) and oxidative-bisulfite-converted cord blood (n = 70) were assayed on Illumina MethylationEPIC BeadChips to measure total DNA methylation (5-mC + 5-hmC) and 5-mC/5-hmC. Correcting for multiple comparisons, beta regressions were used to assess associations between levels of PFAS and total methylation, 5-mC, or 5-hmC. Nonlinear mediation analyses were used to assess the epigenetic meditation effect between PFAS and birth outcomes.

RESULTS

PFAS was significantly associated with total methylation (q < 0.05: PFHxS-12 sites; PFOS-19 sites; PFOA-2 sites; PFNA-3 sites; PFDA-4 sites). In 72 female infants and 69 male infants, there were sex-specific associations between five PFAS and DNA methylation. 5-mC and 5-hmC were each significantly associated with thousands of sites for PFHxS, PFOS, PFNA, PFDA, PFUnDA, and MeFOSAA (q < 0.05). Clusters of 5-mC and 5-hmC sites were significant mediators between PFNA and PFUnDA and decreased gestational age (q < 0.05).

CONCLUSIONS

This study demonstrates the mediation role of specific types of DNA methylation on the relationship between PFAS exposure and birth outcomes. These results suggest that 5-mC and 5-hmC may be more sensitive to the developmental impacts of PFAS than total DNA methylation.

Epigenome-Wide Association Studies of the Fractional Exhaled Nitric Oxide and Bronchodilator Drug Response in Moderate-to-Severe Pediatric Asthma

Asthma is the most prevalent pediatric chronic disease. Bronchodilator drug response (BDR) and fractional exhaled nitric oxide (FeNO) are clinical biomarkers of asthma. Although DNA methylation (DNAm) contributes to asthma pathogenesis, the influence of DNAm on BDR and FeNO is scarcely investigated. This study aims to identify DNAm markers in whole blood associated either with BDR or FeNO in pediatric asthma. We analyzed 121 samples from children with moderate-to-severe asthma. The association of genome-wide DNAm with BDR and FeNO has been assessed using regression models, adjusting for age, sex, ancestry, and tissue heterogeneity. Cross-tissue validation was assessed in 50 nasal samples. Differentially methylated regions (DMRs) and enrichment in traits and biological pathways were assessed. A false discovery rate (FDR) < 0.1 and a genome-wide significance threshold of p < 9 × 10−8 were used to control for false-positive results. The CpG cg12835256 (PLA2G12A) was genome-wide associated with FeNO in blood samples (coefficient= −0.015, p = 2.53 × 10−9) and nominally associated in nasal samples (coefficient = −0.015, p = 0.045). Additionally, three CpGs were suggestively associated with BDR (FDR < 0.1). We identified 12 and four DMRs associated with FeNO and BDR (FDR < 0.05), respectively. An enrichment in allergic and inflammatory processes, smoking, and aging was observed. We reported novel associations of DNAm markers associated with BDR and FeNO enriched in asthma-related processes.

Epigenome-wide DNA methylation in leukocytes and toenail metals: The normative aging study

BACKGROUND

Environmental metal exposures have been associated with multiple deleterious health endpoints. DNA methylation (DNAm) may provide insight into the mechanisms underlying these relationships. Toenail metals are non-invasive biomarkers, reflecting a medium-term time exposure window.

OBJECTIVES

This study examined variation in leukocyte DNAm and toenail arsenic (As), cadmium (Cd), lead (Pb), manganese (Mn), and mercury (Hg) among elderly men in the Normative Aging Study, a longitudinal cohort.

METHODS

We repeatedly collected samples of blood and toenail clippings. We measured DNAm in leukocytes with the Illumina HumanMethylation450 K BeadChip. We first performed median regression to evaluate the effects of each individual toenail metal on DNAm at three levels: individual cytosine-phosphate-guanine (CpG) sites, regions, and pathways. Then, we applied a Bayesian kernel machine regression (BKMR) to assess the joint and individual effects of metal mixtures on DNAm. Significant CpGs were identified using a multiple testing correction based on the independent degrees of freedom approach for correlated outcomes. The approach considers the effective degrees of freedom in the DNAm data using the principal components that explain >95% variation of the data.

RESULTS

We included 564 subjects (754 visits) between 1999 and 2013. The numbers of significantly differentially methylated CpG sites, regions, and pathways varied by metals. For example, we found six significant pathways for As, three for Cd, and one for Mn. The As-associated pathways were associated with cancer (e.g., skin cancer) and cardiovascular disease, whereas the Cd-associated pathways were related to lung cancer. Metal mixtures were also associated with 47 significant CpG sites, as well as pathways, mainly related to cancer and cardiovascular disease.

CONCLUSIONS

This study provides an approach to understanding the potential epigenetic mechanisms underlying observed relations between toenail metals and adverse health endpoints.

Alcohol consumption and epigenetic age acceleration in young adults

Alcohol is a widely consumed substance in the United States, however its effect on aging remains understudied. In this study of young adults, we examined whether cumulative alcohol consumption, i.e., alcohol years of beer, liquor, wine, and total alcohol, and recent binge drinking, were associated with four measures of age-related epigenetic changes via blood DNA methylation. A random subset of study participants in the Coronary Artery Risk Development in Young Adults Study underwent DNA methylation profiling using the Illumina MethylationEPIC Beadchip. Participants with alcohol consumption and methylation data at examination years 15 (n = 1,030) and 20 (n = 945) were included. Liquor and total alcohol consumption were associated with a 0.31-year (P = 0.002) and a 0.12-year (P = 0.013) greater GrimAge acceleration (GAA) per additional five alcohol years, while beer and wine consumption observed marginal (P = 0.075) and no associations (P = 0.359) with GAA, respectively. Any recent binge drinking and the number of days of binge drinking were associated with a 1.38-year (P < 0.001) and a 0.15-year (P < 0.001) higher GAA, respectively. We observed statistical interactions between cumulative beer (P < 0.001) and total alcohol (P = 0.004) consumption with chronological age, with younger participants exhibiting a higher average in GAA compared to older participants. No associations were observed with the other measures of epigenetic aging. These results suggest cumulative liquor and total alcohol consumption and recent binge drinking may alter age-related epigenetic changes as captured by GAA. With the increasing aging population and widespread consumption of alcohol, these findings may have potential implications for lifestyle modification to promote healthy aging.

Peripheral Immune Cell Composition is Altered in Women Before and After a Hypertension Diagnosis

BACKGROUND

The development and consequences of hypertension involve multiple biological systems that may include changes in immune profiles. Whether hypertension is related to peripheral immune cell composition has not been examined in large human cohorts.

METHODS

We estimated circulating proportions of 12 leukocyte subsets from the lymphoid and myeloid lineages by deconvolving cell-type-specific DNA methylation data from 4124 women. Hypertension status at baseline was defined by current use of antihypertensive medication and blood pressure measurements while new incident cases were identified during follow-up via annual health questionnaires.

RESULTS

Among hypertension-free women at baseline, higher B cell and lower naïve CD4+ helper T cell proportions were associated with subsequent increased hazard of hypertension incidence (B cells; adjusted HR: 1.17 [95% CI: 1.02-1.35]; P=0.03; naïve CD4+ T cell, adjusted HR: 0.88 [95% CI: 0.78-0.99]; P=0.04). Blood pressure measurements at baseline were similarly positively associated with B cells and inversely associated with naïve CD4+ helper T cells. Compared to normotensive women, women with hypertension had higher circulating proportions of neutrophils (adjusted odds ratio: 1.18 [95% CI: 1.07-1.31]; P=0.001) and lower proportions of CD4+ helper T cells (adjusted odds ratio: 0.90 [95% CI: 0.81-1.00] P=0.05), natural killers (adjusted odds ratio: 0.82 [95% CI: 0.74-0.91]; P<0.001), and B cells (adjusted odds ratio: 0.84 [95% CI: 0.74-0.96]; P=0.01).

CONCLUSIONS

These observations suggest that shifts in lymphocyte subsets occur before hypertension development, followed by later changes to neutrophils and additional lymphocytes.

Stable methylation loci are associated with systolic blood pressure in a Croatian island population

Background: The objective was to identify stable and dynamic DNA methylation loci associated with cardiometabolic traits among an adult population from the Croatian island of Hvar. Materials & methods: An epigenome-wide association study was conducted using peripheral blood longitudinally collected at two time points 10 years apart via Infinium MethylationEPIC beadarray (n = 112). Stable and dynamic loci were identified using linear mixed models. Associations between cardiometabolic traits and loci were assessed using linear models. Results: 22 CpG loci were significantly associated with systolic blood pressure. Twenty were stable and two were dynamic. Conclusion: Multiple genes may be involved in the determination of systolic blood pressure level via stable epigenetic programming, potentially established earlier in life.

Marijuana use and DNA methylation-based biological age in young adults

BACKGROUND

Marijuana is the third most commonly used drug in the USA and efforts to legalize it for medical and recreational use are growing. Despite the increase in use, marijuana's effect on aging remains understudied and understanding the effects of marijuana on molecular aging may provide novel insights into the role of marijuana in the aging process. We therefore sought to investigate the association between cumulative and recent use of marijuana with epigenetic age acceleration (EAA) as estimated from blood DNA methylation.

RESULTS

A random subset of participants from The Coronary Artery Risk Development in Young Adults (CARDIA) Study with available whole blood at examination years (Y) 15 and Y20 underwent epigenomic profiling. Four EAA estimates (intrinsic epigenetic age acceleration, extrinsic epigenetic age acceleration, PhenoAge acceleration, and GrimAge acceleration) were calculated from DNA methylation levels measured at Y15 and Y20. Ever use and cumulative marijuana-years were calculated from the baseline visit to Y15 and Y20, and recent marijuana use (both any and number of days of use in the last 30 days) were calculated at Y15 and Y20. Ever use of marijuana and each additional marijuana-year were associated with a 6-month (P < 0.001) and a 2.5-month (P < 0.001) higher average in GrimAge acceleration (GAA) using generalized estimating equations, respectively. Recent use and each additional day of recent use were associated with a 20-month (P < 0.001) and a 1-month (P < 0.001) higher GAA, respectively. A statistical interaction between marijuana-years and alcohol consumption on GAA was observed (P = 0.011), with nondrinkers exhibiting a higher GAA (β = 0.21 [95% CI 0.05, 0.36], P = 0.008) compared to heavy drinkers (β = 0.05 [95% CI - 0.09, 0.18], P = 0.500) per each additional marijuana-year. No associations were observed for the remaining EAA estimates.

CONCLUSIONS

These findings suggest cumulative and recent marijuana use are associated with age-related epigenetic changes that are related to lifespan. These observed associations may be modified by alcohol consumption. Given the increase in use and legalization, these findings provide novel insight on the effect of marijuana use on the aging process as captured through blood DNA methylation.

Significant variation in the performance of DNA methylation predictors across data preprocessing and normalization strategies

BACKGROUND

DNA methylation (DNAm)-based predictors hold great promise to serve as clinical tools for health interventions and disease management. While these algorithms often have high prediction accuracy, the consistency of their performance remains to be determined. We therefore conduct a systematic evaluation across 101 different DNAm data preprocessing and normalization strategies and assess how each analytical strategy affects the consistency of 41 DNAm-based predictors.

RESULTS

Our analyses are conducted in a large EPIC DNAm array dataset from the Jackson Heart Study (N = 2053) that included 146 pairs of technical replicate samples. By estimating the average absolute agreement between replicate pairs, we show that 32 out of 41 predictors (78%) demonstrate excellent consistency when appropriate data processing and normalization steps are implemented. Across all pairs of predictors, we find a moderate correlation in performance across analytical strategies (mean rho = 0.40, SD = 0.27), highlighting significant heterogeneity in performance across algorithms. Successful or unsuccessful removal of technical variation furthermore significantly impacts downstream phenotypic association analysis, such as all-cause mortality risk associations.

CONCLUSIONS

We show that DNAm-based algorithms are sensitive to technical variation. The right choice of data processing strategy is important to achieve reproducible estimates and improve prediction accuracy in downstream phenotypic association analyses. For each of the 41 DNAm predictors, we report its degree of consistency and provide the best performing analytical strategy as a guideline for the research community. As DNAm-based predictors become more and more widely used, our work helps improve their performance and standardize their implementation.

Sex-specific and generational effects of alcohol and tobacco use on epigenetic age acceleration in the Michigan longitudinal study

Background:

Methods:

Results:

Conclusions:

Association between prenatal cadmium exposure and cord blood DNA methylation

Prenatal cadmium exposure is known to affect infant growth and organ development. Nonetheless, the role of DNA methylation in cadmium-related health effects has yet to be determined. To this end, we investigated the relationship between prenatal cadmium exposure and cord blood DNA methylation in Korean infants through an epigenome-wide association study. Cadmium concentrations in maternal blood during early and late pregnancy and in cord blood collected from newborns were measured using atomic adsorption spectrometry and DNA methylation analysis was conducted using HumanMethylationEPIC BeadChip kits. After adjusting for infant sex, maternal pregnancy body mass index, smoking status, and estimated leukocyte composition, we analyzed the association between CpG methylation and cadmium concentration in 364 samples. Among 835,252 CpG sites, maternal blood cadmium concentration in early pregnancy was significantly associated with two differentially methylated CpG sites, cg05537752 and cg24904393, which were annotated ATP9A and no gene, respectively. The study findings indicate that prenatal cadmium exposure is significantly associated with methylation statuses of several CpG sites and regions in Korean infants, especially during early pregnancy.

The shared mother-child epigenetic signature of neglect is related to maternal adverse events

Studies of DNA methylation have revealed the biological mechanisms by which life adversity confers risk for later physical and mental health problems. What remains unknown is the “biologically embedding” of maternal adverse experiences resulting in maladaptive parenting and whether these epigenetic effects are transmitted to the next generation. This study focuses on neglectful mothering indexed by a severe disregard for the basic and psychological needs of the child. Using the Illumina Human Methylation EPIC BeadChip in saliva samples, we identified genes with differentially methylated regions (DMRs) in those mothers with (n = 51), versus those without (n = 87), neglectful behavior that present similar DMRs patterns in their children being neglected versus non-neglected (n = 40 vs. 75). Mothers reported the emotional intensity of adverse life events. After covariate adjustment and multiple testing corrections, we identified 69 DMRs in the mother epigenome and 42 DMRs in the child epigenome that were simultaneously above the α = 0.01 threshold. The common set of nine DMRs contained genes related to childhood adversity, neonatal and infant diabetes, child neurobehavioral development and other health problems such as obesity, hypertension, cancer, posttraumatic stress, and the Alzheimer’s disease; four of the genes were associated with maternal life adversity. Identifying a shared epigenetic signature of neglect linked to maternal life adversity is an essential step in breaking the intergenerational transmission of one of the most common forms of childhood maltreatment.

Nicotine dose-dependent epigenomic-wide DNA methylation changes in the mice with long-term electronic cigarette exposure

Epigenomic-wide DNA methylation profiling holds the potential to reflect both electronic cigarette exposure-associated risks and individual poor health outcomes. However, a systemic study in animals or humans is still lacking. Using the Infinium Mouse Methylation BeadChip, we examined the DNA methylation status of white blood cells in male ApoE-/- mice after 14 weeks of electronic cigarette exposure with the InExpose system (2 hr/day, 5 days/week, 50% PG and 50% VG) with low (6 mg/ml) and high (36 mg/ml) nicotine concentrations. Our results indicate that electronic cigarette aerosol inhalation induces significant alteration of 8,985 CpGs in a dose-dependent manner (FDR<0.05); 7,389 (82.2%) of the CpG sites are annotated with known genes. Among the top 6 significant CpG sites (P-value<1e-8), 4 CpG sites are located in the known genes, and most (3/5) of these genes have been related to cigarette smoking. The other two CpGs are close to/associated with the Phc2 gene that was recently linked to smoking in a transcriptome-wide associations study. Furthermore, the gene set enrichment analysis highlights the activation of MAPK and 4 cardiomyocyte/cardiomyopathy-related signaling pathways (including adrenergic signaling in cardiomyocytes and arrhythmogenic right ventricular cardiomyopathy) following repeated electronic cigarette use. The MAPK pathway activation correlates well with our finding of increased cytokine mRNA expression after electronic cigarette exposure in the same batch of mice. Interestingly, two pathways related to mitochondrial activities, namely mitochondrial gene expression and mitochondrial translation, are also activated after electronic cigarette exposure. Elucidating the relationship between these pathways and the increased circulating mitochondrial DNA observed here will provide further insight into the cell-damaging effects of prolonged inhalation of e-cigarette aerosols.

Association of Cardiovascular Health Through Young Adulthood With Genome-Wide DNA Methylation Patterns in Midlife: The CARDIA Study

BACKGROUND

Cardiovascular health (CVH) from young adulthood is strongly associated with an individual's future risk of cardiovascular disease (CVD) and total mortality. Defining epigenomic biomarkers of lifelong CVH exposure and understanding their roles in CVD development may help develop preventive and therapeutic strategies for CVD.

METHODS

In 1085 CARDIA study (Coronary Artery Risk Development in Young Adults) participants, we defined a clinical cumulative CVH score that combines body mass index, blood pressure, total cholesterol, and fasting glucose measured longitudinally from young adulthood through middle age over 20 years (mean age, 25-45). Blood DNA methylation at >840 000 methylation markers was measured twice over 5 years (mean age, 40 and 45). Epigenome-wide association analyses on the cumulative CVH score were performed in CARDIA and compared in the FHS (Framingham Heart Study). We used penalized regression to build a methylation-based risk score to evaluate the risk of incident coronary artery calcification and clinical CVD events.

RESULTS

We identified 45 methylation markers associated with cumulative CVH at false discovery rate <0.01 (P=4.7E-7-5.8E-17) in CARDIA and replicated in FHS. These associations were more pronounced with methylation measured at an older age. CPT1A, ABCG1, and SREBF1 appeared as the most prominent genes. The 45 methylation markers were mostly located in transcriptionally active chromatin and involved lipid metabolism, insulin secretion, and cytokine production pathways. Three methylation markers located in genes SARS1, SOCS3, and LINC-PINT statistically mediated 20.4% of the total effect between CVH and risk of incident coronary artery calcification. The methylation risk score added information and significantly (P=0.004) improved the discrimination capacity of coronary artery calcification status versus CVH score alone and showed association with risk of incident coronary artery calcification 5 to 10 years later independent of cumulative CVH score (odds ratio, 1.87; P=9.66E-09). The methylation risk score was also associated with incident clinical CVD in FHS (hazard ratio, 1.28; P=1.22E-05).

CONCLUSIONS

Cumulative CVH from young adulthood contributes to midlife epigenetic programming over time. Our findings demonstrate the role of epigenetic markers in response to CVH changes and highlight the potential of epigenomic markers for precision CVD prevention, and earlier detection of subclinical CVD, as well.

Low reliability of DNA methylation across Illumina Infinium platforms in cord blood: implications for replication studies and meta-analyses of prenatal exposures

Background

There is an increasing interest in the role of epigenetics in epidemiology, but the emerging research field faces several critical biological and technical challenges. In particular, recent studies have shown poor correlation of measured DNA methylation (DNAm) levels within and across Illumina Infinium platforms in various tissues. In this study, we have investigated concordance between 450 k and EPIC Infinium platforms in cord blood. We could not replicate our previous findings on the association of prenatal paracetamol exposure with cord blood DNAm, which prompted an investigation of cross-platform DNAm differences.

Results

This study is based on two DNAm data sets from cord blood samples selected from the Norwegian Mother, Father and Child Cohort Study (MoBa). DNAm of one data set was measured using the 450 k platform and the other data set was measured using the EPIC platform. Initial analyses of the EPIC data could not replicate any of our previous significant findings in the 450 k data on associations between prenatal paracetamol exposure and cord blood DNAm. A subset of the samples (n = 17) was included in both data sets, which enabled analyses of technical sources potentially contributing to the negative replication. Analyses of these 17 samples with repeated measurements revealed high per-sample correlations (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\stackrel{\mathrm{-}}{\text{R}}\hspace{0.17em}\approx$$\end{document}R-≈ 0.99), but low per-CpG correlations (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\stackrel{\mathrm{-}}{\text{R}}$$\end{document}R- ≈ 0.24) between the platforms. 1.7% of the CpGs exhibited a mean DNAm difference across platforms > 0.1. Furthermore, only 26.7% of the CpGs exhibited a moderate or better cross-platform reliability (intra-class correlation coefficient ≥ 0.5).

Conclusion

The observations of low cross-platform probe correlation and reliability corroborate previous reports in other tissues. Our study cannot determine the origin of the differences between platforms. Nevertheless, it emulates the setting in studies using data from multiple Infinium platforms, often analysed several years apart. Therefore, the findings may have important implications for future epigenome-wide association studies (EWASs), in replication, meta-analyses and longitudinal studies. Cognisance and transparency of the challenges related to cross-platform studies may enhance the interpretation, replicability and validity of EWAS results both in cord blood and other tissues, ultimately improving the clinical relevance of epigenetic epidemiology.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-022-01299-3.

Prenatal lead exposure, telomere length in cord blood, and DNA methylation age in the PROGRESS prenatal cohort

BACKGROUND

Lead is a ubiquitous pollutant with deleterious effects on human health and remains a major current public health concern in developing countries. This heavy metal may interfere with nucleic acids via oxidative stress or epigenetic changes that affect biological markers of aging, e.g., telomere length and DNA methylation (DNAm). Telomere shortening associates with biological age in newborns, and DNA methylation at specific CpG sites can be used to calculate "epigenetic clocks".

OBJECTIVE

The aim of this study was to examine the associations of prenatal lead exposures with telomere length and DNA-methylation-based predictors of age in cord blood.

DESIGN

The study included 507 mother-child pairs from the Programming Research in Obesity, Growth, Environment and Social Stressors (PROGRESS) study, a birth cohort in Mexico City. Maternal blood (second trimester, third trimester and at delivery) and bone lead levels (one month postpartum) were measured using inductively coupled plasma-mass spectrometry and X-ray fluorescence, respectively. Cord blood leukocyte telomere length was measured using quantitative PCR and apparent age by DNA methylation biomarkers, i.e., Horvath's DNA methylation age and the Knight's predictor of gestational age.

RESULTS

Average maternal age was 28.5 ± 5.5 years, and 51.5% reported low socioeconomic status. Children's mean telomere length was 1.2 ± 1.3 relative units, and mean DNA methylation ages using the Horvath's and Knight's clocks were -2.6 ± 0.1 years and 37.9 ± 1.4 weeks (mean ± SD), respectively. No significant associations were found between maternal blood and bone lead concentrations with telomere length and DNAm age in newborns.

CONCLUSION

We found no associations of prenatal lead exposure with telomere length and DNA methylation age biomarkers.

Prenatal Exposures to Common Phthalates and Prevalent Phthalate Alternatives and Infant DNA Methylation at Birth

Phthalates are a diverse group of chemicals used in consumer products. Because they are so widespread, exposure to these compounds is nearly unavoidable. Recently, growing scientific consensus has suggested that phthalates produce health effects in developing infants and children. These effects may be mediated through mechanisms related to the epigenome, the constellation of mitotically heritable chemical marks and small compounds that guide transcription and translation. The present study examined the relationship between prenatal, first-trimester exposure of seven phthalates and epigenetics in two pregnancy cohorts (n = 262) to investigate sex-specific alterations in infant blood DNA methylation at birth (cord blood or neonatal blood spots). Prenatal exposure to several phthalates was suggestive of association with altered DNA methylation at 4 loci in males (all related to ΣDEHP) and 4 loci in females (1 related to ΣDiNP; 2 related to BBzP; and 1 related to MCPP) at a cutoff of q < 0.2. Additionally, a subset of dyads (n = 79) was used to interrogate the relationships between two compounds increasingly used as substitutions for common phthalates (ΣDINCH and ΣDEHTP) and cord blood DNA methylation. ΣDINCH, but not ΣDEHTP, was suggestive of association with DNA methylation (q < 0.2). Together, these results demonstrate that prenatal exposure to both classically used phthalate metabolites and their newer alternatives is associated with sex-specific infant DNA methylation. Research and regulatory actions regarding this chemical class should consider the developmental health effects of these compounds and aim to avoid regrettable substitution scenarios in the present and future.

Navigating the hydroxymethylome: experimental biases and quality control tools for the tandem bisulfite and oxidative bisulfite Illumina microarrays

Aim: Tandem bisulfite (BS) and oxidative bisulfite (oxBS) conversion on DNA followed by hybridization to Infinium HumanMethylation BeadChips allows nucleotide resolution of 5-hydroxymethylcytosine genome-wide. Here, the authors compared data quality acquired from BS-treated and oxBS-treated samples. Materials & methods: Raw BeadArray data from 417 pairs of samples across 12 independent datasets were included in the study. Probe call rates were compared between paired BS and oxBS treatments controlling for technical variables. Results: oxBS-treated samples had a significantly lower call rate. Among technical variables, DNA-specific extraction kits performed better with higher call rates after oxBS conversion. Conclusion: The authors emphasize the importance of quality control during oxBS conversion to minimize information loss and recommend using a DNA-specific extraction kit for DNA extraction and an oxBSQC package for data preprocessing.

GMQN: A Reference-Based Method for Correcting Batch Effects and Probe Bias in HumanMethylation BeadChip

The Illumina HumanMethylation BeadChip is one of the most cost-effective methods to quantify DNA methylation levels at single-base resolution across the human genome, which makes it a routine platform for epigenome-wide association studies. It has accumulated tens of thousands of DNA methylation array samples in public databases, providing great support for data integration and further analysis. However, the majority of public DNA methylation data are deposited as processed data without background probes which are widely used in data normalization. Here, we present Gaussian mixture quantile normalization (GMQN), a reference based method for correcting batch effects as well as probe bias in the HumanMethylation BeadChip. Availability and implementation: https://github.com/MengweiLi-project/gmqn.

The Relationship between Epigenetic Age and Myocardial Infarction/Acute Coronary Syndrome in a Population-Based Nested Case-Control Study

We investigated the relationship between ‘epigenetic age’ (EA) derived from DNA methylation (DNAm) and myocardial infarction (MI)/acute coronary syndrome (ACS). A random population sample was examined in 2003/2005 (n = 9360, 45–69, the HAPIEE project) and followed up for 15 years. From this cohort, incident MI/ACS (cases, n = 129) and age- and sex-stratified controls (n = 177) were selected for a nested case-control study. Baseline EA (Horvath’s, Hannum’s, PhenoAge, Skin and Blood) and the differences between EA and chronological age (CA) were calculated (ΔAHr, ΔAHn, ΔAPh, ΔASB). EAs by Horvath’s, Hannum’s and Skin and Blood were close to CA (median absolute difference, MAD, of 1.08, –1.91 and –2.03 years); PhenoAge had MAD of −9.29 years vs. CA. The adjusted odds ratios (ORs) of MI/ACS per 1–year increments of ΔAHr, ΔAHn, ΔASB and ΔAPh were 1.01 (95% CI 0.95–1.07), 1.01 (95% CI 0.95–1.08), 1.02 (95% CI 0.97–1.06) and 1.01 (0.93–1.09), respectively. When classified into tertiles, only the highest tertile of ΔAPh showed a suggestion of increased risk of MI/ACS with OR 2.09 (1.11–3.94) independent of age and 1.84 (0.99–3.52) in the age- and sex-adjusted model. Metabolic modulation may be the likely mechanism of this association. In conclusion, this case-control study nested in a prospective population-based cohort did not find strong associations between accelerated epigenetic age markers and risk of MI/ACS. Larger cohort studies are needed to re-examine this important research question.

Healthy eating patterns and epigenetic measures of biological age

BACKGROUND

Healthy eating is associated with lower risks of disease and mortality, but the mechanisms underlying these associations are unclear. Age is strongly related to health outcomes, and biological age can be estimated using the blood methylome.

OBJECTIVES

To determine whether healthy eating patterns are associated with methylation-based measures of biological age.

METHODS

Among women in the Sister Study, we calculated scores on 4 recommendation-based healthy eating indexes [Dietary Approaches to Stop Hypertension diet, Healthy Eating Index-2015, Alternative Healthy Eating Index (aHEI-2010), and the Alternative Mediterranean diet] using a validated 110-item Block FFQ completed at enrollment. Genome-wide DNA methylation data were generated using the HumanMethylation450 BeadChip on whole blood samples collected at enrollment from a case-cohort sample of 2694 women and were used to calculate 4 measures of epigenetic age acceleration (Hannum AgeAccel, Horvath AgeAccel, PhenoAgeAccel, and GrimAgeAccel). Linear regression models, adjusted for covariates and cohort sampling weights, were used to examine cross-sectional associations between eating patterns and measures of biological age.

RESULTS

All 4 healthy eating indexes had inverse associations with epigenetic age acceleration, most notably with PhenoAgeAccel and GrimAgeAccel. Of these, the strongest associations were for aHEI-2010 [per 1-SD increase in diet quality, PhenoAgeAccel β = -0.5 y (95% CI: -0.8 to -0.2 y) and GrimAgeAccel β = -0.4 y (95% CI: -0.6 to -0.3 y)]. Although effect modification was not observed for most lifestyle factors, in analyses stratified by physical activity, the benefits of a healthy diet on epigenetic age acceleration were more pronounced among women who did not meet physical activity guidelines (reporting <2.5 h/wk of exercise).

CONCLUSIONS

Higher diet quality is inversely associated with methylation-based measures of biological age. Improving diet could have the most benefits in lowering biological age among women with lower levels of physical activity. This trial was registered at clinicaltrials.gov as NCT00047970.

Short- and intermediate-term exposure to ambient fine particulate elements and leukocyte epigenome-wide DNA methylation in older men: the Normative Aging Study

BACKGROUND

Several epigenome-wide association studies (EWAS) of ambient particulate matter with aerodynamic diameter ≤ 2.5 µm (PM_2.5) have been reported. However, EWAS of PM_2.5 elements (PEs), reflecting different emission sources, are very limited.

OBJECTIVES

We performed EWAS of short- and intermediate-term exposure to PM_2.5 and 13 PEs. We hypothesized that significant changes in DNAm may vary by PM_2.5 mass and its elements.

METHODS

We repeatedly collected blood samples in the Normative Aging Study and measured leukocyte DNA methylation (DNAm) with the Illumina HumanMethylation450K BeadChip. We collected daily PM_2.5 and 13 PEs at a fixed central site. To estimate the associations between each PE and DNAm at individual cytosine-phosphate-guanine (CpG) sites, we incorporated a distributed-lag (0-27 d) term in the setting of median regression with subject-specific intercept and examined cumulative lag associations. We also accounted for selection bias due to loss to follow-up and mortality prior to enrollment. Significantly differentially methylated probes (DMPs) were identified using Bonferroni correction for multiple testing. We further conducted regional and pathway analyses to identify significantly differentially methylated regions (DMRs) and pathways.

RESULTS

We included 695 men with 1,266 visits between 1999 and 2013. The subjects had a mean age of 75 years. The significant DMPs, DMRs, and pathways varied by to PM_2.5 total mass and PEs. For example, PM_2.5 total mass was associated with 2,717 DMPs and 10,470 DMRs whereas Pb was associated with 3,173 DMPs and 637 DMRs. The identified pathways by PM_2.5 mass were mostly involved in mood disorders, neuroplasticity, immunity, and inflammation, whereas the pathways associated with motor vehicles (BC, Cu, Pb, and Zn) were related with cardiovascular disease and cancer (e.g., "PPARs signaling").

CONCLUSIONS

PM_2.5 and PE were associated with methylation changes at multiple probes and along multiple pathways, in ways that varied by particle components.

Iron homeostasis pathway DNA methylation trajectories reveal a role for STEAP3 metalloreductase in patient outcomes after aneurysmal subarachnoid hemorrhage

BACKGROUND

Following aneurysmal subarachnoid hemorrhage (aSAH), the brain is susceptible to ferroptosis, a type of iron-dependent cell death. Therapeutic intervention targeting the iron homeostasis pathway shows promise for mitigating ferroptosis and improving recovery in animal models, but little work has been conducted in humans. DNA methylation (DNAm) plays a key role in gene expression and brain function, plasticity, and injury recovery, making it a potentially useful biomarker of outcomes or therapeutic target for intervention. Therefore, in this longitudinal, observational study, we examined the relationships between trajectories of DNAm in candidate genes related to iron homeostasis and acute (cerebral vasospasm and delayed cerebral ischemia) and long-term (Glasgow Outcome Scale [GOS, unfavorable = 1-3] and death) patient outcomes after aSAH.

RESULTS

Longitudinal, genome-wide DNAm data were generated from DNA extracted from post-aSAH cerebrospinal fluid (n = 260 participants). DNAm trajectories of 637 CpG sites in 36 candidate genes related to iron homeostasis were characterized over 13 days post-aSAH using group-based trajectory analysis, an unsupervised clustering method. Significant associations were identified between inferred DNAm trajectory groups at several CpG sites and acute and long-term outcomes. Among our results, cg25713625 in the STEAP3 metalloreductase gene (STEAP3) stood out. Specifically, in comparing the highest cg25713625 DNAm trajectory group with the lowest, we observed significant associations (i.e., based on p-values less than an empirical significance threshold) with unfavorable GOS at 3 and 12 months (OR = 11.7, p = 0.0006 and OR = 15.6, p = 0.0018, respectively) and death at 3 and 12 months (OR = 19.1, p = 0.0093 and OR = 12.8, p = 0.0041, respectively). These results were replicated in an independent sample (n = 100 participants) observing significant associations with GOS at 3 and 12 months (OR = 8.2, p = 0.001 and OR = 6.3, p = 0.0.0047, respectively) and death at 3 months (OR = 2.3, p = 0.008) and a suggestive association (i.e., p-value < 0.05 not meeting an empirical significance threshold) with death at 12 months (OR = 2.0, p = 0.0272). In both samples, an additive effect of the DNAm trajectory group was observed as the percentage of participants with unfavorable long-term outcomes increased substantially with higher DNAm trajectory groups.

CONCLUSION

Our results support a role for DNAm of cg25713625/STEAP3 in recovery following aSAH. Additional research is needed to further explore the role of DNAm of cg25713625/STEAP3 as a biomarker of unfavorable outcomes, or therapeutic target to improve outcomes, to translate these findings clinically.

The ENmix DNA methylation analysis pipeline for Illumina BeadChip and comparisons with seven other preprocessing pipelines

Background

Illumina DNA methylation arrays are high-throughput platforms for cost-effective genome-wide profiling of individual CpGs. Experimental and technical factors introduce appreciable measurement variation, some of which can be mitigated by careful “preprocessing” of raw data.

Methods

Here we describe the ENmix preprocessing pipeline and compare it to a set of seven published alternative pipelines (ChAMP, Illumina, SWAN, Funnorm, Noob, wateRmelon, and RnBeads). We use two large sets of duplicate sample measurements with 450 K and EPIC arrays, along with mixtures of isogenic methylated and unmethylated cell line DNA to compare raw data and that preprocessed via different pipelines.

Results

Our evaluations show that the ENmix pipeline performs the best with significantly higher correlation and lower absolute difference between duplicate pairs, higher intraclass correlation coefficients (ICC) and smaller deviations from expected methylation level in mixture experiments. In addition to the pipeline function, ENmix software provides an integrated set of functions for reading in raw data files from mouse and human arrays, quality control, data preprocessing, visualization, detection of differentially methylated regions (DMRs), estimation of cell type proportions, and calculation of methylation age clocks. ENmix is computationally efficient, flexible and allows parallel computing. To facilitate further evaluations, we make all datasets and evaluation code publicly available.

Conclusion

Careful selection of robust data preprocessing methods is critical for DNA methylation array studies. ENmix outperformed other pipelines in our evaluations to minimize experimental variation and to improve data quality and study power.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-021-01207-1.

An exploratory study of white blood cell proportions across preeclamptic and normotensive pregnancy by self-identified race in individuals with overweight or obesity

Objective: Examine white blood cell (WBC) proportions across preeclamptic (n = 28 cases) and normotensive (n = 28 controls) pregnancy in individuals with overweight/obesity.Methods: WBC proportions were inferred from genome-wide DNA methylation data and compared by case/control status and self-identified race.Results: In Trimester 1, ean B cell proportions were suggestively lower in cases in the overall sample and significantly lower in White participants but not in Black participants. More significant WBC proportion changes were observed across normotensive than preeclamptic pregnancy.Conclusions: These findings in a small sample demonstrate need for additional studies investigating the relationship between self-identified race and WBCs in pregnancy.

Associations between infant sex and DNA methylation across umbilical cord blood, artery, and placenta samples

DNA methylation (DNAm) is vulnerable to dysregulation by environmental exposures during epigenetic reprogramming that occurs in embryogenesis. Sexual dimorphism in environmentally induced DNAm dysregulation has been identified and therefore it is important to understand sex-specific DNAm patterns. DNAm at several autosomal sites has been consistently associated with sex in cord blood and placental foetal tissues. However, there is limited research comparing sex-specific DNAm across tissues, particularly differentially methylated regions (DMRs). This study leverages DNAm data measured using the Illumina HumanMethylation450 BeadChip in cord blood (N = 179), placenta (N = 229), and umbilical artery samples (N = 229) in the PRogramming of Intergenerational Stress Mechanisms (PRISM) cohort to identify autosomal DMRs and differentially methylated positions (DMPs). A replication analyses was conducted in an independent cohort (GEO Accession GSE129841). We identified 183, 257, and 419 DMRs and 2119, 2281, and 3405 DMPs (p_Bonferroni < 0.05) in cord blood, placenta, and artery samples, respectively. Thirty-nine DMRs overlapped in all three tissues, overlapping with genes involved in spermatogenesis (NKAPL, PIWIL2 and AURKC) and X-inactivation (LRIF1). In replication analysis, 85% of DMRs overlapped with those identified in PRISM. Overall, DMRs and DMPs had higher methylation levels among females in cord blood and artery samples, but higher methylation levels among males in placenta samples. Further research is necessary to understand biological mechanisms that contribute to differences in sex-specific DNAm signatures across tissues, as well as to determine if sexual dimorphism in the epigenome impacts response to environmental stressors.