Fetal DNA Methylation Associates with Early Spontaneous Preterm Birth and Gestational Age

Spontaneous preterm birth (PTB, <37 weeks gestation) is a major public health concern, and children born preterm have a higher risk of morbidity and mortality throughout their lives. Recent studies suggest that fetal DNA methylation of several genes varies across a range of gestational ages (GA), but it is not yet clear if fetal epigenetic changes associate with PTB. The objective of this study is to interrogate methylation patterns across the genome in fetal leukocyte DNA from African Americans with early PTB (24^1/7–34^0/7 weeks; N = 22) or term births (39^0/7–40^6/7weeks; N = 28) and to evaluate the association of each CpG site with PTB and GA. DNA methylation was assessed across the genome with the HumanMethylation450 BeadChip. For each individual sample and CpG site, the proportion of DNA methylation was estimated. The associations between methylation and PTB or GA were evaluated by fitting a separate linear model for each CpG site, adjusting for relevant covariates. Overall, 29 CpG sites associated with PTB (FDR<.05; 5.7×10⁻¹⁰<p<2.9×10⁻⁶) independent of GA. Also, 9637 sites associated with GA (FDR<.05; 9.5×10⁻¹⁶<p<1.0×10⁻³), with 61.8% decreasing in methylation with shorter GA. GA-associated CpG sites were depleted in the CpG islands of their respective genes (p<2.2×10⁻¹⁶). Gene set enrichment analysis (GSEA) supported enrichment of GA-associated CpG sites in genes that play a role in embryonic development as well as the extracellular matrix. Additionally, this study replicated the association of several CpG sites associated with gestational age in other studies (CRHBP, PIK3CD and AVP). Dramatic differences in fetal DNA methylation are evident in fetuses born preterm versus at term, and the patterns established at birth may provide insight into the long-term consequences associated with PTB.     

Genetic associations for two biological age measures point to distinct aging phenotypes

Biological age measures outperform chronological age in predicting various aging outcomes, yet little is known regarding genetic predisposition. We performed genome‐wide association scans of two age‐adjusted biological age measures (PhenoAgeAcceleration and BioAgeAcceleration), estimated from clinical biochemistry markers (Levine et al., 2018; Levine, 2013) in European‐descent participants from UK Biobank. The strongest signals were found in the APOE gene, tagged by the two major protein‐coding SNPs, PhenoAgeAccel—rs429358 (APOE e4 determinant) (p = 1.50 × 10⁻⁷²); BioAgeAccel—rs7412 (APOE e2 determinant) (p = 3.16 × 10⁻⁶⁰). Interestingly, we observed inverse APOE e2 and e4 associations and unique pathway enrichments when comparing the two biological age measures. Genes associated with BioAgeAccel were enriched in lipid related pathways, while genes associated with PhenoAgeAccel showed enrichment for immune system, cell function, and carbohydrate homeostasis pathways, suggesting the two measures capture different aging domains. Our study reaffirms that aging patterns are heterogeneous across individuals, and the manner in which a person ages may be partly attributed to genetic predisposition.     

Clonal hematopoiesis associated with epigenetic aging and clinical outcomes

Clonal hematopoiesis of indeterminate potential (CHIP) is a common precursor state for blood cancers that most frequently occurs due to mutations in the DNA‐methylation modifying enzymes DNMT3A or TET2. We used DNA‐methylation array and whole‐genome sequencing data from four cohorts together comprising 5522 persons to study the association between CHIP, epigenetic clocks, and health outcomes. CHIP was strongly associated with epigenetic age acceleration, defined as the residual after regressing epigenetic clock age on chronological age, in several clocks, ranging from 1.31 years (GrimAge, p < 8.6 × 10⁻⁷) to 3.08 years (EEAA, p < 3.7 × 10⁻¹⁸). Mutations in most CHIP genes except DNA‐damage response genes were associated with increases in several measures of age acceleration. CHIP carriers with mutations in multiple genes had the largest increases in age acceleration and decrease in estimated telomere length. Finally, we found that ~40% of CHIP carriers had acceleration >0 in both Hannum and GrimAge (referred to as AgeAccelHG+). This group was at high risk of all‐cause mortality (hazard ratio 2.90, p < 4.1 × 10⁻⁸) and coronary heart disease (CHD) (hazard ratio 3.24, p < 9.3 × 10⁻⁶) compared to those who were CHIP−/AgeAccelHG−. In contrast, the other ~60% of CHIP carriers who were AgeAccelHG− were not at increased risk of these outcomes. In summary, CHIP is strongly linked to age acceleration in multiple clocks, and the combination of CHIP and epigenetic aging may be used to identify a population at high risk for adverse outcomes and who may be a target for clinical interventions.     

