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

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.

Vault RNAs (vtRNAs): Rediscovered non-coding RNAs with diverse physiological and pathological activities

The physicochemical characteristics of RNA admit non-coding RNAs to perform a different range of biological acts through various mechanisms and are involved in regulating a diversity of fundamental processes. Notably, some reports of pathological conditions have proved abnormal expression of many non-coding RNAs guides the ailment. Vault RNAs are a class of non-coding RNAs containing stem regions or loops with well-conserved sequence patterns that play a fundamental role in the function of vault particles through RNA-ligand, RNA-RNA, or RNA-protein interactions. Taken together, vault RNAs have been proposed to be involved in a variety of functions such as cell proliferation, nucleocytoplasmic transport, intracellular detoxification processes, multidrug resistance, apoptosis, and autophagy, and serve as microRNA precursors and signaling pathways. Despite decades of investigations devoted, the biological function of the vault particle or the vault RNAs is not yet completely cleared. In this review, the current scientific assertions of the vital vault RNAs functions were discussed.

A POLR3B-variant reveals a Pol III transcriptome response dependent on La protein/SSB

RNA polymerase III (Pol III, POLR3) synthesizes tRNAs and other small non-coding RNAs. Human POLR3 pathogenic variants cause a range of developmental disorders, recapitulated in part by mouse models, yet some aspects of POLR3 deficiency have not been explored. We characterized a human POLR3B:c.1625A>G;p.(Asn542Ser) disease variant that was found to cause mis-splicing of POLR3B. Genome-edited POLR3B1625A>G HEK293 cells acquired the mis-splicing with decreases in multiple POLR3 subunits and TFIIIB, although display auto-upregulation of the Pol III termination-reinitiation subunit POLR3E. La protein was increased relative to its abundant pre-tRNA ligands which bind via their U(n)U-3'-termini. Assays for cellular transcription revealed greater deficiencies for tRNA genes bearing terminators comprised of 4Ts than of ≥5Ts. La-knockdown decreased Pol III ncRNA expression unlinked to RNA stability. Consistent with these effects, small-RNAseq showed that POLR3B1625A>G and patient fibroblasts express more tRNA fragments (tRFs) derived from pre-tRNA 3'-trailers (tRF-1) than from mature-tRFs, and higher levels of multiple miRNAs, relative to control cells. The data indicate that decreased levels of Pol III transcripts can lead to functional excess of La protein which reshapes small ncRNA profiles revealing new depth in the Pol III system. Finally, patient cell RNA analysis uncovered a strategy for tRF-1/tRF-3 as POLR3-deficiency biomarkers.

ICF1-Syndrome-Associated DNMT3B Mutations Prevent De Novo Methylation at a Subset of Imprinted Loci during iPSC Reprogramming

Parent-of-origin-dependent gene expression of a few hundred human genes is achieved by differential DNA methylation of both parental alleles. This imprinting is required for normal development, and defects in this process lead to human disease. Induced pluripotent stem cells (iPSCs) serve as a valuable tool for in vitro disease modeling. However, a wave of de novo DNA methylation during reprogramming of iPSCs affects DNA methylation, thus limiting their use. The DNA methyltransferase 3B (DNMT3B) gene is highly expressed in human iPSCs; however, whether the hypermethylation of imprinted loci depends on DNMT3B activity has been poorly investigated. To explore the role of DNMT3B in mediating de novo DNA methylation at imprinted DMRs, we utilized iPSCs generated from patients with immunodeficiency, centromeric instability, facial anomalies type I (ICF1) syndrome that harbor biallelic hypomorphic DNMT3B mutations. Using a whole-genome array-based approach, we observed a gain of methylation at several imprinted loci in control iPSCs but not in ICF1 iPSCs compared to their parental fibroblasts. Moreover, in corrected ICF1 iPSCs, which restore DNMT3B enzymatic activity, imprinted DMRs did not acquire control DNA methylation levels, in contrast to the majority of the hypomethylated CpGs in the genome that were rescued in the corrected iPSC clones. Overall, our study indicates that DNMT3B is responsible for de novo methylation of a subset of imprinted DMRs during iPSC reprogramming and suggests that imprinting is unstable during a specific time window of this process, after which the epigenetic state at these regions becomes resistant to perturbation.

