Genetic and epigenetic regulation of gene expression in fetal and adult human livers

Animal and plant protein intake during infancy and childhood DNA methylation: a meta-analysis in the NutriPROGRAM consortium

BACKGROUND

Higher early-life animal protein intake is associated with a higher childhood obesity risk compared to plant protein intake. Differential DNA methylation may represent an underlying mechanism.

METHODS

We analysed associations of infant animal and plant protein intakes with DNA methylation in early (2-6 years, N = 579) and late (7̄-12 years, N = 604) childhood in two studies. Study-specific robust linear regression models adjusted for relevant confounders were run, and then meta-analysed using a fixed-effects model. We also performed sex-stratified meta-analyses. Follow-up analyses included pathway analysis and eQTM look-up.

RESULTS

Infant animal protein intake was not associated with DNA methylation in early childhood, but was associated with late-childhood DNA methylation at cg21300373 (P = 4.27 × 10¯8, MARCHF1) and cg10633363 (P = 1.09 × 10¯7, HOXB9) after FDR correction. Infant plant protein intake was associated with early-childhood DNA methylation at cg25973293 (P = 2.26 × 10-7, C1orf159) and cg15407373 (P = 2.13 × 10-7, MBP) after FDR correction. There was no overlap between the findings from the animal and plant protein analyses. We did not find enriched functional pathways at either time point using CpGs associated with animal and plant protein. These CpGs were not previously associated with childhood gene expression. Sex-stratified meta-analyses showed sex-specific DNA methylation associations for both animal and plant protein intake.

CONCLUSION

Infant animal protein intake was associated with DNA methylation at two CpGs in late childhood. Infant plant protein intake was associated with DNA methylation in early childhood at two CpGs. A potential mediating role of DNA methylation at these CpGs between infant protein intake and health outcomes requires further investigation.

Epigenome-wide association study on the plasma metabolome suggests self-regulation of the glycine and serine pathway through DNA methylation

BACKGROUND

The plasma metabolome reflects the physiological state of various biological processes and can serve as a proxy for disease risk. Plasma metabolite variation, influenced by genetic and epigenetic mechanisms, can also affect the cellular microenvironment and blood cell epigenetics. The interplay between the plasma metabolome and the blood cell epigenome remains elusive. In this study, we performed an epigenome-wide association study (EWAS) of 1183 plasma metabolites in 693 participants from the LifeLines-DEEP cohort and investigated the causal relationships in DNA methylation-metabolite associations using bidirectional Mendelian randomization and mediation analysis.

RESULTS

After rigorously adjusting for potential confounders, including genetics, we identified five robust associations between two plasma metabolites (L-serine and glycine) and three CpG sites located in two independent genomic regions (cg14476101 and cg16246545 in PHGDH and cg02711608 in SLC1A5) at a false discovery rate of less than 0.05. Further analysis revealed a complex bidirectional relationship between plasma glycine/serine levels and DNA methylation. Moreover, we observed a strong mediating role of DNA methylation in the effect of glycine/serine on the expression of their metabolism/transport genes, with the proportion of the mediated effect ranging from 11.8 to 54.3%. This result was also replicated in an independent population-based cohort, the Rotterdam Study. To validate our findings, we conducted in vitro cell studies which confirmed the mediating role of DNA methylation in the regulation of PHGDH gene expression.

CONCLUSIONS

Our findings reveal a potential feedback mechanism in which glycine and serine regulate gene expression through DNA methylation.

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.

A novel approach to risk exposure and epigenetics-the use of multidimensional context to gain insights into the early origins of cardiometabolic and neurocognitive health

BACKGROUND

Each mother-child dyad represents a unique combination of genetic and environmental factors. This constellation of variables impacts the expression of countless genes. Numerous studies have uncovered changes in DNA methylation (DNAm), a form of epigenetic regulation, in offspring related to maternal risk factors. How these changes work together to link maternal-child risks to childhood cardiometabolic and neurocognitive traits remains unknown. This question is a key research priority as such traits predispose to future non-communicable diseases (NCDs). We propose viewing risk and the genome through a multidimensional lens to identify common DNAm patterns shared among diverse risk profiles.

METHODS

We identified multifactorial Maternal Risk Profiles (MRPs) generated from population-based data (n = 15,454, Avon Longitudinal Study of Parents and Children (ALSPAC)). Using cord blood HumanMethylation450 BeadChip data, we identified genome-wide patterns of DNAm that co-vary with these MRPs. We tested the prospective relation of these DNAm patterns (n = 914) to future outcomes using decision tree analysis. We then tested the reproducibility of these patterns in (1) DNAm data at age 7 and 17 years within the same cohort (n = 973 and 974, respectively) and (2) cord DNAm in an independent cohort, the Generation R Study (n = 686).

RESULTS

We identified twenty MRP-related DNAm patterns at birth in ALSPAC. Four were prospectively related to cardiometabolic and/or neurocognitive childhood outcomes. These patterns were replicated in DNAm data from blood collected at later ages. Three of these patterns were externally validated in cord DNAm data in Generation R. Compared to previous literature, DNAm patterns exhibited novel spatial distribution across the genome that intersects with chromatin functional and tissue-specific signatures.

CONCLUSIONS

To our knowledge, we are the first to leverage multifactorial population-wide data to detect patterns of variability in DNAm. This context-based approach decreases biases stemming from overreliance on specific samples or variables. We discovered molecular patterns demonstrating prospective and replicable relations to complex traits. Moreover, results suggest that patterns harbour a genome-wide organisation specific to chromatin regulation and target tissues. These preliminary findings warrant further investigation to better reflect the reality of human context in molecular studies of NCDs.

Preneoplastic liver colonization by 11p15.5 altered mosaic cells in young children with hepatoblastoma

Pediatric liver tumors are very rare tumors with the most common diagnosis being hepatoblastoma. While hepatoblastomas are predominantly sporadic, around 15% of cases develop as part of predisposition syndromes such as Beckwith-Wiedemann (11p15.5 locus altered). Here, we identify mosaic genetic alterations of 11p15.5 locus in the liver of hepatoblastoma patients without a clinical diagnosis of Beckwith-Wiedemann syndrome. We do not retrieve these alterations in children with other types of pediatric liver tumors. We show that mosaic 11p15.5 alterations in liver FFPE sections of hepatoblastoma patients display IGF2 overexpression and H19 downregulation together with an alteration of the liver zonation. Moreover, mosaic livers' microenvironment is enriched in extracellular matrix and angiogenesis. Spatial transcriptomics and single-nucleus RNAseq analyses identify a 60-gene signature in 11p15.5 altered hepatocytes. These data provide insights for 11p15.5 mosaicism detection and its functional consequences during the early steps of carcinogenesis.

The Biology of Stress Intolerance in Patients with Chronic Pain—State of the Art and Future Directions

Stress has been consistently linked to negative impacts on physical and mental health. More specifically, patients with chronic pain experience stress intolerance, which is an exacerbation or occurrence of symptoms in response to any type of stress. The pathophysiological mechanisms underlying this phenomenon remain unsolved. In this state-of-the-art paper, we summarised the role of the autonomic nervous system (ANS) and hypothalamus-pituitary-adrenal (HPA) axis, the two major stress response systems in stress intolerance. We provided insights into such mechanisms based on evidence from clinical studies in both patients with chronic pain, showing dysregulated stress systems, and healthy controls supported by preclinical studies, highlighting the link between these systems and symptoms of stress intolerance. Furthermore, we explored the possible regulating role for (epi)genetic mechanisms influencing the ANS and HPA axis. The link between stress and chronic pain has become an important area of research as it has the potential to inform the development of interventions to improve the quality of life for individuals living with chronic pain. As stress has become a prevalent concern in modern society, understanding the connection between stress, HPA axis, ANS, and chronic health conditions such as chronic pain is crucial to improve public health and well-being.

Distinct genome-wide DNA methylation and gene expression signatures in classical monocytes from African American patients with systemic sclerosis

BACKGROUND

Systemic sclerosis (SSc) is a multisystem autoimmune disorder that has an unclear etiology and disproportionately affects women and African Americans. Despite this, African Americans are dramatically underrepresented in SSc research. Additionally, monocytes show heightened activation in SSc and in African Americans relative to European Americans. In this study, we sought to investigate DNA methylation and gene expression patterns in classical monocytes in a health disparity population.

METHODS

Classical monocytes (CD14+ + CD16-) were FACS-isolated from 34 self-reported African American women. Samples from 12 SSc patients and 12 healthy controls were hybridized on MethylationEPIC BeadChip array, while RNA-seq was performed on 16 SSc patients and 18 healthy controls. Analyses were computed to identify differentially methylated CpGs (DMCs), differentially expressed genes (DEGs), and CpGs associated with changes in gene expression (eQTM analysis).

RESULTS

