Critical evaluation of the Illumina MethylationEPIC BeadChip microarray for whole-genome DNA methylation profiling

Accelerated epigenetic aging and inflammatory/immunological profile (ipAGE) in patients with chronic kidney disease

Chronic kidney disease (CKD) is defined by a reduced estimated glomerular filtration rate (eGFR). This failure can be related to a phenotype of accelerated aging. In this work, we considered 76 patients with end-stage renal disease (ESRD) and 83 healthy controls. We concomitantly evaluated for the first time two measures that can be informative of the rate of aging, i.e., whole blood DNA methylation using the Illumina Infinium EPIC array and plasma levels of a selection of inflammatory/immunological proteins using multiplex immunoassays. First of all, we demonstrated accelerated aging in terms of the most common epigenetic age estimators in CKD patients. Moreover, we developed a new clock/predictor of age based on the inflammatory/immunological profile (ipAGE) and identified the inflammatory/immunological biomarkers differentially expressed between cases and controls. IpAGE appeared to be more sensitive than epigenetic clocks in quantifying the accelerated aging phenotype of ESRD patients. Interestingly, we did not find any correlation between the age acceleration evaluated according to the epigenetic clocks and ipAGE in either the ESRD group or the control group. On the whole, our data show a consistent accelerated aging phenotype in ESRD patients, which is better appreciated by quantifying the underlying inflammatory processes (inflammaging) by ipAGE than by using epigenetic clocks.

Risks and rewards of big-data in epigenomics research: an interview with Melanie Ehrlich

Melanie Ehrlich, PhD, is a professor in the Tulane Cancer Center, the Tulane Center for Medical Bioinformatics and Genomics and the Hayward Human Genetics Program at Tulane Medical School, New Orleans, LA. She obtained her PhD in molecular biology in 1971 from the State University of New York at Stony Brook and completed postdoctoral research at Albert Einstein College of Medicine in 1972. She has been working on various aspects of epigenetics, starting with DNA methylation, since 1973. Her group made many first findings about DNA methylation (see below). For example, in 1982 and 1983, in collaboration with Charles Gehrke at the University of Missouri, she was the first to report tissue-specific and cancer-specific differences in overall DNA methylation in humans. In 1985, Xian-Yang Zhang and Richard Wang in her lab discovered a class of human DNA sequences specifically hypomethylated in sperm. In 1998, her group was the first to describe extensive losses of DNA methylation in pericentromeric and centromeric DNA repeats in human cancer. Her lab's many publications on the prevalence of both DNA hypermethylation and hypomethylation in the same cancers brought needed balance to our understanding of the epigenetics of cancer and to its clinical implications [1]. Besides working on cancer epigenetics, her research group has helped elucidate cytogenetic and gene expression abnormalities in the immunodeficiency, centromeric and facial anomalies (ICF) syndrome, a rare recessive disease often caused by mutations in DNMT3B. Her group also studied the epigenetics and transcriptomics of facioscapulohumeral muscular dystrophy (FSHD), whose disease locus is a tandem 3.3-kb repeat at subtelomeric 4q (that happens to be hypomethylated in ICF DNA [2]). Her study of FSHD has taken her in the direction of muscle (skeletal muscle, heart and aorta) epigenetics [3-6]. Recently, she has led research that applies epigenetics much more rigorously than usual to the evaluation of genetic variants from genome-wide association studies (GWAS) of osteoporosis and obesity. In continued collaboration with Sriharsa Pradhan at New England Biolabs and Michelle Lacey at Tulane University, she has compared 5-hydroxymethylcytosine and 5-methylcytosine clustering in various human tissues [7] and is studying myoblast methylomes that they generated by a new high-resolution enzymatic technique (enzymatic methyl-seq).

Gestational Perfluoroalkyl Substance Exposure and DNA Methylation at Birth and 12 Years of Age: A Longitudinal Epigenome-Wide Association Study

BACKGROUND

DNA methylation alterations may underlie associations between gestational perfluoroalkyl substances (PFAS) exposure and later-life health outcomes. To the best of our knowledge, no longitudinal studies have examined the associations between gestational PFAS and DNA methylation.

OBJECTIVES

We examined associations of gestational PFAS exposure with longitudinal DNA methylation measures at birth and in adolescence using the Health Outcomes and Measures of the Environment (HOME) Study (2003-2006; Cincinnati, Ohio).

METHODS

We quantified serum concentrations of perfluorooctanoate (PFOA), perfluorooctane sulfonate (PFOS), perfluorononanoate (PFNA), and perfluorohexane sulfonate (PFHxS) in mothers during pregnancy. We measured DNA methylation in cord blood (n=266) and peripheral leukocytes at 12 years of age (n=160) using the Illumina HumanMethylation EPIC BeadChip. We analyzed associations between log2-transformed PFAS concentrations and repeated DNA methylation measures using linear regression with generalized estimating equations. We included interaction terms between children's age and gestational PFAS. We performed Gene Ontology enrichment analysis to identify molecular pathways. We used Project Viva (1999-2002; Boston, Massachusetts) to replicate significant associations.

RESULTS

After adjusting for covariates, 435 cytosine-guanine dinucleotide (CpG) sites were associated with PFAS (false discovery rate, q<0.05). Specifically, we identified 2 CpGs for PFOS, 12 for PFOA, 8 for PFHxS, and 413 for PFNA; none overlapped. Among these, 2 CpGs for PFOA and 4 for PFNA were replicated in Project Viva. Some of the PFAS-associated CpG sites annotated to gene regions related to cancers, cognitive health, cardiovascular disease, and kidney function. We found little evidence that the associations between PFAS and DNA methylation differed by children's age.

DISCUSSION

In these longitudinal data, PFAS biomarkers were associated with differences in several CpGs at birth and at 12 years of age in or near genes linked to some PFAS-associated health outcomes. Future studies should examine whether DNA methylation mediates associations between gestational PFAS exposure and health. https://doi.org/10.1289/EHP10118.

Maternal Glycemic Dysregulation During Pregnancy and Neonatal Blood DNA Methylation: Meta-analyses of Epigenome-Wide Association Studies

OBJECTIVE

Maternal glycemic dysregulation during pregnancy increases the risk of adverse health outcomes in her offspring, a risk thought to be linearly related to maternal hyperglycemia. It is hypothesized that changes in offspring DNA methylation (DNAm) underline these associations.

RESEARCH DESIGN AND METHODS

To address this hypothesis, we conducted fixed-effects meta-analyses of epigenome-wide association study (EWAS) results from eight birth cohorts investigating relationships between cord blood DNAm and fetal exposure to maternal glucose (Nmaximum = 3,503), insulin (Nmaximum = 2,062), and area under the curve of glucose (AUCgluc) following oral glucose tolerance tests (Nmaximum = 1,505). We performed lookup analyses for identified cytosine-guanine dinucleotides (CpGs) in independent observational cohorts to examine associations between DNAm and cardiometabolic traits as well as tissue-specific gene expression.

RESULTS

Greater maternal AUCgluc was associated with lower cord blood DNAm at neighboring CpGs cg26974062 (β [SE] -0.013 [2.1 × 10-3], P value corrected for false discovery rate [PFDR] = 5.1 × 10-3) and cg02988288 (β [SE]-0.013 [2.3 × 10-3], PFDR = 0.031) in TXNIP. These associations were attenuated in women with GDM. Lower blood DNAm at these two CpGs near TXNIP was associated with multiple metabolic traits later in life, including type 2 diabetes. TXNIP DNAm in liver biopsies was associated with hepatic expression of TXNIP. We observed little evidence of associations between either maternal glucose or insulin and cord blood DNAm.

CONCLUSIONS

Maternal hyperglycemia, as reflected by AUCgluc, was associated with lower cord blood DNAm at TXNIP. Associations between DNAm at these CpGs and metabolic traits in subsequent lookup analyses suggest that these may be candidate loci to investigate in future causal and mediation analyses.

Characterization of Early Peripheral Immune Responses in Patients with Sepsis and Septic Shock

(1) Background: Sepsis is a life-threatening condition caused by an abnormal host response to infection that produces altered physiological responses causing tissue damage and can result in organ dysfunction and, in some cases, death. Although sepsis is characterized by a malfunction of the immune system leading to an altered immune response and immunosuppression, the high complexity of the pathophysiology of sepsis requires further investigation to characterize the immune response in sepsis and septic shock. (2) Methods: This study analyzes the immune-related responses occurring during the early stages of sepsis by comparing the amounts of cytokines, immune modulators and other endothelial mediators of a control group and three types of severe patients: critically ill non-septic patients, septic and septic shock patients. (3) Results: We showed that in the early stages of sepsis the innate immune system attempts to counteract infection, probably via neutrophils. Conversely, the adaptive immune system is not yet fully activated, either in septic or in septic shock patients. In addition, immunosuppressive responses and pro-coagulation signals are active in patients with septic shock. (4) Conclusions: The highest levels of IL-6 and pyroptosis-related cytokines (IL-18 and IL-1α) were found in septic shock patients, which correlated with D-dimer. Moreover, endothelial function may be affected as shown by the overexpression of adhesion molecules such as s-ICAM1 and E-Selectin during septic shock.

