LINE-1 methylation in granulocyte DNA and trihalomethane exposure is associated with bladder cancer risk

DNA methylation changes contribute to bladder carcinogenesis. Trihalomethanes (THM), a class of disinfection by-products, are associated with increased urothelial bladder cancer (UBC) risk. THM exposure in animal models produces DNA hypomethylation. We evaluated the relationship of LINE-1 5-methylcytosine levels (LINE-1%5mC) as outcome of long-term THM exposure among controls and as an effect modifier in the association between THM exposure and UBC risk. We used a case-control study of UBC conducted in Spain. We obtained personal lifetime residential THM levels and measured LINE-1%5mC by pyrosequencing in granulocyte DNA from blood samples in 548 incident cases and 559 hospital controls. Two LINE-1%5mC clusters (above and below 64%) were identified through unsupervised hierarchical cluster analysis. The association between THM levels and LINE-1%5mC was evaluated with β regression analyses and logistic regression was used to estimate odds ratios (OR) adjusting for covariables. LINE-1%5mC change between percentiles 75^th and 25^th of THM levels was 1.8% (95% confidence interval (CI): 0.1, 3.4%) among controls. THM levels above vs. below the median (26 μg/L) were associated with increased UBC risk, OR = 1.86 (95% CI: 1.25, 2.75), overall and among subjects with low levels of LINE-1%5mC (n = 975), OR = 2.14 (95% CI: 1.39, 3.30), but not associated with UBC risk among subjects’ high levels of LINE-1%5mC (n = 162), interaction P = 0.03. Results suggest a positive association between LINE-1%5mC and THM levels among controls, and LINE-1%5mC status may modify the association between UBC risk and THM exposure. Because reverse causation and chance cannot be ruled out, confirmation studies are warranted.     

LINE-1 methylation in granulocyte DNA and trihalomethane exposure is associated with bladder cancer risk

DNA methylation changes contribute to bladder carcinogenesis. Trihalomethanes (THM), a class of disinfection by-products, are associated with increased urothelial bladder cancer (UBC) risk. THM exposure in animal models produces DNA hypomethylation. We evaluated the relationship of LINE-1 5-methylcytosine levels (LINE-1%5mC) as outcome of long-term THM exposure among controls and as an effect modifier in the association between THM exposure and UBC risk. We used a case-control study of UBC conducted in Spain. We obtained personal lifetime residential THM levels and measured LINE-1%5mC by pyrosequencing in granulocyte DNA from blood samples in 548 incident cases and 559 hospital controls. Two LINE-1%5mC clusters (above and below 64%) were identified through unsupervised hierarchical cluster analysis. The association between THM levels and LINE-1%5mC was evaluated with β regression analyses and logistic regression was used to estimate odds ratios (OR) adjusting for covariables. LINE-1%5mC change between percentiles 75^th and 25^th of THM levels was 1.8% (95% confidence interval (CI): 0.1, 3.4%) among controls. THM levels above vs. below the median (26 μg/L) were associated with increased UBC risk, OR = 1.86 (95% CI: 1.25, 2.75), overall and among subjects with low levels of LINE-1%5mC (n = 975), OR = 2.14 (95% CI: 1.39, 3.30), but not associated with UBC risk among subjects’ high levels of LINE-1%5mC (n = 162), interaction P = 0.03. Results suggest a positive association between LINE-1%5mC and THM levels among controls, and LINE-1%5mC status may modify the association between UBC risk and THM exposure. Because reverse causation and chance cannot be ruled out, confirmation studies are warranted.     

Imputation of missing covariate values in epigenome-wide analysis of DNA methylation data

DNA methylation is a widely studied epigenetic mechanism and alterations in methylation patterns may be involved in the development of common diseases. Unlike inherited changes in genetic sequence, variation in site-specific methylation varies by tissue, developmental stage, and disease status, and may be impacted by aging and exposure to environmental factors, such as diet or smoking. These non-genetic factors are typically included in epigenome-wide association studies (EWAS) because they may be confounding factors to the association between methylation and disease. However, missing values in these variables can lead to reduced sample size and decrease the statistical power of EWAS. We propose a site selection and multiple imputation (MI) method to impute missing covariate values and to perform association tests in EWAS. Then, we compare this method to an alternative projection-based method. Through simulations, we show that the MI-based method is slightly conservative, but provides consistent estimates for effect size. We also illustrate these methods with data from the Atherosclerosis Risk in Communities (ARIC) study to carry out an EWAS between methylation levels and smoking status, in which missing cell type compositions and white blood cell counts are imputed.     

