Genetic impacts on DNA methylation: research findings and future perspectives

Transgenerational hormesis in healthy aging and antiaging medicine from bench to clinics: Role of food components

Neurodegenerative diseases have multifactorial pathogenesis, mainly involving neuroinflammatory processes. Finding drugs able to treat these diseases, expecially because for most of these diseases there are no effective drugs, and the current drugs cause undesired side effects, represent a crucial point. Most in vivo and in vitro studies have been concentrated on various aspects related to neurons (e.g. neuroprotection), however, there has not been focus on the prevention of early stages involving glial cell activation and neuroinflammation. Recently, it has been demonstrated that nutritional phytochemicals including polyphenols, the main active constituents of the Mediterranean diet, maintain redox balance and neuroprotection through the activation of hormetic vitagene pathway. Recent lipidomics data from our laboratory indicate mushrooms as strong nutritional neuronutrients with strongly activity against neuroinflammation in Meniere' diseaseas, a model of cochleovestibular neural degeneration, as well as in animal model of traumatic brain injury, or rotenone induced parkinson's disease. Moreover, Hidrox®, an aqueous extract of olive containing hydroxytyrosol, and Boswellia, acting as Nrf2 activators, promote resilience by enhancing the redox potential, and thus, regulate through hormetic mechanisms, cellular stress response mechanisms., Thus, modulation of cellular stress pathways, in particular vitagenes system, may be an innovative approach for therapeutic intervention in neurodegenerative disorders.

Genome-Wide Methylation Analysis Reveals a KCNK3-Prominent Causal Cascade on Hypertension

BACKGROUND

Despite advances in understanding hypertension's genetic structure, how noncoding genetic variants influence it remains unclear. Studying their interaction with DNA methylation is crucial to deciphering this complex disease's genetic mechanisms.

METHODS

We investigated the genetic and epigenetic interplay in hypertension using whole-genome bisulfite sequencing. Methylation profiling in 918 males revealed allele-specific methylation and methylation quantitative trait loci. We engineered rs1275988T/C mutant mice using CRISPR (clustered regularly interspaced short palindromic repeats)/Cas9 (CRISPR-associated protein 9), bred them for homozygosity, and subjected them to a high-salt diet. Telemetry captured their cardiovascular metrics. Protein-DNA interactions were elucidated using DNA pull-downs, mass spectrometry, and Western blots. A wire myograph assessed vascular function, and analysis of the Kcnk3 gene methylation highlighted the mutation's role in hypertension.

RESULTS

We discovered that DNA methylation-associated genetic effects, especially in non-cytosine-phosphate-guanine (non-CpG) island and noncoding distal regulatory regions, significantly contribute to hypertension predisposition. We identified distinct methylation quantitative trait locus patterns in the hypertensive population and observed that the onset of hypertension is influenced by the transmission of genetic effects through the demethylation process. By evidence-driven prioritization and in vivo experiments, we unearthed rs1275988 in a cell type-specific enhancer as a notable hypertension causal variant, intensifying hypertension through the modulation of local DNA methylation and consequential alterations in Kcnk3 gene expression and vascular remodeling. When exposed to a high-salt diet, mice with the rs1275988C/C genotype exhibited exacerbated hypertension and significant vascular remodeling, underscored by increased aortic wall thickness. The C allele of rs1275988 was associated with elevated DNA methylation levels, driving down the expression of the Kcnk3 gene by attenuating Nr2f2 (nuclear receptor subfamily 2 group F member 2) binding at the enhancer locus.

CONCLUSIONS

Our research reveals new insights into the complex interplay between genetic variations and DNA methylation in hypertension. We underscore hypomethylation's potential in hypertension onset and identify rs1275988 as a causal variant in vascular remodeling. This work advances our understanding of hypertension's molecular mechanisms and encourages personalized health care strategies.

Enhancing Disease Risk Gene Discovery by Integrating Transcription Factor-Linked Trans-located Variants into Transcriptome-Wide Association Analyses

Transcriptome-wide association studies (TWAS) have been successful in identifying disease susceptibility genes by integrating cis-variants predicted gene expression with genome-wide association studies (GWAS) data. However, trans-located variants for predicting gene expression remain largely unexplored. Here, we introduce transTF-TWAS, which incorporates transcription factor (TF)-linked trans-located variants to enhance model building. Using data from the Genotype-Tissue Expression project, we predict gene expression and alternative splicing and applied these models to large GWAS datasets for breast, prostate, and lung cancers. We demonstrate that transTF-TWAS outperforms other existing TWAS approaches in both constructing gene prediction models and identifying disease-associated genes, as evidenced by simulations and real data analysis. Our transTF-TWAS approach significantly contributes to the discovery of disease risk genes. Findings from this study have shed new light on several genetically driven key regulators and their associated regulatory networks underlying disease susceptibility.

Epigenetic diversity of genes with copy number variations among natural populations of the three-spined stickleback

Duplicated genes provide the opportunity for evolutionary novelty and adaptive divergence. In many cases, having more gene copies increases gene expression, which might facilitate adaptation to stressful or novel environments. Conversely, overexpression or misexpression of duplicated genes can be detrimental and subject to negative selection. In this scenario, newly duplicate genes may evade purifying selection if they are epigenetically silenced, at least temporarily, leading them to persist in populations as copy number variations (CNVs). In animals and plants, younger gene duplicates tend to have higher levels of DNA methylation and lower levels of gene expression, suggesting epigenetic regulation could promote the retention of gene duplications via expression repression or silencing. Here, we test the hypothesis that DNA methylation variation coincides with young duplicate genes that are segregating as CNVs in six populations of the three-spined stickleback that span a salinity gradient from 4 to 30 PSU. Using reduced-representation bisulfite sequencing, we found DNA methylation and CNV differentiation outliers rarely overlapped. Whereas lineage-specific genes and young duplicates were found to be highly methylated, just two gene CNVs showed a significant association between promoter methylation level and copy number, suggesting that DNA methylation might not interact with CNVs in our dataset. If most new duplications are regulated for dosage by epigenetic mechanisms, our results do not support a strong contribution from DNA methylation soon after duplication. Instead, our results are consistent with a preference to duplicate genes that are already highly methylated.

Integrative Genetic Variation, DNA Methylation, and Gene Expression Analysis of Escitalopram and Aripiprazole Treatment Outcomes in Depression: A CAN-BIND-1 Study

INTRODUCTION

Little is known about the interplay between genetics and epigenetics on antidepressant treatment (1) response and remission, (2) side effects, and (3) serum levels. This study explored the relationship among single nucleotide polymorphisms (SNPs), DNA methylation (DNAm), and mRNA levels of four pharmacokinetic genes, CYP2C19, CYP2D6, CYP3A4, and ABCB1, and its effect on these outcomes.

METHODS

The Canadian Biomarker Integration Network for Depression-1 dataset consisted of 177 individuals with major depressive disorder treated for 8 weeks with escitalopram (ESC) followed by 8 weeks with ESC monotherapy or augmentation with aripiprazole. DNAm quantitative trait loci (mQTL), identified by SNP-CpG associations between 20 SNPs and 60 CpG sites in whole blood, were tested for associations with our outcomes, followed by causal inference tests (CITs) to identify methylation-mediated genetic effects.

RESULTS

