Metaboepigenetic regulation of gene expression in obesity and insulin resistance.

INTRODUCTION: Variation in DNA methylation (DNAm) in adipose tissue is associated with the pathogenesis of obesity and insulin resistance. The activity of enzymes involved in altering DNAm levels is dependent on several metabolite cofactors. OBJECTIVES: To understand the role of metabolites as mechanistic regulators of epigenetic marks, we tested the association between selected plasma metabolites and DNAm levels in the adipose tissue of African Americans. METHODS: In the AAGMEx cohort (N = 256), plasma levels of metabolites were measured by untargeted liquid chromatography-mass spectrometry; adipose tissue DNAm and transcript levels were measured by reduced representation bisulfite sequencing, and expression microarray, respectively. RESULTS: Among the 21 one-carbon metabolism pathway metabolites evaluated, six were associated with gluco-metabolic traits (PFDR < 0.05, for BMI, SI, or Matsuda index) in AAGMEx. Methylation levels of 196, 116, and 180 CpG-sites were associated (P < 0.0001) with S-adenosylhomocysteine (SAH), cystine, and hypotaurine, respectively. Cis-expression quantitative trait methylation (cis eQTM) analyses suggested the role of metabolite-level-associated CpG sites in regulating the expression of adipose tissue transcripts, including genes in G-protein coupled receptor signaling pathway. Plasma SAH level-associated CpG sites chr19:3403712 and chr19:3403735 were also associated with the expression of G-protein subunit alpha 15 (GNA15) in adipose. The expression of GNA15 was significantly correlated with BMI (ß = 1.87, P = 1.9 × 10-16) and SI (ß = -1.61, P = 2.49 × 10-5). CONCLUSION: Our study suggests that a subset of metabolites modulates the methylation levels of CpG sites in specific loci and, in turn, regulates the expression of transcripts involved in obesity and insulin resistance.

Torch-eCpG: a fast and scalable eQTM mapper for thousands of molecular phenotypes with graphical processing units.

BACKGROUND: Gene expression may be regulated by the DNA methylation of regulatory elements in cis, distal, and trans regions. One method to evaluate the relationship between DNA methylation and gene expression is the mapping of expression quantitative trait methylation (eQTM) loci (also called expression associated CpG loci, eCpG). However, no open-source tools are available to provide eQTM mapping. In addition, eQTM mapping can involve a large number of comparisons which may prevent the analyses due to limitations of computational resources. Here, we describe Torch-eCpG, an open-source tool to perform eQTM mapping that includes an optimized implementation that can use the graphical processing unit (GPU) to reduce runtime. RESULTS: We demonstrate the analyses using the tool are reproducible, up to 18 × faster using the GPU, and scale linearly with increasing methylation loci. CONCLUSIONS: Torch-eCpG is a fast, reliable, and scalable tool to perform eQTM mapping. Source code for Torch-eCpG is available at https://github.com/kordk/torch-ecpg .

Epigenome-wide DNA methylation association study of circulating IgE levels identifies novel targets for asthma.