Fetal lung and placental methylation is associated with in utero nicotine exposure

In utero smoke exposure has been shown to have detrimental effects on lung function and to be associated with persistent wheezing and asthma in children. One potential mechanism of IUS effects could be alterations in DNA methylation, which may have life-long implications. The goal of this study was to examine the association between DNA methylation and nicotine exposure in fetal lung and placental tissue in early development; nicotine exposure in this analysis represents a likely surrogate for in-utero smoke. We performed an epigenome-wide analysis of DNA methylation in fetal lung tissue (n = 85, 41 smoke exposed (48%), 44 controls) and the corresponding placental tissue samples (n = 80, 39 smoke exposed (49%), 41 controls) using the Illumina HumanMethylation450 BeadChip array. Differential methylation analyses were conducted to evaluate the variation associated with nicotine exposure. The most significant CpG sites in the fetal lung analysis mapped to the PKP3 (P = 2.94 × 10⁻⁰³), ANKRD33B (P = 3.12 × 10⁻⁰³), CNTD2 (P = 4.9 × 10⁻⁰³) and DPP10 (P = 5.43 × 10⁻⁰³) genes. In the placental methylome, the most significant CpG sites mapped to the GTF2H2C and GTF2H2D genes (P = 2.87 × 10⁻⁰⁶ − 3.48 × 10⁻⁰⁵). One hundred and one unique CpG sites with P-values < 0.05 were concordant between lung and placental tissue analyses. Gene Set Enrichment Analysis demonstrated enrichment of specific disorders, such as asthma and immune disorders. Our findings demonstrate an association between in utero nicotine exposure and variable DNA methylation in fetal lung and placental tissues, suggesting a role for DNA methylation variation in the fetal origins of chronic diseases.     

DNA methylation levels and long-term trihalomethane exposure in drinking water: an epigenome-wide association study

Trihalomethanes (THM) are undesired disinfection byproducts (DBPs) formed during water treatment. Mice exposed to DBPs showed global DNA hypomethylation and c-myc and c-jun gene-specific hypomethylation, while evidence of epigenetic effects in humans is scarce. We explored the association between lifetime THM exposure and DNA methylation through an epigenome-wide association study. We selected 138 population-based controls from a case-control study of colorectal cancer conducted in Barcelona, Spain, exposed to average lifetime THM levels ≤85 μg/L vs. >85 μg/L (N = 68 and N = 70, respectively). Mean age of participants was 70 years, and 54% were male. Average lifetime THM level in the exposure groups was 64 and 130 µg/L, respectively. DNA was extracted from whole blood and was bisulphite converted to measure DNA methylation levels using the Illumina HumanMethylation450 BeadChip. Data preprocessing was performed using RnBeads. Methylation was compared between exposure groups using empirical Bayes moderated linear regression for CpG sites and Gaussian kernel for CpG regions. ConsensusPathDB was used for gene set enrichment. Statistically significant differences in methylation between exposure groups was found in 140 CpG sites and 30 gene-related regions, after false discovery rate <0.05 and adjustment for age, sex, methylation first principal component, and blood cell proportion. The annotated genes were localized to several cancer pathways. Among them, 29 CpGs had methylation levels associated with THM levels (|Δβ|≥0.05) located in 11 genes associated with cancer in other studies. Our results suggest that THM exposure may affect DNA methylation in genes related to tumors, including colorectal and bladder cancers. Future confirmation studies are required.     

Key epigenetic changes associated with lung cancer development: Results from dense methylation array profiling

Epigenetic alterations are a common event in lung cancer and their identification can serve to inform on the carcinogenic process and provide clinically relevant biomarkers. Using paired tumor and non-tumor lung tissues from 146 individuals from three independent populations we sought to identify common changes in DNA methylation associated with the development of non-small cell lung cancer. Pathologically normal lung tissue taken at the time of cancer resection was matched to tumorous lung tissue and together were probed for methylation using Illumina GoldenGate arrays in the discovery set (n = 47 pairs) followed by bisulfite pyrosequencing for validation sets (n = 99 pairs). For each matched pair the change in methylation at each CpG was calculated (the odds ratio), and these ratios were averaged across individuals and ranked by magnitude to identify the CpGs with the greatest change in methylation associated with tumor development. We identified the top gene-loci representing an increase in methylation (HOXA9, 10.3-fold and SOX1, 5.9-fold) and decrease in methylation (DDR1, 8.1-fold). In replication testing sets, methylation was higher in tumors for HOXA9 (p < 2.2 × 10⁻¹⁶) and SOX1 (p < 2.2 × 10⁻¹⁶) and lower for DDR1 (p < 2.2 × 10⁻¹⁶). The magnitude and strength of these changes were consistent across squamous cell and adenocarcinoma tumors. Our data indicate that the identified genes consistently have altered methylation in lung tumors. Our identified genes should be included in translational studies that aim to develop screening for early disease detection.     

