Protein-coding mRNAs Research Articles

Transcriptomic profiles of the ovaries from piglets neonatally exposed to 4-tert-octylphenol

The environmental pollutants with hormonal activities may influence steroid-mediated processes in neonatal ovaries and increase the incidence of reproductive disorders. The aim of the current study was to examine effects of 4-tert-octylphenol (OP), a non-ionic surfactant widely used in a variety of industrial applications which has been reported to mimic the 17β-estradiol activity, on the expression of protein-coding (mRNAs) and long non-coding (lncRNAs) transcripts in neonatal ovaries of the pig. By employing RNA-Seq we aimed to gain insights into regulatory networks underlying the OP effects on the follicular development in pigs. Piglets were injected (sc) daily with OP (100 mg/kg bw) or corn oil (controls) between postnatal Days 1 and 10 (n = 3/group). Ovaries were excised from the 11-day-old piglets and total cellular RNA was isolated and sequenced. Two hundred three differentially expressed genes (DEGs; P-adjusted < 0.05 and log2 fold change ≥1.0) and 23 differentially expressed lncRNAs (DELs; P-adjusted < 0.05 and log2 fold change ≥ 1.0) were identified in OP-treated piglet ovaries. The DEGs were assigned to Gene Ontology terms, covering biological processes, molecular functions and cellular components, which linked the DEGs to functions associated with movement of cell or subcellular component, regulation of plasma membrane bounded cell projection assembly as well as hydrolase and endopeptidase activity. In addition, STRING analysis demonstrated the strongest interactions between genes related to negative regulation of endopeptidase activity. Some correlations between DEGs and DELs were also found, revealing that the OP action on the ovary may be partially executed via the changes in the lncRNA expression. These results suggest that neonatal exposure of pigs to OP induces changes in the ovarian transcriptomic profile associated with genes encoding serine protease inhibitors and involved in steroid synthesis as well as genes linked to intracellular and membrane transport. We suggest that the changes in the mRNA and lncRNA expression in the ovaries of OP-treated piglets may disturb ovarian cellular function, including steroidogenesis, proliferation and apoptosis.

Long Non-Coding RNA and mRNA Profiling in Early-Stage Bovine Embryos Treated with Glutathione.

We measured differential expression profiles of genes and long non-coding RNA (lncRNA) using RNA sequencing in bovine embryos with or without glutathione (GSH) treatment. Bovine embryos fertilized in vitro were treated with GSH to blastocyst. Embryos at the 8-16-cell and morula stages were collected, with embryos without GSH treatment as the control. RNA was isolated, amplified, and sequenced. Differentially expressed genes (DEGs) and lncRNAs (DElncRNAs) were identified and bioinformatic analyses carried out. Transcript levels were confirmed using quantitative RT-PCR. A total of 4100 DEGs were identified, of which 3952 were in GSH-treated morulae and 884 in untreated morulae. More gene ontology (GO) terms were associated with GSH treatment than with control conditions. KEGG analysis showed that glutathione metabolism, citrate cycle, and metabolic pathways involving glycine, serine, and threonine were observed only in GSH-treated embryos. Among 4273 DElncRNAs identified, 59 were potentially important in GSH-treated embryo development, including 14 involved in glutathione metabolism. The 59 DElncRNAs co-expressed with protein-coding mRNAs involved similar GO terms and pathways as the DEGs. This appears to be the first comprehensive profiling of DEGs and DElncRNAs in bovine embryos fertilized in vitro with or without GSH, and the first systematic screen of potential lncRNAs in bovine embryos.

Learning-induced mRNA alterations in olfactory bulb mitral cells in neonatal rats