BACKGROUND: Identifying novel epigenetic signatures associated with serum immunoglobulin E (IgE) may improve our understanding of molecular mechanisms underlying asthma and IgE-mediated diseases. METHODS: We performed an epigenome-wide association study using whole blood from Framingham Heart Study (FHS; n = 3,471, 46% females) participants and validated results using the Childhood Asthma Management Program (CAMP; n = 674, 39% females) and the Genetic Epidemiology of Asthma in Costa Rica Study (CRA; n = 787, 41% females). Using the closest gene to each IgE-associated CpG, we highlighted biologically plausible pathways underlying IgE regulation and analyzed the transcription patterns linked to IgE-associated CpGs (expression quantitative trait methylation loci; eQTMs). Using prior UK Biobank summary data from genome-wide association studies of asthma and allergy, we performed Mendelian randomization (MR) for causal inference testing using the IgE-associated CpGs from FHS with methylation quantitative trait loci (mQTLs) as instrumental variables. FINDINGS: We identified 490 statistically significant differentially methylated CpGs associated with IgE in FHS, of which 193 (39.3%) replicated in CAMP and CRA (FDR < 0.05). Gene ontology analysis revealed enrichment in pathways related to transcription factor binding, asthma, and other immunological processes. eQTM analysis identified 124 cis-eQTMs for 106 expressed genes (FDR < 0.05). MR in combination with drug-target analysis revealed CTSB and USP20 as putatively causal regulators of IgE levels (Bonferroni adjusted P < 7.94E-04) that can be explored as potential therapeutic targets. INTERPRETATION: By integrating eQTM and MR analyses in general and clinical asthma populations, our findings provide a deeper understanding of the multidimensional inter-relations of DNA methylation, gene expression, and IgE levels. FUNDING: US NIH/NHLBI grants: P01HL132825, K99HL159234. N01-HC-25195 and HHSN268201500001I.

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.

Cis- and trans-eQTM analysis reveals novel epigenetic and transcriptomic immune markers of atopic asthma in airway epithelium.

BACKGROUND: Expression quantitative trait methylation (eQTM) analyses uncover associations between DNA methylation markers and gene expression. Most eQTM analyses of complex diseases have focused on cis-eQTM pairs (within 1 megabase). OBJECTIVES: This study sought to identify cis- and trans-methylation markers associated with gene expression in airway epithelium from youth with and without atopic asthma. METHODS: In this study, the investigators conducted both cis- and trans-eQTM analyses in nasal (airway) epithelial samples from 158 Puerto Rican youth with atopic asthma and 100 control subjects without atopy or asthma. The investigators then attempted to replicate their findings in nasal epithelial samples from 2 studies of children, while also examining whether their results in nasal epithelium overlap with those from an eQTM analysis in white blood cells from the Puerto Rican subjects. RESULTS: This study identified 9,108 cis-eQTM pairs and 2,131,500 trans-eQTM pairs. Trans-associations were significantly enriched for transcription factor and microRNA target genes. Furthermore, significant cytosine-phosphate-guanine sites (CpGs) were differentially methylated in atopic asthma and significant genes were enriched for genes differentially expressed in atopic asthma. In this study, 50.7% to 62.6% of cis- and trans-eQTM pairs identified in Puerto Rican youth were replicated in 2 smaller cohorts at false discovery rate-adjusted P < .1. Replicated genes in the trans-eQTM analysis included biologically plausible asthma-susceptibility genes (eg, HDC, NLRP3, ITGAE, CDH26, and CST1) and are enriched in immune pathways. CONCLUSIONS: Studying both cis- and trans-epigenetic regulation of airway epithelial gene expression can identify potential causal and regulatory pathways or networks for childhood asthma. Trans-eQTM CpGs may regulate gene expression in airway epithelium through effects on transcription factor and microRNA target genes.

Expression Quantitative Trait Methylation Analysis Identifies Whole Blood Molecular Footprint in Fetal Alcohol Spectrum Disorder (FASD).

Fetal alcohol spectrum disorder (FASD) encompasses neurodevelopmental disabilities and physical birth defects associated with prenatal alcohol exposure. Previously, we attempted to identify epigenetic biomarkers for FASD by investigating the genome-wide DNA methylation (DNAm) profiles of individuals with FASD compared to healthy controls. In this study, we generated additional gene expression profiles in a subset of our previous FASD cohort, encompassing the most severely affected individuals, to examine the functional integrative effects of altered DNAm status on gene expression. We identified six differentially methylated regions (annotated to the SEC61G, REEP3, ZNF577, HNRNPF, MSC, and SDHAF1 genes) associated with changes in gene expression (p-value < 0.05). To the best of our knowledge, this study is the first to assess whole blood gene expression and DNAm-gene expression associations in FASD. Our results present novel insights into the molecular footprint of FASD in whole blood and opens opportunities for future research into multi-omics biomarkers for the diagnosis of FASD.

