Comparative anatomy of chromosomal domains with imprinted and non-imprinted allele-specific DNA methylation

Regulation of nc886 (vtRNA2-1) RNAs is associated with cardiometabolic risk factors and diseases

Non-coding 886 (nc886, vtRNA2-1) is a polymorphically imprinted gene. The methylation status of this locus has been shown to be associated with periconceptional conditions, and both the methylation status and the levels of nc886 RNAs have been shown to associate with later-life health traits. We have previously shown that nc886 RNA levels are associated not only with the methylation status of the locus, but also with a genetic polymorphism upstream from the locus. In this study, we describe the genetic and epigenetic regulators that predict lifelong nc886 RNA levels, as well as their association with cardiometabolic disease (CMD) risk factors and events. We utilised six population cohorts and one CMD cohort comprising 9058 individuals in total. The association of nc886 RNA levels, as predicted by epigenetic and genetic regulators, with CMD phenotypes was analysed using regression models, with a meta-analysis of the results. The meta-analysis showed that individuals with upregulated nc886 RNA levels have higher diastolic blood pressure (β = 0.07, p = 0.008), lower HDL levels (β =  - 0.07, p = 0.006) and an increased incidence of type 2 diabetes (OR = 1.260, p = 0.013). Moreover, CMD patients with upregulated nc886 RNA levels have an increased incidence of stroke (OR = 1.581, p = 0.006) and death (OR = 1.290, p = 0.046). In conclusion, we show that individuals who are predicted to present elevated nc886 RNA levels have poorer cardiovascular health and are at an elevated risk of complications in secondary prevention. This unique mechanism yields metabolic variation in human populations, constituting a CMD risk factor that cannot be modified through lifestyle choices.

Non-coding 886 (nc886/vtRNA2-1), the epigenetic odd duck - implications for future studies

Non-coding 886 (nc886, vtRNA2-1) is the only human polymorphically imprinted gene, in which the methylation status is not determined by genetics. Existing literature regarding the establishment, stability and consequences of the methylation pattern, as well as the nature and function of the nc886 RNAs transcribed from the locus, are contradictory. For example, the methylation status of the locus has been reported to be stable through life and across somatic tissues, but also susceptible to environmental effects. The nature of the produced nc886 RNA(s) has been redefined multiple times, and in carcinogenesis, these RNAs have been reported to have conflicting roles. In addition, due to the bimodal methylation pattern of the nc886 locus, traditional genome-wide methylation analyses can lead to false-positive results, especially in smaller datasets. Herein, we aim to summarize the existing literature regarding nc886, discuss how the characteristics of nc886 give rise to contradictory results, as well as to reinterpret, reanalyse and, where possible, replicate the results presented in the current literature. We also introduce novel findings on how the distribution of the nc886 methylation pattern is associated with the geographical origins of the population and describe the methylation changes in a large variety of human tumours. Through the example of this one peculiar genetic locus and RNA, we aim to highlight issues in the analysis of DNA methylation and non-coding RNAs in general and offer our suggestions for what should be taken into consideration in future analyses.

The mystique of epigenetic regulation: the remarkable case of a human noncoding RNA, nc886

nc886 is a regulatory noncoding RNA that is transcribed by RNA polymerase III (Pol III), is variably expressed in different biological contexts, and plays roles in inflammation and cancer. Epigenetic mechanisms play an intriguing role in regulating nc886 expression. As a maternally imprinted gene and metastable epiallele, nc866 exhibits polymorphic imprinting, with a methylation status that is influenced by environmental and biological factors. Consequently, the promoter DNA methylation status and the different resulting RNA expression levels of nc886 are associated with physiological and pathological conditions. In this review, we summarize the literature and explore the significance in relation to diverse roles of nc886.

Oocyte age and preconceptual alcohol use are highly correlated with epigenetic imprinting of a noncoding RNA (nc886)

Genomic imprinting occurs before fertilization, impacts every cell of the developing child, and may be sensitive to environmental perturbations. The noncoding RNA, nc886 (also called VTRNA2-1) is the only known example of the ∼100 human genes imprinted by DNA methylation, that shows polymorphic imprinting in the population. The nc886 gene is part of an ∼1.6-kb differentially methylated region (DMR) that is methylated in the oocyte and silenced on the maternal allele in about 75% of humans worldwide. Here, we show that the presence or absence of imprinting at the nc886 DMR in an individual is consistent across different tissues, confirming that the imprint is established before cellular differentiation and is maintained into adulthood. We investigated the relationships between the frequency of imprinting in newborns and maternal age, alcohol consumption and cigarette smoking before conception in more than 1,100 mother/child pairs from South Africa. The probability of imprinting in newborns was increased in older mothers and decreased in mothers who drank alcohol before conception. On the other hand, cigarette smoking had no apparent relationship with the frequency of imprinting. These data show an epigenetic change during oocyte maturation which is potentially subject to environmental influence. Much focus has been placed on avoiding alcohol consumption during pregnancy, but our data suggest that drinking before conception may affect the epigenome of the newborn.

Methylation status of nc886 epiallele reflects periconceptional conditions and is associated with glucose metabolism through nc886 RNAs

BACKGROUND

Non-coding RNA 886 (nc886) is coded from a maternally inherited metastable epiallele. We set out to investigate the determinants and dynamics of the methylation pattern at the nc886 epiallele and how this methylation status associates with nc886 RNA expression. Furthermore, we investigated the associations between the nc886 methylation status or the levels of nc886 RNAs and metabolic traits in the YFS and KORA cohorts. The association between nc886 epiallele methylation and RNA expression was also validated in induced pluripotent stem cell (iPSC) lines.

RESULTS

