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.     

BREX is a novel phage resistance system widespread in microbial genomes

The perpetual arms race between bacteria and phage has resulted in the evolution of efficient resistance systems that protect bacteria from phage infection. Such systems, which include the CRISPR‐Cas and restriction‐modification systems, have proven to be invaluable in the biotechnology and dairy industries. Here, we report on a six‐gene cassette in Bacillus cereus which, when integrated into the Bacillus subtilis genome, confers resistance to a broad range of phages, including both virulent and temperate ones. This cassette includes a putative Lon‐like protease, an alkaline phosphatase domain protein, a putative RNA‐binding protein, a DNA methylase, an ATPase‐domain protein, and a protein of unknown function. We denote this novel defense system BREX (Bacteriophage Exclusion) and show that it allows phage adsorption but blocks phage DNA replication. Furthermore, our results suggest that methylation on non‐palindromic TAGGAG motifs in the bacterial genome guides self/non‐self discrimination and is essential for the defensive function of the BREX system. However, unlike restriction‐modification systems, phage DNA does not appear to be cleaved or degraded by BREX, suggesting a novel mechanism of defense. Pan genomic analysis revealed that BREX and BREX‐like systems, including the distantly related Pgl system described in Streptomyces coelicolor, are widely distributed in ~10% of all sequenced microbial genomes and can be divided into six coherent subtypes in which the gene composition and order is conserved. Finally, we detected a phage family that evades the BREX defense, implying that anti‐BREX mechanisms may have evolved in some phages as part of their arms race with bacteria.     

DNA hypermethylation and DNA hypomethylation is present at different loci in chronic kidney disease

Genetic risk factors for chronic kidney disease (CKD) are being identified through international collaborations. By comparison, epigenetic risk factors for CKD have only recently been considered using population-based approaches. DNA methylation is a major epigenetic modification that is associated with complex diseases, so we investigated methylome-wide loci for association with CKD. A total of 485,577 unique features were evaluated in 255 individuals with CKD (cases) and 152 individuals without evidence of renal disease (controls). Following stringent quality control, raw data were quantile normalized and β values calculated to reflect the methylation status at each site. The difference in methylation status was evaluated between cases and controls with resultant P values adjusted for multiple testing. Genes with significantly increased and decreased levels of DNA methylation were considered for biological relevance by functional enrichment analysis using KEGG pathways in Partek Genomics Suite. Twenty-three genes, where more than one CpG per loci was identified with P_adjusted < 10⁻⁸, demonstrated significant methylation changes associated with CKD and additional support for these associated loci was sought from published literature. Strong biological candidates for CKD that showed statistically significant differential methylation include CUX1, ELMO1, FKBP5, INHBA-AS1, PTPRN2, and PRKAG2 genes; several genes are differentially methylated in kidney tissue and RNA-seq supports a functional role for differential methylation in ELMO1 and PRKAG2 genes. This study reports the largest, most comprehensive, genome-wide quantitative evaluation of DNA methylation for association with CKD. Evidence confirming methylation sites influence development of CKD would stimulate research to identify epigenetic therapies that might be clinically useful for CKD.     

Characterization of SpPol4, a unique X-family DNA polymerase in Schizosaccharomyces pombe

As predicted by the amino acid sequence, the purified protein coded by Schizosaccharomyces pombe SPAC2F7.06c is a DNA polymerase (SpPol4) whose biochemical properties resemble those of other X family (PolX) members. Thus, this new PolX is template-dependent, polymerizes in a distributive manner, lacks a detectable 3′→5′ proofreading activity and its preferred substrates are small gaps with a 5′-phosphate group. Similarly to Polμ, SpPol4 can incorporate a ribonucleotide (rNTP) into a primer DNA. However, it is not responsible for the 1–2 rNTPs proposed to be present at the mating-type locus and those necessary for mating-type switching. Unlike Polμ, SpPol4 lacks terminal deoxynucleotidyltransferase activity and realigns the primer terminus to alternative template bases only under certain sequence contexts and, therefore, it is less error-prone than Polμ. Nonetheless, the biochemical properties of this gap-filling DNA polymerase are suitable for a possible role of SpPol4 in non-homologous end-joining. Unexpectedly based on sequence analysis, SpPol4 has deoxyribose phosphate lyase activity like Polβ and Polλ, and unlike Polμ, suggesting also a role of this enzyme in base excision repair. Therefore, SpPol4 is a unique enzyme whose enzymatic properties are hybrid of those described for mammalian Polβ, Polλ and Polμ.     

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.     

