A subset of conserved mammalian long non-coding RNAs are fossils of ancestral protein-coding genes

Syntenic lncRNA locus exhibits DNA regulatory functions with sequence evolution

Syntenic long non-coding RNAs (lncRNAs) often show limited sequence conservation across species, prompting concern in the field. This study delves into functional signatures of syntenic lncRNAs between humans and zebrafish. Syntenic lncRNAs are highly expressed in zebrafish, with ∼90 % located near protein-coding genes, either in sense or antisense orientation. During early zebrafish development and in human embryonic stem cells (H1-hESC), syntenic lncRNA loci are enriched with cis-regulatory repressor signatures, influencing the expression of development-associated genes. In later zebrafish developmental stages and specific human cell lines, these syntenic lncRNA loci function as enhancers or transcription start sites (TSS) for protein-coding genes. Analysis of transposable elements (TEs) in syntenic lncRNA sequences revealed intriguing patterns: human lncRNAs are enriched in simple repeat elements, while their zebrafish counterparts show enrichment in LTR elements. This sequence evolution likely arises from post-rearrangement mutations that enhance DNA elements or cis-regulatory functions. It may also contribute to vertebrate innovation by creating novel transcription factor binding sites within the locus. This study highlights the conserved functionality of syntenic lncRNA loci through DNA elements, emphasizing their conserved roles across species despite sequence divergence.

Enriched atlas of lncRNA and protein-coding genes for the GRCg7b chicken assembly and its functional annotation across 47 tissues

Gene atlases for livestock are steadily improving thanks to new genome assemblies and new expression data improving the gene annotation. However, gene content varies across databases due to differences in RNA sequencing data and bioinformatics pipelines, especially for long non-coding RNAs (lncRNAs) which have higher tissue and developmental specificity and are harder to consistently identify compared to protein coding genes (PCGs). As done previously in 2020 for chicken assemblies galgal5 and GRCg6a, we provide a new gene atlas, lncRNA-enriched, for the latest GRCg7b chicken assembly, integrating "NCBI RefSeq", "EMBL-EBI Ensembl/GENCODE" reference annotations and other resources such as FAANG and NONCODE. As a result, the number of PCGs increases from 18,022 (RefSeq) and 17,007 (Ensembl) to 24,102, and that of lncRNAs from 5789 (RefSeq) and 11,944 (Ensembl) to 44,428. Using 1400 public RNA-seq transcriptome representing 47 tissues, we provided expression evidence for 35,257 (79%) lncRNAs and 22,468 (93%) PCGs, supporting the relevance of this atlas. Further characterization including tissue-specificity, sex-differential expression and gene configurations are provided. We also identified conserved miRNA-hosting genes with human counterparts, suggesting common function. The annotated atlas is available at gega.sigenae.org.

The long non-coding RNA LINDA restrains cellular collapse following DNA damage in Arabidopsis thaliana

The genomic integrity of every organism is endangered by various intrinsic and extrinsic stresses. To maintain genomic integrity, a sophisticated DNA damage response (DDR) network is activated rapidly after DNA damage. Notably, the fundamental DDR mechanisms are conserved in eukaryotes. However, knowledge about many regulatory aspects of the plant DDR is still limited. Important, yet little understood, regulatory factors of the DDR are the long non-coding RNAs (lncRNAs). In humans, 13 lncRNAs functioning in DDR have been characterized to date, whereas no such lncRNAs have been characterized in plants yet. By meta-analysis, we identified the putative long intergenic non-coding RNA induced by DNA damage (LINDA) that responds strongly to various DNA double-strand break-inducing treatments, but not to replication stress induced by mitomycin C. After DNA damage, LINDA is rapidly induced in an ATM- and SOG1-dependent manner. Intriguingly, the transcriptional response of LINDA to DNA damage is similar to that of its flanking hypothetical protein-encoding gene. Phylogenetic analysis of putative Brassicales and Malvales LINDA homologs indicates that LINDA lncRNAs originate from duplication of a flanking small protein-encoding gene followed by pseudogenization. We demonstrate that LINDA is not only needed for the regulation of this flanking gene but also fine-tuning of the DDR after the occurrence of DNA double-strand breaks. Moreover, Δlinda mutant root stem cells are unable to recover from DNA damage, most likely due to hyper-induced cell death.

Multidimensional conservation analysis decodes the expression of conserved long noncoding RNAs

Although long noncoding RNAs (lncRNAs) experience weaker evolutionary constraints and exhibit lower sequence conservation than coding genes, they can still conserve their features in various aspects. Here, we used multiple approaches to systemically evaluate the conservation between human and mouse lncRNAs from various dimensions including sequences, promoter, global synteny, and local synteny, which led to the identification of 1,731 conserved lncRNAs with 427 high-confidence ones meeting multiple criteria. Conserved lncRNAs, compared with non-conserved ones, generally have longer gene bodies, more exons and transcripts, stronger connections with human diseases, and are more abundant and widespread across different tissues. Transcription factor (TF) profile analysis revealed a significant enrichment of TF types and numbers in the promoters of conserved lncRNAs. We further identified a set of TFs that preferentially bind to conserved lncRNAs and exert stronger regulation on conserved than non-conserved lncRNAs. Our study has reconciled some discrepant interpretations of lncRNA conservation and revealed a new set of transcriptional factors ruling the expression of conserved lncRNAs.

Species-specific regulation of XIST by the JPX/FTX orthologs