In the olfactory bulb, a cAMP/PKA/CREB-dependent form of learning occurs in the first week of life that provides a unique mammalian model for defining the epigenetic role of this evolutionarily ancient plasticity cascade. Odor preference learning in the week-old rat pup is rapidly induced by a 10-min pairing of odor and stroking. Memory is demonstrable at 24 h, but not 48 h, posttraining. Using this paradigm, pups that showed peppermint preference 30 min posttraining were sacrificed 20 min later for laser microdissection of odor-encoding mitral cells. Controls were given odor only. Microarray analysis revealed that 13 nonprotein-coding mRNAs linked to mRNA translation and splicing and 11 protein-coding mRNAs linked to transcription differed with odor preference training. MicroRNA23b, a translation inhibitor of multiple plasticity-related mRNAs, was down-regulated. Protein-coding transcription was up-regulated for Sec23b, Clic2, Rpp14, Dcbld1, Magee2, Mstn, Fam229b, RGD1566265, and Mgst2. Gng12 and Srcg1 mRNAs were down-regulated. Increases in Sec23b, Clic2, and Dcbld1 proteins were confirmed in mitral cells in situ at the same time point following training. The protein-coding changes are consistent with extracellular matrix remodeling and ryanodine receptor involvement in odor preference learning. A role for CREB and AP1 as triggers of memory-related mRNA regulation is supported. The small number of gene changes identified in the mitral cell input/output link for 24 h memory will facilitate investigation of the nature, and reversibility, of changes supporting temporally restricted long-term memory.

Splicing conservation signals in plant long noncoding RNAs.

Long noncoding RNAs (lncRNAs) have recently emerged as prominent regulators of gene expression in eukaryotes. LncRNAs often drive the modification and maintenance of gene activation or gene silencing states via chromatin conformation rearrangements. In plants, lncRNAs have been shown to participate in gene regulation, and are essential to processes such as vernalization and photomorphogenesis. Despite their prominent functions, only over a dozen lncRNAs have been experimentally and functionally characterized. Similar to its animal counterparts, the rates of sequence divergence are much higher in plant lncRNAs than in protein coding mRNAs, making it difficult to identify lncRNA conservation using traditional sequence comparison methods. Beyond this, little is known about the evolutionary patterns of lncRNAs in plants. Here, we characterized the splicing conservation of lncRNAs in Brassicaceae. We generated a whole-genome alignment of 16 Brassica species and used it to identify synthenic lncRNA orthologs. Using a scoring system trained on transcriptomes from A. thaliana and B. oleracea, we identified splice sites across the whole alignment and measured their conservation. Our analysis revealed that 17.9% (112/627) of all intergenic lncRNAs display splicing conservation in at least one exon, an estimate that is substantially higher than previous estimates of lncRNA conservation in this group. Our findings agree with similar studies in vertebrates, demonstrating that splicing conservation can be evidence of stabilizing selection. We provide conclusive evidence for the existence of evolutionary deeply conserved lncRNAs in plants and describe a generally applicable computational workflow to identify functional lncRNAs in plants.

Transcriptome analysis reveals the link between lncRNA-mRNA co-expression network and tumor immune microenvironment and overall survival in head and neck squamous cell carcinoma

BackgroundAs the sixth most common cancer worldwide, head and neck squamous cell carcinoma (HNSCC) develops visceral metastases during the advanced stage of the disease and exhibits a low five-year survival rate. The importance of tumor microenvironment (TME) in tumor initiation and metastasis is widely recognized. In addition, accumulating evidence indicates that long non-coding RNA (lncRNA) is involved in crosstalk between TME and tumor cells. However, the lncRNA-associated regulators modulating the HNSCC microenvironment and progression remain largely unknown.MethodsThe publicly available transcriptome data and matched clinical HNSCC data were collected from The Cancer Genome Atlas (TCGA). Immune scores (ISs) and stromal scores (SSs) of HNSCC TME were calculated using ESTIMATE algorithm. Weighted gene co-expression network analysis (WGCNA) was conducted to determine the co-expressed lncRNAs and protein-coding mRNAs.ResultsResults showed that the high IS HNSCC male patient subgroup exhibited improved survival. Additionally, we identified 169 lncRNAs and 825 protein-coding mRNAs that were differentially expressed in high IS HNSCC samples, with the up-regulated mRNAs displaying enrichment in immune-related biological processes. Notably, we identified a high co-expression lncRNA-mRNA module (i.e., purple module) that showed strong correlation with ISs. This module contained 79 lncRNAs and 442 mRNAs, including 26 lncRNAs and 215 mRNAs showing association between expression and male HNSCC survival. Consistently, 207 of the 215 mRNAs were up-regulated in high IS HNSCC group and were enriched in immune-related signaling pathways. Based on bioinformatics analyses and previous functional assays, certain lncRNAs (e.g., AL365361.1 and PCED1B-AS1) in the purple module likely contributed to the modification of tumor immune microenvironment (TIME) in the high IS HNSCC patients, achieved by regulating transcription of abundant immune-related genes (e.g., CCR7 and TLR8).ConclusionsIn summary, we ascertained a HNSCC male patient subgroup that displayed high ISs and good survival probability. We identified hundreds of genes with specific expression patterns in this HNSCC subgroup as well as a highly co-expressed lncRNA-mRNA module with great potential for the modulation of TIME of HNSCC. Our study provides evidence of a link between the lncRNA-associated gene network, TIME, and HNSCC progression, and highlights potential therapeutic targets for this disease.

