Allelic methylation levels of the noncoding VTRNA2-1 located on chromosome 5q31.1 predict outcome in AML

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.

Altered methylation of imprinted genes in neuroblastoma: implications for prognostic refinement

BACKGROUND

Neuroblastoma (NB) is a complex disease, and the current understanding of NB biology is limited. Deregulation in genomic imprinting is a common event in malignancy. Since imprinted genes play crucial roles in early fetal growth and development, their role in NB pathogenesis could be suggested.

METHODS

We examined alterations in DNA methylation patterns of 369 NB tumours at 49 imprinted differentially methylated regions (DMRs) and assessed its association with overall survival probabilities and selected clinical and genomic features of the tumours. In addition, an integrated analysis of DNA methylation and allele-specific copy number alterations (CNAs) was performed, to understand the correlation between the two molecular events.

RESULTS

Several imprinted regions with aberrant methylation patterns in NB were identified. Regions that underwent loss of methylation in > 30% of NB samples were DMRs annotated to the genes NDN, SNRPN, IGF2, MAGEL2 and HTR5A and regions with gain of methylation were NNAT, RB1 and GPR1. Methylation alterations at six of the 49 imprinted DMRs were statistically significantly associated with reduced overall survival: MIR886, RB1, NNAT/BLCAP, MAGEL2, MKRN3 and INPP5F. RB1, NNAT/BLCAP and MKRN3 were further able to stratify low-risk NB tumours i.e. tumours that lacked MYCN amplification and 11q deletion into risk groups. Methylation alterations at NNAT/BLCAP, MAGEL2 and MIR886 predicted risk independently of MYCN amplification or 11q deletion and age at diagnosis. Investigation of the allele-specific CNAs demonstrated that the imprinted regions that displayed most alterations in NB tumours harbor true epigenetic changes and are not result of the underlying CNAs.

CONCLUSIONS

Aberrant methylation in imprinted regions is frequently occurring in NB tumours and several of these regions have independent prognostic value. Thus, these could serve as potentially important clinical epigenetic markers to identify individuals with adverse prognosis. Incorporation of methylation status of these regions together with the established risk predictors may further refine the prognostication of NB patients.

Human Vault RNAs: Exploring Their Potential Role in Cellular Metabolism

Non-coding RNAs have been described as crucial regulators of gene expression and guards of cellular homeostasis. Some recent papers focused on vault RNAs, one of the classes of non-coding RNA, and their role in cell proliferation, tumorigenesis, apoptosis, cancer response to therapy, and autophagy, which makes them potential therapy targets in oncology. In the human genome, four vault RNA paralogues can be distinguished. They are associated with vault complexes, considered the largest ribonucleoprotein complexes. The protein part of these complexes consists of a major vault protein (MVP) and two minor vault proteins (vPARP and TEP1). The name of the complex, as well as vault RNA, comes from the hollow barrel-shaped structure that resembles a vault. Their sequence and structure are highly evolutionarily conserved and show many similarities in comparison with different species, but vault RNAs have various roles. Vaults were discovered in 1986, and their functions remained unclear for many years. Although not much is known about their contribution to cell metabolism, it has become clear that vault RNAs are involved in various processes and pathways associated with cancer progression and modulating cell functioning in normal and pathological stages. In this review, we discuss known functions of human vault RNAs in the context of cellular metabolism, emphasizing processes related to cancer and cancer therapy efficacy.

Crosstalk between vault RNAs and innate immunity

PURPOSE

Vault (vt) RNAs are noncoding (nc) RNAs transcribed by RNA polymerase III (RNA Pol III) with 5'-triphosphate (5'-PPP) termini that play significant roles and are recognized by innate immune sensors, including retinoic acid-inducible protein 1 (RIG-I). In addition, vtRNAs adopt secondary structures that can be targets of interferon-inducible protein kinase R (PKR) and the oligoadenylate synthetase (OAS)/RNase L system, both of which are important for activating antiviral defenses. However, changes in the expression of vtRNAs have been associated with pathological processes that activate proinflammatory pathways, which influence cellular events such as differentiation, aging, autophagy, apoptosis, and drug resistance in cancer cells.

RESULTS

In this review, we summarized the biology of vtRNAs and focused on their interactions with the innate immune system. These findings provide insights into the diverse roles of vtRNAs and their correlation with various cellular processes to improve our understanding of their biological functions.

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.

Metastable epialleles in humans

First identified in isogenic mice, metastable epialleles (MEs) are loci where the extent of DNA methylation (DNAm) is variable between individuals but correlates across tissues derived from different germ layers within a given individual. This property, termed systemic interindividual variation (SIV), is attributed to stochastic methylation establishment before germ layer differentiation. Evidence suggests that some putative human MEs are sensitive to environmental exposures in early development. In this review we introduce key concepts pertaining to human MEs, describe methods used to identify MEs in humans, and review their genomic features. We also highlight studies linking DNAm at putative human MEs to early environmental exposures and postnatal (including disease) phenotypes.

The Pol III transcriptome: Basic features, recurrent patterns, and emerging roles in cancer

