Transcriptomic analyses reveal groups of co-expressed, syntenic lncRNAs in four species of the genus Caenorhabditis

Long non-coding RNAs mediate fish gene expression in response to ocean acidification

The majority of the transcribed genome does not have coding potential but these non-coding transcripts play crucial roles in transcriptional and post-transcriptional regulation of protein-coding genes. Regulation of gene expression is important in shaping an organism's response to environmental changes, ultimately impacting their survival and persistence as population or species face global change. However, the roles of long non-coding RNAs (lncRNAs), when confronted with environmental changes, remain largely unclear. To explore the potential role of lncRNAs in fish exposed to ocean acidification (OA), we analyzed publicly available brain RNA-seq data from a coral reef fish Acanthochromis polyacanthus. We annotated the lncRNAs in its genome and examined the expression changes of intergenic lncRNAs (lincRNAs) between A. polyacanthus samples from a natural CO2 seep and a nearby control site. We identified 4728 lncRNAs, including 3272 lincRNAs in this species. Remarkably, 93.03% of these lincRNAs were species-specific. Among the 125 highly expressed lincRNAs and 403 differentially expressed lincRNAs in response to elevated CO2, we observed that lincRNAs were either neighboring or potentially trans-regulating differentially expressed coding genes associated with pH regulation, neural signal transduction, and ion transport, which are known to be important in the response to OA in fish. In summary, lncRNAs may facilitate fish acclimation and mediate the responses of fish to OA by modulating the expression of crucial coding genes, which offers insight into the regulatory mechanisms underlying fish responses to environmental changes.

Functional evolutionary convergence of long noncoding RNAs involved in embryonic development

Long noncoding RNAs have been identified in most vertebrates, but the functional characterization of these molecules is challenging, mainly due to the lack of linear sequence homology between species. In this work, we aimed to find functional evolutionary convergent lncRNAs involved in development by screening of k-mer content (nonlinear similarity) and secondary structure-based approaches combining in silico, in vitro and in vivo validation analysis. From the Madagascar gecko genes, we have found a non-orthologous lncRNA with a similar k-mer content and structurally concordant with the human lncRNA EVX1AS. Analysis of function-related characteristics together with locus-specific targeting of human EVX1AS and gecko EVX1AS-like (i.e., CRISPR Display) in human neuroepithelial cells and chicken mesencephalon have confirmed that gecko EVX1AS-like lncRNA mimics human EVX1AS function and induces EVX1 expression independently of the target species. Our data shows functional convergence of non-homologous lncRNAs and presents a useful approach for the definition and manipulation of lncRNA function within different model organisms.

A hybrid approach to assess the structural impact of long noncoding RNA mutations uncovers key NEAT1 interactions in colorectal cancer

Long noncoding RNAs (lncRNAs) are emerging players in cancer and they entail potential as prognostic biomarkers or therapeutic targets. Earlier studies have identified somatic mutations in lncRNAs that are associated with tumor relapse after therapy, but the underlying mechanisms behind these associations remain unknown. Given the relevance of secondary structure for the function of some lncRNAs, some of these mutations may have a functional impact through structural disturbance. Here, we examined the potential structural and functional impact of a novel A > G point mutation in NEAT1 that has been recurrently observed in tumors of colorectal cancer patients experiencing relapse after treatment. Here, we used the nextPARS structural probing approach to provide first empirical evidence that this mutation alters NEAT1 structure. We further evaluated the potential effects of this structural alteration using computational tools and found that this mutation likely alters the binding propensities of several NEAT1-interacting miRNAs. Differential expression analysis on these miRNA networks shows upregulation of Vimentin, consistent with previous findings. We propose a hybrid pipeline that can be used to explore the potential functional effects of lncRNA somatic mutations.

The Conservation of Long Intergenic Non-Coding RNAs and Their Response to Verticillium dahliae Infection in Cotton

