Transposable elements are major contributors to the origin, diversification, and regulation of vertebrate long noncoding RNAs

Dynamic Signatures of the Epigenome: Friend or Foe?

Highly dynamic epigenetic signaling is influenced mainly by (micro)environmental stimuli and genetic factors. The exact mechanisms affecting particular epigenomic patterns differ dependently on the context. In the current review, we focus on the causes and effects of the dynamic signatures of the human epigenome as evaluated with the high-throughput profiling data and single-gene approaches. We will discuss three different aspects of phenotypic outcomes occurring as a consequence of epigenetics interplaying with genotype and environment. The first issue is related to the cases of environmental impacts on epigenetic profile, and its adverse and advantageous effects related to human health and evolutionary adaptation. The next topic will present a model of the interwoven co-evolution of genetic and epigenetic patterns exemplified with transposable elements (TEs) and their epigenetic repressors Krüppel-associated box zinc finger proteins (KRAB–ZNFs). The third aspect concentrates on the mitosis-based microevolution that takes place during carcinogenesis, leading to clonal diversity and expansion of tumor cells. The whole picture of epigenome plasticity and its role in distinct biological processes is still incomplete. However, accumulating data define epigenomic dynamics as an essential co-factor driving adaptation at the cellular and inter-species levels with a benefit or disadvantage to the host.

Origins and evolutionary consequences of ancient endogenous retroviruses

Retroviruses infect a broad range of vertebrate hosts that includes amphibians, reptiles, fish, birds and mammals. In addition, a typical vertebrate genome contains thousands of loci composed of ancient retroviral sequences known as endogenous retroviruses (ERVs). ERVs are molecular remnants of ancient retroviruses and proof that the ongoing relationship between retroviruses and their vertebrate hosts began hundreds of millions of years ago. The long-term impact of retroviruses on vertebrate evolution is twofold: first, as with other viruses, retroviruses act as agents of selection, driving the evolution of host genes that block viral infection or that mitigate pathogenesis, and second, through the phenomenon of endogenization, retroviruses contribute an abundance of genetic novelty to host genomes, including unique protein-coding genes and cis-acting regulatory elements. This Review describes ERV origins, their diversity and their relationships to retroviruses and discusses the potential for ERVs to reveal virus-host interactions on evolutionary timescales. It also describes some of the many examples of cellular functions, including protein-coding genes and regulatory elements, that have evolved from ERVs.

Seq'ing identity and function in a repeat-derived noncoding RNA world

Innovations in high-throughout sequencing approaches are being marshaled to both reveal the composition of the abundant and heterogeneous noncoding RNAs that populate cell nuclei and lend insight to the mechanisms by which noncoding RNAs influence chromosome biology and gene expression. This review focuses on some of the recent technological developments that have enabled the isolation of nascent transcripts and chromatin-associated and DNA-interacting RNAs. Coupled with emerging genome assembly and analytical approaches, the field is poised to achieve a comprehensive catalog of nuclear noncoding RNAs, including those derived from repetitive regions within eukaryotic genomes. Herein, particular attention is paid to the challenges and advances in the sequence analyses of repeat and transposable element-derived noncoding RNAs and in ascribing specific function(s) to such RNAs.

Transposable elements have contributed human regulatory regions that are activated upon bacterial infection

Transposable elements (TEs) are increasingly recognized as important contributors to mammalian regulatory systems. For instance, they have been shown to play a role in the human interferon response, but their involvement in other mechanisms of immune cell activation remains poorly understood. Here, we investigated the profile of accessible chromatin enhanced in stimulated human macrophages using ATAC-seq to assess the role of different TE subfamilies in regulating gene expression following an immune response. We found that both previously identified and new repeats belonging to the MER44, THE1, Tigger3 and MLT1 families provide 14 subfamilies that are enriched in differentially accessible chromatin and found near differentially expressed genes. These TEs also harbour binding motifs for several candidate transcription factors, including important immune regulators AP-1 and NF-κB, present in 96% of accessible MER44B and 83% of THE1C instances, respectively. To more directly assess their regulatory potential, we evaluated the presence of these TEs in regions putatively affecting gene expression, as defined by quantitative trait locus (QTL) analysis, and found that repeats are also contributing to accessible elements near QTLs. Together, these results suggest that a number of TE families have contributed to the regulation of gene expression in the context of the immune response to infection in humans. This article is part of a discussion meeting issue 'Crossroads between transposons and gene regulation'.

Comparative Transcriptomics Analyses across Species, Organs, and Developmental Stages Reveal Functionally Constrained lncRNAs

The functionality of long noncoding RNAs (lncRNAs) is disputed. In general, lncRNAs are under weak selective pressures, suggesting that the majority of lncRNAs may be nonfunctional. However, although some surveys showed negligible phenotypic effects upon lncRNA perturbation, key biological roles were demonstrated for individual lncRNAs. Most lncRNAs with proven functions were implicated in gene expression regulation, in pathways related to cellular pluripotency, differentiation, and organ morphogenesis, suggesting that functional lncRNAs may be more abundant in embryonic development, rather than in adult organs. To test this hypothesis, we perform a multidimensional comparative transcriptomics analysis, across five developmental time points (two embryonic stages, newborn, adult, and aged individuals), four organs (brain, kidney, liver, and testes), and three species (mouse, rat, and chicken). We find that, overwhelmingly, lncRNAs are preferentially expressed in adult and aged testes, consistent with the presence of permissive transcription during spermatogenesis. LncRNAs are often differentially expressed among developmental stages and are less abundant in embryos and newborns compared with adult individuals, in agreement with a requirement for tighter expression control and less tolerance for noisy transcription early in development. For differentially expressed lncRNAs, we find that the patterns of expression variation among developmental stages are generally conserved between mouse and rat. Moreover, lncRNAs expressed above noise levels in somatic organs and during development show higher evolutionary conservation, in particular, at their promoter regions. Thus, we show that functionally constrained lncRNA loci are enriched in developing organs, and we suggest that many of these loci may function in an RNA-independent manner.

The Developmental Gene Hypothesis for Punctuated Equilibrium: Combined Roles of Developmental Regulatory Genes and Transposable Elements