The RNA polymerase III (Pol III) transcriptome is universally comprised of short, highly structured noncoding RNA (ncRNA). Through RNA-protein interactions, the Pol III transcriptome actuates functional activities ranging from nuclear gene regulation (7SK), splicing (U6, U6atac), and RNA maturation and stability (RMRP, RPPH1, Y RNA), to cytoplasmic protein targeting (7SL) and translation (tRNA, 5S rRNA). In higher eukaryotes, the Pol III transcriptome has expanded to include additional, recently evolved ncRNA species that effectively broaden the footprint of Pol III transcription to additional cellular activities. Newly evolved ncRNAs function as riboregulators of autophagy (vault), immune signaling cascades (nc886), and translation (Alu, BC200, snaR). Notably, upregulation of Pol III transcription is frequently observed in cancer, and multiple ncRNA species are linked to both cancer progression and poor survival outcomes among cancer patients. In this review, we outline the basic features and functions of the Pol III transcriptome, and the evidence for dysregulation and dysfunction for each ncRNA in cancer. When taken together, recurrent patterns emerge, ranging from shared functional motifs that include molecular scaffolding and protein sequestration, overlapping protein interactions, and immunostimulatory activities, to the biogenesis of analogous small RNA fragments and noncanonical miRNAs, augmenting the function of the Pol III transcriptome and further broadening its role in cancer. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Processing > Processing of Small RNAs RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications.

DNA strand asymmetry generated by CpG hemimethylation has opposing effects on CTCF binding

CpG methylation generally occurs on both DNA strands and is essential for mammalian development and differentiation. Until recently, hemimethylation, in which only one strand is methylated, was considered to be simply a transitory state generated during DNA synthesis. The discovery that a subset of CCCTC-binding factor (CTCF) binding sites is heritably hemimethylated suggests that hemimethylation might have an unknown biological function. Here we show that the binding of CTCF is profoundly altered by which DNA strand is methylated and by the specific CTCF binding motif. CpG methylation on the motif strand can inhibit CTCF binding by up to 7-fold, whereas methylation on the opposite strand can stimulate binding by up to 4-fold. Thus, hemimethylation can alter binding by up to 28-fold in a strand-specific manner. The mechanism for sensing methylation on the opposite strand requires two critical residues, V454 and S364, within CTCF zinc fingers 7 and 4. Similar to methylation, CpG hydroxymethylation on the motif strand can inhibit CTCF binding by up to 4-fold. However, hydroxymethylation on the opposite strand removes the stimulatory effect. Strand-specific methylation states may therefore provide a mechanism to explain the transient and dynamic nature of CTCF-mediated chromatin interactions.

nc886, an RNA Polymerase III-Transcribed Noncoding RNA Whose Expression Is Dynamic and Regulated by Intriguing Mechanisms

nc886 is a medium-sized non-coding RNA that is transcribed by RNA polymerase III (Pol III) and plays diverse roles in tumorigenesis, innate immunity, and other cellular processes. Although Pol III-transcribed ncRNAs were previously thought to be expressed constitutively, this concept is evolving, and nc886 is the most notable example. The transcription of nc886 in a cell, as well as in human individuals, is controlled by multiple mechanisms, including its promoter CpG DNA methylation and transcription factor activity. Additionally, the RNA instability of nc886 contributes to its highly variable steady-state expression levels in a given situation. This comprehensive review discusses nc886's variable expression in physiological and pathological conditions and critically examines the regulatory factors that determine its expression levels.

Epigenetic silencing of ZCCHC10 by the lncRNA SNHG1 promotes progression and venetoclax resistance of acute myeloid leukemia

The gene encoding the tumor suppressor p53 is the most frequently mutated gene in cancers. However, p53 mutation is rare in acute myeloid leukemia (AML), and p53 is inactivated predominantly by aberrant expression of p53 regulators (such as MDM2). A previous study by the authors revealed that the ZCCHC10 protein suppressed MDM2‑mediated degradation of the p53 protein in lung cancer. However, the expression and role of the ZCCHC10 gene in AML have not been investigated. In the present study, it was found that ZCCHC10 expression was downregulated in bone marrow samples of AML patients and that ZCCHC10 expression was significantly and negatively correlated with the expression of the lncRNA SNHG1. Suppression of SNHG1 decreased ZCCHC10 promoter methylation and increased ZCCHC10 expression. Notably, there is a putative binding motif in SNHG1 with full complementarity to five sites surrounding the CpG island in the ZCCHC10 promoter. Overexpression of wild‑type SNHG1 promoted ZCCHC10 methylation, but overexpression of SNHG1 with deletion of the binding motif did not. Further study identified that SNHG1 simultaneously bound to the ZCCHC10 promoter and the DNA methyltransferases DNMT1 and DNMT3B. These results indicated that SNHG1 recruits DNMT1 and DNMT3B to the ZCCHC10 promoter, resulting in hypermethylation of the ZCCHC10 promoter. Kaplan‑Meier survival analysis showed that ZCCHC10 expression was positively associated with overall survival in AML patients. In vitro experiments demonstrated that ZCCHC10 increased p53 expression and suppressed AML cell proliferation and survival. In the xenograft mouse model, ZCCHC10 decreased the proliferation of leukemic cells, improved the survival of leukemic mice, and increased sensitivity to the BCL inhibitor venetoclax. In conclusion, ZCCHC10 expression is suppressed by SNHG1‑induced DNA methylation in AML. Downregulation of ZCCHC10 decreases p53 activation, promotes cell proliferation and survival, and thereby accelerates AML progression and the acquisition of venetoclax resistance. The present study identified a SNHG1/ZCCHC10/p53 signaling axis in AML that may be a therapeutic target in this malignancy.

Small noncoding vault RNA2-1 disrupts gut epithelial barrier function via interaction with HuR