A Genome-Wide Methylation Study on Essential Hypertension in Young African American Males

Objective

There is emerging evidence from animal studies suggesting a key role for methylation in the pathogenesis of essential hypertension. However, to date, very few studies have investigated the role of methylation in the development of human hypertension, and none has taken a genome-wide approach. Based on the recent studies that highlight the involvement of inflammation in the development of hypertension, we hypothesize that changes in DNA methylation of leukocytes are involved in the pathogenesis of hypertension.

Method & Results

We conducted a genome-wide methylation analysis on 8 hypertensive cases and 8 normotensive age-matched controls aged 14–23 years and performed validation of the most significant CpG sites in 2 genes in an independent sample of 36 hypertensive cases and 60 normotensive controls aged 14–30 years. Validation of the CpG sites in the SULF1 gene was further conducted in a second replication sample of 36 hypertensive cases and 34 controls aged 15.8–40 years. A CpG site in the SULF1 gene showed higher methylation levels in cases than in healthy controls in the genome-wide step (p = 6.2×10⁻⁵), which was confirmed in the validation step (p = 0.011) for subjects ≤30 years old but was not significant for subjects of all ages combined (p = 0.095).

Conclusion

The identification of a difference in a blood leukocyte DNA methylation site between hypertensive cases and normotensive controls suggests that changes in DNA methylation may play an important role in the pathogenesis of hypertension. The age dependency of the effect further suggests complexity of epigenetic regulation in this age-related disease.     

Methylation at CPT1A locus is associated with lipoprotein subfraction profiles

Lipoprotein subfractions help discriminate cardiometabolic disease risk. Genetic loci validated as associating with lipoprotein measures do not account for a large proportion of the individual variation in lipoprotein measures. We hypothesized that DNA methylation levels across the genome contribute to interindividual variation in lipoprotein measures. Using data from participants of the Genetics of Lipid Lowering Drugs and Diet Network (n = 663 for discovery and n = 331 for replication stages, respectively), we conducted the first systematic screen of the genome to determine associations between methylation status at ∼470,000 cytosine-guanine dinucleotide (CpG) sites in CD4⁺ T cells and 14 lipoprotein subfraction measures. We modeled associations between methylation at each CpG site and each lipoprotein measure separately using linear mixed models, adjusted for age, sex, study site, cell purity, and family structure. We identified two CpGs, both in the carnitine palmitoyltransferase-1A (CPT1A) gene, which reached significant levels of association with VLDL and LDL subfraction parameters in both discovery and replication phases (P < 1.1 × 10⁻⁷ in the discovery phase, P < .004 in the replication phase, and P < 1.1 × 10⁻¹² in the full sample). CPT1A is regulated by PPARα, a ligand for drugs used to reduce CVD. Our associations between methylation in CPT1A and lipoprotein measures highlight the epigenetic role of this gene in metabolic dysfunction.     

Genome-Wide Investigation of DNA Methylation Marks Associated with FV Leiden Mutation

In order to investigate whether DNA methylation marks could contribute to the incomplete penetrance of the FV Leiden mutation, a major genetic risk factor for venous thrombosis (VT), we measured genome-wide DNA methylation levels in peripheral blood samples of 98 VT patients carrying the mutation and 251 VT patients without the mutation using the dedicated Illumina HumanMethylation450 array. The genome-wide analysis of 388,120 CpG probes identified three sites mapping to the SLC19A2 locus whose DNA methylation levels differed significantly (p<3 10⁻⁸) between carriers and non-carriers. The three sites replicated (p<2 10⁻⁷) in an independent sample of 214 individuals from five large families ascertained on VT and FV Leiden mutation among which 53 were carriers and 161 were non-carriers of the mutation. In both studies, these three CpG sites were also associated (2.33 10⁻¹¹<p<3.02 10⁻⁴) with biomarkers of the Protein C pathway known to be influenced by the FV Leiden mutation. A comprehensive linkage disequilibrium (LD) analysis of the whole locus revealed that the original associations were due to LD between the FV Leiden mutation and a block of single nucleotide polymorphisms (SNP) located in SLC19A2. After adjusting for this block of SNPs, the FV Leiden mutation was no longer associated with any CpG site (p>0.05). In conclusion, our work clearly illustrates some promises and pitfalls of DNA methylation investigations on peripheral blood DNA in large epidemiological cohorts. DNA methylation levels at SLC19A2 are influenced by SNPs in LD with FV Leiden, but these DNA methylation marks do not explain the incomplete penetrance of the FV Leiden mutation.     