DNA methylation and obesity traits: An epigenome-wide association study. The REGICOR study

Obesity is associated with increased risk of several diseases and has become epidemic. Obesity is highly heritable but the genetic variants identified by genome-wide association studies explain only limited variability. Epigenetics could contribute to explain the missing variability. The study aim was to discover differential methylation patterns related to obesity. We designed an epigenome-wide association study with a discovery phase in a subsample of 641 REGICOR study participants, validated by analysis of 2,515 participants in the Framingham Offspring Study. Blood DNA methylation was assessed using Illumina HumanMethylation450 BeadChip. Next, we meta-analyzed the data using the fixed effects method and performed a functional and pathway analysis using the Ingenuity Pathway Analysis software. We were able to validate 94 CpGs associated with body mass index (BMI) and 49 CpGs associated with waist circumference, located in 95 loci. In addition, we newly discovered 70 CpGs associated with BMI and 33 CpGs related to waist circumference. These CpGs explained 25.94% and 29.22% of the variability of BMI and waist circumference, respectively, in the REGICOR sample. We also evaluated 65 of the 95 validated loci in the GIANT genome-wide association data; 10 of them had Tag SNPs associated with BMI. The top-ranked diseases and functions identified in the functional and pathway analysis were neurologic, psychological, endocrine, and metabolic.     

Differential DNA methylation and PM2.5 species in a 450K epigenome-wide association study

Although there is growing evidence that exposure to ambient particulate matter is associated with global DNA methylation and gene-specific methylation, little is known regarding epigenome-wide changes in DNA methylation in relation to particles and, especially, particle components. Using the Illumina Infinium HumanMethylation450 BeadChip, we examined the relationship between one-year moving averages of PM_2.5 species (Al, Ca, Cu, Fe, K, Na, Ni, S, Si, V, and Zn) and DNA methylation at 484,613 CpG probes in a longitudinal cohort that included 646 subjects. Bonferroni correction was applied to adjust for multiple comparisons. Bioinformatics analysis of the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment was also performed. We observed 20 Bonferroni significant (P-value < 9.4× 10^?9) CpGs for Fe, 8 for Ni, and 1 for V. Particularly, methylation at Schlafen Family Member 11 (SLFN11) cg10911913 was positively associated with measured levels of all 3 species. The SLFN11 gene codes for an interferon-induced protein that inhibits retroviruses and sensitizes cancer cells to DNA-damaging agents. Bioinformatics analysis suggests that gene targets may be relevant to pathways including cancers, signal transduction, and cell growth and death. Ours is the first study to examine the epigenome-wide association between ambient particles species and DNA methylation. We found that long-term exposures to specific components of ambient particle pollution, especially particles emitted during oil combustion, were associated with methylation changes in genes relevant to immune responses. Our findings provide insight into potential biologic mechanisms on an epigenetic level.     

In utero arsenic exposure and epigenome-wide associations in placenta,umbilical artery,and human umbilical vein endothelial cells

Exposure to arsenic early in life has been associated with increased risk of several chronic diseases and is believed to alter epigenetic programming in utero. In the present study, we evaluate the epigenome-wide association of arsenic exposure in utero and DNA methylation in placenta (n = 37), umbilical artery (n = 45) and human umbilical vein endothelial cells (HUVEC) (n = 52) in a birth cohort using the Infinium HumanMethylation450 BeadChip array. Unadjusted and cell mixture adjusted associations for each tissue were examined along with enrichment analyses relative to CpG island location and omnibus permutation tests of association among biological pathways. One CpG in artery (cg26587014) and 4 CpGs in placenta (cg12825509; cg20554753; cg23439277; cg21055948) reached a Bonferroni adjusted level of significance. Several CpGs were differentially methylated in artery and placenta when controlling the false discovery rate (q-value<0.05), but none in HUVEC. Enrichment of hypomethylated CpG islands was observed for artery while hypermethylation of open sea regions were present in placenta relative to prenatal arsenic exposure. The melanogenesis pathway was differentially methylated in artery (Max F P < 0.001), placenta (Max F P < 0.001), and HUVEC (Max F P = 0.02). Similarly, the insulin-signaling pathway was differentially methylated in artery (Max F P = 0.02), placenta (Max F P = 0.02), and HUVEC (Max F P = 0.02). Our results show that prenatal arsenic exposure can alter DNA methylation in artery and placenta but not in HUVEC. Further studies are needed to determine if these alterations in DNA methylation mediate the effect of prenatal arsenic exposure and health outcomes later in life.     