Vault RNAs (vtRNAs) are small noncoding RNAs and highly expressed in many eukaryotes. Here, we identified vtRNA2-1 as a novel regulator of the intestinal barrier via interaction with RNA-binding protein HuR. Intestinal mucosal tissues from patients with inflammatory bowel diseases and from mice with colitis or sepsis express increased levels of vtRNAs relative to controls. Ectopically expressed vtRNA2-1 decreases the levels of intercellular junction (IJ) proteins claudin 1, occludin, and E-cadherin and causes intestinal epithelial barrier dysfunction in vitro, whereas vtRNA2-1 silencing promotes barrier function. Increased vtRNA2-1 also decreases IJs in intestinal organoid, inhibits epithelial renewal, and causes Paneth cell defects ex vivo. Elevating the levels of tissue vtRNA2-1 in the intestinal mucosa increases the vulnerability of the gut barrier to septic stress in mice. vtRNA2-1 interacts with HuR and prevents HuR binding to claudin 1 and occludin mRNAs, thus decreasing their translation. These results indicate that vtRNA2-1 impairs intestinal barrier function by repressing HuR-facilitated translation of claudin 1 and occludin.

Hypermethylation of nc886 in HPV-positive oropharyngeal cancer and its clinical implications: An epigenome-wide association study

The incidence of oropharyngeal squamous cell carcinoma (OPSCC) has increased rapidly in the United States, driven by rising human papillomavirus (HPV) infections in the U.S. population. HPV-positive OPSCC patients have a better prognosis than HPV-negative patients. To gain insights into the unique biology of HPV(+) OPSCC that may contribute to its clinical behaviors, we performed a multi-stage epigenome-wide methylation profiling of leukocyte and tumor DNA in OPSCC patients and compared the methylation levels of CpG sites between HPV(+) and HPV(-) OPSCC patients. We identified and validated a significantly differentially methylated region (DMR) of 1,355 bp encompassing non-coding RNA 886 (nc886) gene and its promoter region. Nc886 is hypermethylated in both leukocytes and tumor DNA of HPV(+) OPSCC patients. Homozygous knockout of nc886 by CRISPR-Cas9 in head and neck cell lines was lethal, but nc886 could be knocked out on the background of protein kinase R (PKR) knockout. Our data suggest that HPV induces nc886 hypermethylation, and nc886 acts as both a viral sensor and a tumor sensor in OPSCC patients and contribute to the better prognosis of HPV(+) OPSCC patients. Nc886 may become a therapeutic target in OPSCC.

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.

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.

Small but Powerful: The Human Vault RNAs as Multifaceted Modulators of Pro-Survival Characteristics and Tumorigenesis

Simple Summary

Small non-protein-coding RNAs have been recognized as valuable regulators of gene expression in all three domains of life. Particularly in multicellular organisms, ncRNAs-mediated gene expression control has evolved as a central principle of cellular homeostasis. Thus, it is not surprising that non-coding RNA misregulation has been linked to various diseases. Here, we review the contributions of the four human vault RNAs to cellular proliferation, apoptosis and cancer biology.

Abstract

The importance of non-coding RNAs for regulating gene expression has been uncovered in model systems spanning all three domains of life. More recently, their involvement in modulating signal transduction, cell proliferation, tumorigenesis and cancer progression has also made them promising tools and targets for oncotherapy. Recent studies revealed a class of highly conserved small ncRNAs, namely vault RNAs, as regulators of several cellular homeostasis mechanisms. The human genome encodes four vault RNA paralogs that share significant sequence and structural similarities, yet they seem to possess distinct roles in mammalian cells. The alteration of vault RNA expression levels has frequently been observed in cancer tissues, thus hinting at a putative role in orchestrating pro-survival characteristics. Over the last decade, significant advances have been achieved in clarifying the relationship between vault RNA and cellular mechanisms involved in cancer development. It became increasingly clear that vault RNAs are involved in controlling apoptosis, lysosome biogenesis and function, as well as autophagy in several malignant cell lines, most likely by modulating signaling pathways (e.g., the pro-survival MAPK cascade). In this review, we discuss the identified and known functions of the human vault RNAs in the context of cell proliferation, tumorigenesis and chemotherapy resistance.

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.

Are We Studying Non-Coding RNAs Correctly? Lessons from nc886

Non-coding RNAs (ncRNAs), such as microRNAs or long ncRNAs, have brought about a new paradigm in the regulation of gene expression. Sequencing technologies have detected transcripts with tremendous sensitivity and throughput and revealed that the majority of them lack protein-coding potential. Myriad articles have investigated numerous ncRNAs and many of them claim that ncRNAs play gene-regulatory roles. However, it is questionable whether all these articles draw conclusions through cautious gain- and loss-of function experiments whose design was reasonably based on an ncRNA's correct identity and features. In this review, these issues are discussed with a regulatory ncRNA, nc886, as an example case to represent cautions and guidelines when studying ncRNAs.

A validation study of potential prognostic DNA methylation biomarkers in patients with acute myeloid leukemia using a custom DNA methylation sequencing panel

BACKGROUND

Multiple studies have reported the prognostic impact of DNA methylation changes in acute myeloid leukemia (AML). However, these epigenetic markers have not been thoroughly validated and therefore are still not considered in clinical practice. Hence, we aimed to independently verify results of selected studies describing the relationship between DNA methylation of specific genes and their prognostic potential in predicting overall survival (OS) and event-free survival (EFS).

RESULTS

Fourteen studies (published 2011-2019) comprising of 27 genes were subjected to validation by a custom NGS-based sequencing panel in 178 newly diagnosed non-M3 AML patients treated by 3 + 7 induction regimen. The results were considered as successfully validated, if both the log-rank test and multivariate Cox regression analysis had a p-value ≤ 0.05. The predictive role of DNA methylation was confirmed for three studies comprising of four genes: CEBPA (OS: p = 0.02; EFS: p = 0.03), PBX3 (EFS: p = 0.01), LZTS2 (OS: p = 0.05; EFS: p = 0.0003), and NR6A1 (OS: p = 0.004; EFS: p = 0.0003). For all of these genes, higher methylation was an indicator of longer survival. Concurrent higher methylation of both LZTS2 and NR6A1 was highly significant for survival in cytogenetically normal (CN) AML group (OS: p < 0.0001; EFS: p < 0.0001) as well as for the whole AML cohort (OS: p = 0.01; EFS < 0.0001). In contrast, for two studies reporting the poor prognostic effect of higher GPX3 and DLX4 methylation, we found the exact opposite, again linking higher GPX3 (OS: p = 0.006; EFS: p < 0.0001) and DLX4 (OS: p = 0.03; EFS = 0.03) methylation to a favorable treatment outcome. Individual gene significance levels refer to the outcomes of multivariate Cox regression analysis.

