A comparative analysis of cell-type adjustment methods for epigenome-wide association studies based on simulated and real data sets

Epigenetic age acceleration in follicular fluid and markers of ovarian response among women undergoing IVF

STUDY QUESTION

Are markers of epigenetic age acceleration in follicular fluid associated with outcomes of ovarian stimulation?

SUMMARY ANSWER

Increased epigenetic age acceleration of follicular fluid using the Horvath clock, but not other epigenetic clocks (GrimAge and Granulosa Cell), was associated with lower peak estradiol levels and decreased number of total and mature oocytes.

WHAT IS KNOWN ALREADY

In granulosa cells, there are inconsistent findings between epigenetic age acceleration and ovarian response outcomes.

STUDY DESIGN, SIZE, DURATION

Our study included 61 women undergoing IVF at an academic fertility clinic in the New England area who were part of the Environment and Reproductive Health Study (2006-2016).

PARTICIPANTS/MATERIALS, SETTING, METHODS

Participants provided a follicular fluid sample during oocyte retrieval. DNA methylation of follicular fluid was assessed using a genome-wide methylation screening tool. Three established epigenetic clocks (Horvath, GrimAge, and Granulosa Cell) were used to predict DNA-methylation-based epigenetic age. To calculate the age acceleration, we regressed epigenetic age on chronological age and extracted the residuals. The association between epigenetic age acceleration and ovarian response outcomes (peak estradiol levels, follicle stimulation hormone, number of total, and mature oocytes) was assessed using linear and Poisson regression adjusted for chronological age, three surrogate variables (to account for cellular heterogeneity), race, smoking status, initial infertility diagnosis, and stimulation protocol.

MAIN RESULTS AND ROLE OF CHANCE

Compared to the median chronological age of our participants (34 years), the Horvath clock predicted, on an average, a younger epigenetic age (median: 24.2 years) while the GrimAge (median: 38.6 years) and Granulosa Cell (median: 39.0 years) clocks predicted, on an average, an older epigenetic age. Age acceleration based on the Horvath clock was associated with lower peak estradiol levels (-819.4 unit decrease in peak estradiol levels per standard deviation increase; 95% CI: -1265.7, -373.1) and fewer total (% change in total oocytes retrieved per standard deviation increase: -21.8%; 95% CI: -37.1%, -2.8%) and mature oocytes retrieved (% change in mature oocytes retrieved per standard deviation increase: -23.8%; 95% CI: -39.9%, -3.4%). The age acceleration based on the two other epigenetic clocks was not associated with markers of ovarian response.

LIMITATIONS, REASONS FOR CAUTION

Our sample size was small and we did not specifically isolate granulosa cells from follicular fluid samples so our samples could have included mixed cell types.

WIDER IMPLICATIONS OF THE FINDINGS

Our results highlight that certain epigenetic clocks may be predictive of ovarian stimulation outcomes when applied to follicular fluid; however, the inconsistent findings for specific clocks across studies indicate a need for further research to better understand the clinical utility of epigenetic clocks to improve IVF treatment.

STUDY FUNDING/COMPETING INTEREST(S)

The study was supported by grants ES009718, ES022955, ES000002, and ES026648 from the National Institute of Environmental Health Sciences (NIEHS) and a pilot grant from the NIEHS-funded HERCULES Center at Emory University (P30 ES019776). RBH was supported by the Emory University NIH Training Grant (T32-ES012870).

TRIAL REGISTRATION NUMBER

N/A.

Deciphering the molecular regulatory of RAB32/GPRC5A axis in chronic obstructive pulmonary disease

BACKGROUND

Chronic obstructive pulmonary disease (COPD) is a significant public health problem characterized by persistent airflow limitation. Despite previous research into the pathogenesis of COPD, a comprehensive understanding of the cell-type-specific mechanisms in COPD remains lacking. Recent studies have implicated Rab GTPases in regulating chronic immune response and inflammation via multiple pathways. In this study, the molecular regulating mechanism of RAB32 in COPD was investigated by multiple bioinformatics mining and experimental verification.

METHODS