The influence of early environment and micronutrient availability on developmental epigenetic programming: lessons from the placenta

DNA methylation is the most commonly studied epigenetic mark in humans, as it is well recognised as a stable, heritable mark that can affect genome function and influence gene expression. Somatic DNA methylation patterns that can persist throughout life are established shortly after fertilisation when the majority of epigenetic marks, including DNA methylation, are erased from the pre-implantation embryo. Therefore, the period around conception is potentially critical for influencing DNA methylation, including methylation at imprinted alleles and metastable epialleles (MEs), loci where methylation varies between individuals but is correlated across tissues. Exposures before and during conception can affect pregnancy outcomes and health throughout life. Retrospective studies of the survivors of famines, such as those exposed to the Dutch Hunger Winter of 1944-45, have linked exposures around conception to later disease outcomes, some of which correlate with DNA methylation changes at certain genes. Animal models have shown more directly that DNA methylation can be affected by dietary supplements that act as cofactors in one-carbon metabolism, and in humans, methylation at birth has been associated with peri-conceptional micronutrient supplementation. However, directly showing a role of micronutrients in shaping the epigenome has proven difficult. Recently, the placenta, a tissue with a unique hypomethylated methylome, has been shown to possess great inter-individual variability, which we highlight as a promising target tissue for studying MEs and mixed environmental exposures. The placenta has a critical role shaping the health of the fetus. Placenta-associated pregnancy complications, such as preeclampsia and intrauterine growth restriction, are all associated with aberrant patterns of DNA methylation and expression which are only now being linked to disease risk later in life.

Methylation status of VTRNA2-1/nc886 is stable across populations, monozygotic twin pairs and in majority of tissues

Aims & methods: The aim of this study was to characterize the methylation level of a polymorphically imprinted gene, VTRNA2-1/nc886, in human populations and somatic tissues.48 datasets, consisting of more than 30 tissues and >30,000 individuals, were used. Results: nc886 methylation status is associated with twin status and ethnic background, but the variation between populations is limited. Monozygotic twin pairs present concordant methylation, whereas ∼30% of dizygotic twin pairs present discordant methylation in the nc886 locus. The methylation levels of nc886 are uniform across somatic tissues, except in cerebellum and skeletal muscle. Conclusion: The nc886 imprint may be established in the oocyte, and, after implantation, the methylation status is stable, excluding a few specific tissues.

Intrachromosomal Looping and Histone K27 Methylation Coordinately Regulates the lncRNA H19-Fetal Mitogen IGF2 Imprinting Cluster in the Decidual Microenvironment of Early Pregnancy

Recurrent spontaneous abortion (RSA) is a highly heterogeneous complication of pregnancy with the underlying mechanisms remaining uncharacterized. Dysregulated decidualization is a critical contributor to the phenotypic alterations related to pregnancy complications. To understand the molecular factors underlying RSA, we explored the role of longnoncoding RNAs (lncRNAs) in the decidual microenvironment where the crosstalk at the fetal–maternal interface occurs. By exploring RNA-seq data from RSA patients, we identified H19, a noncoding RNA that exhibits maternal monoallelic expression, as one of the most upregulated lncRNAs associated with RSA. The paternally expressed fetal mitogen IGF2, which is reciprocally coregulated with H19 within the same imprinting cluster, was also upregulated. Notably, both genes underwent loss of imprinting, as H19 and IGF2 were actively transcribed from both parental alleles in some decidual tissues. This loss of imprinting in decidual tissues was associated with the loss of the H3K27m3 repressive histone marker in the IGF2 promoter, CpG hypomethylation at the central CTCF binding site in the imprinting control center (ICR), and the loss of CTCF-mediated intrachromosomal looping. These data suggest that dysregulation of the H19/IGF2 imprinting pathway may be an important epigenetic factor in the decidual microenvironment related to poor decidualization.