We observed modest DNA methylation and gene expression differences between cases and controls. The genes harboring the top DMCs, the top DEGs, as well as the top eQTM loci were enriched for metabolic processes. Genes involved in immune processes and pathways showed a weak upregulation in the transcriptomic analysis. While many genes were newly identified, several other have been previously reported as differentially methylated or expressed in different blood cells from patients with SSc, supporting for their potential dysregulation in SSc.

CONCLUSIONS

While contrasting with results found in other blood cell types in largely European-descent groups, the results of this study support that variation in DNA methylation and gene expression exists among different cell types and individuals of different genetic, clinical, social, and environmental backgrounds. This finding supports the importance of including diverse, well-characterized patients to understand the different roles of DNA methylation and gene expression variability in the dysregulation of classical monocytes in diverse populations, which might help explaining the health disparities.

Promises and challenges in pharmacoepigenetics

Pharmacogenetics, the study of how interindividual genetic differences affect drug response, does not explain all observed heritable variance in drug response. Epigenetic mechanisms, such as DNA methylation, and histone acetylation may account for some of the unexplained variances. Epigenetic mechanisms modulate gene expression and can be suitable drug targets and can impact the action of nonepigenetic drugs. Pharmacoepigenetics is the field that studies the relationship between epigenetic variability and drug response. Much of this research focuses on compounds targeting epigenetic mechanisms, called epigenetic drugs, which are used to treat cancers, immune disorders, and other diseases. Several studies also suggest an epigenetic role in classical drug response; however, we know little about this area. The amount of information correlating epigenetic biomarkers to molecular datasets has recently expanded due to technological advances, and novel computational approaches have emerged to better identify and predict epigenetic interactions. We propose that the relationship between epigenetics and classical drug response may be examined using data already available by (1) finding regions of epigenetic variance, (2) pinpointing key epigenetic biomarkers within these regions, and (3) mapping these biomarkers to a drug-response phenotype. This approach expands on existing knowledge to generate putative pharmacoepigenetic relationships, which can be tested experimentally. Epigenetic modifications are involved in disease and drug response. Therefore, understanding how epigenetic drivers impact the response to classical drugs is important for improving drug design and administration to better treat disease.

Liquid Biopsies or Therapeutic Drug Monitoring for CYP Activity Profile Determination

During pharmacotherapy, knowledge about the actual drug and metabolite concentrations in plasma is often critical. Individual dose adjustments can be performed based on pre-emptive genotyping of certain absorption, distribution, metabolism, and excretion (ADME) genes but also using therapeutic drug monitoring (TDM). Analyses of liquid biopsies for tumor-derived components are well-established and have been found to be a good complement to biopsy examinations. Recently, liquid biopsy-based quantification of cell-free RNA (cfRNA) in plasma exosomes was proposed as a proxy measurement for the expression of different hepatic ADME genes and for the rate of drug metabolism, constituting an alternative to TDM. In this study, we validated these findings by examining the correlation between mRNA expression of eight different CYP genes in liver and the corresponding rate of enzyme-specific drug metabolism in 96 donor-matched liver samples. Analyses of CYP-dependent drug metabolism in liver microsomes in comparison to the level of mRNA expression for the different CYP genes revealed a mean Pearson correlation coefficient of 0.28. The highest correlations (0.33-0.34) were obtained for CYP2D6 and CYP3A4 and the weakest correlations were observed for CYP1A2 and CYP2B6 (0.18-0.21). In all cases, the correlations obtained were too weak to demonstrate a predictive relationship, likely due to different regulatory and post-translational events controlling the rate of enzyme activity. Our results reinforce the notion that, whilst liquid biopsy-based approaches might have utility for prediction of hepatic CYP protein expression, they are not currently an important substitute for TDM.

Methylation status of VTRNA2-1/nc886 is stable across populations, monozygotic twin pairs and in majority of tissues

Aims & methods: The aim of this study was to characterize the methylation level of a polymorphically imprinted gene, VTRNA2-1/nc886, in human populations and somatic tissues.48 datasets, consisting of more than 30 tissues and >30,000 individuals, were used. Results: nc886 methylation status is associated with twin status and ethnic background, but the variation between populations is limited. Monozygotic twin pairs present concordant methylation, whereas ∼30% of dizygotic twin pairs present discordant methylation in the nc886 locus. The methylation levels of nc886 are uniform across somatic tissues, except in cerebellum and skeletal muscle. Conclusion: The nc886 imprint may be established in the oocyte, and, after implantation, the methylation status is stable, excluding a few specific tissues.

Fetoplacental oxygen homeostasis in pregnancies with maternal diabetes mellitus and obesity

Despite improvements in clinical management, pregnancies complicated by pre-existing diabetes mellitus, gestational diabetes mellitus or obesity carry substantial risks for parent and offspring. Some of the endocrine and metabolic changes in parent and fetus in diabetes mellitus and obesity lead to fetal oxygen deficit, mostly due to insulin-induced accelerated fetal metabolism. The human fetus deals with reduced oxygenation through a wide range of adaptive responses that act at various levels in the placenta as well as the fetus. These responses ensure adequate oxygen delivery to the fetus, increase the oxygen transport capacity of fetal blood and redistribute oxygen-rich blood to vital organs such as the brain and heart. The liver has a central role in adapting to reduced oxygenation by increasing its oxygen extraction and stimulating erythropoietin synthesis to increase haematocrit. The type of adaptive response depends on the onset and duration of hypoxia and the severity of the metabolic disturbance. In pregnancies characterized by diabetes mellitus or obesity, these adaptive systems come under additional strain owing to the increased maternal supply of glucose and resultant fetal hyperinsulinaemia, both of which stimulate oxidative metabolism. In the rare situation that the adaptive responses are overwhelmed, stillbirth can ensue.

Spatiotemporal specificity of correlated DNA methylation and gene expression pairs across different human tissues and stages of brain development

DNA methylation (DNAm) that occurs on promoter regions is primarily considered to repress gene expression. Previous studies indicated that DNAm could also show positive correlations with gene expression. Both DNAm and gene expression profiles are known to be tissue- and development-specific. This study aims to investigate how DNAm and gene expression are coordinated across different human tissues and developmental stages, as well as the biological significance of such correlations. By analyzing 2,239 samples with both DNAm and gene expression data in the same human subjects obtained from six published datasets, we evaluated the correlations between gene and CpG pairs (GCPs) at cis-regions and compared significantly correlated GCPs (cGCPs) across different tissues and brains at different age groups. A total of 37,363 cGCPs was identified in the six datasets; approximately 38% of the cGCPs were positively correlated. The majority (>90%) of cGCPs was tissue- or development-specific. We also observed that the correlation direction can be opposite in different tissues and ages. Further analysis highlights the importance of cGCPs for their cellular functions and potential roles in complex traits and human diseases. For instance, the early developmental brain possessed a highly unique set of cGCPs that were associated with neurogenesis and psychiatric disorders. By assessing the epigenetic factors involved in cGCPs, we discovered novel regulatory mechanisms of positive cGCPs distinct from negative cGCPs, which were related to multiple factors, such as H3K27me3, CTCF, and JARD2. The catalogue of cGCPs compiled can be used to guide functional interpretation of genetic and epigenetic studies.

Aquaporin 9 induction in human iPSC-derived hepatocytes facilitates modeling of ornithine transcarbamylase deficiency

BACKGROUND AND AIMS

Patient-derived human-induced pluripotent stem cells (hiPSCs) differentiated into hepatocytes (hiPSC-Heps) have facilitated the study of rare genetic liver diseases. Here, we aimed to establish an in vitro liver disease model of the urea cycle disorder ornithine transcarbamylase deficiency (OTCD) using patient-derived hiPSC-Heps.

APPROACH AND RESULTS

Before modeling OTCD, we addressed the question of why hiPSC-Heps generally secrete less urea than adult primary human hepatocytes (PHHs). Because hiPSC-Heps are not completely differentiated and maintain some characteristics of fetal PHHs, we compared gene-expression levels in human fetal and adult liver tissue to identify genes responsible for reduced urea secretion in hiPSC-Heps. We found lack of aquaporin 9 (AQP9) expression in fetal liver tissue as well as in hiPSC-Heps, and showed that forced expression of AQP9 in hiPSC-Heps restores urea secretion and normalizes the response to ammonia challenge by increasing ureagenesis. Furthermore, we proved functional ureagenesis by challenging AQP9-expressing hiPSC-Heps with ammonium chloride labeled with the stable isotope [¹⁵ N] (¹⁵ NH₄ Cl) and by assessing enrichment of [¹⁵ N]-labeled urea. Finally, using hiPSC-Heps derived from patients with OTCD, we generated a liver disease model that recapitulates the hepatic manifestation of the human disease. Restoring OTC expression-together with AQP9-was effective in fully correcting OTC activity and normalizing ureagenesis as assessed by ¹⁵ NH₄ Cl stable-isotope challenge.

CONCLUSION

Our results identify a critical role for AQP9 in functional urea metabolism and establish the feasibility of in vitro modeling of OTCD with hiPSC-Heps. By facilitating studies of OTCD genotype/phenotype correlation and drug screens, our model has potential for improving the therapy of OTCD.