How to Translate DNA Methylation Biomarkers Into Clinical Practice

Recent advances in sequencing technologies provide unprecedented opportunities for epigenetic biomarker development. Particularly the DNA methylation pattern-which is modified at specific sites in the genome during cellular differentiation, aging, and disease-holds high hopes for a wide variety of diagnostic applications. While many epigenetic biomarkers have been described, only very few of them have so far been successfully translated into clinical practice and almost exclusively in the field of oncology. This discrepancy might be attributed to the different demands of either publishing a new finding or establishing a standardized and approved diagnostic procedure. This is exemplified for epigenetic leukocyte counts and epigenetic age-predictions. To ease later clinical translation, the following hallmarks should already be taken into consideration when designing epigenetic biomarkers: 1) Identification of best genomic regions, 2) pre-analytical processing, 3) accuracy of DNA methylation measurements, 4) identification of confounding parameters, 5) accreditation as diagnostic procedure, 6) standardized data analysis, 7) turnaround time, and 8) costs and customer requirements. While the initial selection of relevant genomic regions is usually performed on genome wide DNA methylation profiles, it might be advantageous to subsequently establish targeted assays that focus on specific genomic regions. Development of an epigenetic biomarker for clinical application is a long and cumbersome process that is only initiated with the identification of an epigenetic signature.

Integrating DNA Methylation Measures of Biological Aging into Social Determinants of Health Research

PURPOSE OF REVIEW

Acceleration of biological processes of aging is hypothesized to drive excess morbidity and mortality in socially disadvantaged populations. DNA methylation measures of biological aging provide tools for testing this hypothesis.

RECENT FINDINGS

Next-generation DNA methylation measures of biological aging developed to predict mortality risk and physiological decline are more predictive of morbidity and mortality than the original epigenetic clocks developed to predict chronological age. These new measures show consistent evidence of more advanced and faster biological aging in people exposed to socioeconomic disadvantage and may be able to record the emergence of socially determined health inequalities as early as childhood. Next-generation DNA methylation measures of biological aging also indicate race/ethnic disparities in biological aging. More research is needed on these measures in samples of non-Western and non-White populations. New DNA methylation measures of biological aging open opportunities for refining inference about the causes of social disparities in health and devising policies to eliminate them. Further refining measures of biological aging by including more diversity in samples used for measurement development is a critical priority for the field.

Gender-affirming hormone therapy induces specific DNA methylation changes in blood

Background

DNA methylation is an epigenetic mark that is influenced by underlying genetic profile, environment, and ageing. In addition to X-linked DNA methylation, sex-specific methylation patterns are widespread across autosomal chromosomes and can be present from birth or arise over time. In individuals where gender identity and sex assigned at birth are markedly incongruent, as in the case of transgender people, feminization or masculinization may be sought through gender-affirming hormone therapy (GAHT). GAHT is a cornerstone of transgender care, yet no studies to date have investigated its effect on genome-wide methylation. We profiled genome-wide DNA methylation in blood of transgender women (n = 13) and transgender men (n = 13) before and during GAHT (6 months and 12 months into feminizing or masculinizing hormone therapy).

Results

We identified several thousand differentially methylated CpG sites (DMPs) (Δβ ≥ 0.02, unadjusted p value < 0.05) and several differentially methylated regions (DMRs) in both people undergoing feminizing and masculinizing GAHT, the vast majority of which were progressive changes over time. X chromosome and sex-specific autosomal DNA methylation patterns established in early development are largely refractory to change in association with GAHT, with only 3% affected (Δβ ≥ 0.02, unadjusted p value < 0.05). The small number of sex-specific DMPs that were affected by GAHT were those that become sex-specific during the lifetime, known as sex-and-age DMPs, including DMRs in PRR4 and VMP1 genes. The GAHT-induced changes at these sex-associated probes consistently demonstrated a shift towards the methylation signature of the GAHT-naïve opposite sex, and we observed enrichment of previously reported adolescence-associated methylation changes.

Conclusion

We provide evidence for GAHT inducing a unique blood methylation signature in transgender people. This study advances our understanding of the complex interplay between sex hormones, sex chromosomes, and DNA methylation in the context of immunity. We highlight the need to broaden the field of ‘sex-specific’ immunity beyond cisgender males and cisgender females, as transgender people on GAHT exhibit a unique molecular profile.

Supplementary Information

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

A mammalian methylation array for profiling methylation levels at conserved sequences

Infinium methylation arrays are not available for the vast majority of non-human mammals. Moreover, even if species-specific arrays were available, probe differences between them would confound cross-species comparisons. To address these challenges, we developed the mammalian methylation array, a single custom array that measures up to 36k CpGs per species that are well conserved across many mammalian species. We designed a set of probes that can tolerate specific cross-species mutations. We annotate the array in over 200 species and report CpG island status and chromatin states in select species. Calibration experiments demonstrate the high fidelity in humans, rats, and mice. The mammalian methylation array has several strengths: it applies to all mammalian species even those that have not yet been sequenced, it provides deep coverage of conserved cytosines facilitating the development of epigenetic biomarkers, and it increases the probability that biological insights gained in one species will translate to others.

The applications of DNA methylation as a biomarker in kidney transplantation: a systematic review

BACKGROUND

Although kidney transplantation improves patient survival and quality of life, long-term results are hampered by both immune- and non-immune-mediated complications. Current biomarkers of post-transplant complications, such as allograft rejection, chronic renal allograft dysfunction, and cutaneous squamous cell carcinoma, have a suboptimal predictive value. DNA methylation is an epigenetic modification that directly affects gene expression and plays an important role in processes such as ischemia/reperfusion injury, fibrosis, and alloreactive immune response. Novel techniques can quickly assess the DNA methylation status of multiple loci in different cell types, allowing a deep and interesting study of cells' activity and function. Therefore, DNA methylation has the potential to become an important biomarker for prediction and monitoring in kidney transplantation.

PURPOSE OF THE STUDY

The aim of this study was to evaluate the role of DNA methylation as a potential biomarker of graft survival and complications development in kidney transplantation. MATERIAL AND METHODS: A systematic review of several databases has been conducted. The Newcastle-Ottawa scale and the Jadad scale have been used to assess the risk of bias for observational and randomized studies, respectively.

RESULTS

Twenty articles reporting on DNA methylation as a biomarker for kidney transplantation were included, all using DNA methylation for prediction and monitoring. DNA methylation pattern alterations in cells isolated from different tissues, such as kidney biopsies, urine, and blood, have been associated with ischemia-reperfusion injury and chronic renal allograft dysfunction. These alterations occurred in different and specific loci. DNA methylation status has also proved to be important for immune response modulation, having a crucial role in regulatory T cell definition and activity. Research also focused on a better understanding of the role of this epigenetic modification assessment for regulatory T cells isolation and expansion for future tolerance induction-oriented therapies.

CONCLUSIONS

Studies included in this review are heterogeneous in study design, biological samples, and outcome. More coordinated investigations are needed to affirm DNA methylation as a clinically relevant biomarker important for prevention, monitoring, and intervention.

Placental genomics mediates genetic associations with complex health traits and disease

As the master regulator in utero, the placenta is core to the Developmental Origins of Health and Disease (DOHaD) hypothesis but is historically understudied. To identify placental gene-trait associations (GTAs) across the life course, we perform distal mediator-enriched transcriptome-wide association studies (TWAS) for 40 traits, integrating placental multi-omics from the Extremely Low Gestational Age Newborn Study. At [Formula: see text], we detect 248 GTAs, mostly for neonatal and metabolic traits, across 176 genes, enriched for cell growth and immunological pathways. In aggregate, genetic effects mediated by placental expression significantly explain 4 early-life traits but no later-in-life traits. 89 GTAs show significant mediation through distal genetic variants, identifying hypotheses for distal regulation of GTAs. Investigation of one hypothesis in human placenta-derived choriocarcinoma cells reveal that knockdown of mediator gene EPS15 upregulates predicted targets SPATA13 and FAM214A, both associated with waist-hip ratio in TWAS, and multiple genes involved in metabolic pathways. These results suggest profound health impacts of placental genomic regulation in developmental programming across the life course.

Does DNA methylation mediate the association of age at puberty with FVC or FEV1?

Background

Age of pubertal onset is associated with lung function in adulthood. However, the underlying role of epigenetics as a mediator of this association remains unknown.