Long intergenic non-coding RNAs (lincRNAs) have been demonstrated to be vital regulators of diverse biological processes in both animals and plants. While many lincRNAs have been identified in cotton, we still know little about the repositories and conservativeness of lincRNAs in different cotton species or about their role in responding to biotic stresses. Here, by using publicly available RNA-seq datasets from diverse sources, including experiments of Verticillium dahliae (Vd) infection, we identified 24,425 and 17,713 lincRNAs, respectively, in Gossypium hirsutum (Ghr) and G. barbadense (Gba), the two cultivated allotetraploid cotton species, and 6933 and 5911 lincRNAs, respectively, in G. arboreum (Gar) and G. raimondii (Gra), the two extant diploid progenitors of the allotetraploid cotton. While closely related subgenomes, such as Ghr_At and Gba_At, tend to have more conserved lincRNAs, most lincRNAs are species-specific. The majority of the synthetic and transcribed lincRNAs (78.2%) have a one-to-one orthologous relationship between different (sub)genomes, although a few of them (0.7%) are retained in all (sub)genomes of the four species. The Vd responsiveness of lincRNAs seems to be positively associated with their conservation level. The major functionalities of the Vd-responsive lincRNAs seem to be largely conserved amongst Gra, Ghr, and Gba. Many Vd-responsive Ghr-lincRNAs overlap with Vd-responsive QTL, and several lincRNAs were predicted to be endogenous target mimicries of miR482/2118, with a pair being highly conserved between Ghr and Gba. On top of the confirmation of the feature characteristics of the lincRNAs previously reported in cotton and other species, our study provided new insights into the conservativeness and divergence of lincRNAs during cotton evolution and into the relationship between the conservativeness and Vd responsiveness of lincRNAs. The study also identified candidate lincRNAs with a potential role in disease response for functional characterization.

Genomic and functional conservation of lncRNAs: lessons from flies

Over the last decade, the increasing interest in long non-coding RNAs (lncRNAs) has led to the discovery of these transcripts in multiple organisms. LncRNAs tend to be specifically, and often lowly, expressed in certain tissues, cell types and biological contexts. Although lncRNAs participate in the regulation of a wide variety of biological processes, including development and disease, most of their functions and mechanisms of action remain unknown. Poor conservation of the DNA sequences encoding for these transcripts makes the identification of lncRNAs orthologues among different species very challenging, especially between evolutionarily distant species such as flies and humans or mice. However, the functions of lncRNAs are unexpectedly preserved among different species supporting the idea that conservation occurs beyond DNA sequences and reinforcing the potential of characterising lncRNAs in animal models. In this review, we describe the features and roles of lncRNAs in the fruit fly Drosophila melanogaster, focusing on genomic and functional comparisons with human and mouse lncRNAs. We also discuss the current state of advances and limitations in the study of lncRNA conservation and future perspectives.

The evolutionary landscape and expression pattern of plant lincRNAs

Long intergenic non-coding RNAs (lincRNAs) are important regulators of cellular processes, including development and stress response. Many lincRNAs have been bioinformatically identified in plants, but their evolutionary dynamics and expression characteristics are still elusive. Here, we systematically identified thousands of lincRNAs in 26 plant species, including 6 non-flowering plants, investigated the conservation of the identified lincRNAs in different levels of plant lineages based on sequence and/or synteny homology and explored characteristics of the conserved lincRNAs during plant evolution and their co-expression relationship with protein-coding genes (PCGs). In addition to confirmation of the features well documented in literature for lincRNAs, such as species-specific, fewer exons, tissue-specific expression patterns and less abundantly expressed, we revealed that histone modification signals and/or binding sites of transcription factors were enriched in the conserved lincRNAs, implying their biological functionalities, as demonstrated by identifying conserved lincRNAs related to flower development in both the Brassicaceae and grass families and ancient lincRNAs potentially functioning in meristem development of non-flowering plants. Compared to PCGs, lincRNAs are more likely to be associated with transposable elements (TEs), but with different characteristics in different evolutionary lineages, for instance, the types of TEs and the variable level of association in lincRNAs with different conservativeness. Together, these results provide a comprehensive view on the evolutionary landscape of plant lincRNAs and shed new insights on the conservation and functionality of plant lincRNAs.

The long non-coding RNA landscape of Candida yeast pathogens

Long non-coding RNAs (lncRNAs) constitute a poorly studied class of transcripts with emerging roles in key cellular processes. Despite efforts to characterize lncRNAs across a wide range of species, these molecules remain largely unexplored in most eukaryotic microbes, including yeast pathogens of the Candida clade. Here, we analyze thousands of publicly available sequencing datasets to infer and characterize the lncRNA repertoires of five major Candida pathogens: Candida albicans, Candida tropicalis, Candida parapsilosis, Candida auris and Candida glabrata. Our results indicate that genomes of these species encode hundreds of lncRNAs that show levels of evolutionary constraint intermediate between those of intergenic genomic regions and protein-coding genes. Despite their low sequence conservation across the studied species, some lncRNAs are syntenic and are enriched in shared sequence motifs. We find co-expression of lncRNAs with certain protein-coding transcripts, hinting at potential functional associations. Finally, we identify lncRNAs that are differentially expressed during infection of human epithelial cells for four of the studied species. Our comprehensive bioinformatic analyses of Candida lncRNAs pave the way for future functional characterization of these transcripts.