Torch-eCpG: A fast and scalable eQTM mapper for thousands of molecular phenotypes with graphical processing units.

Background: Gene expression may be regulated by the DNA methylation of regulatory elements in cis, distal, and trans regions. One method to evaluate the relationship between DNA methylation and gene expression is the mapping of expression quantitative trait methylation (eQTM) loci (also called expression associated CpG loci, eCpG). However, no open-source tools are available to provide eQTM mapping. In addition, eQTM mapping can involve a large number of comparisons which may prevent the analyses due to limitations of computational resources. Here, we describe Torch-eCpG, an open-source tool to perform eQTM mapping that includes an optimized implementation that can use the graphical processing unit (GPU) to reduce runtime. Results: We demonstrate the analyses using the tool are reproducible, up to 18x faster using the GPU, and scale linearly with increasing methylation loci. Conclusions: Torch-eCpG is a fast, reliable, and scalable tool to perform eQTM mapping. Source code for Torch-eCpG is available at https://github.com/kordk/torch-ecpg.

Local diversity of drought resistance and resilience in Populus tomentosa correlates with the variation of DNA methylation.

Little information is known about DNA methylation variation in shaping environment-specific drought resistance and resilience for tree adaptation. In this study, we leveraged RNA sequencing and whole-genome bisulfite sequencing data to dissect the distinction of epigenetic regulation under drought stress and rewater condition of Populus tomentosa accessions from three geographical regions. We demonstrated low resistance and high resilience for accessions from South. Non-CG methylation levels in promoter regions of Southern accessions were lower than accessions from higher latitudes and negatively regulated gene expression. CHH context methylation was more sensitive to drought stress, and the geographical-specific differentially methylated regions were scarcely changed by environmental fluctuation. We identified 60 conserved hub genes within the co-expression networks that correlate with photosynthetic and stomatal morphological traits. Epigenome-wide association studies and genome-wide association studies of these 60 hub genes revealed the interdependency between genetic and epigenetic variation in GATA9 and LECRK-VIII.2, which was associated with stomatal morphology and chlorophyll content. The natural epigenetic variation in GATA9 was also faithfully transmitted to progenies in two family-based F1 populations. This study indicates a functional relationship of DNA methylation diversity with drought resistance and resilience which offers new insights into plants' local adaptation to a stressful environment.

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).

Cells can fine-tune which genes they activate, when and at which levels using a range of chemical marks on the DNA and certain proteins that help to organise the genome. One well-known example of such 'epigenetic tags' is DNA methylation, whereby a methyl group is added onto particular positions in the genome. Many factors ­ including environmental effects such as diet ­ control DNA methylation, allowing an organism to adapt to ever-changing conditions. An expression quantitative trait methylation (eQTM) is a specific position of the genome whose DNA methylation status regulates the activity of a given gene. A catalogue of eQTMs would be useful in helping to reveal how the environment and disease impacts the way cells work. Yet, currently, the relationships between most epigenetic tags and gene activity remains unclear, especially in children. To fill this gap, Ruiz-Arenas et al. studied DNA methylation in blood samples from over 800 healthy children across Europe. Amongst all tested DNA methylation sites, 22,000 (5.7% of total) were associated with the expression of a gene ­ and therefore were eQTMs; reciprocally, 9,000 genes (15.3% of all tested genes) were linked to at least one methylation site, leading to a total of 40,000 pairs of DNA methylation sites and genes. Most often, eQTMs regulated the expression of nearby genes ­ but only half controlled the gene that was the closest to them. Age and the genetic background of the individuals influenced the nature of eQTMs. This catalogue is a useful resource for the scientific community to start understanding the relationship between epigenetics and gene activity. Similar studies are now needed for other tissues and age ranges. Overall, extending our knowledge of eQTMs may help reveal how life events lead to illness, and could inform prevention efforts.

Prenatal NeuN+ neurons of Down syndrome display aberrant integrative DNA methylation and gene expression profiles.

