SSS-test: a novel test for detecting positive selection on RNA secondary structure

Deep Conservation and Unexpected Evolutionary History of Neighboring lncRNAs MALAT1 and NEAT1

Long non-coding RNAs (lncRNAs) have begun to receive overdue attention for their regulatory roles in gene expression and other cellular processes. Although most lncRNAs are lowly expressed and tissue-specific, notable exceptions include MALAT1 and its genomic neighbor NEAT1, two highly and ubiquitously expressed oncogenes with roles in transcriptional regulation and RNA splicing. Previous studies have suggested that NEAT1 is found only in mammals, while MALAT1 is present in all gnathostomes (jawed vertebrates) except birds. Here we show that these assertions are incomplete, likely due to the challenges associated with properly identifying these two lncRNAs. Using phylogenetic analysis and structure-aware annotation of publicly available genomic and RNA-seq coverage data, we show that NEAT1 is a common feature of tetrapod genomes except birds and squamates. Conversely, we identify MALAT1 in representative species of all major gnathostome clades, including birds. Our in-depth examination of MALAT1, NEAT1, and their genomic context in a wide range of vertebrate species allows us to reconstruct the series of events that led to the formation of the locus containing these genes in taxa from cartilaginous fish to mammals. This evolutionary history includes the independent loss of NEAT1 in birds and squamates, since NEAT1 is found in the closest living relatives of both clades (crocodilians and tuataras, respectively). These data clarify the origins and relationships of MALAT1 and NEAT1 and highlight an opportunity to study the change and continuity in lncRNA structure and function over deep evolutionary time.

Comparative RNA Genomics

Over the last quarter of a century it has become clear that RNA is much more than just a boring intermediate in protein expression. Ancient RNAs still appear in the core information metabolism and comprise a surprisingly large component in bacterial gene regulation. A common theme with these types of mostly small RNAs is their reliance of conserved secondary structures. Large-scale sequencing projects, on the other hand, have profoundly changed our understanding of eukaryotic genomes. Pervasively transcribed, they give rise to a plethora of large and evolutionarily extremely flexible non-coding RNAs that exert a vastly diverse array of molecule functions. In this chapter we provide a-necessarily incomplete-overview of the current state of comparative analysis of non-coding RNAs, emphasizing computational approaches as a means to gain a global picture of the modern RNA world.

Computational Prediction of RNA-RNA Interactions between Small RNA Tracks from Betacoronavirus Nonstructural Protein 3 and Neurotrophin Genes during Infection of an Epithelial Lung Cancer Cell Line: Potential Role of Novel Small Regulatory RNA

Whether RNA-RNA interactions of cytoplasmic RNA viruses, such as Betacoronavirus, might end in the biogenesis of putative virus-derived small RNAs as miRNA-like molecules has been controversial. Even more, whether RNA-RNA interactions of wild animal viruses may act as virus-derived small RNAs is unknown. Here, we address these issues in four ways. First, we use conserved RNA structures undergoing negative selection in the genomes of SARS-CoV, MERS-CoV, and SARS-CoV-2 circulating in different bat species, intermediate animals, and human hosts. Second, a systematic literature review was conducted to identify Betacoronavirus-targeting hsa-miRNAs involved in lung cell infection. Third, we employed sophisticated long-range RNA-RNA interactions to refine the seed sequence homology of hsa-miRNAs with conserved RNA structures. Fourth, we used high-throughput RNA sequencing of a Betacoronavirus-infected epithelial lung cancer cell line (Calu-3) to validate the results. We proposed nine potential virus-derived small RNAs: two vsRNAs in SARS-CoV (Bats: SB-vsRNA-ORF1a-3p; SB-vsRNA-S-5p), one vsRNA in MERS-CoV (Bats: MB-vsRNA-ORF1b-3p), and six vsRNAs in SARS-CoV-2 (Bats: S2B-vsRNA-ORF1a-5p; intermediate animals: S2I-vsRNA-ORF1a-5p; and humans: S2H-vsRNA-ORF1a-5p, S2H-vsRNA-ORF1a-3p, S2H-vsRNA-ORF1b-3p, S2H-vsRNA-ORF3a-3p), mainly encoded by nonstructural protein 3. Notably, Betacoronavirus-derived small RNAs targeted 74 differentially expressed genes in infected human cells, of which 55 upregulate the molecular mechanisms underlying acute respiratory distress syndrome (ARDS), and the 19 downregulated genes might be implicated in neurotrophin signaling impairment. These results reveal a novel small RNA-based regulatory mechanism involved in neuropathogenesis that must be further studied to validate its therapeutic use.

Evolutionary Conservation of RNA Secondary Structure

Noncoding RNAs, ncRNAs, naturally fold into structures, which allow them to perform their functions in the cell. Evolutionarily close species share structures and functions. This occurs because of shared selective pressures, resulting in conserved groups. Previous efforts in finding functional RNAs have been made in detecting conserved structures in genomes or alignments. It may occur that, within a conserved group, species-specific structures arise after species split due to positive selection. Detecting positive selection in ncRNAs is a hard problem in biology as well as bioinformatics. To detect positive selection, one should find species-specific structures within a conserved set. This chapter provides protocols to detect and analyze positive selection in ncRNA structures with the SSS-test and other free software.