We confirm that the methylation status of the nc886 epiallele is mostly binomial, with individuals displaying either a non- or hemi-methylated status, but we also describe intermediately and close to fully methylated individuals. We show that an individual's methylation status is associated with the mother's age and socioeconomic status, but not with the individual's own genetics. Once established, the methylation status of the nc886 epiallele remains stable for at least 25 years. This methylation status is strongly associated with the levels of nc886 non-coding RNAs in serum, blood, and iPSC lines. In addition, nc886 methylation status associates with glucose and insulin levels during adolescence but not with the indicators of glucose metabolism or the incidence of type 2 diabetes in adulthood. However, the nc886-3p RNA levels also associate with glucose metabolism in adulthood.

CONCLUSIONS

These results indicate that nc886 metastable epiallele methylation is tuned by the periconceptional conditions and it associates with glucose metabolism through the expression of the ncRNAs coded in the epiallele region.

Methylation status of nc886 epiallele reflects periconceptional conditions and is associated with glucose metabolism through nc886 RNAs

BACKGROUND

Non-coding RNA 886 (nc886) is coded from a maternally inherited metastable epiallele. We set out to investigate the determinants and dynamics of the methylation pattern at the nc886 epiallele and how this methylation status associates with nc886 RNA expression. Furthermore, we investigated the associations between the nc886 methylation status or the levels of nc886 RNAs and metabolic traits in the YFS and KORA cohorts. The association between nc886 epiallele methylation and RNA expression was also validated in induced pluripotent stem cell (iPSC) lines.

RESULTS

We confirm that the methylation status of the nc886 epiallele is mostly binomial, with individuals displaying either a non- or hemi-methylated status, but we also describe intermediately and close to fully methylated individuals. We show that an individual's methylation status is associated with the mother's age and socioeconomic status, but not with the individual's own genetics. Once established, the methylation status of the nc886 epiallele remains stable for at least 25 years. This methylation status is strongly associated with the levels of nc886 non-coding RNAs in serum, blood, and iPSC lines. In addition, nc886 methylation status associates with glucose and insulin levels during adolescence but not with the indicators of glucose metabolism or the incidence of type 2 diabetes in adulthood. However, the nc886-3p RNA levels also associate with glucose metabolism in adulthood.

CONCLUSIONS

These results indicate that nc886 metastable epiallele methylation is tuned by the periconceptional conditions and it associates with glucose metabolism through the expression of the ncRNAs coded in the epiallele region.

Pan-cancer chromatin analysis of the human vtRNA genes uncovers their association with cancer biology

Background: The vault RNAs (vtRNAs) are a class of 84-141-nt eukaryotic non-coding RNAs transcribed by RNA polymerase III, associated to the ribonucleoprotein complex known as vault particle. Of the four human vtRNA genes, vtRNA1-1, vtRNA1-2 and vtRNA1-3, clustered at locus 1, are integral components of the vault particle, while vtRNA2-1 is a more divergent homologue located in a second locus. Gene expression studies of vtRNAs in large cohorts have been hindered by their unsuccessful sequencing using conventional transcriptomic approaches.

Methods: VtRNA expression in The Cancer Genome Atlas (TCGA) Pan-Cancer cohort was estimated using the genome-wide DNA methylation and chromatin accessibility data (ATAC-seq) of their genes as surrogate variables. The association between vtRNA expression and patient clinical outcome, immune subtypes and transcriptionally co-regulated gene programs was analyzed in the dataset.

Results: VtRNAs promoters are enriched in transcription factors related to viral infection. VtRNA2-1 is likely the most independently regulated homologue. VtRNA1-1 has the most accessible chromatin, followed by vtRNA1-2, vtRNA2-1 and vtRNA1-3. VtRNA1-1 and vtRNA1-3 chromatin status does not significantly change in cancer tissues. Meanwhile, vtRNA2-1 and vtRNA1-2 expression is widely deregulated in neoplastic tissues and its alteration is compatible with a broad oncogenic role for vtRNA1-2, and both tumor suppressor and oncogenic functions for vtRNA2-1. Yet, vtRNA1-1, vtRNA1-2 and vtRNA2-1 promoter DNA methylation predicts a shorter patient overall survival cancer-wide. In addition, gene ontology analyses of vtRNAs co-regulated genes identify a chromosome regulatory domain, epithelial differentiation, immune and thyroid cancer gene sets for specific vtRNAs. Furthermore, vtRNA expression patterns are associated with cancer immune subtypes and vtRNA1-2 expression is positively associated with cell proliferation and wound healing.

Conclusions: Our study presents the landscape of vtRNA chromatin status cancer-wide, identifying co-regulated gene networks and ontological pathways associated with the different vtRNA genes that may account for their diverse roles in cancer.

Identification of methylation changes associated with positive and negative growth deviance in Gambian infants using a targeted methyl sequencing approach of genomic DNA

Low birthweight and reduced height gain during infancy (stunting) may arise at least in part from adverse early life environments that trigger epigenetic reprogramming that may favor survival. We examined differential DNA methylation patterns using targeted methyl sequencing of regions regulating gene activity in groups of rural Gambian infants: (a) low and high birthweight (DNA from cord blood (n = 16 and n = 20, respectively), from placental trophoblast tissue (n = 21 and n = 20, respectively), and DNA from peripheral blood collected from infants at 12 months of age (n = 23 and n = 17, respectively)), and, (b) the top 10% showing rapid postnatal length gain (high, n = 20) and the bottom 10% showing slow postnatal length gain (low, n = 20) based on z score change between birth and 12 months of age (LAZ) (DNA from peripheral blood collected from infants at 12 months of age). Using BiSeq analysis to identify significant methylation marks, for birthweight, four differentially methylated regions (DMRs) were identified in trophoblast DNA, compared to 68 DMRs in cord blood DNA, and 54 DMRs in 12‐month peripheral blood DNA. Twenty‐five DMRs were observed to be associated with high and low length for age (LAZ) at 12 months. With the exception of five loci (associated with two different genes), there was no overlap between these groups of methylation marks. Of the 194 CpG methylation marks contained within DMRs, 106 were located to defined gene regulatory elements (promoters, CTCF‐binding sites, transcription factor‐binding sites, and enhancers), 58 to gene bodies (introns or exons), and 30 to intergenic DNA. Distinct methylation patterns associated with birthweight between comparison groups were observed in DNA collected at birth (at the end of intrauterine growth window) compared to those established by 12 months (near the infancy/childhood growth transition). The longitudinal differences in methylation patterns may arise from methylation adjustments, changes in cellular composition of blood or both that continue during the critical postnatal growth period, and in response to early nutritional and infectious environmental exposures with impacts on growth and longer‐term health outcomes.