X chromosome inactivation (XCI) is an essential process, yet it initiates with remarkable diversity in various mammalian species. XIST, the main trigger of XCI, is controlled in the mouse by an interplay of lncRNA genes (LRGs), some of which evolved concomitantly to XIST and have orthologues across all placental mammals. Here, we addressed the functional conservation of human orthologues of two such LRGs, FTX and JPX. By combining analysis of single-cell RNA-seq data from early human embryogenesis with various functional assays in matched human and mouse pluripotent stem- or differentiated post-XCI cells, we demonstrate major functional differences for these orthologues between species, independently of primary sequence conservation. While the function of FTX is not conserved in humans, JPX stands as a major regulator of XIST expression in both species. However, we show that different entities of JPX control the production of XIST at various steps depending on the species. Altogether, our study highlights the functional versatility of LRGs across evolution, and reveals that functional conservation of orthologous LRGs may involve diversified mechanisms of action. These findings represent a striking example of how the evolvability of LRGs can provide adaptative flexibility to constrained gene regulatory networks.

De novo emergence, existence, and demise of a protein-coding gene in murids

BACKGROUND

Genes, principal units of genetic information, vary in complexity and evolutionary history. Less-complex genes (e.g., long non-coding RNA (lncRNA) expressing genes) readily emerge de novo from non-genic sequences and have high evolutionary turnover. Genesis of a gene may be facilitated by adoption of functional genic sequences from retrotransposon insertions. However, protein-coding sequences in extant genomes rarely lack any connection to an ancestral protein-coding sequence.

RESULTS

We describe remarkable evolution of the murine gene D6Ertd527e and its orthologs in the rodent Muroidea superfamily. The D6Ertd527e emerged in a common ancestor of mice and hamsters most likely as a lncRNA-expressing gene. A major contributing factor was a long terminal repeat (LTR) retrotransposon insertion carrying an oocyte-specific promoter and a 5' terminal exon of the gene. The gene survived as an oocyte-specific lncRNA in several extant rodents while in some others the gene or its expression were lost. In the ancestral lineage of Mus musculus, the gene acquired protein-coding capacity where the bulk of the coding sequence formed through CAG (AGC) trinucleotide repeat expansion and duplications. These events generated a cytoplasmic serine-rich maternal protein. Knock-out of D6Ertd527e in mice has a small but detectable effect on fertility and the maternal transcriptome.

CONCLUSIONS

While this evolving gene is not showing a clear function in laboratory mice, its documented evolutionary history in Muroidea during the last ~ 40 million years provides a textbook example of how a several common mutation events can support de novo gene formation, evolution of protein-coding capacity, as well as gene's demise.

Prostate cancer resistance leads to a global deregulation of translation factors and unconventional translation

Emerging evidence associates translation factors and regulators to tumorigenesis. However, our understanding of translational changes in cancer resistance is still limited. Here, we generated an enzalutamide-resistant prostate cancer (PCa) model, which recapitulated key features of clinical enzalutamide-resistant PCa. Using this model and poly(ribo)some profiling, we investigated global translation changes that occur during acquisition of PCa resistance. We found that enzalutamide-resistant cells exhibit an overall decrease in mRNA translation with a specific deregulation in the abundance of proteins involved in mitochondrial processes and in translational regulation. However, several mRNAs escape this translational downregulation and are nonetheless bound to heavy polysomes in enzalutamide-resistant cells suggesting active translation. Moreover, expressing these corresponding genes in enzalutamide-sensitive cells promotes resistance to enzalutamide treatment. We also found increased association of long non-coding RNAs (lncRNAs) with heavy polysomes in enzalutamide-resistant cells, suggesting that some lncRNAs are actively translated during enzalutamide resistance. Consistent with these findings, expressing the predicted coding sequences of known lncRNAs JPX, CRNDE and LINC00467 in enzalutamide-sensitive cells drove resistance to enzalutamide. Taken together, this suggests that aberrant translation of specific mRNAs and lncRNAs is a strong indicator of PCa enzalutamide resistance, which points towards novel therapeutic avenues that may target enzalutamide-resistant PCa.

The Genomics of Auditory Function and Disease

Current estimates suggest that nearly half a billion people worldwide are affected by hearing loss. Because of the major psychological, social, economic, and health ramifications, considerable efforts have been invested in identifying the genes and molecular pathways involved in hearing loss, whether genetic or environmental, to promote prevention, improve rehabilitation, and develop therapeutics. Genomic sequencing technologies have led to the discovery of genes associated with hearing loss. Studies of the transcriptome and epigenome of the inner ear have characterized key regulators and pathways involved in the development of the inner ear and have paved the way for their use in regenerative medicine. In parallel, the immense preclinical success of using viral vectors for gene delivery in animal models of hearing loss has motivated the industry to work on translating such approaches into the clinic. Here, we review the recent advances in the genomics of auditory function and dysfunction, from patient diagnostics to epigenetics and gene therapy.

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.

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.

High-Throughput Evolutionary Comparative Analysis of Long Intergenic Noncoding RNAs in Multiple Organisms

