High‐resolution mapping of function and protein binding in an RNA nuclear enrichment sequence

The functions of long RNAs, including mRNAs and long noncoding RNAs (lncRNAs), critically depend on their subcellular localization. The identity of the sequences that dictate subcellular localization and their high‐resolution anatomy remain largely unknown. We used a suite of massively parallel RNA assays and libraries containing thousands of sequence variants to pinpoint the functional features within the SIRLOIN element, which dictates nuclear enrichment through hnRNPK recruitment. In addition, we profiled the endogenous SIRLOIN RNA‐nucleoprotein complex and identified the nuclear RNA‐binding proteins SLTM and SNRNP70 as novel SIRLOIN binders. Taken together, using massively parallel assays, we identified the features that dictate binding of hnRNPK, SLTM, and SNRNP70 to SIRLOIN and found that these factors are jointly required for SIRLOIN activity. Our study thus provides a roadmap for high‐throughput dissection of functional sequence elements in long RNAs.     

m6A‐mediated alternative splicing coupled with nonsense‐mediated mRNA decay regulates SAM synthetase homeostasis

Alternative splicing of pre‐mRNAs can regulate gene expression levels by coupling with nonsense‐mediated mRNA decay (NMD). In order to elucidate a repertoire of mRNAs regulated by alternative splicing coupled with NMD (AS‐NMD) in an organism, we performed long‐read RNA sequencing of poly(A)⁺ RNAs from an NMD‐deficient mutant strain of Caenorhabditis elegans, and obtained full‐length sequences for mRNA isoforms from 259 high‐confidence AS‐NMD genes. Among them are the S‐adenosyl‐L‐methionine (SAM) synthetase (sams) genes sams‐3 and sams‐4. SAM synthetase activity autoregulates sams gene expression through AS‐NMD in a negative feedback loop. We furthermore find that METT‐10, the orthologue of human U6 snRNA methyltransferase METTL16, is required for the splicing regulation in␣vivo, and specifically methylates the invariant AG dinucleotide at the distal 3′ splice site (3′SS) in␣vitro. Direct RNA sequencing coupled with machine learning confirms m⁶A modification of endogenous sams mRNAs. Overall, these results indicate that homeostasis of SAM synthetase in C. elegans is maintained by alternative splicing regulation through m⁶A modification at the 3′SS of the sams genes.     

Circular RNA VRK1 facilitates pre‐eclampsia progression via sponging miR‐221‐3P to regulate PTEN/Akt

Pre‐eclampsia (PE) is a worldwide pregnancy‐related disorder. It is mainly characterized by defect migration and invasion of trophoblast cells. Recently, circular RNAs (circRNAs) have been believed to play a vital role in PE. The expression patterns and the biological functions of circRNAs in PE remain elusive. Here, we performed a circRNA microarray to identify putative PE‐related circRNAs. Bioinformatics analyses were used to screen the circRNAs which have potential relationships with pre‐eclampsia, and we identified a novel circRNA (circVRK1) that was up‐regulated in PE placenta tissues. By using HTR‐8/SVneo cells, circVRK1 knockdown significantly enhanced cell migration and invasion abilities, as well as epithelial‐mesenchymal transition (EMT). Mechanistically, we found that circVRK1 and PTEN could function as the ceRNAs to miR‐221‐3p. Overexpression of miR‐221‐3p promoted cell migration, invasion and EMT via regulating PTEN. The cotransfection of miR‐221‐3p inhibitor or PTEN reversed the effect from circVRK1 knockdown. Moreover, the circVRK1/miR‐221‐3p/PTEN axis greatly regulated Akt phosphorylation. In general, circVRK1 suppresses trophoblast cell migration, invasion and EMT, by acting as a ceRNA to miR‐221‐3p to regulate PTEN, and further inhibit PI3K/Akt activation. The purpose of this paper is to open wide insights to investigate the onset of PE and provide new potential therapeutic targets in PE.     

Long non‐coding RNA TINCR suppresses metastatic melanoma dissemination by preventing ATF4 translation

