Non-coding Transcriptome Research Articles

Comprehensive analysis of the coding and non-coding RNA transcriptome expression profiles of hippocampus tissue in tx-J animal model of Wilson's disease

Wilson's disease (WD) is an autosomal recessive disorder with a genetic basis. The predominant non-motor symptom of WD is cognitive dysfunction, although the specific genetic regulatory mechanism remains unclear. Tx-J mice, with an 82% sequence homology of the ATP7B gene to the human gene, are considered the most suitable model for WD. This study employs deep sequencing to investigate the differences in RNA transcript profiles, both coding and non-coding, as well as the functional characteristics of the regulatory network involved in WD cognitive impairment. The cognitive function of tx-J mice was evaluated using the Water Maze Test (WMT). Long non-coding RNA (lncRNA), circular RNA (circRNA), and messenger RNA (mRNA) profiles were analyzed in the hippocampal tissue of tx-J mice to identify differentially expressed RNAs (DE-RNAs). Subsequently, the DE-RNAs were used to construct protein–protein interaction (PPI) networks, as well as DE-circRNAs and lncRNAs-associated competing endogenous RNA (ceRNA) expression networks, and coding-noncoding co-expression (CNC) networks. To elucidate their biological functions and pathways, the PPI and ceRNA networks were subjected to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. A total of 361 differentially expressed mRNAs (DE-mRNAs), comprising 193 up-regulated and 168 down-regulated mRNAs, 2627 differentially expressed long non-coding RNAs (DE-lncRNAs), consisting of 1270 up-regulated and 1357 down-regulated lncRNAs, and 99 differentially expressed circular RNAs (DE-circRNAs), consisting of 68 up-regulated and 31 down-regulated circRNAs, were observed in the tx-J mice group when compared to the control mice group. Gene Ontology (GO) and pathway analyses revealed that DE-mRNAs were enriched in cellular processes, calcium signaling pathways, and mRNA surveillance pathways. In contrast, the DE-circRNAs-associated competing endogenous RNA (ceRNA) network was enriched for covalent chromatin modification, histone modification, and axon guidance, whereas the DE-lncRNAs-associated ceRNA network was enriched for dendritic spine, regulation of cell morphogenesis involved in differentiation, and mRNA surveillance pathway. The study presented the expression profiles of lncRNA, circRNA, and mRNA in the hippocampal tissue of tx-J mice. Furthermore, the study constructed PPI, ceRNA, and CNC expression networks. The findings are significant in comprehending the function of regulatory genes in WD associated with cognitive impairment. These results also offer valuable information for the diagnosis and treatment of WD.

Coding and noncoding transcriptomes of NODULIN HOMEOBOX (NDX)-deficient Arabidopsis inflorescence

Arabidopsis NODULIN HOMEOBOX (NDX) is a plant-specific transcriptional regulator whose role in small RNA biogenesis and heterochromatin homeostasis has recently been described. Here we extend our previous transcriptomic analysis to the flowering stage of development. We performed mRNA-seq and small RNA-seq measurements on inflorescence samples of wild-type and ndx1-4 mutant (WiscDsLox344A04) Arabidopsis plants. We identified specific groups of differentially expressed genes and noncoding heterochromatic siRNA (hetsiRNA) loci/regions whose transcriptional activity was significantly changed in the absence of NDX. In addition, data obtained from inflorescence were compared with seedling transcriptomics data, which revealed development-specific changes in gene expression profiles. Overall, we provide a comprehensive data source on the coding and noncoding transcriptomes of NDX-deficient Arabidopsis flowers to serve as a basis for further research on NDX function.

Novel Interactomes of Cardiac Noncoding RNA Transcriptome under Altered Thyroid Hormonal States

