Identifying mislabeled and contaminated DNA methylation microarray data: an extended quality control toolset with examples from GEO

Epigenome-wide association study of pregnancy exposure to green space and placental DNA methylation

Green space exposure during pregnancy has been associated with lower risk of adverse birth outcomes, but the biological mechanisms remain unclear. Epigenetic changes, such as DNA methylation (DNAm), may contribute to this association. The placenta, crucial for foetal development, has been understudied in relation to prenatal green space exposure and DNAm on a genome-wide scale. Here, we aimed to investigate the association between green space exposure during pregnancy and epigenome-wide placental DNAm in 550 mother-child pairs from the Barcelona Life Study Cohort (BiSC) in Spain. Green space exposure was assessed as (i) residential surrounding greenness (satellite-based Normalized Difference Vegetation Index (NDVI) in buffers of 100 m, 300 m and 500 m), (ii) residential distance to the nearest major green space (meters), (iii) use of green space (hours/week), and (iv) visual access to greenery through the home window (≥half of the view). Placental DNAm was measured with the EPIC array. Differentially methylated positions (DMPs) were identified using robust linear regression models adjusted for covariates, while differentially methylated regions (DMRs) were identified using the dmrff method. After Bonferroni correction, cg14852540, annotated to SLC25A10 gene, showed an inverse association with residential greenness within 500 m buffer. Additionally, 101 DMPs were suggestively significant (p-values <1 × 10-5) and annotated to genes involved in glucocorticoid-related pathways, inflammatory response, oxidative stress response, and oocyte maturation. No DMRs were identified. Overall, we identified an association between residential greenness and DNAm levels at one CpG in the SLC25A10 gene. Larger studies are needed to validate these findings and understand the biological pathways.

Discrimination, Coping, and DNAm Accelerated Aging Among African American Mothers of the InterGEN Study

Background: Racial discrimination experiences are associated with the activation of stress biology pathways and signs of accelerated biological aging, including alterations in DNA methylation (DNAm). Coping strategies may mitigate stress from racial discrimination and protect against long-term adverse health outcomes. Methods: We conducted a secondary analysis of data from the Intergenerational Impact of Genetic and Psychological Factors on Blood Pressure cohort, an all-African-American sample, to test the hypothesis that social support can protect against accelerated biological aging associated with experiences of racial discrimination. We measured biological aging from saliva DNAm using six epigenetic clocks. Clock values were residualized on participant age and the estimated proportion of epithelial cells contributing to the DNA sample and standardized to M = 0, SD = 1 within the analysis sample. The primary analysis was focused on the second-generation PhenoAge and GrimAge clocks and the third-generation DunedinPACE "speedometer," which previous studies have linked with racial discrimination. Results: In our sample (n = 234; mean age = 31.9 years; SD = 5.80), we found evidence consistent with our hypothesis in the case of the PhenoAge clock, but not the other clocks. Among mothers who did not seek social support, experiences of racial discrimination were associated with an older PhenoAge (b = 0.26, 95% CI = 0.02-0.50, p = 0.03). However, social-support seeking mitigated this risk; at the highest levels of social support, no adverse consequences of discrimination were observed (interaction b = -0.01, 95% CI = -0.02--0.00, p = 0.03). Conclusions: The replication of results is needed. Future research should also investigate additional adaptive and maladaptive coping strategies utilized by African American women and mothers to identify protective measures that influence health outcomes.

Isolating the effects of HIV infection and HIV exposure on epigenetic profiles in infants using historical data from the Mothers and Infants Cohort Study

BACKGROUND

Epigenetics offers insight into the mechanisms by which early life HIV infection and HIV exposure in utero affects offspring health. However, due to the widespread use of antiretroviral therapy (ART) during pregnancy/infancy, contemporary studies are unable to disentangle effects of HIV from ART exposure on epigenetic profiles.

METHODS

Using historical specimens collected before widespread use of ART (1985-1991), we compared DNA methylation (DNAm) profiles among infants with perinatally-acquired HIV (PHIV), HIV-exposed but uninfected (HEU), and HIV-unexposed uninfected (HUU). DNAm in peripheral blood mononuclear cells collected at 3 and 12 months of age (36 PHIV, 33 HEU, and 33 HUU) was profiled using the Illumina Infinium MethylationEPIC BeadChip. We tested for differentially methylated (DM) CpG sites between groups at 3 and 12 months, adjusting for sex, race/ethnicity, and cell type proportions. Biological pathway enrichment analyses were conducted.

FINDINGS

Comparing PHIV to HEU, there were 2 DM sites at 3 months and 11 at 12 months. Comparing PHIV to HUU, there was 1 DM CpG site at 3 months and 6 at 12 months. Immune-related pathways, including interferon-mediated signalling pathways were enriched. HIV exposure was not associated with any variation in DNA methylation, as no differences were detected between HEU vs. HUU at 3 or 12 months.

INTERPRETATION

HIV infection (in the absence of ART during pregnancy/infancy) was associated with DNA methylation changes at 3 and 12 months of life in infants. Differential methylation in PHIV is related to immune processes and HIV exposure in the absence of infection does not contribute to differential methylation.

FUNDING