Berberine reduces methylation of the MTTP promoter and alleviates fatty liver induced by a high-fat diet in rats

High-calorie food leads to nonalcoholic fatty liver disease (NAFLD) through dysregulation of genes involved in lipid metabolism, but the precise mechanism remains unclear. DNA methylation represents one of the mechanisms that contributes to dysregulation of gene expression via interaction with environmental factors. Berberine can alleviate fatty liver in db/db and ob/ob mice. Here, we investigated whether DNA methylation is involved in the pathogenesis of NAFLD induced by a high-fat diet (HFD) and whether berberine improves NAFLD through influencing the methylation status of promoters of key genes. HFD markedly decreased the mRNA levels encoding CPT-1α, MTTP, and LDLR in the liver. In parallel, DNA methylation levels in the MTTP promoter of rats with NAFLD were elevated in the liver. Interestingly, berberine reversed the downregulated expression of these genes and selectively inhibited HFD-induced increase in the methylation of MTTP. Consistently, berberine increased hepatic triglyceride (TG) export and ameliorated HFD-induced fatty liver. Furthermore, a close negative correlation was observed between the MTTP expression and its DNA methylation (at sites −113 and −20). These data indicate that DNA methylation of the MTTP promoter likely contributes to its downregulation during HFD-induced NAFLD and, further, that berberine can partially counteract the HFD-elicited dysregulation of MTTP by reversing the methylation state of its promoter, leading to reduced hepatic fat content.     

Whole-Genome Saliva and Blood DNA Methylation Profiling in Individuals with a Respiratory Allergy

The etiology of respiratory allergies (RA) can be partly explained by DNA methylation changes caused by adverse environmental and lifestyle factors experienced early in life. Longitudinal, prospective studies can aid in the unravelment of the epigenetic mechanisms involved in the disease development. High compliance rates can be expected in these studies when data is collected using non-invasive and convenient procedures. Saliva is an attractive biofluid to analyze changes in DNA methylation patterns. We investigated in a pilot study the differential methylation in saliva of RA (n = 5) compared to healthy controls (n = 5) using the Illumina Methylation 450K BeadChip platform. We evaluated the results against the results obtained in mononuclear blood cells from the same individuals. Differences in methylation patterns from saliva and mononuclear blood cells were clearly distinguishable (P_Adj<0.001 and |Δβ|>0.2), though the methylation status of about 96% of the cg-sites was comparable between peripheral blood mononuclear cells and saliva. When comparing RA cases with healthy controls, the number of differentially methylated sites (DMS) in saliva and blood were 485 and 437 (P<0.05 and |Δβ|>0.1), respectively, of which 216 were in common. The methylation levels of these sites were significantly correlated between blood and saliva. The absolute levels of methylation in blood and saliva were confirmed for 3 selected DMS in the PM20D1, STK32C, and FGFR2 genes using pyrosequencing analysis. The differential methylation could only be confirmed for DMS in PM20D1 and STK32C genes in saliva. We show that saliva can be used for genome-wide methylation analysis and that it is possible to identify DMS when comparing RA cases and healthy controls. The results were replicated in blood cells of the same individuals and confirmed by pyrosequencing analysis. This study provides proof-of-concept for the applicability of saliva-based whole-genome methylation analysis in the field of respiratory allergy.     