Background: Noncoding RNAs are emerging as key players in cardiovascular diseases. Both Hypothyroidism (HypoTH; low thyroid hormones [THs]) and hyperthyroidism (HyperTH; high THs) can be deleterious to the heart. We recently showed the involvement of long noncoding RNAs (lncRNAs) in HypoTH under caloric restriction. However, a comprehensive view of the cardiac transcriptome and its interactions under altered TH states is lacking. We hypothesized that cardiac noncoding interaction networks are significantly altered in the cardiac-TH regulatory axis. Methods: In this randomized, blinded study (2.5 mon), adult C57 mixed-sex mice (n=8-10/group; 5 females, 3-5 males) were divided into Control, HypoTH (low iodine diet with oral methimazole), HyperTH (1 microgram [ug]/gram Thyroxine, T4; intraperitoneal) and HypoTH+T3 (similar to HypoTH with oral 12 ug/kilogram/day Triiodothyronine, T3). Morphometrics, ultrasound, immunosorbent assay, Quantitative real-time PCR (qPCR) and whole transcriptome sequencing were performed. Statistical analyses were done by Tukey’s analysis of variance (with post hoc), Benjamini & Hochberg and others using Prism, R, etc. RNA quality was validated by Nanodrop, gel electrophoresis and Bioanalyzer. Results: In mixed-sex mice, HypoTH showed significant cardiac atrophy (whole heart, LV and right ventricle), low serum total T3 and T4 levels (p<0.001), low heart rate, (p<0.01) diminished systolic LV posterior (p<0.05) and septal wall thicknesses (p<0.01) and fractional shortening (-34%; p<0.01) compared to controls. T3 therapy showed feedback inhibition of T4 (p<0.001) and restored all other parameters including LV systolic diameter (p<0.05). Among others, HyperTH had cardiac hypertrophy and increased T3 and T4 (p<0.05) compared to controls (p<0.001). Similar significant results were obtained in males and females. qPCR of shorter-term treatment showed significantly altered lncRNAs in inflammatory and immune pathways. Whole transcriptome sequencing showed numerous differentially expressed, significantly altered and enriched ( GO and KEGG) targets across the groups among all RNA species studied (p-adjusted; padj<0.05). These include novel and known lncRNAs, microRNAs (miRNAs) and messenger (m) RNAs that change in opposite directions in low vs high TH groups. These also include 11 novel circular RNAs (circRNAs; CIRI2) decreased with HyperTH (padj<0.05). We also uncovered novel co-expression and interactions of miRNA-mRNA, miRNA-lncRNA, mRNA-lncRNA, circRNA-miRNA and miRNA-lncRNA-mRNA (including TH receptor alpha; miRanda, Cytoscape) and also splicing and base editing events (padj<0.05). Conclusions: This is the first cross-disciplinary study reporting a comprehensive cardiac transcriptome-wide interactome network of circRNA, lncRNA, miRNA and mRNA in a TH-responsive manner. This offers significant translational potential to develop diadiagnostic and therapeutic targets for cardiac-endocrine co-morbid conditions. NYITCOM at A-State, ABI This is the full abstract presented at the American Physiology Summit 2023 meeting and is only available in HTML format. There are no additional versions or additional content available for this abstract. Physiology was not involved in the peer review process.

Dysregulation of murine long noncoding single-cell transcriptome in nonalcoholic steatohepatitis and liver fibrosis.

LncRNAs comprise a heterogeneous class of RNA-encoding genes typified by low expression, nuclear enrichment, high tissue-specificity, and functional diversity, but the vast majority remain uncharacterized. Here, we assembled the mouse liver noncoding transcriptome from >2000 bulk RNA-seq samples and discovered 48,261 liver-expressed lncRNAs, a majority novel. Using these lncRNAs as a single-cell transcriptomic reference set, we elucidated lncRNA dysregulation in mouse models of high fat diet-induced nonalcoholic steatohepatitis and carbon tetrachloride-induced liver fibrosis. Trajectory inference analysis revealed lncRNA zonation patterns across the liver lobule in each major liver cell population. Perturbations in lncRNA expression and zonation were common in several disease-associated liver cell types, including nonalcoholic steatohepatitis-associated macrophages, a hallmark of fatty liver disease progression, and collagen-producing myofibroblasts, a central feature of liver fibrosis. Single-cell-based gene regulatory network analysis using bigSCale2 linked individual lncRNAs to specific biological pathways, and network-essential regulatory lncRNAs with disease-associated functions were identified by their high network centrality metrics. For a subset of these lncRNAs, promoter sequences of the network-defined lncRNA target genes were significantly enriched for lncRNA triplex formation, providing independent mechanistic support for the lncRNA-target gene linkages predicted by the gene regulatory networks. These findings elucidate liver lncRNA cell-type specificities, spatial zonation patterns, associated regulatory networks, and temporal patterns of dysregulation during hepatic disease progression. A subset of the liver disease-associated regulatory lncRNAs identified have human orthologs and are promising candidates for biomarkers and therapeutic targets.