Theories of the genetics underlying punctuated equilibrium (PE) have been vague to date. Here the developmental gene hypothesis is proposed, which states that: 1) developmental regulatory (DevReg) genes are responsible for the orchestration of metazoan morphogenesis and their extreme conservation and mutation intolerance generates the equilibrium or stasis present throughout much of the fossil record and 2) the accumulation of regulatory elements and recombination within these same genes-often derived from transposable elements-drives punctuated bursts of morphological divergence and speciation across metazoa. This two-part hypothesis helps to explain the features that characterize PE, providing a theoretical genetic basis for the once-controversial theory. Also see the video abstract here https://youtu.be/C-fu-ks5yDs.

Compensating for over-production inhibition of the Hsmar1 transposon in Escherichia coli using a series of constitutive promoters

Background

Transposable elements (TEs) are a diverse group of self-mobilizing DNA elements. Transposition has been exploited as a powerful tool for molecular biology and genomics. However, transposition is sometimes limited because of auto-regulatory mechanisms that presumably allow them to cohabit within their hosts without causing excessive genomic damage. The papillation assay provides a powerful visual screen for hyperactive transposases. Transposition is revealed by the activation of a promoter-less lacZ gene when the transposon integrates into a non-essential gene on the host chromosome. Transposition events are detected as small blue speckles, or papillae, on the white background of the main Escherichia coli colony.

Results

We analysed the parameters of the papillation assay including the strength of the transposase transcriptional and translational signals. To overcome certain limitations of inducible promoters, we constructed a set of vectors based on constitutive promoters of different strengths to widen the range of transposase expression. We characterized and validated our expression vectors with Hsmar1, a member of the mariner transposon family. The highest rate of transposition was observed with the weakest promoters. We then took advantage of our approach to investigate how the level of transposition responds to selected point mutations and the effect of joining the transposase monomers into a single-chain dimer.

Conclusions

We generated a set of vectors to provide a wide range of transposase expression which will be useful for screening libraries of transposase mutants. The use of weak promoters should allow screening for truly hyperactive transposases rather than those that are simply resistant to auto-regulatory mechanisms, such as overproduction inhibition (OPI). We also found that mutations in the Hsmar1 dimer interface provide resistance to OPI in bacteria, which could be valuable for improving bacterial transposon mutagenesis techniques.

Transdifferentiation of Mouse Embryonic Fibroblasts into Dopaminergic Neurons Reactivates LINE-1 Repetitive Elements

In mammals, LINE-1 (L1) retrotransposons constitute between 15% and 20% of the genome. Although only a few copies have retained the ability to retrotranspose, evidence in brain and differentiating pluripotent cells indicates that L1 retrotransposition occurs and creates mosaics in normal somatic tissues. The function of de novo insertions remains to be understood. The transdifferentiation of mouse embryonic fibroblasts to dopaminergic neuronal fate provides a suitable model for studying L1 dynamics in a defined genomic and unaltered epigenomic background. We found that L1 elements are specifically re-expressed and mobilized during the initial stages of reprogramming and that their insertions into specific acceptor loci coincides with higher chromatin accessibility and creation of new transcribed units. Those events accompany the maturation of neuronal committed cells. We conclude that L1 retrotransposition is a non-random process correlating with chromatin opening and lncRNA production that accompanies direct somatic cell reprogramming.

•

L1 activation accompanies induced dopaminergic neuron maturation

•

L1 inhibition impairs the transdifferentiation potential of MEFs

•

L1 retrotransposition creates a lineage-specific genetic mosaicism

•

L1 insertions correlates with open chromatin and lncRNA transcription

Beyond classic editing: innovative CRISPR approaches for functional studies of long non-coding RNA

Long non-coding RNAs (lncRNAs) makeup a considerable part of the non-coding human genome and had been well-established as crucial players in an array of biological processes. In spite of their abundance and versatile roles, their functional characteristics remain largely undiscovered mainly due to the lack of suitable genetic manipulation tools. The emerging CRISPR/Cas9 technology has been widely adapted in several studies that aim to screen and identify novel lncRNAs as well as interrogate the functional properties of specific lncRNAs. However, the complexity of lncRNAs genes and the regulatory mechanisms that govern their transcription, as well as their unique functionality pose several limitations the utilization of classic CRISPR methods in lncRNAs functional studies. Here, we overview the unique characteristics of lncRNAs transcription and function and the suitability of the CRISPR toolbox for applications in functional characterization of lncRNAs. We discuss some of the novel variations to the classic CRISPR/Cas9 system that have been tailored and applied previously to study several aspects of lncRNAs functionality. Finally, we share perspectives on the potential applications of various CRISPR systems, including RNA-targeting, in the direct editing and manipulation of lncRNAs.

A primate-specific retroviral enhancer wires the XACT lncRNA into the core pluripotency network in humans

Transposable elements (TEs) have been proposed to play an important role in driving the expansion of gene regulatory networks during mammalian evolution, notably by contributing to the evolution and function of long non-coding RNAs (lncRNAs). XACT is a primate-specific TE-derived lncRNA that coats active X chromosomes in pluripotent cells and may contribute to species-specific regulation of X-chromosome inactivation. Here we explore how different families of TEs have contributed to shaping the XACT locus and coupling its expression to pluripotency. Through a combination of sequence analysis across primates, transcriptional interference, and genome editing, we identify a critical enhancer for the regulation of the XACT locus that evolved from an ancestral group of mammalian endogenous retroviruses (ERVs), prior to the emergence of XACT. This ERV was hijacked by younger hominoid-specific ERVs that gave rise to the promoter of XACT, thus wiring its expression to the pluripotency network. This work illustrates how retroviral-derived sequences may intervene in species-specific regulatory pathways.

Maize transposable elements contribute to long non-coding RNAs that are regulatory hubs for abiotic stress response

Background

Several studies have mined short-read RNA sequencing datasets to identify long non-coding RNAs (lncRNAs), and others have focused on the function of individual lncRNAs in abiotic stress response. However, our understanding of the complement, function and origin of lncRNAs – and especially transposon derived lncRNAs (TE-lncRNAs) - in response to abiotic stress is still in its infancy.

Results

We utilized a dataset of 127 RNA sequencing samples that included total RNA datasets and PacBio fl-cDNA data to discover lncRNAs in maize. Overall, we identified 23,309 candidate lncRNAs from polyA+ and total RNA samples, with a strong discovery bias within total RNA. The majority (65%) of the 23,309 lncRNAs had sequence similarity to transposable elements (TEs). Most had similarity to long-terminal-repeat retrotransposons from the Copia and Gypsy superfamilies, reflecting a high proportion of these elements in the genome. However, DNA transposons were enriched for lncRNAs relative to their genomic representation by ~ 2-fold. By assessing the fraction of lncRNAs that respond to abiotic stresses like heat, cold, salt and drought, we identified 1077 differentially expressed lncRNA transcripts, including 509 TE-lncRNAs. In general, the expression of these lncRNAs was significantly correlated with their nearest gene. By inferring co-expression networks across our large dataset, we found that 39 lncRNAs are as major hubs in co-expression networks that respond to abiotic stress, and 18 appear to be derived from TEs.