Epigenetics in drug disposition & drug therapy: symposium report of the 24th North American meeting of the International Society for the Study of Xenobiotics (ISSX)

The 24th North American International Society for the Study of Xenobiotics (ISSX) meeting, held virtually from September 13 to 17, 2021, embraced the theme of "Broadening Our Horizons." This reinforces a key mission of ISSX: striving to share innovative science related to drug discovery and development. Session speakers and the ISSX New Investigators Group, which supports the scientific and professional development of student and early career ISSX members, elected to highlight the scientific content presented during the captivating session titled, "Epigenetics in Drug Disposition & Drug Therapy." The impact genetic variation has on drug response is well established; however, this session underscored the importance of investigating the role of epigenetics in drug disposition and drug discovery. Session speakers, Drs. Ning, McClay, and Lazarus, detailed mechanisms by which epigenetic players including long non-coding RNA (lncRNAs), microRNA (miRNAs), DNA methylation, and histone acetylation can alter the expression of genes involved in pharmacokinetics, pharmacodynamics, and toxicity. Dr. Ning detailed current knowledge about miRNAs and lncRNAs and the mechanisms by which they can affect the expression of drug metabolizing enzymes (DMEs) and nuclear receptors. Dr. Lazarus discussed the potential role of miRNAs on UDP-glucuronosyltransferase (UGT) expression and activity. Dr. McClay provided evidence that aging alters methylation and acetylation of DMEs in the liver, affecting gene expression and activity. These topics, compiled by the symposium organizers, presenters, and the ISSX New Investigators Group, are herein discussed, along with exciting future perspectives for epigenetics in drug disposition and drug discovery research.

DNA methylation fingerprint of hepatocellular carcinoma from tissue and liquid biopsies

Hepatocellular carcinoma (HCC) is amongst the cancers with highest mortality rates and is the most common malignancy of the liver. Early detection is vital to provide the best treatment possible and liquid biopsies combined with analysis of circulating tumour DNA methylation show great promise as a non-invasive approach for early cancer diagnosis and monitoring with low false negative rates. To identify reliable diagnostic biomarkers of early HCC, we performed a systematic analysis of multiple hepatocellular studies and datasets comprising > 1500 genome-wide DNA methylation arrays, to define a methylation signature predictive of HCC in both tissue and cell-free DNA liquid biopsy samples. Our machine learning pipeline identified differentially methylated regions in HCC, some associated with transcriptional repression of genes related with cancer progression, that benchmarked positively against independent methylation signatures. Combining our signature of 38 DNA methylation regions, we derived a HCC detection score which confirmed the utility of our approach by identifying in an independent dataset 96% of HCC tissue samples with a precision of 98%, and most importantly successfully separated cfDNA of tumour samples from healthy controls. Notably, our risk score could identify cell-free DNA samples from patients with other tumours, including colorectal cancer. Taken together, we propose a comprehensive HCC DNA methylation fingerprint and an associated risk score for detection of HCC from tissue and liquid biopsies.

Innovations in Genomics and Big Data Analytics for Personalized Medicine and Health Care: A Review

Big data in health care is a fast-growing field and a new paradigm that is transforming case-based studies to large-scale, data-driven research. As big data is dependent on the advancement of new data standards, technology, and relevant research, the future development of big data applications holds foreseeable promise in the modern day health care revolution. Enormously large, rapidly growing collections of biomedical omics-data (genomics, proteomics, transcriptomics, metabolomics, glycomics, etc.) and clinical data create major challenges and opportunities for their analysis and interpretation and open new computational gateways to address these issues. The design of new robust algorithms that are most suitable to properly analyze this big data by taking into account individual variability in genes has enabled the creation of precision (personalized) medicine. We reviewed and highlighted the significance of big data analytics for personalized medicine and health care by focusing mostly on machine learning perspectives on personalized medicine, genomic data models with respect to personalized medicine, the application of data mining algorithms for personalized medicine as well as the challenges we are facing right now in big data analytics.

Detecting Blood Methylation Signatures in Response to Childhood Cancer Radiotherapy via Machine Learning Methods

Radiotherapy is a helpful treatment for cancer, but it can also potentially cause changes in many molecules, resulting in adverse effects. Among these changes, the occurrence of abnormal DNA methylation patterns has alarmed scientists. To explore the influence of region-specific radiotherapy on blood DNA methylation, we designed a computational workflow by using machine learning methods that can identify crucial methylation alterations related to treatment exposure. Irrelevant methylation features from the DNA methylation profiles of 2052 childhood cancer survivors were excluded via the Boruta method, and the remaining features were ranked using the minimum redundancy maximum relevance method to generate feature lists. These feature lists were then fed into the incremental feature selection method, which uses a combination of deep forest, k-nearest neighbor, random forest, and decision tree to find the most important methylation signatures and build the best classifiers and classification rules. Several methylation signatures and rules have been discovered and confirmed, allowing for a better understanding of methylation patterns in response to different treatment exposures.

Impact of DNA methylation on ADME gene expression, drug disposition and efficacy

Interindividual differences in drug response have always existed in clinical treatment. Genes involved in drug absorption, distribution, metabolism, and excretion (ADME) play an important role in the process of pharmacokinetics. The effects of genetic polymorphism and nuclear receptors on the expression of drug metabolism enzymes and transporters can only explain some individual differences in clinical treatment. Several key ADME genes have been demonstrated to be regulated by epigenetic mechanisms that can potentially affect interindividual variability in medical treatment. Emerging studies have focused on the importance of DNA methylation for ADME gene expression and for drug response. Among them, the most studied is anti-tumor drugs, and followed by anti-tuberculous and anti-platelet drugs. Therefore, we provide an epigenetics perspective on variability in drug response. The review summarizes the correlation between ADME gene expression and DNA methylation, including the exact methylation locations, and focuses on the corresponding drug disposition and effects to illuminate interindividual differences in clinical medication.

Identification of autosomal cis expression quantitative trait methylation (cis eQTMs) in children's blood

Background

The identification of expression quantitative trait methylation (eQTMs), defined as associations between DNA methylation levels and gene expression, might help the biological interpretation of epigenome-wide association studies (EWAS). We aimed to identify autosomal cis eQTMs in children's blood, using data from 832 children of the Human Early Life Exposome (HELIX) project.

Methods

Blood DNA methylation and gene expression were measured with the Illumina 450K and the Affymetrix HTA v2 arrays, respectively. The relationship between methylation levels and expression of nearby genes (1 Mb window centered at the transcription start site, TSS) was assessed by fitting 13.6 M linear regressions adjusting for sex, age, cohort, and blood cell composition.

Results

We identified 39,749 blood autosomal cis eQTMs, representing 21,966 unique CpGs (eCpGs, 5.7% of total CpGs) and 8,886 unique transcript clusters (eGenes, 15.3% of total transcript clusters, equivalent to genes). In 87.9% of these cis eQTMs, the eCpG was located at <250 kb from eGene's TSS; and 58.8% of all eQTMs showed an inverse relationship between the methylation and expression levels. Only around half of the autosomal cis-eQTMs eGenes could be captured through annotation of the eCpG to the closest gene. eCpGs had less measurement error and were enriched for active blood regulatory regions and for CpGs reported to be associated with environmental exposures or phenotypic traits. In 40.4% of the eQTMs, the CpG and the eGene were both associated with at least one genetic variant. The overlap of autosomal cis eQTMs in children's blood with those described in adults was small (13.8%), and age-shared cis eQTMs tended to be proximal to the TSS and enriched for genetic variants.

Conclusions

This catalogue of autosomal cis eQTMs in children's blood can help the biological interpretation of EWAS findings and is publicly available at https://helixomics.isglobal.org/ and at Dryad (doi:10.5061/dryad.fxpnvx0t0).

Funding

