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.     

Meta-analysis of human methylation data for evidence of sex-specific autosomal patterns

Background

Several individual studies have suggested that autosomal CpG methylation differs by sex both in terms of individual CpG sites and global autosomal CpG methylation. However, these findings have been inconsistent and plagued by spurious associations due to the cross reactivity of CpG probes on commercial microarrays. We collectively analysed 76 published studies (n = 6,795) for sex-associated differences in both autosomal and sex chromosome CpG sites.

Results

Overall autosomal methylation profiles varied substantially by study, and we encountered substantial batch effects. We accounted for these by conducting random effects meta-analysis for individual autosomal CpG methylation associations. After excluding non-specific probes, we found 184 autosomal CpG sites differentially methylated by sex after correction for multiple testing. In line with previous studies, average beta differences were small. Many of the most significantly associated CpG probes were new. Of note was differential CpG methylation in the promoters of genes thought to be involved in spermatogenesis and male fertility, such as SLC9A2, SPESP1, CRISP2, and NUPL1. Pathway analysis revealed overrepresentation of genes differentially methylated by sex in several broad Gene Ontology biological processes, including RNA splicing and DNA repair.

Conclusions

This study represents a comprehensive analysis of sex-specific methylation patterns. We demonstrate the existence of sex-specific methylation profiles and report a large number of novel DNA methylation differences in autosomal CpG sites between sexes.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-981) contains supplementary material, which is available to authorized users.     

Methylated Host Cell Gene Promoters and Human Papillomavirus Type 16 and 18 Predicting Cervical Lesions and Cancer

Change in the host and/or human papillomavirus (HPV) DNA methylation profile is probably one of the main factors responsible for the malignant progression of cervical lesions to cancer. To investigate those changes we studied 173 cervical samples with different grades of cervical lesion, from normal to cervical cancer. The methylation status of nine cellular gene promoters, CCNA1, CDH1, C13ORF18, DAPK1, HIC1, RARβ2, hTERT1, hTERT2 and TWIST1, was investigated by Methylation Specific Polymerase Chain Reaction (MSP). The methylation of HPV18 L1-gene was also investigated by MSP, while the methylated cytosines within four regions, L1, 5’LCR, enhancer, and promoter of the HPV16 genome covering 19 CpG sites were evaluated by bisulfite sequencing. Statistically significant methylation biomarkers distinguishing between cervical precursor lesions from normal cervix were primarily C13ORF18 and secondly CCNA1, and those distinguishing cervical cancer from normal or cervical precursor lesions were CCNA1, C13ORF18, hTERT1, hTERT2 and TWIST1. In addition, the methylation analysis of individual CpG sites of the HPV16 genome in different sample groups, notably the 7455 and 7694 sites, proved to be more important than the overall methylation frequency. The majority of HPV18 positive samples contained both methylated and unmethylated L1 gene, and samples with L1-gene methylated forms alone had better prognosis when correlated with the host cell gene promoters’ methylation profiles. In conclusion, both cellular and viral methylation biomarkers should be used for monitoring cervical lesion progression to prevent invasive cervical cancer.     

Prognostic Value of PLAGL1-Specific CpG Site Methylation in Soft-Tissue Sarcomas

Soft tissue sarcomas (STS) are rare, complex tumors with a poor prognosis. The identification of new prognostic biomarkers is needed to improve patient management. Our aim was to determine the methylation status of the 118 CpG sites in the PLAGL1 tumor-suppressor gene P1 CpG island promoter and study the potential prognostic impact of PLAGL1 promoter methylation CpG sites in STS. Training cohorts constituted of 28 undifferentiated sarcomas (US) and 35 leiomyosarcomas (LMS) were studied. PLAGL1 mRNA expression was investigated by microarray analysis and validated by RT-qPCR. Pyrosequencing was used to analyze quantitative methylation of the PLAGL1 promoter. Associations between global promoter or specific CpG site methylation and mRNA expression were evaluated using Pearson’s product moment correlation coefficient. Cox univariate and multivariate proportional hazard models were used to assess the predictive power of CpG site methylation status. Sixteen CpG sites associated with PLAGL1 mRNA expression were identified in US and 6 in LMS. Statistical analyses revealed an association between CpG107 methylation status and both overall and metastasis-free survival in US, which was confirmed in a validation cohort of 37 US. The exhaustive study of P1 PLAGL1 promoter methylation identified a specific CpG site methylation correlated with mRNA expression, which was predictive for both metastasis-free and overall survival and may constitute the first US-specific biomarker. Such a biomarker may be relevant for identifying patients likely to derive greater benefit from treatment.     

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.     

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.     