Conclusions

Our results show that lncRNAs are enriched in total RNA samples, that most (65%) are derived from TEs, that at least 1077 are differentially expressed during abiotic stress, and that 39 are hubs in co-expression networks, including a small number that are evolutionary conserved. These results suggest that lncRNAs, including TE-lncRNAs, may play key regulatory roles in moderating abiotic responses.

The long non-coding RNA PCA3: an update of its functions and clinical applications as a biomarker in prostate cancer

Prostate cancer antigen 3 (PCA3) is an overexpressed prostate long non-coding RNA (lncRNA), transcribed from an intronic region at the long arm of human chromosome 9q21–22. It has been described that PCA3 modulates prostate cancer (PCa) cell survival through modulating androgen receptor (AR) signaling, besides controlling the expression of several androgen responsive and cancer-related genes, including epithelial–mesenchymal transition (EMT) markers and those regulating gene expression and cell signaling. Also, PCA3 urine levels have been successfully used as a PCa diagnostic biomarker. In this review, we have highlighted recent findings regarding PCA3, addressing its gene structure, putative applications as a biomarker, a proposed origin of this lncRNA, roles in PCa biology and expression patterns. We also updated data regarding PCA3 interactions with cancer-related miRNAs and expression in other tissues and diseases beyond the prostate. Altogether, literature data indicate aberrant expression and dysregulated activity of PCA3, suggesting PCA3 as a promising relevant target that should be even further evaluated on its applicability for PCa detection and management.

High-Quality Genome Assemblies Reveal Long Non-coding RNAs Expressed in Ant Brains

Ants are an emerging model system for neuroepigenetics, as embryos with virtually identical genomes develop into different adult castes that display diverse physiology, morphology, and behavior. Although a number of ant genomes have been sequenced to date, their draft quality is an obstacle to sophisticated analyses of epigenetic gene regulation. We reassembled de novo high-quality genomes for two ant species, Camponotus floridanus and Harpegnathos saltator. Using long reads enabled us to span large repetitive regions and improve genome contiguity, leading to comprehensive and accurate protein-coding annotations that facilitated the identification of a Gp-9-like gene as differentially expressed in Harpegnathos castes. The new assemblies also enabled us to annotate long non-coding RNAs in ants, revealing caste-, brain-, and developmental-stage-specific long non-coding RNAs (lncRNAs) in Harpegnathos. These upgraded genomes, along with the new gene annotations, will aid future efforts to identify epigenetic mechanisms of phenotypic and behavioral plasticity in ants.

Long Noncoding RNAs and Repetitive Elements: Junk or Intimate Evolutionary Partners?

Our recent ability to sequence entire genomes, along with all of their transcribed RNAs, has led to the surprising finding that only ∼1% of the human genome is used to encode proteins. This finding has led to vigorous debate over the functional importance of the transcribed but untranslated portions of the genome. Currently, scientists tend to assume coding genes are functional until proven not to be, while the opposite is true for noncoding genes. This review takes a new look at the evidence for and against widespread noncoding gene functionality. We focus in particular on long noncoding RNA (noncoding RNAs longer than 200 nucleotides) genes and their 'junk' associates, transposable elements, and satellite repeats. Taken together, the suggestion put forward is that more of this junk DNA may be functional than nonfunctional and that noncoding RNAs and transposable elements act symbiotically to drive evolution.

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

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

Continuous salt stress-induced long non-coding RNAs and DNA methylation patterns in soybean roots

BACKGROUND

Environmental stimuli can activate a series of physiological and biochemical responses in plants accompanied by extensive transcriptional reprogramming. Long non-coding RNAs (lncRNAs), as versatile regulators, control gene expression in multiple ways and participate in the adaptation to biotic and abiotic stresses.

RESULTS

In this study, soybean seedlings were continuously cultured for 15 days with high salinity solutions started from seed germination. Strand-specific whole transcriptome sequencing and stringent bioinformatic analysis led to the identification of 3030 long intergenic non-coding RNAs (lincRNAs) and 275 natural antisense transcripts (lncNATs) in soybean roots. In contrast to mRNAs, newly identified lncRNAs exhibited less exons, similar AU content to UTRs, even distribution across the genome and low evolutionary conservation. Remarkably, more than 75% of discovered lncRNAs that were activated or up-regulated by continuous salt stress mainly targeted proteins with binding and catalytic activities. Furthermore, two DNA methylation maps with single-base resolution were generated by using reduced representation bisulfite sequencing, offering a genome-wide perspective and important clues for epigenetic regulation of stress-associated lncRNAs and protein-coding genes.

CONCLUSIONS

Taken together, our findings systematically demonstrated the characteristics of continuous salt stress-induced lncRNAs and extended the knowledge of corresponding methylation profiling, providing valuable evidence for a better understanding of how plants cope with long-term salt stress circumstances.

Interrogation of nonconserved human adipose lincRNAs identifies a regulatory role of linc-ADAL in adipocyte metabolism

Long intergenic noncoding RNAs (lincRNAs) have emerged as important modulators of cellular functions. Most lincRNAs are not conserved among mammals, raising the fundamental question of whether nonconserved adipose-expressed lincRNAs are functional. To address this, we performed deep RNA sequencing of gluteal subcutaneous adipose tissue from 25 healthy humans. We identified 1001 putative lincRNAs expressed in all samples through de novo reconstruction of noncoding transcriptomes and integration with existing lincRNA annotations. One hundred twenty lincRNAs had adipose-enriched expression, and 54 of these exhibited peroxisome proliferator-activated receptor γ (PPARγ) or CCAAT/enhancer binding protein α (C/EBPα) binding at their loci. Most of these adipose-enriched lincRNAs (~85%) were not conserved in mice, yet on average, they showed degrees of expression and binding of PPARγ and C/EBPα similar to those displayed by conserved lincRNAs. Most adipose lincRNAs differentially expressed (n = 53) in patients after bariatric surgery were nonconserved. The most abundant adipose-enriched lincRNA in our subcutaneous adipose data set, linc-ADAL, was nonconserved, up-regulated in adipose depots of obese individuals, and markedly induced during in vitro human adipocyte differentiation. We demonstrated that linc-ADAL interacts with heterogeneous nuclear ribonucleoprotein U (hnRNPU) and insulin-like growth factor 2 mRNA binding protein 2 (IGF2BP2) at distinct subcellular locations to regulate adipocyte differentiation and lipogenesis.