Differential DNA methylation in umbilical cord blood of infants exposed to mercury and arsenic in utero

Mercury and arsenic are known developmental toxicants. Prenatal exposures are associated with adverse childhood health outcomes that could be in part mediated by epigenetic alterations that may also contribute to altered immune profiles. In this study, we examined the association between prenatal mercury exposure on both DNA methylation and white blood cell composition of cord blood, and evaluated the interaction with prenatal arsenic exposure. A total of 138 mother-infant pairs with postpartum maternal toenail mercury, prenatal urinary arsenic concentrations, and newborn cord blood were assessed using the Illumina Infinium Methylation450 array. White blood cell composition was inferred from DNA methylation measurements. A doubling in toenail mercury concentration was associated with a 2.5% decrease (95% CI: 5.0%, 1.0%) in the estimated monocyte proportion. An increase of 3.5% (95% CI: 1.0, 7.0) in B-cell proportion was observed for females only. Among the top 100 CpGs associated with toenail mercury levels (ranked on P-value), there was a significant enrichment of loci located in North shore regions of CpG islands (P = 0.049), and the majority of these loci were hypermethylated (85%). Among the top 100 CpGs for the interaction between arsenic and mercury, there was a greater than expected proportion of loci located in CpG islands (P = 0.045) and in South shore regions (P = 0.009) and all of these loci were hypermethylated. This work supports the hypothesis that mercury may be contributing to epigenetic variability and immune cell proportion changes, and suggests that in utero exposure to mercury and arsenic, even at low levels, may interact to impact the epigenome.     

Improved reproducibility in genome-wide DNA methylation analysis for PAXgene-fixed samples compared with restored formalin-fixed and paraffin-embedded DNA

Formalin fixation has been the standard method for conservation of clinical specimens for decades. However, a major drawback is the high degradation of nucleic acids, which complicates its use in genome-wide analyses. Unbiased identification of biomarkers, however, requires genome-wide studies, precluding the use of the valuable archives of specimens with long-term follow-up data. Therefore, restoration protocols for DNA from formalin-fixed and paraffin-embedded (FFPE) samples have been developed, although they are cost-intensive and time-consuming. An alternative to FFPE and snap-freezing is the PAXgene Tissue System, developed for simultaneous preservation of morphology, proteins, and nucleic acids. In the current study, we compared the performance of DNA from either PAXgene or formalin-fixed tissues to snap-frozen material for genome-wide DNA methylation analysis using the Illumina 450K BeadChip. Quantitative DNA methylation analysis demonstrated that the methylation profile in PAXgene-fixed tissues showed, in comparison with restored FFPE samples, a higher concordance with the profile detected in frozen samples. We demonstrate, for the first time, that DNA from PAXgene conserved tissue performs better compared with restored FFPE DNA in genome-wide DNA methylation analysis. In addition, DNA from PAXgene tissue can be directly used on the array without prior restoration, rendering the analytical process significantly more time- and cost-effective.     

Global DNA methylation in a population with aflatoxin B1 exposure