Differential expression of long noncoding RNAs from dental pulp stem cells in the microenvironment of the angiogenesis

IntroductionAngiogenesis is important in pulp-dentin formation. Among the regulatory factors, long noncoding RNA (LncRNA) is a class of functional RNA molecules that are not translated into protein and involved in regulating multiple physiological processes. The different expression of LncRNA and its target gene in dental pulp stem cells (DPSCs) were explored and may provide a theoretical basis for future regulation of dental pulp angiogenesis. MethodsIn this study, we cultured DPSCs from healthy dental pulp tissues and divided them into two groups: the normal DPSCs and the DPSCs cultured in vascular induction medium. In total, 40,173 LncRNA probes and 20,730 protein coding mRNAs were detected through microarray, which were then verified by the quantitative reverse transcription-polymerase chain reaction (qRT-PCR) method. ResultsThe result of differential expressions measured in LncRNA through microarray showed that 376 LncRNAs increased significantly and 426 were downregulated among the two groups of cells. Moreover, the mRNA microarray in normal cultured DPSCs showed that 629 LncRNAs were significantly upregulated, while 529 of them were downregulated compared with the DPSCs that were cultured in vascular induction medium. Gene ontology (GO) analysis inferred the molecular function of mRNAs. Pathway analysis showed that 52 signaling pathways were involved in the differentiation process of DPSCs. qRT-PCR analysis, conducted for validation, showed results consistent with the microarray analysis. ConclusionsWe found that a number of different regulators are involved in inducing vascular differentiation of DPSCs, which provides a foundation for subsequent experiments.

Abstract P1-05-06: Functional characterization of lncRNA152 in ER+ and triple-negative breast cancers

Abstract Genomic analyses have revealed that many RNA species other than protein-coding mRNAs are transcribed from mammalian genomes and may have profound effects on cellular physiology and pathology. One such class of RNAs are the long noncoding RNAs (lncRNAs). In our previous work using Global Run-On and sequencing (GRO-seq), we identified and annotated ~1,900 lncRNAs in MCF-7 human breast cancer cells, which are differentially expressed across the distinct molecular subtypes of breast cancer. Further meta-analysis identified lncRNA152 as a predictor of increased metastasis-free survival. Thus, lncRNA152 levels may serve as a biomarker to track disease progression in breast cancer patients. However, the precise mechanisms by which lncRNA152 regulates these clinical outcomes remain largely unknown. Interestingly, lncRNA152 is highly expressed in the luminal subtypes of breast cancer, but downregulated in triple-negative breast cancers (TNBC), suggesting that lncRNA152 expression may be associated with more aggressive cancer features. Through multiple complementary experimental approaches, we found that (1) FoxA1 regulates lncRNA152 transcription, (2) knockdown of lncRNA152 inhibits the growth, but promotes cell migration and invasion, of luminal breast cancer cell subtypes, and (3) ectopic overexpression of lncRNA152 inhibits TNBC cell growth, migration, and invasion. Importantly, our xenograft studies provide evidence that lncRNA152 inhibits estrogen-dependent tumor growth of luminal and TNBC cells. In addition, transcriptome analysis of TNBC cell-derived xenografts indicates that lncRNA152 expression downregulates many genes that are involved in biological processes such as cell growth, migration, invasion, and angiogenesis. Collectively, these results suggest that lncRNA152 is a novel tumor suppressor in luminal and basal subtype of breast cancer. This work is supported by grants from the Cancer Prevention and Research Institute of Texas (CPRIT RP160318/RP190235) to W.L.K. Citation Format: Daeseok Kim. Functional characterization of lncRNA152 in ER+ and triple-negative breast cancers [abstract]. In: Proceedings of the 2019 San Antonio Breast Cancer Symposium; 2019 Dec 10-14; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2020;80(4 Suppl):Abstract nr P1-05-06.