Eleven cis-SNP-CpG pairs (q<0.05) constituting four unique SNPs were identified. Although no significant associations were observed between mQTLs and response/remission, CYP2C19 rs4244285 was associated with treatment-related weight gain (q=0.027) and serum concentrations of ESCadj (q<0.001). Between weeks 2-4, 6.7% and 14.9% of those with *1/*1 (normal metabolizers) and *1/*2 (intermediate metabolizers) genotypes, respectively, reported ≥2 lbs of weight gain. In contrast, the *2/*2 genotype (poor metabolizers) did not report weight gain during this period and demonstrated the highest ESCadj concentrations. CITs did not indicate that these effects were epigenetically mediated.

DISCUSSION

These results elucidate functional mechanisms underlying the established associations between CYP2C19 rs4244285 and ESC pharmacokinetics. This mQTL SNP as a marker for antidepressant-related weight gain needs to be further explored.

High-Resolution Melting (HRM) analysis of DNA methylation using semiconductor chip-based digital PCR

BACKGROUND

Digital PCR (dPCR) technology allows absolute quantification and detection of disease-associated rare variants, and thus the use of dPCR technology has been increasing in clinical research and diagnostics. The high-resolution melting curve analysis (HRM) of qPCR is widely used to distinguish true positives from false positives and detect rare variants. In particular, qPCR-HRM is commonly used for methylation assessment in research and diagnostics due to its simplicity and high reproducibility. Most dPCR instruments have limited fluorescence channels available and separate heating and imaging systems. Therefore, it is difficult to perform HRM analysis using dPCR instruments.

OBJECTIVE

A new digital real-time PCR instrument (LOAA) has been recently developed to integrate partitioning, thermocycling, and imaging in a single dPCR instrument. In addition, a new technique to perform HRM analysis is utilized in LOAA. The aim of the present study is to evaluate the efficiency and accuracy of LOAA dPCR on HRM analysis for the detection of methylation.

METHODS

In this study, comprehensive comparison with Bio-Rad qRT-PCR and droplet-based dPCR equipment was performed to verify the HRM analysis-based methylation detection efficiency of the LOAA digital PCR equipment. Here, sodium bisulfite modification method was applied to detect methylated DNA sequences by each PCR method.

RESULTS

Melting curve analysis detected four different Tm values using LOAA and qPCR, and found that LOAA, unlike qPCR, successfully distinguished between different Tm values when the Tm values were very similar. In addition, melting temperatures increased by each methylation were about 0.5℃ for qPCR and about 0.2 ~ 0.6℃ for LOAA. The melting temperature analyses of methylated and unmethylated DNA samples were conducted using LOAA dPCR with TaqMan probes and EvaGreen, and the result found that Tm values of methylated DNA samples are higher than those of unmethylated DNA samples.

CONCLUSION

The present study shows that LOAA dPCR could detect different melting temperatures according to methylation status of target sequences, indicating that LOAA dPCR would be useful for diagnostic applications that require the accurate quantification and assessment of DNA methylation.

Epigenetic alterations in patients with anorexia nervosa-a systematic review

Anorexia nervosa (AN) is a complex metabolic and psychological disorder that is influenced by both heritable genetic components and environmental factors. Exposure to various environmental influences can lead to epigenetically induced changes in gene expression. Epigenetic research in AN is still in its infancy, and studies to date are limited in determining clear, valid links to disease onset and progression are limited. Therefore, the aim of this systematic review was to compile and critically evaluate the available results of epigenetic studies specifically in AN and to provide recommendations for future studies. In accordance with the PRISMA guidelines, a systematic literature search was performed in three different databases (PubMed, Embase, and Web of Science) through May 2023. Twenty-three original papers or conference abstracts on epigenetic studies in AN were collected. Epigenome-wide association studies (EWASs), which analyze DNA methylation across the genome in patients with AN and identify potential disease-relevant changes in promoter/regulatory regions of genes, are the most promising for future research. To date, five EWASs on AN have been published, suggesting a potential reversibility of malnutrition-induced epigenetic changes once patients recover. Hence, determining differential DNA methylation levels could serve as a biomarker for disease status or early diagnosis and might be involved in disease progression or chronification. For future research, EWASs with a larger sample size, longitudinal study design and uniform methods should be performed to contribute to the understanding of the pathophysiology of AN, the development of individual interventions and a better prognosis for affected patients.

Osteoarthritis as an Enhanceropathy: Gene Regulation in Complex Musculoskeletal Disease

PURPOSE OF REVIEW

Osteoarthritis is a complex and highly polygenic disease. Over 100 reported osteoarthritis risk variants fall in non-coding regions of the genome, ostensibly conferring functional effects through the disruption of regulatory elements impacting target gene expression. In this review, we summarise the progress that has advanced our knowledge of gene enhancers both within the field of osteoarthritis and more broadly in complex diseases.

RECENT FINDINGS

Advances in technologies such as ATAC-seq have facilitated our understanding of chromatin states in specific cell types, bolstering the interpretation of GWAS and the identification of effector genes. Their application to osteoarthritis research has revealed enhancers as the principal regulatory element driving disease-associated changes in gene expression. However, tissue-specific effects in gene regulatory mechanisms can contribute added complexity to biological interpretation. Understanding gene enhancers and their altered activity in specific cell and tissue types is the key to unlocking the genetic complexity of osteoarthritis. The use of single-cell technologies in osteoarthritis research is still in its infancy. However, such tools offer great promise in improving our functional interpretation of osteoarthritis GWAS and the identification of druggable targets. Large-scale collaborative efforts will be imperative to understand tissue and cell-type specific molecular mechanisms underlying enhancer function in disease.

Regulatory and evolutionary impact of DNA methylation in two songbird species and their naturally occurring F1 hybrids

BACKGROUND

Regulation of transcription by DNA methylation in 5'-CpG-3' context is a widespread mechanism allowing differential expression of genetically identical cells to persist throughout development. Consequently, differences in DNA methylation can reinforce variation in gene expression among cells, tissues, populations, and species. Despite a surge in studies on DNA methylation, we know little about the importance of DNA methylation in population differentiation and speciation. Here we investigate the regulatory and evolutionary impact of DNA methylation in five tissues of two Ficedula flycatcher species and their naturally occurring F1 hybrids.

RESULTS

We show that the density of CpG in the promoters of genes determines the strength of the association between DNA methylation and gene expression. The impact of DNA methylation on gene expression varies among tissues with the brain showing unique patterns. Differentially expressed genes between parental species are predicted by genetic and methylation differentiation in CpG-rich promoters. However, both these factors fail to predict hybrid misexpression suggesting that promoter mismethylation is not a main determinant of hybrid misexpression in Ficedula flycatchers. Using allele-specific methylation estimates in hybrids, we also determine the genome-wide contribution of cis- and trans effects in DNA methylation differentiation. These distinct mechanisms are roughly balanced in all tissues except the brain, where trans differences predominate.

CONCLUSIONS

Overall, this study provides insight on the regulatory and evolutionary impact of DNA methylation in songbirds.

No association between peripheral serotonin-gene-related DNA methylation and brain serotonin neurotransmission in the healthy and depressed state

BACKGROUND

Methylation of serotonin-related genes has been proposed as a plausible gene-by-environment link which may mediate environmental stress, depressive and anxiety symptoms. DNA methylation is often measured in blood cells, but little is known about the association between this peripheral epigenetic modification and brain serotonergic architecture. Here, we evaluated the association between whole-blood-derived methylation of four CpG sites in the serotonin transporter (SLC6A4) and six CpG sites of the tryptophan hydroxylase 2 (TPH2) gene and in-vivo brain levels of serotonin transporter (5-HTT) and serotonin 4 receptor (5-HT4) in a cohort of healthy individuals (N = 254) and, for 5-HT4, in a cohort of unmedicated patients with depression (N = 90). To do so, we quantified SLC6A4/TPH2 methylation using bisulfite pyrosequencing and estimated brain 5-HT4 and 5-HTT levels using positron emission tomography. In addition, we explored the association between SLC6A4 and TPH2 methylation and measures of early life and recent stress, depressive and anxiety symptoms on 297 healthy individuals.

RESULTS

We found no statistically significant association between peripheral DNA methylation and brain markers of serotonergic neurotransmission in patients with depression or in healthy individuals. In addition, although SLC6A4 CpG2 (chr17:30,236,083) methylation was marginally associated with the parental bonding inventory overprotection score in the healthy cohort, statistical significance did not remain after accounting for blood cell heterogeneity.

CONCLUSIONS

We suggest that findings on peripheral DNA methylation in the context of brain serotonin-related features should be interpreted with caution. More studies are needed to rule out a role of SLC6A4 and TPH2 methylation as biomarkers for environmental stress, depressive or anxiety symptoms.

A review of the role of epigenetic studies for intramuscular fat deposition in beef cattle

Intramuscular fat (IMF) deposition profoundly influences meat quality and economic value in beef cattle production. Meanwhile, contemporary developments in epigenetics have opened new outlooks for understanding the molecular basics of IMF regulation, and it has become a key area of research for world scholars. Therefore, the aim of this paper was to provide insight and synthesis into the intricate relationship between epigenetic mechanisms and IMF deposition in beef cattle. The methodology involves a thorough analysis of existing literature, including pertinent books, academic journals, and online resources, to provide a comprehensive overview of the role of epigenetic studies in IMF deposition in beef cattle. This review summarizes the contemporary studies in epigenetic mechanisms in IMF regulation, high-resolution epigenomic mapping, single-cell epigenomics, multi-omics integration, epigenome editing approaches, longitudinal studies in cattle growth, environmental epigenetics, machine learning in epigenetics, ethical and regulatory considerations, and translation to industry practices from perspectives of IMF deposition in beef cattle. Moreover, this paper highlights DNA methylation, histone modifications, acetylation, phosphorylation, ubiquitylation, non-coding RNAs, DNA hydroxymethylation, epigenetic readers, writers, and erasers, chromatin immunoprecipitation followed by sequencing, whole genome bisulfite sequencing, epigenome-wide association studies, and their profound impact on the expression of crucial genes governing adipogenesis and lipid metabolism. Nutrition and stress also have significant influences on epigenetic modifications and IMF deposition. The key findings underscore the pivotal role of epigenetic studies in understanding and enhancing IMF deposition in beef cattle, with implications for precision livestock farming and ethical livestock management. In conclusion, this review highlights the crucial significance of epigenetic pathways and environmental factors in affecting IMF deposition in beef cattle, providing insightful information for improving the economics and meat quality of cattle production.

Integrated analysis of blood DNA methylation, genetic variants, circulating proteins, microRNAs, and kidney failure in type 1 diabetes

Variation in DNA methylation (DNAmet) in white blood cells and other cells/tissues has been implicated in the etiology of progressive diabetic kidney disease (DKD). However, the specific mechanisms linking DNAmet variation in blood cells with risk of kidney failure (KF) and utility of measuring blood cell DNAmet in personalized medicine are not clear. We measured blood cell DNAmet in 277 individuals with type 1 diabetes and DKD using Illumina EPIC arrays; 51% of the cohort developed KF during 7 to 20 years of follow-up. Our epigenome-wide analysis identified DNAmet at 17 CpGs (5'-cytosine-phosphate-guanine-3' loci) associated with risk of KF independent of major clinical risk factors. DNAmet at these KF-associated CpGs remained stable over a median period of 4.7 years. Furthermore, DNAmet variations at seven KF-associated CpGs were strongly associated with multiple genetic variants at seven genomic regions, suggesting a strong genetic influence on DNAmet. The effects of DNAmet variations at the KF-associated CpGs on risk of KF were partially mediated by multiple KF-associated circulating proteins and KF-associated circulating miRNAs. A prediction model for risk of KF was developed by adding blood cell DNAmet at eight selected KF-associated CpGs to the clinical model. This updated model significantly improved prediction performance (c-statistic = 0.93) versus the clinical model (c-statistic = 0.85) at P = 6.62 × 10-14. In conclusion, our multiomics study provides insights into mechanisms through which variation of DNAmet may affect KF development and shows that blood cell DNAmet at certain CpGs can improve risk prediction for KF in T1D.

Genetic variation and pesticide exposure influence blood DNA methylation signatures in females with early-stage Parkinson's disease

Although sex, genetics, and exposures can individually influence risk for sporadic Parkinson's disease (PD), the joint contributions of these factors to the epigenetic etiology of PD have not been comprehensively assessed. Here, we profiled sex-stratified genome-wide blood DNAm patterns, SNP genotype, and pesticide exposure in agricultural workers (71 early-stage PD cases, 147 controls) and explored replication in three independent samples of varying demographics (n = 218, 222, and 872). Using a region-based approach, we found more associations of blood DNAm with PD in females (69 regions) than in males (2 regions, Δβadj| ≥0.03, padj ≤ 0.05). For 48 regions in females, models including genotype or genotype and pesticide exposure substantially improved in explaining interindividual variation in DNAm (padj ≤ 0.05), and accounting for these variables decreased the estimated effect of PD on DNAm. The results suggested that genotype, and to a lesser degree, genotype-exposure interactions contributed to variation in PD-associated DNAm. Our findings should be further explored in larger study populations and in experimental systems, preferably with precise measures of exposure.

Genome- and epigenome-wide association studies identify susceptibility of CpG sites and regions for metabolic syndrome in a Korean population

BACKGROUND

While multiple studies have investigated the relationship between metabolic syndrome (MetS) and its related traits (fasting glucose, triglyceride, HDL cholesterol, blood pressure, waist circumference) and DNA methylation, our understanding of the epigenetic mechanisms in MetS remains limited. Therefore, we performed an epigenome-wide meta-analysis of blood DNA methylation to identify differentially methylated probes (DMPs) and differentially methylated regions (DMRs) associated with MetS and its components using two independent cohorts comprising a total of 2,334 participants. We also investigated the specific genetic effects on DNA methylation, identified methylation quantitative trait loci (meQTLs) through genome-wide association studies and further utilized Mendelian randomization (MR) to assess how these meQTLs subsequently influence MetS status.

RESULTS

We identified 40 DMPs and 27 DMRs that are significantly associated with MetS. In addition, we identified many novel DMPs and DMRs underlying inflammatory and steroid hormonal processes. The most significant associations were observed in 3 DMPs (cg19693031, cg26974062, cg02988288) and a DMR (chr1:145440444-145441553) at the TXNIP, which are involved in lipid metabolism. These CpG sites were identified as coregulators of DNA methylation in MetS, TG and FAG levels. We identified a total of 144 cis-meQTLs, out of which only 13 were found to be associated with DMPs for MetS. Among these, we confirmed the identified causal mediators of genetic effects at CpG sites cg01881899 at ABCG1 and cg00021659 at the TANK genes for MetS.

CONCLUSIONS

This study observed whether specific CpGs and methylated regions act independently or are influenced by genetic effects for MetS and its components in the Korean population. These associations between the identified DNA methylation and MetS, along with its individual components, may serve as promising targets for the development of preventive interventions for MetS.

Investigating the relationship between DNA methylation, genetic variation, and suicide attempt in bipolar disorder

Individuals with bipolar disorder are at increased risk for suicide, and this can be influenced by a range of biological, clinical, and environmental risk factors. Biological components associated with suicide include DNA modifications that lead to changes in gene expression. Common genetic variation and DNA methylation changes are some of the most frequent types of DNA findings associated with an increased risk for suicidal behavior. Importantly, the interplay between genetic predisposition and DNA methylation patterns is becoming more prevalent in genetic studies. We hypothesized that DNA methylation patterns in specific loci already genetically associated with suicide would be altered in individuals with bipolar disorder and a history of suicide attempt. To test this hypothesis, we searched the literature to identify common genetic variants (N=34) previously associated with suicidal thoughts and behaviors in individuals with bipolar disorder. We then created a customized sequencing panel that covered our chosen genomic loci. We profiled DNA methylation patterns from blood samples collected from bipolar disorder participants with suicidal behavior (N=55) and without suicidal behavior (N=51). We identified seven differentially methylated CpG sites and five differentially methylated regions between the two groups. Additionally, we found that DNA methylation changes in MIF and CACNA1C were associated with lethality or number of suicide attempts. Finally, we identified three meQTLs in SIRT1 , IMPA2 , and INPP1 . This study illustrates that DNA methylation is altered in individuals with bipolar disorder and a history of suicide attempts in regions known to harbor suicide-related variants.

QTL mapping of human retina DNA methylation identifies 87 gene-epigenome interactions in age-related macular degeneration

DNA methylation provides a crucial epigenetic mark linking genetic variations to environmental influence. We have analyzed array-based DNA methylation profiles of 160 human retinas with co-measured RNA-seq and >8 million genetic variants, uncovering sites of genetic regulation in cis (37,453 methylation quantitative trait loci and 12,505 expression quantitative trait loci) and 13,747 DNA methylation loci affecting gene expression, with over one-third specific to the retina. Methylation and expression quantitative trait loci show non-random distribution and enrichment of biological processes related to synapse, mitochondria, and catabolism. Summary data-based Mendelian randomization and colocalization analyses identify 87 target genes where methylation and gene-expression changes likely mediate the genotype effect on age-related macular degeneration. Integrated pathway analysis reveals epigenetic regulation of immune response and metabolism including the glutathione pathway and glycolysis. Our study thus defines key roles of genetic variations driving methylation changes, prioritizes epigenetic control of gene expression, and suggests frameworks for regulation of macular degeneration pathology by genotype-environment interaction in retina.

Is developmental plasticity triggered by DNA methylation changes in the invasive cane toad (Rhinella marina)?

Many organisms can adjust their development according to environmental conditions, including the presence of conspecifics. Although this developmental plasticity is common in amphibians, its underlying molecular mechanisms remain largely unknown. Exposure during development to either 'cannibal cues' from older conspecifics, or 'alarm cues' from injured conspecifics, causes reduced growth and survival in cane toad (Rhinella marina) tadpoles. Epigenetic modifications, such as changes in DNA methylation patterns, are a plausible mechanism underlying these developmental plastic responses. Here we tested this hypothesis, and asked whether cannibal cues and alarm cues trigger the same DNA methylation changes in developing cane toads. We found that exposure to both cannibal cues and alarm cues was associated with local changes in DNA methylation patterns. These DNA methylation changes affected genes putatively involved in developmental processes, but in different genomic regions for different conspecific-derived cues. Genetic background explains most of the epigenetic variation among individuals. Overall, the molecular mechanisms triggered by exposure to cannibal cues seem to differ from those triggered by alarm cues. Studies linking epigenetic modifications to transcriptional activity are needed to clarify the proximate mechanisms that regulate developmental plasticity in cane toads.

Comparing feature selection and machine learning approaches for predicting CYP2D6 methylation from genetic variation

Introduction

Pharmacogenetics currently supports clinical decision-making on the basis of a limited number of variants in a few genes and may benefit paediatric prescribing where there is a need for more precise dosing. Integrating genomic information such as methylation into pharmacogenetic models holds the potential to improve their accuracy and consequently prescribing decisions. Cytochrome P450 2D6 (CYP2D6) is a highly polymorphic gene conventionally associated with the metabolism of commonly used drugs and endogenous substrates. We thus sought to predict epigenetic loci from single nucleotide polymorphisms (SNPs) related to CYP2D6 in children from the GUSTO cohort.

Methods

Buffy coat DNA methylation was quantified using the Illumina Infinium Methylation EPIC beadchip. CpG sites associated with CYP2D6 were used as outcome variables in Linear Regression, Elastic Net and XGBoost models. We compared feature selection of SNPs from GWAS mQTLs, GTEx eQTLs and SNPs within 2 MB of the CYP2D6 gene and the impact of adding demographic data. The samples were split into training (75%) sets and test (25%) sets for validation. In Elastic Net model and XGBoost models, optimal hyperparameter search was done using 10-fold cross validation. Root Mean Square Error and R-squared values were obtained to investigate each models' performance. When GWAS was performed to determine SNPs associated with CpG sites, a total of 15 SNPs were identified where several SNPs appeared to influence multiple CpG sites.

Results

Overall, Elastic Net models of genetic features appeared to perform marginally better than heritability estimates and substantially better than Linear Regression and XGBoost models. The addition of nongenetic features appeared to improve performance for some but not all feature sets and probes. The best feature set and Machine Learning (ML) approach differed substantially between CpG sites and a number of top variables were identified for each model.

Discussion

The development of SNP-based prediction models for CYP2D6 CpG methylation in Singaporean children of varying ethnicities in this study has clinical application. With further validation, they may add to the set of tools available to improve precision medicine and pharmacogenetics-based dosing.

Dissecting the sequence and structural determinants guiding m6A deposition and evolution via inter- and intra-species hybrids

BACKGROUND

N6-methyladenosine (m6A) is the most abundant mRNA modification, and controls mRNA stability. m6A distribution varies considerably between and within species. Yet, it is unclear to what extent this variability is driven by changes in genetic sequences ('cis') or cellular environments ('trans') and via which mechanisms.

RESULTS

Here we dissect the determinants governing RNA methylation via interspecies and intraspecies hybrids in yeast and mammalian systems, coupled with massively parallel reporter assays and m6A-QTL reanalysis. We find that m6A evolution and variability is driven primarily in 'cis', via two mechanisms: (1) variations altering m6A consensus motifs, and (2) variation impacting mRNA secondary structure. We establish that mutations impacting RNA structure - even when distant from an m6A consensus motif - causally dictate methylation propensity. Finally, we demonstrate that allele-specific differences in m6A levels lead to allele-specific changes in gene expression.

CONCLUSIONS

Our findings define the determinants governing m6A evolution and diversity and characterize the consequences thereof on gene expression regulation.

Increased CpG methylation at the CDH1 locus in inflamed ileal mucosa of patients with Crohn disease

BACKGROUND

E-cadherin, a major actor of cell adhesion in the intestinal barrier, is encoded by the CDH1 gene associated with susceptibility to Crohn Disease (CD) and colorectal cancer. Since epigenetic mechanisms are suspected to contribute to the multifactorial pathogenesis of CD, we studied CpG methylation at the CDH1 locus. The methylation of the CpG island (CGI) and of the 1st enhancer, two critical regulatory positions, was quantified in surgical specimens of inflamed ileal mucosa and in peripheral blood mononuclear cells (PBMC) of 21 CD patients. Sixteen patients operated on for a non-inflammatory bowel disease, although not normal controls, provided a macroscopically normal ileal mucosa and PBMC for comparison.

RESULTS

In ileal mucosa, 19/21 (90%) CD patients vs 8/16 control patients (50%) (p < 0.01) had a methylated CDH1 promoter CGI. In PBMC, CD patients with methylated CGI were 11/21 (52%) vs 7/16 controls (44%), respectively. Methylation in the 1st enhancer of CDH1 was also higher in the CD group for each of the studied CpGs and for their average value (45 ± 17% in CD patients vs 36 ± 17% in controls; p < 0.001). Again, methylation was comparable in PBMC. Methylation of CGI and 1st enhancer were not correlated in mucosa or PBMC.

CONCLUSIONS

Methylation of several CpGs at the CDH1 locus was increased in the inflamed ileal mucosa, not in the PBMC, of CD patients, suggesting the association of CDH1 methylation with ileal inflammation. Longitudinal studies will explore if this increased methylation is a risk marker for colorectal cancer.

Global DNA Methylation Levels Viz-a-Viz Genetic and Biochemical Variations in One Carbon Metabolic Pathway: An Exploratory Study from North India

Despite the importance of one carbon metabolic pathway (OCMP) in modulating the DNA methylation process, only a few population-based studies have explored their relationship among healthy individuals. This study aimed to understand the variations in global DNA methylation levels with respect to selected genetic (CBS 844ins68, MTRR A66G, MTR A2756G, and MTHFR C677T polymorphisms) and biochemical (folate, vitamin B12, and homocysteine) markers associated with OCMP among healthy North Indian adults. The study has been conducted among 1095 individuals of either sex (69.5% females), aged 30-75 years. A sample of 5 mL of blood was collected from each participant. Homocysteine, folate, and vitamin B12 levels were determined using the chemiluminescence technique. Restriction digestion was performed for genotyping MTRR A66G, MTR A2756G, and MTHFR C677T polymorphisms and allele-specific PCR amplification for CBS 844ins68 polymorphism. Global DNA methylation levels were analyzed using ELISA-based colorimetric technique. Of the selected genetic and biochemical markers, the mutant MTRR A66G allele was positively associated with global DNA methylation levels. Further, advanced age was inversely associated with methylation levels. MTRR 66GG genotype group was hypermethylated than other genotypes in folate replete and vitamin B12 deficient group (a condition prevalent among vegetarians), suggesting that the G allele may be more efficient than the wild-type allele in such conditions. Global DNA methylation levels appeared to be more influenced by genetic than biochemical factors. MTRR 66G allele may have a selective advantage in vitamin B12 deficient conditions. Further research should be undertaken to understand how genetics affects epigenetic processes.

Potential Role of DNA Methylation as a Driver of Plastic Responses to the Environment Across Cells, Organisms, and Populations

There is great interest in exploring epigenetic modifications as drivers of adaptive organismal responses to environmental change. Extending this hypothesis to populations, epigenetically driven plasticity could influence phenotypic changes across environments. The canonical model posits that epigenetic modifications alter gene regulation and subsequently impact phenotypes. We first discuss origins of epigenetic variation in nature, which may arise from genetic variation, spontaneous epimutations, epigenetic drift, or variation in epigenetic capacitors. We then review and synthesize literature addressing three facets of the aforementioned model: (i) causal effects of epigenetic modifications on phenotypic plasticity at the organismal level, (ii) divergence of epigenetic patterns in natural populations distributed across environmental gradients, and (iii) the relationship between environmentally induced epigenetic changes and gene expression at the molecular level. We focus on DNA methylation, the most extensively studied epigenetic modification. We find support for environmentally associated epigenetic structure in populations and selection on stable epigenetic variants, and that inhibition of epigenetic enzymes frequently bears causal effects on plasticity. However, there are pervasive confounding issues in the literature. Effects of chromatin-modifying enzymes on phenotype may be independent of epigenetic marks, alternatively resulting from functions and protein interactions extrinsic of epigenetics. Associations between environmentally induced changes in DNA methylation and expression are strong in plants and mammals but notably absent in invertebrates and nonmammalian vertebrates. Given these challenges, we describe emerging approaches to better investigate how epigenetic modifications affect gene regulation, phenotypic plasticity, and divergence among populations.

Epigenetic Methylation Changes in Pregnant Women: Bisphenol Exposure and Atopic Dermatitis

Bisphenol is a chemical substance widely used in plastic products and food containers. In this study, we observed a relationship between DNA methylation and atopic dermatitis (AD) in the peripheral blood mononuclear cells (PBMCs) of pregnant women exposed to bisphenol A (BPA) and its alternatives, bisphenol S (BPS) and bisphenol F (BPF). DNA methylation is an epigenetic mechanism that regulates gene expression, which can be altered by environmental factors, and affects the onset and progression of diseases. We found that genes belonging to the JAK-STAT and PI3K-AKT signaling pathways were hypomethylated in the blood of pregnant women exposed to bisphenols. These genes play important roles in skin barrier function and immune responses, and may influence AD. Therefore, we suggest that not only BPA, but also BPS and BPF, which are used as alternatives, can have a negative impact on AD through epigenetic mechanisms.

Metastable epialleles in humans

First identified in isogenic mice, metastable epialleles (MEs) are loci where the extent of DNA methylation (DNAm) is variable between individuals but correlates across tissues derived from different germ layers within a given individual. This property, termed systemic interindividual variation (SIV), is attributed to stochastic methylation establishment before germ layer differentiation. Evidence suggests that some putative human MEs are sensitive to environmental exposures in early development. In this review we introduce key concepts pertaining to human MEs, describe methods used to identify MEs in humans, and review their genomic features. We also highlight studies linking DNAm at putative human MEs to early environmental exposures and postnatal (including disease) phenotypes.

Somatic mutation as an explanation for epigenetic aging

DNA methylation marks have recently been used to build models known as "epigenetic clocks" which predict calendar age. As methylation of cytosine promotes C-to-T mutations, we hypothesized that the methylation changes observed with age should reflect the accrual of somatic mutations, and the two should yield analogous aging estimates. In analysis of multimodal data from 9,331 human individuals, we find that CpG mutations indeed coincide with changes in methylation, not only at the mutated site but also with pervasive remodeling of the methylome out to ±10 kilobases. This one-to-many mapping enables mutation-based predictions of age that agree with epigenetic clocks, including which individuals are aging faster or slower than expected. Moreover, genomic loci where mutations accumulate with age also tend to have methylation patterns that are especially predictive of age. These results suggest a close coupling between the accumulation of sporadic somatic mutations and the widespread changes in methylation observed over the course of life.

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

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

Pleiotropic influence of DNA methylation QTLs on physiological and ageing traits

DNA methylation is influenced by genetic and non-genetic factors. Here, we chart quantitative trait loci (QTLs) that modulate levels of methylation at highly conserved CpGs using liver methylome data from mouse strains belonging to the BXD family. A regulatory hotspot on chromosome 5 had the highest density of trans-acting methylation QTLs (trans-meQTLs) associated with multiple distant CpGs. We refer to this locus as meQTL.5a. Trans-modulated CpGs showed age-dependent changes and were enriched in developmental genes, including several members of the MODY pathway (maturity onset diabetes of the young). The joint modulation by genotype and ageing resulted in a more 'aged methylome' for BXD strains that inherited the DBA/2J parental allele at meQTL.5a. Further, several gene expression traits, body weight, and lipid levels mapped to meQTL.5a, and there was a modest linkage with lifespan. DNA binding motif and protein-protein interaction enrichment analyses identified the hepatic nuclear factor, Hnf1a (MODY3 gene in humans), as a strong candidate. The pleiotropic effects of meQTL.5a could contribute to variations in body size and metabolic traits, and influence CpG methylation and epigenetic ageing that could have an impact on lifespan.

Genetics and Epigenetics in Obesity: What Do We Know so Far?

PURPOSE OF REVIEW

Enormous progress has been made in understanding the genetic architecture of obesity and the correlation of epigenetic marks with obesity and related traits. This review highlights current research and its challenges in genetics and epigenetics of obesity.

RECENT FINDINGS

Recent progress in genetics of polygenic traits, particularly represented by genome-wide association studies, led to the discovery of hundreds of genetic variants associated with obesity, which allows constructing polygenic risk scores (PGS). In addition, epigenome-wide association studies helped identifying novel targets and methylation sites being important in the pathophysiology of obesity and which are essential for the generation of methylation risk scores (MRS). Despite their great potential for predicting the individual risk for obesity, the use of PGS and MRS remains challenging. Future research will likely discover more loci being involved in obesity, which will contribute to better understanding of the complex etiology of human obesity. The ultimate goal from a clinical perspective will be generating highly robust and accurate prediction scores allowing clinicians to predict obesity as well as individual responses to body weight loss-specific life-style interventions.

Methyl-TWAS: A powerful method for in silico transcriptome-wide association studies (TWAS) using long-range DNA methylation

In silico transcriptome-wide association studies (TWAS) are commonly used to test whether expression of specific genes is linked to a complex trait. However, genotype-based in silico TWAS such as PrediXcan, exhibit low prediction accuracy for a majority of genes because genotypic data lack tissue- and disease-specificity and are not affected by the environment. Because methylation is tissue-specific and, like gene expression, can be modified by environment or disease status, methylation should predict gene expression with more accuracy than SNPs. Therefore, we propose Methyl-TWAS, the first approach that utilizes long-range methylation markers to impute gene expression for in silico TWAS through penalized regression. Methyl-TWAS 1) predicts epigenetically regulated/associated expression (eGReX), which incorporates tissue-specific expression and both genetically- (GReX) and environmentally-regulated expression to identify differentially expressed genes (DEGs) that could not be identified by genotype-based methods; and 2) incorporates both cis- and trans- CpGs, including various regulatory regions to identify DEGs that would be missed using cis- methylation only. Methyl-TWAS outperforms PrediXcan and two other methods in imputing gene expression in the nasal epithelium, particularly for immunity-related genes and DEGs in atopic asthma. Methyl-TWAS identified 3,681 (85.2%) of the 4,316 DEGs identified in a previous TWAS of atopic asthma using measured expression, while PrediXcan could not identify any gene. Methyl-TWAS also outperforms PrediXcan for expression imputation as well as in silico TWAS in white blood cells. Methyl-TWAS is a valuable tool for in silico TWAS, leveraging a growing body of publicly available genome-wide DNA methylation data for a variety of human tissues.

DNA methylation sites in early adulthood characterised by pubertal timing and development: a twin study

BACKGROUND

Puberty is a highly heritable and variable trait, with environmental factors having a role in its eventual timing and development. Early and late pubertal onset are both associated with various diseases developing later in life, and epigenetic characterisation of pubertal timing and development could lead to important insights. Blood DNA methylation, reacting to both genotype and environment, has been associated with puberty; however, such studies are relatively scarce. We investigated peripheral blood DNA methylation profiles (using Illumina 450 K and EPIC platforms) of 1539 young adult Finnish twins associated with pubertal development scale (PDS) at ages 12 and 14 as well as pubertal age (PA).

RESULTS

Fixed effect meta-analysis of the two platforms on 347,521 CpGs in common identified 58 CpG sites associated (p < 1 × 10-5) with either PDS or PA. All four CpGs associated with PA and 45 CpGs associated with PDS were sex-specific. Thirteen CpGs had a high heritability (h2: 0.51-0.98), while one CpG site (mapped to GET4) had a high shared environmental component accounting for 68% of the overall variance in methylation at the site. Utilising twin discordance analysis, we found 6 CpG sites (5 associated with PDS and 1 with PA) that had an environmentally driven association with puberty. Furthermore, genes with PDS- or PA-associated CpGs were consistently linked to various developmental processes and diseases such as breast, prostate and ovarian cancer, while methylation quantitative trait loci of associated CpG sites were enriched in immune pathways developing during puberty.

CONCLUSIONS

By identifying puberty-associated DNA methylation sites and examining the effects of sex, environment and genetics, we shed light on the intricate interplay between environment and genetics in the context of puberty. Through our comprehensive analysis, we not only deepen the understanding of the significance of both genetic and environmental factors in the complex processes of puberty and its timing, but also gain insights into potential links with disease risks.

Epigenetics insights from perceived facial aging

Facial aging is the most visible manifestation of aging. People desire to look younger than others of the same chronological age. Hence, perceived age is often used as a visible marker of aging, while biological age, often estimated by methylation markers, is used as an objective measure of age. Multiple epigenetics-based clocks have been developed for accurate estimation of general biological age and the age of specific organs, including the skin. However, it is not clear whether the epigenetic biomarkers (CpGs) used in these clocks are drivers of aging processes or consequences of aging. In this proof-of-concept study, we integrate data from GWAS on perceived facial aging and EWAS on CpGs measured in blood. By running EW Mendelian randomization, we identify hundreds of putative CpGs that are potentially causal to perceived facial aging with similar numbers of damaging markers that causally drive or accelerate facial aging and protective methylation markers that causally slow down or protect from aging. We further demonstrate that while candidate causal CpGs have little overlap with known epigenetics-based clocks, they affect genes or proteins with known functions in skin aging, such as skin pigmentation, elastin, and collagen levels. Overall, our results suggest that blood methylation markers reflect facial aging processes, and thus can be used to quantify skin aging and develop anti-aging solutions that target the root causes of aging.

Investigation of the impact of biogeographic ancestry on DNA methylation based age predictions comparing a Middle East and a Central European population

DNA methylation based age prediction is a new method in the toolbox of forensic genetics. Typically, the method is applied in the course of police investigation e.g. to predict the age of an unknown person that has left a biological trace at a crime scene. The method can also be used to answer other forensic questions, for example to estimate the age of unknown human bodies in the course of the identification process. In the present study, we tested for a potential impact of biogeographic ancestry (BGA) on age predictions using five age dependent methylated CpG sites within the genetic regions of ELOVL2, MIR29B2CHG, FHL2, KLF14 and TRIM59. We collected 102 blood samples each from donors living in Iraq, Middle East (ME) and Germany, Central Europe (EU). Both sample sets were matched in sex and age ranging from 18 to 68 years with exactly one male and female sample per year of age. All samples were analyzed by bisulfite pyrosequencing applying a multiplex pre-amplification strategy based on a single input of 35 ng converted DNA in the PCR. For the CpGs in MIR29B2CHG, FHL2 and KLF14, we observed significantly different methylation levels between the two populations. While we were able to train two highly accurate prediction models for the respective population with mean absolute deviations between predicted and actual ages (MAD) of 3.34 years for the ME model, and 2.72 years for the EU model, we found an absolute prediction difference between the two population specific models of more than 4 years. A combined model for both populations compensated the methylation difference between the two populations, providing MADs of prediction of only 3.81 years for ME and 3.31 years for EU samples. In total, the results of the present study strongly support the benefit of BGA information for more reliable methylation based age predictions.

Parental Exposure to Environmentally Relevant Concentrations of Bisphenol-A Bis(diphenyl phosphate) Impairs Vascular Development in Offspring through DNA/RNA Methylation-Dependent Transmission

Bisphenol-A bis(diphenyl phosphate) (BDP) has been increasingly detected in indoor environmental and human samples. Little is known about its developmental toxicity, particularly the intergenerational effects of parental exposure. In this study, adult zebrafish were exposed to BDP at 30-30,000 ng/L for 28 days, with results showing that exposure did not cause a transfer of BDP or its metabolites to offspring. Vascular morphometric profiling revealed that parental exposure to BDP at 30 and 300 ng/L exerted significant effects on the vascular development of offspring, encompassing diverse alterations in multiple types of blood vessels. N6-Methyladenosine (m6A) methylated RNA immunoprecipitation sequencing of larvae in the 300 ng/L group revealed 378 hypomethylated and 350 hypermethylated m6A peaks that were identified in mRNA transcripts of genes crucial for vascular development, including the Notch/Vegf signaling pathway. Concomitant changes in 5 methylcytosine (m5C) DNA methylation and gene expression of m6A modulators (alkbh5, kiaa1429, and ythdf1) were observed in both parental gonads and offspring exposed to BDP. These results reveal that parental exposure to low concentrations of BDP caused offspring vascular disorders by interfering with DNA and RNA methylation, uncovering a unique DNA-RNA modification pattern in the intergenerational transmission of BDP's developmental toxicity.

Genome-wide DNA methylation landscape of four Chinese populations and epigenetic variation linked to Tibetan high-altitude adaptation

DNA methylation (DNAm) is one of the major epigenetic mechanisms in humans and is important in diverse cellular processes. The variation of DNAm in the human population is related to both genetic and environmental factors. However, the DNAm profiles have not been investigated in the Chinese population of diverse ethnicities. Here, we performed double-strand bisulfite sequencing (DSBS) for 32 Chinese individuals representing four major ethnic groups including Han Chinese, Tibetan, Zhuang, and Mongolian. We identified a total of 604,649 SNPs and quantified DNAm at more than 14 million CpGs in the population. We found global DNAm-based epigenetic structure is different from the genetic structure of the population, and ethnic difference only partially explains the variation of DNAm. Surprisingly, non-ethnic-specific DNAm variations showed stronger correlation with the global genetic divergence than these ethnic-specific DNAm. Differentially methylated regions (DMRs) among these ethnic groups were found around genes in diverse biological processes. Especially, these DMR-genes between Tibetan and non-Tibetans were enriched around high-altitude genes including EPAS1 and EGLN1, suggesting DNAm alteration plays an important role in high-altitude adaptation. Our results provide the first batch of epigenetic maps for Chinese populations and the first evidence of the association of epigenetic changes with Tibetans' high-altitude adaptation.

O’Brien T, Perucca P, Ge Z, Anderson A. Applications for Deep Learning in Epilepsy Genetic Research

Epilepsy is a group of brain disorders characterised by an enduring predisposition to generate unprovoked seizures. Fuelled by advances in sequencing technologies and computational approaches, more than 900 genes have now been implicated in epilepsy. The development and optimisation of tools and methods for analysing the vast quantity of genomic data is a rapidly evolving area of research. Deep learning (DL) is a subset of machine learning (ML) that brings opportunity for novel investigative strategies that can be harnessed to gain new insights into the genomic risk of people with epilepsy. DL is being harnessed to address limitations in accuracy of long-read sequencing technologies, which improve on short-read methods. Tools that predict the functional consequence of genetic variation can represent breaking ground in addressing critical knowledge gaps, while methods that integrate independent but complimentary data enhance the predictive power of genetic data. We provide an overview of these DL tools and discuss how they may be applied to the analysis of genetic data for epilepsy research.

Cis-meQTL for cocaine use-associated DNA methylation in an HIV-positive cohort show pleiotropic effects on multiple traits

BACKGROUND

Cocaine use (CU) is associated with psychiatric and medical diseases. Little is known about the mechanisms of CU-related comorbidities. Findings from preclinical and clinical studies have suggested that CU is associated with aberrant DNA methylation (DNAm) that may be influenced by genetic variants [i.e., methylation quantitative trait loci (meQTLs)]. In this study, we mapped cis-meQTLs for CU-associated DNAm sites (CpGs) in an HIV-positive cohort (Ntotal = 811) and extended the meQTLs to multiple traits.

RESULTS

We conducted cis-meQTL analysis for 224 candidate CpGs selected for their association with CU in blood. We identified 7,101 significant meQTLs [false discovery rate (FDR) < 0.05], which mostly mapped to genes involved in immunological functions and were enriched in immune pathways. We followed up the meQTLs using phenome-wide association study and trait enrichment analyses, which revealed 9 significant traits. We tested for causal effects of CU on these 9 traits using Mendelian Randomization and found evidence that CU plays a causal role in increasing hypertension (p-value = 2.35E-08) and decreasing heel bone mineral density (p-value = 1.92E-19).

CONCLUSIONS

These findings suggest that genetic variants for CU-associated DNAm have pleiotropic effects on other relevant traits and provide new insights into the causal relationships between cocaine use and these complex traits.

PCSK9 Involvement in Autism Etiology: Sequence Variations, Protein Concentration, and Promoter Methylation

BACKGROUND

Besides its main role in the control of blood cholesterol, PCSK9 has a role in the regulation of neuronal development and apoptosis. We suggest, for the first time, the possible involvement of PCSK9 in autism.

METHOD

In this case-control study, Sanger sequencing was used to analyze sequence variations in the PCSK9 gene exons and their flanking intronic sequences. ELISA assay was used to determine the plasma concentration of PCSK9. The methylation percentage of the PCSK9 gene promoter was assessed by methylation-specific PCR (MSP).

RESULTS

Forty-three variants were found; out of them, seven showed differential frequency between patients and controls. rs.45448095, rs.45613943, rs.630431, rs.529500286, and rs.45439391 are risk factors for autism, while rs.11800231 and rs.483462 are protective variants. The concentration of plasma PCSK9 protein was significantly elevated and the methylation percentage of PCSK9 gene promoter was significantly lower in cases than in controls (p <0.001 and = 0.002, respectively). ROC curve analysis identified an area under the curve (AUC) of 0.915 for plasma protein concentration and 0.693 for percent gene promoter methylation. In addition, two new variants were identified (g.23809C>T in intron 11 and g.24071T>G in 3' UTR).

CONCLUSION

This is the first study to investigate the correlation between PCSK9 protein and autism and suggests the potential involvement of PCSK9 as one of the susceptibility genes for autism. Further studies with a larger number of subjects are recommended.

Whole blood DNA methylation changes are associated with anti-TNF drug concentration in patients with Crohn's disease

BACKGROUND AND AIMS

Anti-TNF treatment failure in patients with inflammatory bowel disease (IBD) is common and frequently related to low drug concentrations. In order to identify patients who may benefit from dose optimisation at the outset of anti-TNF therapy, we sought to define epigenetic biomarkers in whole blood at baseline associated with anti-TNF drug concentrations at week 14.

METHODS

DNA methylation from 1,104 whole blood samples from 385 patients in the Personalised Anti-TNF Therapy in Crohn's disease (PANTS) study were assessed using the Illumina EPIC Beadchip (v1.0) at baseline, weeks 14, 30 and 54. We compared DNA methylation profiles in anti-TNF-treated patients who experienced primary non-response at week 14 and if they were assessed at subsequent time points, were not in remission at week 30 or 54 (infliximab n = 99, adalimumab n = 94), with patients who responded at week 14 and when assessed at subsequent time points, were in remission at week 30 or 54 (infliximab n = 99, adalimumab n = 93).

RESULTS

Overall, between baseline and week 14, we observed 4,999 differentially methylated probes (DMPs) annotated to 2376 genes following anti-TNF treatment. Pathway analysis identified 108 significant gene ontology terms enriched in biological processes related to immune system processes and responses.Epigenome-wide association (EWAS) analysis identified 323 DMPs annotated to 210 genes at baseline associated with higher anti-TNF drug concentrations at week 14. Of these, 125 DMPs demonstrated shared associations with other common traits (proportion of shared CpGs compared to DMPs) including body mass index (23.2%), followed by CRP (11.5%), smoking (7.4%), alcohol consumption per day (7.1%) and IBD type (6.8%). EWAS of primary non-response to anti-TNF identified 20 DMPs that were associated with both anti-TNF drug concentration and primary non-response to anti-TNF with a strong correlation of the coefficients (Spearman's rho = -0.94, p < 0.001).

CONCLUSION

Baseline DNA methylation profiles may be used as a predictor for anti-TNF drug concentration at week 14 to identify patients who may benefit from dose optimisation at the outset of anti-TNF therapy.

Genetic impacts on DNA methylation help elucidate regulatory genomic processes

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Stability selection enhances feature selection and enables accurate prediction of gestational age using only five DNA methylation sites

BACKGROUND

DNA methylation (DNAm) is robustly associated with chronological age in children and adults, and gestational age (GA) in newborns. This property has enabled the development of several epigenetic clocks that can accurately predict chronological age and GA. However, the lack of overlap in predictive CpGs across different epigenetic clocks remains elusive. Our main aim was therefore to identify and characterize CpGs that are stably predictive of GA.

RESULTS

We applied a statistical approach called 'stability selection' to DNAm data from 2138 newborns in the Norwegian Mother, Father, and Child Cohort study. Stability selection combines subsampling with variable selection to restrict the number of false discoveries in the set of selected variables. Twenty-four CpGs were identified as being stably predictive of GA. Intriguingly, only up to 10% of the CpGs in previous GA clocks were found to be stably selected. Based on these results, we used generalized additive model regression to develop a new GA clock consisting of only five CpGs, which showed a similar predictive performance as previous GA clocks (R2 = 0.674, median absolute deviation = 4.4 days). These CpGs were in or near genes and regulatory regions involved in immune responses, metabolism, and developmental processes. Furthermore, accounting for nonlinear associations improved prediction performance in preterm newborns.

CONCLUSION

We present a methodological framework for feature selection that is broadly applicable to any trait that can be predicted from DNAm data. We demonstrate its utility by identifying CpGs that are highly predictive of GA and present a new and highly performant GA clock based on only five CpGs that is more amenable to a clinical setting.

QTL mapping of human retina DNA methylation identifies 87 gene-epigenome interactions in age-related macular degeneration

DNA methylation (DNAm) provides a crucial epigenetic mark linking genetic variations to environmental influence. We analyzed array-based DNAm profiles of 160 human retinas with co-measured RNA-seq and > 8 million genetic variants, uncovering sites of genetic regulation in cis (37,453 mQTLs and 12,505 eQTLs) and 13,747 eQTMs (DNAm loci affecting gene expression), with over one-third specific to the retina. mQTLs and eQTMs show non-random distribution and enrichment of biological processes related to synapse, mitochondria, and catabolism. Summary data-based Mendelian randomization and colocalization analyses identify 87 target genes where methylation and gene-expression changes likely mediate the genotype effect on age-related macular degeneration (AMD). Integrated pathway analysis reveals epigenetic regulation of immune response and metabolism including the glutathione pathway and glycolysis. Our study thus defines key roles of genetic variations driving methylation changes, prioritizes epigenetic control of gene expression, and suggests frameworks for regulation of AMD pathology by genotype-environment interaction in retina.

Underexplored reciprocity between genome-wide methylation status and long non-coding RNA expression reflected in breast cancer research: potential impacts for the disease management in the framework of 3P medicine

Breast cancer (BC) is the most common female malignancy reaching a pandemic scale worldwide. A comprehensive interplay between genetic alterations and shifted epigenetic regions synergistically leads to disease development and progression into metastatic BC. DNA and histones methylations, as the most studied epigenetic modifications, represent frequent and early events in the process of carcinogenesis. To this end, long non-coding RNAs (lncRNAs) are recognized as potent epigenetic modulators in pathomechanisms of BC by contributing to the regulation of DNA, RNA, and histones' methylation. In turn, the methylation status of DNA, RNA, and histones can affect the level of lncRNAs expression demonstrating the reciprocity of mechanisms involved. Furthermore, lncRNAs might undergo methylation in response to actual medical conditions such as tumor development and treated malignancies. The reciprocity between genome-wide methylation status and long non-coding RNA expression levels in BC remains largely unexplored. Since the bio/medical research in the area is, per evidence, strongly fragmented, the relevance of this reciprocity for BC development and progression has not yet been systematically analyzed. Contextually, the article aims at:consolidating the accumulated knowledge on both-the genome-wide methylation status and corresponding lncRNA expression patterns in BC andhighlighting the potential benefits of this consolidated multi-professional approach for advanced BC management. Based on a big data analysis and machine learning for individualized data interpretation, the proposed approach demonstrates a great potential to promote predictive diagnostics and targeted prevention in the cost-effective primary healthcare (sub-optimal health conditions and protection against the health-to-disease transition) as well as advanced treatment algorithms tailored to the individualized patient profiles in secondary BC care (effective protection against metastatic disease). Clinically relevant examples are provided, including mitochondrial health control and epigenetic regulatory mechanisms involved.

Variation in histone configurations correlates with gene expression across nine inbred strains of mice

The Diversity Outbred (DO) mice and their inbred founders are widely used models of human disease. However, although the genetic diversity of these mice has been well documented, their epigenetic diversity has not. Epigenetic modifications, such as histone modifications and DNA methylation, are important regulators of gene expression and, as such, are a critical mechanistic link between genotype and phenotype. Therefore, creating a map of epigenetic modifications in the DO mice and their founders is an important step toward understanding mechanisms of gene regulation and the link to disease in this widely used resource. To this end, we performed a strain survey of epigenetic modifications in hepatocytes of the DO founders. We surveyed four histone modifications (H3K4me1, H3K4me3, H3K27me3, and H3K27ac), as well as DNA methylation. We used ChromHMM to identify 14 chromatin states, each of which represents a distinct combination of the four histone modifications. We found that the epigenetic landscape is highly variable across the DO founders and is associated with variation in gene expression across strains. We found that epigenetic state imputed into a population of DO mice recapitulated the association with gene expression seen in the founders, suggesting that both histone modifications and DNA methylation are highly heritable mechanisms of gene expression regulation. We illustrate how DO gene expression can be aligned with inbred epigenetic states to identify putative cis-regulatory regions. Finally, we provide a data resource that documents strain-specific variation in the chromatin state and DNA methylation in hepatocytes across nine widely used strains of laboratory mice.

Melatonin: A Potential Regulator of DNA Methylation

The pineal gland-derived indoleamine hormone, melatonin, regulates multiple cellular processes, ranging from chronobiology, proliferation, apoptosis, and oxidative damage to pigmentation, immune regulation, and mitochondrial metabolism. While melatonin is best known as a master regulator of the circadian rhythm, previous studies also have revealed connections between circadian cycle disruption and genomic instability, including epigenetic changes in the pattern of DNA methylation. For example, melatonin secretion is associated with differential circadian gene methylation in night shift workers and the regulation of genomic methylation during embryonic development, and there is accumulating evidence that melatonin can modify DNA methylation. Since the latter one impacts cancer initiation, and also, non-malignant diseases development, and that targeting DNA methylation has become a novel intervention target in clinical therapy, this review discusses the potential role of melatonin as an under-investigated candidate epigenetic regulator, namely by modulating DNA methylation via changes in mRNA and the protein expression of DNA methyltransferases (DNMTs) and ten-eleven translocation (TET) proteins. Furthermore, since melatonin may impact changes in the DNA methylation pattern, the authors of the review suggest its possible use in combination therapy with epigenetic drugs as a new anticancer strategy.

Osteoarthritis genetic risk acting on the galactosyltransferase gene COLGALT2 has opposing functional effects in articulating joint tissues

BACKGROUND

Investigation of cartilage and chondrocytes has revealed that the osteoarthritis risk marked by the independent DNA variants rs11583641 and rs1046934 mediate  their effects by decreasing the methylation status of CpG dinucleotides in enhancers and increasing the expression of shared target gene COLGALT2. We set out to investigate if these functional effects operate in a non-cartilaginous joint tissue.

METHODS

Nucleic acids were extracted from the synovium of osteoarthritis patients. Samples were genotyped, and DNA methylation was quantified by pyrosequencing at CpGs within the COLGALT2 enhancers. CpGs were tested for enhancer effects using a synovial cell line and a reporter gene assay. DNA methylation was altered using epigenetic editing, with the impact on gene expression determined using quantitative polymerase chain reaction. In silico analysis complemented laboratory experiments.

RESULTS

The rs1046934 genotype did not associate with DNA methylation or COLGALT2 expression in the synovium, whereas the rs11583641 genotype did. Surprisingly, the effects for rs11583641 were opposite to those previously observed in cartilage. Epigenetic editing in synovial cells revealed that enhancer methylation is causally linked to COLGALT2 expression.

CONCLUSIONS

This is the first direct demonstration for osteoarthritis genetic risk of a functional link between DNA methylation and gene expression operating in opposite directions between articular joint tissues. It highlights pleiotropy in the action of osteoarthritis risk and provides a cautionary note in the application of future genetically based osteoarthritis therapies: an intervention that decreases the detrimental effect of a risk allele in one joint tissue may inadvertently increase its detrimental effect in another joint tissue.

meQTL mapping in the GENOA study reveals genetic determinants of DNA methylation in African Americans

Identifying genetic variants that are associated with variation in DNA methylation, an analysis commonly referred to as methylation quantitative trait locus (meQTL) mapping, is an important first step towards understanding the genetic architecture underlying epigenetic variation. Most existing meQTL mapping studies have focused on individuals of European ancestry and are underrepresented in other populations, with a particular absence of large studies in populations with African ancestry. We fill this critical knowledge gap by performing a large-scale cis-meQTL mapping study in 961 African Americans from the Genetic Epidemiology Network of Arteriopathy (GENOA) study. We identify a total of 4,565,687 cis-acting meQTLs in 320,965 meCpGs. We find that 45% of meCpGs harbor multiple independent meQTLs, suggesting potential polygenic genetic architecture underlying methylation variation. A large percentage of the cis-meQTLs also colocalize with cis-expression QTLs (eQTLs) in the same population. Importantly, the identified cis-meQTLs explain a substantial proportion (median = 24.6%) of methylation variation. In addition, the cis-meQTL associated CpG sites mediate a substantial proportion (median = 24.9%) of SNP effects underlying gene expression. Overall, our results represent an important step toward revealing the co-regulation of methylation and gene expression, facilitating the functional interpretation of epigenetic and gene regulation underlying common diseases in African Americans.

The Epigenetics of Animal Personality

Animal personality, consistent individual differences in behaviour, is an important concept for understanding how individuals vary in how they cope with environmental challenges. In order to understand the evolutionary significance of animal personality, it is crucial to understand the underlying regulatory mechanisms. Epigenetic marks such as DNA methylation are hypothesised to play a major role in explaining variation in phenotypic changes in response to environmental alterations. Several characteristics of DNA methylation also align well with the concept of animal personality. In this review paper, we summarise the current literature on the role that molecular epigenetic mechanisms may have in explaining personality variation. We elaborate on the potential for epigenetic mechanisms to explain behavioural variation, behavioural development and temporal consistency in behaviour. We then suggest future routes for this emerging field and point to potential pitfalls that may be encountered. We conclude that a more inclusive approach is needed for studying the epigenetics of animal personality and that epigenetic mechanisms cannot be studied without considering the genetic background.

Analysis of Y-STR diversity and DNA methylation variation among Black and Indian males from KwaZulu-Natal, South Africa

Y-chromosome short tandem repeats (Y-STRs) are essential in understanding genetic structure and diversity of human populations and, most importantly, in identification of male perpetrators in criminal investigations. DNA methylation differences have been reported in human populations and methylation pattern at the CpG sites found within or flanking the Y-STR sites could also aid in human identification. Studies based on DNA methylation (DNAm) at Y-STRs are currently limited. The current study aimed to analyze the Y-STR diversity in South African Black and Indian individuals living in KwaZulu-Natal, Durban, South Africa, with the Yfiler™ Plus Kit and to analyze DNAm patterns in Y-STR markers CpG sites. DNA from 247 stored saliva samples were isolated and quantified. Across the 27 Y-STR loci in the Yfiler™ Plus Kit, 253 alleles were observed in 113 South African Black and Indian males, 112 unique haplotypes were observed, and one haplotype appeared twice (two Black individuals). No statistically significant differences were observed in the genetic diversity between the two population groups (Fst = 0.028, p-value ≥ 0.05). The kit showed a high discrimination capacity (DC) of 0.9912 and an overall haplotype diversity (HD) = 0.9995 among the sampled population groups. DYS438 and DYS448 markers displayed 2 and 3 CpG sites, respectively. Based on the two-tailed Fisher's Exact test, there were no statistically significant differences in the DNAm levels at DYS438 CpGs of Black and Indian males (p > 0.05). The Yfiler™ Plus Kit can be considered highly discriminatory among South African Black and Indian males. Studies on the South African population using Yfiler™ Plus Kit are scarce. Hence, accumulating Y-STR data on the diverse South African population will enhance the representation of South Africa in STR databases. Knowing which Y-STR markers are significantly informative for South Africa is essential for developing Y-STR kits better suited for the different ethnic groups. And to the best of our knowledge, DNA methylation analysis in Y-STR for different ethnic groups has never been done before. Complementing Y-STR data with methylation knowledge could provide population-specific information for forensic identification.

Maternal Immune Activation Induces Adolescent Cognitive Deficits Preceded by Developmental Perturbations in Cortical Reelin Signalling

Exposure to maternal immune activation (MIA) in utero significantly elevates the risk of developing schizophrenia and other neurodevelopmental disorders. To understand the biological mechanisms underlying the link between MIA and increased risk, preclinical animal models have focussed on specific signalling pathways in the brain that mediate symptoms associated with neurodevelopmental disorders such as cognitive dysfunction. Reelin signalling in multiple brain regions is involved in neuronal migration, synaptic plasticity and long-term potentiation, and has been implicated in cognitive deficits. However, how regulation of Reelin expression is affected by MIA across cortical development and associated cognitive functions remains largely unclear. Using a MIA rat model, here we demonstrate cognitive deficits in adolescent object-location memory in MIA offspring and reductions in Reln expression prenatally and in the adult prefrontal cortex. Further, developmental disturbances in gene/protein expression and DNA methylation of downstream signalling components occurred subsequent to MIA-induced Reelin dysregulation and prior to cognitive deficits. We propose that MIA-induced dysregulation of Reelin signalling contributes to the emergence of prefrontal cortex-mediated cognitive deficits through altered NMDA receptor function, resulting in inefficient long-term potentiation. Our data suggest a developmental window during which attenuation of Reelin signalling may provide a possible therapeutic target.

Aging Hallmarks and the Role of Oxidative Stress

Aging is a complex biological process accompanied by a progressive decline in the physical function of the organism and an increased risk of age-related chronic diseases such as cardiovascular diseases, cancer, and neurodegenerative diseases. Studies have established that there exist nine hallmarks of the aging process, including (i) telomere shortening, (ii) genomic instability, (iii) epigenetic modifications, (iv) mitochondrial dysfunction, (v) loss of proteostasis, (vi) dysregulated nutrient sensing, (vii) stem cell exhaustion, (viii) cellular senescence, and (ix) altered cellular communication. All these alterations have been linked to sustained systemic inflammation, and these mechanisms contribute to the aging process in timing not clearly determined yet. Nevertheless, mitochondrial dysfunction is one of the most important mechanisms contributing to the aging process. Mitochondria is the primary endogenous source of reactive oxygen species (ROS). During the aging process, there is a decline in ATP production and elevated ROS production together with a decline in the antioxidant defense. Elevated ROS levels can cause oxidative stress and severe damage to the cell, organelle membranes, DNA, lipids, and proteins. This damage contributes to the aging phenotype. In this review, we summarize recent advances in the mechanisms of aging with an emphasis on mitochondrial dysfunction and ROS production.