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

Long and small noncoding RNAs during oocyte-to-embryo transition in mammals

Oocyte-to-embryo transition is a process during which an oocyte ovulates, is fertilized, and becomes a developing embryo. It involves the first major genome reprogramming event in life of an organism where gene expression, which gave rise to a differentiated oocyte, is remodeled in order to establish totipotency in blastomeres of an early embryo. This remodeling involves replacement of maternal RNAs with zygotic RNAs through maternal RNA degradation and zygotic genome activation. This review is focused on expression and function of long noncoding RNAs (lncRNAs) and small RNAs during oocyte-to-embryo transition in mammals. LncRNAs are an assorted rapidly evolving collection of RNAs, which have no apparent protein-coding capacity. Their biogenesis is similar to mRNAs including transcriptional control and post-transcriptional processing. Diverse molecular and biological roles were assigned to lncRNAs although most of them probably did not acquire a detectable biological role. Since some lncRNAs serve as precursors for small noncoding regulatory RNAs in RNA silencing pathways, both types of noncoding RNA are reviewed together.

Long noncoding RNAs in the model species Brachypodium distachyon

Eukaryotic genomes are pervasively transcribed and only a small portion of the transcribed sequences belongs to protein coding genes. High-throughput sequencing technology contributed to consolidate this perspective, allowing the identification of numerous noncoding RNAs with key roles in biological processes. Long noncoding RNAs (lncRNAs) are transcripts longer than 200 nt with limited phylogenetic conservation, expressed at low levels and characterized by tissue/organ specific expression profiles. Although a large set of lncRNAs has been identified, the functional roles of lncRNAs are only beginning to be recognized and the molecular mechanism of lncRNA-mediated gene regulation remains largely unexplored, particularly in plants where their annotation and characterization are still incomplete. Using public and proprietary poly-(A)+ RNA-seq data as well as a collection of full length ESTs from several organs, developmental stages and stress conditions in three Brachypodium distachyon inbred lines, we describe the identification and the main features of thousands lncRNAs. Here we provide a genome-wide characterization of lncRNAs, highlighting their intraspecies conservation and describing their expression patterns among several organs/tissues and stress conditions. This work represents a fundamental resource to deepen our knowledge on long noncoding RNAs in C3 cereals, allowing the Brachypodium community to exploit these results in future research programs.

The functions of long noncoding RNAs in development and stem cells

Eukaryotic genomes are pervasively transcribed, with tens of thousands of RNAs emanating from uni- and bi-directional promoters and from active enhancers. In vertebrates, thousands of loci in each species produce a class of transcripts called long noncoding RNAs (lncRNAs) that are typically expressed at low levels and do not appear to give rise to functional proteins. Substantial numbers of lncRNAs are expressed at specific stages of embryonic development, in many cases from regions flanking key developmental regulators. Here, we review the known biological functions of such lncRNAs and the emerging paradigms of their modes of action. We also provide an overview of the growing arsenal of methods for lncRNA identification, perturbation and functional characterization.

Transposable elements are the primary source of novelty in primate gene regulation

Gene regulation shapes the evolution of phenotypic diversity. We investigated the evolution of liver promoters and enhancers in six primate species using ChIP-seq (H3K27ac and H3K4me1) to profile cis-regulatory elements (CREs) and using RNA-seq to characterize gene expression in the same individuals. To quantify regulatory divergence, we compared CRE activity across species by testing differential ChIP-seq read depths directly measured for orthologous sequences. We show that the primate regulatory landscape is largely conserved across the lineage, with 63% of the tested human liver CREs showing similar activity across species. Conserved CRE function is associated with sequence conservation, proximity to coding genes, cell-type specificity, and transcription factor binding. Newly evolved CREs are enriched in immune response and neurodevelopmental functions. We further demonstrate that conserved CREs bind master regulators, suggesting that while CREs contribute to species adaptation to the environment, core functions remain intact. Newly evolved CREs are enriched in young transposable elements (TEs), including Long-Terminal-Repeats (LTRs) and SINE-VNTR-Alus (SVAs), that significantly affect gene expression. Conversely, only 16% of conserved CREs overlap TEs. We tested the cis-regulatory activity of 69 TE subfamilies by luciferase reporter assays, spanning all major TE classes, and showed that 95.6% of tested TEs can function as either transcriptional activators or repressors. In conclusion, we demonstrated the critical role of TEs in primate gene regulation and illustrated potential mechanisms underlying evolutionary divergence among the primate species through the noncoding genome.

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

Background

Only a small portion of human long non-coding RNAs (lncRNAs) appear to be conserved outside of mammals, but the events underlying the birth of new lncRNAs in mammals remain largely unknown. One potential source is remnants of protein-coding genes that transitioned into lncRNAs.

Results

We systematically compare lncRNA and protein-coding loci across vertebrates, and estimate that up to 5% of conserved mammalian lncRNAs are derived from lost protein-coding genes. These lncRNAs have specific characteristics, such as broader expression domains, that set them apart from other lncRNAs. Fourteen lncRNAs have sequence similarity with the loci of the contemporary homologs of the lost protein-coding genes. We propose that selection acting on enhancer sequences is mostly responsible for retention of these regions. As an example of an RNA element from a protein-coding ancestor that was retained in the lncRNA, we describe in detail a short translated ORF in the JPX lncRNA that was derived from an upstream ORF in a protein-coding gene and retains some of its functionality.

Conclusions

We estimate that ~ 55 annotated conserved human lncRNAs are derived from parts of ancestral protein-coding genes, and loss of coding potential is thus a non-negligible source of new lncRNAs. Some lncRNAs inherited regulatory elements influencing transcription and translation from their protein-coding ancestors and those elements can influence the expression breadth and functionality of these lncRNAs.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1293-0) contains supplementary material, which is available to authorized users.

Expression dynamics and relations with nearby genes of rat transposable elements across 11 organs, 4 developmental stages and both sexes

BACKGROUND

TEs pervade mammalian genomes. However, compared with mice, fewer studies have focused on the TE expression patterns in rat, particularly the comparisons across different organs, developmental stages and sexes. In addition, TEs can influence the expression of nearby genes. The temporal and spatial influences of TEs remain unclear yet.

RESULTS

To evaluate the TEs transcription patterns, we profiled their transcript levels in 11 organs for both sexes across four developmental stages of rat. The results show that most short interspersed elements (SINEs) are commonly expressed in all conditions, which are also the major TE types with commonly expression patterns. In contrast, long terminal repeats (LTRs) are more likely to exhibit specific expression patterns. The expression tendency of TEs and genes are similar in most cases. For example, few specific genes and TEs are in the liver, muscle and heart. However, TEs perform superior over genes on classing organ, which imply their higher organ specificity than genes. By associating the TEs with the closest genes in genome, we find their expression levels are correlated, independent of their distance in some cases.

CONCLUSIONS

TEs sex-dependently associate with nearest genes. A gene would be associated with more than one TE. Our works can help to functionally annotate the genome and further understand the role of TEs in gene regulation.