Downregulation of nc886 contributes to prostate cancer cell invasion and TGFβ1-induced EMT

Epithelial-to-mesenchymal transition (EMT) activation is important in cancer progression and metastasis. Evidence indicates that nc886 is a representative Pol III gene that processes microRNA products via Dicer and further downregulates its target gene transforming growth factor- β1 (TGF-β1), which is the most prominent inducer of EMT in prostate cancer (PC). Consistent with the previous literature, we found that nc886 downregulation was strongly associated with metastatic behavior and showed worse outcomes in PC patients. However, little is known about the association between nc886 and the EMT signaling pathway. We developed a PC cell model with stable overexpression of nc886 and found that nc886 changed cellular morphology and drove MET. The underlying mechanism may be related to its promotion of SNAIL protein degradation via ubiquitination, but not to its neighboring genes, TGFβ-induced protein (TGFBI) and SMAD5, which are Pol II-transcribed. TGF-β1 also override nc886 promotion of MET via transient suppression the transcription of nc886, promotion of TGFBI or increase in SMAD5 phosphorylation. Both nc886 inhibition and TGFBI activation occur regardless of their methylation status. The literature suggests that MYC inhibition by TGF-β1 is attributed to nc886 downregulation. We incidentally identified MYC-associated zinc finger protein (MAZ) as a suppressive transcription factor of TGFBI, which is controlled by TGF-β1. We elucidate a new mechanism of TGF-β1 differential control of Pol II and the transcription of its neighboring Pol III gene and identify a new EMT unit consisting of nc886 and its neighboring genes.

Methylation pattern of polymorphically imprinted nc886 is not conserved across mammalia

Background

In humans, the nc886 locus is a polymorphically imprinted metastable epiallele. Periconceptional conditions have an effect on the methylation status of nc886, and further, this methylation status is associated with health outcomes in later life, in line with the Developmental Origins of Health and Disease (DOHaD) hypothesis. Animal models would offer opportunities to study the associations between periconceptional conditions, nc886 methylation status and metabolic phenotypes further. Thus, we set out to investigate the methylation pattern of the nc886 locus in non-human mammals.

Data

We obtained DNA methylation data from the data repository GEO for mammals, whose nc886 gene included all three major parts of nc886 and had sequency similarity of over 80% with the human nc886. Our final sample set consisted of DNA methylation data from humans, chimpanzees, bonobos, gorillas, orangutangs, baboons, macaques, vervets, marmosets and guinea pigs.

Results

In human data sets the methylation pattern of nc886 locus followed the expected bimodal distribution, indicative of polymorphic imprinting. In great apes, we identified a unimodal DNA methylation pattern with 50% methylation level in all individuals and in all subspecies. In Old World monkeys, the between individual variation was greater and methylation on average was close to 60%. In guinea pigs the region around the nc886 homologue was non-methylated. Results obtained from the sequence comparison of the CTCF binding sites flanking the nc886 gene support the results on the DNA methylation data.

Conclusions

Our results indicate that unlike in humans, nc886 is not a polymorphically imprinted metastable epiallele in non-human primates or in guinea pigs, thus implying that animal models are not applicable for nc886 research. The obtained data suggests that the nc886 region may be classically imprinted in great apes, and potentially also in Old World monkeys, but not in guinea pigs.

Genome-wide DNA methylation profiling in nonagenarians suggests an effect of PM20D1 in late onset Alzheimer's disease

Background

The aim of this study is to identify differentially methylated regions (DMRs) in the genomes of a sample of cognitively healthy individuals and a sample of individuals with LOAD, all of them nonagenarians from Costa Rica.