This study was supported by funding from the National Institutes of Health (R21HD104558 to SS, K01DA053157 to SS, P30ES019776 to CJM, and P30ES005022 to ESB.

Applying blood-derived epigenetic algorithms to saliva: cross-tissue similarity of DNA-methylation indices of aging, physiology, and cognition

BACKGROUND

Epigenetic algorithms of aging, health, and cognition, based on DNA-methylation (DNAm) patterns, are prominent tools for measuring biological age and have been linked to age-related diseases, cognitive decline, and mortality. While most of these methylation profile scores (MPSs) are developed in blood tissue, there is growing interest in using less invasive tissues like saliva. The aim of the current study is to probe the cross-tissue intraclass correlation coefficients (ICCs) of MPSs developed in blood applied to saliva DNAm from the same people. While our primary focus is on MPSs that were previously found to be robustly correlated with social determinants of health, including second- and third-generation clocks and MPSs of physiology and cognition, we also report ICC values for first-generation clocks to enable comparison across metrics. We pooled three publicly available datasets that had both saliva and blood DNAm from the same individuals (total n = 107, aged 5-74 years), corrected MPSs for cell composition within each tissue, and computed the cross-tissue ICCs.

RESULTS

We found that after correcting for cell composition, saliva-blood cross-tissue ICCs were moderate for second- and third-generation indices of aging and MPSs of physiology and cognition. Specifically, PCGrimAge had the highest ICC (0.76), followed by PCPhenoAge (0.72), a measure of cognitive performance (Epigenetic-g, 0.69), DunedinPACE (0.68), PCGrimAge Acceleration (0.67), PCPhenoAge Acceleration (0.66), an MPS of hs-CRP (0.58), and BMI (0.54). These ICCs appear lower than previous reports on within-tissue ICCs (saliva ICCs range from 0.67 to 0.85, blood ICCs range from 0.73 to 0.93). Cross-tissue ICCs values for first-generation biological age acceleration measures were poor, ranging from 0.19 to 0.25.

CONCLUSIONS

Our findings suggest that applying second- and third-generation MPSs of biological age acceleration and related phenotypes developed in blood to saliva DNAm results in moderate cross-tissue similarity and the precise cross-tissue correspondence differs by measure. While the degree of cross-tissue similarity of several MPSs may suffice for some research settings, it may not be suitable in clinical or commercial applications. Collection of both blood and saliva DNAm samples is necessary to validate existing algorithms and to customize MPSs in saliva DNAm.

Epigenetic age acceleration predicts subject-specific white matter degeneration in the human brain

Epigenetic clocks provide powerful tools for estimating health and lifespan but their ability to predict brain degeneration and neuronal damage during the aging process is unknown. In this study, we use GrimAge, an epigenetic clock correlated to several blood plasma proteins, to longitudinally investigate brain cellular microstructure in axonal white matter from a cohort of healthy aging individuals. A specific focus was made on white matter hyperintensities, a visible neurological manifestation of small vessel disease, and the axonal pathways throughout each individual's brain affected by their unique white matter hyperintensity location and volume. 98 subjects over 55 years of age were scanned at baseline with 41 returning for a follow-up scan 2 years later. Using diffusion MRI lesionometry, we reconstructed subject-specific networks of affected axonal tracts and examined the diffusion cellular microstructure composition of these areas, both at baseline and longitudinally, for evidence of cellular degeneration. A chronological age-adjusted version of GrimAge was significantly correlated with baseline WMH volume and markers of neuronal decline, indicated by increased extracellular free water, increased intracellular signal, and decreased axonal signal within WMH. By isolating subject-specific axonal regions "lesioned" by crossing through a WMH, age-adjusted GrimAge was also able to predict longitudinal development of similar patterns of neuronal decline throughout the brain. This study is the first to demonstrate WMH lesionometry as a subject-specific precision imaging technique to study degeneration in aging and the first to establish a relationship between accelerated epigenetic GrimAge and brain cellular microstructure in humans.

Comprehensive guide for epigenetics and transcriptomics data quality control

Host response to environmental exposures such as pathogens and chemicals can include modifications to the epigenome and transcriptome. Improved signature discovery, including the identification of the agent and timing of exposure, has been enabled by advancements in assaying techniques to detect RNA expression, DNA base modifications, histone modifications, and chromatin accessibility. The interrogation of the epigenome and transcriptome cascade requires analyzing disparate datasets from multiple assay types, often at single-cell resolution, derived from the same biospecimen. However, there remains a paucity of rigorous quality control standards of those datasets that reflect quality assurance of the underlying assay. This guide outlines a comprehensive suite of metrics that can be used to ensure quality from 11 different epigenetics and transcriptomics assays. Recommended mitigative actions to address failed metrics are provided. The workflow presented aims to improve benchwork protocols and dataset quality to enable accurate discovery of exposure signatures.

Maternal epigenetic index links early neglect to later neglectful care and other psychopathological, cognitive, and bonding effects

BACKGROUND

Past experiences of maltreatment and life adversity induce DNA methylation changes in adults, but less is known about their impact on mothers' maladaptive neglectful parenting and its negative effects. We performed an epigenome-wide association study to investigate the role of DNA methylation levels in mothers with neglectful care, who were exposed to childhood maltreatment and neglect, and their current negative effects. Saliva DNA methylation was determined with the Illumina Human Methylation EPIC BeadChip v1. The individual epigenome was the input to a machine learning algorithm for trajectory inference, which assigned a specific state to each mother in the progression from healthy controls to the extreme neglect condition. A compound epigenetic maternal neglect score (EMN) was derived from 138 mothers (n = 51 in the neglectful group; n = 87 in the control non-neglectful group) having young children. Differential methylation between groups was utilized to derive the EMNs adjusted for education level, age, experimental variables, and blood cell types in saliva samples.

RESULTS

Structural equation modeling: X2 (29) = 37.81; p = 0.127; RMSEA = 0.048, confirmed that EMNs link their early experience of physical neglect to current reports of psychopathological symptoms, lower cognitive status, and observed poor mother-child emotional availability. A third of the genes annotated to the CpGs that affect EMNs are related to cognitive impairment and neurodegenerative and psychopathological disorders.

CONCLUSIONS

EMNs are a novel index to assess the contribution of DNA methylations as a neglected girl to later neglectful caregiving behavior and other negative effects. The evidence provided expands the possibilities for earlier interventions on the neglect condition to prevent and ameliorate the direct or indirect epigenetic impact of maternal adversities on mother-child care, helping to break the cycle of maltreatment.

Placental PFAS concentrations are associated with perturbations of placental DNA methylation

The placenta is crucial for fetal development, is affected by PFAS toxicity, and evidence is accumulating that gestational PFAS perturb the epigenetic activity of the placenta. Gestational PFAS exposure can adversely affect offspring, yet individual and cumulative impacts of PFAS on the placental epigenome remain underexplored. Here, we conducted an epigenome-wide association study (EWAS) to examine the relationships between placental PFAS levels and DNA methylation in a cohort of mother-infant dyads in Arkansas (N = 151). We measured 17 PFAS in human placental tissues and quantified placental DNA methylation levels via the Illumina EPIC Microarray. We tested for differential DNA methylation with individual PFAS, and with mixtures of multiple PFAS. Our results demonstrated that numerous epigenetic loci were perturbed by PFAS, with PFHxS exhibiting the most abundant effects. Mixture analyses suggested cumulative effects of PFOA and PFOS, while PFHxS may act more independently. We additionally explored whether sex-specific effects may be present and concluded that future large studies should explicitly test for sex-specific effects. The genes that are annotated to our PFAS-associated epigenetic loci are primarily involved in growth processes and cardiometabolic health, while some genes are involved in neurodevelopment. These findings shed light on how prenatal PFAS exposures affect birth outcomes and children's health, emphasizing the importance of understanding PFAS mechanisms in the in-utero environment.

Epigenetic signatures of intergenerational exposure to violence in three generations of Syrian refugees

Maternal trauma influences infant and adult health outcomes and may impact future generations through epigenetic modifications such as DNA methylation (DNAm). Research in humans on the intergenerational epigenetic transmission of trauma effects is limited. In this study, we assessed DNAm signatures of war-related violence by comparing germline, prenatal, and direct exposures to violence across three generations of Syrian refugees. We compared families in which a pregnant grandmother versus a pregnant mother was exposed to violence and included a control group with no exposure to war. We collected buccal swab samples and survey data from mothers and 1-2 children in each of 48 families (n = 131 participants). Based on an epigenome-wide association study (EWAS), we identified differentially methylated regions (DMPs): 14 were associated with germline and 21 with direct exposure to violence. Most DMPs showed the same directionality in DNAm change across germline, prenatal, and direct exposures, suggesting a common epigenetic response to violence. Additionally, we identified epigenetic age acceleration in association with prenatal exposure to violence in children, highlighting the critical period of in utero development. This is the first report of an intergenerational epigenetic signature of violence, which has important implications for understanding the inheritance of trauma.

Epigenome-wide association study of cerebrospinal fluid-based biomarkers of Alzheimer's disease in cognitively normal individuals

INTRODUCTION

Cerebrospinal fluid (CSF) biomarkers of Alzheimer's disease (AD) are reliable predictors of future AD risk. We investigated whether pre-clinical changes in AD CSF biomarkers are reflected in blood DNA methylation (DNAm) levels in cognitively normal participants.

METHODS

We profiled blood-based DNAm with the EPIC array in participants without a diagnosis of cognitive impairment in the Emory Healthy Brain Study (EHBS; N=495) and ADNI (N=122). Their CSF Aβ42, tTau, and pTau levels were quantified using Elecsys immunoassays. We conducted epigenome-wide association studies to assess associations between DNAm and CSF biomarkers of AD.

RESULTS

In EHBS, no loci were Bonferroni-significant after adjusting for confounding factors. In ADNI, two loci were significant, but they were not replicated in EHBS. There was little agreement between the top loci from EHBS and ADNI.

DISCUSSION

Our study showed little evidence of an association between differential blood-based DNAm and pre-clinical AD CSF biomarkers.

Placental and immune cell DNA methylation reference panel for bulk tissue cell composition estimation in epidemiological studies

To distinguish DNA methylation (DNAm) from cell proportion changes in whole placental villous tissue research, we developed a robust cell type-specific DNAm reference to estimate cell composition. We collated new and existing cell type DNAm profiles quantified via Illumina EPIC or 450k microarrays. To estimate cell composition, we deconvoluted whole placental samples (n = 36) with robust partial correlation based on the top 30 hyper- and hypomethylated sites identified per cell type. To test deconvolution performance, we evaluated root mean square error in predicting principal components of DNAm variation in 204 external placental samples. We analyzed DNAm profiles (n = 368,435 sites) from 12 cell types: cytotrophoblasts (n = 18), endothelial cells (n = 19), Hofbauer cells (n = 26), stromal cells (n = 21), syncytiotrophoblasts (n = 4), six lymphocyte types (n = 36), and nucleated red blood cells (n = 11). Median cell composition was consistent with placental biology: 60.9% syncytiotrophoblast, 17.3% stromal, 8.8% endothelial, 3.7% cytotrophoblast, 3.7% Hofbauer, 1.7% nucleated red blood cells, and 1.2% neutrophils. Our expanded reference outperformed an existing reference in predicting DNAm variation (PC1, 15.4% variance explained, IQR = 21.61) with cell composition estimates (mean square error of prediction: 8.62 vs. 10.79, p-value < 0.001). This cell type reference can robustly estimate cell composition from whole placental DNAm data to detect important cell types, reveal biological mechanisms, and improve causal inference.

Parental arsenic exposure and tissue-specific DNA methylation in Bangladeshi infants with spina bifida

An emerging hypothesis linking arsenic toxicity involves altered epigenetic mechanisms, such as DNA methylation. In this study, we examined the relationship between parents' arsenic exposure and DNA methylation in tissues obtained from 28 infants with spina bifida from Bangladesh. We analyzed arsenic in parents' toenails using inductively coupled plasma mass spectrometry (ICP-MS). DNA methylation was measured in infants' dural tissue, buccal swabs, and whole blood using the Illumina Infinium MethylationEPIC BeadChip. We performed epigenome-wide association analyses (EWAS) and tested differentially methylated regions (DMRs). In EWAS, DNA methylation at cg24039697 in dural tissue was positively associated (β = 0.59, p = 7.6 × 10-9) with father's toenail arsenic concentrations, adjusting for covariates. We did not identify any CpG sites related to father's arsenic exposure in the other tissues, or any CpG sites related to mother's arsenic exposure. Gene ontology analysis identified many biological pathways of interest, including the Wnt signaling pathways. We identified several DMRs across the tissues related to arsenic exposure that included probes mapping to genes that have previously been identified in studies of neural tube defects. This study emphasizes the potential impact of arsenic exposure in fathers, often understudied in epidemiological studies, on DNA methylation in a unique neurological tissue specific to spina bifida.

Technical variability across the 450K, EPICv1, and EPICv2 DNA methylation arrays: lessons learned for clinical and longitudinal studies

DNA methylation (DNAm) is the most commonly measured epigenetic mechanism in human populations, with most studies using Illumina arrays to assess DNAm levels. In 2023, Illumina updated their DNAm arrays to the EPIC version 2 (EPICv2), building on prior iterations, namely the EPIC version 1 (EPICv1) and 450K arrays. Whether DNAm measurements are stable across these three generations of arrays has yet not been investigated, limiting the ability of researchers-especially those with longitudinal data-to compare and replicate results across arrays. Here, we present results from a study of 30 child participants (15 male; 15 female) from the Drakenstein Child Health Study, who had DNAm measured on all three of the latest arrays: 450K, EPICv1, and EPICv2. Using these data, we created an annotation of probe quality across arrays, which includes the intraclass correlations, interquartile ranges, correlations, and array bias (i.e., the extent to which DNAm levels were explained by array type) of all CpGs. We also present results from an analysis of sex differences, where we found that CpGs with lower replicability across arrays had higher array-based variance, suggesting this variance metric help guide replication efforts. We also showed that epigenetic age estimates across arrays were more stable when using the principal component versions of epigenetic clocks. Ultimately, this collection of results provides a framework for investigating the replicability and longitudinal stability of epigenetic changes across multiple versions of Illumina DNAm arrays.

Associations between DNA methylation and cognitive function in early-stage hormone receptor-positive breast cancer patients

Background

Approximately one-third of breast cancer (BC) patients show poorer cognitive function (CF) before receiving adjuvant therapy compared with age-matched healthy controls. However, the biological mechanisms driving CF variation in the context of BC remain unclear. In this study, we aimed to identify genes and biological pathways associated with CF in postmenopausal women with early-stage hormone receptor-positive (HR+) BC using DNA methylation (DNAm) data, a dynamic regulator of gene activity.

Methods

Epigenome-wide association studies (EWAS) and differentially methylated region analyses were performed for each CF phenotype (seven objective domains and one subjective phenotype) using DNAm data from whole blood samples (n=109) taken at time of enrollment. Post-EWAS functional analyses were performed to enhance the understanding of the CF-related cytosine-phosphate-guanine (CpG) sites.

Results

When adjusting for age, verbal IQ scores, and global DNAm signature, cg10331779 near CTNND2 (p-value= 9.65 × 10 -9 ) and cg25906741 in MLIP (p-value= 2.01 × 10 -8 ) were associated with processing speed and subjective CF, respectively, while regions in/near SLC6A11 , PRKG1/CSTF2T , and FAM3B for processing speed, and regions in/near PI4KB and SGCE/PEG10 for mental flexibility were differentially methylated. In addition, beta-estradiol was identified as a common upstream regulator for all the CF phenotypes, suggesting an essential role of estrogen in explaining variation in CF of HR+ BC patients.

Conclusions

In our EWAS of 8 CF phenotypes, we found two epigenome-wide significant signals, one at cg10331779 near CTNND2 with processing speed and the other at cg25906741 in MLIP with subjective CF. We also found three differentially methylated regions associated with processing speed and two associated with mental flexibility. These findings need replication in larger cohorts.

Evaluating the association between placenta DNA methylation and cognitive functions in the offspring

The placenta plays a crucial role in protecting the fetus from environmental harm and supports the development of its brain. In fact, compromised placental function could predispose an individual to neurodevelopmental disorders. Placental epigenetic modifications, including DNA methylation, could be considered a proxy of placental function and thus plausible mediators of the association between intrauterine environmental exposures and genetics, and childhood and adult mental health. Although neurodevelopmental disorders such as autism spectrum disorder have been investigated in relation to placenta DNA methylation, no studies have addressed the association between placenta DNA methylation and child's cognitive functions. Thus, our goal here was to investigate whether the placental DNA methylation profile measured using the Illumina EPIC array is associated with three different cognitive domains (namely verbal score, perceptive performance score, and general cognitive score) assessed by the McCarthy Scales of Children's functions in childhood at age 4. To this end, we conducted epigenome-wide association analyses, including data from 255 mother-child pairs within the INMA project, and performed a follow-up functional analysis to help the interpretation of the findings. After multiple-testing correction, we found that methylation at 4 CpGs (cg1548200, cg02986379, cg00866476, and cg14113931) was significantly associated with the general cognitive score, and 2 distinct differentially methylated regions (DMRs) (including 27 CpGs) were significantly associated with each cognitive dimension. Interestingly, the genes annotated to these CpGs, such as DAB2, CEP76, PSMG2, or MECOM, are involved in placenta, fetal, and brain development. Moreover, functional enrichment analyses of suggestive CpGs (p < 1 × 10-4) revealed gene sets involved in placenta development, fetus formation, and brain growth. These findings suggest that placental DNA methylation could be a mechanism contributing to the alteration of important pathways in the placenta that have a consequence on the offspring's brain development and cognitive function.

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.

DNA methylation-based estimators of telomere length show low correspondence with paternal age at conception and other measures of external validity of telomere length

In humans, DNA methylation (DNAm) based estimators of telomere length (TL) have been shown to better predict TL-associated variables (e.g., age, sex, and mortality) than TL itself. The biological significance of DNAm-based estimators of TL (DNAmTL) is unclear. In vitro DNAmTL shortens with cell replications, even when telomerase is maintaining TL. Telomerase is typically suppressed in humans, except in testes. Accordingly, sperm TL increases with age, and offspring with greater paternal age at conception (PAC) have longer TL. Thus, we expect that PAC associations with DNAmTL can shed light on whether in vivo cell replications in the presence of high telomerase activity (production of sperm) shorten DNAmTL or if PAC-lengthened TL causes lengthened DNAmTL. In a pre-registered analysis, using data from 1733 blood samples from the Philippines, we examined the association between paternal age at conception (PAC) and offspring DNAmTL. We did not find an association between PAC and DNAmTL but found a positive association of paternal grandfather's age at father's conception predicting grandchild's DNAmTL. In post hoc analyses, we examined how DNAmTL versus qPCR-measured TL (qPCR-TL) correlated with measures typically associated with TL. Contrary to previous findings, on almost all measures of external validity (correlations with parental TLs, southern blot TL, and age), qPCR-TL outperformed DNAmTL. The "kilobase" units of DNAm-based estimators of TL showed considerable deviations from southern blot-derived kilobase measures. Our findings suggest that DNAmTL is not a reliable index of inherited aspects of TL and underscores uncertainty about the biological meaning of DNAmTL.

Epigenetic link between Agent Orange exposure and type 2 diabetes in Korean veterans

Conflicting findings have been reported regarding the association between Agent Orange (AO) exposure and type 2 diabetes. This study aimed to examine whether AO exposure is associated with the development of type 2 diabetes and to verify the causal relationship between AO exposure and type 2 diabetes by combining DNA methylation with DNA genotype analyses. An epigenome-wide association study and DNA genotype analyses of the blood of AO-exposed and AO-unexposed individuals with type 2 diabetes and that of healthy controls were performed. Methylation quantitative trait locus and Mendelian randomisation analyses were performed to evaluate the causal effect of AO-exposure-identified CpGs on type 2 diabetes. AO-exposed individuals with type 2 diabetes were associated with six hypermethylated CpG sites (cg20075319, cg21757266, cg05203217, cg20102280, cg26081717, and cg21878650) and one hypo-methylated CpG site (cg07553761). Methylation quantitative trait locus analysis showed the methylation levels of some CpG sites (cg20075319, cg20102280, and cg26081717) to be significantly different. Mendelian randomisation analysis showed that CpG sites that were differentially methylated in AO-exposed individuals were causally associated with type 2 diabetes; the reverse causal effect was not significant. These findings reflect the need for further epigenetic studies on the causal relationship between AO exposure and type 2 diabetes.

DNA methylation-based telomere length is associated with HIV infection, physical frailty, cancer, and all-cause mortality

Telomere length (TL) is an important indicator of cellular aging. Shorter TL is associated with several age-related diseases including coronary heart disease, heart failure, diabetes, osteoporosis, and cancer. Recently, a DNA methylation-based TL (DNAmTL) estimator has been developed as an alternative method for directly measuring TL. In this study, we examined the association of DNAmTL with cancer prevalence and mortality risk among people with and without HIV in the Veterans Aging Cohort Study Biomarker Cohort (VACS, N = 1917) and Women's Interagency HIV Study Cohort (WIHS, N = 481). We profiled DNAm in whole blood (VACS) or in peripheral blood mononuclear cells (WIHS) using an array-based method. Cancer prevalence was estimated from electronic medical records and cancer registry data. The VACS Index was used as a measure of physiologic frailty. Models were adjusted for self-reported race and ethnicity, batch, smoking status, alcohol consumption, and five cell types (CD4, CD8, NK, B cell, and monocyte). We found that people with HIV had shorter average DNAmTL than those without HIV infection [beta = -0.25, 95% confidence interval (-0.32, -0.18), p = 1.48E-12]. Greater value of VACS Index [beta = -0.002 (-0.003, -0.001), p = 2.82E-05] and higher cancer prevalence [beta = -0.07 (-0.10, -0.03), p = 1.37E-04 without adjusting age] were associated with shortened DNAmTL. In addition, one kilobase decrease in DNAmTL was associated with a 40% increase in mortality risk [hazard ratio: 0.60 (0.44, 0.82), p = 1.42E-03]. In summary, HIV infection, physiologic frailty, and cancer are associated with shortening DNAmTL, contributing to an increased risk of all-cause mortality.

Epigenetic Aging and Musculoskeletal Outcomes in a Cohort of Women Living With HIV

BACKGROUND

The relationship between accelerated epigenetic aging and musculoskeletal outcomes in women with HIV (WWH) has not been studied.

METHODS

We measured DNA methylation age using the Infinium MethylationEPIC BeadChip in a cohort from the Women's Interagency HIV Study (n = 190) with measures of bone mineral density (BMD) and physical function. We estimated 6 biomarkers of epigenetic aging-epigenetic age acceleration (EAA), extrinsic EAA, intrinsic EAA, GrimAge, PhenoAge, and DNA methylation-estimated telomere length-and evaluated associations of epigenetic aging measures with BMD and physical function. We also performed epigenome-wide association studies to examine associations of DNA methylation signatures with BMD and physical function.

RESULTS

This study included 118 WWH (mean age, 49.7 years; 69% Black) and 72 without HIV (mean age, 48.9 years; 69% Black). WWH had higher EAA (mean ± SD, 1.44 ± 5.36 vs -1.88 ± 5.07; P < .001) and lower DNA methylation-estimated telomere length (7.13 ± 0.31 vs 7.34 ± 0.23, P < .001) than women without HIV. There were no significant associations between accelerated epigenetic aging and BMD. Rather, measures of accelerated epigenetic aging were associated with lower physical function.

CONCLUSIONS

Accelerated epigenetic aging was observed in WWH as compared with women without HIV and was associated with lower physical function in both groups.

Microarray-Based DNA Methylation Profiling: Validation Considerations for Clinical Testing

Microarray-based methylation profiling has emerged as a valuable tool for refining diagnoses and revealing novel tumor subtypes, particularly in central nervous system tumors. Despite the increasing adoption of this technique in clinical genomic laboratories, no technical standards have been published in establishing minimum criteria for test validation. A working group with experience and expertise in DNA-based methylation profiling tests on central nervous system tumors collaborated to develop practical discussion points and focus on important considerations for validating this test in clinical laboratory settings. The experience in validating this methodology in a clinical setting is summarized. Specifically, the advantages and challenges associated with utilizing an in-house classifier compared with a third-party classifier are highlighted. Additionally, experiences in demonstrating the assay's sensitivity and specificity, establishing minimum sample criteria, and implementing quality control metrics are described. As methylation profiling for tumor classification expands to other tumor types and continues to evolve for various other applications, the critical considerations described here are expected to serve as a guidance for future efforts in establishing professional guidelines for this assay.

Prediction of disease-free survival for precision medicine using cooperative learning on multi-omic data

In precision medicine, both predicting the disease susceptibility of an individual and forecasting its disease-free survival are areas of key research. Besides the classical epidemiological predictor variables, data from multiple (omic) platforms are increasingly available. To integrate this wealth of information, we propose new methodology to combine both cooperative learning, a recent approach to leverage the predictive power of several datasets, and polygenic hazard score models. Polygenic hazard score models provide a practitioner with a more differentiated view of the predicted disease-free survival than the one given by merely a point estimate, for instance computed with a polygenic risk score. Our aim is to leverage the advantages of cooperative learning for the computation of polygenic hazard score models via Cox's proportional hazard model, thereby improving the prediction of the disease-free survival. In our experimental study, we apply our methodology to forecast the disease-free survival for Alzheimer's disease (AD) using three layers of data. One layer contains epidemiological variables such as sex, APOE (apolipoprotein E, a genetic risk factor for AD) status and 10 leading principal components. Another layer contains selected genomic loci, and the last layer contains methylation data for selected CpG sites. We demonstrate that the survival curves computed via cooperative learning yield an AUC of around $0.7$, above the state-of-the-art performance of its competitors. Importantly, the proposed methodology returns (1) a linear score that can be easily interpreted (in contrast to machine learning approaches), and (2) a weighting of the predictive power of the involved data layers, allowing for an assessment of the importance of each omic (or other) platform. Similarly to polygenic hazard score models, our methodology also allows one to compute individual survival curves for each patient.

Placental and Immune Cell DNA Methylation Reference Panel for Bulk Tissue Cell Composition Estimation in Epidemiological Studies

To distinguish DNA methylation (DNAm) from cell proportion changes in whole placental tissue research, we developed a robust cell type-specific DNAm reference to estimate cell composition. We collated newly collected and existing cell type DNAm profiles quantified via Illumina EPIC or 450k microarrays. To estimate cell composition, we deconvoluted whole placental samples (n=36) with robust partial correlation based on the top 50 hyper- and hypomethylated sites per cell type. To test deconvolution performance, we evaluated RMSE in predicting principal component one of DNAm variation in 204 external placental samples. We analyzed DNAm profiles (n=368,435 sites) from 12 cell types: cytotrophoblasts (n=18), endothelial cells (n=19), Hofbauer cells (n=26), stromal cells (n=21), syncytiotrophoblasts (n=4), six lymphocyte types (n=36), and nucleated red blood cells (n=11). Median cell composition was consistent with placental biology: 60.4% syncytiotrophoblast, 17.1% stromal, 8.8% endothelial, 4.5% cytotrophoblast, 3.9% Hofbauer, 1.7% nucleated red blood cells, and 1.2% neutrophils. Our expanded reference outperformed an existing reference in predicting DNAm variation (15.4% variance explained, IQR=21.61) with cell composition estimates (RMSE:10.51 vs. 11.43, p-value<0.001). This cell type reference can robustly estimate cell composition from whole placental DNAm data to detect important cell types, reveal biological mechanisms, and improve casual inference.

Accelerated epigenetic age is associated with whole-brain functional connectivity and impaired cognitive performance in older adults

While chronological age is a strong predictor for health-related risk factors, it is an incomplete metric that fails to fully characterize the unique aging process of individuals with different genetic makeup, neurodevelopment, and environmental experiences. Recent advances in epigenomic array technologies have made it possible to generate DNA methylation-based biomarkers of biological aging, which may be useful in predicting a myriad of cognitive abilities and functional brain network organization across older individuals. It is currently unclear which cognitive domains are negatively correlated with epigenetic age above and beyond chronological age, and it is unknown if functional brain organization is an important mechanism for explaining these associations. In this study, individuals with accelerated epigenetic age (i.e. AgeAccelGrim) performed worse on tasks that spanned a wide variety of cognitive faculties including both fluid and crystallized intelligence (N = 103, average age = 68.98 years, 73 females, 30 males). Additionally, fMRI connectome-based predictive models suggested a mediating mechanism of functional connectivity on epigenetic age acceleration-cognition associations primarily in medial temporal lobe and limbic structures. This research highlights the important role of epigenetic aging processes on the development and maintenance of healthy cognitive capacities and function of the aging brain.

The association of cigarette smoking with DNA methylation and gene expression in human tissue samples

Cigarette smoking adversely affects many aspects of human health, and epigenetic responses to smoking may reflect mechanisms that mediate or defend against these effects. Prior studies of smoking and DNA methylation (DNAm), typically measured in leukocytes, have identified numerous smoking-associated regions (e.g., AHRR). To identify smoking-associated DNAm features in typically inaccessible tissues, we generated array-based DNAm data for 916 tissue samples from the GTEx (Genotype-Tissue Expression) project representing 9 tissue types (lung, colon, ovary, prostate, blood, breast, testis, kidney, and muscle). We identified 6,350 smoking-associated CpGs in lung tissue (n = 212) and 2,735 in colon tissue (n = 210), most not reported previously. For all 7 other tissue types (sample sizes 38-153), no clear associations were observed (false discovery rate 0.05), but some tissues showed enrichment for smoking-associated CpGs reported previously. For 1,646 loci (in lung) and 22 (in colon), smoking was associated with both DNAm and local gene expression. For loci detected in both lung and colon (e.g., AHRR, CYP1B1, CYP1A1), top CpGs often differed between tissues, but similar clusters of hyper- or hypomethylated CpGs were observed, with hypomethylation at regulatory elements corresponding to increased expression. For lung tissue, 17 hallmark gene sets were enriched for smoking-associated CpGs, including xenobiotic- and cancer-related gene sets. At least four smoking-associated regions in lung were impacted by lung methylation quantitative trait loci (QTLs) that co-localize with genome-wide association study (GWAS) signals for lung function (FEV1/FVC), suggesting epigenetic alterations can mediate the effects of smoking on lung health. Our multi-tissue approach has identified smoking-associated regions in disease-relevant tissues, including effects that are shared across tissue types.

Mother adversity and co-residence time impact mother-child similarity in genome-wide and gene-specific methylation profiles

BACKGROUND

The effects of adverse life events on physical and psychological health, with DNA methylation (DNAm) as a critical underlying mechanism, have been extensively studied. However, the epigenetic resemblance between mother and child in the context of neglectful caregiving, and whether it may be shaped by the emotional impact of maternal stressful events and the duration of co-residence (indexed by child age), remains unknown. The present study examined mother-child similarity in methylation profiles, considering the potential effect of mother adversity, mother empathy, neglect-control group, child age (an index of years of mother-child co-residence), and mother age. Using Illumina Epic arrays, we quantified DNAm in 115 mother-child saliva samples. We obtained a methylation similarity index by computing correlation coefficients between methylation profiles within dyads, for the entire epigenome, and five specific genes related to stress and empathy: NR3C1, FKPB5, OXTR, SCL6A4, and BDNF.

RESULTS

The methylation profiles of the mother-child familial pairs significantly correlated as compared to mother-child random pairs for the entire epigenome and NR3C1, FKBP5, OXTR and BDNF genes. Next, multiple linear regression models observed associations of mother adversity, child age, and neglect-control group on mother-child methylation similarity, only significant in mother-child familial pairs, after correcting for multiple comparisons. Higher mother adversity was associated with lower mother-child methylation similarity for the epigenome-wide analysis, for the BDNF gene, and in the neglect-control group for the OXTR gene. In turn, being an older child (longer co-residence) was associated with higher mother-child methylation similarity.

CONCLUSIONS

Mother adversity and co-residence time are modulating factors in the intergenerational methylation process that offer a window into development-dependent adaptations that can be affected by both hereditary and environmental factors, significantly observed only in biological dyads. A twofold implication for child well-being emerges, one is positive in that children of mothers exposed to life adversity or neglect did not necessarily inherit their methylation patterns. The other is concerning due to the influence of time spent living together, which affects similarity with the mother and potentially increases the risk of inheriting an epigenetic profile associated with future dysfunctional parenting patterns. This underscores the importance of the 'the earlier, the better' recommendation by the Child Protection System, which is not always followed.

The Epigenetic Evolution of Glioma Is Determined by the IDH1 Mutation Status and Treatment Regimen

Tumor adaptation or selection is thought to underlie therapy resistance in glioma. To investigate longitudinal epigenetic evolution of gliomas in response to therapeutic pressure, we performed an epigenomic analysis of 132 matched initial and recurrent tumors from patients with IDH-wildtype (IDHwt) and IDH-mutant (IDHmut) glioma. IDHwt gliomas showed a stable epigenome over time with relatively low levels of global methylation. The epigenome of IDHmut gliomas showed initial high levels of genome-wide DNA methylation that was progressively reduced to levels similar to those of IDHwt tumors. Integration of epigenomics, gene expression, and functional genomics identified HOXD13 as a master regulator of IDHmut astrocytoma evolution. Furthermore, relapse of IDHmut tumors was accompanied by histologic progression that was associated with survival, as validated in an independent cohort. Finally, the initial cell composition of the tumor microenvironment varied between IDHwt and IDHmut tumors and changed differentially following treatment, suggesting increased neoangiogenesis and T-cell infiltration upon treatment of IDHmut gliomas. This study provides one of the largest cohorts of paired longitudinal glioma samples with epigenomic, transcriptomic, and genomic profiling and suggests that treatment of IDHmut glioma is associated with epigenomic evolution toward an IDHwt-like phenotype.

SIGNIFICANCE

Standard treatments are related to loss of DNA methylation in IDHmut glioma, resulting in epigenetic activation of genes associated with tumor progression and alterations in the microenvironment that resemble treatment-naïve IDHwt glioma.

Self-control is associated with health-relevant disparities in buccal DNA-methylation measures of biological aging in older adults

Self-control is a personality dimension that is associated with better physical health and a longer lifespan. Here, we examined (1) whether self-control is associated with buccal and saliva DNA-methylation (DNAm) measures of biological aging quantified in children, adolescents, and adults, and (2) whether biological aging measured in buccal DNAm is associated with self-reported health. Following preregistered analyses, we computed two DNAm measures of advanced biological age (principal-component PhenoAge and GrimAge Acceleration) and a DNAm measure of pace of aging (DunedinPACE) in buccal samples from the German Socioeconomic Panel Study (SOEP-G[ene], n = 1058, age range 0-72, Mage = 42.65) and saliva samples from the Texas Twin Project (TTP, n = 1327, age range 8-20, Mage = 13.50). We found that lower self-control was associated with advanced biological age in older adults (PhenoAge Acceleration β = - .34, [- .51, - .17], p < .001; GrimAge Acceleration β = - .34, [- .49, - .19], p < .001), but not young adults, adolescents or children. These associations remained statistically robust even after correcting for possible confounders such as socioeconomic contexts, BMI, or genetic correlates of low self-control. Moreover, a faster pace of aging and advanced biological age measured in buccal DNAm were associated with self-reported disease (PhenoAge Acceleration: β = .13 [.06, .19], p < .001; GrimAge Acceleration: β = .19 [.12, .26], p < .001; DunedinPACE: β = .09 [.02, .17], p = .01). However, effect sizes were weaker than observations in blood, suggesting that customization of DNAm aging measures to buccal and saliva tissues may be necessary. Our findings are consistent with the hypothesis that self-control is associated with health via pathways that accelerate biological aging in older adults.

Immunological shifts during early-stage Parkinson's disease identified with DNA methylation data on longitudinally collected blood samples

Parkinson's disease (PD) is the second most common neurodegenerative disease in the United States. Decades before motor symptoms manifest, non-motor symptoms such as hyposmia and rapid eye movement (REM) sleep behavior disorder are highly predictive of PD. Previous immune profiling studies have identified alterations to the proportions of immune cells in the blood of clinically defined PD patients. However, it remains unclear if these phenotypes manifest before the clinical diagnosis of PD. We utilized longitudinal DNA methylation (DNAm) microarray data from the Parkinson's Progression Marker's Initiative (PPMI) to perform immune profiling in clinically defined PD and prodromal PD patients (Prod). We identified previously reported changes in neutrophil, monocyte, and T cell numbers in PD patients. Additionally, we noted previously unrecognized decreases in the naive B cell compartment in the defined PD and Prod patient group. Over time, we observed the proportion of innate immune cells in PD blood increased, but the proportion of adaptive immune cells decreased. We identified decreases in T and B cell subsets associated with REM sleep disturbances and early cognitive decline. Lastly, we identified increases in B memory cells associated with both genetic (LRRK2 genotype) and infectious (cytomegalovirus seropositivity) risk factors of PD. Our analysis shows that the peripheral immune system is dynamic as the disease progresses. The study provides a platform to understand how and when peripheral immune alterations occur in PD and whether intervention at particular stages may be therapeutically advantageous.

Socioeconomic status is negatively associated with immunosenescence but positively associated with inflammation among middle-aged women in Cebu, Philippines

BACKGROUND

Socioeconomic status (SES) gradients in health are well-documented, and while biological pathways are incompletely understood, chronic inflammation and accelerated immune aging (immunosenescence) among lower SES individuals have been implicated. However, previous findings have come from samples in higher income countries, and it is unclear how generalizable they are to lower- and middle-income countries (LMIC) with different infectious exposures and where adiposity-an important contributor to chronic inflammation-might show different SES patterning. To address this gap, we explored associations between SES and inflammation and immunosenescence in a sample of women in Cebu, Philippines.

METHODS

Data came from the mothers of the Cebu Longitudinal Health and Nutrition Survey birth cohort (mean age: 47.7, range: 35-69 years). SES was measured as a combination of annual household income, education level, and assets. Chronic inflammation was measured using C-reactive protein (CRP) in plasma samples from 1,834 women. Immunosenescence was measured by the abundance of exhausted CD8T (CD8 + CD28-CD45RA-) and naïve CD8T and CD4T cells, estimated from DNA methylation in whole blood in a random subsample of 1,028. Possible mediators included waist circumference and a collection of proxy measures of pathogen exposure.

RESULTS

SES was negatively associated with the measures of immunosenescence, with slight evidence for mediation by a proxy measure for pathogen exposure from the household's drinking water source. In contrast, SES was positively associated with CRP, which was explained by the positive association with waist circumference.

CONCLUSIONS

Similar to higher income populations, in Cebu there is an SES-gradient in pathogen exposures and immunosenescence. However, lifestyle changes occurring more rapidly among higher SES individuals is contributing to a positive association between SES and adiposity and inflammation. Our results suggest more studies are needed to clarify the relationship between SES and inflammation and immunosenescence across LMIC.

Sex-specific DNA methylation in saliva from the multi-ethnic Future of Families and Child Wellbeing Study

The prevalence and severity of many diseases differs by sex, potentially due to sex-specific patterns in DNA methylation. Autosomal sex-specific differences in DNA methylation have been observed in cord blood and placental tissue but are not well studied in saliva or in diverse populations. We sought to characterize sex-specific DNA methylation on autosomal chromosomes in saliva samples from children in the Future of Families and Child Wellbeing Study, a multi-ethnic prospective birth cohort containing an oversampling of Black, Hispanic and low-income families. DNA methylation from saliva samples was analysed on 796 children (50.6% male) at both ages 9 and 15 with DNA methylation measured using the Illumina HumanMethylation 450k array. An epigenome-wide association analysis of the age 9 samples identified 8,430 sex-differentiated autosomal DNA methylation sites (P < 2.4 × 10-7), of which 76.2% had higher DNA methylation in female children. The strongest sex-difference was in the cg26921482 probe, in the AMDHD2 gene, with 30.6% higher DNA methylation in female compared to male children (P < 1 × 10-300). Treating the age 15 samples as an internal replication set, we observed highly consistent results between the ages 9 and 15 measurements, indicating stable and replicable sex-differentiation. Further, we directly compared our results to previously published DNA methylation sex differences in both cord blood and saliva and again found strong consistency. Our findings support widespread and robust sex-differential DNA methylation across age, human tissues, and populations. These findings help inform our understanding of potential biological processes contributing to sex differences in human physiology and disease.

Prenatal maternal stress is associated with site-specific and age acceleration changes in maternal and newborn DNA methylation

Prenatal maternal stress has a negative impact on child health but the mechanisms through which maternal stress affects child health are unclear. Epigenetic variation, such as DNA methylation, is a likely mechanistic candidate as DNA methylation is sensitive to environmental insults and can regulate long-term changes in gene expression. We recruited 155 mother-newborn dyads in the Democratic Republic of Congo to investigate the effects of maternal stress on DNA methylation in mothers and newborns. We used four measures of maternal stress to capture a range of stressful experiences: general trauma, sexual trauma, war trauma, and chronic stress. We identified differentially methylated positions (DMPs) associated with general trauma, sexual trauma, and war trauma in both mothers and newborns. No DMPs were associated with chronic stress. Sexual trauma was positively associated with epigenetic age acceleration across several epigenetic clocks in mothers. General trauma and war trauma were positively associated with newborn epigenetic age acceleration using the extrinsic epigenetic age clock. We tested the top DMPs for enrichment of DNase I hypersensitive sites (DHS) and found no enrichment in mothers. In newborns, top DMPs associated with war trauma were enriched for DHS in embryonic and foetal cell types. Finally, one of the top DMPs associated with war trauma in newborns also predicted birthweight, completing the cycle from maternal stress to DNA methylation to newborn health outcome. Our results indicate that maternal stress is associated with site-specific changes in DNAm and epigenetic age acceleration in both mothers and newborns.

Stress and DNA Methylation of Blood Leukocytes among Pregnant Latina Women

Latinas experience physical and psychological stressors in pregnancy leading to increased morbidity and higher risk for adverse birth outcomes. Epigenetic changes, including DNA methylation (DNAm), have been proposed as markers to create more refined risk stratification, yet few of these studies have examined these changes in Latinas. We conducted a secondary analysis of stored blood leukocytes of Latina women (n = 58) enrolled in a larger National Institutes of Health funded R01 project (2011-2016). We examined DNAm on eight candidate stress genes to compare physically and psychologically stressed participants to healthy (low stress) participants. We found unique CpGs that were differentially methylated in stressed women early- and mid-pregnancy compared to the healthy group, though none remained significant after FDR correction. Both physical and psychological stress were associated with hypomethylation at two consecutive CpG sites on NR3C1 in early pregnancy and one CpG site on NR3C1 in mid-pregnancy before adjustment. Stress was also associated with hypomethylation at two CpG sites on FKBP5 in early and mid-pregnancy but were no longer significant after FDR adjustment. Though we did not find statistically significant differences in DNAm during pregnancy between stressed and healthy women in this sample, signals were consistent with previous findings. Future work in larger samples should further examine the associations between stress and DNAm in pregnancy as this mechanism may explain underlying perinatal health inequities.

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.

Negative parenting, epigenetic age, and psychological problems: prospective associations from adolescence to young adulthood

BACKGROUND

Epigenetic clocks are based on DNA methylation levels of several genomic loci and have been developed as indices of biological aging. Studies examining the effects of stressful environmental exposures have shown that stress is associated with differences between epigenetic age and chronological age (i.e., Epigenetic Age acceleration, EA). This pre-registered longitudinal study examined the long-term effects of negative parenting and psychological problems throughout adolescence (ages 13-17 years) on EA in late adolescence (age 17 years) and EA changes from late adolescence to young adulthood (age 25 years). Further, it examined how (change in) EA is related to changes in psychological problems from adolescence to young adulthood.

METHODS

We used data from a sample of 434 participants followed from age 13 to age 25, with saliva collected at ages 17 and 25. We estimated EA using four commonly used epigenetic clocks and analyzed the data using Structural Equation Modeling.

RESULTS

While negative parenting was not related to EA nor change in EA, (change in) EA was related to developmental indices such as externalizing problems and self-concept clarity.

CONCLUSIONS

Declining psychological well-being during young adulthood was preceded by EA.

Multi-omic atlas of the parahippocampal gyrus in Alzheimer's disease

Alzheimer's disease (AD) is the most common form of dementia worldwide, with a projection of 151 million cases by 2050. Previous genetic studies have identified three main genes associated with early-onset familial Alzheimer's disease, however this subtype accounts for less than 5% of total cases. Next-generation sequencing has been well established and holds great promise to assist in the development of novel therapeutics as well as biomarkers to prevent or slow the progression of this devastating disease. Here we present a public resource of functional genomic data from the parahippocampal gyrus of 201 postmortem control, mild cognitively impaired (MCI) and AD individuals from the Mount Sinai brain bank, of which whole-genome sequencing (WGS), and bulk RNA sequencing (RNA-seq) were previously published. The genomic data include bulk proteomics and DNA methylation, as well as cell-type-specific RNA-seq and assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) data. We have performed extensive preprocessing and quality control, allowing the research community to access and utilize this public resource available on the Synapse platform at https://doi.org/10.7303/syn51180043.2 .