We previously reported that global DNA hypomethylation, measured as Sat2 methylation in white blood cells (WBC), and aflatoxin B₁ (AFB₁) exposure were associated with increased hepatocellular carcinoma risk. In this study, we assessed the association between AFB₁ exposure and global DNA methylation. We measured LINE-1 and Sat2 methylation in WBC DNA samples from 1140 cancer free participants of the Cancer Screening Program (CSP) cohort. Blood and urine samples were used to determine the level of AFB₁-albumin (AFB₁-Alb) adducts and urinary AFB₁ metabolites. In continuous models, we found reverse associations of urinary AFB₁ with LINE-1 and Sat2 methylation. The odds ratio (OR) per 1 unit decrease were 1.12 (95%CI = 1.03–1.22) for LINE-1 and 1.48 (95%CI = 1.10–2.00) for Sat2 methylation. When compared with subjects in the highest quartile of LINE-1, we found that individuals in the 2nd and 3rd quartiles were less likely to have detectable AFB₁-Alb adducts, with ORs (95%CI) of 0.61 (0.40–0.93), 0.61 (0.40-.94), and 1.09 (0.69–1.72), respectively. The OR for detectable AFB₁-Alb was 1.81 (95%CI = 1.15–2.85) for subjects in the lowest quartile of Sat2 methylation. The OR for detection of urinary AFB₁ for those with LINE-1 methylation in the lowest quartile compared with those in the highest quartile was 1.87 (95%CI = 1.15–3.04). The corresponding OR was 1.75 (95%CI = 1.08–2.82) for subjects in the lowest quartile of Sat2 methylation. The association between AFB₁ exposure and global DNA methylation may have implications for the epigenetic effect of AFB₁ on hepatocellular carcinoma development and also suggests that changes in DNA methylation may represent an epigenetic biomarker of dietary AFB₁ exposure.     

Identification of a new locus and validation of previously reported loci showing differential methylation associated with smoking. The REGICOR study

Smoking increases the risk of many diseases and could act through changes in DNA methylation patterns. The aims of this study were to determine the association between smoking and DNA methylation throughout the genome at cytosine-phosphate-guanine (CpG) site level and genomic regions. A discovery cross-sectional epigenome-wide association study nested in the follow-up of the REGICOR cohort was designed and included 645 individuals. Blood DNA methylation was assessed using the Illumina HumanMethylation450 BeadChip. Smoking status was self-reported using a standardized questionnaire. We identified 66 differentially methylated CpG sites associated with smoking, located in 38 genes. In most of these CpG sites, we observed a trend among those quitting smoking to recover methylation levels typical of never smokers. A CpG site located in a novel smoking-associated gene (cg06394460 in LNX2) was hypomethylated in current smokers. Moreover, we validated two previously reported CpG sites (cg05886626 in THBS1, and cg24838345 in MTSS1) for their potential relation to atherosclerosis and cancer diseases, using several different approaches: CpG site methylation, gene expression, and plasma protein level determinations. Smoking was also associated with higher THBS1 gene expression but with lower levels of thrombospondin-1 in plasma. Finally, we identified differential methylation regions in 13 genes and in four non-coding RNAs. In summary, this study replicated previous findings and identified and validated a new CpG site located in LNX2 associated with smoking.     

Quality control and statistical modeling for environmental epigenetics: A study on in utero lead exposure and DNA methylation at birth

DNA methylation data assayed using pyrosequencing techniques are increasingly being used in human cohort studies to investigate associations between epigenetic modifications at candidate genes and exposures to environmental toxicants and to examine environmentally-induced epigenetic alterations as a mechanism underlying observed toxicant-health outcome associations. For instance, in utero lead (Pb) exposure is a neurodevelopmental toxicant of global concern that has also been linked to altered growth in human epidemiological cohorts; a potential mechanism of this association is through alteration of DNA methylation (e.g., at growth-related genes). However, because the associations between toxicants and DNA methylation might be weak, using appropriate quality control and statistical methods is important to increase reliability and power of such studies. Using a simulation study, we compared potential approaches to estimate toxicant-DNA methylation associations that varied by how methylation data were analyzed (repeated measures vs. averaging all CpG sites) and by method to adjust for batch effects (batch controls vs. random effects). We demonstrate that correcting for batch effects using plate controls yields unbiased associations, and that explicitly modeling the CpG site-specific variances and correlations among CpG sites increases statistical power. Using the recommended approaches, we examined the association between DNA methylation (in LINE-1 and growth related genes IGF2, H19 and HSD11B2) and 3 biomarkers of Pb exposure (Pb concentrations in umbilical cord blood, maternal tibia, and maternal patella), among mother-infant pairs of the Early Life Exposures in Mexico to Environmental Toxicants (ELEMENT) cohort (n = 247). Those with 10 μg/g higher patella Pb had, on average, 0.61% higher IGF2 methylation (P = 0.05). Sex-specific trends between Pb and DNA methylation (P < 0.1) were observed among girls including a 0.23% increase in HSD11B2 methylation with 10 μg/g higher patella Pb.     