Transition from proliferative‐to‐invasive phenotypes promotes metastasis and therapy resistance in melanoma. Reversion of the invasive phenotype, however, is challenged by the poor understanding of mechanisms underlying its maintenance. Here, we report that the lncRNA TINCR is down‐regulated in metastatic melanoma and its silencing increases the expression levels of invasive markers, in vitro migration, in vivo tumor growth, and resistance to BRAF and MEK inhibitors. The critical mediator is ATF4, a central player of the integrated stress response (ISR), which is activated in TINCR‐depleted cells in the absence of starvation and eIF2α phosphorylation. TINCR depletion increases global protein synthesis and induces translational reprogramming, leading to increased translation of mRNAs encoding ATF4 and other ISR proteins. Strikingly, re‐expression of TINCR in metastatic melanoma suppresses the invasive phenotype, reduces numbers of tumor‐initiating cells and metastasis formation, and increases drug sensitivity. Mechanistically, TINCR interacts with mRNAs associated with the invasive phenotype, including ATF4, preventing their binding to ribosomes. Thus, TINCR is a suppressor of the melanoma invasive phenotype, which functions in nutrient‐rich conditions by repressing translation of selected ISR RNAs.     

Regulation and roles of RNA modifications in aging‐related diseases

With the aging of the global population, accumulating interest is focused on manipulating the fundamental aging‐related signaling pathways to delay the physiological aging process and eventually slow or prevent the appearance or severity of multiple aging‐related diseases. Recently, emerging evidence has shown that RNA modifications, which were historically considered infrastructural features of cellular RNAs, are dynamically regulated across most of the RNA species in cells and thereby critically involved in major biological processes, including cellular senescence and aging. In this review, we summarize the current knowledge about RNA modifications and provide a catalog of RNA modifications on different RNA species, including mRNAs, miRNAs, lncRNA, tRNAs, and rRNAs. Most importantly, we focus on the regulation and roles of these RNA modifications in aging‐related diseases, including neurodegenerative diseases, cardiovascular diseases, cataracts, osteoporosis, and fertility decline. This would be an important step toward a better understanding of fundamental aging mechanisms and thereby facilitating the development of novel diagnostics and therapeutics for aging‐related diseases.     

Hsa_circ_0003945 promotes progression of hepatocellular carcinoma by mediating miR‐34c‐5p/LGR4/β‐catenin axis activity

Accumulating evidence suggests that circular RNAs (circRNAs) play essential roles in regulating cancer progression, but many circRNAs in hepatocellular carcinoma (HCC) remain unknown. Dysregulated circRNAs in HCC were identified through bioinformatics analysis of Gene Expression Omnibus data sets. Quantitative real‐time PCR (qRT‐PCR), Sanger sequencing, RNase R digestion and actinomycin D treatment were conducted to confirm the characterization of circRNAs. CCK‐8, wound‐healing and Transwell assays were performed to assess the functional roles of Hsa_circ_0003945 (Circ_0003945) in HCC cell lines. Subcellular fractionation and fluorescence in situ hybridization (FISH) were performed to locate Circ_0003945 in HCC cells. Dual‐luciferase reporter assay was executed to verify the binding of Circ_0003945 to microRNAs (miRNAs) or the miRNAs to their target genes. In this study, we found that Circ_0003945 was upregulated in HCC tissue, and higher Circ_0003945 expression was positively correlated with tumour size and tumour stage. Furthermore, high plasma levels of circulating Circ_0003945 were confirmed in HCC patients compared with those in non‐HCC groups. The functional experiments revealed that overexpression or knockdown of Circ_0003945 promoted or attenuated tumour growth and migration, respectively. Mechanistically, Circ_0003945 might exert as a miR‐34c‐5p sponge to upregulate the expression of leucine‐rich repeat‐containing G protein‐coupled receptor 4 (LGR4), activating the β‐catenin pathway, and finally facilitating HCC progression. Additionally, a β‐catenin activator could reverse the effect of Circ_0003945 knockdown. In conclusion, Circ_0003945 exerts a tumour‐promoting role in HCC cells by regulating the miR‐34c‐5p/LGR4/β‐catenin axis, which may be a potential target for HCC therapy.     

Prediction of a competing endogenous RNA co‐expression network as a prognostic marker in glioblastoma