The 3' Untranslated Region Protects the Heart from Angiotensin II-Induced Cardiac Dysfunction via AGGF1 Expression.

The messenger RNA (mRNA) 3' untranslated regions (3' UTRs), as cis-regulated elements bound by microRNAs (miRNAs), affect their gene translation. However, the role of the trans-regulation of 3' UTRs during heart dysfunction remains elusive. Compared with administration of angiogenic factor with G-patch and forkhead-associate domains 1 (Aggf1), ectopic expression of Aggf1 with its 3' UTR significantly suppressed cardiac dysfunction in angiotensin II-infused mice, with upregulated expression of both Aggf1 and myeloid cell leukemia 1 (Mcl1). Along their 3' UTRs, Mcl1 and Aggf1 mRNAs share binding sites for the same miRNAs, including miR-105, miR-101, and miR-93. We demonstrated that the protein-coding Mcl1 and Aggf1 mRNAs communicate and co-regulate each other's expression through competition for these three miRNAs that target both transcripts via their 3' UTRs. Our results indicate that Aggf1 3' UTR, as a trans-regulatory element, accelerates the cardioprotective role of Aggf1 in response to hypertensive conditions by elevating Mcl1 expression. Our work broadens the scope of gene therapy targets and provides a new insight into gene therapy strategies involving 3' UTRs.

A CRISPR/Cas13-based approach demonstrates biological relevance of vlinc class of long non-coding RNAs in anticancer drug response

Long non-coding (lnc) RNAs represent a fascinating class of transcripts that remains highly controversial mainly due to ambiguity surrounding overall biological relevance of these RNAs. Multitude of reverse genetics studies showing functionality of lncRNAs are unfortunately based on assays that are either plagued by non-specific effects and/or cannot unambiguously assign observed phenotypes to the transcript per se. Here, we show application of the novel CRISPR/Cas13 RNA knockdown system that has superior specificity compared to other transcript-targeting knockdown methods like RNAi. We applied this method to a novel widespread subclass of nuclear lncRNAs — very long intergenic non-coding (vlinc) RNAs — in a high-throughput phenotypic assay based on survival challenge in response to anticancer drug treatments. We used multiple layers of controls including mismatch control for each targeting gRNA to ensure uncovering true phenotype-transcript relationships. We found evidence supporting importance for cellular survival for up to 60% of the tested protein-coding mRNAs and, importantly, 64% of vlincRNAs. Overall, this study demonstrates utility of CRISPR/Cas13 as a highly sensitive and specific tool for reverse genetics study of both protein-coding genes and lncRNAs. Furthermore, importantly, this approach provides evidence supporting biological significance of the latter transcripts in anticancer drug response.

A MAFG-lncRNA axis links systemic nutrient abundance to hepatic glucose metabolism

Obesity and type 2 diabetes mellitus are global emergencies and long noncoding RNAs (lncRNAs) are regulatory transcripts with elusive functions in metabolism. Here we show that a high fraction of lncRNAs, but not protein-coding mRNAs, are repressed during diet-induced obesity (DIO) and refeeding, whilst nutrient deprivation induced lncRNAs in mouse liver. Similarly, lncRNAs are lost in diabetic humans. LncRNA promoter analyses, global cistrome and gain-of-function analyses confirm that increased MAFG signaling during DIO curbs lncRNA expression. Silencing Mafg in mouse hepatocytes and obese mice elicits a fasting-like gene expression profile, improves glucose metabolism, de-represses lncRNAs and impairs mammalian target of rapamycin (mTOR) activation. We find that obesity-repressed LincIRS2 is controlled by MAFG and observe that genetic and RNAi-mediated LincIRS2 loss causes elevated blood glucose, insulin resistance and aberrant glucose output in lean mice. Taken together, we identify a MAFG-lncRNA axis controlling hepatic glucose metabolism in health and metabolic disease.