Quality control and statistical modeling for environmental epigenetics: A study on in utero lead exposure and DNA methylation at birth

DNA methylation data assayed using pyrosequencing techniques are increasingly being used in human cohort studies to investigate associations between epigenetic modifications at candidate genes and exposures to environmental toxicants and to examine environmentally-induced epigenetic alterations as a mechanism underlying observed toxicant-health outcome associations. For instance, in utero lead (Pb) exposure is a neurodevelopmental toxicant of global concern that has also been linked to altered growth in human epidemiological cohorts; a potential mechanism of this association is through alteration of DNA methylation (e.g., at growth-related genes). However, because the associations between toxicants and DNA methylation might be weak, using appropriate quality control and statistical methods is important to increase reliability and power of such studies. Using a simulation study, we compared potential approaches to estimate toxicant-DNA methylation associations that varied by how methylation data were analyzed (repeated measures vs. averaging all CpG sites) and by method to adjust for batch effects (batch controls vs. random effects). We demonstrate that correcting for batch effects using plate controls yields unbiased associations, and that explicitly modeling the CpG site-specific variances and correlations among CpG sites increases statistical power. Using the recommended approaches, we examined the association between DNA methylation (in LINE-1 and growth related genes IGF2, H19 and HSD11B2) and 3 biomarkers of Pb exposure (Pb concentrations in umbilical cord blood, maternal tibia, and maternal patella), among mother-infant pairs of the Early Life Exposures in Mexico to Environmental Toxicants (ELEMENT) cohort (n = 247). Those with 10 μg/g higher patella Pb had, on average, 0.61% higher IGF2 methylation (P = 0.05). Sex-specific trends between Pb and DNA methylation (P < 0.1) were observed among girls including a 0.23% increase in HSD11B2 methylation with 10 μg/g higher patella Pb.     

Correlation of DNA methylation levels in blood and saliva DNA in young girls of the LEGACY Girls study

Many epidemiologic studies of environmental exposures and disease susceptibility measure DNA methylation in white blood cells (WBC). Some studies are also starting to use saliva DNA as it is usually more readily available in large epidemiologic studies. However, little is known about the correlation of methylation between WBC and saliva DNA. We examined DNA methylation in three repetitive elements, Sat2, Alu, and LINE-1, and in four CpG sites, including AHRR (cg23576855, cg05575921), cg05951221 at 2q37.1, and cg11924019 at CYP1A1, in 57 girls aged 6–15 years with blood and saliva collected on the same day. We measured all DNA methylation markers by bisulfite-pyrosequencing, except for Sat2 and Alu, which were measured by the MethyLight assay. Methylation levels measured in saliva DNA were lower than those in WBC DNA, with differences ranging from 2.8% for Alu to 14.1% for cg05575921. Methylation levels for the three repetitive elements measured in saliva DNA were all positively correlated with those in WBC DNA. However, there was a wide range in the Spearman correlations, with the smallest correlation found for Alu (0.24) and the strongest correlation found for LINE-1 (0.73). Spearman correlations for cg05575921, cg05951221, and cg11924019 were 0.33, 0.42, and 0.79, respectively. If these findings are replicated in larger studies, they suggest that, for selected methylation markers (e.g., LINE-1), methylation levels may be highly correlated between blood and saliva, while for others methylation markers, the levels may be more tissue specific. Thus, in studies that differ by DNA source, each interrogated site should be separately examined in order to evaluate the correlation in DNA methylation levels across DNA sources.     

Cooperative activity of DNA methyltransferases for maintenance of symmetrical and non-symmetrical cytosine methylation in Arabidopsis thaliana