Evidence for Existence of Multiple Functional Human Small RNAs Derived from Transcripts of Protein-Coding Genes.

The human genome encodes a multitude of different noncoding transcripts that have been traditionally separated on the basis of their lengths into long (>200 nt) or small (<200 nt) noncoding RNAs. The functions, mechanisms of action, and biological relevance of the vast majority of both long and short noncoding transcripts remain unknown. However, according to the functional understanding of the known classes of long and small noncoding RNAs (sncRNAs) that have been shown to play crucial roles in multiple biological processes, it is generally assumed that many unannotated long and small transcripts participate in important cellular functions as well. Nevertheless, direct evidence of functionality is lacking for most noncoding transcripts, especially for sncRNAs that are often dismissed as stable degradation products of longer RNAs. Here, we developed a high-throughput assay to test the functionality of sncRNAs by overexpressing them in human cells. Surprisingly, we found that a significant fraction (>40%) of unannotated sncRNAs appear to have biological relevance. Furthermore, contrary to the expectation, the potentially functional transcripts are not highly abundant and can be derived from protein-coding mRNAs. These results strongly suggest that the small noncoding transcriptome can harbor multiple functional transcripts that warrant future studies.

Ataxin-1 controls the expression of specific noncoding RNAs in B cells upon autoimmune demyelination.

B cells play a key mechanistic role in the pathogenesis of multiple sclerosis (MS), a chronic neurological disease of the central nervous system with an autoimmune etiology. B cells contribute to disease initiation and progression by acting as professional antigen-presenting cells as well as via secreting autoantibodies and proinflammatory cytokines. We have recently shown that the polyglutamine protein ataxin-1, which was first linked to the movement disorder spinocerebellar ataxia type 1, also acts as a master regulator of B-cell functions in the context of central nervous system autoimmunity. In fact, ataxin-1-deficient mice display an aggravated manifestation of the MS disease model experimental autoimmune encephalomyelitis along with aberrant B-cell functions. Consistent with this scenario, transcriptomic analysis of Atxn1-null Bcells highlighted distinct genetic signatures involved in cell activation, proliferation and antigen presentation. To further characterize the role of ataxin-1, we profiled the noncoding transcriptome controlled by ataxin-1 in the B-cell compartment upon an encephalitogenic challenge. We show that two specific classes of noncoding RNAs, namely, processed pseudogenes and intergenic long noncoding RNAs, are differentially regulated along disease. Furthermore, pathway and protein network analyses on their putative protein-coding gene targets found a significant enrichment in ontologies related to cell mitosis, together with molecular processes relevant to MS such as chitin metabolism. Altogether, these findings shed light on the possible contribution of noncoding RNAs to B-cell biology and MS pathogenesis, and further establish the immunomodulatory role of ataxin-1 in autoimmune demyelination.

Integrative analysis of the expression profiles of whole coding and non-coding RNA transcriptomes and construction of the competing endogenous RNA networks for chronic obstructive pulmonary disease.