Long Noncoding RNAs as Diagnostic and Therapeutic Targets in Type 2 Diabetes and Related Complications

Protein-coding genes represent only a small fraction of the human genome. In the past, the majority of the genomic sequence has been considered transcriptionally silent, but recent large-scale studies have uncovered an array of functionally significant elements, including non-protein-coding transcripts, within these noncoding regions of the human genome. Long noncoding RNAs (lncRNAs), a class of noncoding transcripts with lengths >200 nucleotides, are pervasively transcribed in the genome and function as signals, decoys, guides, or scaffolds to regulate gene expression. More than 200 diseases have been associated with dysregulated or dysfunctional lncRNAs, and new associations continue to accumulate in the literature. The role of lncRNAs in the pathogenesis of type 2 diabetes mellitus and related complications has only recently been recognized, but there is already evidence for their involvement in many of the pathophysiological mechanisms underlying the disease. In this review, we summarize the current knowledge of the functions and underlying mechanisms of lncRNA activity with a focus on type 2 diabetes mellitus and related renal and retinal complications of the disease. We also discuss the potential of lncRNAs to serve as therapeutic targets for drug development and diagnostic markers for clinical applications in the management of diabetes.

Cis- and trans-acting lncRNAs in pluripotency and reprogramming

Pervasive transcription in mammalian genomes produces thousands of long noncoding RNA (lncRNA) transcripts. Although they have been implicated in diverse biological processes, the functional relevance of most lncRNAs remains unknown. Recent studies reveal the prevalence of lncRNA-mediated cis regulation on nearby transcription. In this review, we summarize cis- and trans-acting lncRNAs involved in stem cell pluripotency and reprogramming, highlighting the role of regulatory lncRNAs in providing an additional layer of complexity to the regulation of genes that govern cell fate during development.

Biological action in Read-Write genome evolution

Many of the most important evolutionary variations that generated phenotypic adaptations and originated novel taxa resulted from complex cellular activities affecting genome content and expression. These activities included (i) the symbiogenetic cell merger that produced the mitochondrion-bearing ancestor of all extant eukaryotes, (ii) symbiogenetic cell mergers that produced chloroplast-bearing ancestors of photosynthetic eukaryotes, and (iii) interspecific hybridizations and genome doublings that generated new species and adaptive radiations of higher plants and animals. Adaptive variations also involved horizontal DNA transfers and natural genetic engineering by mobile DNA elements to rewire regulatory networks, such as those essential to viviparous reproduction in mammals. In the most highly evolved multicellular organisms, biological complexity scales with 'non-coding' DNA content rather than with protein-coding capacity in the genome. Coincidentally, 'non-coding' RNAs rich in repetitive mobile DNA sequences function as key regulators of complex adaptive phenotypes, such as stem cell pluripotency. The intersections of cell fusion activities, horizontal DNA transfers and natural genetic engineering of Read-Write genomes provide a rich molecular and biological foundation for understanding how ecological disruptions can stimulate productive, often abrupt, evolutionary transformations.

Global and cell-type specific properties of lincRNAs with ribosome occupancy

Advances in transcriptomics have led to the discovery of a large number of long intergenic non-coding RNAs (lincRNAs), which are now recognized as important regulators of diverse cellular processes. Although originally thought to be non-coding, recent studies have revealed that many lincRNAs are bound by ribosomes, with a few lincRNAs even having ability to generate micropeptides. The question arises: how widespread the translation of lincRNAs may be and whether such translation is likely to be functional. To better understand biological relevance of lincRNA translation, we systematically characterized lincRNAs with ribosome occupancy by the expression, structural, sequence, evolutionary and functional features for eight human cell lines, revealed that lincRNAs with ribosome occupancy have remarkably distinctive properties compared with those without ribosome occupancy, indicating that translation has important biological implication in categorizing and annotating lincRNAs. Further analysis revealed lincRNAs exhibit remarkable cell-type specificity with differential translational repertoires and substantial discordance in functionality. Collectively, our analyses provide the first attempt to characterize global and cell-type specific properties of translation of lincRNAs in human cells, highlighting that translation of lincRNAs has clear molecular, evolutionary and functional implications. This study will facilitate better understanding of the diverse functions of lincRNAs.

Human Retrotransposon Insertion Polymorphisms Are Associated with Health and Disease via Gene Regulatory Phenotypes

The human genome hosts several active families of transposable elements (TEs), including the Alu, LINE-1, and SVA retrotransposons that are mobilized via reverse transcription of RNA intermediates. We evaluated how insertion polymorphisms generated by human retrotransposon activity may be related to common health and disease phenotypes that have been previously interrogated through genome-wide association studies (GWAS). To address this question, we performed a genome-wide screen for retrotransposon polymorphism disease associations that are linked to TE induced gene regulatory changes. Our screen first identified polymorphic retrotransposon insertions found in linkage disequilibrium (LD) with single nucleotide polymorphisms that were previously associated with common complex diseases by GWAS. We further narrowed this set of candidate disease associated retrotransposon polymorphisms by identifying insertions that are located within tissue-specific enhancer elements. We then performed expression quantitative trait loci analysis on the remaining set of candidates in order to identify polymorphic retrotransposon insertions that are associated with gene expression changes in B-cells of the human immune system. This progressive and stringent screen yielded a list of six retrotransposon insertions as the strongest candidates for TE polymorphisms that lead to disease via enhancer-mediated changes in gene regulation. For example, we found an SVA insertion within a cell-type specific enhancer located in the second intron of the B4GALT1 gene. B4GALT1 encodes a glycosyltransferase that functions in the glycosylation of the Immunoglobulin G (IgG) antibody in such a way as to convert its activity from pro- to anti-inflammatory. The disruption of the B4GALT1 enhancer by the SVA insertion is associated with down-regulation of the gene in B-cells, which would serve to keep the IgG molecule in a pro-inflammatory state. Consistent with this idea, the B4GALT1 enhancer SVA insertion is linked to a genomic region implicated by GWAS in both inflammatory conditions and autoimmune diseases, such as systemic lupus erythematosus and Crohn’s disease. We explore this example and the other cases uncovered by our genome-wide screen in an effort to illuminate how retrotransposon insertion polymorphisms can impact human health and disease by causing changes in gene expression.

New Insights into Early Human Development: Lessons for Stem Cell Derivation and Differentiation

Pathways underlying mouse embryonic development have always informed efforts to derive, maintain, and drive differentiation of human pluripotent stem cells. However, direct application of mouse embryology to the human system has not always been successful because of fundamental developmental differences between species. The naive pluripotent state of mouse embryonic stem cells (ESCs), in particular, has been difficult to capture in human ESCs, and appears to be transitory in the human embryo itself. Further studies of human and non-human primate embryo development are needed to untangle the complexities of pluripotency networks across mammalian species.

The transformative potential of an integrative approach to pregnancy

BACKGROUND