The study has received funding from the European Community's Seventh Framework Programme (FP7/2007-206) under grant agreement no 308333 (HELIX project); the H2020-EU.3.1.2. - Preventing Disease Programme under grant agreement no 874583 (ATHLETE project); from the European Union's Horizon 2020 research and innovation programme under grant agreement no 733206 (LIFECYCLE project), and from the European Joint Programming Initiative "A Healthy Diet for a Healthy Life" (JPI HDHL and Instituto de Salud Carlos III) under the grant agreement no AC18/00006 (NutriPROGRAM project). The genotyping was supported by the projects PI17/01225 and PI17/01935, funded by the Instituto de Salud Carlos III and co-funded by European Union (ERDF, "A way to make Europe") and the Centro Nacional de Genotipado-CEGEN (PRB2-ISCIII). BiB received core infrastructure funding from the Wellcome Trust (WT101597MA) and a joint grant from the UK Medical Research Council (MRC) and Economic and Social Science Research Council (ESRC) (MR/N024397/1). INMA data collections were supported by grants from the Instituto de Salud Carlos III, CIBERESP, and the Generalitat de Catalunya-CIRIT. KANC was funded by the grant of the Lithuanian Agency for Science Innovation and Technology (6-04-2014_31V-66). The Norwegian Mother, Father and Child Cohort Study is supported by the Norwegian Ministry of Health and Care Services and the Ministry of Education and Research. The Rhea project was financially supported by European projects (EU FP6-2003-Food-3-NewGeneris, EU FP6. STREP Hiwate, EU FP7 ENV.2007.1.2.2.2. Project No 211250 Escape, EU FP7-2008-ENV-1.2.1.4 Envirogenomarkers, EU FP7-HEALTH-2009- single stage CHICOS, EU FP7 ENV.2008.1.2.1.6. Proposal No 226285 ENRIECO, EU- FP7- HEALTH-2012 Proposal No 308333 HELIX), and the Greek Ministry of Health (Program of Prevention of obesity and neurodevelopmental disorders in preschool children, in Heraklion district, Crete, Greece: 2011-2014; "Rhea Plus": Primary Prevention Program of Environmental Risk Factors for Reproductive Health, and Child Health: 2012-15). We acknowledge support from the Spanish Ministry of Science and Innovation through the "Centro de Excelencia Severo Ochoa 2019-2023" Program (CEX2018-000806-S), and support from the Generalitat de Catalunya through the CERCA Program. MV-U and CR-A were supported by a FI fellowship from the Catalan Government (FI-DGR 2015 and #016FI_B 00272). MC received funding from Instituto Carlos III (Ministry of Economy and Competitiveness) (CD12/00563 and MS16/00128).

Epigenome-wide association study of alcohol consumption in N = 8161 individuals and relevance to alcohol use disorder pathophysiology: identification of the cystine/glutamate transporter SLC7A11 as a top target

Alcohol misuse is common in many societies worldwide and is associated with extensive morbidity and mortality, often leading to alcohol use disorders (AUD) and alcohol-related end-organ damage. The underlying mechanisms contributing to the development of AUD are largely unknown; however, growing evidence suggests that alcohol consumption is strongly associated with alterations in DNA methylation. Identification of alcohol-associated methylomic variation might provide novel insights into pathophysiology and novel treatment targets for AUD. Here we performed the largest single-cohort epigenome-wide association study (EWAS) of alcohol consumption to date (N = 8161) and cross-validated findings in AUD populations with relevant endophenotypes, as well as alcohol-related animal models. Results showed 2504 CpGs significantly associated with alcohol consumption (Bonferroni p value < 6.8 × 10-8) with the five leading probes located in SLC7A11 (p = 7.75 × 10-108), JDP2 (p = 1.44 × 10-56), GAS5 (p = 2.71 × 10-47), TRA2B (p = 3.54 × 10-42), and SLC43A1 (p = 1.18 × 10-40). Genes annotated to associated CpG sites are implicated in liver and brain function, the cellular response to alcohol and alcohol-associated diseases, including hypertension and Alzheimer's disease. Two-sample Mendelian randomization confirmed the causal relationship of consumption on AUD risk (inverse variance weighted (IVW) p = 5.37 × 10-09). A methylation-based predictor of alcohol consumption was able to discriminate AUD cases in two independent cohorts (p = 6.32 × 10-38 and p = 5.41 × 10-14). The top EWAS probe cg06690548, located in the cystine/glutamate transporter SLC7A11, was replicated in an independent cohort of AUD and control participants (N = 615) and showed strong hypomethylation in AUD (p < 10-17). Decreased CpG methylation at this probe was consistently associated with clinical measures including increased heavy drinking days (p < 10-4), increased liver function enzymes (GGT (p = 1.03 × 10-21), ALT (p = 1.29 × 10-6), and AST (p = 1.97 × 10-8)) in individuals with AUD. Postmortem brain analyses documented increased SLC7A11 expression in the frontal cortex of individuals with AUD and animal models showed marked increased expression in liver, suggesting a mechanism by which alcohol leads to hypomethylation-induced overexpression of SLC7A11. Taken together, our EWAS discovery sample and subsequent validation of the top probe in AUD suggest a strong role of abnormal glutamate signaling mediated by methylomic variation in SLC7A11. Our data are intriguing given the prominent role of glutamate signaling in brain and liver and might provide an important target for therapeutic intervention.

Effect of host genetics on the gut microbiome in 7,738 participants of the Dutch Microbiome Project

Host genetics are known to influence the gut microbiome, yet their role remains poorly understood. To robustly characterize these effects, we performed a genome-wide association study of 207 taxa and 205 pathways representing microbial composition and function in 7,738 participants of the Dutch Microbiome Project. Two robust, study-wide significant (P < 1.89 × 10-10) signals near the LCT and ABO genes were found to be associated with multiple microbial taxa and pathways and were replicated in two independent cohorts. The LCT locus associations seemed modulated by lactose intake, whereas those at ABO could be explained by participant secretor status determined by their FUT2 genotype. Twenty-two other loci showed suggestive evidence (P < 5 × 10-8) of association with microbial taxa and pathways. At a more lenient threshold, the number of loci we identified strongly correlated with trait heritability, suggesting that much larger sample sizes are needed to elucidate the remaining effects of host genetics on the gut microbiome.

Genetic variation influencing DNA methylation provides insights into molecular mechanisms regulating genomic function

We determined the relationships between DNA sequence variation and DNA methylation using blood samples from 3,799 Europeans and 3,195 South Asians. We identify 11,165,559 SNP-CpG associations (methylation quantitative trait loci (meQTL), P < 10-14), including 467,915 meQTL that operate in trans. The meQTL are enriched for functionally relevant characteristics, including shared chromatin state, High-throuhgput chromosome conformation interaction, and association with gene expression, metabolic variation and clinical traits. We use molecular interaction and colocalization analyses to identify multiple nuclear regulatory pathways linking meQTL loci to phenotypic variation, including UBASH3B (body mass index), NFKBIE (rheumatoid arthritis), MGA (blood pressure) and COMMD7 (white cell counts). For rs6511961 , chromatin immunoprecipitation followed by sequencing (ChIP-seq) validates zinc finger protein (ZNF)333 as the likely trans acting effector protein. Finally, we used interaction analyses to identify population- and lineage-specific meQTL, including rs174548 in FADS1, with the strongest effect in CD8+ T cells, thus linking fatty acid metabolism with immune dysregulation and asthma. Our study advances understanding of the potential pathways linking genetic variation to human phenotype.

Identifying causal models between genetically regulated methylation patterns and gene expression in healthy colon tissue

BACKGROUND

DNA methylation is involved in the regulation of gene expression and phenotypic variation, but the inter-relationship between genetic variation, DNA methylation and gene expression remains poorly understood. Here we combine the analysis of genetic variants related to methylation markers (methylation quantitative trait loci: mQTLs) and gene expression (expression quantitative trait loci: eQTLs) with methylation markers related to gene expression (expression quantitative trait methylation: eQTMs), to provide novel insights into the genetic/epigenetic architecture of colocalizing molecular markers.

RESULTS

Normal mucosa from 100 patients with colon cancer and 50 healthy donors included in the Colonomics project have been analyzed. Linear models have been used to find mQTLs and eQTMs within 1 Mb of the target gene. From 32,446 eQTLs previously detected, we found a total of 6850 SNPs, 114 CpGs and 52 genes interrelated, generating 13,987 significant combinations of co-occurring associations (meQTLs) after Bonferromi correction. Non-redundant meQTLs were 54, enriched in genes involved in metabolism of glucose and xenobiotics and immune system. SNPs in meQTLs were enriched in regulatory elements (enhancers and promoters) compared to random SNPs within 1 Mb of genes. Three colorectal cancer GWAS SNPs were related to methylation changes, and four SNPs were related to chemerin levels. Bayesian networks have been used to identify putative causal relationships among associated SNPs, CpG and gene expression triads. We identified that most of these combinations showed the canonical pathway of methylation markers causes gene expression variation (60.1%) or non-causal relationship between methylation and gene expression (33.9%); however, in up to 6% of these combinations, gene expression was causing variation in methylation markers.

CONCLUSIONS

In this study we provided a characterization of the regulation between genetic variants and inter-dependent methylation markers and gene expression in a set of 150 healthy colon tissue samples. This is an important finding for the understanding of molecular susceptibility on colon-related complex diseases.

A targeted multi-omics approach reveals paraoxonase-1 as a determinant of obesity-associated fatty liver disease

Background

The multifactorial nature of non-alcoholic fatty liver disease cannot be explained solely by genetic factors. Recent evidence revealed that DNA methylation changes take place at proximal promoters within susceptibility genes. This emphasizes the need for integrating multiple data types to provide a better understanding of the disease’s pathogenesis. One such candidate gene is paraoxonase-1 (PON1). Substantial interindividual differences in PON1 are apparent and could influence disease risk later in life. The aim of this study was therefore to determine the different regulatory aspects of PON1 variability and to examine them in relation to the predisposition to obesity-associated fatty liver disease.

Results

A targeted multi-omics approach was applied to investigate the interplay between PON1 genetic variants, promoter methylation, expression profile and enzymatic activity in an adult patient cohort with extensive metabolic and hepatic characterisation including liver biopsy. Alterations in PON1 status were shown to correlate with waist-to-hip ratio and relevant features of liver pathology. Particularly, the regulatory polymorphism rs705379:C > T was strongly associated with more severe liver disease. Multivariable data analysis furthermore indicated a significant association of combined genetic and epigenetic PON1 regulation. This identified relationship postulates a role for DNA methylation as a mediator between PON1 genetics and expression, which is believed to further influence liver disease progression via modifications in PON1 catalytic efficiency.

Conclusions

Our findings demonstrate that vertical data-integration of genetic and epigenetic regulatory mechanisms generated a more in-depth understanding of the molecular basis underlying the development of obesity-associated fatty liver disease. We gained novel insights into how NAFLD classification and outcome are orchestrated, which could not have been obtained by exclusively considering genetic variation.

Supplementary Information

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

Epigenetic modifications precede molecular alterations and drive human hepatocarcinogenesis

Development of primary liver cancer is a multi-stage process. Detailed understanding of sequential epigenetic alterations is largely missing. Here, we performed Infinium Human Methylation 450k BeadChips and RNA sequencing analyses for genome-wide methylome and transcriptome profiling of cirrhotic liver (n=7), low- (n=4) and high-grade (n=9) dysplastic lesions, early (n=5) and progressed (n=3) hepatocellular carcinomas (HCC) synchronously detected in eight HCC patients with chronic hepatitis B infection. Integrative analyses of epigenetically driven molecular changes were identified and validated in two independent cohorts comprising 887 HCC. Mitochondrial DNA sequencing was further employed for clonality analyses and indicates multi-clonal origin in the majority of investigated HCC. Alterations in DNA methylation progressively increased from CL to dysplastic lesions and reached a maximum in early HCC. Associated early alterations identified by IPA pathway analyses involved apoptosis, immune regulation and stemness pathways, while late changes centered on cell survival, proliferation and invasion. We further validated putative 23 epi-drivers with concomitant expression changes and associated with overall survival. Functionally, Striatin 4 (STRN4) was demonstrated to be epigenetically regulated and inhibition of STRN4 significantly suppressed tumorigenicity of HCC cell lines.Overall, application of integrative genomic analyses defines epigenetic driver alterations and provides promising targets for novel therapeutic approaches.

Hard wiring of normal tissue-specific chromosome-wide gene expression levels is an additional factor driving cancer type-specific aneuploidies

Background

Many carcinomas have recurrent chromosomal aneuploidies specific to the tissue of tumor origin. The reason for this specificity is not completely understood.

Methods

In this study, we looked at the frequency of chromosomal arm gains and losses in different cancer types from the The Cancer Genome Atlas (TCGA) and compared them to the mean gene expression of each chromosome arm in corresponding normal tissues of origin from the Genotype-Tissue Expression (GTEx) database, in addition to the distribution of tissue-specific oncogenes and tumor suppressors on different chromosome arms.

Results

This analysis revealed a complex picture of factors driving tumor karyotype evolution in which some recurrent chromosomal copy number reflect the chromosome arm-wide gene expression levels of the their normal tissue of tumor origin.

Conclusions

We conclude that the cancer type-specific distribution of chromosomal arm gains and losses is potentially “hardwiring” gene expression levels characteristic of the normal tissue of tumor origin, in addition to broadly modulating the expression of tissue-specific tumor driver genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-021-00905-y.

Human methylome variation across Infinium 450K data on the Gene Expression Omnibus

While DNA methylation (DNAm) is the most-studied epigenetic mark, few recent studies probe the breadth of publicly available DNAm array samples. We collectively analyzed 35 360 Illumina Infinium HumanMethylation450K DNAm array samples published on the Gene Expression Omnibus. We learned a controlled vocabulary of sample labels by applying regular expressions to metadata and used existing models to predict various sample properties including epigenetic age. We found approximately two-thirds of samples were from blood, one-quarter were from brain and one-third were from cancer patients. About 19% of samples failed at least one of Illumina's 17 prescribed quality assessments; signal distributions across samples suggest modifying manufacturer-recommended thresholds for failure would make these assessments more informative. We further analyzed DNAm variances in seven tissues (adipose, nasal, blood, brain, buccal, sperm and liver) and characterized specific probes distinguishing them. Finally, we compiled DNAm array data and metadata, including our learned and predicted sample labels, into database files accessible via the recountmethylation R/Bioconductor companion package. Its vignettes walk the user through some analyses contained in this paper.

Prenatal lead exposure and cord blood DNA methylation in the Korean Exposome Study

BACKGROUND

Prenatal lead exposure has been reported to affect infant growth and nervous system development, as well as to influence DNA methylation. We conducted an epigenome-wide association study to identify associations between prenatal lead exposure and cord blood DNA methylation in Korean infants.

METHODS

Cord blood samples were assayed with the Illumina HumanMethylationEPIC BeadChip kits, and maternal blood lead levels during early and late pregnancy, as well as cord blood lead level, were measured. The association between CpG methylation and lead level was analyzed using the limma package, with adjusting for infant sex, maternal pre-pregnancy body mass index, and estimated leukocyte composition.

RESULTS

Among 364 blood samples (182 males and 182 females), those for which maternal and cord blood lead concentrations during early and later pregnancy was known were used for analysis. Maternal lead concentration in blood during early pregnancy was significantly associated with the methylation status of specific positions. After data stratification by infant sex, we found that, in males, the level of maternal blood lead was associated with 18 CpG sites during early pregnancy, and with one CpG site near the NBAS gene, during late pregnancy. In female samples, there was no significant association between DNA methylation and lead concentrations.

CONCLUSIONS

Prenatal lead exposure was associated with altered, gender-specific patterns of DNA methylation in Korean infants.

DNA methylation of blood cells is associated with prevalent type 2 diabetes in a meta-analysis of four European cohorts

BACKGROUND

Type 2 diabetes (T2D) is a heterogeneous disease with well-known genetic and environmental risk factors contributing to its prevalence. Epigenetic mechanisms related to changes in DNA methylation (DNAm), may also contribute to T2D risk, but larger studies are required to discover novel markers, and to confirm existing ones.

RESULTS

We performed a large meta-analysis of individual epigenome-wide association studies (EWAS) of prevalent T2D conducted in four European studies using peripheral blood DNAm. Analysis of differentially methylated regions (DMR) was also undertaken, based on the meta-analysis results. We found three novel CpGs associated with prevalent T2D in Europeans at cg00144180 (HDAC4), cg16765088 (near SYNM) and cg24704287 (near MIR23A) and confirmed three CpGs previously identified (mapping to TXNIP, ABCG1 and CPT1A). We also identified 77 T2D associated DMRs, most of them hypomethylated in T2D cases versus controls. In adjusted regressions among diabetic-free participants in ALSPAC, we found that all six CpGs identified in the meta-EWAS were associated with white cell-types. We estimated that these six CpGs captured 11% of the variation in T2D, which was similar to the variation explained by the model including only the common risk factors of BMI, sex, age and smoking (R2 = 10.6%).

CONCLUSIONS

This study identifies novel loci associated with T2D in Europeans. We also demonstrate associations of the same loci with other traits. Future studies should investigate if our findings are generalizable in non-European populations, and potential roles of these epigenetic markers in T2D etiology or in determining long term consequences of T2D.

Differential DNA Methylation Landscape in Skin Fibroblasts from African Americans with Systemic Sclerosis

The etiology and reasons underlying the ethnic disparities in systemic sclerosis (SSc) remain unknown. African Americans are disproportionally affected by SSc and yet are underrepresented in research. The aim of this study was to comprehensively investigate the association of DNA methylation levels with SSc in dermal fibroblasts from patients of African ancestry. Reduced representation bisulfite sequencing (RRBS) was performed on primary dermal fibroblasts from 15 SSc patients and 15 controls of African ancestry, and over 3.8 million CpG sites were tested for differential methylation patterns between cases and controls. The dermal fibroblasts from African American patients exhibited widespread reduced DNA methylation. Differentially methylated CpG sites were most enriched in introns and intergenic regions while depleted in 5' UTR, promoters, and CpG islands. Seventeen genes and eleven promoters showed significant differential methylation, mostly in non-coding RNA genes and pseudogenes. Gene set enrichment analysis (GSEA) and gene ontology (GO) analyses revealed an enrichment of pathways related to interferon signaling and mesenchymal differentiation. The hypomethylation of DLX5 and TMEM140 was accompanied by these genes' overexpression in patients but underexpression for lncRNA MGC12916. These data show that differential methylation occurs in dermal fibroblasts from African American patients with SSc and identifies novel coding and non-coding genes.

DNA methylation perturbations may link altered development and aging in the lung

Fetal perturbations in DNA methylation during lung development may reveal insights into the enduring impacts on adult lung health and disease during aging that have not been explored altogether before. We studied the association between genome-wide DNA-methylation and post-conception age in fetal-lung (n=78, 42 exposed to in-utero-smoke (IUS)) tissue and chronological age in adult-lung tissue (n=160, 114 with Chronic Obstructive Pulmonary Disease) using multi-variate linear regression models with covariate adjustment and tested for effect modification by phenotypes. Overlapping age-associations were evaluated for functional and tissue-specific enrichment using the Genotype-Tissue-Expression (GTEx) project. We identified 244 age-associated differentially methylated positions and 878 regions overlapping between fetal and adult-lung tissues. Hyper-methylated CpGs (96%) were enriched in transcription factor activity (FDR adjusted P=2x10-33) and implicated in developmental processes including embryonic organ morphogenesis, neurogenesis and growth delay. Hypo-methylated CpGs (2%) were enriched in oxido-reductase activity and VEGFA-VEGFR2 Signaling. Twenty-one age-by-sex and eleven age-by-pack-years interactions were statistically significant (FDR<0.05) in adult-lung tissue. DNA methylation in transcription factors during development in fetal lung recapitulates in adult-lung tissue with aging. These findings reveal molecular mechanisms and pathways that may link disrupted development in early-life and age-associated lung diseases.

Epigenome-wide association study of maternal hemoglobin A1c in pregnancy and cord blood DNA methylation

Background: Gestational hyperglycemia is associated with adverse perinatal outcomes and long-term offspring metabolic programming, likely through dysregulation of DNA methylation (DNAm). Materials & methods: We tested associations between maternal HbA1c and cord blood DNAm among 412 mother-child pairs in the genetics of glucose regulation in gestation and growth (Gen3G) and implemented Mendelian randomization to infer causality. We sought replication in an independent sample from Healthy Start. Results: Higher second trimester HbA1c levels were associated with lower DNAm at cg21645848 (p = 3.9 × 10^-11) near URGCP. Mendelian randomization and replication analyses showed same direction of effect between HbA1c and DNAm at cg21645848, but did not reach statistical significance. Conclusion: We found that higher maternal glycemia reflected by HbA1c is associated with cord blood DNAm at URGCP, a gene related with inflammatory pathways.

DNA methylation and lipid metabolism: an EWAS of 226 metabolic measures

BACKGROUND

The discovery of robust and trans-ethnically replicated DNA methylation markers of metabolic phenotypes, has hinted at a potential role of epigenetic mechanisms in lipid metabolism. However, DNA methylation and the lipid compositions and lipid concentrations of lipoprotein sizes have been scarcely studied. Here, we present an epigenome-wide association study (EWAS) (N = 5414 total) of mostly lipid-related metabolic measures, including a fine profiling of lipoproteins. As lipoproteins are the main players in the different stages of lipid metabolism, examination of epigenetic markers of detailed lipoprotein features might improve the diagnosis, prognosis, and treatment of metabolic disturbances.

RESULTS

We conducted an EWAS of leukocyte DNA methylation and 226 metabolic measurements determined by nuclear magnetic resonance spectroscopy in the population-based KORA F4 study (N = 1662) and replicated the results in the LOLIPOP, NFBC1966, and YFS cohorts (N = 3752). Follow-up analyses in the discovery cohort included investigations into gene transcripts, metabolic-measure ratios for pathway analysis, and disease endpoints. We identified 161 associations (p value < 4.7 × 10-10), covering 16 CpG sites at 11 loci and 57 metabolic measures. Identified metabolic measures were primarily medium and small lipoproteins, and fatty acids. For apolipoprotein B-containing lipoproteins, the associations mainly involved triglyceride composition and concentrations of cholesterol esters, triglycerides, free cholesterol, and phospholipids. All associations for HDL lipoproteins involved triglyceride measures only. Associated metabolic measure ratios, proxies of enzymatic activity, highlight amino acid, glucose, and lipid pathways as being potentially epigenetically implicated. Five CpG sites in four genes were associated with differential expression of transcripts in blood or adipose tissue. CpG sites in ABCG1 and PHGDH showed associations with metabolic measures, gene transcription, and metabolic measure ratios and were additionally linked to obesity or previous myocardial infarction, extending previously reported observations.

CONCLUSION

Our study provides evidence of a link between DNA methylation and the lipid compositions and lipid concentrations of different lipoprotein size subclasses, thus offering in-depth insights into well-known associations of DNA methylation with total serum lipids. The results support detailed profiling of lipid metabolism to improve the molecular understanding of dyslipidemia and related disease mechanisms.

Incorporation of DNA methylation into eQTL mapping in African Americans

Epigenetics is a reversible molecular mechanism that plays a critical role in many developmental, adaptive, and disease processes. DNA methylation has been shown to regulate gene expression and the advent of high throughput technologies has made genome-wide DNA methylation analysis possible. We investigated the effect of DNA methylation on eQTL mapping (methylation-adjusted eQTLs), by incorporating DNA methylation as a SNP-based covariate in eQTL mapping in African American derived hepatocytes. We found that the addition of DNA methylation uncovered new eQTLs and eGenes. Previously discovered eQTLs were significantly altered by the addition of DNA methylation data suggesting that methylation may modulate the association of SNPs to gene expression. We found that methylation-adjusted eQTLs that were less significant compared to PC-adjusted eQTLs were enriched in lipoprotein measurements (FDR=0.0040), immune system disorders (FDR = 0.0042), and liver enzyme measurements (FDR=0.047), suggesting that DNA methylation modulates the genetic regulation of these phenotypes. Our methylation-adjusted eQTL analysis also uncovered novel SNP-gene pairs. For example, we found that the SNP, rs1332018, was associated to GSTM3. GSTM3 expression has been linked to Hepatitis B which African Americans suffer from disproportionately. Our methylation-adjusted method adds new understanding to the genetic basis of complex diseases that disproportionally affect African Americans.

Hepatoblastomas exhibit marked NNMT downregulation driven by promoter DNA hypermethylation

Hepatoblastomas exhibit the lowest mutational burden among pediatric tumors. We previously showed that epigenetic disruption is crucial for hepatoblastoma carcinogenesis. Our data revealed hypermethylation of nicotinamide N-methyltransferase, a highly expressed gene in adipocytes and hepatocytes. The expression pattern and the role of nicotinamide N-methyltransferase in pediatric liver tumors have not yet been explored, and this study aimed to evaluate the effect of nicotinamide N-methyltransferase hypermethylation in hepatoblastomas. We evaluated 45 hepatoblastomas and 26 non-tumoral liver samples. We examined in hepatoblastomas if the observed nicotinamide N-methyltransferase promoter hypermethylation could lead to dysregulation of expression by measuring mRNA and protein levels by real-time quantitative polymerase chain reaction, immunohistochemistry, and Western blot assays. The potential impact of nicotinamide N-methyltransferase changes was evaluated on the metabolic profile by high-resolution magic angle spinning nuclear magnetic resonance spectroscopy. Significant nicotinamide N-methyltransferase downregulation was revealed in hepatoblastomas, with two orders of magnitude lower nicotinamide N-methyltransferase expression in tumor samples and hepatoblastoma cell lines than in hepatocellular carcinoma cell lines. A specific TSS1500 CpG site (cg02094283) of nicotinamide N-methyltransferase was hypermethylated in tumors, with an inverse correlation between its methylation level and nicotinamide N-methyltransferase expression. A marked global reduction of the nicotinamide N-methyltransferase protein was validated in tumors, with strong correlation between gene and protein expression. Of note, higher nicotinamide N-methyltransferase expression was statistically associated with late hepatoblastoma diagnosis, a known clinical variable of worse prognosis. In addition, untargeted metabolomics analysis detected aberrant lipid metabolism in hepatoblastomas. Data presented here showed the first evidence that nicotinamide N-methyltransferase reduction occurs in hepatoblastomas, providing further support that the nicotinamide N-methyltransferase downregulation is a wide phenomenon in liver cancer. Furthermore, this study unraveled the role of DNA methylation in the regulation of nicotinamide N-methyltransferase expression in hepatoblastomas, in addition to evaluate the potential effect of nicotinamide N-methyltransferase reduction in the metabolism of these tumors. These preliminary findings also suggested that nicotinamide N-methyltransferase level may be a potential prognostic biomarker for hepatoblastoma.

Differences in cytochrome p450 enzyme expression and activity in fetal and adult tissues

INTRODUCTION

Human cytochrome p450 (CYP) enzyme expression and activity is lower in the fetus as compared to the adult; however, limited quantitative data exists regarding the specific differences in magnitude or the degree of inducibility due to environmental factors.

METHODS

We utilized a combination of in silico- and molecular-based approaches to profile and compare CYP expression/activity in human adult liver and fetal tissues. Using public datasets, we evaluated human CYP expression between: 1) placenta vs. adult livers; 2) fetal vs. adult livers; or 3) five compartments of the human placenta. We generated new experimental data, characterizing expression levels of nine CYPs in placenta/fetal liver vs. adult liver. In a subset of samples, we evaluated CYP3A4 activity. Finally, we summarized evidence of human fetal CYP expression/activity and environmental exposures during pregnancy.