Due to its high proliferation capacity and rapid intracranial spread, glioblastoma (GBM) has become one of the least curable malignant cancers. Recently, the competing endogenous RNAs (ceRNAs) hypothesis has become a focus in the researches of molecular biological mechanisms of cancer occurrence and progression. However, there is a lack of correlation studies on GBM, as well as a lack of comprehensive analyses of GBM molecular mechanisms based on high‐throughput sequencing and large‐scale sample sizes. We obtained RNA‐seq data from The Cancer Genome Atlas (TCGA) and Genotype‐Tissue Expression (GTEx) databases. Further, differentially expressed mRNAs were identified from normal brain tissue and GBM tissue. The similarities between the mRNA modules with clinical traits were subjected to weighted correlation network analysis (WGCNA). With the mRNAs from clinical‐related modules, a survival model was constructed by univariate and multivariate Cox proportional hazard regression analyses. Thereafter, we carried out Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Finally, we predicted interactions between lncRNAs, miRNAs and mRNAs by TargetScan, miRDB, miRTarBase and starBase. We identified 2 lncRNAs (NORAD, XIST), 5 miRNAs (hsa‐miR‐3613, hsa‐miR‐371, hsa‐miR‐373, hsa‐miR‐32, hsa‐miR‐92) and 2 mRNAs (LYZ, PIK3AP1) for the construction of a ceRNA network, which might act as a prognostic biomarker of GBM. Combined with previous studies and our enrichment analysis results, we hypothesized that this ceRNA network affects immune activities and tumour microenvironment variations. Our research provides novel aspects to study GBM development and treatment.     

Characterization and Complexity of Wheat Developing Endosperm mRNAs

Free and membrane-bound (MB) polysomes and the corresponding polyadenylated RNAs (polyA⁺ RNAs) have been isolated from developing wheat endosperm (Triticum aestivum L.) Free and MB poly(A)⁺ RNAs, analyzed on isokinetic sucrose gradient with [³H]polyuridylic acid [poly(U)] hybridization detection, appear to be 11S to 12S in size with a 7% poly(A) tail for MB RNAs. cDNAs synthesized using both of these mRNA populations in presence of a potent RNase inhibitor (RNasin), have been used for hybridization kinetics experiments. The mean square fitting analysis of the hybridization kinetics between MB cDNA and its template reveals the presence of two abundance classes representing roughly ⅔ and ⅓ of the MB poly(A)⁺ RNAs and containing the information for approximately 75 superabundant species (21,000 copies per cell) and 750 intermediate species (530 copies per cell), respectively. The mRNA population extracted from free polysomes is divided into three abundance classes. The first one is composed of superabundant sequences which would correspond to the MB superabundant mRNAs. The free mRNAs consist of about 11,000 diverse sequences, most of them being rare sequences. Heterologous hybridizations of MB cDNAs to free mRNAs have shown that some mRNAs are common to both populations. This could be explained either by a partial contamination or by free polysomes en route to their membrane destination. Contrary to the low number of diverse mRNAs corresponding to the legume seed storage proteins, the wheat endosperm superabundant mRNAs consist of about 75 different sequences which would encode most of the seed storage proteins, especially gliadins.     

Physiologically relevant miRNAs in mammalian oocytes are rare and highly abundant

miRNAs, ~22nt small RNAs associated with Argonaute (AGO) proteins, are important negative regulators of gene expression in mammalian cells. However, mammalian maternal miRNAs show negligible repressive activity and the miRNA pathway is dispensable for oocytes and maternal‐to‐zygotic transition. The stoichiometric hypothesis proposed that this is caused by dilution of maternal miRNAs during oocyte growth. As the dilution affects miRNAs but not mRNAs, it creates unfavorable miRNA:mRNA stoichiometry for efficient repression of cognate mRNAs. Here, we report that porcine ssc‐miR‐205 and bovine bta‐miR‐10b are exceptional miRNAs, which resist the diluting effect of oocyte growth and can efficiently suppress gene expression. Additional analysis of ssc‐miR‐205 shows that it has higher stability, reduces expression of endogenous targets, and contributes to the porcine oocyte‐to‐embryo transition. Consistent with the stoichiometric hypothesis, our results show that the endogenous miRNA pathway in mammalian oocytes is intact and that maternal miRNAs can efficiently suppress gene expression when a favorable miRNA:mRNA stoichiometry is established.     

Cap‐independent translation of GPLD1 enhances markers of brain health in long‐lived mutant and drug‐treated mice