We collected lung tissue surgical specimens from Zhongshan Hospital, Fudan University, and RT-qPCR and western blotting were used to detect the expression of Rabs in COPD lung tissues. Four COPD microarray datasets from the Gene Expression Omnibus (GEO) were analyzed. COPD-related epithelial cell scRNA-seq data was obtained from the GSE173896 dataset. Weighted gene co-expression network analysis (WGCNA), mfuzz cluster, and Spearman correlation analysis were combined to obtain the regulatory network of RAB32 in COPD. The slingshot algorithm was used to identify the regulatory molecule, and the co-localization of RAB32 and GPRC5A was observed with immunofluorescence.

RESULTS

WGCNA identified 771 key module genes significantly associated with the occurrence of COPD, including five Rab genes. RAB32 was up-regulated in lung tissues from subjects with COPD as contrast to those without COPD on both mRNA and protein levels. Integrating the results of WGCNA, Mfuzz clusters, and Spearman analysis, nine potential interacting genes with RAB32 were identified. Among these genes, GPRC5A exhibited a similar molecular expression pattern to RAB32. Co-expression density analysis at the cell level demonstrated that the co-expression density of RAB32 and GPRC5A was higher in type I alveolar epithelial cells (AT1s) than in type II alveolar epithelial cells (AT2s). The immunofluorescence also confirmed the co-localization of RAB32 and GPRC5A, and the Pearson correlation analysis found the relationship between RAB32 and GPRC5A was significantly stronger in the COPD lungs (r = 0.65) compared to the non-COPD lungs (r = 0.33).

CONCLUSIONS

Our study marked endeavor to delineate the molecular regulatory axis of RAB32 in COPD by employing diverse methods and identifying GPRC5A as a potential interacting molecule with RAB32. These findings offered novel perspectives on the mechanism of COPD.

Challenges and perspectives in computational deconvolution of genomics data

Deciphering cell-type heterogeneity is crucial for systematically understanding tissue homeostasis and its dysregulation in diseases. Computational deconvolution is an efficient approach for estimating cell-type abundances from a variety of omics data. Despite substantial methodological progress in computational deconvolution in recent years, challenges are still outstanding. Here we enlist four important challenges related to computational deconvolution: the quality of the reference data, generation of ground truth data, limitations of computational methodologies, and benchmarking design and implementation. Finally, we make recommendations on reference data generation, new directions of computational methodologies, and strategies to promote rigorous benchmarking.

Dynamics and consequences of nutrition-related microbial dysbiosis in early life: study protocol of the VITERBI GUT project

Introduction

Early life under- and overnutrition (jointly termed malnutrition) is increasingly recognized as an important risk factor for adult obesity and metabolic syndrome, a diet-related cluster of conditions including high blood sugar, fat and cholesterol. Nevertheless, the exact factors linking early life malnutrition with metabolic syndrome remain poorly characterized. We hypothesize that the microbiota plays a crucial role in this trajectory and that the pathophysiological mechanisms underlying under- and overnutrition are, to some extent, shared. We further hypothesize that a "dysbiotic seed microbiota" is transmitted to children during the birth process, altering the children's microbiota composition and metabolic health. The overall objective of this project is to understand the precise causes and biological mechanisms linking prenatal or early life under- or overnutrition with the predisposition to develop overnutrition and/or metabolic disease in later life, as well as to investigate the possibility of a dysbiotic seed microbiota inheritance in the context of maternal malnutrition.

Methods/design

VITERBI GUT is a prospective birth cohort allowing to study the link between early life malnutrition, the microbiota and metabolic health. VITERBI GUT will include 100 undernourished, 100 normally nourished and 100 overnourished pregnant women living in Vientiane, Lao People's Democratic Republic (PDR). Women will be recruited during their third trimester of pregnancy and followed with their child until its second birthday. Anthropometric, clinical, metabolic and nutritional data are collected from both the mother and the child. The microbiota composition of maternal and child's fecal and oral samples as well as maternal vaginal and breast milk samples will be determined using amplicon and shotgun metagenomic sequencing. Epigenetic modifications and lipid profiles will be assessed in the child's blood at 2 years of age. We will investigate for possible associations between metabolic health, epigenetics, and microbial changes.