Methods

DNA methylation (DNAm) in peripheral blood was measured at age 18 years in the Isle of Wight birth cohort (IOWBC) along with data on age of pubertal events, forced vital capacity (FVC) and forced expiratory volume in 1 s (FEV₁) at 26 years. Structural equation models were applied to examine mediation effects of DNAm on the association of age at pubertal events with FVC and FEV₁. Findings were further tested in the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort.

Results

In the IOWBC, for females, 21 cytosine-phosphate-guanine sites (CpGs) were shown to mediate the association of age at puberty with FVC or FEV₁ at 26 years (p<0.05). In males, DNAm at 20 CpGs was found to mediate the association of age at puberty with FVC (p<0.05). At almost all these CpGs, indirect effects (effects of age at pubertal events on FVC or FEV₁via DNAm) contributed a smaller portion to the total effects compared to direct effects (e.g. at cg08680129, ∼22% of the estimated total effect of age at menarche on FVC at age 26 was contributed by an indirect effect). Among the IOWBC-discovered CpGs available in ALSPAC, none of them was replicated in ALSPAC (p>0.05).

Conclusions

Our findings suggest that post-adolescence DNAm in peripheral blood is likely not to mediate the association of age at pubertal onset with young adulthood FVC or FEV₁.

The association between age at pubertal onset and lung function parameters FVC or FEV₁ in young adulthood is not likely to be mediated by DNA methylation in peripheral bloodhttps://bit.ly/31G8hDi

The WID-BC-index identifies women with primary poor prognostic breast cancer based on DNA methylation in cervical samples

Genetic and non-genetic factors contribute to breast cancer development. An epigenome-based signature capturing these components in easily accessible samples could identify women at risk. Here, we analyse the DNA methylome in 2,818 cervical, 357 and 227 matched buccal and blood samples respectively, and 42 breast tissue samples from women with and without breast cancer. Utilising cervical liquid-based cytology samples, we develop the DNA methylation-based Women’s risk IDentification for Breast Cancer index (WID-BC-index) that identifies women with breast cancer with an AUROC (Area Under the Receiver Operator Characteristic) of 0.84 (95% CI: 0.80–0.88) and 0.81 (95% CI: 0.76–0.86) in internal and external validation sets, respectively. CpGs at progesterone receptor binding sites hypomethylated in normal breast tissue of women with breast cancer or in BRCA mutation carriers are also hypomethylated in cervical samples of women with poor prognostic breast cancer. Our data indicate that a systemic epigenetic programming defect is highly prevalent in women who develop breast cancer. Further studies validating the WID-BC-index may enable clinical implementation for monitoring breast cancer risk.

Breast cancer is most commonly diagnosed via a needle biopsy. In this study, the authors show that cervical samples from women with breast cancer have a methylation signature different to that of healthy controls.

The DNA methylome of cervical cells can predict the presence of ovarian cancer