Methods

In this study, we compared whole blood DNA methylation profiles of 32 individuals: 21 cognitively healthy and 11 with LOAD, using the Infinium MethylationEPIC BeadChip. First, we calculated the epigenetic age of the participants based on Horvath’s epigenetic clock. DMRcate and Bumphunter were used to identify DMRs. After in silico and knowledge-based filtering of the DMRs, we performed a methylation quantitative loci (mQTL) analysis (rs708727 and rs960603).

Results

On average, the epigenetic age was 73 years in both groups, which represents a difference of over 20 years between epigenetic and chronological age in both affected and unaffected individuals. Methylation analysis revealed 11 DMRs between groups, which contain six genes and two pseudogenes. These genes are involved in cell cycle regulation, embryogenesis, synthesis of ceramides, and migration of interneurons to the cerebral cortex. One of the six genes is PM20D1, for which altered expression has been reported in LOAD. After genotyping previously reported mQTL SNPs for the gene, we found that average methylation in the PM20D1 DMR differs between genotypes for rs708727, but not for rs960603.

Conclusions

This work supports the possible role of PM20D1 in protection against AD, by showing differential methylation in blood of affected and unaffected nonagenarians. Our results also support the influence of genetic factors on PM20D1 methylation levels.

The nuclear and cytoplasmic activities of RNA polymerase III, and an evolving transcriptome for surveillance

A 1969 report that described biochemical and activity properties of the three eukaryotic RNA polymerases revealed Pol III as highly distinguishable, even before its transcripts were identified. Now known to be the most complex, Pol III contains several stably-associated subunits referred to as built-in transcription factors (BITFs) that enable highly efficient RNA synthesis by a unique termination-associated recycling process. In vertebrates, subunit RPC7(α/β) can be of two forms, encoded by POLR3G or POLR3GL, with differential activity. Here we review promoter-dependent transcription by Pol III as an evolutionary perspective of eukaryotic tRNA expression. Pol III also provides nonconventional functions reportedly by promoter-independent transcription, one of which is RNA synthesis from DNA 3'-ends during repair. Another is synthesis of 5'ppp-RNA signaling molecules from cytoplasmic viral DNA in a pathway of interferon activation that is dysfunctional in immunocompromised patients with mutations in Pol III subunits. These unconventional functions are also reviewed, including evidence that link them to the BITF subunits. We also review data on a fraction of the human Pol III transcriptome that evolved to include vault RNAs and snaRs with activities related to differentiation, and in innate immune and tumor surveillance. The Pol III of higher eukaryotes does considerably more than housekeeping.

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.

VTRNA2-1: Genetic Variation, Heritable Methylation and Disease Association

VTRNA2-1 is a metastable epiallele with accumulating evidence that methylation at this region is heritable, modifiable and associated with disease including risk and progression of cancer. This study investigated the influence of genetic variation and other factors such as age and adult lifestyle on blood DNA methylation in this region. We first sequenced the VTRNA2-1 gene region in multiple-case breast cancer families in which VTRNA2-1 methylation was identified as heritable and associated with breast cancer risk. Methylation quantitative trait loci (mQTL) were investigated using a prospective cohort study (4500 participants with genotyping and methylation data). The cis-mQTL analysis (334 variants ± 50 kb of the most heritable CpG site) identified 43 variants associated with VTRNA2-1 methylation (p < 1.5 × 10⁻⁴); however, these explained little of the methylation variation (R² < 0.5% for each of these variants). No genetic variants elsewhere in the genome were found to strongly influence VTRNA2-1 methylation. SNP-based heritability estimates were consistent with the mQTL findings (h² = 0, 95%CI: −0.14 to 0.14). We found no evidence that age, sex, country of birth, smoking, body mass index, alcohol consumption or diet influenced blood DNA methylation at VTRNA2-1. Genetic factors and adult lifestyle play a minimal role in explaining methylation variability at the heritable VTRNA2-1 cluster.

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.