Allele-specific DNA methylation is increased in cancers and its dense mapping in normal plus neoplastic cells increases the yield of disease-associated regulatory SNPs

BACKGROUND

Mapping of allele-specific DNA methylation (ASM) can be a post-GWAS strategy for localizing regulatory sequence polymorphisms (rSNPs). The advantages of this approach, and the mechanisms underlying ASM in normal and neoplastic cells, remain to be clarified.

RESULTS

We perform whole genome methyl-seq on diverse normal cells and tissues and three cancer types. After excluding imprinting, the data pinpoint 15,112 high-confidence ASM differentially methylated regions, of which 1838 contain SNPs in strong linkage disequilibrium or coinciding with GWAS peaks. ASM frequencies are increased in cancers versus matched normal tissues, due to widespread allele-specific hypomethylation and focal allele-specific hypermethylation in poised chromatin. Cancer cells show increased allele switching at ASM loci, but disruptive SNPs in specific classes of CTCF and transcription factor binding motifs are similarly correlated with ASM in cancer and non-cancer. Rare somatic mutations affecting these same motif classes track with de novo ASM. Allele-specific transcription factor binding from ChIP-seq is enriched among ASM loci, but most ASM differentially methylated regions lack such annotations, and some are found in otherwise uninformative "chromatin deserts."

CONCLUSIONS

ASM is increased in cancers but occurs by a shared mechanism involving disruptive SNPs in CTCF and transcription factor binding sites in both normal and neoplastic cells. Dense ASM mapping in normal plus cancer samples reveals candidate rSNPs that are difficult to find by other approaches. Together with GWAS data, these rSNPs can nominate specific transcriptional pathways in susceptibility to autoimmune, cardiometabolic, neuropsychiatric, and neoplastic diseases.

Analyzing DNA methylation patterns in subjects diagnosed with schizophrenia using machine learning methods

Schizophrenia is a common mental disorder with high heritability. It is genetically complex and to date more than a hundred risk loci have been identified. Association of environmental factors and schizophrenia has also been reported, while epigenetic analyses have yielded ambiguous and sometimes conflicting results. Here, we analyzed fresh frozen post-mortem brain tissue from a cohort of 73 subjects diagnosed with schizophrenia and 52 control samples, using the Illumina Infinium HumanMethylation450 Bead Chip, to investigate genome-wide DNA methylation patterns in the two groups. Analysis of differential methylation was performed with the Bioconductor Minfi package and modern machine-learning and visualization techniques, which were shown previously to be successful in detecting and highlighting differentially methylated patterns in case-control studies. In this dataset, however, these methods did not uncover any significant signals discerning the patient group and healthy controls, suggesting that if there are methylation changes associated with schizophrenia, they are heterogeneous and complex with small effect.

Protein kinase R and its cellular regulators in cancer: An active player or a surveillant?

Protein kinase R (PKR), originally known as an antiviral protein, senses various stresses as well as pathogen-driven double-stranded RNAs. Thereby activated PKR provokes diverse downstream events, including eIF2α phosphorylation and nuclear factor kappa-light-chain-enhancer of activated B cells activation. Consequently, PKR induces apoptosis and inflammation, both of which are highly important in cancer as much as its original antiviral role. Therefore, cellular proteins and RNAs should tightly control PKR activity. PKR and its regulators are often dysregulated in cancer and it is undoubted that such dysregulation contributes to tumorigenesis. However, PKR's precise role in cancer is still in debate, due to incomprehensible and even contradictory data. In this review, we introduce important cellular PKR regulators and discuss about their roles in cancer. Among them, we pay particular attention to nc886, a PKR repressor noncoding RNA that has been identified relatively recently, because its expression pattern in cancer can explain interesting yet obscure oncologic aspects of PKR. Based on nc886 and its regulation of PKR, we have proposed a tumor surveillance model, which reconciles contradictory data about PKR in cancer. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications.

Mother-child transmission of epigenetic information by tunable polymorphic imprinting

First, our work provides critical biological interpretation of intermediate DNA methylation readouts at the nc886 differentially methylated region (DMR). nc886 was identified in multiple large-scale epigenome-wide association studies (EWAS) that did not recognize that this region acts as a contiguous DMR imposed by genomic imprinting, highlighting the need to reexamine several 450k data sets. Second, strict control of genomic imprinting was thought to be required for organismal viability. Reports of polymorphic imprinting are limited to specific tissue types such as placenta and brain. In blood and somatic tissues, we show nc886 imprinting is mosaic in the population and influenced by maternal environment.

Genomic imprinting mediated by DNA methylation restricts gene expression to a single allele determined by parental origin and is not generally considered to be under genetic or environmental influence. Here, we focused on a differentially methylated region (DMR) of approximately 1.9 kb that includes a 101-bp noncoding RNA gene (nc886/VTRNA2-1), which is maternally imprinted in ∼75% of humans. This is unlike other imprinted genes, which demonstrate monoallelic methylation in 100% of individuals. The DMR includes a CTCF binding site on the centromeric side defining the DMR boundary and is flanked by a CTCF binding site on the telomeric side. The centromeric CTCF binding site contains an A/C polymorphism (rs2346018); the C allele is associated with less imprinting. The frequency of imprinting of the nc886 DMR in infants was linked to at least two nongenetic factors, maternal age at delivery and season of conception. In a separate cohort, nc886 imprinting was associated with lower body mass index in children at 5 y of age. Thus, we propose that the imprinting status of the nc886 DMR is “tunable” in that it is associated with maternal haplotype and prenatal environment. This provides a potential mechanism for transmitting information, with phenotypic consequences, from mother to child.