RNA structure-altering mutations underlying positive selection on Spike protein reveal novel putative signatures to trace crossing host-species barriers in Betacoronavirus

Similar to other RNA viruses, the emergence of Betacoronavirus relies on cross-species viral transmission, which requires careful health surveillance monitoring of protein-coding information as well as genome-wide analysis. Although the evolutionary jump from natural reservoirs to humans may be mainly traced-back by studying the effect that hotspot mutations have on viral proteins, it is largely unexplored if other impacts might emerge on the structured RNA genome of Betacoronavirus. In this survey, the protein-coding and viral genome architecture were simultaneously studied to uncover novel insights into cross-species horizontal transmission events. We analysed 1,252,952 viral genomes of SARS-CoV, MERS-CoV, and SARS-CoV-2 distributed across the world in bats, intermediate animals, and humans to build a new landscape of changes in the RNA viral genome. Phylogenetic analyses suggest that bat viruses are the most closely related to the time of most recent common ancestor of Betacoronavirus, and missense mutations in viral proteins, mainly in the S protein S1 subunit: SARS-CoV (G > T; A577S); MERS-CoV (C > T; S746R and C > T; N762A); and SARS-CoV-2 (A > G; D614G) appear to have driven viral diversification. We also found that codon sites under positive selection on S protein overlap with non-compensatory mutations that disrupt secondary RNA structures in the RNA genome complement. These findings provide pivotal factors that might be underlying the eventual jumping the species barrier from bats to intermediate hosts. Lastly, we discovered that nearly half of the Betacoronavirus genomes carry highly conserved RNA structures, and more than 90% of these RNA structures show negative selection signals, suggesting essential functions in the biology of Betacoronavirus that have not been investigated to date. Further research is needed on negatively selected RNA structures to scan for emerging functions like the potential of coding virus-derived small RNAs and to develop new candidate antiviral therapeutic strategies.

Knockdown of PDX1 enhances the osteogenic differentiation of ADSCs partly via activation of the PI3K/Akt signaling pathway

Background

Osteoporosis (OP) is a systemic bone disease manifested as low bone mass, destruction of bone microstructure, increased bone fragility and fracture risk. The purpose of this study was to explore the role and mechanism of PDX1 for osteogenic differentiation of adipose derived stem cells (ADSCs).

Methods

GSE37329 dataset was retrieved from NCBI Gene Expression Omnibus (GEO) database and performed bioinformatic analyses. ADSCs were incubated with normal medium, osteogenic induction medium (OIM) and OIM+si-PDX1. Then, alkaline phosphatase (ALP) staining and Alizarin Red Staining (ARS) were performed to assess the role of PDX1 for osteogenesis of ADSCs. PI3K inhibitor, LY294002 was then added to further explore the mechanism of PDX1 for osteogenic differentiation of ADSCs. Western blot assay was used to assess the osteogenic-related markers. Graphpad software was used to perform statistically analysis.

Results

A total of 285 DEGs were obtained from analysis of the dataset GSE37329, of which 145 were upregulated and 140 were downregulated genes. These differentially expressed genes mainly enriched in cell differentiation and PI3K/Akt signaling pathway. Moreover, PDX1 was decreased in osteogenic induced ADSCs. Knockdown of PDX1 significantly increased osteogenic differentiation capacity and p-PI3K and p-Akt protein levels. Administration with LY294002 could partially reversed the promotion effects of si-PDX1.

Conclusion

In conclusion, knockdown of PDX1 promotes osteogenic differentiation of ADSCs through the PI3K/Akt signaling pathway.

OUP accepted manuscript

Accelerated evolution of any portion of the genome is of significant interest, potentially signaling positive selection of phenotypic traits and adaptation. Accelerated evolution remains understudied for structured RNAs, despite the fact that an RNA’s structure is often key to its function. RNA structures are typically characterized by compensatory (structure-preserving) basepair changes that are unexpected given the underlying sequence variation, i.e., they have evolved through negative selection on structure. We address the question of how fast the primary sequence of an RNA can change through evolution while conserving its structure. Specifically, we consider predicted and known structures in vertebrate genomes. After careful control of false discovery rates, we obtain 13 de novo structures (and three known Rfam structures) that we predict to have rapidly evolving sequences—defined as structures where the primary sequences of human and mouse have diverged at least twice as fast (1.5 times for Rfam) as nearby neutrally evolving sequences. Two of the three known structures function in translation inhibition related to infection and immune response. We conclude that rapid sequence divergence does not preclude RNA structure conservation in vertebrates, although these events are relatively rare.

LINC02288 promotes chondrocyte apoptosis and inflammation through miR-374a-3p targeting RTN3

BACKGROUND