CONCLUSIONS

Out of twenty-seven genes subjected to DNA methylation validation, a prognostic role was observed for six genes. Therefore, independent validation studies are necessary to reveal truly prognostic DNA methylation changes and to enable the introduction of these promising epigenetic markers into clinical practice.

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.

Elevated methylation of the vault RNA2-1 promoter in maternal blood is associated with preterm birth

Background

Preterm birth, defined as parturition before 37 completed weeks of gestation, is associated with an increased risk of neonatal complications and death, as well as poor health and disease later in life. Epigenetics could contribute to the mechanism underlying preterm birth.

Results

Genome-wide DNA methylation analysis of whole blood cells from 10 women (5 term and 5 preterm deliveries) was performed using an Illumina Infinium HumanMethylation450 BeadChips array. We identified 1,581 differentially methylated CpG sites in promoter regions between term and preterm birth. Although the differences were not significant after correcting for multiple tests, seven CpGs on the genomically imprinted vault RNA2-1 (VTRNA2-1; also known as non-coding RNA, nc886 or miR-886) showed the largest differences (range: 26–39 %). Pyrosequencing verification was performed with blood samples from pregnant women recruited additionally (39 term and 43 preterm deliveries). In total, 28 (34.1 %) samples showed hypomethylation of the VTRNA2-1 promoter (< 13 % methylation), while 54 (65.9 %) samples showed elevated methylation levels between 30 and 60 %. Elevated methylation of VTRNA2-1 promoter was associated with an increased risk of preterm birth after adjusting for maternal age, season of delivery, parity and white blood cell count. The mRNA expression of VTRNA2-1 was 0.51-fold lower in women with preterm deliveries (n = 20) compared with women with term deliveries (n = 20).

Conclusions

VTRNA2-1 is a noncoding transcript to environmentally responsive epialleles. Our results suggest that elevated methylation of the VTRNA2-1 promoter may result in increased risk of PTB caused by the pro-inflammatory cytokines. Further studies are needed to confirm the association of VTRNA2-1 methylation with preterm birth in a large population, and to elucidate the underlying mechanism.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07865-y.

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.

Vault RNAs: hidden gems in RNA and protein regulation

Non-coding RNAs are important regulators of differentiation during embryogenesis as well as key players in the fine-tuning of transcription and furthermore, they control the post-transcriptional regulation of mRNAs under physiological conditions. Deregulated expression of non-coding RNAs is often identified as one major contribution in a number of pathological conditions. Non-coding RNAs are a heterogenous group of RNAs and they represent the majority of nuclear transcripts in eukaryotes. An evolutionary highly conserved sub-group of non-coding RNAs is represented by vault RNAs, named since firstly discovered as component of the largest known ribonucleoprotein complexes called "vault". Although they have been initially described 30 years ago, vault RNAs are largely unknown and their molecular role is still under investigation. In this review we will summarize the known functions of vault RNAs and their involvement in cellular mechanisms.

MiR-410-3p activates the NF-κB pathway by targeting ZCCHC10 to promote migration, invasion and EMT of colorectal cancer

BACKGROUND

Colorectal cancer (CRC) is a prevalent malignancy of the digestive tract. miR-410-3p is involved in oncogenesis and development of CRC, but the specific regulation mechanism is still not known clearly.

METHODS

The expression of miR-410-3p and zinc finger CCHC-type containing 10 (ZCCHC10) in CRC cells was detected by qRT-PCR and western blot method, respectively. The dual-luciferase reporter gene detection was applied for determination of interaction between miR-410-3p and ZCCHC10. The wound healing assay and transwell assay were carried out to measure cell migration and invasive ability, respectively.

RESULTS

The miR-410-3p expression levels were markedly increased, but ZCCHC10 levels were reduced in CRC cells and tissues. Dual-luciferase reporter gene detection indicated that miR-410-3p targeted ZCCHC10 directly. Functionally knockdown of ZCCHC10 or overexpression of miR-410-3p activated nuclear factor-κB (NF-κB) signaling pathway, promoted epithelial-mesenchymal transition (EMT) process, as well as cell migration and invasion of CRC cells. After adding NF-κB inhibitor BAY 11-708 to inhibit NF-κB pathway, the promoting effects of si-ZCCHC10 on cell migration, invasion and EMT of HT29 and SW480 cells were suppressed. Meanwhile, overexpression of ZCCHC10 inhibited the effects of miR-410-3p on cell migration, invasion and EMT of HT29 and SW480.

CONCLUSION

miR-410-3p-mediated ZCCHC10 suppression regulates NF-κB activation, thereby promoting EMT process, cell migration and invasion of CRC cells. This study provides a new insight into the specific mechanism by which miR-410-3p mediates CRC progression.

Detection of haplotype-dependent allele-specific DNA methylation in WGBS data