DNA Methylation Profile Distinguishes Clear Cell Sarcoma of the Kidney from Other Pediatric Renal Tumors

A number of specific, distinct neoplastic entities occur in the pediatric kidney, including Wilms’ tumor, clear cell sarcoma of the kidney (CCSK), congenital mesoblastic nephroma (CMN), rhabdoid tumor of the kidney (RTK), and the Ewing’s sarcoma family of tumors (ESFT). By employing DNA methylation profiling using Illumina Infinium HumanMethylation27, we analyzed the epigenetic characteristics of the sarcomas including CCSK, RTK, and ESFT in comparison with those of the non-neoplastic kidney (NK), and these tumors exhibited distinct DNA methylation profiles in a tumor-type-specific manner. CCSK is the most frequently hypermethylated, but least frequently hypomethylated, at CpG sites among these sarcomas, and exhibited 490 hypermethylated and 46 hypomethylated CpG sites in compared with NK. We further validated the results by MassARRAY, and revealed that a combination of four genes was sufficient for the DNA methylation profile-based differentiation of these tumors by clustering analysis. Furthermore, THBS1 CpG sites were found to be specifically hypermethylated in CCSK and, thus, the DNA methylation status of these THBS1 sites alone was sufficient for the distinction of CCSK from other pediatric renal tumors, including Wilms’ tumor and CMN. Moreover, combined bisulfite restriction analysis could be applied for the detection of hypermethylation of a THBS1 CpG site. Besides the biological significance in the pathogenesis, the DNA methylation profile should be useful for the differential diagnosis of pediatric renal tumors.     

Genetic Control of Individual Differences in Gene-Specific Methylation in Human Brain

We have observed extensive interindividual differences in DNA methylation of 8590 CpG sites of 6229 genes in 153 human adult cerebellum samples, enriched in CpG island “shores” and at further distances from CpG islands. To search for genetic factors that regulate this variation, we performed a genome-wide association study (GWAS) mapping of methylation quantitative trait loci (mQTLs) for the 8590 testable CpG sites. cis association refers to correlation of methylation with SNPs within 1 Mb of a CpG site. 736 CpG sites showed phenotype-wide significant cis association with 2878 SNPs (after permutation correction for all tested markers and methylation phenotypes). In trans analysis of methylation, which tests for distant regulation effects, associations of 12 CpG sites and 38 SNPs remained significant after phenotype-wide correction. To examine the functional effects of mQTLs, we analyzed 85 genes that were with genetically regulated methylation we observed and for which we had quality gene expression data. Ten genes showed SNP-methylation-expression three-way associations—the same SNP simultaneously showed significant association with both DNA methylation and gene expression, while DNA methylation was significantly correlated with gene expression. Thus, we demonstrated that DNA methylation is frequently a heritable continuous quantitatively variable trait in human brain. Unlike allele-specific methylation, genetic polymorphisms mark both cis- and trans-regulatory genetic sites at measurable distances from their CpG sites. Some of the genetically regulated DNA methylation is directly connected with genetically regulated gene expression variation.     

DNA methylation regulates sclerostin (SOST) expression in osteoarthritic chondrocytes by bone morphogenetic protein 2 (BMP-2) induced changes in Smads binding affinity to the CpG region of SOST promoter