RESULTS

In silico, CYPs were predominately expressed at higher levels in the adult liver vs. fetal tissues, with a few noted exceptions. Sixty percent of CYP enzymes were expressed at nominal levels in the placenta. In wet-lab analyses, we observed significant CYP-specific differences in expression/activity between adult and fetal tissues; CYP2E1 and -3A4 were expressed significantly lower in fetal vs. adult livers, while CYP2J2 levels were similar.

DISCUSSION

We provide a qualitative review of the expression of the CYP enzyme family in critical sites of xenobiotic distribution during human pregnancy and novel quantitative data regarding fetal CYP expression and activity during mid-gestation. Data outputs may be a resource for modeling predictions of chemical distribution and sensitivity.

Epigenetics of paraoxonases

PURPOSE OF REVIEW

Studies have shown the three-member paraoxonase (PON) multigene family to be involved in the development of a large variety of diseases with an inflammatory component. Environmental factors such as lifestyle-related factors differ widely between populations and it is important to consider that their impacts may be exerted through the epigenetic mechanisms, which connect genes, the environment and disease development and are a potential therapeutic avenue.

RECENT FINDINGS

In the review period, very little was published on epigenetics of PON2 or PON3, mostly on their diagnostic value in cancer by measuring methylation levels of these genes. However, the picture is more promising with PON1. Here, several studies have linked the epigenetic regulation of PON1 to various metabolic processes and particularly to the development of several diseases, including stroke, heart disease, aortic valve stenosis and chronic obstructive pulmonary disease.