M6A RNA modification modulates PI3K/Akt/mTOR signal pathway in Gastrointestinal Cancer.

Rationale: Methylation at the N6 position of adenosine (m6A) is the most prevalent RNA modification within protein-coding mRNAs in mammals, and it is a reversible modification with various important biological functions. The formation and function of m6A are regulated by methyltransferases (writers), demethylases (erasers), and special binding proteins (readers) as key factors. However, the underlying modification mechanisms of m6A in gastrointestinal (GI) cancer remain unclear. Here, we performed comprehensive molecular profiling of the nine known m6A writer, eraser, and reader proteins in GI cancer.Methods: Data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) were used. Gene alteration and pathway analysis were done in cBioportal. The protein network of m6A regulators and its related pathway members was analyzed in STRING online platform. Phylogenetic tree was constructed in MEGA7. m6A modification sites were predicted by SRAMP. m6A related SNPs were analyzed by m6ASNP. The modulation of m6A on its related pathway members was validated by m6A-seq, real-time PCR and phosphor-MAPK array.Results: We found that m6A regulators were mostly upregulated in GI cancer and their differential expression significantly influenced the overall survival of patients with GI cancer. The phosphatidylinositol-3-kinase (PI3K)/Akt and mammalian target of rapamycin (mTOR) signaling pathways were found to be potentially affected by m6A modification in most human cancers, including GI cancer, which was further verified by m6A-Seq and phospho-MAPK array.Conclusions: Our findings suggest that m6A RNA modification has a fundamental role in the regulation of PI3K/Akt and mTOR signaling pathway function in cancer.

Subcellular mRNA Localization Regulates Ribosome Biogenesis&amp;nbsp;In Migrating Cells

Translation of Ribosomal Protein-coding mRNAs (RP-mRNAs) constitutes a key step in ribosome biogenesis, but the mechanisms which modulate RP-mRNA translation in coordination with other cellular processes are poorly defined. Here we show that the subcellular localization of RP-mRNAs acts as a key regulator of their translation during cell migration. As cells migrate into their surroundings, RP-mRNAs localize to the actin-rich protrusive fronts of the cells. This localization is mediated by La-related protein 6 (LARP6), an RNA binding protein that is enriched in protrusions. Protrusions act as hotspots of translation for RP-mRNAs, resulting in enhancement of RP synthesis, ribosome biogenesis, and overall protein synthetic capacity of migratory cells, required for supporting efficient migration and proliferation. In human breast carcinomas, Epithelial to Mesenchymal Transition (EMT) upregulates LARP6 expression to enhance protein synthesis and support invasive growth. Our findings reveal LARP6 mediated mRNA localization as a key regulator of ribosome biogenesis during cell migration, and demonstrate a role for this process in cancer progression downstream of EMT.

Electrophoretic Mobility Shift Assay of in vitro Phosphorylated RNA Polymerase II Carboxyl-terminal Domain Substrates.

Eukaryotic RNA polymerase II transcribes all protein-coding mRNAs and is highly regulated. A key mechanism directing RNA polymerase II and facilitating the co-transcriptional processing of mRNAs is the phosphorylation of its highly repetitive carboxyl-terminal domain (CTD) of its largest subunit, RPB1, at specific residues. A variety of techniques exist to identify and quantify the degree of CTD phosphorylation, including phosphorylation-specific antibodies and mass spectrometry. Electrophoretic mobility shift assays (EMSAs) have been utilized since the discovery of CTD phosphorylation and continue to represent a simple, direct, and widely applicable approach for qualitatively monitoring CTD phosphorylation. We present a standardized method for EMSA analysis of recombinant GST-CTD substrates phosphorylated by a variety of CTD kinases. Strategies to analyze samples under both denatured/reduced and semi-native conditions are provided. This method represents a simple, direct, and reproducible means to monitor CTD phosphorylation in recombinant substrates utilizing equipment common to molecular biology labs and readily applicable to downstream analyses including immunoblotting and mass spectrometry.