Aim: To detect expression quantitative trait methylation (eQTM) loci within the cerebrum of prenatal Down syndrome (DS) and controls. Material & methods: DNA methylation gene expression profiles were acquired from NeuN+ nuclei, obtained from cerebrum sections of DS and controls. Linear regression models were applied to both datasets and were subsequently applied in an integrative analysis model to detect DS-associated eQTM loci. Results & conclusion: Widespread aberrant DNA methylation and gene expression were observed in DS. A substantial number of differentially methylated loci were replicated according to a previously reported study. Subsequent integrative analyses (eQTM) yielded numerous associated DS loci. the authors associated DNA methylation, gene expression and eQTM loci with DS that may underlie particular DS phenotypical characteristics.

Down syndrome (DS) is a common (1 of 1000 live births) autosomal aneuploidy in humans. Epigenetic programming regulates gene expression, defines cell fate and differentiation and drives early development. The authors aimed to detect DNA loci that are linked with the early developing brain of DS by analyzing DNA methylation and gene expression in prenatal DS neuronal samples. Numerous differential DNA methylated and expressed loci were found to be linked with DS. These findings may underlie particular DS characteristics, yet follow-up confirmation is required.

PTSD is associated with increased DNA methylation across regions of HLA-DPB1 and SPATC1L.

Posttraumatic stress disorder (PTSD) is characterized by intrusive thoughts, avoidance, negative alterations in cognitions and mood, and arousal symptoms that adversely affect mental and physical health. Recent evidence links changes in DNA methylation of CpG cites to PTSD. Since clusters of proximal CpGs share similar methylation signatures, identification of PTSD-associated differentially methylated regions (DMRs) may elucidate the pathways defining differential risk and resilience of PTSD. Here we aimed to identify epigenetic differences associated with PTSD. DNA methylation data profiled from blood samples using the MethylationEPIC BeadChip were used to perform a DMR analysis in 187 PTSD cases and 367 trauma-exposed controls from the Grady Trauma Project (GTP). DMRs were assessed with R package bumphunter. We identified two regions that associate with PTSD after multiple test correction. These regions were in the gene body of HLA-DPB1 and in the promoter of SPATC1L. The DMR in HLA-DPB1 was associated with PTSD in an independent cohort. Both DMRs included CpGs whose methylation associated with nearby sequence variation (meQTL) and that associated with expression of their respective genes (eQTM). This study supports an emerging literature linking PTSD risk to genetic and epigenetic variation in the HLA region.

Expression Quantitative Trait Methylation Analysis Reveals Methylomic Associations With Gene Expression in Childhood Asthma.

BACKGROUND: Nasal (airway) epithelial methylation profiles have been associated with asthma, but the effects of such profiles on expression of distant cis-genes are largely unknown. RESEARCH QUESTION: To identify genes whose expression is associated with proximal and distal CpG probes (within 1 Mb), and to assess whether and how such genes are differentially expressed in atopic asthma. STUDY DESIGN AND METHODS: Genome-wide expression quantitative trait methylation (eQTM) analysis in nasal epithelium from Puerto Rican subjects (aged 9-20 years) with (n = 219) and without (n = 236) asthma. After the eQTM analysis, a Gene Ontology Enrichment analysis was conducted for the top 500 eQTM genes, and mediation analyses were performed to identify paths from DNA methylation to atopic asthma through gene expression. Asthma was defined as physician-diagnosed asthma and wheeze in the previous year, and atopy was defined as at least one positive IgE to allergens. Atopic asthma was defined as the presence of both atopy and asthma. RESULTS: We identified 16,867 significant methylation-gene expression pairs (false-discovery rate-adjusted P < .01) in nasal epithelium from study participants. Most eQTM methylation probes were distant (average distance, ∼378 kb) from their target genes, and also more likely to be located in enhancer regions of their target genes in lung tissue than control probes. The top 500 eQTM genes were enriched in pathways for immune processes and epithelial integrity and were more likely to have been previously identified as differentially expressed in atopic asthma. In a mediation analysis, we identified 5,934 paths through which methylation markers could affect atopic asthma through gene expression in nasal epithelium. INTERPRETATION: Previous epigenome-wide association studies of asthma have estimated the effects of DNA methylation markers on expression of nearby genes in airway epithelium. Our findings suggest that distant epigenetic regulation of gene expression in airway epithelium plays a role in atopic asthma.