In heterozygous genomes, allele-specific measurements can reveal biologically significant differences in DNA methylation between homologous alleles associated with local changes in genetic sequence. Current approaches for detecting such events from whole-genome bisulfite sequencing (WGBS) data perform statistically independent marginal analysis at individual cytosine-phosphate-guanine (CpG) sites, thus ignoring correlations in the methylation state, or carry-out a joint statistical analysis of methylation patterns at four CpG sites producing unreliable statistical evidence. Here, we employ the one-dimensional Ising model of statistical physics and develop a method for detecting allele-specific methylation (ASM) events within segments of DNA containing clusters of linked single-nucleotide polymorphisms (SNPs), called haplotypes. Comparisons with existing approaches using simulated and real WGBS data show that our method provides an improved fit to data, especially when considering large haplotypes. Importantly, the method employs robust hypothesis testing for detecting statistically significant imbalances in mean methylation level and methylation entropy, as well as for identifying haplotypes for which the genetic variant carries significant information about the methylation state. As such, our ASM analysis approach can potentially lead to biological discoveries with important implications for the genetics of complex human diseases.

MicroRNA expression profiles in the esophagus of children with caustic stenosis: A pathway towards esophageal cancer?

BACKGROUND

Eighty percent of caustic ingestions occur in children and esophageal neoplasms may develop as a late complication of such injury. The identification of biomarkers is a promising strategy to improve early diagnosis of esophageal cancer or caustic lesions that are at an increased risk of progression.

STUDY DESIGN/AIMS

This study aimed at identifying global microRNA (miRNA) expression changes in esophageal mucosa from children with caustic stenosis. The study included 27 biopsy samples from 15 patients. Samples were divided into two groups, according to the time elapsed after injury (N = 15 in Group A, with less than five years of follow-up and N = 12 in Group B, with more than five years of follow-up). miRNA expression profiles were determined in each lesion, compared with normal esophageal tissues from control group. We used the TaqMan Human MicroRNA Arrays (Thermo Fisher) platform. Furthermore, bioinformatic algorithms were used to identify miRNA target genes and biological pathways including miRNAs and their target genes potentially associated with esophageal disease.

RESULTS

Thirteen miRNAs were significantly deregulated (9 over- and 4 underexpressed) in patients from Group A. In patients from Group B, two miRNAs were over- and two were underexpressed. Of note, miR-374 and miR-574 were deregulated in Group B patients and have been linked to esophageal tumorigenesis. We identified signal transduction and transcription factor networks with genes strongly related to development and progression of esophageal cancer.

CONCLUSION

miRNAs identified here contribute to a better understanding of pathways associated with malignant transformation from caustic stenosis to neoplastic lesions. This study may serve as a basis for validation of miRNAs, including miR-374 and miR-574, as potential biomarkers of early cancer detection.

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.

A Regulatory Noncoding RNA, nc886, Suppresses Esophageal Cancer by Inhibiting the AKT Pathway and Cell Cycle Progression

nc886 is a regulatory non-coding RNA (ncRNA) whose expression is frequently silenced in malignancies. In the case of esophageal squamous cell carcinoma (ESCC), nc886 silencing is associated with shorter survival of patients, suggesting nc886’s tumor suppressor role in ESCC. However, this observation has not been complemented by an in-detail study about nc886’s impact on gene expression and cellular phenotypes. Here we have shown that nc886 inhibits AKT, a key protein in a renowned pro-survival pathway in cancer. nc886-silenced cells (nc886⁻ cells) have activated AKT and altered expression of cell cycle genes. nc886⁻ cells tend to have lower expression of CDKN2A and CDKN2C, both of which are inhibitors for cyclin-dependent kinase (CDK), and higher expression of CDK4 than nc886-expressing cells. As a result, nc886⁻ cells are hyperactive in the progression of the G1 to S cell cycle phase, proliferate faster, and are more sensitive to palbociclib, which is a cancer therapeutic drug that targets CDK4/6. Experimentally by nc886 expression and knockdown, we have determined the AKT target genes and cell cycle genes that are controlled by nc886 (nc886-associated gene sets). These gene sets, in combination with pathologic staging and nc886 expression levels, are a vastly superior predictor for the survival of 108 ESCC patients. In summary, our study has elucidated in ESCC how nc886 inhibits cell proliferation to explain its tumor suppressor role and identified gene sets that are of future clinical utility, by predicting patient survival and responsiveness to a therapeutic drug.

vtRNA2-1/nc886 Produces a Small RNA That Contributes to Its Tumor Suppression Action through the microRNA Pathway in Prostate Cancer

vtRNA2-1 is a vault RNA initially classified as microRNA precursor hsa-mir-886 and recently proposed as "nc886", a new type of non-coding RNA involved in cancer progression acting as an oncogene and tumor suppressor gene in different tissues. We have shown that vtRNA2-1/nc886 is epigenetically repressed in neoplastic cells, increasing cell proliferation and invasion in prostate tissue. Here we investigate the ability of vtRNA2-1/nc886 to produce small-RNAs and their biological effect in prostate cells. The interrogation of public small-RNA transcriptomes of prostate and other tissues uncovered two small RNAs, snc886-3p and snc886-5p, derived from vtRNA2-1/nc886 (previously hsa-miR-886-3p and hsa-miR-886-5p). Re-analysis of PAR-CLIP and knockout of microRNA biogenesis enzymes data showed that these small RNAs are products of DICER, independent of DROSHA, and associate with Argonaute proteins, satisfying microRNA attributes. In addition, the overexpression of snc886-3p provokes the downregulation of mRNAs bearing sequences complementary to its "seed" in their 3'-UTRs. Microarray and in vitro functional assays in DU145, LNCaP and PC3 cell lines revealed that snc886-3p reduced cell cycle progression and increases apoptosis, like its precursor vtRNA2-1/nc886. Finally, we found a list of direct candidate targets genes of snc886-3p upregulated and associated with disease condition and progression in PRAD-TCGA data. Overall, our findings suggest that vtRNA2-1/nc886 and its processed product snc886-3p are synthesized in prostate cells, exerting a tumor suppressor action.