The pathogenesis of Chronic Obstructive Pulmonary Disease (COPD) is implicated in airway inflammation, oxidative stress, protease/anti-protease and emphysema. Abnormally expressed non-coding RNAs (ncRNAs) play a vital role in regulation of COPD occurrence and progression. The regulatory mechanisms of the circRNA/lncRNA-miRNA-mRNA (competing endogenous RNA, ceRNA) networks might facilitate our cognition of RNA interactions in COPD. This study aimed to identified novel RNA transcripts and constructed the potential ceRNA networks of COPD patients. Total transcriptome sequencing of the tissues from patients with COPD (COPD) (n = 7) and non-COPD control subjects (Normal) (n = 6) was performed, and the expression profiles of differentially expressed genes (DEGs), including mRNAs, lncRNAs, circRNAs, and miRNAs, were analyzed. The ceRNA network was established based on the miRcode and miRanda databases. Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene Ontology (GO), Gene Set Enrichment Analysis (GSEA), and Gene set variation analysis (GSVA) were implemented for functional enrichment analysis of DEGs. Finally, CIBERSORTx was extracted to analyze the relevance between hub genes and various immune cells.The Starbase and JASPAR databases were used to construct hub-RNA binding proteins (RBPs) and lncRNA-transcription factor (TF) interaction networks. A total of 1,796 mRNAs, 2,207 lncRNAs, and 11 miRNAs showed differentially expression between the lung tissue samples from the normal and COPD groups. Based on these DEGs, lncRNA/circRNA-miRNA-mRNA ceRNA networks were constructed respectively. In addition, ten hub genes were identified. Among them, RPS11, RPL32, RPL5, and RPL27A were associated with the proliferation, differentiation, and apoptosis of the lung tissue. The biological function revealed that TNF-α via NF-kB and IL6/JAK/STAT3 signaling pathways were involved in COPD. Our research constructed the lncRNA/circRNA-miRNA-mRNA ceRNA networks, filtrated ten hub genes may regulate the TNF-α/NF-κB, IL6/JAK/STAT3 signally pathways, which indirectly elucidated the post-transcriptional regulation mechanism of COPD and lay the foundation for excavating the novel targets of diagnosis and treatment in COPD.

Multiplexed target enrichment of coding and non-coding transcriptomes enables studying Candida spp. infections from human derived samples.

The study of transcriptomic interactions between host and pathogens in in vivo conditions is challenged by the low relative amounts of the pathogen RNA. Yeast opportunistic pathogens of the genus Candida can cause life-threatening systemic infections in immunocompromised patients, and are of growing medical concern. Four phylogenetically diverse species account for over 90% of Candida infections, and their specific interactions with various human tissues are still poorly understood. To enable in vivo transcriptomic analysis in these species, we designed and validated pan-Candida target capture probes to enrich protein-coding and non-coding transcriptomes. The probe-based enrichment approach outperformed enrichment based on differential lysis of host cells, and showed similar enrichment performance as an existing capture design, yet achieving better fidelity of expression levels, enabling species multiplexing and capturing of lncRNAs. In addition, we show that our probe-based enrichment strategy allows robust genotype-based identification of the infecting strain present in the sample.

Nano3P-seq: transcriptome-wide analysis of gene expression and tail dynamics using end-capture nanopore cDNA sequencing

RNA polyadenylation plays a central role in RNA maturation, fate, and stability. In response to developmental cues, polyA tail lengths can vary, affecting the translation efficiency and stability of mRNAs. Here we develop Nanopore 3′ end-capture sequencing (Nano3P-seq), a method that relies on nanopore cDNA sequencing to simultaneously quantify RNA abundance, tail composition, and tail length dynamics at per-read resolution. By employing a template-switching-based sequencing protocol, Nano3P-seq can sequence RNA molecule from its 3′ end, regardless of its polyadenylation status, without the need for PCR amplification or ligation of RNA adapters. We demonstrate that Nano3P-seq provides quantitative estimates of RNA abundance and tail lengths, and captures a wide diversity of RNA biotypes. We find that, in addition to mRNA and long non-coding RNA, polyA tails can be identified in 16S mitochondrial ribosomal RNA in both mouse and zebrafish models. Moreover, we show that mRNA tail lengths are dynamically regulated during vertebrate embryogenesis at an isoform-specific level, correlating with mRNA decay. Finally, we demonstrate the ability of Nano3P-seq in capturing non-A bases within polyA tails of various lengths, and reveal their distribution during vertebrate embryogenesis. Overall, Nano3P-seq is a simple and robust method for accurately estimating transcript levels, tail lengths, and tail composition heterogeneity in individual reads, with minimal library preparation biases, both in the coding and non-coding transcriptome.