Maintenance of cytosine methylation in plants is controlled by three DNA methyltransferases. MET1 maintains CG methylation, and DRM1/2 and CMT3 act redundantly to enforce non-CG methylation. RPS, a repetitive hypermethylated DNA fragment from Petunia hybrida, attracts DNA methylation when transferred into Petunia or other species. In Arabidopsis thaliana, which does not contain any RPS homologues, RPS transgenes are efficiently methylated in all sequence contexts. To test which DNA methylation pathways regulate RPS methylation, we examined maintenance of RPS methylation in various mutant backgrounds. Surprisingly, CG methylation was lost in a drm1/2/cmt3 mutant, and non-CG methylation was almost completely eliminated in a met1 mutant. An unusual cooperative activity of all three DNA methyltransferases is therefore required for maintenance of both CG and non-CG methylation in RPS. Other unusual features of RPS methylation are the independence of its non-CG methylation from the RNA-directed DNA methylation (RdDM) pathway and the exceptional maintenance of methylation at a CC(m)TGG site in some epigenetic mutants. This is indicative of activity of a methylation system in plants that may have evolved from the DCM methylation system that controls CC(m)WGG methylation in bacteria. Our data suggest that strict separation of CG and non-CG methylation pathways does not apply to all target regions, and that caution is required in generalizing methylation data obtained for individual genomic regions.     

Altering TET dioxygenase levels within physiological range affects DNA methylation dynamics of HEK293 cells

The TET family of dioxygenases (TET1/2/3) can convert 5-methylcytosine (5mC) into 5-hydroxymethylcytosine (5hmC) and has been shown to be involved in active and passive DNA demethylation. Here, we demonstrate that altering TET dioxygenase levels within physiological range can affect DNA methylation dynamics of HEK293 cells. Overexpression of TET1 increased global 5hmC levels and was accompanied by mild DNA demethylation of promoters, gene bodies and CpG islands. Conversely, the simultaneous knockdown of TET1, TET2, and TET3 led to decreased global 5hmC levels and mild DNA hypermethylation of above-mentioned regions. The methylation changes observed in the overexpression and knockdown studies were mostly non-reciprocal and occurred with different preference depending on endogenous methylation and gene expression levels. Single-nucleotide 5hmC profiling performed on a genome-wide scale revealed that TET1 overexpression induced 5mC oxidation without a distribution bias among genetic elements and structures. Detailed analysis showed that this oxidation was related to endogenous 5hmC levels. In addition, our results support the notion that the effects of TET1 overexpression on gene expression are generally unrelated to its catalytic activity.     

Correlation of DNA methylation levels in blood and saliva DNA in young girls of the LEGACY Girls study

Many epidemiologic studies of environmental exposures and disease susceptibility measure DNA methylation in white blood cells (WBC). Some studies are also starting to use saliva DNA as it is usually more readily available in large epidemiologic studies. However, little is known about the correlation of methylation between WBC and saliva DNA. We examined DNA methylation in three repetitive elements, Sat2, Alu, and LINE-1, and in four CpG sites, including AHRR (cg23576855, cg05575921), cg05951221 at 2q37.1, and cg11924019 at CYP1A1, in 57 girls aged 6–15 years with blood and saliva collected on the same day. We measured all DNA methylation markers by bisulfite-pyrosequencing, except for Sat2 and Alu, which were measured by the MethyLight assay. Methylation levels measured in saliva DNA were lower than those in WBC DNA, with differences ranging from 2.8% for Alu to 14.1% for cg05575921. Methylation levels for the three repetitive elements measured in saliva DNA were all positively correlated with those in WBC DNA. However, there was a wide range in the Spearman correlations, with the smallest correlation found for Alu (0.24) and the strongest correlation found for LINE-1 (0.73). Spearman correlations for cg05575921, cg05951221, and cg11924019 were 0.33, 0.42, and 0.79, respectively. If these findings are replicated in larger studies, they suggest that, for selected methylation markers (e.g., LINE-1), methylation levels may be highly correlated between blood and saliva, while for others methylation markers, the levels may be more tissue specific. Thus, in studies that differ by DNA source, each interrogated site should be separately examined in order to evaluate the correlation in DNA methylation levels across DNA sources.