Heritable methylation marks associated with breast and prostate cancer risk

BACKGROUND

DNA methylation can mimic the effects of germline mutations in cancer predisposition genes. Recently, we identified twenty-four heritable methylation marks associated with breast cancer risk. As breast and prostate cancer share genetic risk factors, including rare, high-risk mutations (eg, in BRCA2), we hypothesized that some of these heritable methylation marks might also be associated with the risk of prostate cancer.

METHODS

We studied 869 incident prostate cancers (430 aggressive and 439 non-aggressive) and 869 matched controls nested within a prospective cohort study. DNA methylation was measured in pre-diagnostic blood samples using the Illumina Infinium HM450K BeadChip. Conditional logistic regression models, adjusted for prostate cancer risk factors and blood cell composition, were used to estimate odds ratios and 95% confidence intervals for the association between the 24 methylation marks and the risk of prostate cancer.

RESULTS

Five methylation marks within the VTRNA2-1 promoter region (cg06536614, cg00124993, cg26328633, cg25340688, and cg26896946), and one in the body of CLGN (cg22901919) were associated with the risk of prostate cancer. In stratified analyses, the five VTRNA2-1 marks were associated with the risk of aggressive prostate cancer.

CONCLUSIONS

This work highlights a potentially important new area of investigation for prostate cancer susceptibility and adds to our knowledge about shared risk factors for breast and prostate cancer.

Allele-specific DNA methylation of disease susceptibility genes in Japanese patients with inflammatory bowel disease

BACKGROUND

Inflammatory bowel disease (IBD) has an unknown etiology; however, accumulating evidence suggests that IBD is a multifactorial disease influenced by a combination of genetic and environmental factors. The influence of genetic variants on DNA methylation in cis and cis effects on expression have been demonstrated. We hypothesized that IBD susceptibility single-nucleotide polymorphisms (SNPs) regulate susceptibility gene expressions in cis by regulating DNA methylation around SNPs. For this, we determined cis-regulated allele-specific DNA methylation (ASM) around IBD susceptibility genes in CD4+ effector/memory T cells (Tem) in lamina propria mononuclear cells (LPMCs) in patients with IBD and examined the association between the ASM SNP genotype and neighboring susceptibility gene expressions.

METHODS