Glycosylphosphatidylinositol‐specific phospholipase D1 (GPLD1) hydrolyzes inositol phosphate linkages in proteins anchored to the cell membrane. Mice overexpressing GPLD1 show enhanced neurogenesis and cognition. Snell dwarf (DW) and growth hormone receptor knockout (GKO) mice show delays in age‐dependent cognitive decline. We hypothesized that augmented GPLD1 might contribute to retained cognitive function in these mice. We report that DW and GKO show higher GPLD1 levels in the liver and plasma. These mice also have elevated levels of hippocampal brain‐derived neurotrophic factor (BDNF) and of doublecortin (DCX), suggesting a mechanism for maintenance of cognitive function at older ages. GPLD1 was not increased in the hippocampus of DW or GKO mice, suggesting that plasma GPLD1 increases elevated these brain proteins. Alteration of the liver and plasma GPLD1 was unaltered in mice with liver‐specific GHR deletion, suggesting that the GH effect was not intrinsic to the liver. GPLD1 was also induced by caloric restriction and by each of four drugs that extend lifespan. The proteome of DW and GKO mice is molded by selective translation of mRNAs, involving cap‐independent translation (CIT) of mRNAs marked by N⁶ methyladenosine. Because GPLD1 protein increases were independent of the mRNA level, we tested the idea that GPLD1 might be regulated by CIT. 4EGI‐1, which enhances CIT, increased GPLD1 protein without changes in GPLD1 mRNA in cultured fibroblasts and mice. Furthermore, transgenic overexpression of YTHDF1, which promotes CIT by reading m6A signals, also led to increased GPLD1 protein, showing that elevation of GPLD1 reflects selective mRNA translation.     

m6A modification of HSATIII lncRNAs regulates temperature‐dependent splicing

Nuclear stress bodies (nSBs) are nuclear membraneless organelles formed around stress‐inducible HSATIII architectural long noncoding RNAs (lncRNAs). nSBs repress splicing of hundreds of introns during thermal stress recovery, which are partly regulated by CLK1 kinase phosphorylation of temperature‐dependent Ser/Arg‐rich splicing factors (SRSFs). Here, we report a distinct mechanism for this splicing repression through protein sequestration by nSBs. Comprehensive identification of RNA‐binding proteins revealed HSATIII association with proteins related to N⁶‐methyladenosine (m⁶A) RNA modification. 11% of the first adenosine in the repetitive HSATIII sequence were m⁶A‐modified. nSBs sequester the m⁶A writer complex to methylate HSATIII, leading to subsequent sequestration of the nuclear m⁶A reader, YTHDC1. Sequestration of these factors from the nucleoplasm represses m⁶A modification of pre‐mRNAs, leading to repression of m⁶A‐dependent splicing during stress recovery phase. Thus, nSBs serve as a common platform for regulation of temperature‐dependent splicing through dual mechanisms employing two distinct ribonucleoprotein modules with partially m⁶A‐modified architectural lncRNAs.     

Circular RNA circ‐RCCD promotes cardiomyocyte differentiation in mouse embryo development via recruiting YY1 to the promoter of MyD88

Congenital heart disease (CHD) is the most common birth defect, affecting approximately 1% of live births. Genetic and environmental factors are leading factors to CHD, but the mechanism of CHD pathogenesis remains unclear. Circular RNAs (circRNAs) are kinds of endogenous non‐coding RNAs (ncRNAs) involved in a variety of physiological and pathological processes, especially in heart diseases. In this study, three significant differently expressed circRNA between maternal embryonic day (E) E13 and E17 was found by microarray assay. Among them, the content of circ‐RCCD increases with the development of heart and was enriched in primary cardiomyocytes of different species, which arouses our attention. Functional experiments revealed that inhibition of circ‐RCCD dramatically suppressed the formation of beating cell clusters, the fluorescence intensity of cardiac differentiation marker MF20, and the expression of the myocardial‐specific markers CTnT, Mef2c, and GATA4. Next, we found that circ‐RCCD was involved in cardiomyocyte differentiation through negative regulation of MyD88 expression. Further experiments proved that circ‐RCCD inhibited MyD88 levels by recruiting YY1 to the promoter of MyD88; circ‐RCCD inhibited nuclear translocation of YY1. These results reported that circ‐RCCD promoted cardiomyocyte differentiation by recruiting YY1 to the promoter of MyD88. And, this study provided a potential role and molecular mechanism of circ‐RCCD as a target for the treatment of CHD.