SUMMARY

Studies into the epigenetic regulation of the PON family are in their infancy. However, recent studies linking epigenetic regulation of PON1 to disease development will encourage further research and open up the possibility for new potential therapeutic interventions.

Integrated multiomics analysis of hepatoblastoma unravels its heterogeneity and provides novel druggable targets

Although hepatoblastoma is the most common pediatric liver cancer, its genetic heterogeneity and therapeutic targets are not well elucidated. Therefore, we conducted a multiomics analysis, including mutatome, DNA methylome, and transcriptome analyses, of 59 hepatoblastoma samples. Based on DNA methylation patterns, hepatoblastoma was classified into three clusters exhibiting remarkable correlation with clinical, histological, and genetic features. Cluster F was largely composed of cases with fetal histology and good outcomes, whereas clusters E1 and E2 corresponded primarily to embryonal/combined histology and poor outcomes. E1 and E2, albeit distinguishable by different patient age distributions, were genetically characterized by hypermethylation of the HNF4A/CEBPA-binding regions, fetal liver-like expression patterns, upregulation of the cell cycle pathway, and overexpression of NQO1 and ODC1. Inhibition of NQO1 and ODC1 in hepatoblastoma cells induced chemosensitization and growth suppression, respectively. Our results provide a comprehensive description of the molecular basis of hepatoblastoma and rational therapeutic strategies for high-risk cases.