CD4+ effector/memory T cells (Tem) were isolated from LPMCs in 15 Japanese IBD patients (ten Crohn's disease [CD] and five ulcerative colitis [UC] patients). ASM analysis was performed by methylation-sensitive SNP array analysis. We defined ASM as a changing average relative allele score ([Formula: see text]) >0.1 after digestion by methylation-sensitive restriction enzymes. Among SNPs showing [Formula: see text] >0.1, we extracted the probes located on tag-SNPs of 200 IBD susceptibility loci and around IBD susceptibility genes as candidate ASM SNPs. To validate ASM, bisulfite-pyrosequencing was performed. Transcriptome analysis was examined in 11 IBD patients (seven CD and four UC patients). The relation between rs36221701 genotype and neighboring gene expressions were analyzed.

RESULTS

We extracted six candidate ASM SNPs around IBD susceptibility genes. The top of [Formula: see text] (0.23) was rs1130368 located on HLA-DQB1. ASM around rs36221701 ([Formula: see text] = 0.14) located near SMAD3 was validated using bisulfite pyrosequencing. The SMAD3 expression was significantly associated with the rs36221701 genotype (p = 0.016).

CONCLUSIONS

We confirmed the existence of cis-regulated ASM around IBD susceptibility genes and the association between ASM SNP (rs36221701) genotype and SMAD3 expression, a susceptibility gene for IBD. These results give us supporting evidence that DNA methylation mediates genetic effects on disease susceptibility.

Maternal 5mCpG Imprints at the PARD6G-AS1 and GCSAML Differentially Methylated Regions Are Decoupled From Parent-of-Origin Expression Effects in Multiple Human Tissues

A hallmark of imprinted genes in mammals is the occurrence of parent-of-origin-dependent asymmetry of DNA cytosine methylation (5^mC) of alleles at CpG islands (CGIs) in their promoter regions. This 5^mCpG asymmetry between the parental alleles creates allele-specific imprinted differentially methylated regions (iDMRs). iDMRs are often coupled to the transcriptional repression of the methylated allele and the activation of the unmethylated allele in a tissue-specific, developmental-stage-specific and/or isoform-specific fashion. iDMRs function as regulatory platforms, built through the recruitment of chemical modifications to histones to achieve differential, parent-of-origin-dependent chromatin segmentation states. Here, we used a comparative computational data mining approach to identify 125 novel constitutive candidate iDMRs that integrate the maximal number of allele-specific methylation region records overlapping CGIs in human methylomes. Twenty-nine candidate iDMRs display gametic 5^mCpG asymmetry, and another 96 are candidate secondary iDMRs. We established the maternal origin of the 5^mCpG imprints of one gametic (PARD6G-AS1) and one secondary (GCSAML) iDMRs. We also found a constitutively hemimethylated, nonimprinted domain at the PWWP2AP1 promoter CGI with oocyte-derived methylation asymmetry. Given that the 5^mCpG level at the iDMRs is not a sufficient criterion to predict active or silent locus states and that iDMRs can regulate genes from a distance of more than 1 Mb, we used RNA-Seq experiments from the Genotype-Tissue Expression project and public archives to assess the transcriptional expression profiles of SNPs across 4.6 Mb spans around the novel maternal iDMRs. We showed that PARD6G-AS1 and GCSAML are expressed biallelically in multiple tissues. We found evidence of tissue-specific monoallelic expression of ZNF124 and OR2L13, located 363 kb upstream and 419 kb downstream, respectively, of the GCSAML iDMR. We hypothesize that the GCSAML iDMR regulates the tissue-specific, monoallelic expression of ZNF124 but not of OR2L13. We annotated the non-coding epigenomic marks in the two maternal iDMRs using data from the Roadmap Epigenomics project and showed that the PARD6G-AS1 and GCSAML iDMRs achieve contrasting activation and repression chromatin segmentations. Lastly, we found that the maternal 5^mCpG imprints are perturbed in several hematopoietic cancers. We conclude that the maternal 5^mCpG imprints at PARD6G-AS1 and GCSAML iDMRs are decoupled from parent-of-origin transcriptional expression effects in multiple tissues.

Heritable DNA methylation marks associated with susceptibility to breast cancer

Mendelian-like inheritance of germline DNA methylation in cancer susceptibility genes has been previously reported. We aimed to scan the genome for heritable methylation marks associated with breast cancer susceptibility by studying 25 Australian multiple-case breast cancer families. Here we report genome-wide DNA methylation measured in 210 peripheral blood DNA samples provided by family members using the Infinium HumanMethylation450. We develop and apply a new statistical method to identify heritable methylation marks based on complex segregation analysis. We estimate carrier probabilities for the 1000 most heritable methylation marks based on family structure, and we use Cox proportional hazards survival analysis to identify 24 methylation marks with corresponding carrier probabilities significantly associated with breast cancer. We replicate an association with breast cancer risk for four of the 24 marks using an independent nested case-control study. Here, we report a novel approach for identifying heritable DNA methylation marks associated with breast cancer risk.

Nc886 is epigenetically repressed in prostate cancer and acts as a tumor suppressor through the inhibition of cell growth

Background

Nc886 is a 102 bp non-coding RNA transcript initially classified as a microRNA precursor (Pre-miR-886), later as a divergent homologue of the vault RNAs (vtRNA 2–1) and more recently as a novel type of RNA (nc886). Although nc886/vtRNA2–1/Pre-miR-886 identity is still controversial, it was shown to be epigenetically controlled, presenting both tumor suppressor and oncogenic function in different cancers. Here, we study for the first time the role of nc886 in prostate cancer.

Methods

Nc886 promoter methylation status and its correlation with patient clinical parameters or DNMTs levels were evaluated in TCGA and specific GEO prostate tissue datasets. Nc886 level was measured by RT-qPCR to compare normal/neoplastic prostate cells from radical prostatectomies and cell lines, and to assess nc886 response to demethylating agents. The effect of nc886 recovery in cell proliferation (in vitro and in vivo) and invasion (in vitro) was evaluated using lentiviral transduced DU145 and LNCaP cell lines. The association between the expression of nc886 and selected genes was analyzed in the TCGA-PRAD cohort.

Results

Nc886 promoter methylation increases in tumor vs. normal prostate tissue, as well as in metastatic vs. normal prostate tissue. Additionally, nc886 promoter methylation correlates with prostate cancer clinical staging, including biochemical recurrence, Clinical T-value and Gleason score. Nc886 transcript is downregulated in tumor vs. normal tissue -in agreement with its promoter methylation status- and increases upon demethylating treatment. In functional studies, the overexpression of nc886 in the LNCaP and DU145 cell line leads to a decreased in vitro cell proliferation and invasion, as well as a reduced in vivo cell growth in NUDE-mice tumor xenografts. Finally, nc886 expression associates with the prostate cancer cell cycle progression gene signature in TCGA-PRAD.

Conclusions

Our data suggest a tumor suppressor role for nc886 in the prostate, whose expression is epigenetically silenced in cancer leading to an increase in cell proliferation and invasion. Nc886 might hold clinical value in prostate cancer due to its association with clinical parameters and with a clinically validated gene signature.

Electronic supplementary material

The online version of this article (10.1186/s12885-018-4049-7) contains supplementary material, which is available to authorized users.

Regions of common inter-individual DNA methylation differences in human monocytes: genetic basis and potential function

Background

There is increasing evidence for inter-individual methylation differences at CpG dinucleotides in the human genome, but the regional extent and function of these differences have not yet been studied in detail. For identifying regions of common methylation differences, we used whole genome bisulfite sequencing data of monocytes from five donors and a novel bioinformatic strategy.

Results

We identified 157 differentially methylated regions (DMRs) with four or more CpGs, almost none of which has been described before. The DMRs fall into different chromatin states, where methylation is inversely correlated with active, but not repressive histone marks. However, methylation is not correlated with the expression of associated genes. High-resolution single nucleotide polymorphism (SNP) genotyping of the five donors revealed evidence for a role of cis-acting genetic variation in establishing methylation patterns. To validate this finding in a larger cohort, we performed genome-wide association studies (GWAS) using SNP genotypes and 450k array methylation data from blood samples of 1128 individuals. Only 30/157 (19%) DMRs include at least one 450k CpG, which shows that these arrays miss a large proportion of DNA methylation variation. In most cases, the GWAS peak overlapped the CpG position, and these regions are enriched for CREB group, NF-1, Sp100 and CTCF binding motifs. In two cases, there was tentative evidence for a trans-effect by KRAB zinc finger proteins.

Conclusions

Allele-specific DNA methylation occurs in discrete chromosomal regions and is driven by genetic variation in cis and trans, but in general has little effect on gene expression.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-017-0144-2) contains supplementary material, which is available to authorized users.

The methyltransferase SETDB1 regulates a large neuron-specific topological chromatin domain