Human-specific CpG “beacons” identify loci associated with human-specific traits and disease

Regulatory change has long been hypothesized to drive the delineation of the human phenotype from other closely related primates. Here we provide evidence that CpG dinucleotides play a special role in this process. CpGs enable epigenome variability via DNA methylation, and this epigenetic mark functions as a regulatory mechanism. Therefore, species-specific CpGs may influence species-specific regulation. We report non-polymorphic species-specific CpG dinucleotides (termed “CpG beacons”) as a distinct genomic feature associated with CpG island (CGI) evolution, human traits and disease. Using an inter-primate comparison, we identified 21 extreme CpG beacon clusters (≥ 20/kb peaks, empirical p < 1.0 × 10⁻³) in humans, which include associations with four monogenic developmental and neurological disease related genes (Benjamini-Hochberg corrected p = 6.03 × 10⁻³). We also demonstrate that beacon-mediated CpG density gain in CGIs correlates with reduced methylation in these species in orthologous CGIs over time, via human, chimpanzee and macaque MeDIP-seq. Therefore mapping into both the genomic and epigenomic space the identified CpG beacon clusters define points of intersection where a substantial two-way interaction between genetic sequence and epigenetic state has occurred. Taken together, our data support a model for CpG beacons to contribute to CGI evolution from genesis to tissue-specific to constitutively active CGIs.     

ELOVL2 gene polymorphisms are associated with increases in plasma eicosapentaenoic and docosahexaenoic acid proportions after fish oil supplement

Fish oil supplementation provides an inconsistent degree of protection from cardiovascular disease (CVD), which may be attributed to genetic variation. Single nucleotide polymorphisms (SNPs) in the elongation-of-very-long-chain-fatty-acids-2 (ELOVL2) gene have been strongly associated with plasma proportions of n-3 long-chain polyunsaturated fatty acids (LC-PUFA). We investigated the effect of genotype interaction with fish oil dosage on plasma n-3 LC-PUFA proportions in a parallel double-blind controlled trial, involving 367 subjects randomised to treatment with 0.45, 0.9 and 1.8 g/day eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) (1.51:1) or olive oil placebo for 6 months. We genotyped 310 subjects for ELOVL2 gene SNPs rs3734398, rs2236212 and rs953413. At baseline, carriers of all minor alleles had lower proportions of plasma DHA than non-carriers (P = 0.021–0.030). Interaction between genotype and treatment was a significant determinant of plasma EPA (P < 0.0001) and DHA (P = 0.004–0.032). After the 1.8 g/day dose, carriers of ELOVL2 SNP minor alleles had approximately 30 % higher proportions of EPA (P = 0.002–0.004) and 9 % higher DHA (P = 0.013–0.017) than non-carriers. Minor allele carriers could therefore particularly benefit from a high intake of EPA and DHA in maintaining high levels of plasma n-3 PUFA conducive to protection from CVD.

Electronic supplementary material

The online version of this article (doi:10.1007/s12263-013-0362-6) contains supplementary material, which is available to authorized users.     

The impact of genetic variation and cigarette smoke on DNA methylation in current and former smokers from the COPDGene study

DNA methylation can be affected by systemic exposures, such as cigarette smoking and genetic sequence variation; however, the relative impact of each on the epigenome is unknown. We aimed to assess if cigarette smoking and genetic variation are associated with overlapping or distinct sets of DNA methylation marks and pathways. We selected 85 Caucasian current and former smokers with genome-wide single nucleotide polymorphism (SNP) genotyping available from the COPDGene study.  Genome-wide methylation was obtained on DNA from whole blood using the Illumina HumanMethylation27 platform. To determine the impact of local sequence variation on DNA methylation (mQTL), we examined the association between methylation and SNPs within 50 kb of each CpG site.  To examine the impact of cigarette smoking on DNA methylation, we examined the differences in methylation by current cigarette smoking status. We detected 770 CpG sites annotated to 708 genes associated at an FDR < 0.05 in the cis-mQTL analysis and 1,287 CpG sites annotated to 1,242 genes, which were nominally associated in the smoking-CpG association analysis (P_unadjusted < 0.05). Forty-three CpG sites annotated to 40 genes were associated with both SNP variation and current smoking; this overlap was not greater than that expected by chance. Our results suggest that cigarette smoking and genetic variants impact distinct sets of DNA methylation marks, the further elucidation of which may partially explain the variable susceptibility to the health effects of cigarette smoking. Ascertaining how genetic variation and systemic exposures differentially impact the human epigenome has relevance for both biomarker identification and therapeutic target development for smoking-related diseases.     