Electroacupuncture Alleviates Experimental Chronic Inflammatory Pain by Inhibiting Calcium Voltage-Gated Channel-Mediated Inflammation.

Background Both experimental and clinical studies have shown that electroacupuncture (EA) administration ameliorates chronic inflammatory pain (CIP). However, the multifaceted mechanism underlying the effects of EA on CIP is poorly understood. In this study, the mRNA transcriptome was used to study various therapeutic targets of EA. Methods Using RNA-sequencing, protein-coding mRNA expression profiles of the L4-L5 dorsal root ganglion (DRG) were examined in the control (CN), complete Freund's adjuvant- (CFA-) induced CIP, and EA-treated CIP groups. A series of bioinformatics analyses was performed; “EA-reversed upregulated genes with CIP” (up-DEGs) and “EA-reversed downregulated genes with CIP” (down-DEGs) were identified. Thereafter, based on up-DEGs and down-DEGs, biological functions and signaling pathways were enriched using gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) pathway analyses. Results In total, 189 DEGs were identified, including 134 up- and 55 down-DEGs, which were enriched in arachidonic acid metabolism (rno00590), glutamatergic synapse (rno04724), serotonergic synapse (rno04726), FoxO signaling pathway (rno04068), insulin signaling pathway (rno04910), amyotrophic lateral sclerosis (rno05014), cholinergic synapse (rno04725), ECM-receptor interaction (rno04512), and choline metabolism in cancer (rno05231). Conclusion We identified a few GOs, pathways, and genes that could play key roles in the amelioration of CIP by EA. Hence, this study may provide a theoretical basis for CIP amelioration by EA.

Long Noncoding RNA (lncRNA)-Mediated Competing Endogenous RNA Networks Provide Novel Potential Biomarkers and Therapeutic Targets for Colorectal Cancer.

Colorectal cancer (CRC) is the third most common cancer and has a high metastasis and reoccurrence rate. Long noncoding RNAs (lncRNAs) play an important role in CRC growth and metastasis. Recent studies revealed that lncRNAs participate in CRC progression by coordinating with microRNAs (miRNAs) and protein-coding mRNAs. LncRNAs function as competitive endogenous RNAs (ceRNAs) by competitively occupying the shared binding sequences of miRNAs, thus sequestering the miRNAs and changing the expression of their downstream target genes. Such ceRNA networks formed by lncRNA/miRNA/mRNA interactions have been found in a broad spectrum of biological processes in CRC, including liver metastasis, epithelial to mesenchymal transition (EMT), inflammation formation, and chemo-/radioresistance. In this review, we summarize typical paradigms of lncRNA-associated ceRNA networks, which are involved in the underlying molecular mechanisms of CRC initiation and progression. We comprehensively discuss the competitive crosstalk among RNA transcripts and the novel targets for CRC prognosis and therapy.

The RNA-binding protein ILF3 binds to transposable element sequences in SINEUP lncRNAs.

Transposable elements (TEs) compose about half of the mammalian genome and, as embedded sequences, up to 40% of long noncoding RNA (lncRNA) transcripts. Embedded TEs may represent functional domains within lncRNAs, providing a structured RNA platform for protein interaction. Here we show the interactome profile of the mouse inverted short interspersed nuclear element (SINE) of subfamily B2 (invSINEB2) alone and embedded in antisense (AS) ubiquitin C-terminal hydrolase L1 (Uchl1), an lncRNA that is AS to Uchl1 gene. AS Uchl1 is the representative member of a functional class of AS lncRNAs, named SINEUPs, in which the invSINEB2 acts as effector domain (ED)–enhancing translation of sense protein-coding mRNAs. By using RNA-interacting domainome technology, we identify the IL enhancer-binding factor 3 (ILF3) as a protein partner of AS Uchl1 RNA. We determine that this interaction is mediated by the RNA-binding motif 2 of ILF3 and the invSINEB2. Furthermore, we show that ILF3 is able to bind a free right Arthrobacter luteus (Alu) monomer sequence, the embedded TE acting as ED in human SINEUPs. Bioinformatic analysis of Encyclopedia of DNA Elements–enhanced cross-linking immunoprecipitation data reveals that ILF3 binds transcribed human SINE sequences at transcriptome-wide levels. We then demonstrate that the embedded TEs modulate AS Uchl1 RNA nuclear localization to an extent moderately influenced by ILF3. This work unveils the existence of a specific interaction between embedded TEs and an RNA-binding protein, strengthening the model of TEs as functional modules in lncRNAs.—Fasolo, F., Patrucco, L., Volpe, M., Bon, C., Peano, C., Mignone, F., Carninci, P., Persichetti, F., Santoro, C., Zucchelli, S., Sblattero, D., Sanges, R., Cotella, D., Gustincich, S. The RNA-binding protein ILF3 binds to transposable element sequences in SINEUP lncRNAs.