Complex traits typically involve diverse biological pathways and are shaped by numerous genetic and environmental factors. Pregnancy-associated traits and pathologies are further complicated by extensive communication across multiple tissues in two individuals, interactions between two genomes-maternal and fetal-that obscure causal variants and lead to genetic conflict, and rapid evolution of pregnancy-associated traits across mammals and in the human lineage. Given the multi-faceted complexity of human pregnancy, integrative approaches that synthesize diverse data types and analyses harbor tremendous promise to identify the genetic architecture and environmental influences underlying pregnancy-associated traits and pathologies.

METHODS

We review current research that addresses the extreme complexities of traits and pathologies associated with human pregnancy.

RESULTS

We find that successful efforts to address the many complexities of pregnancy-associated traits and pathologies often harness the power of many and diverse types of data, including genome-wide association studies, evolutionary analyses, multi-tissue transcriptomic profiles, and environmental conditions.

CONCLUSION

We propose that understanding of pregnancy and its pathologies will be accelerated by computational platforms that provide easy access to integrated data and analyses. By simplifying the integration of diverse data, such platforms will provide a comprehensive synthesis that transcends many of the inherent challenges present in studies of pregnancy.

Systematic identification and characterization of regulatory elements derived from human endogenous retroviruses

Human endogenous retroviruses (HERVs) and other long terminal repeat (LTR)-type retrotransposons (HERV/LTRs) have regulatory elements that possibly influence the transcription of host genes. We systematically identified and characterized these regulatory elements based on publicly available datasets of ChIP-Seq of 97 transcription factors (TFs) provided by ENCODE and Roadmap Epigenomics projects. We determined transcription factor-binding sites (TFBSs) using the ChIP-Seq datasets and identified TFBSs observed on HERV/LTR sequences (HERV-TFBSs). Overall, 794,972 HERV-TFBSs were identified. Subsequently, we identified "HERV/LTR-shared regulatory element (HSRE)," defined as a TF-binding motif in HERV-TFBSs, shared within a substantial fraction of a HERV/LTR type. HSREs could be an indication that the regulatory elements of HERV/LTRs are present before their insertions. We identified 2,201 HSREs, comprising specific associations of 354 HERV/LTRs and 84 TFs. Clustering analysis showed that HERV/LTRs can be grouped according to the TF binding patterns; HERV/LTR groups bounded to pluripotent TFs (e.g., SOX2, POU5F1, and NANOG), embryonic endoderm/mesendoderm TFs (e.g., GATA4/6, SOX17, and FOXA1/2), hematopoietic TFs (e.g., SPI1 (PU1), GATA1/2, and TAL1), and CTCF were identified. Regulatory elements of HERV/LTRs tended to locate nearby and/or interact three-dimensionally with the genes involved in immune responses, indicating that the regulatory elements play an important role in controlling the immune regulatory network. Further, we demonstrated subgroup-specific TF binding within LTR7, LTR5B, and LTR5_Hs, indicating that gains or losses of the regulatory elements occurred during genomic invasions of the HERV/LTRs. Finally, we constructed dbHERV-REs, an interactive database of HERV/LTR regulatory elements (http://herv-tfbs.com/). This study provides fundamental information in understanding the impact of HERV/LTRs on host transcription, and offers insights into the transcriptional modulation systems of HERV/LTRs and ancestral HERVs.

G Quadruplex in Plants: A Ubiquitous Regulatory Element and Its Biological Relevance

G quadruplexes (G4) are higher-order DNA and RNA secondary structures formed by G-rich sequences that are built around tetrads of hydrogen-bonded guanine bases. Potential G4 quadruplex sequences have been identified in G-rich eukaryotic non-telomeric and telomeric genomic regions. Upon function, G4 formation is known to involve in chromatin remodeling, gene regulation and has been associated with genomic instability, genetic diseases and cancer progression. The natural role and biological validation of G4 structures is starting to be explored, and is of particular interest for the therapeutic interventions for human diseases. However, the existence and physiological role of G4 DNA and G4 RNA in plants species have not been much investigated yet and therefore, is of great interest for the development of improved crop varieties for sustainable agriculture. In this context, several recent studies suggests that these highly diverse G4 structures in plants can be employed to regulate expression of genes involved in several pathophysiological conditions including stress response to biotic and abiotic stresses as well as DNA damage. In the current review, we summarize the recent findings regarding the emerging functional significance of G4 structures in plants and discuss their potential value in the development of improved crop varieties.

Emerging roles of macrosatellite repeats in genome organization and disease development

Abundant repetitive DNA sequences are an enigmatic part of the human genome. Despite increasing evidence on the functionality of DNA repeats, their biologic role is still elusive and under frequent debate. Macrosatellites are the largest of the tandem DNA repeats, located on one or multiple chromosomes. The contribution of macrosatellites to genome regulation and human health was demonstrated for the D4Z4 macrosatellite repeat array on chromosome 4q35. Reduced copy number of D4Z4 repeats is associated with local euchromatinization and the onset of facioscapulohumeral muscular dystrophy. Although the role other macrosatellite families may play remains rather obscure, their diverse functionalities within the genome are being gradually revealed. In this review, we will outline structural and functional features of coding and noncoding macrosatellite repeats, and highlight recent findings that bring these sequences into the spotlight of genome organization and disease development.

Comparative analysis of lincRNA in insect species

BACKGROUND

The ever increasing availability of genomes makes it possible to investigate and compare not only the genomic complements of genes and proteins, but also of RNAs. One class of RNAs, the long noncoding RNAs (lncRNAs) and, in particular, their subclass of long intergenic noncoding RNAs (lincRNAs) have recently gained much attention because of their roles in regulation of important biological processes such as immune response or cell differentiation and as possible evolutionary precursors for protein coding genes. lincRNAs seem to be poorly conserved at the sequence level but at least some lincRNAs have conserved structural elements and syntenic genomic positions. Previous studies showed that transposable elements are a main contribution to the evolution of lincRNAs in mammals. In contrast, plant lincRNA emergence and evolution has been linked with local duplication events. However, little is known about their evolutionary dynamics in general and in insect genomes in particular.

RESULTS

Here we compared lincRNAs between seven insect genomes and investigated possible evolutionary changes and functional roles. We find very low sequence conservation between different species and that similarities within a species are mostly due to their association with transposable elements (TE) and simple repeats. Furthermore, we find that TEs are less frequent in lincRNA exons than in their introns, indicating that TEs may have been removed by selection. When we analysed the predicted thermodynamic stabilities of lincRNAs we found that they are more stable than their randomized controls which might indicate some selection pressure to maintain certain structural elements. We list several of the most stable lincRNAs which could serve as prime candidates for future functional studies. We also discuss the possibility of de novo protein coding genes emerging from lincRNAs. This is because lincRNAs with high GC content and potentially with longer open reading frames (ORF) are candidate loci where de novo gene emergence might occur.

CONCLUSION

The processes responsible for the emergence and diversification of lincRNAs in insects remain unclear. Both duplication and transposable elements may be important for the creation of new lincRNAs in insects.

Regulatory Roles of Long Non-Coding RNAs in the Central Nervous System and Associated Neurodegenerative Diseases

Accumulating studies have revealed that the human genome encodes tens of thousands of long non-coding RNAs (lncRNAs), which participate in multiple biological networks modulating gene expression via transcriptional, post-transcriptional and epigenetic regulation. Strikingly, a large fraction of tissue-specific lncRNAs are expressed in the Central Nervous System (CNS) with precisely regulated temporal and spatial expression patterns. These brain-specific lncRNAs are also featured with the cell-type specificity, the highest signals of evolutionary conservation, and their preferential location adjacent to brain-expressed protein-coding genes. Mounting evidence has indicated dysregulation or mutations in lncRNA gene loci are associated with a variety of CNS-associated neurodegenerative disorders, such as Alzheimer's, Parkinson's, Huntington's diseases, Amyotrophic Lateral Sclerosis and others. However, how lncRNAs contribute to these disorders remains to be further explored and studied. In this review article, we systematically and comprehensively summarize the current studies of lncRNAs, demonstrate the specificity of lncRNAs expressed in the brain, their functions during neural development and expression profiles in major cell types of the CNS, highlight the regulatory mechanisms of several studied lncRNAs that may play essential roles in the pathophysiology of neurodegenerative diseases, and discuss the current challenges and future perspectives of lncRNA studies involved in neurodegenerative and other diseases.

Chromatin-enriched lncRNAs can act as cell-type specific activators of proximal gene transcription

We recently described a new class of long noncoding RNA defined by especially tight chromatin association, whose presence is strongly correlated with expression of nearby genes in HEK293 cells. Here we critically examine the generality and cis-enhancer mechanism of this class of chromatin enriched RNA (cheRNA). CheRNA are largely cell-type specific, and remain the most effective chromatin signature for predicting cis-gene transcription in all cell types examined. Targeted depletion of three cheRNAs decreases gene expression of their neighbors, indicating potential co-activator function. Single-molecule FISH of one cheRNA-distal target gene pair suggests spatial overlap consistent with a role in chromosome looping. In another example, the cheRNA HIDALGO stimulates the fetal hemoglobin HBG1 gene during erythroid differentiation by promoting contacts to a downstream enhancer. Our results suggest that many cheRNAs activate proximal, lineage-specific gene transcription.

Primate-specific Long Non-coding RNAs and MicroRNAs

Non-coding RNAs (ncRNAs) are critical regulators of gene expression in essentially all life forms. Long ncRNAs (lncRNAs) and microRNAs (miRNAs) are two important RNA classes possessing regulatory functions. Up to date, many primate-specific ncRNAs have been identified and investigated. Their expression specificity to primate lineage suggests primate-specific roles. It is thus critical to elucidate the biological significance of primate or even human-specific ncRNAs, and to develop potential ncRNA-based therapeutics. Here, we have summarized the studies regarding regulatory roles of some key primate-specific lncRNAs and miRNAs.

Transposon-driven transcription is a conserved feature of vertebrate spermatogenesis and transcript evolution

Spermatogenesis is associated with major and unique changes to chromosomes and chromatin. Here, we sought to understand the impact of these changes on spermatogenic transcriptomes. We show that long terminal repeats (LTRs) of specific mouse endogenous retroviruses (ERVs) drive the expression of many long non‐coding transcripts (lncRNA). This process occurs post‐mitotically predominantly in spermatocytes and round spermatids. We demonstrate that this transposon‐driven lncRNA expression is a conserved feature of vertebrate spermatogenesis. We propose that transposon promoters are a mechanism by which the genome can explore novel transcriptional substrates, increasing evolutionary plasticity and allowing for the genesis of novel coding and non‐coding genes. Accordingly, we show that a small fraction of these novel ERV‐driven transcripts encode short open reading frames that produce detectable peptides. Finally, we find that distinct ERV elements from the same subfamilies act as differentially activated promoters in a tissue‐specific context. In summary, we demonstrate that LTRs can act as tissue‐specific promoters and contribute to post‐mitotic spermatogenic transcriptome diversity.

Evolinc: A Tool for the Identification and Evolutionary Comparison of Long Intergenic Non-coding RNAs

Long intergenic non-coding RNAs (lincRNAs) are an abundant and functionally diverse class of eukaryotic transcripts. Reported lincRNA repertoires in mammals vary, but are commonly in the thousands to tens of thousands of transcripts, covering ~90% of the genome. In addition to elucidating function, there is particular interest in understanding the origin and evolution of lincRNAs. Aside from mammals, lincRNA populations have been sparsely sampled, precluding evolutionary analyses focused on their emergence and persistence. Here we present Evolinc, a two-module pipeline designed to facilitate lincRNA discovery and characterize aspects of lincRNA evolution. The first module (Evolinc-I) is a lincRNA identification workflow that also facilitates downstream differential expression analysis and genome browser visualization of identified lincRNAs. The second module (Evolinc-II) is a genomic and transcriptomic comparative analysis workflow that determines the phylogenetic depth to which a lincRNA locus is conserved within a user-defined group of related species. Here we validate lincRNA catalogs generated with Evolinc-I against previously annotated Arabidopsis and human lincRNA data. Evolinc-I recapitulated earlier findings and uncovered an additional 70 Arabidopsis and 43 human lincRNAs. We demonstrate the usefulness of Evolinc-II by examining the evolutionary histories of a public dataset of 5,361 Arabidopsis lincRNAs. We used Evolinc-II to winnow this dataset to 40 lincRNAs conserved across species in Brassicaceae. Finally, we show how Evolinc-II can be used to recover the evolutionary history of a known lincRNA, the human telomerase RNA (TERC). These latter analyses revealed unexpected duplication events as well as the loss and subsequent acquisition of a novel TERC locus in the lineage leading to mice and rats. The Evolinc pipeline is currently integrated in CyVerse's Discovery Environment and is free for use by researchers.

EnHERV: Enrichment analysis of specific human endogenous retrovirus patterns and their neighboring genes

Human endogenous retroviruses (HERVs) are flanked by long terminal repeats (LTRs), which contain the regulation part of the retrovirus. Remaining HERVs constitute 7% to 8% of the present day human genome, and most have been identified as solo LTRs. The HERV sequences have been associated with several molecular functions as well as certain diseases in human, but their roles in human diseases are yet to be established. We designed EnHERV to make accessible the identified endogenous retrovirus repetitive sequences from Repbase Update (a database of eukaryotic repetitive elements) that are present in the human genome. Defragmentation process was done to improve the RepeatMasker annotation output. The defragmented elements were used as core database in EnHERV. EnHERV is available at http://sysbio.chula.ac.th/enherv and can be searched using either gene lists of user interest or HERV characteristics. Besides the search function, EnHERV also provides an enrichment analysis function that allows users to perform enrichment analysis between selected HERV characteristics and user-input gene lists, especially genes with the expression profile of a certain disease. EnHERV will facilitate exploratory studies of specific HERV characteristics that control gene expression patterns related to various disease conditions. Here we analyzed 25 selected HERV groups/names from all four HERV superfamilies, using the sense and anti-sense directions of the HERV and gene expression profiles from 49 specific tissue and disease conditions. We found that intragenic HERVs were associated with down-regulated genes in most cancer conditions and in psoriatic skin tissues and associated with up-regulated genes in immune cells particularly from systemic lupus erythematosus (SLE) patients. EnHERV allowed the analysis of how different types of LTRs were differentially associated with specific gene expression profiles in particular disease conditions for further studies into their mechanisms and functions.

Long noncoding RNAs in the pathogenesis of diabetic kidney disease: implications for novel therapeutic strategies

Diabetic kidney disease is a progressive disorder that develops secondary to diabetes. Current strategies for the clinical management of the disease can delay its onset and prevent progression, yet a significant proportion of patients still develop renal failure. The need for more advanced pharmaceuticals is therefore critical for improved treatment strategies. Recent studies support a role for long noncoding RNAs (lncRNAs) in the pathogenesis of human disease. Here we review recent experimental results linking lncRNAs with diabetic kidney disease. A better understanding of the regulatory role that lncRNAs play in the development of diabetic kidney disease may lead to identification of novel targets for therapeutic intervention.

De novo RNA sequence assembly during in vivo inflammatory stress reveals hundreds of unannotated lincRNAs in human blood CD14+ monocytes and in adipose tissue

Long intergenic noncoding RNAs (lincRNAs) have emerged as key regulators of cellular functions and physiology. Yet functional lincRNAs often have low, context-specific and tissue-specific expression. We hypothesized that many human monocyte and adipose lincRNAs would be absent in current public annotations due to lincRNA tissue specificity, modest sequencing depth in public data, limitations of transcriptome assembly algorithms, and lack of dynamic physiological contexts. Deep RNA sequencing (RNA-Seq) was performed in peripheral blood CD14⁺ monocytes (monocytes; average ~247 million reads per sample) and adipose tissue (average ~378 million reads per sample) collected before and after human experimental endotoxemia, an in vivo inflammatory stress, to identify tissue-specific and clinically relevant lincRNAs. Using a stringent filtering pipeline, we identified 109 unannotated lincRNAs in monocytes and 270 unannotated lincRNAs in adipose. Most unannotated lincRNAs are not conserved in rodents and are tissue specific, while many have features of regulated expression and are enriched in transposable elements. Specific subsets have enhancer RNA characteristics or are expressed only during inflammatory stress. A subset of unannotated lincRNAs was validated and replicated for their presence and inflammatory induction in independent human samples and for their monocyte and adipocyte origins. Through interrogation of public genome-wide association data, we also found evidence of specific disease association for selective unannotated lincRNAs. Our findings highlight the critical need to perform deep RNA-Seq in a cell-, tissue-, and context-specific manner to annotate the full repertoire of human lincRNAs for a complete understanding of lincRNA roles in dynamic cell functions and in human disease.

Transposable elements (TEs) contribute to stress-related long intergenic noncoding RNAs in plants

Noncoding RNAs have been extensively described in plant and animal transcriptomes by using high-throughput sequencing technology. Of these noncoding RNAs, a growing number of long intergenic noncoding RNAs (lincRNAs) have been described in multicellular organisms, however the origins and functions of many lincRNAs remain to be explored. In many eukaryotic genomes, transposable elements (TEs) are widely distributed and often account for large fractions of plant and animal genomes yet the contribution of TEs to lincRNAs is largely unknown. By using strand-specific RNA-sequencing, we profiled the expression patterns of lincRNAs in Arabidopsis, rice and maize, and identified 47 611 and 398 TE-associated lincRNAs (TE-lincRNAs), respectively. TE-lincRNAs were more often derived from retrotransposons than DNA transposons and as retrotransposon copy number in both rice and maize genomes so did TE-lincRNAs. We validated the expression of these TE-lincRNAs by strand-specific RT-PCR and also demonstrated tissue-specific transcription and stress-induced TE-lincRNAs either after salt, abscisic acid (ABA) or cold treatments. For Arabidopsis TE-lincRNA11195, mutants had reduced sensitivity to ABA as demonstrated by longer roots and higher shoot biomass when compared to wild-type. Finally, by altering the chromatin state in the Arabidopsis chromatin remodelling mutant ddm1, unique lincRNAs including TE-lincRNAs were generated from the preceding untranscribed regions and interestingly inherited in a wild-type background in subsequent generations. Our findings not only demonstrate that TE-associated lincRNAs play important roles in plant abiotic stress responses but lincRNAs and TE-lincRNAs might act as an adaptive reservoir in eukaryotes.

Human population-specific gene expression and transcriptional network modification with polymorphic transposable elements

Transposable element (TE) derived sequences are known to contribute to the regulation of the human genome. The majority of known TE-derived regulatory sequences correspond to relatively ancient insertions, which are fixed across human populations. The extent to which human genetic variation caused by recent TE activity leads to regulatory polymorphisms among populations has yet to be thoroughly explored. In this study, we searched for associations between polymorphic TE (polyTE) loci and human gene expression levels using an expression quantitative trait loci (eQTL) approach. We compared locus-specific polyTE insertion genotypes to B cell gene expression levels among 445 individuals from 5 human populations. Numerous human polyTE loci correspond to both cis and trans eQTL, and their regulatory effects are directly related to cell type-specific function in the immune system. PolyTE loci are associated with differences in expression between European and African population groups, and a single polyTE loci is indirectly associated with the expression of numerous genes via the regulation of the B cell-specific transcription factor PAX5. The polyTE-gene expression associations we found indicate that human TE genetic variation can have important phenotypic consequences. Our results reveal that TE-eQTL are involved in population-specific gene regulation as well as transcriptional network modification.

An atlas of human long non-coding RNAs with accurate 5' ends

Long non-coding RNAs (lncRNAs) are largely heterogeneous and functionally uncharacterized. Here, using FANTOM5 cap analysis of gene expression (CAGE) data, we integrate multiple transcript collections to generate a comprehensive atlas of 27,919 human lncRNA genes with high-confidence 5' ends and expression profiles across 1,829 samples from the major human primary cell types and tissues. Genomic and epigenomic classification of these lncRNAs reveals that most intergenic lncRNAs originate from enhancers rather than from promoters. Incorporating genetic and expression data, we show that lncRNAs overlapping trait-associated single nucleotide polymorphisms are specifically expressed in cell types relevant to the traits, implicating these lncRNAs in multiple diseases. We further demonstrate that lncRNAs overlapping expression quantitative trait loci (eQTL)-associated single nucleotide polymorphisms of messenger RNAs are co-expressed with the corresponding messenger RNAs, suggesting their potential roles in transcriptional regulation. Combining these findings with conservation data, we identify 19,175 potentially functional lncRNAs in the human genome.

Conserved expression of transposon-derived non-coding transcripts in primate stem cells

BACKGROUND

A significant portion of expressed non-coding RNAs in human cells is derived from transposable elements (TEs). Moreover, it has been shown that various long non-coding RNAs (lncRNAs), which come from the human endogenous retrovirus subfamily H (HERVH), are not only expressed but required for pluripotency in human embryonic stem cells (hESCs).

RESULTS

To identify additional TE-derived functional non-coding transcripts, we generated RNA-seq data from induced pluripotent stem cells (iPSCs) of four primate species (human, chimpanzee, gorilla, and rhesus) and searched for transcripts whose expression was conserved. We observed that about 30% of TE instances expressed in human iPSCs had orthologous TE instances that were also expressed in chimpanzee and gorilla. Notably, our analysis revealed a number of repeat families with highly conserved expression profiles including HERVH but also MER53, which is known to be the source of a placental-specific family of microRNAs (miRNAs). We also identified a number of repeat families from all classes of TEs, including MLT1-type and Tigger families, that contributed a significant amount of sequence to primate lncRNAs whose expression was conserved.

CONCLUSIONS

Together, these results describe TE families and TE-derived lncRNAs whose conserved expression patterns can be used to identify what are likely functional TE-derived non-coding transcripts in primate iPSCs.

Identification and evolutionary analysis of long non-coding RNAs in zebra finch

Background

Long non-coding RNAs (lncRNAs) are important in various biological processes, but very few studies on lncRNA have been conducted in birds. To identify IncRNAs expressed during feather development, we analyzed single-stranded RNA-seq (ssRNA-seq) data from the anterior and posterior dorsal regions during zebra finch (Taeniopygia guttata) embryonic development. Using published transcriptomic data, we further analyzed the evolutionary conservation of IncRNAs in birds and amniotes.

Results

A total of 1,081 lncRNAs, including 965 intergenic lncRNAs (lincRNAs), 59 intronic lncRNAs, and 57 antisense lncRNAs (lncNATs), were identified using our newly developed pipeline. These avian IncRNAs share similar characteristics with lncRNAs in mammals, such as shorter transcript length, lower exon number, lower average expression level and less sequence conservation than mRNAs. However, the proportion of lncRNAs overlapping with transposable elements in birds is much lower than that in mammals. We predicted the functions of IncRNAs based on the enriched functions of co-expressed protein-coding genes. Clusters of lncRNAs associated with natal down development were identified. The sequences and expression levels of candidate lncRNAs that shared conserved sequences among birds were validated by qPCR in both zebra finch and chicken. Finally, we identified three highly conserved lncRNAs that may be associated with natal down development.

Conclusions

Our study provides the first systematical identification of avian lncRNAs using ssRNA-seq analysis and offers a resource of embryonically expressed lncRNAs in zebra finch. We also predicted the biological function of identified lncRNAs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3506-z) contains supplementary material, which is available to authorized users.