MLH1 Region Polymorphisms Show a Significant Association with CpG Island Shore Methylation in a Large Cohort of Healthy Individuals

Single nucleotide polymorphisms (SNPs) are the most common form of genetic variation. We previously demonstrated that SNPs (rs1800734, rs749072, and rs13098279) in the MLH1 gene region are associated with MLH1 promoter island methylation, loss of MLH1 protein expression, and microsatellite instability (MSI) in colorectal cancer (CRC) patients. Recent studies have identified less CpG-dense “shore” regions flanking many CpG islands. These shores often exhibit distinct methylation profiles between different tissues and matched normal versus tumor cells of patients. To date, most epigenetic studies have focused on somatic methylation events occurring within solid tumors; less is known of the contributions of peripheral blood cell (PBC) methylation to processes such as aging and tumorigenesis. To address whether MLH1 methylation in PBCs is correlated with tumorigenesis we utilized the Illumina 450 K microarrays to measure methylation in PBC DNA of 846 healthy controls and 252 CRC patients from Ontario, Canada. Analysis of a region of chromosome 3p21 spanning the MLH1 locus in healthy controls revealed that a CpG island shore 1 kb upstream of the MLH1 gene exhibits different methylation profiles when stratified by SNP genotypes (rs1800734, rs749072, and rs13098279). Individuals with wild-type genotypes incur significantly higher PBC shore methylation than heterozygous or homozygous variant carriers (p<1.1×10⁻⁶; ANOVA). This trend is also seen in CRC cases (p<0.096; ANOVA). Shore methylation also decreases significantly with increasing age in cases and controls. This is the first study of its kind to integrate PBC methylation at a CpG island shore with SNP genotype status in CRC cases and controls. These results indicate that CpG island shore methylation in PBCs may be influenced by genotype as well as the normal aging process.     

Genetic variants modify the effect of age on APOE methylation in the Genetics of Lipid Lowering Drugs and Diet Network study

Although apolipoprotein E (APOE) variants are associated with age-related diseases, the underlying mechanism is unknown and DNA methylation may be a potential one. With methylation data, measured by the Infinium Human Methylation 450 array, from 993 participants (age ranging from 18 to 87 years) in the Genetics of Lipid Lowering Drugs and Diet Network (GOLDN) study, and from Encyclopedia of DNA Elements (ENCODE) consortium, combined with published methylation datasets, we described the methylation pattern of 13 CpG sites within APOE locus, their correlations with gene expression across cell types, and their relationships with age, plasma lipids, and sequence variants. Based on methylation levels and the genetic regions, we categorized the 13 APOE CpG sites into three groups: Group 1 showed hypermethylation (> 50%) and were located in the promoter region, Group 2 exhibited hypomethylation (< 50%) and were located in the first two exons and introns, and Group 3 showed hypermethylation (> 50%) and were located in the exon 4. APOE methylation was negatively correlated with gene expression (minimum r = −0.66, P = 0.004). APOE methylation was significantly associated with age (minimum P = 2.06E-08) and plasma total cholesterol (minimum P = 3.53E-03). Finally, APOE methylation patterns differed across APOE ε variants (minimum P = 3.51E-05) and the promoter variant rs405509 (minimum P = 0.01), which further showed a significant interaction with age (P = 0.03). These findings suggest that methylation may be a potential mechanistic explanation for APOE functions related to aging and call for further molecular mechanistic studies.     

Increased P16 DNA Methylation in Mouse Thymic Lymphoma Induced by Irradiation

DNA methylation is an important part of epigenetics. In this study, we examined the methylation state of two CpG islands in the promoter of the p16 gene in radiation-induced thymic lymphoma samples. The mRNA and protein levels of P16 were significantly reduced in radiation-induced thymic lymphoma tissue samples. Twenty-three CpG sites of the CpG islands in the p16 promoter region were detected, and the methylation percentages of −71, −63, −239, −29, −38, −40, −23, 46 CpG sites were significantly higher in radiation-induced thymic lymphoma tissue samples than those in matched non-irradiated thymus tissue samples. This study provides new evidence for the methylation state of p16 in the radiation-induced thymic lymphoma samples, which suggests that the methylation of these CpG sites in the p16 promoter may reduce its expression in the thymic lymphoma after irradiation.