The control of gene expression and cell identity by H3K9 trimethylation

Histone 3 lysine 9 trimethylation (H3K9me3) is a conserved histone modification that is best known for its role in constitutive heterochromatin formation and the repression of repetitive DNA elements. More recently, it has become evident that H3K9me3 is also deposited at certain loci in a tissue-specific manner and plays important roles in regulating cell identity. Notably, H3K9me3 can repress genes encoding silencing factors, pointing to a fundamental principle of repressive chromatin auto-regulation. Interestingly, recent studies have shown that H3K9me3 deposition requires protein SUMOylation in different contexts, suggesting that the SUMO pathway functions as an important module in gene silencing and heterochromatin formation. In this Review, we discuss the role of H3K9me3 in gene regulation in various systems and the molecular mechanisms that guide the silencing machinery to target loci.

Transcriptome analyses of tumor-adjacent somatic tissues reveal genes co-expressed with transposable elements

BACKGROUND

Despite the long-held assumption that transposons are normally only expressed in the germ-line, recent evidence shows that transcripts of transposable element (TE) sequences are frequently found in the somatic cells. However, the extent of variation in TE transcript levels across different tissues and different individuals are unknown, and the co-expression between TEs and host gene mRNAs have not been examined.

RESULTS

Here we report the variation in TE derived transcript levels across tissues and between individuals observed in the non-tumorous tissues collected for The Cancer Genome Atlas. We found core TE co-expression modules consisting mainly of transposons, showing correlated expression across broad classes of TEs. Despite this co-expression within tissues, there are individual TE loci that exhibit tissue-specific expression patterns, when compared across tissues. The core TE modules were negatively correlated with other gene modules that consisted of immune response genes in interferon signaling. KRAB Zinc Finger Proteins (KZFPs) were over-represented gene members of the TE modules, showing positive correlation across multiple tissues. But we did not find overlap between TE-KZFP pairs that are co-expressed and TE-KZFP pairs that are bound in published ChIP-seq studies.

CONCLUSIONS

We find unexpected variation in TE derived transcripts, within and across non-tumorous tissues. We describe a broad view of the RNA state for non-tumorous tissues exhibiting higher level of TE transcripts. Tissues with higher level of TE transcripts have a broad range of TEs co-expressed, with high expression of a large number of KZFPs, and lower RNA levels of immune genes.

Host-transposon interactions: conflict, cooperation, and cooption

Transposable elements (TEs) are mobile DNA sequences that colonize genomes and threaten genome integrity. As a result, several mechanisms appear to have emerged during eukaryotic evolution to suppress TE activity. However, TEs are ubiquitous and account for a prominent fraction of most eukaryotic genomes. We argue that the evolutionary success of TEs cannot be explained solely by evasion from host control mechanisms. Rather, some TEs have evolved commensal and even mutualistic strategies that mitigate the cost of their propagation. These coevolutionary processes promote the emergence of complex cellular activities, which in turn pave the way for cooption of TE sequences for organismal function.

Human SETMAR is a DNA sequence-specific histone-methylase with a broad effect on the transcriptome

Transposons impart dynamism to the genomes they inhabit and their movements frequently rewire the control of nearby genes. Occasionally, their proteins are domesticated when they evolve a new function. SETMAR is a protein methylase with a sequence-specific DNA binding domain. It began to evolve about 50 million years ago when an Hsmar1 transposon integrated downstream of a SET-domain methylase gene. Here we show that the DNA-binding domain of the transposase targets the enzyme to transposon-end remnants and that this is capable of regulating gene expression, dependent on the methylase activity. When SETMAR was modestly overexpressed in human cells, almost 1500 genes changed expression by more than 2-fold (65% up- and 35% down-regulated). These genes were enriched for the KEGG Pathways in Cancer and include several transcription factors important for development and differentiation. Expression of a similar level of a methylase-deficient SETMAR changed the expression of many fewer genes, 77% of which were down-regulated with no significant enrichment of KEGG Pathways. Our data is consistent with a model in which SETMAR is part of an anthropoid primate-specific regulatory network centered on the subset of genes containing a transposon end.

Parallels between artificial reprogramming and the biogenesis of cancer stem cells: Involvement of lncRNAs

Cellular identity is established and maintained by the interplay of cell type-specific transcription factors and epigenetic regulation of the genome. During development in vivo and differentiation in vitro, transitions from one cell type to the next are triggered by cell signaling events culminating in modifications of chromatin that render genes accessible or inaccessible to the transcriptional apparatus. In recent years it has become apparent that cellular identity is plastic, and technological reprogramming methods such as somatic cell nuclear transfer and induced pluripotency can yield reprogrammed cells that have been restored to a state of developmental potency. Long noncoding RNAs (lncRNAs) are untranslated functional RNA molecules that are intimately involved in the regulation of the chromatin of protein-coding genes. In fact, recent evidence shows that there are more lncRNA species in the cell than mRNA species and that most protein-coding genes are likely to be under epigenetic regulation mediated by lncRNAs. This review examines lncRNA function in reprogrammed pluripotent cells and cancer stem cells. Because cancer stem cells arise from normal cells, their biogenesis can be viewed as a reprogramming process that occurs in vivo, and parallels between artificial reprogramming and cancer stem cell biogenesis are discussed.

Prediction of regulatory long intergenic non-coding RNAs acting in trans through base-pairing interactions

Background

Long intergenic non-coding RNAs (lincRNAs) can act as regulators of expression of protein-coding genes. Trans-natural antisense transcripts (trans-NATs) are a type of lincRNAs that contain sequence complementary to mRNA from other loci. The regulatory potential of trans-NATs has been poorly studied in eukaryotes and no example of trans-NATs regulating gene expression in plants are reported. The goal of this study was to identify lincRNAs, and particularly trans-NATs, in Arabidopsis thaliana that have a potential to regulate expression of target genes in trans at the transcriptional or translational level.

Results