Differential hypermethylation of the VTRNA2-1 promoter in hepatocellular carcinoma as a prognostic factor: Tumor marker prognostic study

BACKGROUND

Vault RNA 2-1 (VTRNA2-1, also called nc886) is a 108-nucleotide noncoding transcript that is epigenetically controlled via 18 CpG dinucleotide modifications of its promoter, and can exert either tumor suppressor or oncogenic functions depending on cell types of cancers. In hepatocellular carcinoma (HCC), the role of VTRNA2-1 in prognosis of patients remains unexplored. Here, we analysed the methylation status of the VTRNA2-1 promoter and its correlation with clinical parameters in patients with HCC.

PATIENTS AND METHODS

A total of 92 patients with HCC were enrolled, genomic DNA of tumor versus normal tissues were extracted and bisulfite modified. VTRNA2-1 promoter regions chr5: 135416381 (cg06536614), 135416388, 135416394 (cg26328633), and 135416398 (cg25340688) were PCR amplified and pyrosequenced. The methylation status of VTRNA2-1 in patients was further analysed with other clinical parameters via univariate and multivariate analysis.

RESULTS

The differential hypermethylation status (tumor- normal) of the VTRNA2-1 promoter in HCC correlated well with the presence of large tumor size (p = 0.001), pathological vascular invasion (p = 0.036), tumor recurrence (p = 0.007) and more advanced tumor stage (stage III AJCC) in patients (p = 0.03). In addition, the methylation of the VTRNA2-1 promoter increased in stage III HCC tumor compared with stage I & II tumor (64.7% versus 36.0%, p = 0.030). Furthermore, the differential hypermethylation status of the VTRNA2-1 promoter was an independent factor for patient outcome after partial hepatectomy using multivariate Cox regression analysis (p = 0.011, HR = 2.305). Using another public dataset (GSE89852), we found that the differential hypermethylation of the VTRNA2-1 promoter was also significantly associated with tumor recurrence.

CONCLUSIONS

Patients had unfavourable outcomes when the VTRNA2-1 promoter was differentially hypermethylated in tumor tissues compared to its adjacent normal tissues. These findings suggest that such patients should receive intensive follow-up care or possible adjuvant therapy after liver resection.

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.

Insights into imprinting from parent-of-origin phased methylomes and transcriptomes

Imprinting is the preferential expression of one parental allele over the other. It is controlled primarily through differential methylation of cytosine at CpG dinucleotides. Here we combine 285 methylomes and 11,617 transcriptomes from peripheral blood samples with parent-of-origin phased haplotypes, to produce a new map of imprinted methylation and gene expression patterns across the human genome. We demonstrate how imprinted methylation is a continuous rather than a binary characteristic. We describe at high resolution the parent-of-origin methylation pattern at the 15q11.2 Prader-Willi/Angelman syndrome locus, with nearly confluent stochastic paternal methylation punctuated by 'spikes' of maternal methylation. We find examples of polymorphic imprinted methylation unrelated (at VTRNA2-1 and PARD6G) or related (at CHRNE) to nearby SNP genotypes. We observe RNA isoform-specific imprinted expression patterns suggestive of a methylation-sensitive transcriptional elongation block. Finally, we gain new insights into parent-of-origin-specific effects on phenotypes at the DLK1/MEG3 and GNAS loci.

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.

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.

Epigenetic supersimilarity of monozygotic twin pairs

BACKGROUND

Monozygotic twins have long been studied to estimate heritability and explore epigenetic influences on phenotypic variation. The phenotypic and epigenetic similarities of monozygotic twins have been assumed to be largely due to their genetic identity.

RESULTS