Single-cell analysis of the long noncoding RNA transcriptome in cardiac fibroblasts identifies FIXER, a novel therapeutic target for fibrosis in the damaged heart

Single-cell analysis of the long noncoding RNA transcriptome in cardiac fibroblasts identifies FIXER, a novel therapeutic target for fibrosis in the damaged heart

TMIC-76. HUMANIN RELEASE FROM TUMOR ASSOCIATED MYELOID CELLS PROMOTES GLIOMA CHEMORESISTANCE

Abstract Transcriptomic screens of brain tumor (glioblastoma; GBM) parenchymal cells indicated tumor supporting traits of the GBM microenvironment, but largely excluded mitochondrially enriched gene-sets. Here, we show that a mitochondrial transcript of GBM parenchymal cells contributes to therapy resistance. We inspected the non-coding transcriptome of human GBM associated myeloid cells (GAM) and observed an upregulation of the mitochondrial ribosomal subunit MT-RNR2, which contains an open reading frame for the signaling peptide humanin. Immunohistology disclosed that humanin was preponderant in GAM. In vitro assays with a range of human stem-like GBM cells (GSCs) revealed that nanomolar concentrations of humanin can drive tumor cell expansion and chemoresistance to temozolomide. A series of genetic, pharmacological and Western blotting experiments showed that extracellular humanin increased GSC viability via stimulation of the GP130/IL6ST receptor and MAPK (Erk) activation. Humanin responsiveness was subject to inter-individual heterogeneity and predominated in GSC with high expression levels of the IL6ST subunit. Mechanistically, humanin promoted the ATR-dependent DNA-repair machinery in GSCs via induction of the DNA-clamp component HUS1. A slice culture model containing human microglia and GSCs confirmed these observations. Humanin is absent from the mouse brain. Exogenous delivery of humanin into orthotopic GBM mouse models or over-expression specifically of a secreted humanin variant recapitulated the in vitro findings. Blockade of the trimeric interleukin receptor with the brain permeant, FDA-approved drug bazedoxifene blocked humanin-mediated chemoresistance in vitro and in vivo. Overall, we identified a clinically applicable compound together with a predictive marker to prevent chemoresistance in a human-specific GBM model.

The Immunological Epigenetic Landscape of the Human Life Trajectory.

Adaptive immunity changes over an individual's lifetime, maturing by adulthood and diminishing with old age. Epigenetic mechanisms involving DNA and histone methylation form the molecular basis of immunological memory during lymphocyte development. Monocytes alter their function to convey immune tolerance, yet the epigenetic influences at play remain to be fully understood in the context of lifespan. This study of a healthy genetically homogenous cohort of children, adults and seniors sought to decipher the epigenetic dynamics in B-lymphocytes and monocytes. Variable global cytosine methylation within retro-transposable LINE-1 repeats was noted in monocytes compared to B-lymphocytes across age groups. The expression of the human leukocyte antigen (HLA)-DQ&amp;nbsp;alpha chain gene HLA-DQA1*01 revealed significantly reduced levels in monocytes in all ages relative to B-lymphocytes, as well as between lifespan groups. High melting point analysis and bisulfite sequencing of the HLA-DQA1*01 promoter in monocytes highlighted variable cytosine methylation in children and seniors but greater stability at this locus in adults. Further epigenetic evaluation revealed higher histone lysine 27 trimethylation in monocytes from this adult group. Chromatin immunoprecipitation and RNA pulldown demonstrated association with a novel lncRNA TINA with structurally conserved similarities to the previously recognized epigenetic modifier PARTICLE. Seeking to interpret the epigenetic immunological landscape across three representative age groups, this study focused on HLA-DQA1*01 to expose cytosine and histone methylation alterations and their association with the non-coding transcriptome. Such insights unveil previously unknown complex epigenetic layers, orchestrating the strength and weakening of adaptive immunity with the progression of life.