We report locus-specific disintegration of megabase-scale chromosomal conformations in brain after neuronal ablation of Kmt1e/Setdb1 histone H3-lysine 9 methyltransferase, including a large topologically associated 1.2Mb domain conserved in human and mouse and encompassing >70 genes at the clustered Protocadherin (cPcdh) locus. TAD^cPcdh in mutant neurons showed abnormal accumulations of CTCF transcriptional regulator and 3D genome organizer at cryptic binding sites, converted into permissive state with DNA cytosine hypomethylation and histone hyperacetylation. Broadly upregulated expression across cPcdh included defective S-type Protocadherin single-cell stochastic constraint. Setdb1-dependent loop formations, bypassing 0.2–1Mb of linear genome, radiated from TAD^Pcdh fringes towards cPcdh cis-regulatory sequences, counterbalanced shorter-range facilitative promoter-enhancer contacts and carried loop-bound polymorphisms associated with genetic risk for schizophrenia. We show that KRAB-zinc finger Setdb1 repressor complex, shielding neuronal 3D genomes from excess CTCF binding, is critically required for structural maintenance of TAD^cPcdh.

Genetic-epigenetic interactions in cis: a major focus in the post-GWAS era

Studies on genetic-epigenetic interactions, including the mapping of methylation quantitative trait loci (mQTLs) and haplotype-dependent allele-specific DNA methylation (hap-ASM), have become a major focus in the post-genome-wide-association-study (GWAS) era. Such maps can nominate regulatory sequence variants that underlie GWAS signals for common diseases, ranging from neuropsychiatric disorders to cancers. Conversely, mQTLs need to be filtered out when searching for non-genetic effects in epigenome-wide association studies (EWAS). Sequence variants in CCCTC-binding factor (CTCF) and transcription factor binding sites have been mechanistically linked to mQTLs and hap-ASM. Identifying these sites can point to disease-associated transcriptional pathways, with implications for targeted treatment and prevention.

DNA Methylation Signatures of Early Childhood Malnutrition Associated With Impairments in Attention and Cognition

BACKGROUND

Early childhood malnutrition affects 113 million children worldwide, impacting health and increasing vulnerability for cognitive and behavioral disorders later in life. Molecular signatures after childhood malnutrition, including the potential for intergenerational transmission, remain unexplored.

METHODS

We surveyed blood DNA methylomes (~483,000 individual CpG sites) in 168 subjects across two generations, including 50 generation 1 individuals hospitalized during the first year of life for moderate to severe protein-energy malnutrition, then followed up to 48 years in the Barbados Nutrition Study. Attention deficits and cognitive performance were evaluated with the Connors Adult Attention Rating Scale and Wechsler Abbreviated Scale of Intelligence. Expression of nutrition-sensitive genes was explored by quantitative reverse transcriptase polymerase chain reaction in rat prefrontal cortex.

RESULTS

We identified 134 nutrition-sensitive, differentially methylated genomic regions, with most (87%) specific for generation 1. Multiple neuropsychiatric risk genes, including COMT, IFNG, MIR200B, SYNGAP1, and VIPR2 showed associations of specific methyl-CpGs with attention and IQ. IFNG expression was decreased in prefrontal cortex of rats showing attention deficits after developmental malnutrition.

CONCLUSIONS

Early childhood malnutrition entails long-lasting epigenetic signatures associated with liability for attention and cognition, and limited potential for intergenerational transmission.

Mechanisms and Disease Associations of Haplotype-Dependent Allele-Specific DNA Methylation

Haplotype-dependent allele-specific methylation (hap-ASM) can impact disease susceptibility, but maps of this phenomenon using stringent criteria in disease-relevant tissues remain sparse. Here we apply array-based and Methyl-Seq approaches to multiple human tissues and cell types, including brain, purified neurons and glia, T lymphocytes, and placenta, and identify 795 hap-ASM differentially methylated regions (DMRs) and 3,082 strong methylation quantitative trait loci (mQTLs), most not previously reported. More than half of these DMRs have cell type-restricted ASM, and among them are 188 hap-ASM DMRs and 933 mQTLs located near GWAS signals for immune and neurological disorders. Targeted bis-seq confirmed hap-ASM in 12/13 loci tested, including CCDC155, CD69, FRMD1, IRF1, KBTBD11, and S100A(∗)-ILF2, associated with immune phenotypes, MYT1L, PTPRN2, CMTM8 and CELF2, associated with neurological disorders, NGFR and HLA-DRB6, associated with both immunological and brain disorders, and ZFP57, a trans-acting regulator of genomic imprinting. Polymorphic CTCF and transcription factor (TF) binding sites were over-represented among hap-ASM DMRs and mQTLs, and analysis of the human data, supplemented by cross-species comparisons to macaques, indicated that CTCF and TF binding likelihood predicts the strength and direction of the allelic methylation asymmetry. These results show that hap-ASM is highly tissue specific; an important trans-acting regulator of genomic imprinting is regulated by this phenomenon; and variation in CTCF and TF binding sites is an underlying mechanism, and maps of hap-ASM and mQTLs reveal regulatory sequences underlying supra- and sub-threshold GWAS peaks in immunological and neurological disorders.

Tissue-specific patterns of allelically-skewed DNA methylation

While DNA methylation is usually thought to be symmetrical across both alleles, there are some notable exceptions. Genomic imprinting and X chromosome inactivation are two well-studied sources of allele-specific methylation (ASM), but recent research has indicated a more complex pattern in which genotypic variation can be associated with allelically-skewed DNA methylation in cis. Given the known heterogeneity of DNA methylation across tissues and cell types we explored inter- and intra-individual variation in ASM across several regions of the human brain and whole blood from multiple individuals. Consistent with previous studies, we find widespread ASM with > 4% of the ∼220,000 loci interrogated showing evidence of allelically-skewed DNA methylation. We identify ASM flanking known imprinted regions, and show that ASM sites are enriched in DNase I hypersensitivity sites and often located in an extended genomic context of intermediate DNA methylation. We also detect examples of genotype-driven ASM, some of which are tissue-specific. These findings contribute to our understanding of the nature of differential DNA methylation across tissues and have important implications for genetic studies of complex disease. As a resource to the community, ASM patterns across each of the tissues studied are available in a searchable online database: http://epigenetics.essex.ac.uk/ASMBrainBlood.