Here, by analyzing data from a genome-scale study of DNA methylation in monozygotic and dizygotic twins, we identified genomic regions at which the epigenetic similarity of monozygotic twins is substantially greater than can be explained by their genetic identity. This "epigenetic supersimilarity" apparently results from locus-specific establishment of epigenotype prior to embryo cleavage during twinning. Epigenetically supersimilar loci exhibit systemic interindividual epigenetic variation and plasticity to periconceptional environment and are enriched in sub-telomeric regions. In case-control studies nested in a prospective cohort, blood DNA methylation at these loci years before diagnosis is associated with risk of developing several types of cancer.

CONCLUSIONS

These results establish a link between early embryonic epigenetic development and adult disease. More broadly, epigenetic supersimilarity is a previously unrecognized phenomenon that may contribute to the phenotypic similarity of monozygotic twins.

DNA methylation in blood from neonatal screening cards and the association with BMI and insulin sensitivity in early childhood

BACKGROUND/OBJECTIVES

There is increasing evidence that metabolic diseases originate in early life, and epigenetic changes have been implicated as key drivers of this early life programming. This led to the hypothesis that epigenetic marks present at birth may predict an individual's future risk of obesity and type 2 diabetes. In this study, we assessed whether epigenetic marks in blood of newborn children were associated with body mass index (BMI) and insulin sensitivity later in childhood.

SUBJECTS/METHODS

DNA methylation was measured in neonatal blood spot samples of 438 children using the Illumina Infinium 450 k BeadChip. Associations were assessed between DNA methylation at birth and BMI z-scores, body fat mass, fasting plasma glucose, insulin and homeostatic model assessment of insulin resistance (HOMA-IR) at age 5 years, as well as birth weight, maternal BMI and smoking status.

RESULTS

No individual methylation sites at birth were associated with obesity or insulin sensitivity measures at 5 years. DNA methylation in 69 genomic regions at birth was associated with BMI z-scores at age 5 years, and in 63 regions with HOMA-IR. The methylation changes were generally small (<5%), except for a region near the non-coding RNA nc886 (VTRNA2-1) where a clear link between methylation status at birth and BMI in childhood was observed (P=0.001). Associations were also found between DNA methylation, maternal smoking and birth weight.

CONCLUSIONS

We identified a number of DNA methylation regions at birth that were associated with obesity or insulin sensitivity measurements in childhood. These findings support the mounting evidence on the role of epigenetics in programming of metabolic health. Whether many of these small changes in DNA methylation are causally related to the health outcomes, and of clinical relevance, remains to be determined, but the nc886 region represents a promising obesity risk marker that warrants further investigation.

Epigenetic regulation of noncoding RNA transcription by mammalian RNA polymerase III

RNA polymerase III (Pol III) synthesizes a range of medium-sized noncoding RNAs (collectively 'Pol III genes') whose early established biological roles were so essential that they were considered 'housekeeping genes'. Besides these fundamental functions, diverse unconventional roles of mammalian Pol III genes have recently been recognized and their expression must be exquisitely controlled. In this review, we summarize the epigenetic regulation of Pol III genes by chromatin structure, histone modification and CpG DNA methylation. We also recapitulate the association between dysregulation of Pol III genes and diseases such as cancer and neurological disorders. Additionally, we will discuss why in-depth molecular studies of Pol III genes have not been attempted and how nc886, a Pol III gene, may resolve this issue.

Epigenetic regulation of RNA polymerase III transcription in early breast tumorigenesis

RNA polymerase III (Pol III) transcribes medium-sized non-coding RNAs (collectively termed Pol III genes). Emerging diverse roles of Pol III genes suggest that individual Pol III genes are exquisitely regulated by transcription and epigenetic factors. Here we report global Pol III expression/methylation profiles and molecular mechanisms of Pol III regulation that have not been as extensively studied, using nc886 as a representative Pol III gene. In a human mammary epithelial cell system that recapitulates early breast tumorigenesis, the fraction of actively transcribed Pol III genes increases reaching a plateau during immortalization. Hyper-methylation of Pol III genes inhibits Pol III binding to DNA via inducing repressed chromatin and is a determinant for the Pol III repertoire. When Pol III genes are hypo-methylated, MYC amplifies their transcription, regardless of its recognition DNA motif. Thus, Pol III expression during tumorigenesis is delineated by methylation and magnified by MYC.