Discussion

We expect the VITERBI GUT project to contribute to the emerging literature linking the early life microbiota, epigenetic changes and growth/metabolic health. We also expect this project to give new (molecular) insights into the mechanisms linking malnutrition-induced early life dysbiosis and metabolic health in later life, opening new avenues for microbiota-engineering using microbiota-targeted interventions.

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

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

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Construction and Verification of a Glycolysis-Associated Gene Signature for the Prediction of Overall Survival in Low Grade Glioma

The overall survival of patients with lower grade glioma (LGG) that might develop into high-grade malignant glioma shows marked heterogeneity. The currently used clinical evaluation index is not sufficient to predict precise prognostic outcomes accurately. To optimize survival risk stratification and the personalized management of patients with LGG, there is an urgent need to develop an accurate risk prediction model. The TCGA-LGG dataset, downloaded from The Cancer Genome Atlas (TCGA) portal, was used as a training cohort, and the Chinese Glioma Genome Atlas (CGGA) dataset and Rembrandt dataset were used as validation cohorts. The levels of various cancer hallmarks were quantified, which identified glycolysis as the primary overall survival-related risk factor in LGGs. Furthermore, using various bioinformatic and statistical methods, we developed a strong glycolysis-related gene signature to predict prognosis. Gene set enrichment analysis showed that in our model, high-risk glioma correlated with the chemoradiotherapy resistance and poor survival. Moreover, based on established risk model and other clinical features, a decision tree and a nomogram were built, which could serve as useful tools in the diagnosis and treatment of LGGs. This study indicates that the glycolysis-related gene signature could distinguish high-risk and low‐risk patients precisely, and thus can be used as an independent clinical feature.

What can we learn from DNA methylation studies in lupus?

During the past twenty years, a wide range of studies have established the existence of epigenetic alterations, particularly DNA methylation changes, in lupus. Epigenetic changes might have different contributions in children-onset versus adult-onset lupus. DNA methylation alterations have been identified and characterized in relation to disease activity and damage, different lupus subtypes and responses to drugs. However, to date there has been no practical application of these findings in the clinical milieu. In this article, we provide a review of key studies showing the relationship between DNA methylation and the many clinical aspects related to lupus. We also propose several options, in relation to the range of methodological developments and experimental design, that could optimize these findings and make them amenable for use in clinical practice.

A transdisciplinary approach to understand the epigenetic basis of race/ethnicity health disparities

Health disparities correspond to differences in disease burden and mortality among socially defined population groups. Such disparities may emerge according to race/ethnicity, socioeconomic status and a variety of other social contexts, and are documented for a wide range of diseases. Here, we provide a transdisciplinary perspective on the contribution of epigenetics to the understanding of health disparities, with a special emphasis on disparities across socially defined racial/ethnic groups. Scientists in the fields of biological anthropology, bioinformatics and molecular epidemiology provide a summary of theoretical, statistical and practical considerations for conducting epigenetic health disparities research, and provide examples of successful applications from cancer research using this approach.

Epigenome-wide analysis of DNA methylation and coronary heart disease: a nested case-control study

Background

Identifying environmentally responsive genetic loci where DNA methylation is associated with coronary heart disease (CHD) may reveal novel pathways or therapeutic targets for CHD. We conducted the first prospective epigenome-wide analysis of DNA methylation in relation to incident CHD in the Asian population.

Methods

We did a nested case-control study comprising incident CHD cases and 1:1 matched controls who were identified from the 10 year follow-up of the China Kadoorie Biobank. Methylation level of baseline blood leukocyte DNA was measured by Infinium Methylation EPIC BeadChip. We performed the single cytosine-phosphate-guanine (CpG) site association analysis and network approach to identify CHD-associated CpG sites and co-methylation gene module.

Results

