Aberrant Methylation of Gene Associated CpG Sites Occurs in Borderline Personality Disorder

Borderline personality disorder (BPD) is a complex psychiatric disease with an increased impact in the last years. While the diagnosis and therapy are well established, little is known on the pathogenesis of borderline personality disorder. Previously, a significant increase in DNA methylation of relevant neuropsychiatric genes in BPD patients has been reported. In our study we performed genome wide methylation analysis and revealed specific CpG sites that exhibited increased methylation in 24 female BPD patients compared to 11 female healthy controls. Bead chip technology and quantitative bisulfite pyrosequencing showed a significantly increased methylation at CpG sites of APBA2 (1.1 fold) and APBA3 (1.1 fold), KCNQ1 (1.5 fold), MCF2 (1.1 fold) and NINJ2 (1.2 fold) in BPD patients. For the CpG sites of GATA4 and HLCS an increase in DNA methylation was observed, but was only significant in the bead chip assay. Moreover genome wide methylation levels of blood samples of BPD patients and control samples are similar. In summary, our results show a significant 1.26 fold average increase in methylation at the analyzed gene associated CpG sites in the blood of BPD patients compared to controls samples (p<0.001). This data may provide new insights into epigenetic mechanisms underlying the pathogenesis of BPD.     

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.     

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.     

Prenatal antidepressant exposure associated with CYP2E1 DNA methylation change in neonates

Some but not all neonates are affected by prenatal exposure to serotonin reuptake inhibitor antidepressants (SRI) and maternal mood disturbances. Distinguishing the impact of these 2 exposures is challenging and raises critical questions about whether pharmacological, genetic, or epigenetic factors can explain the spectrum of reported outcomes. Using unbiased DNA methylation array measurements followed by a detailed candidate gene approach, we examined whether prenatal SRI exposure was associated with neonatal DNA methylation changes and whether such changes were associated with differences in birth outcomes. Prenatal SRI exposure was first associated with increased DNA methylation status primarily at CYP2E1(β_Non-exposed = 0.06, β_SRI-exposed = 0.30, FDR = 0); however, this finding could not be distinguished from the potential impact of prenatal maternal depressed mood. Then, using pyrosequencing of CYP2E1 regulatory regions in an expanded cohort, higher DNA methylation status—both the mean across 16 CpG sites (P < 0.01) and at each specific CpG site (P < 0.05)—was associated with exposure to lower 3rd trimester maternal depressed mood symptoms only in the SRI-exposed neonates, indicating a maternal mood x SRI exposure interaction. In addition, higher DNA methylation levels at CpG2 (P = 0.04), CpG9 (P = 0.04) and CpG10 (P = 0.02), in the interrogated CYP2E1 region, were associated with increased birth weight independently of prenatal maternal mood, SRI drug exposure, or gestational age at birth. Prenatal SRI antidepressant exposure and maternal depressed mood were associated with altered neonatal CYP2E1 DNA methylation status, which, in turn, appeared to be associated with birth weight.     

DNA methylation in neonates born to women receiving psychiatric care

Prenatal exposure both to maternal psychiatric illness and psychiatric medication has been linked with adverse child outcomes that affect physiological, emotional and psychiatric development. Studies suggest that epigenetic mechanisms, such as DNA methylation, may facilitate these effects. In this report, we explore the association between maternal psychiatric illness and treatment during pregnancy and neonatal DNA methylation patterns in a prospectively-characterized clinical cohort of 201 dyads. Associations between the percent of umbilical cord blood DNA methylated at 27,578 CpG sites and maternal psychiatric diagnosis, symptoms and antidepressant use were evaluated by fitting a separate linear mixed effects model for each CpG site. There were no significant changes in neonatal DNA methylation attributable to maternal psychiatric diagnosis or depressive symptoms during pregnancy. Exposure to an antidepressant medication was associated with differential methylation of CpG sites in TNFRSF21 and CHRNA2 (false discovery rate < 0.05), but the average difference in methylation for both CpG sites was less than 3% between each group. The results were not specific to type of antidepressant or duration of the exposure. This study suggests that there are no large effects of maternal psychiatric illness, depressive symptoms or prenatal exposure to antidepressants on neonatal DNA methylation. Delineation of the influence of maternal psychiatric illness and pharmacological exposures on the developing fetuses has critical implications for clinical care during pregnancy.     

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.     

Polycomb group genes are targets of aberrant DNA methylation in renal cell carcinoma

The combined effects of genetic and epigenetic aberrations are well recognized as causal in tumorigenesis. Here, we defined profiles of DNA methylation in primary renal cell carcinomas (RCC) and assessed the association of these profiles with the expression of genes required for the establishment and maintenance of epigenetic marks. A bead-based methylation array platform was used to measure methylation of 1,413 CpG loci in ∼800 cancer-associated genes and three methylation classes were derived by unsupervised clustering of tumors using recursively partitioned mixture modeling (RPMM). Quantitative RT-PCR was performed on all tumor samples to determine the expression of DNMT1, DNMT3B, VEZF1 and EZH2. Additionally, methylation at LINE-1 and AluYb8 repetitive elements was measured using bisulfite pyrosequencing. Associations between methylation class and tumor stage (p = 0.05), LINE-1 (p < 0.0001) and AluYb8 (p < 0.0001) methylation, as well as EZH2 expression (p < 0.0001) were noted following univariate analyses. A multinomial logistic regression model controlling for potential confounders revealed that AluYb8 (p < 0.003) methylation and EZH2 expression (p < 0.008) were significantly associated with methylation class membership. Because EZH2 is a member of the Polycomb repressive complex 2 (PRC2), we next analyzed the distribution of Polycomb group (PcG) targets among methylation classes derived by clustering the 1,413 array CpG loci using RPMM. PcG target genes were significantly enriched (p < 0.0001) in methylation classes with greater differential methylation between RCC and non-diseased kidney tissue. This work contributes to our understanding of how repressive marks on DNA and chromatin are dysregulated in carcinogenesis, knowledge that might aid the development of therapies or preventive strategies for human malignancies.Key words: EZH2, DNA methylation, renal cell carcinoma, polycomb, microarray     

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.     