The Construction and Analysis of ceRNA Network and Patterns of Immune Infiltration in Mesothelioma With Bone Metastasis.

Background: Mesothelioma is a rare and aggressive tumor. Bone metastasis often occurs in the later stages of this disease along with poor quality of life. Thus, it is important to explore the tumorigenesis and bone metastasis mechanism of invasive mesothelioma. For this purpose, we established two nomograms based on tumor-infiltrating immune cells and ceRNA networks to describe the molecular immunity and the clinical prediction of mesothelioma patients with bone metastasis.Method: The expression profiles of mRNAs, lncRNAs, and miRNAs of 87 primary mesotheliomas were obtained from the TCGA database; there were four patients with bone metastasis and 83 patients without. We constructed a ceRNAs network based on the differentially expressed RNAs between mesothelioma and bone metastasis. CIBERSORT was used to distinguish 22 immune cell types from the tumor transcriptomes. Kaplan–Meier survival analysis and the Cox proportional hazards model were used to evaluate the prognostic value of each factor. Prognosis-associated immune cells and ceRNAs were applied to establish prediction nomograms. The receiver operating characteristic curves (ROC) and calibration curves were utilized to assess the discrimination and accuracy of the nomogram.Results: Differential analysis revealed that 20 lncRNAs, 15 miRNAs, and 230 mRNAs were significantly different in mesothelioma samples vs. bone metastasis samples. We constructed the ceRNA network to include 10 protein-coding mRNAs, 8 lncRNAs, and 10 miRNAs. Nine of 28 ceRNAs were found to be significant in the Kaplan–Meier analysis. Out of the 22 cell types, the fraction of dendritic cells resting (P = 0.018) was significantly different between the bone metastasis group and the non-bone metastasis group. The ROC and the calibration curves, based on ceRNA networks and tumor-infiltrating immune cells, respectively, suggested acceptable accuracy (AUC of 3-year survival: 0.827, AUC of 5-year survival: 0.840; AUC of 3-year survival: 0.730; AUC of 5-year survival: 0.753). Notably, based on the co-expression patterns between ceRNAs and Immune cells, we found that the hsa-miR-582-5p, CASP9, dendritic cells resting, ANIX2, T cells CD8, and T cells CD4 memory resting might be associated with the mesothelioma bone metastasis.Conclusion: Based on ceRNA networks and patterns of immune infiltration, our study provided a valid bioinformatics basis in order to explore the molecular mechanism and predict the possibility of mesothelioma bone metastasis.

Transcriptome signature of cellular senescence.

Cellular senescence, an integral component of aging and cancer, arises in response to diverse triggers, including telomere attrition, macromolecular damage and signaling from activated oncogenes. At present, senescent cells are identified by the combined presence of multiple traits, such as senescence-associated protein expression and secretion, DNA damage and β-galactosidase activity; unfortunately, these traits are neither exclusively nor universally present in senescent cells. To identify robust shared markers of senescence, we have performed RNA-sequencing analysis across eight diverse models of senescence triggered in human diploid fibroblasts (WI-38, IMR-90) and endothelial cells (HUVEC, HAEC) by replicative exhaustion, exposure to ionizing radiation or doxorubicin, and expression of the oncogene HRASG12V. The intersection of the altered transcriptomes revealed 50 RNAs consistently elevated and 18 RNAs consistently reduced across all senescence models, including many protein-coding mRNAs and some non-coding RNAs. We propose that these shared transcriptome profiles will enable the identification of senescent cells in vivo, the investigation of their roles in aging and malignancy and the development of strategies to target senescent cells therapeutically.