Efficient Transformation System for Propionibacterium freudenreichii Based on a Novel Vector

A 3.6-kb endogenous plasmid was isolated from a Propionibacterium freudenreichii strain and sequenced completely. Based on homologies with plasmids from other bacteria, notably a plasmid from Mycobacterium, a region harboring putative replicative functions was defined. Outside this region two restriction enzyme recognition sites were used for insertion of an Escherichia coli-specific replicon and an erythromycin resistance gene for selection in Propionibacterium. Hybrid vectors obtained in this way replicated in both E. coli and P. freudenreichii. Whereas electroporation of P. freudenreichii with vector DNA isolated from an E. coli transformant yielded 10 to 30 colonies per μg of DNA, use of vector DNA reisolated from a Propionibacterium transformant dramatically increased the efficiency of transformation (≥10⁸ colonies per μg of DNA). It could be shown that restriction-modification was responsible for this effect. The high efficiency of the system described here permitted successful transformation of Propionibacterium with DNA ligation mixtures.     

Selection against PUMA gene expression in Myc-driven B-cell lymphomagenesis

The p53 tumor suppressor pathway limits oncogenesis by inducing cell cycle arrest or apoptosis. A key p53 target gene is PUMA, which encodes a BH3-only proapoptotic protein. Here we demonstrate that Puma deletion in the Eμ-Myc mouse model of Burkitt lymphoma accelerates lymphomagenesis and that ~75% of Eμ-Myc lymphomas naturally select against Puma protein expression. Furthermore, approximately 40% of primary human Burkitt lymphomas fail to express detectable levels of PUMA and in some tumors this is associated with DNA methylation. Burkitt lymphoma cell lines phenocopy the primary tumors with respect to DNA methylation and diminished PUMA expression, which can be reactivated following inhibition of DNA methyltransferases. These findings establish that PUMA is silenced in human malignancies, and they suggest PUMA as a target for the development of novel chemotherapeutics.     

Higher Urinary Bisphenol A Concentration Is Associated with Unexplained Recurrent Miscarriage Risk: Evidence from a Case-Control Study in Eastern China