Silencing hepatic MCJ attenuates non-alcoholic fatty liver disease (NAFLD) by increasing mitochondrial fatty acid oxidation

Nonalcoholic fatty liver disease (NAFLD) is considered the next major health epidemic with an estimated 25% worldwide prevalence. No drugs have yet been approved and NAFLD remains a major unmet need. Here, we identify MCJ (Methylation-Controlled J protein) as a target for non-alcoholic steatohepatitis (NASH), an advanced phase of NAFLD. MCJ is an endogenous negative regulator of the respiratory chain Complex I that acts to restrain mitochondrial respiration. We show that therapeutic targeting of MCJ in the liver with nanoparticle- and GalNAc-formulated siRNA efficiently reduces liver lipid accumulation and fibrosis in multiple NASH mouse models. Decreasing MCJ expression enhances the capacity of hepatocytes to mediate β-oxidation of fatty acids and minimizes lipid accumulation, which results in reduced hepatocyte damage and fibrosis. Moreover, MCJ levels in the liver of NAFLD patients are elevated relative to healthy subjects. Thus, inhibition of MCJ emerges as an alternative approach to treat NAFLD.

Non-alcoholic fatty liver (NAFLD) disease causes degeneration of the liver, affects about 25% of people globally, and has no approved treatment. Here, the authors show that the therapeutic siRNA-driven silencing of MCJ in the liver is an effective and safe treatment for NAFLD in multiple mouse models.

Identification, Heritability, and Relation With Gene Expression of Novel DNA Methylation Loci for Blood Pressure

We conducted an epigenome-wide association study meta-analysis on blood pressure (BP) in 4820 individuals of European and African ancestry aged 14 to 69. Genome-wide DNA methylation data from peripheral leukocytes were obtained using the Infinium Human Methylation 450k BeadChip. The epigenome-wide association study meta-analysis identified 39 BP-related CpG sites with P<1×10-5. In silico replication in the CHARGE consortium of 17 010 individuals validated 16 of these CpG sites. Out of the 16 CpG sites, 13 showed novel association with BP. Conversely, out of the 126 CpG sites identified as being associated (P<1×10-7) with BP in the CHARGE consortium, 21 were replicated in the current study. Methylation levels of all the 34 CpG sites that were cross-validated by the current study and the CHARGE consortium were heritable and 6 showed association with gene expression. Furthermore, 9 CpG sites also showed association with BP with P<0.05 and consistent direction of the effect in the meta-analysis of the Finnish Twin Cohort (199 twin pairs and 4 singletons; 61% monozygous) and the Netherlands Twin Register (266 twin pairs and 62 singletons; 84% monozygous). Bivariate quantitative genetic modeling of the twin data showed that a majority of the phenotypic correlations between methylation levels of these CpG sites and BP could be explained by shared unique environmental rather than genetic factors, with 100% of the correlations of systolic BP with cg19693031 (TXNIP) and cg00716257 (JDP2) determined by environmental effects acting on both systolic BP and methylation levels.

Association of common variation in ADD3 and GPC1 with biliary atresia susceptibility

Biliary atresia (BA) is an idiopathic neonatal cholestatic disease. Recent genome-wide association study (GWAS) revealed that common variation of ADD3, GPC1, ARF6, and EFEMP1 gene was associated with BA susceptibility. We aimed to evaluate the association of these genes with BA in Chinese population. Twenty single nucleotide polymorphisms (SNPs) in these four genes were genotyped in 340 BA patients and 1,665 controls. Three SNPs in ADD3 were significantly associated with BA, and rs17095355 was the top SNP (P_Allele = 3.23×10^-6). Meta-analysis of published data and current data indicated that rs17095355 was associated with BA susceptibility in Asians and Caucasians. Three associated SNPs were expression quantitative trait loci (eQTL) for ADD3. Two GPC1 SNPs in high linkage disequilibrium (LD) showed nominal association with BA susceptibility (P_Allele = 0.03 for rs6707262 and P_Allele = 0.04 for rs6750380), and were eQTL of GPC1. Haplotype harboring these two SNPs almost reached the study-wide significance (P = 0.0035). No association for ARF6 and EFEMP1 was found with BA risk in the current population. Our study validated associations of ADD3 and GPC1 SNPs with BA risk in Chinese population and provided evidence of epistatic contributions of genetic factors to BA susceptibility.

DNA methylation and histone acetylation changes to cytochrome P450 2E1 regulation in normal aging and impact on rates of drug metabolism in the liver

Aging is associated with reduced liver function that may increase the risk for adverse drug reactions in older adults. We hypothesized that age-related changes to epigenetic regulation of genes involved in drug metabolism may contribute to this effect. We reviewed published epigenome-wide studies of human blood and identified the cytochrome P450 2E1 (CYP2E1) gene as a top locus exhibiting epigenetic changes with age. To investigate potential functional changes with age in the liver, the primary organ of drug metabolism, we obtained liver tissue from mice aged 4-32 months from the National Institute on Aging. We assayed global DNA methylation (5-methylcytosine, 5mC), hydroxymethylation (5-hydroxymethylcytosine, 5hmC), and locus-specific 5mC and histone acetylation changes around mouse Cyp2e1. The mouse livers exhibit significant global decreases in 5mC and 5hmC with age. Furthermore, 5mC significantly increased with age at two regulatory regions of Cyp2e1 in tandem with decreases in its gene and protein expressions. H3K9ac levels also changed with age at both regulatory regions of Cyp2e1 investigated, while H3K27ac did not. To test if these epigenetic changes are associated with varying rates of drug metabolism, we assayed clearance of the CYP2E1-specific probe drug chlorzoxazone in microsome extracts from the same livers. CYP2E1 intrinsic clearance is associated with DNA methylation and H3K9ac levels at the Cyp2e1 locus but not with chronological age. This suggests that age-related epigenetic changes may influence rates of hepatic drug metabolism. In the future, epigenetic biomarkers could prove useful to guide dosing regimens in older adults.