Weighted Gene Co-Expression Analyses Point to Long Non-Coding RNA Hub Genes at Different Schistosoma mansoni Life-Cycle Stages

Long non-coding RNAs (lncRNAs) (>200 nt) are expressed at levels lower than those of the protein-coding mRNAs, and in all eukaryotic model species where they have been characterized, they are transcribed from thousands of different genomic loci. In humans, some four dozen lncRNAs have been studied in detail, and they have been shown to play important roles in transcriptional regulation, acting in conjunction with transcription factors and epigenetic marks to modulate the tissue-type specific programs of transcriptional gene activation and repression. In Schistosoma mansoni, around 10,000 lncRNAs have been identified in previous works. However, the limited number of RNA-sequencing (RNA-seq) libraries that had been previously assessed, together with the use of old and incomplete versions of the S. mansoni genome and protein-coding transcriptome annotations, have hampered the identification of all lncRNAs expressed in the parasite. Here we have used 633 publicly available S. mansoni RNA-seq libraries from whole worms at different stages (n = 121), from isolated tissues (n = 24), from cell-populations (n = 81), and from single-cells (n = 407). We have assembled a set of 16,583 lncRNA transcripts originated from 10,024 genes, of which 11,022 are novel S. mansoni lncRNA transcripts, whereas the remaining 5,561 transcripts comprise 120 lncRNAs that are identical to and 5,441 lncRNAs that have gene overlap with S. mansoni lncRNAs already reported in previous works. Most importantly, our more stringent assembly and filtering pipeline has identified and removed a set of 4,293 lncRNA transcripts from previous publications that were in fact derived from partially processed mRNAs with intron retention. We have used weighted gene co-expression network analyses and identified 15 different gene co-expression modules. Each parasite life-cycle stage has at least one highly correlated gene co-expression module, and each module is comprised of hundreds to thousands lncRNAs and mRNAs having correlated co-expression patterns at different stages. Inspection of the top most highly connected genes within the modules’ networks has shown that different lncRNAs are hub genes at different life-cycle stages, being among the most promising candidate lncRNAs to be further explored for functional characterization.

Spatial transcriptome profiling by MERFISH reveals subcellular RNA compartmentalization and cell cycle-dependent gene expression

The expression profiles and spatial distributions of RNAs regulate many cellular functions. Image-based transcriptomic approaches provide powerful means to measure both expression and spatial information of RNAs in individual cells within their native environment. Among these approaches, multiplexed error-robust fluorescence in situ hybridization (MERFISH) has achieved spatially resolved RNA quantification at transcriptome scale by massively multiplexing single-molecule FISH measurements. Here, we increased the gene throughput of MERFISH and demonstrated simultaneous measurements of RNA transcripts from ∼10,000 genes in individual cells with ∼80% detection efficiency and ∼4% misidentification rate. We combined MERFISH with cellular structure imaging to determine subcellular compartmentalization of RNAs. We validated this approach by showing enrichment of secretome transcripts at the endoplasmic reticulum, and further revealed enrichment of long noncoding RNAs, RNAs with retained introns, and a subgroup of protein-coding mRNAs in the cell nucleus. Leveraging spatially resolved RNA profiling, we developed an approach to determine RNA velocity in situ using the balance of nuclear versus cytoplasmic RNA counts. We applied this approach to infer pseudotime ordering of cells and identified cells at different cell-cycle states, revealing ∼1,600 genes with putative cell cycle-dependent expression and a gradual transcription profile change as cells progress through cell-cycle stages. Our analysis further revealed cell cycle-dependent and cell cycle-independent spatial heterogeneity of transcriptionally distinct cells. We envision that the ability to perform spatially resolved, genome-wide RNA profiling with high detection efficiency and accuracy by MERFISH could help address a wide array of questions ranging from the regulation of gene expression in cells to the development of cell fate and organization in tissues.