Human noncoding RNA 886 (nc886) adopts two structurally distinct conformers that are functionally opposing regulators of PKR

The double-stranded RNA (dsRNA)-activated protein kinase (PKR) senses dsRNA produced during viral infection and halts cellular protein synthesis to block viral replication. How basal PKR activity is controlled in the absence of infection was unclear until the recent identification of a potential endogenous regulator, the cellular noncoding RNA 886 (nc886). However, nc886 adopts two distinct conformations for which the structural details and potential functional differences remain unclear. Here, we isolated and separately dissected the function of each form of nc886 to more clearly define the molecular mechanism of nc886-mediated PKR inhibition. We show that nc886 adopts two stable, noninterconverting RNA conformers that are functionally nonequivalent using complementary RNA structure probing and mutational analyses combined with PKR binding and activity assays. One conformer acts as a potent inhibitor, while the other is a pseudoinhibitor capable of weakly activating the kinase. We mapped the nc886 region necessary for high affinity binding and potent inhibition of PKR to an apical stem-loop structure present in only one conformer of the RNA. This structural feature is not only critical for inhibiting PKR autophosphorylation, but also the phosphorylation of its cellular substrate, the eukaryotic translation initiation factor 2α subunit. The identification of different activities of the nc886 conformers suggests a potential mechanism for producing a gradient of PKR regulation within the cell and reveals a way by which a cellular noncoding RNA can mask or present a structural feature to PKR for inhibition.

Human Oocyte-Derived Methylation Differences Persist in the Placenta Revealing Widespread Transient Imprinting

Thousands of regions in gametes have opposing methylation profiles that are largely resolved during the post-fertilization epigenetic reprogramming. However some specific sequences associated with imprinted loci survive this demethylation process. Here we present the data describing the fate of germline-derived methylation in humans. With the exception of a few known paternally methylated germline differentially methylated regions (DMRs) associated with known imprinted domains, we demonstrate that sperm-derived methylation is reprogrammed by the blastocyst stage of development. In contrast a large number of oocyte-derived methylation differences survive to the blastocyst stage and uniquely persist as transiently methylated DMRs only in the placenta. Furthermore, we demonstrate that this phenomenon is exclusive to primates, since no placenta-specific maternal methylation was observed in mouse. Utilizing single cell RNA-seq datasets from human preimplantation embryos we show that following embryonic genome activation the maternally methylated transient DMRs can orchestrate imprinted expression. However despite showing widespread imprinted expression of genes in placenta, allele-specific transcriptional profiling revealed that not all placenta-specific DMRs coordinate imprinted expression and that this maternal methylation may be absent in a minority of samples, suggestive of polymorphic imprinted methylation.

Differences in gamete DNA methylation is subject to genome-wide reprogramming during preimplantation development to establish an embryo with an epigenetic state compatible with totipotency. DNA sequences associated with imprinted differentially methylated regions (DMRs) are largely protected from this process, retaining their parent-of-origin epigenetic marks. By comparing the methylation profiles of human oocytes, sperm, blastocysts and various somatic tissues including placenta, we observe hundreds of CpG island sequences that maintain methylation on their maternal allele in blastocysts and placenta indicative of incomplete reprogramming. In some cases this maternal methylation influence transcription of nearby genes, revealing transient imprinting in embryos after genome-activation and in placenta. Strikingly, these placenta-specific DMRs are polymorphic between placenta samples with a minority of samples being robustly unmethylated on both alleles.

Whole-exome sequencing and genome-wide methylation analyses identify novel disease associated mutations and methylation patterns in idiopathic hypereosinophilic syndrome