We identified 1001 lincRNAs using an RNAseq dataset from total polyA⁺ and polysome-associated RNA of seedlings grown under high and low phosphate, or shoots and roots treated with different phytohormones, of which 550 were differentially regulated. Approximately 30% of lincRNAs showed conservation amongst Brassicaceae and 25% harbored transposon element (TE) sequences. Gene co-expression network analysis highlighted a group of lincRNAs associated with the response of roots to low phosphate. A total of 129 trans-NATs were predicted, of which 88 were significantly differentially expressed under at least one pairwise comparison. Five trans-NATs showed a positive correlation between their expression and target mRNA steady-state levels, and three showed a negative correlation. Expression of four trans-NATs positively correlated with a change in target mRNA polysome association. The regulatory potential of these trans-NATs did not implicate miRNA mimics nor siRNAs. We also looked for lincRNAs that could regulate gene expression in trans by Watson-Crick DNA:RNA base pairing with target protein-encoding loci. We identified 100 and 81 with a positive or negative correlation, respectively, with steady-state level of their predicted target. The regulatory potential of one such candidate lincRNA harboring a SINE TE sequence was validated in a protoplast assay on three distinct genes containing homologous TE sequence in their promoters. Construction of networks highlighted other putative lincRNAs with multiple predicted target loci for which expression was positively correlated with target gene expression.

Conclusions

This study identified lincRNAs in Arabidopsis with potential in regulating target gene expression in trans by both RNA:RNA and RNA:DNA base pairing and highlights lincRNAs harboring TE sequences in such activity.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5946-0) contains supplementary material, which is available to authorized users.

The role of transposable elements in the ecological morphogenesis under the influence of stress

In natural selection, insertional mutagenesis is an important source of genome variability. Transposons are sensors of environmental stress effects, which contribute to adaptation and speciation. These effects are due to changes in the mechanisms of morphogenesis, since transposons contain regulatory sequences that have cis and trans effects on specific protein-coding genes. In variability of genomes, the horizontal transfer of transposons plays an important role, because it contributes to changing the composition of transposons and the acquisition of new properties. Transposons are capable of site-specific transpositions, which lead to the activation of stress response genes. Transposons are sources of non-coding RNA, transcription factors binding sites and protein-coding genes due to domestication, exonization, and duplication. These genes contain nucleotide sequences that interact with non-coding RNAs processed from transposons transcripts, and therefore they are under the control of epigenetic regulatory networks involving transposons. Therefore, inherited features of the location and composition of transposons, along with a change in the phenotype, play an important role in the characteristics of responding to a variety of environmental stressors. This is the basis for the selection and survival of organisms with a specific composition and arrangement of transposons that contribute to adaptation under certain environmental conditions. In evolution, the capability to transpose into specific genome sites, regulate gene expression, and interact with transcription factors, along with the ability to respond to stressors, is the basis for rapid variability and speciation by altering the regulation of ontogenesis. The review presents evidence of tissue-specific and stage-specific features of transposon activation and their role in the regulation of cell differentiation to confirm their role in ecological morphogenesis.

Long-read direct RNA sequencing by 5'-Cap capturing reveals the impact of Piwi on the widespread exonization of transposable elements in locusts

The large genome of the migratory locust (Locusta migratoria) genome accumulates massive amount of accumulated transposable elements (TEs), which show intrinsic transcriptional activities. Hampering the ability to precisely determine full-length RNA transcript sequences are exonized TEs, which produce numerous highly similar fragments that are difficult to resolve using short-read sequencing technology. Here, we applied a 5'-Cap capturing method using Nanopore long-read direct RNA sequencing to characterize full-length transcripts in their native RNA form and to analyze the TE exonization pattern in the locust transcriptome. Our results revealed the widespread establishment of TE exonization and a substantial contribution of TEs to RNA splicing in the locust transcriptome. The results of the transcriptomic spectrum influenced by Piwi expression indicated that TE-derived sequences were the main targets of Piwi-mediated repression. Furthermore, our study showed that Piwi expression regulates the length of RNA transcripts containing TE-derived sequences, creating an alternative UTR usage. Overall, our results reveal the transcriptomic characteristics of TE exonization in the species characterized by large and repetitive genomes.

Global Positioning System: Understanding Long Noncoding RNAs through Subcellular Localization

The localization of long noncoding RNAs (lncRNAs) within the cell is the primary determinant of their molecular functions. LncRNAs are often thought of as chromatin-restricted regulators of gene transcription and chromatin structure. However, a rich population of cytoplasmic lncRNAs has come to light, with diverse roles including translational regulation, signaling, and respiration. RNA maps of increasing resolution and scope are revealing a subcellular world of highly specific localization patterns and hint at sequence-based address codes specifying lncRNA fates. We propose a new framework for analyzing sequencing-based data, which suggests that numbers of cytoplasmic lncRNA molecules rival those in the nucleus. New techniques promise to create high-resolution, transcriptome-wide maps associated with all organelles of the mammalian cell. Given its intimate link to molecular roles, subcellular localization provides a means of unlocking the mystery of lncRNA functions.

Novel Bioinformatics Approach Identifies Transcriptional Profiles of Lineage-Specific Transposable Elements at Distinct Loci in the Human Dorsolateral Prefrontal Cortex

Expression of transposable elements (TE) is transiently activated during human preimplantation embryogenesis in a developmental stage- and cell type-specific manner and TE-mediated epigenetic regulation is intrinsically wired in developmental genetic networks in human embryos and embryonic stem cells. However, there are no systematic studies devoted to a comprehensive analysis of the TE transcriptome in human adult organs and tissues, including human neural tissues. To investigate TE expression in the human Dorsolateral Prefrontal Cortex (DLPFC), we developed and validated a straightforward analytical approach to chart quantitative genome-wide expression profiles of all annotated TE loci based on unambiguous mapping of discrete TE-encoded transcripts using a de novo assembly strategy. To initially evaluate the potential regulatory impact of DLPFC-expressed TE, we adopted a comparative evolutionary genomics approach across humans, primates, and rodents to document conservation patterns, lineage-specificity, and colocalizations with transcription factor binding sites mapped within primate- and human-specific TE. We identified 654,665 transcripts expressed from 477,507 distinct loci of different TE classes and families, the majority of which appear to have originated from primate-specific sequences. We discovered 4,687 human-specific and transcriptionally active TEs in DLPFC, of which the prominent majority (80.2%) appears spliced. Our analyses revealed significant associations of DLPFC-expressed TE with primate- and human-specific transcription factor binding sites, suggesting potential cross-talks of concordant regulatory functions. We identified 1,689 TEs differentially expressed in the DLPFC of Schizophrenia patients, a majority of which is located within introns of 1,137 protein-coding genes. Our findings imply that identified DLPFC-expressed TEs may affect human brain structures and functions following different evolutionary trajectories. On one side, hundreds of thousands of TEs maintained a remarkably high conservation for ∼8 My of primates’ evolution, suggesting that they are likely conveying evolutionary-constrained primate-specific regulatory functions. In parallel, thousands of transcriptionally active human-specific TE loci emerged more recently, suggesting that they could be relevant for human-specific behavioral or cognitive functions.