Non-coding transcriptomic profiles in the sheep mammary gland during different lactation periods.

Sheep milk production is a dynamic and multifactorial trait regulated by diverse biological mechanisms. To improve the quality and production of sheep milk, it is necessary to understand the underlying non-coding transcriptomic mechanisms. In this study, ribonucleic acid-sequencing (RNA-seq) was used to profile the expression of microRNAs (miRNAs) and circular RNAs (circRNAs) in the sheep mammary gland at three key lactation time points (perinatal period, PP; early lactation, EL; and peak lactation, PL). A total of 2,369 novel circRNAs and 272 miRNAs were profiled, of which 348, 373, and 36 differentially expressed (DE) circRNAs and 30, 34, and 7 DE miRNAs were detected in the comparison of EL vs. PP, PL vs. PP, and PL vs. EL, respectively. A series of bioinformatics analyses including functional enrichment, machine learning prediction, and competing endogenous RNA (ceRNA) network analyses were conducted to identify subsets of the potential candidate miRNAs (e.g., oar_miR_148a, oar_miR_362, and oar_miR_432) and circRNAs (e.g., novel_circ_0011066, novel_circ_0010460, and novel_circ_0006589) involved in sheep mammary gland development. Taken together, this study offers a window into the dynamics of non-coding transcriptomes that occur during sheep lactation and may provide further insights into miRNA and circRNA that influence sheep mammary gland development.

Long noncoding RNA transcriptome analysis reveals novel lncRNAs in Morus alba 'Yu-711' response to drought stress.

Drought stress has been a key environmental factor affecting plant growth and development. The plant genome is capable of producing long noncoding RNAs (lncRNAs). To better understand white mulberry (Morus alba L.) drought response mechanism, we conducted a comparative transcriptome study comparing two treatments: drought-stressed (EG) and well-watered (CK) plants. A total of 674 differentially expressed lncRNAs (DElncRNAs) were identified. In addition, 782 differentially expressed messenger RNAs (DEmRNAs) were identified. We conducted Gene Ontology (GO) and KEGG enrichment analyses focusing on the differential lncRNAs cis-target genes. The target genes of the DElncRNAs were most significantly involved in the biosynthesis of secondary metabolites. Gene regulatory networks of the target genes involving DElncRNAs-mRNAs-DEmRNAs and DElncRNA-miRNA-DEmRNA were constructed. In the DElncRNAs-DEmRNAs network, 30 DEmRNAs involved in the biosynthesis of secondary metabolites are collocated with 46 DElncRNAs. The interaction between DElncRNAs and candidate genes was identified using LncTar. In summary, quantitative real-time polymerase chain reaction (qRT-PCR) validated nine candidate genes and seven target lncRNAs including those identified by LncTar. We predicted that the DElncRNAs-DEmRNAs might recruit microRNAs (miRNAs) to interact with gene regulatory networks under the drought stress response in mulberry. The findings will contribute to our understanding of the regulatory functions of lncRNAs under drought stress and will shed new light on the mulberry-drought stress interactions.

Regulation of lipid droplet (LD) formation in hepatocytes via regulation of SREBP1c by non-coding RNAs.