'High vault-age': non-coding RNA control of autophagy

RNA-binding proteins typically change the fate of RNA, such as stability, translation or processing. Conversely, we recently uncovered that the small non-coding vault RNA 1-1 (vtRNA1-1) directly binds to the autophagic receptor p62/SQSTM1 and changes the protein's function. We refer to this process as 'riboregulation'. Here, we discuss this newly uncovered vault RNA function against the background of three decades of vault RNA research. We highlight the vtRNA1-1-p62 interaction as an example of riboregulation of a key cellular process.

Uncovering Tumour Heterogeneity through PKR and nc886 Analysis in Metastatic Colon Cancer Patients Treated with 5-FU-Based Chemotherapy

Colorectal cancer treatment has advanced over the past decade. The drug 5-fluorouracil is still used with a wide percentage of patients who do not respond. Therefore, a challenge is the identification of predictive biomarkers. The protein kinase R (PKR also called EIF2AK2) and its regulator, the non-coding pre-mir-nc886, have multiple effects on cells in response to numerous types of stress, including chemotherapy. In this work, we performed an ambispective study with 197 metastatic colon cancer patients with unresectable metastases to determine the relative expression levels of both nc886 and PKR by qPCR, as well as the location of PKR by immunohistochemistry in tumour samples and healthy tissues (plasma and colon epithelium). As primary end point, the expression levels were related to the objective response to first-line chemotherapy following the response evaluation criteria in solid tumours (RECIST) and, as the second end point, with survival at 18 and 36 months. Hierarchical agglomerative clustering was performed to accommodate the heterogeneity and complexity of oncological patients’ data. High expression levels of nc886 were related to the response to treatment and allowed to identify clusters of patients. Although the PKR mRNA expression was not associated with chemotherapy response, the absence of PKR location in the nucleolus was correlated with first-line chemotherapy response. Moreover, a relationship between survival and the expression of both PKR and nc886 in healthy tissues was found. Therefore, this work evaluated the best way to analyse the potential biomarkers PKR and nc886 in order to establish clusters of patients depending on the cancer outcomes using algorithms for complex and heterogeneous data.

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.

ZCCHC10 suppresses lung cancer progression and cisplatin resistance by attenuating MDM2-mediated p53 ubiquitination and degradation

The activation of p53 tumor suppressor is essential for preventing abnormal cell proliferation and carcinogenesis. ZCCHC10 was previously identified as a potential p53-interacting partner in a yeast two-hybrid screen, but the interaction in cells and its subsequent influence on p53 activity and cancer development have not been investigated. In this paper, we demonstrate that ZCCHC10 expression levels are statistically lower in lung adenocarcinoma tissues than the corresponding adjacent noncancerous tissues, and decreased expression of ZCCHC10 mRNA predicts poorer survival of the patients. Ectopic expression of ZCCHC10 in lung cancer cells harboring wild-type p53 dramatically suppresses cell proliferation, colony formation, migration, invasion and cisplatin resistance in vitro, as well as tumor growth and metastasis in vivo. Conversely, knockdown of ZCCHC10 exerts opposite effects in the normal lung cell Beas-2b. However, ZCCHC10 has no influence on the biological behaviors of p53-null (H358) or p53-mutant (H1437) lung cancer cells. Mechanistically, ZCCHC10 binds and stabilizes p53 by disrupting the interaction between p53 and MDM2. The p53 inhibitor pifithrin-α attenuated the influences of ZCCHC10 overexpression on p53 pathway, cell cycle, apoptosis, and epithelial-mesenchymal transition, whereas the p53 activator Nutlin3 could reverse the effects of ZCCHC10 knockdown. Collectively, our results indicate that ZCCHC10 exerts its tumor-suppressive effects by stabilizing the p53 protein and can be used a potential prognostic marker and therapeutic target in lung adenocarcinoma.

Role of microRNAs, circRNAs and long noncoding RNAs in acute myeloid leukemia

Acute myeloid leukemia (AML) is a malignant tumor of the immature myeloid hematopoietic cells in the bone marrow (BM). It is a highly heterogeneous disease, with rising morbidity and mortality in older patients. Although researches over the past decades have improved our understanding of AML, its pathogenesis has not yet been fully elucidated. Long noncoding RNAs (lncRNAs), microRNAs (miRNAs), and circular RNAs (circRNAs) are three noncoding RNA (ncRNA) molecules that regulate DNA transcription and translation. With the development of RNA-Seq technology, more and more ncRNAs that are closely related to AML leukemogenesis have been discovered. Numerous studies have found that these ncRNAs play an important role in leukemia cell proliferation, differentiation, and apoptosis. Some may potentially be used as prognostic biomarkers. In this systematic review, we briefly described the characteristics and molecular functions of three groups of ncRNAs, including lncRNAs, miRNAs, and circRNAs, and discussed their relationships with AML in detail.

Mechanism mediated by a noncoding RNA, nc886, in the cytotoxicity of a DNA-reactive compound

DNA-reactive compounds target actively proliferating cells. Therefore, they are presupposed to kill cancer cells selectively, and many of them are used as chemotherapeutic agents. In this study, we have discovered a cell death pathway involving nc886 and PKR as another mechanism for the cytotoxicity. Our study provides an insight how a proapoptotic protein responds to a DNA-reactive compound via a regulatory noncoding RNA (ncRNA) as a molecular signal. Since the nc886/PKR pathway operates in most normal cells including nonproliferating ones, our finding may answer to a conundrum why a DNA-damaging compound harms quiescent cells and is of future clinical utility by considering nc886/PKR when designing a chemotherapeutic regimen with minimal side effects on normal cells.