DNA methylation loci in placenta associated with birthweight and expression of genes relevant for early development and adult diseases.

BACKGROUND: Birthweight marks an important milestone of health across the lifespan, including cardiometabolic disease risk in later life. The placenta, a transient organ at the maternal-fetal interface, regulates fetal growth. Identifying genetic loci where DNA methylation in placenta is associated with birthweight can unravel genomic pathways that are dysregulated in aberrant fetal growth and cardiometabolic diseases in later life. RESULTS: We performed placental epigenome-wide association study (EWAS) of birthweight in an ethnic diverse cohort of pregnant women (n = 301). Methylation at 15 cytosine-(phosphate)-guanine sites (CpGs) was associated with birthweight (false discovery rate (FDR) < 0.05). Methylation at four (26.7%) CpG sites was associated with placental transcript levels of 15 genes (FDR < 0.05), including genes known to be associated with adult lipid traits, inflammation and oxidative stress. Increased methylation at cg06155341 was associated with higher birthweight and lower FOSL1 expression, and lower FOSL1 expression was correlated with higher birthweight. Given the role of the FOSL1 transcription factor in regulating developmental processes at the maternal-fetal interface, epigenetic mechanisms at this locus may regulate fetal development. We demonstrated trans-tissue portability of methylation at four genes (MLLT1, PDE9A, ASAP2, and SLC20A2) implicated in birthweight by a previous study in cord blood. We also found that methylation changes known to be related to maternal underweight, preeclampsia and adult type 2 diabetes were associated with lower birthweight in placenta. CONCLUSION: We identified novel placental DNA methylation changes associated with birthweight. Placental epigenetic mechanisms may underlie dysregulated fetal development and early origins of adult cardiometabolic diseases. CLINICAL TRIAL REGISTRATION: ClinicalTrials.gov, NCT00912132.

Epigenome-wide association study of DNA methylation and microRNA expression highlights novel pathways for human complex traits.

DNA methylation (DNAm) and microRNAs (miRNAs) have been implicated in a wide-range of human diseases. While often studied in isolation, DNAm and miRNAs are not independent. We analyzed associations of expression of 283 miRNAs with DNAm at >400K CpG sites in whole blood obtained from 3565 individuals and identified 227 CpGs at which differential methylation was associated with the expression of 40 nearby miRNAs (cis-miR-eQTMs) at FDR<0.01, including 91 independent CpG sites at r2 < 0.2. cis-miR-eQTMs were enriched for CpGs in promoter and polycomb-repressed state regions, and 60% were inversely associated with miRNA expression. Bidirectional Mendelian randomization (MR) analysis further identified 58 cis-miR-eQTMCpG-miRNA pairs where DNAm changes appeared to drive miRNA expression changes and opposite directional effects were unlikely. Integration of genetic variants in joint analyses revealed an average partial between cis-miR-eQTM CpGs and miRNAs of 2% after conditioning on site-specific genetic variation, suggesting that DNAm is an important epigenetic regulator of miRNA expression. Finally, two-step MR analysis was performed to identify putatively causal CpGs driving miRNA expression in relation to human complex traits. We found that an imprinted region on 14q32 that was previously identified in relation to age at menarche is enriched with cis-miR-eQTMs. Nine CpGs and three miRNAs at this locus tested causal for age at menarche, reflecting novel epigenetic-driven molecular pathways underlying this complex trait. Our study sheds light on the joint genetic and epigenetic regulation of miRNA expression and provides insights into the relations of miRNAs to their targets and to complex phenotypes.

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 F1 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.