BackgroundEvidence about the association between Bisphenol A (BPA) and the risk of recurrent miscarriage (RM) in human being is still limited.ObjectiveWe evaluated the association of urinary BPA concentrations with RM in human being.MethodsA hospital-based 1:2 matched case-control study on RM was carried out in Suzhou and Kunshan in Jiangsu Province in China between August 2008 and November 2011. Total urinary BPA concentrations in 264 eligible urine samples (102 RM patients and 162 controls) were measured using liquid chromatography-tandem mass spectrometry (LC-MS/MS). The Wilcoxon test and conditional logistic regression were used to estimate the differences between the groups and odds ratios (OR) with 95% confidence intervals (CI), respectively.ResultsThe median ± IQR (interquartile range) (P₇₅-P₂₅) values of non-creatinine-adjusted total urinary BPA levels in the RM patients and the controls were 1.66±3.69ng/ml and 0.58±1.07ng/ml, respectively (0.98±2.67μg/g Cr (creatinine) and 0.40±0.77μg/g Cr. The adjusted BPA level was significantly higher in the RM patients than in the controls (Wilcoxon test, Z = 4.476, P<0.001). Higher level of urinary BPA was significantly associated with an increased risk of RM (P-trend <0.001). Compared to the groups with urinary BPA levels less than 0.16μg/g Cr, the women with levels of 0.40–0.93μg/g Cr and 0.93μg/g Cr or above had a significantly higher risk of RM (OR = 3.91, 95%CI: 1.23–12.45 and OR = 9.34, 95%CI: 3.06–28.44) that persisted after adjusting for confounding factors. The time from recently RM date to recruitment does not significantly influence the urinary BPA level (P = 0.090).ConclusionExposure to BPA may be associated with RM risk.     

TGF-β induces global changes in DNA methylation during the epithelial-to-mesenchymal transition in ovarian cancer cells

A key step in the process of metastasis is the epithelial-to-mesenchymal transition (EMT). We hypothesized that epigenetic mechanisms play a key role in EMT and to test this hypothesis we analyzed global and gene-specific changes in DNA methylation during TGF-β-induced EMT in ovarian cancer cells. Epigenetic profiling using the Infinium HumanMethylation450 BeadChip (HM450) revealed extensive (P < 0.01) methylation changes after TGF-β stimulation (468 and 390 CpG sites altered at 48 and 120 h post cytokine treatment, respectively). The majority of gene-specific TGF-β-induced methylation changes occurred in CpG islands located in or near promoters (193 and 494 genes hypermethylated at 48 and 120 h after TGF-β stimulation, respectively). Furthermore, methylation changes were sustained for the duration of TGF-β treatment and reversible after the cytokine removal. Pathway analysis of the hypermethylated loci identified functional networks strongly associated with EMT and cancer progression, including cellular movement, cell cycle, organ morphology, cellular development, and cell death and survival. Altered methylation and corresponding expression of specific genes during TGF-β-induced EMT included CDH1 (E-cadherin) and COL1A1 (collagen 1A1). Furthermore, TGF-β induced both expression and activity of DNA methyltransferases (DNMT) -1, -3A, and -3B, and treatment with the DNMT inhibitor SGI-110 prevented TGF-β-induced EMT. These results demonstrate that dynamic changes in the DNA methylome are implicated in TGF-β-induced EMT and metastasis. We suggest that targeting DNMTs may inhibit this process by reversing the EMT genes silenced by DNA methylation in cancer.     

Esp1396I restriction-modification system: structural organization and mode of regulation

Esp1396I restriction–modification (RM) system recognizes an interrupted palindromic DNA sequ ence 5′-CCA(N)₅TGG-3′. The Esp1396I RM system was found to reside on pEsp1396, a 5.6 kb plasmid naturally occurring in Enterobacter sp. strain RFL1396. The nucleotide sequence of the entire 5622 bp pEsp1396 plasmid was determined on both strands. Identified genes for DNA methyltransferase (esp1396IM) and restriction endonuclease (esp1396IR) are transcribed convergently. The restriction endonuclease gene is preceded by the small ORF (esp1396IC) that possesses a strong helix-turn-helix motif and resembles regulatory proteins found in PvuII, BamHI and few other RM systems. Gene regulation studies revealed that C.Esp1396I acts as both a repressor of methylase expression and an activator of regulatory protein and restriction endonuclease expression. Our data indicate that C protein from Esp1396I RM system activates the expression of the Enase gene, which is co-transcribed from the promoter of regulatory gene, by the mechanism of coupled translation.     

Chromosomal Gene Inactivation in the Green Sulfur Bacterium Chlorobium tepidum by Natural Transformation

Conditions for inactivating chromosomal genes of Chlorobium tepidum by natural transformation and homologous recombination were established. As a model, mutants unable to perform nitrogen fixation were constructed by interrupting nifD with various antibiotic resistance markers. Growth of wild-type C. tepidum at 40°C on agar plates could be completely inhibited by 100 μg of gentamicin ml⁻¹, 2 μg of erythromycin ml⁻¹, 30 μg of chloramphenicol ml⁻¹, or 1 μg of tetracycline ml⁻¹ or a combination of 300 μg of streptomycin ml⁻¹ and 150 μg of spectinomycin ml⁻¹. Transformation was performed by spotting cells and DNA on an agar plate for 10 to 20 h. Transformation frequencies on the order of 10⁻⁷ were observed with gentamicin and erythromycin markers, and transformation frequencies on the order of 10⁻³ were observed with a streptomycin-spectinomycin marker. The frequency of spontaneous mutants resistant to gentamicin, erythromycin, or spectinomycin-streptomycin was undetectable or significantly lower than the transformation frequency. Transformation with the gentamicin marker was observed when the transforming DNA contained 1 or 3 kb of total homologous flanking sequence but not when the transforming DNA contained only 0.3 kb of homologous sequence. Linearized plasmids transformed at least an order of magnitude better than circular plasmids. This work forms a foundation for the systematic targeted inactivation of genes in C. tepidum, whose 2.15-Mb genome has recently been completely sequenced.     

Isolation of a Multiheme Protein with Features of a Hydrazine-Oxidizing Enzyme from an Anaerobic Ammonium-Oxidizing Enrichment Culture

A multiheme protein having hydrazine-oxidizing activity was purified from enriched culture from a reactor in which an anammox bacterium, strain KSU-1, was dominant. The enzyme has oxidizing activity toward hydrazine but not hydroxylamine and is a 130-kDa homodimer composed of a 62-kDa polypeptide containing eight hemes. It was therefore named hydrazine-oxidizing enzyme (HZO). With cytochrome c as an electron acceptor, the V_max and K_m for hydrazine are 6.2 ± 0.3 μmol/min · mg and 5.5 ± 0.6 μM, respectively. Hydrazine (25 μM) induced an increase in the proportion of reduced form in the spectrum, whereas hydroxylamine (500 μM) did not. Two genes coding for HZO, hzoA and hzoB, were identified within the metagenomic DNA from the culture. The genes encode the same amino acid sequence except for two residues. The sequences deduced from these genes showed low-level identities (<30%) to those of all of the hydroxylamine oxidoreductases reported but are highly homologous to two hao genes found by sequencing the genome of “Candidatus Kuenenia stuttgartiensis” (88% and 89% identities). The purified enzyme might therefore be a novel hydrazine-oxidizing enzyme having a critical role in anaerobic ammonium oxidation.     

Genome-wide differentially methylated genes in prostate cancer tissues from African-American and Caucasian men

Increasing evidence suggests that aberrant DNA methylation changes may contribute to prostate cancer (PCa) ethnic disparity. To comprehensively identify DNA methylation alterations in PCa disparity, we used the Illumina 450K methylation platform to interrogate the methylation status of 485,577 CpG sites focusing on gene-associated regions of the human genome. Genomic DNA from African-American (AA; 7 normal and 3 cancers) and Caucasian (Cau; 8 normal and 3 cancers) was used in the analysis. Hierarchical clustering analysis identified probe-sets unique to AA and Cau samples, as well as common to both. We selected 25 promoter-associated novel CpG sites most differentially methylated by race (fold change > 1.5-fold; adjusted P < 0.05) and compared the β-value of these sites provided by the Illumina, Inc. array with quantitative methylation obtained by pyrosequencing in 7 prostate cell lines. We found very good concordance of the methylation levels between β-value and pyrosequencing. Gene expression analysis using qRT-PCR in a subset of 8 genes after treatment with 5-aza-2′-deoxycytidine and/or trichostatin showed up-regulation of gene expression in PCa cells. Quantitative analysis of 4 genes, SNRPN, SHANK2, MST1R, and ABCG5, in matched normal and PCa tissues derived from AA and Cau PCa patients demonstrated differential promoter methylation and concomitant differences in mRNA expression in prostate tissues from AA vs. Cau. Regression analysis in normal and PCa tissues as a function of race showed significantly higher methylation prevalence for SNRPN (P = 0.012), MST1R (P = 0.038), and ABCG5 (P < 0.0002) for AA vs. Cau samples. We selected the ABCG5 and SNRPN genes and verified their biological functions by Western blot analysis and siRNA gene knockout effects on cell proliferation and invasion in 4 PCa cell lines (2 AA and 2 Cau patients-derived lines). Knockdown of either ABCG5 or SNRPN resulted in a significant decrease in both invasion and proliferation in Cau PCa cell lines but we did not observe these remarkable loss-of-function effects in AA PCa cell lines. Our study demonstrates how differential genome-wide DNA methylation levels influence gene expression and biological functions in AA and Cau PCa.