DNA-reactive compounds are harnessed for cancer chemotherapy. Their genotoxic effects are considered to be the main mechanism for the cytotoxicity to date. Because this mechanism preferentially affects actively proliferating cells, it is postulated that the cytotoxicity is specific to cancer cells. Nonetheless, they do harm normal quiescent cells, suggesting that there are other cytotoxic mechanisms to be uncovered. By employing doxorubicin as a representative DNA-reactive compound, we have discovered a cytotoxic mechanism that involves a cellular noncoding RNA (ncRNA) nc886 and protein kinase R (PKR) that is a proapoptotic protein. nc886 is transcribed by RNA polymerase III (Pol III), binds to PKR, and prevents it from aberrant activation in most normal cells. We have shown here that doxorubicin evicts Pol III from DNA and, thereby, shuts down nc886 transcription. Consequently, the instantaneous depletion of nc886 provokes PKR and leads to apoptosis. In a short-pulse treatment of doxorubicin, these events are the main cause of cytotoxicity preceding the DNA damage response in a 3D culture system as well as the monolayer cultures. By identifying nc886 as a molecular signal for PKR to sense doxorubicin, we have provided an explanation for the conundrum why DNA-damaging drugs can be cytotoxic to quiescent cells that have the competent nc886/PKR pathway.

nc886, a non-coding RNA, inhibits UVB-induced MMP-9 and COX-2 expression via the PKR pathway in human keratinocytes

nc886, a long non-coding RNA (ncRNA) of 101 nucleotides in length, is known as a vault RNA or microRNA precursor. Despite the recent discovery that ncRNAs in the nucleus play a crucial role in regulating chromosomal transformation and transcription, only a few studies have focused on the function of ncRNAs in the cytoplasm, such as nc886. Several studies have investigated the function of nc886 as a suppressor of carcinogenesis and inflammation in different cancer cell types; however, its role in the skin has yet to be clearly elucidated. The two RNA binding sites for protein kinase RNA-activated (PKR) are located in the central region of the stable structure of nc886, which competes with other double-stranded RNA species. Successful binding results in decreased PKR activity. Among changes in skin cells induced by ultraviolet B (UVB) radiation, nc886 expression decreases, whereas PKR phosphorylation via mitogen-activated protein kinases (MAPKs) increases. Reduced nc886 expression leads to uncontrolled PKR activity and increases in the expression of inflammatory cytokines, matrix metalloproteinase-9 (MMP-9), type IV collagenase, and cyclooxygenase (COX-2), which ultimately accelerate inflammatory responses and skin aging. The present study investigated the regulatory mechanism associated with PKR activity and nc886-PKR binding in skin cell aging and inflammation. These results suggest a role for nc886 in controlling photoaging and inflammation in skin cells.

The Noncoding RNA nc886 Regulates PKR Signaling and Cytokine Production in Human Cells

Protein kinase RNA-activated (PKR) is a cytoplasmic receptor for dsRNA, and as such is involved in detection of viral infection. On binding dsRNA, PKR dimerizes, autophosphorylates, and then phosphorylates its substrate, eukaryotic translation initiation factor 2 subunit α (eIF2α), causing inhibition of mRNA translation and shutdown of viral protein production. However, active PKR has also been found to be involved in the NF-κB signaling pathway by inducing phosphorylation of IκBα. PKR is regulated by the noncoding RNA nc886, which has altered expression in cancer. We have found that expression of nc886 is highly upregulated during activation of human CD4⁺ T cells. As has been described in other cell types, nc886 bound to PKR in human T cell lysates, preventing PKR phosphorylation by polyinosinic:polycytidylic acid or HIV trans-activation response element RNA in lysates of T cell lines or primary human CD4⁺ T cells. Using clonal human T cell lines, we found that nc886 expression was strictly required for IFN-γ and IL-2 expression and secretion after T cell activation but did not affect proliferation or activation-induced cell death. In stimulated human PBMCs, nc886 expression strongly correlated with IFN-γ expression. Although nc886 inhibited PKR activation by dsRNA, it was required for PKR phosphorylation during T cell stimulation, with subsequent NF-κB signaling and CREB phosphorylation. nc886 also regulated PKR phosphorylation during human monocyte-derived macrophage activation. We have therefore identified nc886 as a noncoding RNA marker of T cell activation and regulator of PKR-dependent signaling.

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.

Prenatal arsenic exposure is associated with increased plasma IGFBP3 concentrations in 9-year-old children partly via changes in DNA methylation

Exposure to inorganic arsenic (As), a carcinogen and epigenetic toxicant, has been associated with lower circulating levels of insulin-like growth factor 1 (IGF1) and impaired growth in children of pre-school age. The aim of this study was to assess the potential impact of exposure to As on IGF1 and insulin-like growth factor-binding protein 3 (IGFBP3) as well as DNA methylation changes in 9-year-old children. To this end, we studied 9-year-old children from a longitudinal mother-child cohort in rural Bangladesh (n = 551). Prenatal and concurrent exposure to As was assessed via concentrations in maternal urine at gestational week 8 and in child urine at 9 years, measured by HPLC-HG-ICPMS. Plasma IGF1 and IGFBP3 concentrations were quantified with immunoassays. DNA methylation was measured in blood mononuclear cells at 9 years in a sub-sample (n = 113) using the Infinium HumanMethylation450K BeadChip. In multivariable-adjusted linear regression models, prenatal As (natural log-transformed), but not children's concurrent urinary As, was positively associated with IGFBP3 concentrations (β = 76, 95% CI 19, 133). As concentrations were not associated with IGF1. DNA methylation analysis revealed CpGs associated with both prenatal As and IGFBP3. Mediation analysis suggested that methylation of 12 CpG sites for all children was mediator of effect for the association between prenatal As and IGFBP3. We also found differentially methylated regions, generally hypermethylated, that were associated with both prenatal As and IGFBP3. In all, our study revealed that prenatal exposure to As was positively associated with IGFBP3 concentrations in children at 9 years, independent of IGF1, and this association may, at least in part, be epigenetically mediated.