After quality control, 982 participants (mean age 50.1 years) were retained. Methylation level at 25 CpG sites across the genome was associated with incident CHD (genome-wide false discovery rate [FDR] < 0.05 or module-specific FDR < 0.01). One SD increase in methylation level of identified CpGs was associated with differences in CHD risk, ranging from a 47 % decrease to a 118 % increase. Mediation analyses revealed 28.5 % of the excessed CHD risk associated with smoking was mediated by methylation level at the promoter region of ANKS1A gene (P for mediation effect = 0.036). Methylation level at the promoter region of SNX30 was associated with blood pressure and subsequent risk of CHD, with the mediating proportion to be 7.7 % (P = 0.003) via systolic blood pressure and 6.4 % (P = 0.006) via diastolic blood pressure. Network analysis revealed a co-methylation module associated with CHD.

Conclusions

We identified novel blood methylation alterations associated with incident CHD in the Asian population and provided evidence of the possible role of epigenetic regulations in the smoking- and blood pressure-related pathways to CHD risk.

Funding

This work was supported by National Natural Science Foundation of China (81390544 and 91846303). The CKB baseline survey and the first re-survey were supported by a grant from the Kadoorie Charitable Foundation in Hong Kong. The long-term follow-up is supported by grants from the UK Wellcome Trust (202922/Z/16/Z, 088158/Z/09/Z, 104085/Z/14/Z), grant (2016YFC0900500, 2016YFC0900501, 2016YFC0900504, 2016YFC1303904) from the National Key R&D Program of China, and Chinese Ministry of Science and Technology (2011BAI09B01).

Longitudinal analysis of individual cfDNA methylome patterns in metastatic prostate cancer

Background

Disease progression and therapeutic resistance are hallmarks of advanced stage prostate cancer (PCa), which remains a major cause of cancer-related mortality around the world. Longitudinal studies, coupled with the use of liquid biopsies, offer a potentially new and minimally invasive platform to study the dynamics of tumour progression. Our aim was to investigate the dynamics of personal DNA methylomic profiles of metastatic PCa (mPCa) patients, during disease progression and therapy administration.

Results

Forty-eight plasma samples from 9 mPCa patients were collected, longitudinally, over 13–21 months. After circulating cell-free DNA (cfDNA) isolation, DNA methylation was profiled using the Infinium MethylationEPIC BeadChip. The top 5% most variable probes across time, within each individual, were utilised to study dynamic methylation patterns during disease progression and therapeutic response. Statistical testing was carried out to identify differentially methylated genes (DMGs) in cfDNA, which were subsequently validated in two independent mPCa (cfDNA and FFPE tissue) cohorts. Individual cfDNA global methylation patterns were temporally stable throughout the disease course. However, a proportion of CpG sites presented a dynamic temporal pattern that was consistent with clinical events, including different therapies, and were prominently associated with genes linked to immune response pathways. Additionally, study of the tumour fraction of cfDNA identified > 2000 DMGs with dynamic methylation patterns.

Conclusions

Longitudinal assessment of cfDNA methylation in mPCa patients unveiled dynamic patterns associated with disease progression and therapy administration, thus highlighting the potential of using liquid biopsies to study PCa evolution at a methylomic level.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-021-01155-w.

Epigenetic inactivation of the autophagy-lysosomal system in appendix in Parkinson's disease

The gastrointestinal tract may be a site of origin for α-synuclein pathology in idiopathic Parkinson's disease (PD). Disruption of the autophagy-lysosome pathway (ALP) may contribute to α-synuclein aggregation. Here we examined epigenetic alterations in the ALP in the appendix by deep sequencing DNA methylation at 521 ALP genes. We identified aberrant methylation at 928 cytosines affecting 326 ALP genes in the appendix of individuals with PD and widespread hypermethylation that is also seen in the brain of individuals with PD. In mice, we find that DNA methylation changes at ALP genes induced by chronic gut inflammation are greatly exacerbated by α-synuclein pathology. DNA methylation changes at ALP genes induced by synucleinopathy are associated with the ALP abnormalities observed in the appendix of individuals with PD specifically involving lysosomal genes. Our work identifies epigenetic dysregulation of the ALP which may suggest a potential mechanism for accumulation of α-synuclein pathology in idiopathic PD.