eoPred: predicting the placental phenotype of early-onset preeclampsia using public DNA methylation data

Background: A growing body of literature has reported molecular and histological changes in the human placenta in association with preeclampsia (PE). Placental DNA methylation (DNAme) and transcriptomic patterns have revealed molecular subgroups of PE that are associated with placental histopathology and clinical phenotypes of the disease. However, the clinical and molecular heterogeneity of PE both across and within subtypes complicates the study of this disease. PE is most strongly associated with placental pathology and adverse fetal and maternal outcomes when it develops early in pregnancy. We focused on placentae from pregnancies affected by preeclampsia that were delivered before 34 weeks of gestation to develop eoPred, a predictor of the DNAme signature associated with the placental phenotype of early-onset preeclampsia (EOPE). Results: Public data from 83 placental samples (HM450K), consisting of 42 EOPE and 41 normotensive preterm birth (nPTB) cases, was used to develop eoPred-a supervised model that relies on a highly discriminative 45 CpG DNAme signature of EOPE in the placenta. The performance of eoPred was assessed using cross-validation (AUC = 0.95) and tested in an independent validation cohort (n = 49, AUC = 0.725). A subset of fetal growth restriction (FGR) and late-PE cases showed a similar DNAme profile at the 45 predictive CpGs, consistent with the overlap in placental pathology between these conditions. The relationship between the EOPE probability generated by eoPred and various phenotypic variables was also assessed, revealing that it is associated with gestational age, and it is not driven by cell composition differences. Conclusion: eoPred relies on a 45-CpG DNAme signature to predict a homogeneous placental phenotype of EOPE in a discrete or continuous manner. Using this classifier should 1) aid in the study of placental insufficiency and improve the consistency of future placental DNAme studies of PE, 2) facilitate identifying the placental phenotype of EOPE in public data sets and 3) importantly, standardize the placental diagnosis of EOPE to allow better cross-cohort comparisons. Lastly, classification of cases with eoPred will be useful for investigating the relationship between placental pathology and genetic or environmental variables.