Emerging concepts on the role of epigenetics in the relationships between nutrition and health

Understanding the physiological and metabolic underpinnings that confer individual differences in responses to diet and diet-related chronic disease is essential to advance the field of nutrition. This includes elucidating the differences in gene expression that are mediated through programming of the genome through epigenetic chromatin modifications. Epigenetic landscapes are influenced by age, genetics, toxins and other environmental factors, including dietary exposures and nutritional status. Epigenetic modifications influence transcription and genome stability are established during development with life-long consequences. They can be inherited from one generation to the next. The covalent modifications of chromatin, which include methylation and acetylation, on DNA nucleotide bases, histone proteins and RNA are derived from intermediates of one-carbon metabolism and central metabolism. They influence key physiological processes throughout life, and together with inherited DNA primary sequence, contribute to responsiveness to environmental stresses, diet and risk for age-related chronic disease. Revealing diet-epigenetic relationships has the potential to transform nutrition science by increasing our fundamental understanding of: (i) the role of nutrients in biological systems, (ii) the resilience of living organisms in responding to environmental perturbations, and (iii) the development of dietary patterns that programme physiology for life-long health. Epigenetics may also enable the classification of individuals with chronic disease for specific dietary management and/or for efficacious diet-pharmaceutical combination therapies. These new emerging concepts at the interface of nutrition and epigenetics were discussed, and future research needs identified by leading experts at the 26th Marabou Symposium entitled 'Nutrition, Epigenetics, Genetics: Impact on Health and Disease'. For a compilation of the general discussion at the marabou symposium, click here http://www.marabousymposium.org/.

Deoxyribonucleic Acid Methylation in Systemic Lupus Erythematosus: Implications for Future Clinical Practice

Differential deoxyribonucleic acid (DNA) methylation has emerged as a critical feature of systemic lupus erythematosus (SLE). Genome-wide DNA methylation studies have revealed methylation patterns characteristic of SLE—in particular, robust hypomethylation of interferon-regulated genes is a prominent finding in all cells of the immune system studied to date. These patterns reliably distinguish individuals with SLE from healthy controls and from individuals with other autoimmune diseases. For example, hypomethylation within IFI44L is both highly sensitive and highly specific for SLE, superior to currently available biomarkers. Furthermore, methylation status of other genetic loci has been associated with clinically relevant features of SLE including disease severity and organ-specific manifestations. Finally, DNA methylation studies have provided important insights into the pathophysiology of SLE. Most recently, there is a growing body of evidence that the transcription factor enhancer of zeste homolog 2 (EZH2) plays an important role in triggering SLE disease activity via epigenetic mechanisms, and that EZH2 blockade may be a future treatment option in SLE. In this short review, we discuss the DNA methylation patterns associated with SLE, their relationship to clinically significant features of SLE, and their implications in the development of novel diagnostic and therapeutic approaches to this complex disease.

nc886 is induced by TGF-β and suppresses the microRNA pathway in ovarian cancer

Transforming growth factor-β (TGF-β) signaling and microRNAs (miRNAs) are important gene regulatory components in cancer. Usually in advanced malignant stages, TGF-β signaling is elevated but global miRNA expression is suppressed. Such a gene expression signature is well illustrated in a fibrosis (or mesenchymal) subtype of ovarian cancer (OC) that is of poor prognosis. However, the interplay between the two pathways in the OC subtype has not yet been elucidated. nc886 is a recently identified non-coding RNA implicated in several malignancies. The high expression of nc886 is associated with poor prognosis in 285 OC patients. Herein, we find that in OC nc886 expression is induced by TGF-β and that nc886 binds to Dicer to inhibit miRNA maturation. By preventing the miRNA pathway, nc886 emulates TGF-β in gene expression patterns and potentiates cell adhesion, migration, invasion, and drug resistance. Here we report nc886 to be a molecular link between the TGF-β and miRNA pathways.

Occupational exposure to pesticides is associated with differential DNA methylation

Objectives

Occupational pesticide exposure is associated with a wide range of diseases, including lung diseases, but it is largely unknown how pesticides influence airway disease pathogenesis. A potential mechanism might be through epigenetic mechanisms, like DNA methylation. Therefore, we assessed associations between occupational exposure to pesticides and genome-wide DNA methylation sites.

Methods

1561 subjects of LifeLines were included with either no (n=1392), low (n=108) or high (n=61) exposure to any type of pesticides (estimated based on current or last held job). Blood DNA methylation levels were measured using Illumina 450K arrays. Associations between pesticide exposure and 420 938 methylation sites (CpGs) were assessed using robust linear regression adjusted for appropriate confounders. In addition, we performed genome-wide stratified and interaction analyses by gender, smoking and airway obstruction status, and assessed associations between gene expression and methylation for genome-wide significant CpGs (n=2802).

Results

In total for all analyses, high pesticide exposure was genome-wide significantly (false discovery rate P<0.05) associated with differential DNA methylation of 31 CpGs annotated to 29 genes. Twenty of these CpGs were found in subjects with airway obstruction. Several of the identified genes, for example, RYR1, ALLC, PTPRN2, LRRC3B, PAX2 and VTRNA2-1, are genes previously linked to either pesticide exposure or lung-related diseases. Seven out of 31 CpGs were associated with gene expression levels.

Conclusions

We show for the first time that occupational exposure to pesticides is genome-wide associated with differential DNA methylation. Further research should reveal whether this differential methylation plays a role in the airway disease pathogenesis induced by pesticides.