In silico prediction of housekeeping long intergenic non-coding RNAs reveals HKlincR1 as an essential player in lung cancer cell survival

Prioritising long intergenic noncoding RNAs (lincRNAs) for functional characterisation is a significant challenge. Here we applied computational approaches to discover lincRNAs expected to play a critical housekeeping (HK) role within the cell. Using the Illumina Human BodyMap RNA sequencing dataset as a starting point, we first identified lincRNAs ubiquitously expressed across a panel of human tissues. This list was then further refined by reference to conservation score, secondary structure and promoter DNA methylation status. Finally, we used tumour expression and copy number data to identify lincRNAs rarely downregulated or deleted in multiple tumour types. The resulting list of candidate essential lincRNAs was then subjected to co-expression analyses using independent data from ENCODE and The Cancer Genome Atlas (TCGA). This identified a substantial subset with a predicted role in DNA replication and cell cycle regulation. One of these, HKlincR1, was selected for further characterisation. Depletion of HKlincR1 affected cell growth in multiple lung cancer cell lines, and led to disruption of genes involved in cell growth and viability. In addition, HKlincR1 expression was correlated with overall survival in lung adenocarcinoma patients. Our in silico studies therefore reveal a set of housekeeping noncoding RNAs of interest both in terms of their role in normal homeostasis, and their relevance in tumour growth and maintenance.

Retrotransposons evolution and impact on lncRNA and protein coding genes in pigs

BACKGROUND

Retrotransposons are the major determinants of genome sizes and they have shaped both genes and genomes in mammalian organisms, but their overall activity, diversity, and evolution dynamics, particularly their impact on protein coding and lncRNA genes in pigs remain largely unknown.

RESULTS

In the present study, we performed de novo detection of retrotransposons in pigs by using multiple pipelines, four distinct families of pig-specific L1 s classified into 51 distinct subfamilies and representing four evolution models and three expansion waves of pig-specific SINEs represented by three distinct families were identified. ERVs were classified into 18 families and found two most "modern" subfamilies in the pig genome. The transposition activity of pig L1 was verified by experiment, the sense and antisense promoter activities of young L1 5'UTRs and ERV LTRs and expression profiles of young retrotransposons in multiple tissues and cell lines were also validated. Furthermore, retrotransposons had an extensive impact on lncRNA and protein coding genes at both the genomic and transcriptomic levels. Most protein coding and lncRNA (> 80%) genes contained retrotransposon insertions, and about half of protein coding genes (44.30%) and one-fourth (24.13%) of lncRNA genes contained the youngest retrotransposon insertions. Nearly half of protein coding genes (43.78%) could generate chimeric transcripts with retrotransposons. Significant distribution bias of retrotransposon composition, location, and orientation in lncRNA and protein coding genes, and their transcripts, were observed.

CONCLUSIONS

In the current study, we characterized the classification and evolution profile of retrotransposons in pigs, experimentally proved the transposition activity of the young pig L1 subfamily, characterized the sense and antisense expression profiles and promoter activities of young retrotransposons, and investigated their impact on lncRNA and protein coding genes by defining the mobilome landscapes at the genomic and transcriptomic levels. These findings help provide a better understanding of retrotransposon evolution in mammal and their impact on the genome and transcriptome.

Distribution, Characteristics, and Regulatory Potential of Long Noncoding RNAs in Brown-Rot Fungi

Long noncoding RNAs have been thoroughly studied in plants, animals, and yeasts, where they play important roles as regulators of transcription. Nevertheless, almost nothing is known about their presence and characteristics in filamentous fungi, especially in basidiomycetes. In the present study, we have carried out an exhaustive annotation and characterization of lncRNAs in two lignin degrader basidiomycetes, Coniophora puteana and Serpula lacrymans. We identified 2,712 putative lncRNAs in the former and 2,242 in the latter, mainly originating from intergenic locations of transposon-sparse genomic regions. The lncRNA length, GC content, expression levels, and stability of the secondary structure differ from coding transcripts but are similar in these two species and resemble that of other eukaryotes. Nevertheless, they lack sequence conservation. Also, we found that lncRNAs are transcriptionally regulated in the same proportion as genes when the fungus actively decomposes soil organic matter. Finally, up to 7% of the upstream gene regions of Coniophora puteana and Serpula lacrymans are transcribed and produce lncRNAs. The study of expression trends in these gene-lncRNA pairs uncovered groups with similar and opposite transcriptional profiles which may be the result of cis-transcriptional regulation.

New lncRNA annotation reveals extensive functional divergence of the transcriptome in maize

Long non-coding RNAs (lncRNAs), whose sequences are approximately 200 bp or longer and unlikely to encode proteins, may play an important role in eukaryotic gene regulation. Although the latest maize (Zea mays L.) reference genome provides an essential genomic resource, genome-wide annotations of maize lncRNAs have not been updated. Here, we report on a large transcriptomic dataset collected from 749 RNA sequencing experiments across different tissues and stages of the maize reference inbred B73 line and 60 from its wild relative teosinte. We identified 18,165 high-confidence lncRNAs in maize, of which 6,873 are conserved between maize and teosinte. We uncovered distinct genomic characteristics of conserved lncRNAs, non-conserved lncRNAs, and protein-coding transcripts. Intriguingly, Shannon entropy analysis showed that conserved lncRNAs are likely to be expressed similarly to protein-coding transcripts. Co-expression network analysis revealed significant variation in the degree of co-expression. Furthermore, selection analysis indicated that conserved lncRNAs are more likely than non-conserved lncRNAs to be located in regions subject to recent selection, indicating evolutionary differentiation. Our results provide the latest genome-wide annotation and analysis of maize lncRNAs and uncover potential functional divergence between protein-coding, conserved lncRNA, and non-conserved lncRNA genes, demonstrating the high complexity of the maize transcriptome.