Self-control is associated with health-relevant disparities in buccal DNA-methylation measures of biological aging in older adults

Self-control is a personality dimension that is associated with better physical health and a longer lifespan. Here we examined (1) whether self-control is associated with buccal and saliva DNA-methylation (DNAm) measures of biological aging quantified in children, adolescents, and adults, and (2) whether biological aging measured in buccal DNAm is associated with self-reported health. Following preregistered analyses, we computed two DNAm measures of advanced biological age (PhenoAge and GrimAge Acceleration) and a DNAm measure of pace of aging (DunedinPACE) in buccal samples from the German Socioeconomic Panel Study (SOEP-G[ene], n = 1058, age range 0-72, Mage = 42.65) and saliva samples from the Texas Twin Project (TTP, n = 1327, age range 8-20, Mage = 13.50). We found that lower self-control was associated with advanced biological age in older adults (β =-.34), but not young adults, adolescents or children. This association was not accounted for by statistical correction for socioeconomic contexts, BMI, or genetic correlates of low self-control. Moreover, a faster pace of aging and advanced biological age measured in buccal DNAm were associated with worse self-reported health (β =.13 to β = .19). But, effect sizes were weaker than observations in blood, thus customization of DNAm aging measures to buccal and saliva tissues may be necessary. Our findings are consistent with the hypothesis that self-control is associated with health via pathways that accelerate biological aging in older adults.

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.