Population and clinical genetics of human transposable elements in the (post) genomic era

Recent technological developments-in genomics, bioinformatics and high-throughput experimental techniques-are providing opportunities to study ongoing human transposable element (TE) activity at an unprecedented level of detail. It is now possible to characterize genome-wide collections of TE insertion sites for multiple human individuals, within and between populations, and for a variety of tissue types. Comparison of TE insertion site profiles between individuals captures the germline activity of TEs and reveals insertion site variants that segregate as polymorphisms among human populations, whereas comparison among tissue types ascertains somatic TE activity that generates cellular heterogeneity. In this review, we provide an overview of these new technologies and explore their implications for population and clinical genetic studies of human TEs. We cover both recent published results on human TE insertion activity as well as the prospects for future TE studies related to human evolution and health.

Immunobiology of Long Noncoding RNAs

The discovery of long noncoding RNAs (lncRNA) has provided a new perspective on gene regulation in diverse biological contexts. lncRNAs are remarkably versatile molecules that interact with RNA, DNA, or proteins to promote or restrain the expression of protein-coding genes. Activation of immune cells is associated with dynamic changes in expression of genes, the products of which combat infectious microorganisms, initiate repair, and resolve inflammatory responses in cells and tissues. Recent evidence indicates that lncRNAs play important roles in directing the development of diverse immune cells and controlling the dynamic transcriptional programs that are a hallmark of immune cell activation. The importance of these molecules is underscored by their newly recognized roles in inflammatory diseases. In this review, we discuss the contribution of lncRNAs in the development and activation of immune cells and their roles in immune-related diseases. We also discuss challenges faced in identifying biological functions for this large and complex class of genes.