Comparative genomic and transcriptomic analyses can help prioritize and facilitate the functional analysis of long noncoding RNAs (lncRNAs). Evolinc-II is a bioinformatic pipeline that automates comparative analyses, searching for sequence and structural conservation for thousands of lncRNAs at once. In addition, Evolinc-II takes a phylogenetic approach to infer key evolutionary events that may have occurred during the emergence of each query lncRNA. Here, we describe how to use command line or GUI (CyVerse's Discovery Environment) versions of Evolinc-II to identify lncRNA homologs and prioritize them for functional analysis.

The World of Pseudogenes: New Diagnostic and Therapeutic Targets in Cancers or Still Mystery Molecules?

Pseudogenes were once considered as “junk DNA”, due to loss of their functions as a result of the accumulation of mutations, such as frameshift and presence of premature stop-codons and relocation of genes to inactive heterochromatin regions of the genome. Pseudogenes are divided into two large groups, processed and unprocessed, according to their primary structure and origin. Only 10% of all pseudogenes are transcribed into RNAs and participate in the regulation of parental gene expression at both transcriptional and translational levels through senseRNA (sRNA) and antisense RNA (asRNA). In this review, about 150 pseudogenes in the different types of cancers were analyzed. Part of these pseudogenes seem to be useful in molecular diagnostics and can be detected in various types of biological material including tissue as well as biological fluids (liquid biopsy) using different detection methods. The number of pseudogenes, as well as their function in the human genome, is still unknown. However, thanks to the development of various technologies and bioinformatic tools, it was revealed so far that pseudogenes are involved in the development and progression of certain diseases, especially in cancer.

LncRNAs in domesticated animals: from dog to livestock species

Animal genomes are pervasively transcribed into multiple RNA molecules, of which many will not be translated into proteins. One major component of this transcribed non-coding genome is the long non-coding RNAs (lncRNAs), which are defined as transcripts longer than 200 nucleotides with low coding-potential capabilities. Domestic animals constitute a unique resource for studying the genetic and epigenetic basis of phenotypic variations involving protein-coding and non-coding RNAs, such as lncRNAs. This review presents the current knowledge regarding transcriptome-based catalogues of lncRNAs in major domesticated animals (pets and livestock species), covering a broad phylogenetic scale (from dogs to chicken), and in comparison with human and mouse lncRNA catalogues. Furthermore, we describe different methods to extract known or discover novel lncRNAs and explore comparative genomics approaches to strengthen the annotation of lncRNAs. We then detail different strategies contributing to a better understanding of lncRNA functions, from genetic studies such as GWAS to molecular biology experiments and give some case examples in domestic animals. Finally, we discuss the limitations of current lncRNA annotations and suggest research directions to improve them and their functional characterisation.

HAR1: an insight into lncRNA genetic evolution

Long noncoding RNAs (lncRNAs) have a wide range of functions in health and disease, but many remain uncharacterized because of their complex expression patterns and structures. The genetic loci encoding lncRNAs can be subject to accelerated evolutionary changes within the human lineage. HAR1 is a region that has a significantly altered sequence compared to other primates and is a component of two overlapping lncRNA loci, HAR1A and HAR1B. Although the functions of these lncRNAs are unknown, they have been associated with neurological disorders and cancer. Here, we explore the current state of understanding of evolution in human lncRNA genes, using the HAR1 locus as the case study.

A vast pool of lineage-specific microproteins encoded by long non-coding RNAs in plants

Pervasive transcription of eukaryotic genomes results in expression of long non-coding RNAs (lncRNAs) most of which are poorly conserved in evolution and appear to be non-functional. However, some lncRNAs have been shown to perform specific functions, in particular, transcription regulation. Thousands of small open reading frames (smORFs, <100 codons) located on lncRNAs potentially might be translated into peptides or microproteins. We report a comprehensive analysis of the conservation and evolutionary trajectories of lncRNAs-smORFs from the moss Physcomitrium patens across transcriptomes of 479 plant species. Although thousands of smORFs are subject to substantial purifying selection, the majority of the smORFs appear to be evolutionary young and could represent a major pool for functional innovation. Using nanopore RNA sequencing, we show that, on average, the transcriptional level of conserved smORFs is higher than that of non-conserved smORFs. Proteomic analysis confirmed translation of 82 novel species-specific smORFs. Numerous conserved smORFs containing low complexity regions (LCRs) or transmembrane domains were identified, the biological functions of a selected LCR-smORF were demonstrated experimentally. Thus, microproteins encoded by smORFs are a major, functionally diverse component of the plant proteome.

Long noncoding RNAs as tumorigenic factors and therapeutic targets for renal cell carcinoma

Approximately 338,000 patients are diagnosed with kidney cancer worldwide each year, and renal cell carcinoma (RCC), which is derived from renal epithelium, accounts for more than ninety percent of the malignancy. Next generation RNA sequencing has enabled the identification of novel long noncoding RNAs (lncRNAs) in the past 10 years. Recent studies have provided extensive evidence that lncRNAs bind to chromatin modification proteins, transcription factors, RNA-binding proteins and microRNAs, and thereby modulate gene expression through regulating chromatin status, gene transcription, pre-mRNA splicing, mRNA decay and stability, protein translation and stability. In vitro and in vivo studies have demonstrated that over-expression of oncogenic lncRNAs and silencing of tumor suppressive lncRNAs are a common feature of human RCC, and that aberrant lncRNA expression is a marker for poor patient prognosis, and is essential for the initiation and progression of RCC. Because lncRNAs, compared with mRNAs, are expressed in a tissue-specific manner, aberrantly expressed lncRNAs can be better targeted for the treatment of RCC through screening small molecule compounds which block the interaction between lncRNAs and their binding proteins or microRNAs.

Long non-coding RNAs in development and disease: conservation to mechanisms

Our genomes contain the blueprint of what makes us human and many indications as to why we develop disease. Until the last 10 years, most studies had focussed on protein-coding genes, more specifically DNA sequences coding for proteins. However, this represents less than 5% of our genomes. The other 95% is referred to as the 'dark matter' of our genomes, our understanding of which is extremely limited. Part of this 'dark matter' includes regions that give rise to RNAs that do not code for proteins. A subset of these non-coding RNAs are long non-coding RNAs (lncRNAs), which in particular are beginning to be dissected and their importance to human health revealed. To improve our understanding and treatment of disease it is vital that we understand the molecular and cellular function of lncRNAs, and how their misregulation can contribute to disease. It is not yet clear what proportion of lncRNAs is actually functional; conservation during evolution is being used to understand the biological importance of lncRNA. Here, we present key themes within the field of lncRNAs, emphasising the importance of their roles in both the nucleus and the cytoplasm of cells, as well as patterns in their modes of action. We discuss their potential functions in development and disease using examples where we have the greatest understanding. Finally, we emphasise why lncRNAs can serve as biomarkers and discuss their emerging potential for therapy. © 2020 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

Formation of human long intergenic non-coding RNA genes, pseudogenes, and protein genes: Ancestral sequences are key players

Pathways leading to formation of non-coding RNA and protein genes are varied and complex. We report finding a conserved repeat sequence present in human and chimpanzee genomes that appears to have originated from a common primate ancestor. This sequence is repeatedly copied in human chromosome 22 (chr22) low copy repeats (LCR22) or segmental duplications and forms twenty-one different genes, which include the human long intergenic non-coding RNA (lincRNA) family FAM230, a newly discovered lincRNA gene family termed conserved long intergenic non-coding RNAs (clincRNA), pseudogene families, as well as the gamma-glutamyltransferase (GGT) protein gene family and the RNA pseudogenes that originate from GGT sequences. Of particular interest are the GGT5 and USP18 protein genes that appear to have formed from an homologous repeat sequence that also forms the clincRNA gene family. The data point to ancestral DNA sequences, conserved through evolution and duplicated in humans by chromosomal repeat sequences that may serve as functional genomic elements in the development of diverse genes.

Evolution of new proteins from translated sORFs in long non-coding RNAs

High throughput RNA sequencing techniques have revealed that a large fraction of the genome is transcribed into long non-coding RNAs (lncRNAs). Unlike canonical protein-coding genes, lncRNAs do not contain long open reading frames (ORFs) and tend to be poorly conserved across species. However, many of them contain small ORFs (sORFs) that exhibit translation signatures according to ribosome profiling or proteomics data. These sORFs are a source of putative novel proteins; some of them may confer a selective advantage and be maintained over time, a process known as de novo gene birth. Here we review the mechanisms by which randomly occurring sORFs in lncRNAs can become new functional proteins.

Functional Conservation of LncRNA JPX Despite Sequence and Structural Divergence

Long noncoding RNAs (lncRNAs) have been identified in all eukaryotes and are most abundant in the human genome. However, the functional importance and mechanisms of action for human lncRNAs are largely unknown. Using comparative sequence, structural, and functional analyses, we characterize the evolution and molecular function of human lncRNA JPX. We find that human JPX and its mouse homolog, lncRNA Jpx, have deep divergence in their nucleotide sequences and RNA secondary structures. Despite such differences, both lncRNAs demonstrate robust binding to CTCF, a protein that is central to Jpx's role in X chromosome inactivation. In addition, our functional rescue experiment using Jpx-deletion mutant cells shows that human JPX can functionally complement the loss of Jpx in mouse embryonic stem cells. Our findings support a model for functional conservation of lncRNAs independent from sequence and structural divergence. This study provides mechanistic insight into the evolution of lncRNA function.

Identification of transcribed protein coding sequence remnants within lincRNAs

Long intergenic non-coding RNAs (lincRNAs) are non-coding transcripts >200 nucleotides long that do not overlap protein-coding sequences. Importantly, such elements are known to be tissue-specifically expressed and to play a widespread role in gene regulation across thousands of genomic loci. However, very little is known of the mechanisms for the evolutionary biogenesis of these RNA elements, especially given their poor conservation across species. It has been proposed that lincRNAs might arise from pseudogenes. To test this systematically, we developed a novel method that searches for remnants of protein-coding sequences within lincRNA transcripts; the hypothesis is that we can trace back their biogenesis from protein-coding genes or posterior transposon/retrotransposon insertions. Applying this method, we found 203 human lincRNA genes with regions significantly similar to protein-coding sequences. Our method provides a visualization tool to trace the evolutionary biogenesis of lincRNAs with respect to protein-coding genes by sequence divergence. Subsequently, we show the expression correlation between lincRNAs and their identified parental protein-coding genes using public RNA-seq repositories, hinting at novel gene regulatory relationships. In summary, we developed a novel computational methodology to study non-coding gene sequences, which can be applied to identify the evolutionary biogenesis and function of lincRNAs.

Developmental dynamics of lncRNAs across mammalian organs and species

While many long noncoding RNAs (lncRNAs) have been identified in human and other mammalian genomes, there has been limited systematic functional characterization. In particular, the contribution of lncRNAs to organ development remains largely unexplored. Here we analyze the expression patterns of lncRNAs across developmental timepoints in seven major organs, from early organogenesis to adulthood, across seven species (human, macaque, mouse, rat, rabbit, opossum, and chicken). Our analyses identified ~15,000-35,000 candidate lncRNAs in each species, most of which show species specificity. We characterized expression patterns of lncRNAs across developmental stages, and found many with dynamic expression patterns across time that show signatures of enrichment for functionality. During development, there is a transition from broadly expressed and conserved lncRNAs towards an increasing number of lineage- and organ-specific lncRNAs. Our study provides a resource of candidate lncRNAs and their patterns of expression and evolutionary conservation across mammalian organ development.

Computational Methods for Mapping, Assembly and Quantification for Coding and Non-coding Transcripts

The measurement of gene expression has long provided significant insight into biological functions. The development of high-throughput short-read sequencing technology has revealed transcriptional complexity at an unprecedented scale, and informed almost all areas of biology. However, as researchers have sought to gather more insights from the data, these new technologies have also increased the computational analysis burden. In this review, we describe typical computational pipelines for RNA-Seq analysis and discuss their strengths and weaknesses for the assembly, quantification and analysis of coding and non-coding RNAs. We also discuss the assembly of transposable elements into transcripts, and the difficulty these repetitive elements pose. In summary, RNA-Seq is a powerful technology that is likely to remain a key asset in the biologist's toolkit.

Improved annotation of the domestic pig genome through integration of Iso-Seq and RNA-seq data

BACKGROUND

Our understanding of the pig transcriptome is limited. RNA transcript diversity among nine tissues was assessed using poly(A) selected single-molecule long-read isoform sequencing (Iso-seq) and Illumina RNA sequencing (RNA-seq) from a single White cross-bred pig.

RESULTS

Across tissues, a total of 67,746 unique transcripts were observed, including 60.5% predicted protein-coding, 36.2% long non-coding RNA and 3.3% nonsense-mediated decay transcripts. On average, 90% of the splice junctions were supported by RNA-seq within tissue. A large proportion (80%) represented novel transcripts, mostly produced by known protein-coding genes (70%), while 17% corresponded to novel genes. On average, four transcripts per known gene (tpg) were identified; an increase over current EBI (1.9 tpg) and NCBI (2.9 tpg) annotations and closer to the number reported in human genome (4.2 tpg). Our new pig genome annotation extended more than 6000 known gene borders (5' end extension, 3' end extension, or both) compared to EBI or NCBI annotations. We validated a large proportion of these extensions by independent pig poly(A) selected 3'-RNA-seq data, or human FANTOM5 Cap Analysis of Gene Expression data. Further, we detected 10,465 novel genes (81% non-coding) not reported in current pig genome annotations. More than 80% of these novel genes had transcripts detected in > 1 tissue. In addition, more than 80% of novel intergenic genes with at least one transcript detected in liver tissue had H3K4me3 or H3K36me3 peaks mapping to their promoter and gene body, respectively, in independent liver chromatin immunoprecipitation data.

CONCLUSIONS

These validated results show significant improvement over current pig genome annotations.

The long noncoding RNA ROCKI regulates inflammatory gene expression

Long noncoding RNAs (lncRNAs) can regulate target gene expression by acting in cis (locally) or in trans (non-locally). Here, we performed genome-wide expression analysis of Toll-like receptor (TLR)-stimulated human macrophages to identify pairs of cis-acting lncRNAs and protein-coding genes involved in innate immunity. A total of 229 gene pairs were identified, many of which were commonly regulated by signaling through multiple TLRs and were involved in the cytokine responses to infection by group B Streptococcus We focused on elucidating the function of one lncRNA, named lnc-MARCKS or ROCKI (Regulator of Cytokines and Inflammation), which was induced by multiple TLR stimuli and acted as a master regulator of inflammatory responses. ROCKI interacted with APEX1 (apurinic/apyrimidinic endodeoxyribonuclease 1) to form a ribonucleoprotein complex at the MARCKS promoter. In turn, ROCKI-APEX1 recruited the histone deacetylase HDAC1, which removed the H3K27ac modification from the promoter, thus reducing MARCKS transcription and subsequent Ca2+ signaling and inflammatory gene expression. Finally, genetic variants affecting ROCKI expression were linked to a reduced risk of certain inflammatory and infectious disease in humans, including inflammatory bowel disease and tuberculosis. Collectively, these data highlight the importance of cis-acting lncRNAs in TLR signaling, innate immunity, and pathophysiological inflammation.

Long non-coding RNA in stem cell pluripotency and lineage commitment: functions and evolutionary conservation

LncRNAs have recently emerged as new and fundamental transcriptional and post-transcriptional regulators acting at multiple levels of gene expression. Indeed, lncRNAs participate in a wide variety of stem cell and developmental processes, acting in cis and/or in trans in the nuclear and/or in the cytoplasmic compartments, and generating an intricate network of interactions with RNAs, enhancers, and chromatin-modifier complexes. Given the versatility of these molecules to operate in different subcellular compartments, via different modes of action and with different target specificity, the interest in this research field is rapidly growing. Here, we review recent progress in defining the functional role of lncRNAs in stem cell biology with a specific focus on the underlying mechanisms. We also discuss recent findings on a new family of evolutionary conserved lncRNAs transcribed from ultraconserved elements, which show perfect conservation between human, mouse, and rat genomes, and that are emerging as new player in this complex scenario.

Ancient exapted transposable elements promote nuclear enrichment of human long noncoding RNAs

The sequence domains underlying long noncoding RNA (lncRNA) activities, including their characteristic nuclear enrichment, remain largely unknown. It has been proposed that these domains can originate from neofunctionalized fragments of transposable elements (TEs), otherwise known as RIDLs (repeat insertion domains of lncRNA), although just a handful have been identified. It is challenging to distinguish functional RIDL instances against a numerous genomic background of neutrally evolving TEs. We here show evidence that a subset of TE types experience evolutionary selection in the context of lncRNA exons. Together these comprise an enrichment group of 5374 TE fragments in 3566 loci. Their host lncRNAs tend to be functionally validated and associated with disease. This RIDL group was used to explore the relationship between TEs and lncRNA subcellular localization. By using global localization data from 10 human cell lines, we uncover a dose-dependent relationship between nuclear/cytoplasmic distribution and evolutionarily conserved L2b, MIRb, and MIRc elements. This is observed in multiple cell types and is unaffected by confounders of transcript length or expression. Experimental validation with engineered transgenes shows that these TEs drive nuclear enrichment in a natural sequence context. Together these data reveal a role for TEs in regulating the subcellular localization of lncRNAs.

Evolutionary Patterns of Non-Coding RNA in Cardiovascular Biology

Cardiovascular diseases (CVDs) affect the heart and the vascular system with a high prevalence and place a huge burden on society as well as the healthcare system. These complex diseases are often the result of multiple genetic and environmental risk factors and pose a great challenge to understanding their etiology and consequences. With the advent of next generation sequencing, many non-coding RNA transcripts, especially long non-coding RNAs (lncRNAs), have been linked to the pathogenesis of CVD. Despite increasing evidence, the proper functional characterization of most of these molecules is still lacking. The exploration of conservation of sequences across related species has been used to functionally annotate protein coding genes. In contrast, the rapid evolutionary turnover and weak sequence conservation of lncRNAs make it difficult to characterize functional homologs for these sequences. Recent studies have tried to explore other dimensions of interspecies conservation to elucidate the functional role of these novel transcripts. In this review, we summarize various methodologies adopted to explore the evolutionary conservation of cardiovascular non-coding RNAs at sequence, secondary structure, syntenic, and expression level.

Oncogenicity of lncRNA FOXD2-AS1 and its molecular mechanisms in human cancers

OBJECTIVES

Long non-coding RNAs (lncRNAs) are a group of noncoding RNAs with length larger than 200 nucleotides. LncRNAs have limited or no protein-coding capacity because of lack of obvious open reading frame. An increasing number of researches have shown that lncRNAs participate in the complex regulation network of cancer and play an important role in tumourigenesis and progression such as proliferation, migration and invasion. LncRNA FOXD2 adjacent opposite strand RNA 1 (FOXD2-AS1), located on chromosome 1p33 and with a transcript length of 2527 nucleotides, is a novel cancer-related lncRNA. FOXD2-AS1 was recently found to exhibit aberrant expression in various malignancies, including gastric, lung, bladder, colorectal, nasopharyngeal, esophageal, hepatocellular, thyroid and skin cancer, and its deregulation might be related to survival and prognosis of cancer patients. Pertinent to clinical practice, FOXD2-AS1 might act as a feasible biomarker or therapeutic target in human cancers. In this paper, we made a summary on the current findings concerning the biological functions and molecular mechanisms of FOXD2-AS1 in tumor progression.

MATERIALS AND METHODS

In this paper, we summarized and figured out recent studies about the expression and molecular biological mechanisms of FOXD2-AS1 in tumor progression. Existing relevant studies were obtained through a systematic search from PubMed, Embase, BioMedNet, GEO database and Cochrane Library.

RESULTS

FOXD2-AS1 was a valuable tumor-associated lncRNA. Its expression level was up-regulation in various malignancies, including gastric, lung, bladder, colorectal, nasopharyngeal, esophageal, hepatocellular, thyroid and skin cancer. In addition, the aberrant expressions of FOXD2-AS1 have shown to contribute to proliferation, migration and invasion of cancer cells, and its deregulation is related to carcinogensis, overall survival, disease free survival, prognosis and tumor progression.

CONCLUSIONS

LncRNA FOXD2-AS1 is an oncogene and probably represents a feasible biomarker or therapeutic target in human cancers.

Translation of Small Open Reading Frames: Roles in Regulation and Evolutionary Innovation

The translatome can be defined as the sum of the RNA sequences that are translated into proteins in the cell by the ribosomal machinery. Until recently, it was generally assumed that the translatome was essentially restricted to evolutionary conserved proteins encoded by the set of annotated protein-coding genes. However, it has become increasingly clear that it also includes small regulatory open reading frames (ORFs), functional micropeptides, de novo proteins, and the pervasive translation of likely nonfunctional proteins. Many of these ORFs have been discovered thanks to the development of ribosome profiling, a technique to sequence ribosome-protected RNA fragments. To fully capture the diversity of translated ORFs, we propose a comprehensive classification that includes the new types of translated ORFs in addition to standard proteins.

Ancient exapted transposable elements promote nuclear enrichment of human long noncoding RNAs

The sequence domains underlying long noncoding RNA (lncRNA) activities, including their characteristic nuclear enrichment, remain largely unknown. It has been proposed that these domains can originate from neofunctionalized fragments of transposable elements (TEs), otherwise known as RIDLs (repeat insertion domains of lncRNA), although just a handful have been identified. It is challenging to distinguish functional RIDL instances against a numerous genomic background of neutrally evolving TEs. We here show evidence that a subset of TE types experience evolutionary selection in the context of lncRNA exons. Together these comprise an enrichment group of 5374 TE fragments in 3566 loci. Their host lncRNAs tend to be functionally validated and associated with disease. This RIDL group was used to explore the relationship between TEs and lncRNA subcellular localization. By using global localization data from 10 human cell lines, we uncover a dose-dependent relationship between nuclear/cytoplasmic distribution and evolutionarily conserved L2b, MIRb, and MIRc elements. This is observed in multiple cell types and is unaffected by confounders of transcript length or expression. Experimental validation with engineered transgenes shows that these TEs drive nuclear enrichment in a natural sequence context. Together these data reveal a role for TEs in regulating the subcellular localization of lncRNAs.

Exaptation at the molecular genetic level

The realization that body parts of animals and plants can be recruited or coopted for novel functions dates back to, or even predates the observations of Darwin. S.J. Gould and E.S. Vrba recognized a mode of evolution of characters that differs from adaptation. The umbrella term aptation was supplemented with the concept of exaptation. Unlike adaptations, which are restricted to features built by selection for their current role, exaptations are features that currently enhance fitness, even though their present role was not a result of natural selection. Exaptations can also arise from nonaptations; these are characters which had previously been evolving neutrally. All nonaptations are potential exaptations. The concept of exaptation was expanded to the molecular genetic level which aided greatly in understanding the enormous potential of neutrally evolving repetitive DNA-including transposed elements, formerly considered junk DNA-for the evolution of genes and genomes. The distinction between adaptations and exaptations is outlined in this review and examples are given. Also elaborated on is the fact that such distinctions are sometimes more difficult to determine; this is a widespread phenomenon in biology, where continua abound and clear borders between states and definitions are rare.

Pseudogenes as Functionally Significant Elements of the Genome

Pseudogene is a gene copy that has lost its original function. For a long time, pseudogenes have been considered as "junk DNA" that inevitably arises as a result of ongoing evolutionary process. However, experimental data obtained during recent years indicate this understanding of the nature of pseudogenes is not entirely correct, and many pseudogenes perform important genetic functions. In the review, we have addressed classification of pseudogenes, methods of their detection in the genome, and the problem of their evolutionary conservatism and prevalence among species belonging to different taxonomic groups in the light of modern data. The mechanisms of gene expression regulation by pseudogenes and the role of pseudogenes in pathogenesis of various human diseases are discussed.

Production of Spliced Long Noncoding RNAs Specifies Regions with Increased Enhancer Activity

Active enhancers in mammals produce enhancer RNAs (eRNAs) that are bidirectionally transcribed, unspliced, and unstable. Enhancer regions are also enriched with long noncoding RNA (lncRNA) transcripts, which are typically spliced and substantially more stable. In order to explore the relationship between these two classes of RNAs, we analyzed DNase hypersensitive sites with evidence of bidirectional transcription, which we termed eRNA-producing centers (EPCs). EPCs found very close to transcription start sites of lncRNAs exhibit attributes of both enhancers and promoters, including distinctive DNA motifs and a characteristic chromatin landscape. These EPCs are associated with higher enhancer activity, driven at least in part by the presence of conserved, directional splicing signals that promote lncRNA production, pointing at a causal role of lncRNA processing in enhancer activity. Together, our results suggest that the conserved ability of some enhancers to produce lncRNAs augments their activity in a manner likely mediated through lncRNA maturation.

Copaifera langsdorffii Novel Putative Long Non-Coding RNAs: Interspecies Conservation Analysis in Adaptive Response to Different Biomes

Long non-coding RNAs (lncRNAs) are involved in multiple regulatory pathways and its versatile form of action has disclosed a new layer in gene regulation. LncRNAs have their expression levels modulated during plant development, and in response to stresses with tissue-specific functions. In this study, we analyzed lncRNA from leaf samples collected from the legume Copaifera langsdorffii Desf. (copaíba) present in two divergent ecosystems: Cerrado (CER; Ecological Station of Botanical Garden in Brasília, Brazil) and Atlantic Rain Forest (ARF; Rio de Janeiro, Brazil). We identified 8020 novel lncRNAs, and they were compared to seven Fabaceae genomes and transcriptomes, to which 1747 and 2194 copaíba lncRNAs were mapped, respectively, to at least one species. The secondary structures of the lncRNAs that were conserved and differentially expressed between the populations were predicted using in silico methods. A few selected lncRNA were confirmed by RT-qPCR in the samples from both biomes; Additionally, the analysis of the lncRNA sequences predicted that some might act as microRNA (miRNA) targets or decoys. The emerging studies involving lncRNAs function and conservation have shown their involvement in several types of biotic and abiotic stresses. Thus, the conservation of lncRNAs among Fabaceae species considering their rapid turnover, suggests they are likely to have been under functional conservation pressure. Our results indicate the potential involvement of lncRNAs in the adaptation of C. langsdorffii in two different biomes.

The role of the long non-coding RNA HOXA11-AS in promoting proliferation and metastasis of malignant tumors

Long non-coding RNA (lncRNA) is one of the focuses and hotspots of biological research in recent years. At the same time, tumors have become the main disease that endangers human health. In recent years, a large number of researchers have explored the relationship between lncRNA and tumors. HOXA11-AS is one of these lncRNAs. The long non-coding RNA HOXA11 antisense RNA (HOXA11-AS) is a novel lncRNA recently discovered. It is found in a variety of tumors such as ovarian cancer, glioma, and gastric cancer (GC) and so on, and is defined as an oncogene. It promotes tumor proliferation and metastasis by interacting with proteins such as miRNA and EZH2. In this paper, we review the mechanism of interaction between HOXA11-AS and various tumors in recent years, and believe that it can be a potential tumor marker and therapeutic target in the future prevention and treatment of tumors.

Long noncoding RNAs are associated with metabolic and cellular processes in the jejunum mucosa of pre-weaning calves in response to different diets

Long noncoding RNAs (lncRNAs) emerged as important regulatory component of mechanisms involved in gene expression, chromatin modification and epigenetic processes, but they are rarely annotated in the bovine genome.

Our study monitored the jejunum transcriptome of German Holstein calves fed two different milk diets using transcriptome sequencing (RNA-seq). To identify potential lncRNAs within the pool of unknown transcripts, four bioinformatic lncRNA prediction tools were applied. The intersection of the alignment-free lncRNA prediction tools (CNCI, PLEK and FEELnc) predicted 1,812 lncRNA transcripts concordantly comprising a catalogue of 1,042 putative lncRNA loci expressed in the calves’ intestinal mucosa.

Nine lncRNA loci were differentially expressed (DE lncRNAs) between both calf groups. To elucidate their biological function, we applied a systems biology approach that combines weighted gene co-expression network analysis with functional enrichment and biological pathway analysis. Four DE lncRNAs were found to be strongly correlated with a gene network module (GNM) enriched for genes from canonical pathways of remodeling of epithelial adherens junction, tight junction and integrin signaling. Another DE lncRNA was strongly correlated with a GNM enriched for genes associated with energy metabolism and maintaining of cellular homeostasis with a focus on mitochondrial processes.

Our data suggest that these DE lncRNAs may play potential regulatory roles in modulating biological processes associated with energy metabolism pathways and cellular signaling processes affecting the barrier function of intestinal epithelial cells of calves in response to different feeding regimens in the pre-weaning period.

Paraquat and MPTP induce alteration in the expression profile of long noncoding RNAs in the substantia nigra of mice: Role of the transcription factor Nrf2

Parkinson's disease (PD) is a common age-related degenerative disease of the central nervous system caused mainly by hereditary, pesticides, metals, and polychlorinated biphenyls. Paraquat (PQ), a widely used herbicide, causes PD. Long noncoding RNAs (lncRNAs) are nonprotein-coding transcripts, expressed in the brain and play irreplaceable roles in neurodegenerative diseases. NF-E2-related factor-2 (Nrf2) is an important genetic transcription regulator in oxidative stress. We aimed to discover novel PQ or 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP)-Nrf2-related lncRNAs and explore their association with PD. 17157 lncRNAs and 13707 mRNAs (fold change ≥2, P < 0.05) were identified by Microarray. And the expressions of six lncRNAs were confirmed by using qRT-PCR and two by FISH. Coding-noncoding analysis and qRT-PCR were applied to discover the functions of lncRNAs and predict the targeted genes. In mice, PQ and MPTP exposure caused alteration of the lncRNA expression profile, suggesting lncRNAs may be involved in PQ- and MPTP-induced neurotoxicity. The changes in their lncRNA expression were distinct but related. PQ caused lncRNA expression profiling alteration in the substantia nigra (SN) through an interaction with Nrf2, thus changing the NR_027648/Zc3h14/Cybb and NR_030777/Zfp326/Cpne5 mRNA pathways. Similarly, MPTP caused lncRNA expression profiling alteration in SN through an interaction with Nrf2. Nrf2 may be involved in the development of neurodegeneration induced by PQ and MPTP via interaction with lncRNAs as the molecular mechanism. Our findings indicate the potential roles of lncRNAs in the development of PD by PQ or MPTP and provide positive insights into future mechanism studies.

Long noncoding RNAs: Undeciphered cellular codes encrypting keys of colorectal cancer pathogenesis

Long noncoding RNAs are non-protein coding transcripts longer than 200 nucleotides in length. By the advance in genetic and bioinformatic technologies, the new genomic landscape including noncoding transcripts has been revealed. Despite their non-capacity to be translated into proteins, lncRNAs have a versatile functions through various mechanisms interacting with other cellular molecules including DNA, protein, and RNA. Recent research interest and endeavor have identified the functional role of lncRNAs in various diseases including cancer. Colorectal cancer (CRC) is not only one of the most frequent cancer but also one of the cancer types with remarkable achievements in lncRNA research. Of the numerous notable lncRNAs identified and characterized in CRC, we will focus on key lncRNAs with the high potential as CRC-specific biomarkers in this review.

ZFLNC: a comprehensive and well-annotated database for zebrafish lncRNA

There is emerging evidence showing that lncRNAs can be involved in various critical biological processes. Zebrafish is a fully developed model system being used in a variety of basic research and biomedical studies. Hence, it is an ideal model organism to study the functions and mechanisms of lncRNAs. Here, we constructed ZFLNC—a comprehensive database of zebrafish lncRNA that is dedicated to providing a zebrafish-based platform for deep exploration of zebrafish lncRNAs and their mammalian counterparts to the relevant academic communities. The main data resources of lncRNAs in this database come from the NCBI, Ensembl, NONCODE, zflncRNApedia and literature. We also obtained lncRNAs as a supplement by analysing RNA-Seq datasets from SRA database. With these IncRNAs, we further carried out expression profiling, co-expression network prediction, Gene Ontology (GO)/Kyoto Encyclopedia of Genes and Genomes (KEGG)/Online Mendelian Inheritance in Man (OMIM) annotation and conservation analysis. As far as we know, ZFLNC is the most comprehensive and well-annotated database for zebrafish lncRNA.