Association of immune cell composition with the risk factors and incidence of acute coronary syndrome

BACKGROUND

Although immune cells are involved in acute coronary syndrome (ACS), few studies have explored the association of incident ACS with the relative immune cell proportions. We aimed to investigate the association of immune cell proportions with the incidence and risk factors of ACS in the Dongfeng-Tongji cohort.

METHODS

We conducted the analyses with 38,295 subjects from the first follow-up of the Dongfeng-Tongji cohort, including DNA methylation profiles for 1570 individuals. The proportions of immune cell types were observed from routine blood tests or estimated from DNA methylation profiles. For both observed and estimated immune cell proportions, we tested their associations with risk factors of ACS by multivariable linear regression models. In addition, the association of each immune cell proportion with incident ACS was assessed by the Cox regression model and conditional logistic regression model, respectively, adjusting for the risk factors of ACS.

FINDINGS

The proportions of lymphocytes, monocytes, and neutrophils showed strong associations with sex, followed by diabetes. Moreover, sex and current smoking were the two factors with strongest association with the proportions of lymphocyte subtypes. The hazard ratio (HR) and 95% confidence interval (CI) of incident ACS per standard deviation (SD) increase in proportions of lymphocytes and neutrophils were 0.91 (0.85-0.96) and 1.10 (1.03-1.16), respectively. Furthermore, the OR (95% CI) of incident ACS per SD increase in proportions of NK cells, CD4⁺ T cells, and B cells were 0.88 (0.78-0.99), 1.15 (1.03-1.30), and 1.13 (1.00-1.26), respectively.