Long non-coding RNAs (lncRNAs) play roles in key cellular processes, but remain largely unexplored in fungal pathogens such as Candida. Here, Hovhannisyan and Gabaldón analyze thousands of sequencing datasets to infer and characterize the lncRNA repertoires of five Candida species, paving the way for their future functional characterization.

Comparative genomics in the search for conserved long noncoding RNAs

Long noncoding RNAs (lncRNAs) have emerged as prominent regulators of gene expression in eukaryotes. The identification of lncRNA orthologs is essential in efforts to decipher their roles across model organisms, as homologous genes tend to have similar molecular and biological functions. The relatively high sequence plasticity of lncRNA genes compared with protein-coding genes, makes the identification of their orthologs a challenging task. This is why comparative genomics of lncRNAs requires the development of specific and, sometimes, complex approaches. Here, we briefly review current advancements and challenges associated with four levels of lncRNA conservation: genomic sequences, splicing signals, secondary structures and syntenic transcription.

Ten Years of Collaborative Progress in the Quest for Orthologs

Accurate determination of the evolutionary relationships between genes is a foundational challenge in biology. Homology-evolutionary relatedness-is in many cases readily determined based on sequence similarity analysis. By contrast, whether or not two genes directly descended from a common ancestor by a speciation event (orthologs) or duplication event (paralogs) is more challenging, yet provides critical information on the history of a gene. Since 2009, this task has been the focus of the Quest for Orthologs (QFO) Consortium. The sixth QFO meeting took place in Okazaki, Japan in conjunction with the 67th National Institute for Basic Biology conference. Here, we report recent advances, applications, and oncoming challenges that were discussed during the conference. Steady progress has been made toward standardization and scalability of new and existing tools. A feature of the conference was the presentation of a panel of accessible tools for phylogenetic profiling and several developments to bring orthology beyond the gene unit-from domains to networks. This meeting brought into light several challenges to come: leveraging orthology computations to get the most of the incoming avalanche of genomic data, integrating orthology from domain to biological network levels, building better gene models, and adapting orthology approaches to the broad evolutionary and genomic diversity recognized in different forms of life and viruses.

The landscape of lncRNAs in Cydia pomonella provides insights into their signatures and potential roles in transcriptional regulation

Background

Long noncoding RNAs (lncRNAs) have emerged as an important class of transcriptional regulators in cellular processes. The past decades have witnessed great progress in lncRNA studies in a variety of organisms. The codling moth (Cydia pomonella L.) is an important invasive insect in China. However, the functional impact of lncRNAs in this insect remains unclear. In this study, an atlas of codling moth lncRNAs was constructed based on publicly available RNA-seq datasets.

Results

In total, 9875 lncRNA transcripts encoded by 9161 loci were identified in the codling moth. As expected, the lncRNAs exhibited shorter transcript lengths, lower GC contents, and lower expression levels than protein-coding genes (PCGs). Additionally, the lncRNAs were more likely to show tissue-specific expression patterns than PCGs. Interestingly, a substantial fraction of the lncRNAs showed a testis-biased expression pattern. Additionally, conservation analysis indicated that lncRNA sequences were weakly conserved across insect species, though additional lncRNAs with homologous relationships could be identified based on synteny, suggesting that synteny could be a more reliable approach for the cross-species comparison of lncRNAs. Furthermore, the correlation analysis of lncRNAs with neighbouring PCGs indicated a stronger correlation between them, suggesting potential cis-acting roles of these lncRNAs in the regulation of gene expression.

Conclusions

Taken together, our work provides a valuable resource for the comparative and functional study of lncRNAs, which will facilitate the understanding of their mechanistic roles in transcriptional regulation.

Single-molecule real-time sequencing of the full-length transcriptome of loquat under low-temperature stress