Early-pregnancy maternal body mass index is associated with common DNA methylation markers in cord blood and placenta: a paired-tissue epigenome-wide association study

Women entering pregnancy with elevated body mass index (BMI) face greater risk of adverse outcomes during pregnancy, delivery, and for their offspring later in life, potentially via epigenetics. If epigenetic programming occurs early during in utero development, the differential marks should be detectable in multiple tissues despite the known unique epigenetic profile in each.We used early-pregnancy BMI as reflection of maternal metabolic milieu exposure in peri-conception and early-pregnancy period. We analysed DNA methylation in paired cord blood and placenta samples among 437 newborns from Gen3G, a pre-birth prospective cohort of primarily European descent. We measured DNA methylation in both tissues across the genome in >720,000 CpG sites using the Illumina MethylationEPIC array. At each site, we used linear mixed models (LMMs) with an unstructured variance-covariance matrix to test for an association between maternal early-pregnancy BMI and DNA methylation in both tissues (modelled as M-values). We adjusted for tissue-specific covariates, offspring sex, gestational age at delivery, and maternal smoking and age.Women had a mean (SD) BMI of 25.4 (5.7) kg/m² measured at first trimester visit (mean=9.9 weeks). Early-pregnancy BMI was associated with differential DNA methylation levels in paired-tissue analyses at two sites: cg10593758 (β=0.0126, SE=0.0025; P=4.07e-7), annotated to CRHBP, and cg0762168 (β=-0.0094, SE=0.0018; P=2.78e-7), annotated to CCDC97.Application of LMMs in DNA methylation data from distinct fetal-origin tissues allowed us to identify CpG sites at which early-pregnancy BMI may have an epigenetic 'programming' effect on overall fetus development. One site (CRHBP) may play a role in hypothalamic-pituitary-adrenal axis regulation.

A HML6 endogenous retrovirus on chromosome 3 is upregulated in amyotrophic lateral sclerosis motor cortex

There is increasing evidence that endogenous retroviruses (ERVs) play a significant role in central nervous system diseases, including amyotrophic lateral sclerosis (ALS). Studies of ALS have consistently identified retroviral enzyme reverse transcriptase activity in patients. Evidence indicates that ERVs are the cause of reverse transcriptase activity in ALS, but it is currently unclear whether this is due to a specific ERV locus or a family of ERVs. We employed a combination of bioinformatic methods to identify whether specific ERVs or ERV families are associated with ALS. Using the largest post-mortem RNA-sequence datasets available we selectively identified ERVs that closely resembled full-length proviruses. In the discovery dataset there was one ERV locus (HML6_3p21.31c) that showed significant increased expression in post-mortem motor cortex tissue after multiple-testing correction. Using six replication post-mortem datasets we found HML6_3p21.31c was consistently upregulated in ALS in motor cortex and cerebellum tissue. In addition, HML6_3p21.31c showed significant co-expression with cytokine binding and genes involved in EBV, HTLV-1 and HIV type-1 infections. There were no significant differences in ERV family expression between ALS and controls. Our results support the hypothesis that specific ERV loci are involved in ALS pathology.

Long intergenic noncoding RNA 299 methylation in peripheral blood is a biomarker for triple-negative breast cancer

Aim: To identify DNA methylation biomarkers in peripheral blood samples from triple-negative breast cancer (TNBC) patients. Materials & methods: We conducted an epigenome-wide association study (EWAS): the most promising markers were identified in 233 TNBC case–control pairs (discovery set) and subsequently validated in an independent validation set (57 TNBC patients and 124 controls). Results: cg06588802 (LINC00299/ID2) showed a higher methylation in TNBC patients compared with controls (discovery set: 3% increase, p-value = 0.0009; validation set: 2% increase, p-value = 0.01). Consistent results at four neighboring methylation probes and the strong negative correlation (rho = -0.93) with LINC00299 expression add plausibility to this result. Conclusion: Hypermethylation of LINC00299 in peripheral blood may constitute a useful circulating biomarker for TNBC.