Structural and Functional Annotation of Long Noncoding RNAs

Protein-coding RNAs represent only a small fraction of the transcriptional output in higher eukaryotes. The remaining RNA species encompass a broad range of molecular functions and regulatory roles, a consequence of the structural polyvalence of RNA polymers. Albeit several classes of small noncoding RNAs are relatively well characterized, the accessibility of affordable high-throughput sequencing is generating a wealth of novel, unannotated transcripts, especially long noncoding RNAs (lncRNAs) that are derived from genomic regions that are antisense, intronic, intergenic, and overlapping protein-coding loci. Parsing and characterizing the functions of noncoding RNAs-lncRNAs in particular-is one of the great challenges of modern genome biology. Here we discuss concepts and computational methods for the identification of structural domains in lncRNAs from genomic and transcriptomic data. In the first part, we briefly review how to identify RNA structural motifs in individual lncRNAs. In the second part, we describe how to leverage the evolutionary dynamics of structured RNAs in a computationally efficient screen to detect putative functional lncRNA motifs using comparative genomics.

History, Discovery, and Classification of lncRNAs

The RNA World Hypothesis suggests that prebiotic life revolved around RNA instead of DNA and proteins. Although modern cells have changed significantly in 4 billion years, RNA has maintained its central role in cell biology. Since the discovery of DNA at the end of the nineteenth century, RNA has been extensively studied. Many discoveries such as housekeeping RNAs (rRNA, tRNA, etc.) supported the messenger RNA model that is the pillar of the central dogma of molecular biology, which was first devised in the late 1950s. Thirty years later, the first regulatory non-coding RNAs (ncRNAs) were initially identified in bacteria and then in most eukaryotic organisms. A few long ncRNAs (lncRNAs) such as H19 and Xist were characterized in the pre-genomic era but remained exceptions until the early 2000s. Indeed, when the sequence of the human genome was published in 2001, studies showed that only about 1.2% encodes proteins, the rest being deemed "non-coding." It was later shown that the genome is pervasively transcribed into many ncRNAs, but their functionality remained controversial. Since then, regulatory lncRNAs have been characterized in many species and were shown to be involved in processes such as development and pathologies, revealing a new layer of regulation in eukaryotic cells. This newly found focus on lncRNAs, together with the advent of high-throughput sequencing, was accompanied by the rapid discovery of many novel transcripts which were further characterized and classified according to specific transcript traits.In this review, we will discuss the many discoveries that led to the study of lncRNAs, from Friedrich Miescher's "nuclein" in 1869 to the elucidation of the human genome and transcriptome in the early 2000s. We will then focus on the biological relevance during lncRNA evolution and describe their basic features as genes and transcripts. Finally, we will present a non-exhaustive catalogue of lncRNA classes, thus illustrating the vast complexity of eukaryotic transcriptomes.