Heterogeneity in White Blood Cells Has Potential to Confound DNA Methylation Measurements

Epigenetic studies are commonly conducted on DNA from tissue samples. However, tissues are ensembles of cells that may each have their own epigenetic profile, and therefore inter-individual cellular heterogeneity may compromise these studies. Here, we explore the potential for such confounding on DNA methylation measurement outcomes when using DNA from whole blood. DNA methylation was measured using pyrosequencing-based methodology in whole blood (n = 50–179) and in two white blood cell fractions (n = 20), isolated using density gradient centrifugation, in four CGIs (CpG Islands) located in genes HHEX (10 CpG sites assayed), KCNJ11 (8 CpGs), KCNQ1 (4 CpGs) and PM20D1 (7 CpGs). Cellular heterogeneity (variation in proportional white blood cell counts of neutrophils, lymphocytes, monocytes, eosinophils and basophils, counted by an automated cell counter) explained up to 40% (p<0.0001) of the inter-individual variation in whole blood DNA methylation levels in the HHEX CGI, but not a significant proportion of the variation in the other three CGIs tested. DNA methylation levels in the two cell fractions, polymorphonuclear and mononuclear cells, differed significantly in the HHEX CGI; specifically the average absolute difference ranged between 3.4–15.7 percentage points per CpG site. In the other three CGIs tested, methylation levels in the two fractions did not differ significantly, and/or the difference was more moderate. In the examined CGIs, methylation levels were highly correlated between cell fractions. In summary, our analysis detects region-specific differential DNA methylation between white blood cell subtypes, which can confound the outcome of whole blood DNA methylation measurements. Finally, by demonstrating the high correlation between methylation levels in cell fractions, our results suggest a possibility to use a proportional number of a single white blood cell type to correct for this confounding effect in analyses.     

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.     

Identification of a DNA Methylome Profile of Esophageal Squamous Cell Carcinoma and Potential Plasma Epigenetic Biomarkers for Early Diagnosis

DNA methylation is a critical epigenetic mechanism involved in key cellular processes. Its deregulation has been linked to many human cancers including esophageal squamous cell carcinoma (ESCC). This study was designed to explore the whole methylation status of ESCC and to identify potential plasma biomarkers for early diagnosis. We used Infinium Methylation 450k array to analyze ESCC tissues (n = 4), paired normal surrounding tissues (n = 4) and normal mucosa from healthy individuals (n = 4), and combined these with gene expression data from the GEO database. One hundred and sixty eight genes had differentially methylated CpG sites in their promoter region and a gene expression pattern inverse to the direction of change in DNA methylation. These genes were involved in several cancer-related pathways. Three genes were validated in additional 42 ESCC tissues and paired normal surrounding tissues. The methylation frequency of EPB41L3, GPX3, and COL14A1 were higher in tumor tissues than in normal surrounding tissues (P<0.017). The higher methylation frequency of EPB41l3 was correlated with large tumor size (P = 0.044) and advanced pT tumor stage (P = 0.001). The higher methylation frequency of GPX3 and COL14A1 were correlated with advanced pN tumor stage (P = 0.001 and P<0.001). The methylation of EPB41L3, GPX3, and COL14A1 genes were only found in ESCC patients'' plasma, but not in normal individuals upon testing 42 ESCC patients and 50 healthy individuals. Diagnostic sensitivity was increased when methylation of any of the 3 genes were counted (64.3% sensitivity and 100% specificity). These differentially methylated genes in plasma may be used as biomarkers for early diagnosis of ESCC.     

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.     

Genome-wide methylation analysis identifies differentially methylated CpG loci associated with severe obesity in childhood

Childhood obesity is a major public health issue. Here we investigated whether differential DNA methylation was associated with childhood obesity. We studied DNA methylation profiles in whole blood from 78 obese children (mean BMI Z-score: 2.6) and 71 age- and sex-matched controls (mean BMI Z-score: 0.1). DNA samples from obese and control groups were pooled and analyzed using the Infinium HumanMethylation450 BeadChip array. Comparison of the methylation profiles between obese and control subjects revealed 129 differentially methylated CpG (DMCpG) loci associated with 80 unique genes that had a greater than 10% difference in methylation (P-value < 0.05). The top pathways enriched among the DMCpGs included developmental processes, immune system regulation, regulation of cell signaling, and small GTPase-mediated signal transduction. The associations between the methylation of selected DMCpGs with childhood obesity were validated using sodium bisulfite pyrosequencing across loci within the FYN, PIWIL4, and TAOK3 genes in individual subjects. Three CpG loci within FYN were hypermethylated in obese individuals (all P < 0.01), while obesity was associated with lower methylation of CpG loci within PIWIL4 (P = 0.003) and TAOK3 (P = 0.001). After building logistic regression models, we determined that a 1% increase in methylation in TAOK3, multiplicatively decreased the odds of being obese by 0.91 (95% CI: 0.86 – 0.97), and an increase of 1% methylation in FYN CpG3, multiplicatively increased the odds of being obese by 1.03 (95% CI: 0.99 – 1.07). In conclusion, these findings provide evidence that childhood obesity is associated with specific DNA methylation changes in whole blood, which may have utility as biomarkers of obesity risk.