INTERPRETATION

The proportions of immune cells were associated with several risk factors of ACS, including sex, diabetes, and current smoking. In addition, proportion of neutrophils had a risk effect, while proportion of lymphocytes had a protective effect on the incidence of ACS. The protective effect of lymphocytes was probably driven by NK cells.

An epigenetic clock for gestational age based on human umbilical vein endothelial cells from a diverse population of newborns

Background

Epigenetic clocks are emerging as a useful tool in many areas of research. Many epigenetic clocks have been developed for adults; however, there are fewer clocks focused on newborns and most are trained using blood from European ancestry populations. In this study, we built an epigenetic clock based on primary human umbilical vein endothelial cells from a racially and ethnically diverse population.

Results

Using human umbilical vein endothelial cell [HUVEC]-derived DNA, we calculated epigenetic gestational age using 83 CpG sites selected through elastic net regression. In this study with newborns from different racial/ethnic identities, epigenetic gestational age and clinical gestational age were more highly correlated (r = 0.85), than epigenetic clocks built from adult and other pediatric populations. The correlation was also higher than clocks based on blood samples from newborns with European ancestry. We also found that birth weight was positively associated with epigenetic gestational age acceleration (EGAA), while NICU admission was associated with lower EGAA. Newborns self-identified as Hispanic or non-Hispanic Black had lower EGAA than self-identified as non-Hispanic White.