A thorough understanding of the idiopathic hypereosinophilic syndrome (IHES) and further optimization of diagnostic work-up procedures are warranted. We analyzed purified eosinophils from patients with IHES by next-generation whole-exome sequencing and compared DNA methylation profiles from reactive eosinophilic conditions to known clonal and suspected clonal eosinophilia. Somatic missense mutations in cancer-related genes were detected in three IHES patients. These included the spliceosome gene PUF60 and the cadherin gene CDH17. Furthermore, reactive eosinophilia samples could be differentiated from known- and suspected clonal eosinophilia samples based on 285 differentially methylated CpG sites corresponding to 128 differentially methylated genes. Using Ingenuity pathway analysis, we found that differentially methylated genes were highly enriched in functional pathways such as cancer, cell death and survival, and hematological disease. Our data show that a subset of IHES may be of clonal origin not related to the classical molecular aberrations of FGFR, PDGFRA/B, or T-cells, and that the initiating hits could be point mutations in a variety of genes, including spliceosome mutations or hypermethylated tumor suppressor genes. In addition, we identified a DNA methylation signature that is relevant for distinguishing clonal and suspected clonal eosinophilia from reactive eosinophilia per se, which may be useful in daily clinical work.

DNA Methylation Profiling of Uniparental Disomy Subjects Provides a Map of Parental Epigenetic Bias in the Human Genome

Genomic imprinting is a mechanism in which gene expression varies depending on parental origin. Imprinting occurs through differential epigenetic marks on the two parental alleles, with most imprinted loci marked by the presence of differentially methylated regions (DMRs). To identify sites of parental epigenetic bias, here we have profiled DNA methylation patterns in a cohort of 57 individuals with uniparental disomy (UPD) for 19 different chromosomes, defining imprinted DMRs as sites where the maternal and paternal methylation levels diverge significantly from the biparental mean. Using this approach we identified 77 DMRs, including nearly all those described in previous studies, in addition to 34 DMRs not previously reported. These include a DMR at TUBGCP5 within the recurrent 15q11.2 microdeletion region, suggesting potential parent-of-origin effects associated with this genomic disorder. We also observed a modest parental bias in DNA methylation levels at every CpG analyzed across ∼1.9 Mb of the 15q11-q13 Prader-Willi/Angelman syndrome region, demonstrating that the influence of imprinting is not limited to individual regulatory elements such as CpG islands, but can extend across entire chromosomal domains. Using RNA-seq data, we detected signatures consistent with imprinted expression associated with nine novel DMRs. Finally, using a population sample of 4,004 blood methylomes, we define patterns of epigenetic variation at DMRs, identifying rare individuals with global gain or loss of methylation across multiple imprinted loci. Our data provide a detailed map of parental epigenetic bias in the human genome, providing insights into potential parent-of-origin effects.

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.

Hypermethylation of the VTRNA1-3 Promoter is Associated with Poor Outcome in Lower Risk Myelodysplastic Syndrome Patients

Myelodysplastic syndrome (MDS) is a heterogeneous group of clonal hematopoietic disorders. MDS is frequently associated with deletions on chromosome 5q as well as aberrant DNA methylation patterns including hypermethylation of key tumor suppressors. We have previously shown that hypermethylation and silencing of the non-coding RNA VTRNA2-1 are correlated with poor outcomes in acute myeloid leukemia patients. In this study, we find that VTRNA1-2 and VTRNA1-3, both located on chromosome 5q, can be regulated and silenced by promoter DNA methylation, and that the hypomethylating agent 5-aza-2-deoxycytidine causes reactivation these genes. In normal hematopoiesis, we find that vault RNAs (vtRNAs) show differential methylation between various hematopoietic cell populations, indicating that allele-specific methylation events may occur during hematopoiesis. In addition, we show that VTRNA1-3 promoter hypermethylation is frequent in lower risk MDS patients and is associated with a decreased overall survival.

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.