In this study, third-generation full-length (FL) transcriptome sequencing was performed of loquat using single-molecule real-time(SMRT) sequencing from the pooled cDNA of embryos of young loquat fruit under different low temperatures (three biological replicates for treatments of 1°C, -1°C, and -3°C, for 12 h or 24 h) and the control group(three biological replicates for treatments of room temperature), Illumina sequencing was used to correct FL transcriptome sequences. A total of 3 PacBio Iso-Seq libraries (1–2 kb, 2–3 kb and 3–6 kb) and 21 Illumina transcriptome libraries were constructed, a total of 13.41 Gb of clean reads were generated, which included 215,636 reads of insert (ROIs) and 121,654 FL, non-chimaric (FLNC) reads. Transcript clustering analysis of the FLNC reads revealed 76,586 consensus isoforms, and a total of 12,520 high-quality transcript sequences corrected with non-FL sequences were used for subsequent analysis. After the redundant reads were removed, 38,435 transcripts were obtained. A total of 27,905 coding DNA sequences (CDSs) were identified, and 407 long non-coding RNAs (lncRNAs) were ultimately predicted. Additionally, 24,832 simple sequence repeats (SSRs) were identified, and a total of 1,295 alternative splicing (AS) events were predicted. Furthermore, 37,993 transcripts were annotated in eight functional databases. This is the first study to perform SMRT sequencing of the FL transcriptome of loquat. The obtained transcriptomic data are conducive for further exploration of the mechanism of loquat freezing injury and thus serve as an important theoretical basis for generating new loquat material and for identifying new ways to improve loquat cold resistance.

Dynamic Expression of Long Non-Coding RNAs Throughout Parasite Sexual and Neural Maturation in Schistosoma japonicum

Schistosoma japonicum is a flatworm that causes schistosomiasis, a neglected tropical disease. S. japonicum RNA-Seq analyses has been previously reported in the literature on females and males obtained during sexual maturation from 14 to 28 days post-infection in mouse, resulting in the identification of protein-coding genes and pathways, whose expression levels were related to sexual development. However, this work did not include an analysis of long non-coding RNAs (lncRNAs). Here, we applied a pipeline to identify and annotate lncRNAs in 66 S. japonicum RNA-Seq publicly available libraries, from different life-cycle stages. We also performed co-expression analyses to find stage-specific lncRNAs possibly related to sexual maturation. We identified 12,291 S. japonicum expressed lncRNAs. Sequence similarity search and synteny conservation indicated that some 14% of S. japonicum intergenic lncRNAs have synteny conservation with S. mansoni intergenic lncRNAs. Co-expression analyses showed that lncRNAs and protein-coding genes in S. japonicum males and females have a dynamic co-expression throughout sexual maturation, showing differential expression between the sexes; the protein-coding genes were related to the nervous system development, lipid and drug metabolism, and overall parasite survival. Co-expression pattern suggests that lncRNAs possibly regulate these processes or are regulated by the same activation program as that of protein-coding genes.

Genomic and Transcriptomic Analysis Identified Novel Putative Cassava lncRNAs Involved in Cold and Drought Stress

Long non-coding RNAs (lncRNAs) play important roles in the regulation of complex cellular processes, including transcriptional and post-transcriptional regulation of gene expression relevant for development and stress response, among others. Compared to other important crops, there is limited knowledge of cassava lncRNAs and their roles in abiotic stress adaptation. In this study, we performed a genome-wide study of ncRNAs in cassava, integrating genomics- and transcriptomics-based approaches. In total, 56,840 putative ncRNAs were identified, and approximately half the number were verified using expression data or previously known ncRNAs. Among these were 2229 potential novel lncRNA transcripts with unmatched sequences, 250 of which were differentially expressed in cold or drought conditions, relative to controls. We showed that lncRNAs might be involved in post-transcriptional regulation of stress-induced transcription factors (TFs) such as zinc-finger, WRKY, and nuclear factor Y gene families. These findings deepened our knowledge of cassava lncRNAs and shed light on their stress-responsive roles.

Microsyntenic Clusters Reveal Conservation of lncRNAs in Chordates Despite Absence of Sequence Conservation

Homologous long non-coding RNAs (lncRNAs) are elusive to identify by sequence similarity due to their fast-evolutionary rate. Here we develop LincOFinder, a pipeline that finds conserved intergenic lncRNAs (lincRNAs) between distant related species by means of microsynteny analyses. Using this tool, we have identified 16 bona fide homologous lincRNAs between the amphioxus and human genomes. We characterized and compared in amphioxus and Xenopus the expression domain of one of them, Hotairm1, located in the anterior part of the Hox cluster. In addition, we analyzed the function of this lincRNA in Xenopus, showing that its disruption produces a severe headless phenotype, most probably by interfering with the regulation of the Hox cluster. Our results strongly suggest that this lincRNA has probably been regulating the Hox cluster since the early origin of chordates. Our work pioneers the use of syntenic searches to identify non-coding genes over long evolutionary distances and helps to further understand lncRNA evolution.