The Role of Long Noncoding RNAs in Plant Stress Tolerance

Plants must adapt to multiple biotic and abiotic stresses ; thus, sensing and responding to environmental signals is imperative for their survival. Moreover, understanding these responses is imperative for efforts to improve plant yield and consistency. Regulation of transcript levels is a key aspect of the plant response to environmental signals. Long noncoding RNAs (lncRNAs) have gained widespread attention in recent years with the advance of high-throughput sequencing technologies. As important biological regulators, lncRNAs have been implicated in a wide range of developmental processes and diseases in animals. However, knowledge of the role that lncRNAs play in plant stress tolerance remains limited. Here, we review recent studies on the identification, characteristics, classification, and biological functions of lncRNAs in response to various stresses, including bacterial pathogens, excess light, drought, salinity, hypoxia, extreme temperatures, and nitrogen/phosphate deficiency. We also discuss possible directions for future research.

Study of Transposable Elements and Their Genomic Impact

Transposable elements (TEs) have been considered traditionally as junk DNA, i.e., DNA sequences that despite representing a high proportion of genomes had no evident cellular functions. However, over the last decades, it has become undeniable that not only TE-derived DNA sequences have (and had) a fundamental role during genome evolution, but also TEs have important implications in the origin and evolution of many genomic disorders. This concise review provides a brief overview of the different types of TEs that can be found in genomes, as well as a list of techniques and methods used to study their impact and mobilization. Some of these techniques will be covered in detail in this Method Book.