The Genome of Cucurbita argyrosperma (Silver-Seed Gourd) Reveals Faster Rates of Protein-Coding Gene and Long Noncoding RNA Turnover and Neofunctionalization within Cucurbita

Whole-genome duplications are an important source of evolutionary novelties that change the mode and tempo at which genetic elements evolve within a genome. The Cucurbita genus experienced a whole-genome duplication around 30 million years ago, although the evolutionary dynamics of the coding and noncoding genes in this genus have not yet been scrutinized. Here, we analyzed the genomes of four Cucurbita species, including a newly assembled genome of Cucurbita argyrosperma, and compared the gene contents of these species with those of five other members of the Cucurbitaceae family to assess the evolutionary dynamics of protein-coding and long intergenic noncoding RNA (lincRNA) genes after the genome duplication. We report that Cucurbita genomes have a higher protein-coding gene birth-death rate compared with the genomes of the other members of the Cucurbitaceae family. C. argyrosperma gene families associated with pollination and transmembrane transport had significantly faster evolutionary rates. lincRNA families showed high levels of gene turnover throughout the phylogeny, and 67.7% of the lincRNA families in Cucurbita showed evidence of birth from the neofunctionalization of previously existing protein-coding genes. Collectively, our results suggest that the whole-genome duplication in Cucurbita resulted in faster rates of gene family evolution through the neofunctionalization of duplicated genes.

That is life: communicating RNA networks from viruses and cells in continuous interaction

All the conserved detailed results of evolution stored in DNA must be read, transcribed, and translated via an RNA-mediated process. This is required for the development and growth of each individual cell. Thus, all known living organisms fundamentally depend on these RNA-mediated processes. In most cases, they are interconnected with other RNAs and their associated protein complexes and function in a strictly coordinated hierarchy of temporal and spatial steps (i.e., an RNA network). Clearly, all cellular life as we know it could not function without these key agents of DNA replication, namely rRNA, tRNA, and mRNA. Thus, any definition of life that lacks RNA functions and their networks misses an essential requirement for RNA agents that inherently regulate and coordinate (communicate to) cells, tissues, organs, and organisms. The precellular evolution of RNAs occurred at the core of the emergence of cellular life and the question remained of how both precellular and cellular levels are interconnected historically and functionally. RNA networks and RNA communication can interconnect these levels. With the reemergence of virology in evolution, it became clear that communicating viruses and subviral infectious genetic parasites are bridging these two levels by invading, integrating, coadapting, exapting, and recombining constituent parts in host genomes for cellular requirements in gene regulation and coordination aims. Therefore, a 21st century understanding of life is of an inherently social process based on communicating RNA networks, in which viruses and cells continuously interact.

Evolutionary Origins of Pseudogenes and Their Association with Regulatory Sequences in Plants

Pseudogenes (Ψs), nonfunctional relatives of functional genes, form by duplication or retrotransposition, and loss of gene function by disabling mutations. Evolutionary analysis provides clues to Ψ origins and effects on gene regulation. However, few systematic studies of plant Ψs have been conducted, hampering comparative analyses. Here, we examined the origin, evolution, and expression patterns of Ψs and their relationships with noncoding sequences in seven angiosperm plants. We identified ∼250,000 Ψs, most of which are more lineage specific than protein-coding genes. The distribution of Ψs on the chromosome indicates that genome recombination may contribute to Ψ elimination. Most Ψs evolve rapidly in terms of sequence and expression levels, showing tissue- or stage-specific expression patterns. We found that a surprisingly large fraction of nontransposable element regulatory noncoding RNAs (microRNAs and long noncoding RNAs) originate from transcription of Ψ proximal upstream regions. We also found that transcription factor binding sites preferentially occur in putative Ψ proximal upstream regions compared with random intergenic regions, suggesting that Ψs have conditioned genome evolution by providing transcription factor binding sites that serve as promoters and enhancers. We therefore propose that rapid rewiring of Ψ transcriptional regulatory regions is a major mechanism driving the origin of novel regulatory modules.

RNAi-dependent Polycomb repression controls transposable elements in Tetrahymena

RNAi and Polycomb repression play evolutionarily conserved and often coordinated roles in transcriptional silencing. Here, we show that, in the protozoan Tetrahymena thermophila, germline-specific internally eliminated sequences (IESs)-many related to transposable elements (TEs)-become transcriptionally activated in mutants deficient in the RNAi-dependent Polycomb repression pathway. Germline TE mobilization also dramatically increases in these mutants. The transition from noncoding RNA (ncRNA) to mRNA production accompanies transcriptional activation of TE-related sequences and vice versa for transcriptional silencing. The balance between ncRNA and mRNA production is potentially affected by cotranscriptional processing as well as RNAi and Polycomb repression. We posit that interplay between RNAi and Polycomb repression is a widely conserved phenomenon, whose ancestral role is epigenetic silencing of TEs.

LIONS: analysis suite for detecting and quantifying transposable element initiated transcription from RNA-seq

Summary

Transposable elements (TEs) influence the evolution of novel transcriptional networks yet the specific and meaningful interpretation of how TE-derived transcriptional initiation contributes to the transcriptome has been marred by computational and methodological deficiencies. We developed LIONS for the analysis of RNA-seq data to specifically detect and quantify TE-initiated transcripts.

Availability and implementation

Source code, container, test data and instruction manual are freely available at www.github.com/ababaian/LIONS.

Supplementary information

Supplementary data are available at Bioinformatics online.

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.

Deciphering DMET genetic data: comprehensive assessment of Northwestern Han, Tibetan, Uyghur populations and their comparison to eleven 1000 genome populations

We investigated the allele frequencies of drug absorption, distribution, metabolism and elimination (ADME)-related drug-metabolizing enzymes and transporters (DMET) genes in the Northwestern Han, Tibetan and Uyghur populations and compared the related genes in these three populations with those in eleven 1000 Genome populations. We examined 1936 single nucleotide polymorphisms of 225 DMET genes involved in ADME processes and found 732, 679 and 804 sites were polymorphic in Han, Tibetan and Uyghur. Tibetan differed from Han in only four sites (p < .05), whereas Uyghur differed from Han and Tibetan in 24 and 21 sites, respectively (p < .05). The distributions of 1058 genotyping data of 245 individuals from Han, Tibetan and Uyghur were compared with 1207 other individuals from the eleven 1000 Genomes populations. The top four populations in Han that exhibited the smallest pairwise Fst values were CHB, Tibetan, CHD and JPT; those in Tibetan were Han, CHB, Uyghur and CHD; and those in Uyghur were Han, Tibetan, GIH and CEU. MEGA results revealed that CHB, CHD, JPT, Han, Tibetan and Uyghur were grouped in cluster 1. GIH, MEX, CEU and TSI were grouped in cluster 2. MKK, ASW, LWK and YRI were grouped in cluster 3.