IntroductionSclerostin (SOST), a soluble antagonist of Wnt signaling, is expressed in chondrocytes and contributes to chondrocytes’ hypertrophic differentiation; however its role in osteoarthritis (OA) pathogenesis is not well known. Based on our previous findings on the interaction between Wnt/β-catenin pathway and BMP-2 in OA, we aimed to investigate the role of DNA methylation and BMP-2 on SOST’s expression in OA chondrocytes.MethodsSOST mRNA and protein expression levels were investigated using real-time polymerase chain reaction (PCR) and Western blot, respectively. The methylation status of SOST promoter was analysed using methylation-specific PCR (MSP), quantitative methylation-specific PCR (qMSP) and bisulfite sequencing analysis. The effect of BMP-2 and 5’-Aza-2-deoxycytidine (5-AzadC) on SOST’s expression levels were investigated and Smad1/5/8 binding to SOST promoter was assessed by Chromatin Immunoprecipitation (ChΙP).ResultsWe observed that SOST’s expression was upregulated in OA chondrocytes compared to normal. Moreover, we found that the CpG region of SOST promoter was hypomethylated in OA chondrocytes and 5-AzadC treatment in normal chondrocytes resulted in decreased SOST methylation, whereas its expression was upregulated. BMP-2 treatment in 5-AzadC-treated normal chondrocytes resulted in SOST upregulation, which was mediated through Smad 1/5/8 binding on the CpG region of SOST promoter.ConclusionsWe report novel findings that DNA methylation regulates SOST’s expression in OA, by changing Smad 1/5/8 binding affinity to SOST promoter, providing evidence that changes in DNA methylation pattern could underlie changes in genes’ expression observed in OA.     

Genome-Wide Analysis of DNA Methylation Differences in Muscle and Fat from Monozygotic Twins Discordant for Type 2 Diabetes

Background

Monozygotic twins discordant for type 2 diabetes constitute an ideal model to study environmental contributions to type 2 diabetic traits. We aimed to examine whether global DNA methylation differences exist in major glucose metabolic tissues from these twins.

Methodology/Principal Findings

Skeletal muscle (n = 11 pairs) and subcutaneous adipose tissue (n = 5 pairs) biopsies were collected from 53–80 year-old monozygotic twin pairs discordant for type 2 diabetes. DNA methylation was measured by microarrays at 26,850 cytosine-guanine dinucleotide (CpG) sites in the promoters of 14,279 genes. Bisulfite sequencing was applied to validate array data and to quantify methylation of intergenic repetitive DNA sequences. The overall intra-pair variation in DNA methylation was large in repetitive (LINE1, D4Z4 and NBL2) regions compared to gene promoters (standard deviation of intra-pair differences: 10% points vs. 4% points, P<0.001). Increased variation of LINE1 sequence methylation was associated with more phenotypic dissimilarity measured as body mass index (r = 0.77, P = 0.007) and 2-hour plasma glucose (r = 0.66, P = 0.03) whereas the variation in promoter methylation did not associate with phenotypic differences. Validated methylation changes were identified in the promoters of known type 2 diabetes-related genes, including PPARGC1A in muscle (13.9±6.2% vs. 9.0±4.5%, P = 0.03) and HNF4A in adipose tissue (75.2±3.8% vs. 70.5±3.7%, P<0.001) which had increased methylation in type 2 diabetic individuals. A hypothesis-free genome-wide exploration of differential methylation without correction for multiple testing identified 789 and 1,458 CpG sites in skeletal muscle and adipose tissue, respectively. These methylation changes only reached some percentage points, and few sites passed correction for multiple testing.

Conclusions/Significance

Our study suggests that likely acquired DNA methylation changes in skeletal muscle or adipose tissue gene promoters are quantitatively small between type 2 diabetic and non-diabetic twins. The importance of methylation changes in candidate genes such as PPARGC1A and HNF4A should be examined further by replication in larger samples.     

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.     

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.