The vast majority of epithelial ovarian cancer arises from tissues that are embryologically derived from the Müllerian Duct. Here, we demonstrate that a DNA methylation signature in easy-to-access Müllerian Duct-derived cervical cells from women with and without ovarian cancer (i.e. referred to as the Women's risk IDentification for Ovarian Cancer index or WID-OC-index) is capable of identifying women with an ovarian cancer in the absence of tumour DNA with an AUC of 0.76 and women with an endometrial cancer with an AUC of 0.81. This and the observation that the cervical cell WID-OC-index mimics the epigenetic program of those cells at risk of becoming cancerous in BRCA1/2 germline mutation carriers (i.e. mammary epithelium, fallopian tube fimbriae, prostate) further suggest that the epigenetic misprogramming of cervical cells is an indicator for cancer predisposition. This concept has the potential to advance the field of risk-stratified cancer screening and prevention.

DNA methylation changes in glial cells of the normal-appearing white matter in Multiple Sclerosis patients

Multiple Sclerosis (MS), the leading cause of non-traumatic neurological disability in young adults, is a chronic inflammatory and neurodegenerative disease of the central nervous system (CNS). Due to the poor accessibility to the target organ, CNS-confined processes underpinning the later progressive form of MS remain elusive thereby limiting treatment options. We aimed to examine DNA methylation, a stable epigenetic mark of genome activity, in glial cells to capture relevant molecular changes underlying MS neuropathology. We profiled DNA methylation in nuclei of non-neuronal cells, isolated from 38 post-mortem normal-appearing white matter (NAWM) specimens of MS patients (n = 8) in comparison to white matter of control individuals (n = 14), using Infinium MethylationEPIC BeadChip. We identified 1,226 significant (genome-wide adjusted P-value < 0.05) differentially methylated positions (DMPs) between MS patients and controls. Functional annotation of the altered DMP-genes uncovered alterations of processes related to cellular motility, cytoskeleton dynamics, metabolic processes, synaptic support, neuroinflammation and signaling, such as Wnt and TGF-β pathways. A fraction of the affected genes displayed transcriptional differences in the brain of MS patients, as reported by publically available transcriptomic data. Cell type-restricted annotation of DMP-genes attributed alterations of cytoskeleton rearrangement and extracellular matrix remodelling to all glial cell types, while some processes, including ion transport, Wnt/TGF-β signaling and immune processes were more specifically linked to oligodendrocytes, astrocytes and microglial cells, respectively. Our findings strongly suggest that NAWM glial cells are highly altered, even in the absence of lesional insult, collectively exhibiting a multicellular reaction in response to diffuse inflammation.

An in vitro chronic damage model impairs inflammatory and regenerative responses in human colonoid monolayers

Acute damage to the intestinal epithelium can be repaired via de-differentiation of mature intestinal epithelial cells (IECs) to a stem-like state, but there is a lack of knowledge on how intestinal stem cells function after chronic injury, such as in inflammatory bowel disease (IBD). We developed a chronic-injury model in human colonoid monolayers by repeated rounds of air-liquid interface and submerged culture. We use this model to understand how chronic intestinal damage affects the ability of IECs to (1) respond to microbial stimulation, using the Toll-like receptor 5 (TLR5) agonist FliC and (2) regenerate and protect the epithelium from further damage. Repeated rounds of damage impair the ability of IECs to regrow and respond to TLR stimulation. We also identify mRNA expression and DNA methylation changes in genes associated with IBD and colon cancer. This methodology results in a human model of recurrent IEC injury like that which occurs in IBD.

Novel diagnostic DNA methylation episignatures expand and refine the epigenetic landscapes of Mendelian disorders

Overlapping clinical phenotypes and an expanding breadth and complexity of genomic associations are a growing challenge in the diagnosis and clinical management of Mendelian disorders. The functional consequences and clinical impacts of genomic variation may involve unique, disorder-specific, genomic DNA methylation episignatures. In this study, we describe 19 novel episignature disorders and compare the findings alongside 38 previously established episignatures for a total of 57 episignatures associated with 65 genetic syndromes. We demonstrate increasing resolution and specificity ranging from protein complex, gene, sub-gene, protein domain, and even single nucleotide-level Mendelian episignatures. We show the power of multiclass modeling to develop highly accurate and disease-specific diagnostic classifiers. This study significantly expands the number and spectrum of disorders with detectable DNA methylation episignatures, improves the clinical diagnostic capabilities through the resolution of unsolved cases and the reclassification of variants of unknown clinical significance, and provides further insight into the molecular etiology of Mendelian conditions.

CpG Island Methylator Phenotype-A Hope for the Future or a Road to Nowhere?

The CpG island methylator phenotype (CIMP) can be regarded as the most notable emanation of epigenetic instability in cancer. Since its discovery in the late 1990s, CIMP has been extensively studied, mainly in colorectal cancers (CRC) and gliomas. Consequently, knowledge on molecular and pathological characteristics of CIMP in CRC and other tumour types has rapidly expanded. Concordant and widespread hypermethylation of multiple CpG islands observed in CIMP in multiple cancers raised hopes for future epigenetically based diagnostics and treatments of solid tumours. However, studies on CIMP in solid tumours were hampered by a lack of generalisability and reproducibility of epigenetic markers. Moreover, CIMP was not a satisfactory marker in predicting clinical outcomes. The idea of targeting epigenetic abnormalities such as CIMP for cancer therapy has not been implemented for solid tumours, either. Twenty-one years after its discovery, we aim to cover both the fundamental and new aspects of CIMP and its future application as a diagnostic marker and target in anticancer therapies.

Assessment of the Molecular Heterogeneity of E-Cadherin Expression in Invasive Lobular Breast Cancer

Simple Summary

Invasive lobular breast cancers (ILCs) are histologically classified by their discohesive growth pattern, due to loss of the cell adhesion glycoprotein E-cadherin (CDH1), which arises via mutation in CDH1 in around half of these tumours. A subset of these tumours, however, show mixed levels of E-cadherin expression. Here, we sought to address whether the distinct parts of individual tumours showing heterogeneous E-cadherin expression harbour distinct driver alterations. Using whole genome sequencing and methylation profiling of nine such cases, we identified that these tumours are clonally related, suggesting that they are part of the spectrum of ILC tumours. CDH1 mutant tumours showed a higher mutational burden indicative of APOBEC-mediated mutagenesis. In some cases, known clinically actionable driver mutations, such as PIK3CA, were exclusive to one component. Together, these results highlight the heterogeneity underpinning this special histological breast cancer.

Abstract

Mutations and loss of E-cadherin protein expression define the vast majority of invasive lobular carcinomas. In a subset of these cases, the heterogeneous expression of E-cadherin is observed either as wild-type (strong membranous) expression or aberrant expression (cytoplasmic expression). However, it is unclear as to whether the two components would be driven by distinct genetic or epigenetic alterations. Here, we used whole genome DNA sequencing and methylation array profiling of two separately dissected components of nine invasive lobular carcinomas with heterogeneous E-cadherin expression. E-cadherin negative and aberrant/positive components of E-cadherin heterogeneous tumours showed a similar mutational, copy number and promoter methylation repertoire, suggesting they arise from a common ancestor, as opposed to the collision of two independent tumours. We found that the majority of E-cadherin heterogeneous tumours harboured CDH1 mutations in both the E-cadherin negative and aberrant/positive components together with somatic mutations in additional driver genes known to be enriched in both pure invasive carcinomas of no special type and invasive lobular breast cancers, whereas these were less commonly observed in CDH1 wild-type tumours. CDH1 mutant tumours also exhibited a higher mutation burden as well as increased presence of APOBEC-dependent mutational signatures 2 and 13 compared to CDH1 wild-type tumours. Together, our results suggest that regardless of E-cadherin protein expression, tumours showing heterogeneous expression of E-cadherin should be considered as part of the spectrum of invasive lobular breast cancers.

Genome-Wide Placental Gene Methylations in Gestational Diabetes Mellitus, Fetal Growth and Metabolic Health Biomarkers in Cord Blood

Gestational diabetes mellitus (GDM) "program" an elevated risk of metabolic syndrome in the offspring. Epigenetic alterations are a suspected mechanism. GDM has been associated with placental DNA methylation changes in some epigenome-wide association studies. It remains unclear which genes or pathways are affected, and whether any placental differential gene methylations are correlated to fetal growth or circulating metabolic health biomarkers. In an epigenome-wide association study using the Infinium MethylationEPIC Beadchip, we sought to identify genome-wide placental differentially methylated genes and enriched pathways in GDM, and to assess the correlations with fetal growth and metabolic health biomarkers in cord blood. The study samples were 30 pairs of term placentas in GDM vs. euglycemic pregnancies (controls) matched by infant sex and gestational age at delivery in the Shanghai Birth Cohort. Cord blood metabolic health biomarkers included insulin, C-peptide, proinsulin, IGF-I, IGF-II, leptin and adiponectin. Adjusting for maternal age, pre-pregnancy BMI, parity, mode of delivery and placental cell type heterogeneity, 256 differentially methylated positions (DMPs,130 hypermethylated and 126 hypomethylated) were detected between GDM and control groups accounting for multiple tests with false discovery rate <0.05 and beta-value difference >0.05. WSCD2 was identified as a differentially methylated gene in both site- and region-level analyses. We validated 7 hypermethylated (CYP1A2, GFRA1, HDAC4, LIMS2, NAV3, PAX6, UPK1B) and 10 hypomethylated (DPP10, CPLX1, CSMD2, GPR133, NRXN1, PCSK9, PENK, PRDM16, PTPRN2, TNXB) genes reported in previous epigenome-wide association studies. We did not find any enriched pathway accounting for multiple tests. DMPs in 11 genes (CYP2D7P1, PCDHB15, ERG, SIRPB1, DKK2, RAPGEF5, CACNA2D4, PCSK9, TSNARE1, CADM2, KCNAB2) were correlated with birth weight (z score) accounting for multiple tests. There were no significant correlations between placental gene methylations and cord blood biomarkers. In conclusions, GDM was associated with DNA methylation changes in a number of placental genes, but these placental gene methylations were uncorrelated to the observed metabolic health biomarkers (fetal growth factors, leptin and adiponectin) in cord blood. We validated 17 differentially methylated placental genes in GDM, and identified 11 differentially methylated genes relevant to fetal growth.

Future Prospects for Epigenetics in Autism Spectrum Disorder

Despite decades of investigation into the genetics of autism spectrum disorder (ASD), a current consensus in the field persists that ASD risk is too heterogeneous to be diagnosed by a single set of genetic variants. As such, ASD research has broadened to include assessment of other molecular biomarkers implicated in the condition that may be reflective of environmental exposures or gene by environment interactions. Epigenetic variance, and specifically differential DNA methylation, have emerged as areas of particularly high interest to ASD, as the epigenetic markers from specific chromatin loci collectively can reflect influences of multiple genetic and environmental factors and can also result in differential gene expression patterns. This review examines recent studies of the ASD epigenome, detailing common gene pathways found to be differentially methylated in people with ASD, and considers how these discoveries may inform our understanding of ASD etiology. We also consider future applications of epigenetics in ASD research and clinical practice, focusing on substratification, biomarker development, and experimental preclinical models of ASD that test causality. In combination with other -omics approaches, epigenomics allows an improved conceptualization of the multifactorial nature of ASD, and opens future lines of inquiry for both basic research and clinical practice.

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

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

DNA Methylation Profiles of Ovarian Clear Cell Carcinoma

BACKGROUND

Ovarian clear cell carcinoma (OCCC) is a rare ovarian cancer histotype that tends to be resistant to standard platinum-based chemotherapeutics. We sought to better understand the role of DNA methylation in clinical and biological subclassification of OCCC.

METHODS

We interrogated genome-wide methylation using DNA from fresh frozen tumors from 271 cases, applied nonsmooth nonnegative matrix factorization (nsNMF) clustering, and evaluated clinical associations and biological pathways.

RESULTS

Two approximately equally sized clusters that associated with several clinical features were identified. Compared with Cluster 2 (N = 137), Cluster 1 cases (N = 134) presented at a more advanced stage, were less likely to be of Asian ancestry, and tended to have poorer outcomes including macroscopic residual disease following primary debulking surgery (P < 0.10). Subset analyses of targeted tumor sequencing and IHC data revealed that Cluster 1 tumors showed TP53 mutation and abnormal p53 expression, and Cluster 2 tumors showed aneuploidy and ARID1A/PIK3CA mutation (P < 0.05). Cluster-defining CpGs included 1,388 CpGs residing within 200 bp of the transcription start sites of 977 genes; 38% of these genes (N = 369 genes) were differentially expressed across cluster in transcriptomic subset analysis (P < 10-4). Differentially expressed genes were enriched for six immune-related pathways, including IFNα and IFNγ responses (P < 10-6).

CONCLUSIONS

DNA methylation clusters in OCCC correlate with disease features and gene expression patterns among immune pathways.

IMPACT

This work serves as a foundation for integrative analyses that better understand the complex biology of OCCC in an effort to improve potential for development of targeted therapeutics.

Culture-Associated DNA Methylation Changes Impact on Cellular Function of Human Intestinal Organoids

BACKGROUND & AIMS

Human intestinal epithelial organoids (IEOs) are a powerful tool to model major aspects of intestinal development, health, and diseases because patient-derived cultures retain many features found in vivo. A necessary aspect of the organoid model is the requirement to expand cultures in vitro through several rounds of passaging. This is of concern because the passaging of cells has been shown to affect cell morphology, ploidy, and function.

METHODS

Here, we analyzed 173 human IEO lines derived from the small and large bowel and examined the effect of culture duration on DNA methylation (DNAm). Furthermore, we tested the potential impact of DNAm changes on gene expression and cellular function.

RESULTS

Our analyses show a reproducible effect of culture duration on DNAm in a large discovery cohort as well as 2 publicly available validation cohorts generated in different laboratories. Although methylation changes were seen in only approximately 8% of tested cytosine-phosphate-guanine dinucleotides (CpGs) and global cellular function remained stable, a subset of methylation changes correlated with altered gene expression at baseline as well as in response to inflammatory cytokine exposure and withdrawal of Wnt agonists. Importantly, epigenetic changes were found to be enriched in genomic regions associated with colonic cancer and distant to the site of replication, indicating similarities to malignant transformation.

CONCLUSIONS

Our study shows distinct culture-associated epigenetic changes in mucosa-derived human IEOs, some of which appear to impact gene transcriptomic and cellular function. These findings highlight the need for future studies in this area and the importance of considering passage number as a potentially confounding factor.

Data Analysis of DNA Methylation Epigenome-Wide Association Studies (EWAS): A Guide to the Principles of Best Practice

Array-based EWAS have become an increasingly popular technique to identify population epigenetic effects, particularly in humans. With the arrival of nonhuman species arrays, such as the mouse, this is likely to become an even more widely used technology. This chapter provides the less experienced researcher a guide to the analysis of data from the most widely used platform, the Illumina Infinium Methylation assay. This includes an overview of quality filtering, data normalization, analysis options, and techniques to improve the interpretation of results.

DNA Methylation Profiles of Purified Cell Types in Bronchoalveolar Lavage: Applications for Mixed Cell Paediatric Pulmonary Studies

In epigenome-wide association studies analysing DNA methylation from samples containing multiple cell types, it is essential to adjust the analysis for cell type composition. One well established strategy for achieving this is reference-based cell type deconvolution, which relies on knowledge of the DNA methylation profiles of purified constituent cell types. These are then used to estimate the cell type proportions of each sample, which can then be incorporated to adjust the association analysis. Bronchoalveolar lavage is commonly used to sample the lung in clinical practice and contains a mixture of different cell types that can vary in proportion across samples, affecting the overall methylation profile. A current barrier to the use of bronchoalveolar lavage in DNA methylation-based research is the lack of reference DNA methylation profiles for each of the constituent cell types, thus making reference-based cell composition estimation difficult. Herein, we use bronchoalveolar lavage samples collected from children with cystic fibrosis to define DNA methylation profiles for the four most common and clinically relevant cell types: alveolar macrophages, granulocytes, lymphocytes and alveolar epithelial cells. We then demonstrate the use of these methylation profiles in conjunction with an established reference-based methylation deconvolution method to estimate the cell type composition of two different tissue types; a publicly available dataset derived from artificial blood-based cell mixtures and further bronchoalveolar lavage samples. The reference DNA methylation profiles developed in this work can be used for future reference-based cell type composition estimation of bronchoalveolar lavage. This will facilitate the use of this tissue in studies examining the role of DNA methylation in lung health and disease.

A Bayesian hierarchical model to estimate DNA methylation conservation in colorectal tumors

MOTIVATION

Conservation is broadly used to identify biologically important (epi)genomic regions. In the case of tumor growth, preferential conservation of DNA methylation can be used to identify areas of particular functional importance to the tumor. However, reliable assessment of methylation conservation based on multiple tissue samples per patient requires the decomposition of methylation variation at multiple levels.

RESULTS

We developed a Bayesian hierarchical model that allows for variance decomposition of methylation on three levels: between-patient normal tissue variation, between-patient tumor-effect variation and within-patient tumor variation. We then defined a model-based conservation score to identify loci of reduced within-tumor methylation variation relative to between-patient variation. We fit the model to multi-sample methylation array data from 21 colorectal cancer (CRC) patients using a Monte Carlo Markov Chain algorithm (Stan). Sets of genes implicated in CRC tumorigenesis exhibited preferential conservation, demonstrating the model's ability to identify functionally relevant genes based on methylation conservation. A pathway analysis of preferentially conserved genes implicated several CRC relevant pathways and pathways related to neoantigen presentation and immune evasion. Our findings suggest that preferential methylation conservation may be used to identify novel gene targets that are not consistently mutated in CRC. The flexible structure makes the model amenable to the analysis of more complex multi-sample data structures.

AVAILABILITY AND IMPLEMENTATION

The data underlying this article are available in the NCBI GEO Database, under accession code GSE166212. The R analysis code is available at https://github.com/kevin-murgas/DNAmethylation-hierarchicalmodel.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Reproductive history and blood cell DNA methylation later in life: the Young Finns Study

Background

Women with a history of complications of pregnancy, including hypertensive disorders, gestational diabetes or an infant fetal growth restriction or preterm birth, are at higher risk for cardiovascular disease later in life. We aimed to examine differences in maternal DNA methylation following pregnancy complications.

Methods

Data on women participating in the Young Finns study (n = 836) were linked to the national birth registry. DNA methylation in whole blood was assessed using the Infinium Methylation EPIC BeadChip. Epigenome-wide analysis was conducted on differential CpG methylation at 850 K sites. Reproductive history was also modeled as a predictor of four epigenetic age indices.

Results

Fourteen significant differentially methylated sites were found associated with both history of pre-eclampsia and overall hypertensive disorders of pregnancy. No associations were found between reproductive history and any epigenetic age acceleration measure.

Conclusions

Differences in epigenetic methylation profiles could represent pre-existing risk factors, or changes that occurred as a result of experiencing these complications.

Supplementary Information

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

Identification of DNA methylation biomarkers with potential to predict response to neoadjuvant chemotherapy in triple-negative breast cancer

Neoadjuvant chemotherapy (NAC) is used to treat triple-negative breast cancer (TNBC) prior to resection. Biomarkers that accurately predict a patient's response to NAC are needed to individualise therapy and avoid chemotoxicity from unnecessary chemotherapy. We performed whole-genome DNA methylation profiling on diagnostic TNBC biopsy samples from the Sequential Evaluation of Tumours Undergoing Preoperative (SETUP) NAC study. We found 9 significantly differentially methylated regions (DMRs) at diagnosis which were associated with response to NAC. We show that 4 of these DMRs are associated with TNBC overall survival (P < 0.05). Our results highlight the potential of DNA methylation biomarkers for predicting NAC response in TNBC.

An epigenetic aging analysis of randomized metformin and weight loss interventions in overweight postmenopausal breast cancer survivors

Metformin and weight loss relationships with epigenetic age measures-biological aging biomarkers-remain understudied. We performed a post-hoc analysis of a randomized controlled trial among overweight/obese breast cancer survivors (N = 192) assigned to metformin, placebo, weight loss with metformin, or weight loss with placebo interventions for 6 months. Epigenetic age was correlated with chronological age (r = 0.20-0.86; P < 0.005). However, no significant epigenetic aging associations were observed by intervention arms. Consistent with published reports in non-cancer patients, 6 months of metformin therapy may be inadequate to observe expected epigenetic age deceleration. Longer duration studies are needed to better characterize these relationships.Trial Registration: Registry Name: ClincialTrials.Gov.Registration Number: NCT01302379.Date of Registration: February 2011.URL: https://clinicaltrials.gov/ct2/show/NCT01302379.

Epigenetic Age Acceleration Is Not Associated with Age-Related Macular Degeneration

DNA methylation age (DNAm age) estimation is a powerful biomarker of human ageing. To date, epigenetic clocks have not been evaluated in age-related macular degeneration (AMD). Here, we perform genome-wide DNA methylation analyses in blood of AMD patients with a documented smoking history (14 AMD, 16 Normal), identifying loci of differential methylation (DML) with a relaxed p-value criterion (p ≤ 10⁻⁴). We conduct DNAm age analyses using the Horvath-multi tissue, Hannum and Skin & Blood epigenetic clocks in both blood and retinal pigment epithelium (RPE). We perform Ingenuity Pathway Analysis Causal Network Analysis (IPA CNA) on the topmost significantly differentially methylated CpG probes in blood and RPE. Results show poor performance of epigenetic clocks in RPE. Epigenetic age acceleration (EAA) was not observed in AMD. However, we observe positive EAA in blood of smokers, and in smokers with AMD. DML analysis revealed hypomethylation at cg04953735 within RPTOR (p = 6.51 × 10⁻⁵; Δβ = −11.95%). IPA CNA in the RPE also identified RPTOR as the putative master regulator, predicted to be inhibited in AMD. In conclusion, this is the first study evaluating an association of epigenetic ageing in AMD. We posit a role for RPTOR as a common master regulator of methylation changes in the RPE in AMD.

Reduced NCOR2 expression accelerates androgen deprivation therapy failure in prostate cancer

This study addresses the roles of nuclear receptor corepressor 2 (NCOR2) in prostate cancer (PC) progression in response to androgen deprivation therapy (ADT). Reduced NCOR2 expression significantly associates with shorter disease-free survival in patients with PC receiving adjuvant ADT. Utilizing the CWR22 xenograft model, we demonstrate that stably reduced NCOR2 expression accelerates disease recurrence following ADT, associates with gene expression patterns that include neuroendocrine features, and induces DNA hypermethylation. Stably reduced NCOR2 expression in isogenic LNCaP (androgen-sensitive) and LNCaP-C4–2 (androgen-independent) cells revealed that NCOR2 reduction phenocopies the impact of androgen treatment and induces global DNA hypermethylation patterns. NCOR2 genomic binding is greatest in LNCaP-C4–2 cells and most clearly associates with forkhead box (FOX) transcription factor FOXA1 binding. NCOR2 binding significantly associates with transcriptional regulation most when in active enhancer regions. These studies reveal robust roles for NCOR2 in regulating the PC transcriptome and epigenome and underscore recent mutational studies linking NCOR2 loss of function to PC disease progression.

Long et al. show that reduced levels of NCOR2 lead to accelerated prostate cancer recurrence during androgen withdrawal in a patient-derived xenograft model. NCOR2 reduction is characterized by incomplete response to androgen withdrawal, and recurrent tumors show increased neuroendocrine traits. These phenotypic changes are associated with hypermethylated enhancers.

Locus-Specific Methylation of GSTP1, RNF219, and KIAA1539 Genes with Single Molecule Resolution in Cell-Free DNA from Healthy Donors and Prostate Tumor Patients: Application in Diagnostics

Simple Summary

Prostate cancer (PCa) is the second most commonly diagnosed cancer in men, which is constantly accompanied by benign prostate hyperplasia (BPH). To reach a 100% 5-year survival rate in PCa, which is characteristic for PCa if it is diagnosed in early stages, efficient PCa diagnostics against the background of BPH are demanded. The article describes a liquid biopsy approach to differential PCa diagnostics based on the data on locus-specific methylation of the three genes (GSTP1, RNF219, and KIAA1539) obtained with NGS of cell-free DNA from blood plasma of PCa, BPH, and healthy individuals. We offered a diagnostic approach including the analysis of simultaneous methylation status of two CpGs in one cell-free DNA molecule, allowing the discrimination of all patients with absolute sensitivity and specificity.

Abstract

The locus-specific methylation of three genes (GSTP1, RNF219, and KIAA1539 (also known as FAM214B)) in the blood plasma cell-free DNA (cfDNA) of 20 patients with prostate cancer (PCa), 18 healthy donors (HDs), and 17 patients with benign prostatic hyperplasia (BPH) was studied via the MiSeq platform. The methylation status of two CpGs within the same loci were used as the diagnostic feature for discriminating the patient groups. Many variables had good diagnostic characteristics, e.g., each of the variables GSTP1.C3.C9, GSTP1.C9, and GSTP1.C9.T17 demonstrated an 80% sensitivity at a 100% specificity for PCa patients vs. the others comparison. The analysis of RNF219 gene loci methylation allowed discriminating BPH patients with absolute sensitivity and specificity. The data on the methylation of the genes GSTP1 and RNF219 allowed discriminating PCa patients, as well as HDs, with absolute sensitivity and specificity. Thus, the data on the locus-specific methylation of cfDNA (with single-molecule resolution) combined with a diagnostic approach considering the simultaneous methylation of several CpGs in one locus enabled the discrimination of HD, BPH, and PCa patients.

Transcriptomic Remodelling of Fetal Endothelial Cells During Establishment of Inflammatory Memory

Inflammatory memory involves the molecular and cellular ‘reprogramming’ of innate immune cells following exogenous stimuli, leading to non-specific protection against subsequent pathogen exposure. This phenomenon has now also been described in non-hematopoietic cells, such as human fetal and adult endothelial cells. In this study we mapped the cell-specific DNA methylation profile and the transcriptomic remodelling during the establishment of inflammatory memory in two distinct fetal endothelial cell types – a progenitor cell (ECFC) and a differentiated cell (HUVEC) population. We show that both cell types have a core transcriptional response to an initial exposure to a viral-like ligand, Poly(I:C), characterised by interferon responsive genes. There was also an ECFC specific response, marked by the transcription factor ELF1, suggesting a non-canonical viral response pathway in progenitor endothelial cells. Next, we show that both ECFCs and HUVECs establish memory in response to an initial viral exposure, resulting in an altered subsequent response to lipopolysaccharide. While the capacity to train or tolerize the induction of specific sets of genes was similar between the two cell types, the progenitor ECFCs show a higher capacity to establish memory. Among tolerized cellular pathways are those involved in endothelial barrier establishment and leukocyte migration, both important for regulating systemic immune-endothelial cell interactions. These findings suggest that the capacity for inflammatory memory may be a common trait across different endothelial cell types but also indicate that the specific downstream targets may vary by developmental stage.

Strong Expansion of Human Regulatory T Cells for Adoptive Cell Therapy Results in Epigenetic Changes Which May Impact Their Survival and Function

Adoptive transfer of regulatory T cells (Treg) is a promising new therapeutic option to treat detrimental inflammatory conditions after transplantation and during autoimmune disease. To reach sufficient cell yield for treatment, ex vivo isolated autologous or allogenic Tregs need to be expanded extensively in vitro during manufacturing of the Treg product. However, repetitive cycles of restimulation and prolonged culture have been shown to impact T cell phenotypes, functionality and fitness. It is therefore critical to scrutinize the molecular changes which occur during T cell product generation, and reexamine current manufacturing practices. We performed genome-wide DNA methylation profiling of cells throughout the manufacturing process of a polyclonal Treg product that has proven safety and hints of therapeutic efficacy in kidney transplant patients. We found progressive DNA methylation changes over the duration of culture, which were donor-independent and reproducible between manufacturing runs. Differentially methylated regions (DMRs) in the final products were significantly enriched at promoters and enhancers of genes implicated in T cell activation. Additionally, significant hypomethylation did also occur in promoters of genes implicated in functional exhaustion in conventional T cells, some of which, however, have been reported to strengthen immunosuppressive effector function in Tregs. At the same time, a set of reported Treg-specific demethylated regions increased methylation levels with culture, indicating a possible destabilization of Treg identity during manufacturing, which was independent of the purity of the starting material. Together, our results indicate that the repetitive TCR-mediated stimulation lead to epigenetic changes that might impact functionality of Treg products in multiple ways, by possibly shifting to an effector Treg phenotype with enhanced functional activity or by risking destabilization of Treg identity and impaired TCR activation. Our analyses also illustrate the value of epigenetic profiling for the evaluation of T cell product manufacturing pipelines, which might open new avenues for the improvement of current adoptive Treg therapies with relevance for conventional effector T cell products.

The SEQC2 epigenomics quality control (EpiQC) study

BACKGROUND

Cytosine modifications in DNA such as 5-methylcytosine (5mC) underlie a broad range of developmental processes, maintain cellular lineage specification, and can define or stratify types of cancer and other diseases. However, the wide variety of approaches available to interrogate these modifications has created a need for harmonized materials, methods, and rigorous benchmarking to improve genome-wide methylome sequencing applications in clinical and basic research. Here, we present a multi-platform assessment and cross-validated resource for epigenetics research from the FDA's Epigenomics Quality Control Group.

RESULTS

Each sample is processed in multiple replicates by three whole-genome bisulfite sequencing (WGBS) protocols (TruSeq DNA methylation, Accel-NGS MethylSeq, and SPLAT), oxidative bisulfite sequencing (TrueMethyl), enzymatic deamination method (EMSeq), targeted methylation sequencing (Illumina Methyl Capture EPIC), single-molecule long-read nanopore sequencing from Oxford Nanopore Technologies, and 850k Illumina methylation arrays. After rigorous quality assessment and comparison to Illumina EPIC methylation microarrays and testing on a range of algorithms (Bismark, BitmapperBS, bwa-meth, and BitMapperBS), we find overall high concordance between assays, but also differences in efficiency of read mapping, CpG capture, coverage, and platform performance, and variable performance across 26 microarray normalization algorithms.

CONCLUSIONS

The data provided herein can guide the use of these DNA reference materials in epigenomics research, as well as provide best practices for experimental design in future studies. By leveraging seven human cell lines that are designated as publicly available reference materials, these data can be used as a baseline to advance epigenomics research.

Epigenome-wide association studies: current knowledge, strategies and recommendations

The aetiology and pathophysiology of complex diseases are driven by the interaction between genetic and environmental factors. The variability in risk and outcomes in these diseases are incompletely explained by genetics or environmental risk factors individually. Therefore, researchers are now exploring the epigenome, a biological interface at which genetics and the environment can interact. There is a growing body of evidence supporting the role of epigenetic mechanisms in complex disease pathophysiology. Epigenome-wide association studies (EWASes) investigate the association between a phenotype and epigenetic variants, most commonly DNA methylation. The decreasing cost of measuring epigenome-wide methylation and the increasing accessibility of bioinformatic pipelines have contributed to the rise in EWASes published in recent years. Here, we review the current literature on these EWASes and provide further recommendations and strategies for successfully conducting them. We have constrained our review to studies using methylation data as this is the most studied epigenetic mechanism; microarray-based data as whole-genome bisulphite sequencing remains prohibitively expensive for most laboratories; and blood-based studies due to the non-invasiveness of peripheral blood collection and availability of archived DNA, as well as the accessibility of publicly available blood-cell-based methylation data. Further, we address multiple novel areas of EWAS analysis that have not been covered in previous reviews: (1) longitudinal study designs, (2) the chip analysis methylation pipeline (ChAMP), (3) differentially methylated region (DMR) identification paradigms, (4) methylation quantitative trait loci (methQTL) analysis, (5) methylation age analysis and (6) identifying cell-specific differential methylation from mixed cell data using statistical deconvolution.

Epigenome-Wide Association Study of Prostate Cancer in African Americans Identifies DNA Methylation Biomarkers for Aggressive Disease

DNA methylation plays important roles in prostate cancer (PCa) development and progression. African American men have higher incidence and mortality rates of PCa than other racial groups in U.S. The goal of this study was to identify differentially methylated CpG sites and genes between clinically defined aggressive and nonaggressive PCa in African Americans. We performed genome-wide DNA methylation profiling in leukocyte DNA from 280 African American PCa patients using Illumina MethylationEPIC array that contains about 860K CpG sties. There was a slight increase of overall methylation level (mean β value) with the increasing Gleason scores (GS = 6, GS = 7, GS ≥ 8, P for trend = 0.002). There were 78 differentially methylated CpG sites with P < 10-4 and 9 sites with P < 10-5 in the trend test. We also found 77 differentially methylated regions/genes (DMRs), including 10 homeobox genes and six zinc finger protein genes. A gene ontology (GO) molecular pathway enrichment analysis of these 77 DMRs found that the main enriched pathway was DNA-binding transcriptional factor activity. A few representative DMRs include HOXD8, SOX11, ZNF-471, and ZNF-577. Our study suggests that leukocyte DNA methylation may be valuable biomarkers for aggressive PCa and the identified differentially methylated genes provide biological insights into the modulation of immune response by aggressive PCa.

Common and rare variant association analyses in amyotrophic lateral sclerosis identify 15 risk loci with distinct genetic architectures and neuron-specific biology

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease with a lifetime risk of one in 350 people and an unmet need for disease-modifying therapies. We conducted a cross-ancestry genome-wide association study (GWAS) including 29,612 patients with ALS and 122,656 controls, which identified 15 risk loci. When combined with 8,953 individuals with whole-genome sequencing (6,538 patients, 2,415 controls) and a large cortex-derived expression quantitative trait locus (eQTL) dataset (MetaBrain), analyses revealed locus-specific genetic architectures in which we prioritized genes either through rare variants, short tandem repeats or regulatory effects. ALS-associated risk loci were shared with multiple traits within the neurodegenerative spectrum but with distinct enrichment patterns across brain regions and cell types. Of the environmental and lifestyle risk factors obtained from the literature, Mendelian randomization analyses indicated a causal role for high cholesterol levels. The combination of all ALS-associated signals reveals a role for perturbations in vesicle-mediated transport and autophagy and provides evidence for cell-autonomous disease initiation in glutamatergic neurons.

Exploring the impact of night shift work on methylation of circadian genes

Night shift work is associated with increased breast cancer risk, but the molecular mechanisms are not well-understood. The objective of this study was to explore the relationship between night shift work parameters (current status, duration/years, and intensity) and methylation in circadian genes as a potential mechanism underlying the carcinogenic effects of night shift work. A cross-sectional study was conducted among 74 female healthcare employees (n = 38 day workers, n = 36 night shift workers). The Illumina Infinium MethylationEPIC beadchip was applied to DNA extracted from blood samples to measure methylation using a candidate gene approach at 1150 CpG loci across 22 circadian genes. Linear regression models were used to examine the association between night shift work parameters and continuous methylation measurements (β-values) for each CpG site. The false-discovery rate (q = 0.2) was used to account for multiple comparisons. Compared to day workers, current night shift workers demonstrated hypermethylation in the 5'UTR region of CSNK1E (q = 0.15). Individuals that worked night shifts for ≥10 years exhibited hypomethylation in the gene body of NR1D1 (q = 0.08) compared to those that worked <10 years. Hypermethylation in the gene body of ARNTL was also apparent in those who worked ≥3 consecutive night shifts a week (q = 0.18). These findings suggest that night shift work is associated with differential methylation in core circadian genes, including CSNK1E, NR1D1 and ARNTL. Future, larger-scale studies with long-term follow-up and detailed night shift work assessment are needed to confirm and expand on these findings.

Genome-wide differentially methylated genes associated with posttraumatic stress disorder and longitudinal change in methylation in rape survivors

Rape is associated with a high risk for posttraumatic stress disorder (PTSD). DNA methylation changes may confer risk or protection for PTSD following rape by regulating the expression of genes implicated in pathways affected by PTSD. We aimed to: (1) identify epigenome-wide differences in methylation profiles between rape-exposed women with and without PTSD at 3-months post-rape, in a demographically and ethnically similar group, drawn from a low-income setting; (2) validate and replicate the findings of the epigenome-wide analysis in selected genes (BRSK2 and ADCYAP1); and (3) investigate baseline and longitudinal changes in BRSK2 and ADCYAP1 methylation over six months in relation to change in PTSD symptom scores over 6 months, in the combined discovery/validation and replication samples (n = 96). Rape-exposed women (n = 852) were recruited from rape clinics in the Rape Impact Cohort Evaluation (RICE) umbrella study. Epigenome-wide differentially methylated CpG sites between rape-exposed women with (n = 24) and without (n = 24) PTSD at 3-months post-rape were investigated using the Illumina EPIC BeadChip in a discovery cohort (n = 48). Validation (n = 47) and replication (n = 49) of BRSK2 and ADCYAP1 methylation findings were investigated using EpiTYPER technology. Longitudinal change in BRSK2 and ADCYAP1 was also investigated using EpiTYPER technology in the combined sample (n = 96). In the discovery sample, after adjustment for multiple comparisons, one differentially methylated CpG site (chr10: 61385771/ cg01700569, p = 0.049) and thirty-four differentially methylated regions were associated with PTSD status at 3-months post-rape. Decreased BRSK2 and ADCYAP1 methylation at 3-months and 6-months post-rape were associated with increased PTSD scores at the same time points, but these findings did not remain significant in adjusted models. In conclusion, decreased methylation of BRSK2 may result in abnormal neuronal polarization, synaptic development, vesicle formation, and disrupted neurotransmission in individuals with PTSD. PTSD symptoms may also be mediated by differential methylation of the ADCYAP1 gene which is involved in stress regulation. Replication of these findings is required to determine whether ADCYAP1 and BRSK2 are biomarkers of PTSD and potential therapeutic targets.

Mitochondrial genome-wide analysis of nuclear DNA methylation quantitative trait loci

Mitochondria have a complex communication network with the surrounding cell and can alter nuclear DNA methylation (DNAm). Variation in the mitochondrial DNA (mtDNA) has also been linked to differential DNAm. Genome-wide association studies have identified numerous DNAm quantitative trait loci, but these studies have not examined the mitochondrial genome. Herein, we quantified nuclear DNAm from blood and conducted a mitochondrial genome-wide association study of DNAm, with an additional emphasis on sex- and prediabetes-specific heterogeneity. We used the Young Finns Study (n = 926) with sequenced mtDNA genotypes as a discovery sample and sought replication in the Ludwigshafen Risk and Cardiovascular Health study (n = 2317). We identified numerous significant associations in the discovery phase (P &lt; 10−9), but they were not replicated when accounting for multiple testing. In total, 27 associations were nominally replicated with a P &lt; 0.05. The replication analysis presented no evidence of sex- or prediabetes-specific heterogeneity. The 27 associations were included in a joint meta-analysis of the two cohorts, and 19 DNAm sites associated with mtDNA variants, while four other sites showed haplogroup associations. An expression quantitative trait methylation analysis was performed for the identified DNAm sites, pinpointing two statistically significant associations. This study provides evidence of a mitochondrial genetic control of nuclear DNAm with little evidence found for sex- and prediabetes-specific effects. The lack of a comparable mtDNA data set for replication is a limitation in our study and further studies are needed to validate our results.

Machine Learning-Based Single Cell and Integrative Analysis Reveals That Baseline mDC Predisposition Correlates With Hepatitis B Vaccine Antibody Response

Vaccination to prevent infectious disease is one of the most successful public health interventions ever developed. And yet, variability in individual vaccine effectiveness suggests that a better mechanistic understanding of vaccine-induced immune responses could improve vaccine design and efficacy. We have previously shown that protective antibody levels could be elicited in a subset of recipients with only a single dose of the hepatitis B virus (HBV) vaccine and that a wide range of antibody levels were elicited after three doses. The immune mechanisms responsible for this vaccine response variability is unclear. Using single cell RNA sequencing of sorted innate immune cell subsets, we identified two distinct myeloid dendritic cell subsets (NDRG1-expressing mDC2 and CDKN1C-expressing mDC4), the ratio of which at baseline (pre-vaccination) correlated with the immune response to a single dose of HBV vaccine. Our results suggest that the participants in our vaccine study were in one of two different dendritic cell dispositional states at baseline – an NDRG2-mDC2 state in which the vaccine elicited an antibody response after a single immunization or a CDKN1C-mDC4 state in which the vaccine required two or three doses for induction of antibody responses. To explore this correlation further, genes expressed in these mDC subsets were used for feature selection prior to the construction of predictive models using supervised canonical correlation machine learning. The resulting models showed an improved correlation with serum antibody titers in response to full vaccination. Taken together, these results suggest that the propensity of circulating dendritic cells toward either activation or suppression, their “dispositional endotype” at pre-vaccination baseline, could dictate response to vaccination.

DNA methylation signatures reveal that distinct combinations of transcription factors specify human immune cell epigenetic identity

Epigenetic reprogramming underlies specification of immune cell lineages, but patterns that uniquely define immune cell types and the mechanisms by which they are established remain unclear. Here, we identified lineage-specific DNA methylation signatures of six immune cell types from human peripheral blood and determined their relationship to other epigenetic and transcriptomic patterns. Sites of lineage-specific hypomethylation were associated with distinct combinations of transcription factors in each cell type. By contrast, sites of lineage-specific hypermethylation were restricted mostly to adaptive immune cells. PU.1 binding sites were associated with lineage-specific hypo- and hypermethylation in different cell types, suggesting that it regulates DNA methylation in a context-dependent manner. These observations indicate that innate and adaptive immune lineages are specified by distinct epigenetic mechanisms via combinatorial and context-dependent use of key transcription factors. The cell-specific epigenomics and transcriptional patterns identified serve as a foundation for future studies on immune dysregulation in diseases and aging.

DNA Methylation Patterns in CD8+ T Cells Discern Psoriasis From Psoriatic Arthritis and Correlate With Cutaneous Disease Activity

Background: Psoriasis is a T cell-mediated chronic autoimmune/inflammatory disease. While some patients experience disease limited to the skin (skin psoriasis), others develop joint involvement (psoriatic arthritis; PsA). In the absence of disease- and/or outcome-specific biomarkers, and as arthritis can precede skin manifestations, diagnostic and therapeutic delays are common and contribute to disease burden and damage accrual.

Objective: Altered epigenetic marks, including DNA methylation, contribute to effector T cell phenotypes and altered cytokine expression in autoimmune/inflammatory diseases. This project aimed at the identification of disease-/outcome-specific DNA methylation signatures in CD8⁺ T cells from patients with psoriasis and PsA as compared to healthy controls.

Method: Peripheral blood CD8⁺ T cells from nine healthy controls, 10 psoriasis, and seven PsA patients were collected to analyze DNA methylation marks using Illumina Human Methylation EPIC BeadChips (>850,000 CpGs per sample). Bioinformatic analysis was performed using R (minfi, limma, ChAMP, and DMRcate packages).

Results: DNA methylation profiles in CD8⁺ T cells differentiate healthy controls from psoriasis patients [397 Differentially Methylated Positions (DMPs); 9 Differentially Methylated Regions (DMRs) when ≥CpGs per DMR were considered; 2 DMRs for ≥10 CpGs]. Furthermore, patients with skin psoriasis can be discriminated from PsA patients [1,861 DMPs, 20 DMRs (≥5 CpGs per region), 4 DMRs (≥10 CpGs per region)]. Gene ontology (GO) analyses considering genes with ≥1 DMP in their promoter delivered methylation defects in skin psoriasis and PsA primarily affecting the BMP signaling pathway and endopeptidase regulator activity, respectively. GO analysis of genes associated with DMRs between skin psoriasis and PsA demonstrated an enrichment of GABAergic neuron and cortex neuron development pathways. Treatment with cytokine blockers associated with DNA methylation changes [2,372 DMPs; 1,907 DMPs within promoters, 7 DMRs (≥5 CpG per regions)] affecting transforming growth factor beta receptor and transmembrane receptor protein serine/threonine kinase signaling pathways. Lastly, a methylation score including TNF and IL-17 pathway associated DMPs inverse correlates with skin disease activity scores (PASI).

Conclusion: Patients with skin psoriasis exhibit DNA methylation patterns in CD8⁺ T cells that allow differentiation from PsA patients and healthy individuals, and reflect clinical activity of skin disease. Thus, DNA methylation profiling promises potential as diagnostic and prognostic tool to be used for molecular patient stratification toward individualized treatment.

DNA methylation profiling identifies epigenetic differences between early versus late stages of diabetic chronic kidney disease

BACKGROUND

We investigated a cross-sectional epigenome-wide association study of patients with early and late diabetes-associated chronic kidney disease (CKD) to identify possible epigenetic differences between the two groups as well as changes in methylation across all stages of diabetic CKD. We also evaluated the potential of using a panel of identified 5'-C-phosphate-G-3' (CpG) sites from this cohort to predict the progression of diabetic CKD.

METHODS

This cross-sectional study recruited 119 adults. DNA was extracted from blood using the Qiagen QIAampDNA Mini Spin Kit. Genome-wide methylation analysis was performed using Illumina Infinium MethylationEPIC BeadChips (HM850K). Intensity data files were processed and analysed using the minfi and MissMethyl packages for R. We examined the degree of methylation of CpG sites in early versus late diabetic CKD patients for CpG sites with an unadjusted P-value <0.01 and an absolute change in methylation of 5% (n = 239 CpG sites).

RESULTS

Hierarchical clustering of the 239 CpG sites largely separated the two groups. A heat map for all 239 CpG sites demonstrated distinct methylation patterns in the early versus late groups, with CpG sites showing evidence of progressive change. Based on our differentially methylated region (DMR) analysis of the 239 CpG sites, we highlighted two DMRs, namely the cysteine-rich secretory protein 2 (CRISP2) and piwi-like RNA-mediated gene silencing 1 (PIWIL1) genes. The best predictability for the two groups involved a receiver operating characteristics curve of eight CpG sites alone and achieved an area under the curve of 0.976.

CONCLUSIONS

We have identified distinct DNA methylation patterns between early and late diabetic CKD patients as well as demonstrated novel findings of potential progressive methylation changes across all stages (1-5) of diabetic CKD at specific CpG sites. We have also identified associated genes CRISP2 and PIWIL1, which may have the potential to act as stage-specific diabetes-associated CKD markers, and showed that the use of a panel of eight identified CpG sites alone helps to increase the predictability for the two groups.

Effect of exercise and diet intervention in NAFLD and NASH via GAB2 methylation

BACKGROUND

Nonalcoholic fatty liver disease (NAFLD) is a disorder that extends from simple hepatic steatosis to nonalcoholic steatohepatitis (NASH), which is effectively alleviated by lifestyle intervention. Nevertheless, DNA methylation mechanism underling the effect of environmental factors on NAFLD and NASH is still obscure. The aim of this study was to investigate the effect of exercise and diet intervention in NAFLD and NASH via DNA methylation of GAB2.

METHODS

Methylation of genomic DNA in human NAFLD was quantified using Infinium Methylation EPIC BeadChip assay after exercise (Ex), low carbohydrate diet (LCD) and exercise plus low carbohydrate diet (ELCD) intervention. The output Idat files were processed using ChAMP package. False discovery rate on genome-wide analysis of DNA methylation (q < 0.05), and cytosine-guanine dinucleotides (CpGs) which are located in promoters were used for subsequent analysis (|Δβ|≥ 0.1). K-means clustering was used to cluster differentially methylated genes according to 3D genome information from Human embryonic stem cell. To quantify DNA methylation and mRNA expression of GRB2 associated binding protein 2 (GAB2) in NASH mice after Ex, low fat diet (LFD) and exercise plus low fat diet (ELFD), MassARRAY EpiTYPER and quantitative reverse transcription polymerase chain reaction were used.

RESULTS

Both LCD and ELCD intervention on human NAFLD can induce same DNA methylation alterations at critical genes in blood, e.g., GAB2, which was also validated in liver and adipose of NASH mice after LFD and ELFD intervention. Moreover, methylation of CpG units (i.e., CpG_10.11.12) inversely correlated with mRNA expression GAB2 in adipose tissue of NASH mice after ELFD intervention.

CONCLUSIONS

We highlighted the susceptibility of DNA methylation in GAB2 to ELFD intervention, through which exercise and diet can protect against the progression of NAFLD and NASH on the genome level, and demonstrated that the DNA methylation variation in blood could mirror epigenetic signatures in target tissues of important biological function, i.e., liver and adipose tissue. Trial registration International Standard Randomized Controlled Trial Number Register (ISRCTN 42622771).