Human Pluripotent Stem Cell-Derived Hepatocytes Show Higher Transcriptional Correlation with Adult Liver Tissue than with Fetal Liver Tissue

Human pluripotent stem cell-derived hepatocytes (hPSC-HEP) display many properties of mature hepatocytes, including expression of important genes of the drug metabolizing machinery, glycogen storage, and production of multiple serum proteins. To this date, hPSC-HEP do not, however, fully recapitulate the complete functionality of in vivo mature hepatocytes. In this study, we applied versatile bioinformatic algorithms, including functional annotation and pathway enrichment analyses, transcription factor binding-site enrichment, and similarity and correlation analyses, to datasets collected from different stages during hPSC-HEP differentiation and compared these to developmental stages and tissues from fetal and adult human liver. Our results demonstrate a high level of similarity between the in vitro differentiation of hPSC-HEP and in vivo hepatogenesis. Importantly, the transcriptional correlation of hPSC-HEP with adult liver (AL) tissues was higher than with fetal liver (FL) tissues (0.83 and 0.70, respectively). Functional data revealed mature features of hPSC-HEP including cytochrome P450 enzymes activities and albumin secretion. Moreover, hPSC-HEP showed expression of many genes involved in drug absorption, distribution, metabolism, and excretion. Despite the high similarities observed, we identified differences of specific pathways and regulatory players by analyzing the gene expression between hPSC-HEP and AL. These findings will aid future intervention and improvement of in vitro hepatocyte differentiation protocol in order to generate hepatocytes displaying the complete functionality of mature hepatocytes. Finally, on the transcriptional level, our results show stronger correlation and higher similarity of hPSC-HEP to AL than to FL. In addition, potential targets for further functional improvement of hPSC-HEP were also identified.

Whole-genome re-sequencing and transcriptome reveal cadmium tolerance related genes and pathways in Chlamydomonas reinhardtii

Cadmium (Cd), a common environmental toxic contaminant, is easily accumulated in living organisms, leading to numerous harmful effects. Chlamydomonas reinhardtii, a unicellular eukaryotic green algae strain, is a very suitable candidate for bioremediation of Cd-contaminated water. However, for the poor resistance to Cd, application of C. reinhardtii was restricted and genes mediating Cd tolerance in C. reinhardtii remain unclear. In this paper, adaptive laboratory evolution was performed with algae constant exposure to Cd over 420-day at environmentally relevant concentrations to select C. reinhardtii strains with high tolerance to Cd. Physiological indicators, such as cell proliferation, photosynthetic pigment contents and photosynthetic activity of photosystem were detected to evaluate the Cd tolerance of selected algae strain ALE0.5. Then, whole-genome re-sequencing and transcriptome were applied to identify the genes related to Cd tolerance. Genes involved in photosynthesis (PSBP1), glutathione metabolism (CHLREDRAFT_167073, GPX5) and calcium transport (CHLREDRAFT_189266, CHLREDRAFT_191203, CHLREDRAFT_187187, CSE1) were related to Cd tolerance in C. reinhardtii. This study provides a basis for obtaining transgenic C. reinhardtii strains with high Cd tolerance used for bioremediation of Cd pollution in the future.

Linkage-linkage disequilibrium dissection of the epigenetic quantitative trait loci (epiQTLs) underlying growth and wood properties in Populus

Increasing evidence indicates that DNA methylation is heritable and serves as an essential marker contributing to phenotypic variation. Linkage-linkage disequilibrium mapping was used to decipher the epigenetic architecture underlying nine growth and wood property traits in a linkage population (550 F₁ progeny) and a natural population (435 unrelated individuals) of Populus using methylation-sensitive amplification polymorphism (MSAP)-based analysis. The interactions between genetic and epigenetic variants in the causative genes was further unveiled using expression quantitative trait methylation (eQTM) and nucleotide (eQTN) mapping strategies. A total of 163 epigenetic quantitative trait loci (epiQTLs; LOD ≥ 3.0), explaining 1.7-44.5% of phenotypic variations, were mapped to a high-resolution epigenetic map with 19 linkage groups, which was supported by the significant MSAP associations (P < 0.001) in the two populations. There were 23 causal genes involved in growth regulation and wood formation, whose markers were located in epiQTLs and associated with the same traits in both populations. Further eQTN and eQTM mapping showed that causal genetic and epigenetic variants within the 23 candidate genes may interact more in trans in gene expression and phenotype. The present study provides strategies for investigating epigenetic architecture and the interaction between genetic and epigenetic variants modulating complex traits in forest trees.

The effects of bariatric surgery on clinical profile, DNA methylation, and ageing in severely obese patients

BACKGROUND

Severe obesity is a growing, worldwide burden and conventional therapies including radical change of diet and/or increased physical activity have limited results. Bariatric surgery has been proposed as an alternative therapy showing promising results. It leads to substantial weight loss and improvement of comorbidities such as type 2 diabetes. Increased adiposity is associated with changes in epigenetic profile, including DNA methylation. We investigated the effect of bariatric surgery on clinical profile, DNA methylation, and biological age estimated using Horvath's epigenetic clock.

RESULTS

To determine the impact of bariatric surgery and subsequent weight loss on clinical traits, a cohort of 40 severely obese individuals (BMI = 30-73 kg/m²) was examined at the time of surgery and at three follow-up visits, i.e., 3, 6, and 12 months after surgery. The majority of the individuals were women (65%) and the mean age at surgery was 45.1 ± 8.1 years. We observed a significant decrease over time in BMI, fasting glucose, HbA1c, HOMA-IR, insulin, total cholesterol, triglycerides, LDL and free fatty acids levels, and a significant small increase in HDL levels (all p values < 0.05). Epigenome-wide association analysis revealed 4857 differentially methylated CpG sites 12 months after surgery (at Bonferroni-corrected p value < 1.09 × 10^-7). Including BMI change in the model decreased the number of significantly differentially methylated CpG sites by 51%. Gene set enrichment analysis identified overrepresentation of multiple processes including regulation of transcription, RNA metabolic, and biosynthetic processes in the cell. Bariatric surgery in severely obese patients resulted in a decrease in both biological age and epigenetic age acceleration (EAA) (mean = - 0.92, p value = 0.039).

CONCLUSIONS

Our study shows that bariatric surgery leads to substantial BMI decrease and improvement of clinical outcomes observed 12 months after surgery. These changes explained part of the association between bariatric surgery and DNA methylation. We also observed a small, but significant improvement of biological age. These epigenetic changes may be modifiable by environmental lifestyle factors and could be used as potential biomarkers for obesity and in the future for obesity related comorbidities.

Leveraging cell-specific differentially methylated regions to identify leukocyte infiltration in adipose tissue

Obesity is understood to be an inflammatory condition characterized in part by changes in resident immune cell populations in adipose tissue. However, much of this knowledge has been obtained through experimental animal models. Epigenetic mechanisms, such as DNA methylation may be useful tools for characterizing the changes in immune cell populations in human subjects. In this study, we introduce a simple and intuitive method for assessing cellular infiltration by blood into other heterogeneous, admixed tissues such as adipose tissue, and apply this approach in a large human cohort study. Associations between higher leukocyte infiltration, measured by evaluating a distance measure between the methylation signatures of leukocytes and adipose tissue, and increasing body mass index (BMI) or android fat mass (AFM) were identified and validated in independent replication samples for CD4 (pBMI  = 0.009, pAFM  = 0.020), monocytes (pBMI  = 0.001, pAFM  = 4.3 × 10-4 ), and dendritic cells (pBMI  = 0.571, pAFM  = 0.012). Patterns of depletion with increasing adiposity were observed for plasma B (pBMI  = 0.430, pAFM  = 0.004) and immature B (pBMI  = 0.022, pAFM  = 0.042) cells. CD4, dendritic, monocytes, immature B, and plasma B cells may be important agents in the inflammatory process. Finally, the method used to assess leukocyte infiltration in this study is straightforwardly extended to other cell types and tissues in which infiltration might be of interest.

Molecular Aging of Human Liver: An Epigenetic/Transcriptomic Signature

The feasibility of liver transplantation from old healthy donors suggests that this organ is able to preserve its functionality during aging. To explore the biological basis of this phenomenon, we characterized the epigenetic profile of liver biopsies collected from 45 healthy liver donors ranging from 13 to 90 years old using the Infinium HumanMethylation450 BeadChip. The analysis indicates that a large remodeling in DNA methylation patterns occurs, with 8,823 age-associated differentially methylated CpG probes. Notably, these age-associated changes tended to level off after the age of 60, as confirmed by Horvath's clock. Using stringent selection criteria, we further identified a DNA methylation signature of aging liver including 75 genomic regions. We demonstrated that this signature is specific for liver compared to other tissues and that it is able to detect biological age-acceleration effects associated with obesity. Finally, we combined DNA methylation measurements with available expression data. Although the intersection between the two omic characterizations was low, both approaches suggested a previously unappreciated role of epithelial-mesenchymal transition and Wnt-signaling pathways in the aging of human liver.