Dysregulation of long non-coding RNAs (lncRNAs) is related to the occurrence of osteoarthritis (OA). In the present study, we explored the role of LINC02288 and its regulatory mechanism in OA development.

METHODS

GSE113825 was obtained from Gene Expression Omnibus (GEO) database and analyzed to identify the differentially expressed lncRNAs in OA. Gene enrichment analyses and Kyoto Encyclopedia of Genes and Genomes biological process analysis were performed through Metascape (http://metascape.org/gp). The interactions among LINC02288, miR-374a-3p and RTN3 were determined using RNA immunoprecipitation (RIP) assays and dual luciferase reporter assays. Chondrocyte apoptosis was examined using flow cytometry. Western blot assays were conducted to assess the pro-apoptotic and anti-apoptotic markers.

RESULTS

We identified a total of 4,491 differentially expressed lncRNAs. We focused on LINC02288 as the top-ranked up-regulated lncRNA in OA as indicated by a significant p-value. LINC02288 was significantly up-regulated, which was further verified by a real-time polymerase chain reaction. Down-regulation of LINC02288 significantly reduced the apoptosis of OA chondrocytes induced by interleukin-1β and the production of pro-inflammatory cytokines. These effects were further verified in an OA rat model. An RIP assay and dual luciferase assay further confirmed that LINC02288 served as a sponge of miR-374a-3p. Moreover, the overexpression of RTN3 could partially reverse the effects of LINC02288 knockdown, mediating inhibitory effects on chondrocyte apoptosis and the inflammatory response. Down-regulation of LINC02288 alleviated OA development in an in vivo OA animal model.

CONCLUSIONS

Our findings indicate that LINC02288 contributes to OA progression by targeting the miR-374a-3p/RTN3 axis, which might provide a promising molecular therapy strategy for OA.

Analysis of Long Non-Coding RNA in Cryptosporidium parvum Reveals Significant Stage-Specific Antisense Transcription

Cryptosporidium is a protist parasite that has been identified as the second leading cause of moderate to severe diarrhea in children younger than two and a significant cause of mortality worldwide. Cryptosporidium has a complex, obligate, intracellular but extra cytoplasmic lifecycle in a single host. How genes are regulated in this parasite remains largely unknown. Long non-coding RNAs (lncRNAs) play critical regulatory roles, including gene expression across a broad range of organisms. Cryptosporidium lncRNAs have been reported to enter the host cell nucleus and affect the host response. However, no systematic study of lncRNAs in Cryptosporidium has been conducted to identify additional lncRNAs. In this study, we analyzed a C. parvum in vitro strand-specific RNA-seq developmental time series covering both asexual and sexual stages to identify lncRNAs associated with parasite development. In total, we identified 396 novel lncRNAs, mostly antisense, with 86% being differentially expressed. Surprisingly, nearly 10% of annotated mRNAs have an antisense transcript. lncRNAs occur most often at the 3' end of their corresponding sense mRNA. Putative lncRNA regulatory regions were identified and many appear to encode bidirectional promoters. A positive correlation between lncRNA and upstream mRNA expression was observed. Evolutionary conservation and expression of lncRNA candidates was observed between C. parvum, C. hominis and C. baileyi. Ten C. parvum protein-encoding genes with antisense transcripts have P. falciparum orthologs that also have antisense transcripts. Three C. parvum lncRNAs with exceptional properties (e.g., intron splicing) were experimentally validated using RT-PCR and RT-qPCR. This initial characterization of the C. parvum non-coding transcriptome facilitates further investigations into the roles of lncRNAs in parasite development and host-pathogen interactions.

The role of lncRNAs in innate immunity and inflammation

The innate immune system relies on a germ-line-encoded repertoire of pattern recognition receptors (PRRs), activated by deeply conserved pathogen signatures, such as bacterial cell wall components or foreign nucleic acids. To enable effective defence against invading pathogens and prevent from deleterious inflammation, PRR-driven immune responses are tightly controlled by a dense network of nuclear and cytoplasmic regulators. Long non-coding RNAs (lncRNAs) are increasingly recognized as important components of these regulatory circuitries, providing positive and negative control of PRR-induced innate immune responses. The present review provides an overview of the presently known roles of lncRNAs in human and murine innate antiviral and antibacterial immunity. The emerging roles in host defence and inflammation suggest that further mechanistic insights into the cellular functions of lncRNAs will decisively advance our molecular understanding of immune-associated diseases and open new avenues for therapeutic intervention.

Selection Pressures on RNA Sequences and Structures

With the discovery of increasingly more functional noncoding RNAs (ncRNAs), it becomes eminent to more strongly consider them as important players during species evolution. Although tests for negative selection of ncRNAs already exist since the beginning of this century, the SSS-test is the first one for also investigating positive selection. When analyzing selection in ncRNAs, it should be taken into account that selection pressures can independently act on sequence and structure. We applied the SSS-test to explore the evolution of ncRNAs in primates and identified more than 100 long noncoding RNAs (lncRNAs) that might evolve under positive selection in humans. With this test, it is now possible to more thoroughly include ncRNAs into evolutionary studies.