Endogenous retroviral promoter exaptation in human cancer

Cancer arises from a series of genetic and epigenetic changes, which result in abnormal expression or mutational activation of oncogenes, as well as suppression/inactivation of tumor suppressor genes. Aberrant expression of coding genes or long non-coding RNAs (lncRNAs) with oncogenic properties can be caused by translocations, gene amplifications, point mutations or other less characterized mechanisms. One such mechanism is the inappropriate usage of normally dormant, tissue-restricted or cryptic enhancers or promoters that serve to drive oncogenic gene expression. Dispersed across the human genome, endogenous retroviruses (ERVs) provide an enormous reservoir of autonomous gene regulatory modules, some of which have been co-opted by the host during evolution to play important roles in normal regulation of genes and gene networks. This review focuses on the “dark side” of such ERV regulatory capacity. Specifically, we discuss a growing number of examples of normally dormant or epigenetically repressed ERVs that have been harnessed to drive oncogenes in human cancer, a process we term onco-exaptation, and we propose potential mechanisms that may underlie this phenomenon.

Regulatory activities of transposable elements: from conflicts to benefits

Transposable elements (TEs) are a prolific source of tightly regulated, biochemically active non-coding elements, such as transcription factor-binding sites and non-coding RNAs. Many recent studies reinvigorate the idea that these elements are pervasively co-opted for the regulation of host genes. We argue that the inherent genetic properties of TEs and the conflicting relationships with their hosts facilitate their recruitment for regulatory functions in diverse genomes. We review recent findings supporting the long-standing hypothesis that the waves of TE invasions endured by organisms for eons have catalysed the evolution of gene-regulatory networks. We also discuss the challenges of dissecting and interpreting the phenotypic effect of regulatory activities encoded by TEs in health and disease.

The impact of transposable elements on mammalian development

Despite often being classified as selfish or junk DNA, transposable elements (TEs) are a group of abundant genetic sequences that have a significant impact on mammalian development and genome regulation. In recent years, our understanding of how pre-existing TEs affect genome architecture, gene regulatory networks and protein function during mammalian embryogenesis has dramatically expanded. In addition, the mobilization of active TEs in selected cell types has been shown to generate genetic variation during development and in fully differentiated tissues. Importantly, the ongoing domestication and evolution of TEs appears to provide a rich source of regulatory elements, functional modules and genetic variation that fuels the evolution of mammalian developmental processes. Here, we review the functional impact that TEs exert on mammalian developmental processes and discuss how the somatic activity of TEs can influence gene regulatory networks.

LncRNAs: emerging players in gene regulation and disease pathogenesis

The advent of next-generation sequencing has demonstrated that eukaryotic genomes are extremely complex than what were previously thought. Recent studies revealed that in addition to protein-coding genes, nonprotein-coding genes have allocated a large fraction of the genome. Long noncoding RNA (lncRNA) genes are classified as nonprotein-coding genes, serving as a molecular signal, decoy, guide and scaffold. They were suggested to play important roles in chromatin states, epigenetic and posttranscriptional regulation of genes. Aberrant expression of lncRNAs and changes in their structure are associated with a wide spectrum of diseases ranging from different types of cancer and neurodegeneration to ?-thalassaemia. The purpose of this study was to summarize the current progress in understanding the genomic bases and origin of lncRNAs. Moreover, this study focusses on the diverse functions of lncRNAs in normal cells as well as various types of disease to illustrate the potential impacts of lncRNAs on diverse biological processes and their therapeutic significance.

Inverted repeat Alu elements in the human lincRNA-p21 adopt a conserved secondary structure that regulates RNA function

LincRNA-p21 is a long intergenic non-coding RNA (lincRNA) involved in the p53-mediated stress response. We sequenced the human lincRNA-p21 (hLincRNA-p21) and found that it has a single exon that includes inverted repeat Alu elements (IRAlus). Sense and antisense Alu elements fold independently of one another into a secondary structure that is conserved in lincRNA-p21 among primates. Moreover, the structures formed by IRAlus are involved in the localization of hLincRNA-p21 in the nucleus, where hLincRNA-p21 colocalizes with paraspeckles. Our results underscore the importance of IRAlus structures for the function of hLincRNA-p21 during the stress response.

A lincRNA connected to cell mortality and epigenetically-silenced in most common human cancers

Immortality is an essential characteristic of human carcinoma cells. We recently developed an efficient, reproducible method that immortalizes human mammary epithelial cells (HMEC) in the absence of gross genomic changes by targeting 2 critical senescence barriers. Consistent transcriptomic changes associated with immortality were identified using microarray analysis of isogenic normal finite pre-stasis, abnormal finite post-stasis, and immortal HMECs from 4 individuals. A total of 277 genes consistently changed in cells that transitioned from post-stasis to immortal. Gene ontology analysis of affected genes revealed biological processes significantly altered in the immortalization process. These immortalization-associated changes showed striking similarity to the gene expression changes seen in The Cancer Genome Atlas (TCGA) clinical breast cancer data. The most dramatic change in gene expression seen during the immortalization step was the downregulation of an unnamed, incompletely annotated transcript that we called MORT, for mortality, since its expression was closely associated with the mortal, finite lifespan phenotype. We show here that MORT (ZNF667-AS1) is expressed in all normal finite lifespan human cells examined to date and is lost in immortalized HMEC. MORT gene silencing at the mortal/immortal boundary was due to DNA hypermethylation of its CpG island promoter. This epigenetic silencing is also seen in human breast cancer cell lines and in a majority of human breast tumor tissues. The functional importance of DNA hypermethylation in MORT gene silencing is supported by the ability of 5-aza-2'-deoxycytidine to reactivate MORT expression. Analysis of TCGA data revealed deregulation of MORT expression due to DNA hypermethylation in 15 out of the 17 most common human cancers. The epigenetic silencing of MORT in a large majority of the common human cancers suggests a potential fundamental role in cellular immortalization during human carcinogenesis.