Conclusions

Epigenetic gestational age can be used to estimate clinical gestational age and may help index neonatal development. Caution should be exercised when using epigenetic clocks built from adults with children, especially newborns. We highlight the importance of cell type-specific epigenetic clocks and general pan tissue epigenetic clocks derived from a large racially and ethnically diverse population.

Adolescent peer struggles predict accelerated epigenetic aging in midlife

This study examined struggles to establish autonomy and relatedness with peers in adolescence and early adulthood as predictors of advanced epigenetic aging assessed at age 30. Participants (N = 154; 67 male and 87 female) were observed repeatedly, along with close friends and romantic partners, from ages 13 through 29. Observed difficulty establishing close friendships characterized by mutual autonomy and relatedness from ages 13 to 18, an interview-assessed attachment state of mind lacking autonomy and valuing of attachment at 24, and self-reported difficulties in social integration across adolescence and adulthood were all linked to greater epigenetic age at 30, after accounting for chronological age, gender, race, and income. Analyses assessing the unique and combined effects of these factors, along with lifetime history of cigarette smoking, indicated that each of these factors, except for adult social integration, contributed uniquely to explaining epigenetic age acceleration. Results are interpreted as evidence that the adolescent preoccupation with peer relationships may be highly functional given the relevance of such relationships to long-term physical outcomes.

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.