Implication of a Chromosome 15q15.2 Locus in Regulating UBR1 and Predisposing Smokers to MGMT Methylation in Lung

O(6)-Methylguanine-DNA methyltransferase (MGMT) is a DNA repair enzyme that protects cells from carcinogenic effects of alkylating agents; however, MGMT is silenced by promoter hypermethylation during carcinogenesis. A single-nucleotide polymorphism (SNP) in an enhancer in the MGMT promoter was previously identified to be highly significantly associated with risk for MGMT methylation in lung cancer and sputum from smokers. To further genetic investigations, a genome-wide association and replication study was conducted in two smoker cohorts to identify novel loci for MGMT methylation in sputum that were independent of the MGMT enhancer polymorphism. Two novel trans-acting loci (15q15.2 and 17q24.3) that were identified acted together with the enhancer SNP to empower risk prediction for MGMT methylation. We found that the predisposition to MGMT methylation arising from the 15q15.2 locus involved regulation of the ubiquitin protein ligase E3 component UBR1. UBR1 attenuation reduced turnover of MGMT protein and increased repair of O6-methylguanine in nitrosomethylurea-treated human bronchial epithelial cells, while also reducing MGMT promoter activity and abolishing MGMT induction. Overall, our results substantiate reduced gene transcription as a major mechanism for predisposition to MGMT methylation in the lungs of smokers, and support the importance of UBR1 in regulating MGMT homeostasis and DNA repair of alkylated DNA adducts in cells.

Independent genomewide screens identify the tumor suppressor VTRNA2-1 as a human epiallele responsive to periconceptional environment

BACKGROUND

Interindividual epigenetic variation that occurs systemically must be established prior to gastrulation in the very early embryo and, because it is systemic, can be assessed in easily biopsiable tissues. We employ two independent genome-wide approaches to search for such variants.

RESULTS

First, we screen for metastable epialleles by performing genomewide bisulfite sequencing in peripheral blood lymphocyte (PBL) and hair follicle DNA from two Caucasian adults. Second, we conduct a genomewide screen for genomic regions at which PBL DNA methylation is affected by season of conception in rural Gambia. Remarkably, both approaches identify the genomically imprinted VTRNA2-1 as a top environmentally responsive epiallele. We demonstrate systemic and stochastic interindividual variation in DNA methylation at the VTRNA2-1 differentially methylated region in healthy Caucasian and Asian adults and show, in rural Gambians, that periconceptional environment affects offspring VTRNA2-1 epigenotype, which is stable over at least 10 years. This unbiased screen also identifies over 100 additional candidate metastable epialleles, and shows that these are associated with cis genomic features including transposable elements.

CONCLUSIONS

The non-coding VTRNA2-1 transcript (also called nc886) is a putative tumor suppressor and modulator of innate immunity. Thus, these data indicating environmentally induced loss of imprinting at VTRNA2-1 constitute a plausible causal pathway linking early embryonic environment, epigenetic alteration, and human disease. More broadly, the list of candidate metastable epialleles provides a resource for future studies of epigenetic variation and human disease.

Multiplex bisulfite PCR resequencing of clinical FFPE DNA

BACKGROUND

The clinical utility of DNA methylation as a predictive or prognostic biomarker requires scalable resequencing protocols for bisulfite-converted DNA. Key features of any validation method should be adaptability for low- or high-throughput needs and high reproducibility, and should only require minimal amounts of precious clinical sample as input material. Critically, this method should also deliver robust results when working with bisulfite-converted DNA extracted from formalin-fixed, paraffin-embedded (FFPE) blocks.

RESULTS

We report here for the first time on comparison studies between the Fluidigm Access Array system and multiplex assays for multiplex bisulfite PCR resequencing. The requirement of the Fluidigm Access Array system for high template amounts and its sensitivity to variations in template quality rendered it unsuitable for bisulfite PCR applications utilizing FFPE DNA. In response to this limitation, we established a multiplex bisulfite PCR assay capable of delivering robust methylation data using minimal amounts of FFPE clinical DNA. To evaluate the parameters and reproducibility of this assay, 57 amplicons were used to prepare sequencing libraries in triplicate for 13 FFPE tumour samples, as well as a series of 5 methylated controls (0%, 25%, 50%, 75%, and 100%). Analysis of this data demonstrated that this multiplex assay had high reproducibility (mean standard deviation of 1.4% for methylation values), was low cost, required low sample input (50 ng of DNA or less), and could be scaled for both low- and high-throughput needs. Notably, ExoSAP-IT (exonuclease I) treatment to remove residual primers in bisulfite resequencing libraries appeared to degrade the library and generate a high-molecular weight smear which may impact on the degree of methylation assessed.

CONCLUSIONS

Multiplex bisulfite PCR assays represent a convenient and scalable method for validation and screening of methylated DNA regions from archival FFPE DNA. Moreover, the library construction process detailed here can be rapidly optimized and implemented with a minimal amount of work, can be performed using the standard equipment found in any molecular biology laboratory, and can be easily adapted for use on both genomic DNA and bisulfite DNA applications. However, in preparing bisulfite libraries for sequencing, the use of ExoSAP-IT is not recommended due to potential off-target nuclease effects which may impact downstream methylation analysis.

Genome-wide DNA methylation profile of leukocytes from melanoma patients with and without CDKN2A mutations

Melanoma is a highly aggressive cancer, accounting for up to 75% of skin cancer deaths. A small proportion of melanoma cases can be ascribed to the presence of highly penetrant germline mutations, and approximately 40% of hereditary melanoma cases are caused by CDKN2A mutations. The current study sought to investigate whether the presence of germline CDKN2A mutations or the occurrence of cutaneous melanoma would result in constitutive genome-wide DNA methylation changes. The leukocyte methylomes of two groups of melanoma patients (those with germline CDKN2A mutations and those without CDKN2A mutations) were analyzed together with the profile of a control group of individuals. A pattern of DNA hypomethylation was detected in the CDKN2A-negative patients relative to both CDKN2A-mutated patients and controls. Additionally, we delineated a panel of 90 CpG sites that were differentially methylated in CDKN2A-mutated patients relative to controls. Although we identified a possible constitutive epigenetic signature in CDKN2A-mutated patients, the occurrence of reported SNPs at the detected CpG sites complicated the data interpretation. Thus, further studies are required to elucidate the impact of these findings on melanoma predisposition and their possible effect on the penetrance of CDKN2A mutations.