Developmental piRNA profiles of the invasive vector mosquito Aedes albopictus

Background

In eukaryotic organisms, Piwi-interacting RNAs (piRNAs) control the activities of mobile genetic elements and ensure genome maintenance. Recent evidence indicates that piRNAs are involved in multiple biological pathways, including transcriptional regulation of protein-coding genes, sex determination and even interactions between host and pathogens. Aedes albopictus is a major invasive species that transmits a number of viral diseases in humans. Ae. albopictus has the largest genome and the highest abundance of repetitive sequences when compared with members that belong to Culicidae with a published genome. Analysis of piRNA profiles will provide a developmental and evolutionary perspective on piRNAs in Ae. albopictus.

Methods

piRNAs were identified and characterized during the development of Ae. albopictus, and piRNA expression patterns in adult males and females as well as sugar-fed females and blood-fed females were compared.

Results

Our results reveal that, despite the large genome size of Ae. albopictus, the piRNA pool of Ae. albopictus (1.2 × 10⁷) is smaller than those of Aedes aegypti (1.7 × 10⁷) and Drosophila melanogaster (1.6 × 10⁷). In Ae. albopictus, piRNAs displayed the highest abundance at the embryo stage and the lowest abundance at the pupal stage. Approximately 50 % of the piRNAs mapped to intergenic regions with no known functions. Approximately 30 % of the piRNAs mapped to repetitive elements, and 77.69 % of these repeat-derived piRNAs mapped to Class I TEs; 45.42 % of the observed piRNA reads originated from piRNA clusters, and most of the top 10 highest expressed piRNA clusters and 100 highest expressed piRNAs from each stage displayed biased expression patterns across the developmental stages. All anti-sense-derived piRNAs displayed a preference for uridine at the 5′ end; however, the sense-derived piRNAs showed adenine bias at the tenth nucleotide position and a typical ping-pong signature, suggesting that the biogenesis of piRNAs was conserved throughout development. Our results also show that 962 piRNAs displayed sex-biased expression, and 522 piRNAs showed higher expression in the blood-fed females than in the sugar-fed females.

Conclusions

Our results suggest that piRNAs, aside from silencing transposable elements in Ae. albopictus, may have a role in other biological pathways.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1815-8) contains supplementary material, which is available to authorized users.

The Tetraodon nigroviridis reference transcriptome: developmental transition, length retention and microsynteny of long non-coding RNAs in a compact vertebrate genome

Pufferfish such as fugu and tetraodon carry the smallest genomes among all vertebrates and are ideal for studying genome evolution. However, comparative genomics using these species is hindered by the poor annotation of their genomes. We performed RNA sequencing during key stages of maternal to zygotic transition of Tetraodon nigroviridis and report its first developmental transcriptome. We assembled 61,033 transcripts (23,837 loci) representing 80% of the annotated gene models and 3816 novel coding transcripts from 2667 loci. We demonstrate the similarities of gene expression profiles between pufferfish and zebrafish during maternal to zygotic transition and annotated 1120 long non-coding RNAs (lncRNAs) many of which differentially expressed during development. The promoters for 60% of the assembled transcripts result validated by CAGE-seq. Despite the extreme compaction of the tetraodon genome and the dramatic loss of transposons, the length of lncRNA exons remain comparable to that of other vertebrates and a small set of lncRNAs appears enriched for transposable elements suggesting a selective pressure acting on lncRNAs length and composition. Finally, a set of lncRNAs are microsyntenic between teleost and vertebrates, which indicates potential regulatory interactions between lncRNAs and their flanking coding genes. Our work provides a fundamental molecular resource for vertebrate comparative genomics and embryogenesis studies.

BmncRNAdb: a comprehensive database of non-coding RNAs in the silkworm, Bombyx mori

BACKGROUND

Long non-coding RNAs (lncRNAs) may play critical roles in a wide range of developmental processes of higher organisms. Recently, lncRNAs have been widely identified across eukaryotes and many databases of lncRNAs have been developed for human, mouse, fruit fly, etc. However, there is rare information about them in the only completely domesticated insect, silkworm (Bombyx mori).

DESCRIPTION

In this study, we systematically scanned lncRNAs using the available silkworm RNA-seq data and public unigenes. Finally, we identified and collected 6281 lncRNAs in the silkworm. Besides, we also collected 1986 microRNAs (miRNAs) from previous studies. Then, we organized them into a comprehensive and web-based database, BmncRNAdb. This database offers a user-friendly interface for data browse and online analysis as well as the three online tools for users to predict the target genes of lncRNA or miRNA.

CONCLUSIONS

We have systematically identified and collected the silkworm lncRNAs and constructed a comprehensive database of the silkworm lncRNAs and miRNAs. This work gives a glimpse into lncRNAs of the silkworm and lays foundations for the ncRNAs study of the silkworm and other insects in the future. The BmncRNAdb is freely available at http://gene.cqu.edu.cn/BmncRNAdb/index.php .

Transposable elements and G-quadruplexes

A significant part of eukaryotic genomes is formed by transposable elements (TEs) containing not only genes but also regulatory sequences. Some of the regulatory sequences located within TEs can form secondary structures like hairpins or three-stranded (triplex DNA) and four-stranded (quadruplex DNA) conformations. This review focuses on recent evidence showing that G-quadruplex-forming sequences in particular are often present in specific parts of TEs in plants and humans. We discuss the potential role of these structures in the TE life cycle as well as the impact of G-quadruplexes on replication, transcription, translation, chromatin status, and recombination. The aim of this review is to emphasize that TEs may serve as vehicles for the genomic spread of G-quadruplexes. These non-canonical DNA structures and their conformational switches may constitute another regulatory system that, together with small and long non-coding RNA molecules and proteins, contribute to the complex cellular network resulting in the large diversity of eukaryotes.

Evolutionary impact of transposable elements on genomic diversity and lineage-specific innovation in vertebrates

Since their discovery, a growing body of evidence has emerged demonstrating that transposable elements are important drivers of species diversity. These mobile elements exhibit a great variety in structure, size and mechanisms of transposition, making them important putative actors in organism evolution. The vertebrates represent a highly diverse and successful lineage that has adapted to a wide range of different environments. These animals also possess a rich repertoire of transposable elements, with highly diverse content between lineages and even between species. Here, we review how transposable elements are driving genomic diversity and lineage-specific innovation within vertebrates. We discuss the large differences in TE content between different vertebrate groups and then go on to look at how they affect organisms at a variety of levels: from the structure of chromosomes to their involvement in the regulation of gene expression, as well as in the formation and evolution of non-coding RNAs and protein-coding genes. In the process of doing this, we highlight how transposable elements have been involved in the evolution of some of the key innovations observed within the vertebrate lineage, driving the group's diversity and success.