Associations of socioeconomic disparities with buccal DNA-methylation measures of biological aging

BACKGROUND

Individuals who are socioeconomically disadvantaged are at increased risk for aging-related diseases and perform less well on tests of cognitive function. The weathering hypothesis proposes that these disparities in physical and cognitive health arise from an acceleration of biological processes of aging. Theories of how life adversity is biologically embedded identify epigenetic alterations, including DNA methylation (DNAm), as a mechanistic interface between the environment and health. Consistent with the weathering hypothesis and theories of biological embedding, recently developed DNAm algorithms have revealed profiles reflective of more advanced aging and lower cognitive function among socioeconomically-at-risk groups. These DNAm algorithms were developed using blood-DNA, but social and behavioral science research commonly collect saliva or cheek-swab DNA. This discrepancy is a potential barrier to research to elucidate mechanisms through which socioeconomic disadvantage affects aging and cognition. We therefore tested if social gradients observed in blood DNAm measures could be reproduced using buccal-cell DNA obtained from cheek swabs.

RESULTS

We analyzed three DNAm measures of biological aging and one DNAm measure of cognitive performance, all of which showed socioeconomic gradients in previous studies: the PhenoAge and GrimAge DNAm clocks, DunedinPACE, and Epigenetic-g. We first computed blood-buccal cross-tissue correlations in n = 21 adults (GEO111165). Cross-tissue correlations were low-to-moderate (r = .25 to r = .48). We next conducted analyses of socioeconomic gradients using buccal DNAm data from SOEP-G (n = 1128, 57% female; age mean = 42 yrs, SD = 21.56, range 0-72). Associations of socioeconomic status with DNAm measures of aging were in the expected direction, but were smaller as compared to reports from blood DNAm datasets (r = - .08 to r = - .13).

CONCLUSIONS

Our findings are consistent with the hypothesis that socioeconomic disadvantage is associated with DNAm indicators of worse physical health. However, relatively low cross-tissue correlations and attenuated effect sizes for socioeconomic gradients in buccal DNAm compared with reports from analysis of blood DNAm suggest that in order to take full advantage of buccal DNA samples, DNAm algorithms customized to buccal DNAm are needed.

DNA Methylation Research in Autologous Hematopoietic Stem Cell Transplant Population

Despite increased sophistication in DNA methylation (DNAm) measurement and methods, conducting studies in specific populations such as the hematopoietic stem cell transplant (HCT) population, presents unique challenges and study design considerations. In this article, we explain the motivation for investigating DNAm in the HCT population, highlighting important study design features and key findings in a longitudinal prospective pilot study of DNAm in 32 patients undergoing autologous HCT in Central Virginia, USA. We also discuss limitations and challenges to generating robust results. We observed that HCT does not prevent high-quality DNA from being extracted from whole blood for DNAm research and that longitudinal prospective studies that span pre- and 2-months post-HCT are feasible. Critically, we did not observe significant impacts of cancer diagnosis, time since transplant, age, or chromosomal sex on overall DNAm data dimensionality. These observations demonstrate that while extreme care is required to ensure generalizable, accurate, and interpretable results, researchers should not avoid HCT-DNAm research simply for fear that the transplant procedure or presence of a cancer diagnosis will prevent meaningful conclusions from being drawn. DNAm is an attractive biomarker that is understudied in patients undergoing HCT and needs to expand to improve precise prediction of HCT outcomes.

Risperidone response in patients with schizophrenia drives DNA methylation changes in immune and neuronal systems

Background: The choice of efficient antipsychotic therapy for schizophrenia relies on a time-consuming trial-and-error approach, whereas the social and economic burdens of the disease call for faster alternatives. Material & methods: In a search for predictive biomarkers of antipsychotic response, blood methylomes of 28 patients were analyzed before and 4 weeks into risperidone therapy. Results: Several CpGs exhibiting response-specific temporal dynamics were identified in otherwise temporally stable methylomes and noticeable global response-related differences were observed between good and bad responders. These were associated with genes involved in immunity, neurotransmission and neuronal development. Polymorphisms in many of these genes were previously linked with schizophrenia etiology and antipsychotic response. Conclusion: Antipsychotic response seems to be shaped by both stable and medication-induced methylation differences.

Genetic and epigenetic signatures associated with plasma oxytocin levels in children and adolescents with autism spectrum disorder

Oxytocin (OT), the brain's most abundant neuropeptide, plays an important role in social salience and motivation. Clinical trials of the efficacy of OT in autism spectrum disorder (ASD) have reported mixed results due in part to ASD's complex etiology. We investigated whether genetic and epigenetic variation contribute to variable endogenous OT levels that modulate sensitivity to OT therapy. To carry out this analysis, we integrated genome-wide profiles of DNA-methylation, transcriptional activity, and genetic variation with plasma OT levels in 290 participants with ASD enrolled in a randomized controlled trial of OT. Our analysis identified genetic variants with novel association with plasma OT, several of which reside in known ASD risk genes. We also show subtle but statistically significant association of plasma OT levels with peripheral transcriptional activity and DNA-methylation profiles across several annotated gene sets. These findings broaden our understanding of the effects of the peripheral oxytocin system and provide novel genetic candidates for future studies to decode the complex etiology of ASD and its interaction with OT signaling and OT-based interventions. LAY SUMMARY: Oxytocin (OT) is an abundant chemical produced by neurons that plays an important role in social interaction and motivation. We investigated whether genetic and epigenetic factors contribute to variable OT levels in the blood. To this, we integrated genetic, gene expression, and non-DNA regulated (epigenetic) signatures with blood OT levels in 290 participants with autism enrolled in an OT clinical trial. We identified genetic association with plasma OT, several of which reside in known autism risk genes. We also show statistically significant association of plasma OT levels with gene expression and epigenetic across several gene pathways. These findings broaden our understanding of the factors that influence OT levels in the blood for future studies to decode the complex presentation of autism and its interaction with OT and OT-based treatment.

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.

Oxytocin receptor DNA methylation is associated with exogenous oxytocin needs during parturition and postpartum hemorrhage

BACKGROUND

The oxytocin receptor gene (OXTR) is regulated, in part, by DNA methylation. This mechanism has implications for uterine contractility during labor and for prevention or treatment of postpartum hemorrhage, an important contributor to global maternal morbidity and mortality.

METHODS

We measured and compared the level of OXTR DNA methylation between matched blood and uterine myometrium to evaluate blood as an indicator of uterine methylation status using targeted pyrosequencing and sites from the Illumina EPIC Array. Next, we tested for OXTR DNA methylation differences in blood between individuals who experienced a postpartum hemorrhage arising from uterine atony and matched controls following vaginal birth. Bivariate statistical tests, generalized linear modeling and Poisson regression were used in the analyses.

RESULTS

Here we show a significant positive correlation between blood and uterine DNA methylation levels at several OXTR loci. Females with higher OXTR DNA methylation in blood had required significantly more exogenous oxytocin during parturition. With higher DNA methylation, those who had oxytocin administered during labor had significantly greater relative risk for postpartum hemorrhage (IRR 2.95, 95% CI 1.53-5.71).

CONCLUSIONS

We provide evidence that epigenetic variability in OXTR is associated with the amount of oxytocin administered during parturition and moderates subsequent postpartum hemorrhage. Methylation can be measured using a peripheral tissue, suggesting potential use in identifying individuals susceptible to postpartum hemorrhage. Future studies are needed to quantify myometrial gene expression in connection with OXTR methylation.