Differentially methylated CpG island within human XIST mediates alternative P2 transcription and YY1 binding

Background

X-chromosome inactivation silences one X chromosome in females to achieve dosage compensation with the single X chromosome in males. While most genes are silenced on the inactive X chromosome, the gene for the long non-coding RNA XIST is silenced on the active X chromosome and expressed from the inactive X chromosome with which the XIST RNA associates, triggering silencing of the chromosome. In mouse, an alternative Xist promoter, P2 is also the site of YY1 binding, which has been shown to serve as a tether between the Xist RNA and the DNA of the chromosome. In humans there are many differences from the initial events of mouse Xist activation, including absence of a functional antisense regulator Tsix, and absence of strictly paternal inactivation in extraembryonic tissues, prompting us to examine regulatory regions for the human XIST gene.

Results

We demonstrate that the female-specific DNase hypersensitivity site within XIST is specific to the inactive X chromosome and correlates with transcription from an internal P2 promoter. P2 is located within a CpG island that is differentially methylated between males and females and overlaps conserved YY1 binding sites that are only bound on the inactive X chromosome where the sites are unmethylated. However, YY1 binding is insufficient to drive P2 expression or establish the DHS, which may require a development-specific factor. Furthermore, reduction of YY1 reduces XIST transcription in addition to causing delocalization of XIST.

Conclusions

The differentially methylated DNase hypersensitive site within XIST marks the location of an alternative promoter, P2, that generates a transcript of unknown function as it lacks the A repeats that are critical for silencing. In addition, this region binds YY1 on the unmethylated inactive X chromosome, and depletion of YY1 untethers the XIST RNA as well as decreasing transcription of XIST.

Genetic variants within the second intron of the KCNQ1 gene affect CTCF binding and confer a risk of Beckwith-Wiedemann syndrome upon maternal transmission

BACKGROUND

Disruption of 11p15 imprinting results in two fetal growth disorders with opposite phenotypes: the Beckwith-Wiedemann (BWS; MIM 130650) and the Silver-Russell (SRS; MIM 180860) syndromes. DNA methylation defects account for 60% of BWS and SRS cases and, in most cases, occur without any identified mutation in a cis-acting regulatory sequence or a trans-acting factor.

METHODS

We investigated whether 11p15 cis-acting sequence variants account for primary DNA methylation defects in patients with SRS and BWS with loss of DNA methylation at ICR1 and ICR2, respectively.

RESULTS

We identified a 4.5 kb haplotype that, upon maternal transmission, is associated with a risk of ICR2 loss of DNA methylation in patients with BWS. This novel region is located within the second intron of the KCNQ1 gene, 170 kb upstream of the ICR2 imprinting centre and encompasses two CTCF binding sites. We showed that, within the 4.5 kb region, two SNPs (rs11823023 and rs179436) affect CTCF occupancy at DNA motifs flanking the CTCF 20 bp core motif.

CONCLUSIONS

This study shows that genetic variants confer a risk of DNA methylation defect with a parent-of-origin effect and highlights the crucial role of CTCF for the regulation of genomic imprinting of the CDKN1C/KCNQ1 domain.

GeMes, clusters of DNA methylation under genetic control, can inform genetic and epigenetic analysis of disease

Epigenetic marks such as DNA methylation have generated great interest in the study of human disease. However, studies of DNA methylation have not established population-epigenetics principles to guide design, efficient statistics, or interpretation. Here, we show that the clustering of correlated DNA methylation at CpGs was similar to that of linkage-disequilibrium (LD) correlation in genetic SNP variation but for much shorter distances. Some clustering of methylated CpGs appeared to be genetically driven. Further, a set of correlated methylated CpGs related to a single SNP-based LD block was not always physically contiguous-segments of uncorrelated methylation as long as 300 kb could be interspersed in the cluster. Thus, we denoted these sets of correlated CpGs as GeMes, defined as potentially noncontiguous methylation clusters under the control of one or more methylation quantitative trait loci. This type of correlated methylation structure has implications for both biological functions of DNA methylation and for the design, analysis, and interpretation of epigenome-wide association studies.

Variable maternal methylation overlapping the nc886/vtRNA2-1 locus is locked between hypermethylated repeats and is frequently altered in cancer

Cancer is as much an epigenetic disease as a genetic one; however, the interplay between these two processes is unclear. Recently, it has been shown that a large proportion of DNA methylation variability can be explained by allele-specific methylation (ASM), either at classical imprinted loci or those regulated by underlying genetic variants. During a recent screen for imprinted differentially methylated regions, we identified the genomic interval overlapping the non-coding nc886 RNA (previously known as vtRNA2-1) as an atypical ASM that shows variable levels of methylation, predominantly on the maternal allele in many tissues. Here we show that the nc886 interval is the first example of a polymorphic imprinted DMR in humans. Further analysis of the region suggests that the interval subjected to ASM is approximately 2 kb in size and somatically acquired. An in depth analysis of this region in primary cancer samples with matching normal adjacent tissue from the Cancer Genome Atlas revealed that aberrant methylation in bladder, breast, colon and lung tumors occurred in approximately 27% of cases. Hypermethylation occurred more frequently than hypomethylation. Using additional normal-tumor paired samples we show that on rare occasions the aberrant methylation profile is due to loss-of-heterozygosity. This work therefore suggests that the nc886 locus is subject to variable allelic methylation that undergoes cancer-associated epigenetic changes in solid tumors.