IntroductionIncreased de novo lipogenesis (DNL) is one of the key factors contributing to fat accumulation and non-alcoholic fatty liver disease (NAFLD). Among the critical transcription factors (TFs) regulating DNL is mTOR and its downstream lipogenic TF, SREBP1c. In recent years, it has been established that non-coding RNAs (ncRNAs) play role in both biological processes and disease pathogenesis. Our group has previously characterized microRNAs that can target and regulate the expression of both mTOR and SREBP1c. Accordingly, this study aimed to broaden our understanding of the role of ncRNAs in regulating the mTOR/SREBP1c axis to elucidate the role of the non-coding transcriptome in DNL and lipid droplet (LD) formation. Hence, short ncRNA, miR-615-5p, and long non-coding RNA (lncRNA), H19, were chosen as they were previously proven to target mTOR by our group and in the published literature, respectively.MethodologyHuh-7 cells were treated with 800 μM oleic acid (OA) to promote LD formation. Transfection of miR-615-5p mimics or H19 over-expression vectors was performed, followed by the measurement of their downstream targets, mTOR and SREBP, on the mRNA level by quantitative real-time PCR (qRT-PCR), and on the protein level by Western blot. To determine the functional impact of miR-615-5p and H19 on LD formation and triglyceride (TG) accumulation, post-transfection LDs were stained, imaged, and characterized, and TGs were extracted and quantified.ResultsmiR-615-5p was able to reduce mTOR and SREBP1c significantly on both the mRNA and protein levels compared to control cells, while H19 caused a reduction of both targets on the protein level only. Both miR-615-5p and H19 were able to significantly reduce the LD count and total area, as well as TG levels compared to control cells.ConclusionTo conclude, this study shows, for the first time, the impact of miR-615-5p and H19 on the mTOR/SREBP1c axis, and thus, their functional impact on LDs and TG accumulation. These findings might pave the way for using ncRNAs as potential therapeutic targets in the management of fatty liver.

Functional high-throughput screen identifies microRNAs that promote butyrate-induced death in colorectal cancer cells

The gut fermentation product butyrate displays anti-cancer properties in the human proximal colon, including the ability to inhibit proliferation and induce apoptosis in colorectal cancer (CRC) cells. A natural histone deacetylase inhibitor (HDACi), butyrate can alter histone acetylation patterns in CRC cells, and thereby regulate global gene expression, including the non-coding transcriptome and microRNAs (miRNAs). Dysregulated miRNA expression affects CRC development and progression; however, the interplay between miRNA activity and butyrate response remains to be elucidated. A high-throughput functional screen was employed to identify miRNAs that can act as enhancers of the anti-cancer properties of butyrate. Validation studies confirmed that several miRNAs, including miR-125b, miR-181a, miR-593, and miR-1227, enhanced apoptosis, decreased proliferation, and promoted cell-cycle arrest in the presence of butyrate. Pathway analyses of predicted miRNA target genes highlighted their likely involvement in critical cancer-related growth pathways, including WNT and PI3K signaling. Several cancer-associated miRNA targets, including TRIM29, COX2, PIK3R3, CCND1, MET, EEF2K, DVL3, and NUP62 were synergistically regulated by the combination of cognate miRNAs and butyrate. Overall, this study has exposed the potential of miRNAs to act as enhancers of the anti-cancer effects of HDAC inhibition and identifies specific miRNAs that might be exploited for therapeutic benefit.

The Role of Non-Coding RNAs in Glioma.

For decades, research in cancer biology has been focused on the protein-coding fraction of the human genome. However, with the discovery of non-coding RNAs (ncRNAs), it has become known that these entities not only function in numerous fundamental life processes such as growth, differentiation, and development, but also play critical roles in a wide spectrum of human diseases, including cancer. Dysregulated ncRNA expression is found to affect cancer initiation, progression, and therapy resistance, through transcriptional, post-transcriptional, or epigenetic processes in the cell. In this review, we focus on the recent development and advances in ncRNA biology that are pertinent to their role in glioma tumorigenesis and therapy response. Gliomas are common, and are the most aggressive type of primary tumors, which account for ~30% of central nervous system (CNS) tumors. Of these, glioblastoma (GBM), which are grade IV tumors, are the most lethal brain tumors. Only 5% of GBM patients survive beyond five years upon diagnosis. Hence, a deeper understanding of the cellular non-coding transcriptome might help identify biomarkers and therapeutic agents for a better treatment of glioma. Here, we delve into the functional roles of microRNA (miRNA), long non-coding RNA (lncRNA), and circular RNA (circRNA) in glioma tumorigenesis, discuss the function of their extracellular counterparts, and highlight their potential as biomarkers and therapeutic agents in glioma.

Molecular regulation of hypoxia through the lenses of noncoding RNAs and epitranscriptome.