Target of Rapamycin (TOR) Regulates the Expression of lncRNAs in Response to Abiotic Stresses in Cotton

TOR (Target of Rapamycin) kinase is an evolutionarily conserved protein kinase, which integrates stress-related cues with growth and metabolic outputs. Long non-coding RNAs (lncRNAs) play a vital role in the regulation of eukaryotic genes. However, little is known about TOR's function in regulating the expression of lncRNAs in plants. In this study, four putative homologous genes encoding the TOR protein were identified by utilizing the recently completed cotton genome. Pharmacological experiments with TOR inhibitor AZD8055 and on silencing GhTOR genes resulted in obvious cotton growth retardation, indicating the conserved role of TOR in plant growth. The expression pattern analyses in different tissues reveal that TOR may play a role in root development, and the transcript levels of TOR genes were changed under different stress conditions. Importantly, we found TOR may be a key player in regulating the expression of long non-coding RNAs (lncRNAs). A total of 10,315 lncRNAs were discovered in cotton seedlings, 90.7% of which were long intergenic ncRNAs. Moreover, we identified the differentially expressed lncRNAs, of which 296 were significantly upregulated and 105 were downregulated in TOR inactivated plants. GO and KEGG analyses of differentially expressed lncRNA neighboring genes reveal that these differentially expressed lncRNA-targeted genes are involved in many life processes, including stress response, glutathione, and ribosomes in cotton. A series of differentially expressed lncRNAs potentially involved in plant stress response was identified under TOR inhibition. Collectively, these results suggest that cotton TOR proteins may directly modulate the expression of putative stress-related lncRNAs and eventually play a potential role in the cotton stress response.

The long noncoding RNA CHROME regulates cholesterol homeostasis in primate

The human genome encodes thousands of long non-coding RNAs (lncRNAs), the majority of which are poorly conserved and uncharacterized. Here we identify a primate-specific lncRNA (CHROME), elevated in the plasma and atherosclerotic plaques of individuals with coronary artery disease, that regulates cellular and systemic cholesterol homeostasis. LncRNA CHROME expression is influenced by dietary and cellular cholesterol via the sterol-activated liver X receptor transcription factors, which control genes mediating responses to cholesterol overload. Using gain- and loss-of-function approaches, we show that CHROME promotes cholesterol efflux and HDL biogenesis by curbing the actions of a set of functionally related microRNAs that repress genes in those pathways. CHROME knockdown in human hepatocytes and macrophages increases levels of miR-27b, miR-33a, miR-33b and miR-128, thereby reducing expression of their overlapping target gene networks and associated biologic functions. In particular, cells lacking CHROME show reduced expression of ABCA1, which regulates cholesterol efflux and nascent HDL particle formation. Collectively, our findings identify CHROME as a central component of the non-coding RNA circuitry controlling cholesterol homeostasis in humans.

Diverging RNPs: Toward Understanding lncRNA-Protein Interactions and Functions

RNA-protein interactions are essential to a variety of biological processes. The realization that mammalian genomes are pervasively transcribed brought a tidal wave of tens of thousands of newly identified long noncoding RNAs (lncRNAs) and raised questions about their purpose in cells. The vast majority of lncRNAs have yet to be studied, and it remains to be determined to how many of these transcripts a function can be ascribed. However, results gleaned from studying a handful of these macromolecules have started to reveal common themes of biological function and mechanism of action involving intricate RNA-protein interactions. Some lncRNAs were shown to regulate the chromatin and transcription of distant and neighboring genes in the nucleus, while others regulate the translation or localization of proteins in the cytoplasm. Some lncRNAs were found to be crucial during development, while mutations and aberrant expression of others have been associated with several types of cancer and a plethora of diseases. Over the last few years, the establishment of new technologies has been key in providing the tools to decode the rules governing lncRNA-protein interactions and functions. This chapter will highlight the general characteristics of lncRNAs, their function, and their mode of action, with a special focus on protein interactions. It will also describe the methods at the disposition of scientists to help them cross this next frontier in our understanding of lncRNA biology.

CRlncRC: a machine learning-based method for cancer-related long noncoding RNA identification using integrated features

BACKGROUND

Long noncoding RNAs (lncRNAs) are widely involved in the initiation and development of cancer. Although some computational methods have been proposed to identify cancer-related lncRNAs, there is still a demanding to improve the prediction accuracy and efficiency. In addition, the quick-update data of cancer, as well as the discovery of new mechanism, also underlay the possibility of improvement of cancer-related lncRNA prediction algorithm. In this study, we introduced CRlncRC, a novel Cancer-Related lncRNA Classifier by integrating manifold features with five machine-learning techniques.

RESULTS

CRlncRC was built on the integration of genomic, expression, epigenetic and network, totally in four categories of features. Five learning techniques were exploited to develop the effective classification model including Random Forest (RF), Naïve bayes (NB), Support Vector Machine (SVM), Logistic Regression (LR) and K-Nearest Neighbors (KNN). Using ten-fold cross-validation, we showed that RF is the best model for classifying cancer-related lncRNAs (AUC = 0.82). The feature importance analysis indicated that epigenetic and network features play key roles in the classification. In addition, compared with other existing classifiers, CRlncRC exhibited a better performance both in sensitivity and specificity. We further applied CRlncRC to lncRNAs from the TANRIC (The Atlas of non-coding RNA in Cancer) dataset, and identified 121 cancer-related lncRNA candidates. These potential cancer-related lncRNAs showed a certain kind of cancer-related indications, and many of them could find convincing literature supports.

CONCLUSIONS

Our results indicate that CRlncRC is a powerful method for identifying cancer-related lncRNAs. Machine-learning-based integration of multiple features, especially epigenetic and network features, had a great contribution to the cancer-related lncRNA prediction. RF outperforms other learning techniques on measurement of model sensitivity and specificity. In addition, using CRlncRC method, we predicted a set of cancer-related lncRNAs, all of which displayed a strong relevance to cancer as a valuable conception for the further cancer-related lncRNA function studies.

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.