Cells maintain homeostasis in response to environmental stress through specific cell stress responses. Hypoxic stress, well known to be associated with diverse solid tumors, is one of the main reasons for cancer-related mortality. Although cells can balance themselves well during hypoxic stress, the underlying molecular mechanisms are not well understood. The enhanced appreciation of diverse roles played by noncoding transcriptome and epigenome in recent years has brought to light the involvement of noncoding RNAs and epigenetic modifiers in hypoxic regulation. The emergence of techniques like deep sequencing has facilitated the identification of large numbers of long noncoding RNAs (lncRNAs) that are differentially regulated in various cancers. Similarly, proteomic studies have identified diverse epigenetic modifiers such as HATs, HDACs, DNMTs, polycomb groups of proteins, and their possible roles in the regulation of hypoxia. The crosstalk between lncRNAs and epigenetic modifiers play a pivotal role in hypoxia-induced cancer initiation and progression. Besides the lncRNAs, several other noncoding RNAs like circular RNAs, miRNAs, and so forth are also expressed during hypoxic conditions. Hypoxia has a profound effect on the expression of noncoding RNAs and epigenetic modifiers. Conversely, noncoding RNAs/epigenetic modifies can regulate the hypoxia signaling axis by modulating the stability of the hypoxia-inducible factors (HIFs). The focus of this review is to illustrate the molecular orchestration underlying hypoxia biology, especially in cancers, which can help in identifying promising therapeutic targets in hypoxia-induced cancers. This article is categorized under: RNA Turnover and Surveillance > Regulation of RNA Stability RNA in Disease and Development > RNA in Disease RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry.

Mutant KRAS regulates transposable element RNA and innate immunity via KRAB zinc-finger genes.

SUMMARYRAS genes are the most frequently mutated oncogenes in cancer, yet the effects of oncogenic RAS signaling on the noncoding transcriptome remain unclear. We analyzed the transcriptomes of human airway and bronchial epithelial cells transformed with mutant KRAS to define the landscape of KRAS-regulated noncoding RNAs. We find that oncogenic KRAS signaling upregulates noncoding transcripts throughout the genome, many of which arise from transposable elements (TEs). These TE RNAs exhibit differential expression, are preferentially released in extracellular vesicles, and are regulated by KRAB zinc-finger (KZNF) genes, which are broadly downregulated in mutant KRAS cells and lung adenocarcinomas in vivo. Moreover, mutant KRAS induces an intrinsic IFN-stimulated gene (ISG) signature that is often seen across many different cancers. Our results indicate that mutant KRAS remodels the repetitive noncoding transcriptome, demonstrating the broad scope of intracellular and extracellular RNAs regulated by this oncogenic signaling pathway.

Bioinformatics Analysis of Key Genes and miRNAs Associated with Small-cell Lung Cancer

Small-cell lung cancer (SCLC) is a type of lung cancer, accounting for approximately 15% of lung cancers. SCLC is significantly associated with early recurrence, metastasis, and poor prognosis, with a 2-year survival rate after therapy, which is &lt;5%. The molecular mechanisms underlying the development of SCLC are not clear. In this study, we identified SCLC-specific biomarkers by evaluating the differential expression of mRNA and microRNA (miRNA) profiles in SCLC tissue compared with normal lung tissue. A non-coding RNA sequence dataset (GSE19945) and transcriptome sequencing dataset (GSE6044) were downloaded from the GEO database. We identified 445 DEGs (differentially expressed genes) and 128 DE-miRNAs (differentially expressed miRNAs) using the GEO2R tool of the GEO and R limma software package. Furthermore, using the KEGG database, we identified 15 enrichment pathways, mostly associated with DNA replication, cell cycle, and oocyte meiosis mismatch repair, and the GO function was considerably enriched for 26 items. To investigate the molecular processes of key signaling pathways and cellular activity in SCLC, we used Cytoscape software to construct protein-protein interaction (PPI) networks. Using miRNAWalk, we identified 598 target genes of the 1380 miRNAs and constructed miRNA target networks. In addition, we identified eighteen overlapping genes that are regulated by 28 different miRNAs. The identified hub genes are important because they may be used as biomarkers for prognosis, diagnosis, and therapeutic target for SCLC.