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

The persistent contributions of RNA to eukaryotic gen(om)e architecture and cellular function

Currently, the best scenario for earliest forms of life is based on RNA molecules as they have the proven ability to catalyze enzymatic reactions and harbor genetic information. Evolutionary principles valid today become apparent in such models already. Furthermore, many features of eukaryotic genome architecture might have their origins in an RNA or RNA/protein (RNP) world, including the onset of a further transition, when DNA replaced RNA as the genetic bookkeeper of the cell. Chromosome maintenance, splicing, and regulatory function via RNA may be deeply rooted in the RNA/RNP worlds. Mostly in eukaryotes, conversion from RNA to DNA is still ongoing, which greatly impacts the plasticity of extant genomes. Raw material for novel genes encoding protein or RNA, or parts of genes including regulatory elements that selection can act on, continues to enter the evolutionary lottery.

Retrotransposons in pluripotent cells: Impact and new roles in cellular plasticity

Transposable Elements are pieces of DNA able to mobilize from one location to another within genomes. Although they constitute more than 50% of the human genome, they have been classified as selfish DNA, with the only mission to spread within genomes and generate more copies of themselves that will ensure their presence over generations. Despite their remarkable prevalence, only a minor group of transposable elements remain active in the human genome and can sporadically be associated with the generation of a genetic disorder due to their ongoing mobility. Most of the transposable elements identified in the human genome corresponded to fixed insertions that no longer move in genomes. As selfish DNA, transposable element insertions accumulate in cell types where genetic information can be passed to the next generation. Indeed, work from different laboratories has demonstrated that the main heritable load of TE accumulation in humans occurs during early embryogenesis. Thus, active transposable elements have a clear impact on our pluripotent genome. However, recent findings suggest that the main proportion of fixed non-mobile transposable elements might also have emerging roles in cellular plasticity. In this concise review, we provide an overview of the impact of currently active transposable elements in our pluripotent genome and further discuss new roles of transposable elements (active or not) in regulating pluripotency. This article is part of a Special Issue entitled: Stress as a fundamental theme in cell plasticity.

Identification and characterization of long intergenic non-coding RNAs related to mouse liver development

Long non-coding RNAs (lncRNAs) have been studied extensively over the last few years. Liver is an important organ that plays a crucial role in glucose metabolism and homeostasis; however, there are few reports of the identification and functional characterization of lncRNAs with important roles in liver development. Therefore, it is necessary to systematically identify lncRNAs that are involved in liver development. In this paper, we assembled the transcriptome using published RNA-seq data across three mouse liver developmental stages and identified 4,882 putative long intergenic non-coding RNAs (lincRNAs) expressed in at least one of the investigated stages. Combining these with Ensembl lincRNAs, we established a reference catalog of 6,602 transcribed lincRNAs in the mouse liver. We then analyzed all the lincRNAs in this reference catalog systematically and revealed that liver lincRNAs carry different genomic signatures from protein-coding genes, while the putative lincRNAs are generally comparable with known Ensembl lincRNAs. In addition, putative lincRNAs are functionally associated with essential biological processes, including RNA splicing, protein localization and fatty acid metabolic process, implying that they may play an important role in regulating liver development. The validation of selected lincRNAs that are specifically expressed in developing liver tissues further suggested the effectiveness of our approach. Our study shows that lincRNAs that are differentially expressed during three liver developmental stages could have important regulatory roles in liver development. The identified putative lincRNAs are a valuable resource for further functional studies.

A novel reannotation strategy for dissecting DNA methylation patterns of human long intergenic non-coding RNAs in cancers

Despite growing consensus that long intergenic non-coding ribonucleic acids (lincRNAs) are modulators of cancer, the knowledge about the deoxyribonucleic acid (DNA) methylation patterns of lincRNAs in cancers remains limited. In this study, we constructed DNA methylation profiles for 4629 tumors and 705 normal tissue samples from 20 different types of human cancer by reannotating data of DNA methylation arrays. We found that lincRNAs had different promoter methylation patterns in cancers. We classified 2461 lincRNAs into two categories and three subcategories, according to their promoter methylation patterns in tumors. LincRNAs with resistant methylation patterns in tumors had conserved transcriptional regulation regions and were ubiquitously expressed across normal tissues. By integrating cancer subtype data and patient clinical information, we identified lincRNAs with promoter methylation patterns that were associated with cancer status, subtype or prognosis for several cancers. Network analysis of aberrantly methylated lincRNAs in cancers showed that lincRNAs with aberrant methylation patterns might be involved in cancer development and progression. The methylated and demethylated lincRNAs identified in this study provide novel insights for developing cancer biomarkers and potential therapeutic targets.

A comprehensive analysis of piRNAs from adult human testis and their relationship with genes and mobile elements

Background

Piwi-interacting RNAs (piRNAs) are a recently discovered class of small non-coding RNAs whose best-understood function is to repress mobile element (ME) activity in animal germline. To date, nearly all piRNA studies have been conducted in model organisms and little is known about piRNA diversity, target specificity and biological function in human.

Results

Here we performed high-throughput sequencing of piRNAs from three human adult testis samples. We found that more than 81% of the ~17 million putative piRNAs mapped to ~6,000 piRNA-producing genomic clusters using a relaxed definition of clusters. A set of human protein-coding genes produces a relatively large amount of putative piRNAs from their 3’UTRs, and are significantly enriched for certain biological processes, suggestive of non-random sampling by the piRNA biogenesis machinery. Up to 16% of putative piRNAs mapped to a few hundred annotated long non-coding RNA (lncRNA) genes, suggesting that some lncRNA genes can act as piRNA precursors. Among major ME families, young families of LTR and endogenous retroviruses have a greater association with putative piRNAs than other MEs. In addition, piRNAs preferentially mapped to specific regions in the consensus sequences of several ME (sub)families and some piRNA mapping peaks showed patterns consistent with the “ping-pong” cycle of piRNA targeting and amplification.

Conclusions

Overall our data provide a comprehensive analysis and improved annotation of human piRNAs in adult human testes and shed new light into the relationship of piRNAs with protein-coding genes, lncRNAs, and mobile genetic elements in human.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-545) contains supplementary material, which is available to authorized users.

The RIDL hypothesis: transposable elements as functional domains of long noncoding RNAs

Our genome contains tens of thousands of long noncoding RNAs (lncRNAs), many of which are likely to have genetic regulatory functions. It has been proposed that lncRNA are organized into combinations of discrete functional domains, but the nature of these and their identification remain elusive. One class of sequence elements that is enriched in lncRNA is represented by transposable elements (TEs), repetitive mobile genetic sequences that have contributed widely to genome evolution through a process termed exaptation. Here, we link these two concepts by proposing that exonic TEs act as RNA domains that are essential for lncRNA function. We term such elements Repeat Insertion Domains of LncRNAs (RIDLs). A growing number of RIDLs have been experimentally defined, where TE-derived fragments of lncRNA act as RNA-, DNA-, and protein-binding domains. We propose that these reflect a more general phenomenon of exaptation during lncRNA evolution, where inserted TE sequences are repurposed as recognition sites for both protein and nucleic acids. We discuss a series of genomic screens that may be used in the future to systematically discover RIDLs. The RIDL hypothesis has the potential to explain how functional evolution can keep pace with the rapid gene evolution observed in lncRNA. More practically, TE maps may in the future be used to predict lncRNA function.

Co-option of Sox3 as the male-determining factor on the Y chromosome in the fish Oryzias dancena

Sex chromosomes harbour a primary sex-determining signal that triggers sexual development of the organism. However, diverse sex chromosome systems have been evolved in vertebrates. Here we use positional cloning to identify the sex-determining locus of a medaka-related fish, Oryzias dancena, and find that the locus on the Y chromosome contains a cis-regulatory element that upregulates neighbouring Sox3 expression in developing gonad. Sex-reversed phenotypes in Sox3(Y) transgenic fish, and Sox3(Y) loss-of-function mutants all point to its critical role in sex determination. Furthermore, we demonstrate that Sox3 initiates testicular differentiation by upregulating expression of downstream Gsdf, which is highly conserved in fish sex differentiation pathways. Our results not only provide strong evidence for the independent recruitment of Sox3 to male determination in distantly related vertebrates, but also provide direct evidence that a novel sex determination pathway has evolved through co-option of a transcriptional regulator potentially interacted with a conserved downstream component.

A novel RNA motif mediates the strict nuclear localization of a long noncoding RNA

The ubiquitous presence of long noncoding RNAs (lncRNAs) in eukaryotes points to the importance of understanding how their sequences impact function. As many lncRNAs regulate nuclear events and thus must localize to nuclei, we analyzed the sequence requirements for nuclear localization in an intergenic lncRNA named BORG (BMP2-OP1-responsive gene), which is both spliced and polyadenylated but is strictly localized in nuclei. Subcellular localization of BORG was not dependent on the context or level of its expression or decay but rather depended on the sequence of the mature, spliced transcript. Mutational analyses indicated that nuclear localization of BORG was mediated through a novel RNA motif consisting of the pentamer sequence AGCCC with sequence restrictions at positions -8 (T or A) and -3 (G or C) relative to the first nucleotide of the pentamer. Mutation of the motif to a scrambled sequence resulted in complete loss of nuclear localization, while addition of even a single copy of the motif to a cytoplasmically localized RNA was sufficient to impart nuclear localization. Further, the presence of this motif in other cellular RNAs showed a direct correlation with nuclear localization, suggesting that the motif may act as a general nuclear localization signal for cellular RNAs.

Intergenic lncRNAs and the evolution of gene expression

Eukaryote genomes encode a surprisingly large number of noncoding transcripts. Around two-thirds of human transcribed loci do not encode protein, and many are intergenic and produce long (>200 nucleotides) noncoding RNAs (lncRNAs). Extensive analyses using comparative genomics and transcriptomics approaches have established that lncRNA sequence and transcription tend to turn over rapidly during evolution. Our appreciation of the biological roles of lncRNAs, based only on a handful of transcripts with well-characterized functions, is that lncRNAs have diverse roles in regulating gene expression. These proposed roles together with their rapid rates of evolution suggest that lncRNAs could contribute to the divergent expression patterns observed among species and potentially to the origin of new traits.

Blood from 'junk': the LTR chimeric transcript Pu.2 promotes erythropoiesis

Transposable elements (TEs) are a prominent feature of most eukaryotic genomes. Despite rapidly accumulating evidence for the role of TE-driven insertional mutagenesis and structural variation in genome evolution, few clear examples of individual TEs impacting biology via perturbed gene regulation are available. A recent report describes the discovery of an alternative promoter for the murine erythroid transcription factor Pu.1. This promoter is located in an ORR1A0 long terminal repeat (LTR) retrotransposon intronic to Pu.1 and is regulated by the Krüppel-like factors KLF1 and KLF3. Expression of the resultant chimeric transcript, called Pu.2, spontaneously induces erythroid differentiation in vitro. These experiments illustrate how transcription factor binding sites spread by retrotransposition have the potential to impact networks encoding key biological processes in the host genome.

The rise of regulatory RNA

Discoveries over the past decade portend a paradigm shift in molecular biology. Evidence suggests that RNA is not only functional as a messenger between DNA and protein but also involved in the regulation of genome organization and gene expression, which is increasingly elaborate in complex organisms. Regulatory RNA seems to operate at many levels; in particular, it plays an important part in the epigenetic processes that control differentiation and development. These discoveries suggest a central role for RNA in human evolution and ontogeny. Here, we review the emergence of the previously unsuspected world of regulatory RNA from a historical perspective.

Epigenetic control of mobile DNA as an interface between experience and genome change

Mobile DNA in the genome is subject to RNA-targeted epigenetic control. This control regulates the activity of transposons, retrotransposons and genomic proviruses. Many different life history experiences alter the activities of mobile DNA and the expression of genetic loci regulated by nearby insertions. The same experiences induce alterations in epigenetic formatting and lead to trans-generational modifications of genome expression and stability. These observations lead to the hypothesis that epigenetic formatting directed by non-coding RNA provides a molecular interface between life history events and genome alteration.

Long noncoding RNAs, emerging players in muscle differentiation and disease

The vast majority of the mammalian genome is transcribed giving rise to many different types of noncoding RNAs. Among them, long noncoding RNAs are the most numerous and functionally versatile class. Indeed, the lncRNA repertoire might be as rich as the proteome. LncRNAs have emerged as key regulators of gene expression at multiple levels. They play important roles in the regulation of development, differentiation and maintenance of cell identity and they also contribute to disease. In this review, we present recent advances in the biology of lncRNAs in muscle development and differentiation. We will also discuss the contribution of lncRNAs to muscle disease with a particular focus on Duchenne and facioscapulohumeral muscular dystrophies.

The retrovirus HERVH is a long noncoding RNA required for human embryonic stem cell identity

Human endogenous retrovirus subfamily H (HERVH) is a class of transposable elements expressed preferentially in human embryonic stem cells (hESCs). Here, we report that the long terminal repeats of HERVH function as enhancers and that HERVH is a nuclear long noncoding RNA required to maintain hESC identity. Furthermore, HERVH is associated with OCT4, coactivators and Mediator subunits. Together, these results uncover a new role of species-specific transposable elements in hESCs.

Genetic variation in the non-coding genome: Involvement of micro-RNAs and long non-coding RNAs in disease

It has been found that the majority of disease-associated genetic variants identified by genome-wide association studies are located outside of protein-coding regions, where they seem to affect regions that control transcription (promoters, enhancers) and non-coding RNAs that also can influence gene expression. In this review, we focus on two classes of non-coding RNAs that are currently a major focus of interest: micro-RNAs and long non-coding RNAs. We describe their biogenesis, suggested mechanism of action, and discuss how these non-coding RNAs might be affected by disease-associated genetic alterations. The discovery of these alterations has already contributed to a better understanding of the etiopathology of human diseases and yielded insight into the function of these non-coding RNAs. We also provide an overview of available databases, bioinformatics tools, and high-throughput techniques that can be used to study the mechanism of action of individual non-coding RNAs. This article is part of a Special Issue entitled: From Genome to Function.

Stable C0T-1 repeat RNA is abundant and is associated with euchromatic interphase chromosomes

Recent studies recognize a vast diversity of noncoding RNAs with largely unknown functions, but few have examined interspersed repeat sequences, which constitute almost half our genome. RNA hybridization in situ using C0T-1 (highly repeated) DNA probes detects surprisingly abundant euchromatin-associated RNA comprised predominantly of repeat sequences (C0T-1 RNA), including LINE-1. C0T-1-hybridizing RNA strictly localizes to the interphase chromosome territory in cis and remains stably associated with the chromosome territory following prolonged transcriptional inhibition. The C0T-1 RNA territory resists mechanical disruption and fractionates with the nonchromatin scaffold but can be experimentally released. Loss of repeat-rich, stable nuclear RNAs from euchromatin corresponds to aberrant chromatin distribution and condensation. C0T-1 RNA has several properties similar to XIST chromosomal RNA but is excluded from chromatin condensed by XIST. These findings impact two "black boxes" of genome science: the poorly understood diversity of noncoding RNA and the unexplained abundance of repetitive elements.

Genome-wide discovery and characterization of maize long non-coding RNAs

BACKGROUND

Long non-coding RNAs (lncRNAs) are transcripts that are 200 bp or longer, do not encode proteins, and potentially play important roles in eukaryotic gene regulation. However, the number, characteristics and expression inheritance pattern of lncRNAs in maize are still largely unknown.

RESULTS

By exploiting available public EST databases, maize whole genome sequence annotation and RNA-seq datasets from 30 different experiments, we identified 20,163 putative lncRNAs. Of these lncRNAs, more than 90% are predicted to be the precursors of small RNAs, while 1,704 are considered to be high-confidence lncRNAs. High confidence lncRNAs have an average transcript length of 463 bp and genes encoding them contain fewer exons than annotated genes. By analyzing the expression pattern of these lncRNAs in 13 distinct tissues and 105 maize recombinant inbred lines, we show that more than 50% of the high confidence lncRNAs are expressed in a tissue-specific manner, a result that is supported by epigenetic marks. Intriguingly, the inheritance of lncRNA expression patterns in 105 recombinant inbred lines reveals apparent transgressive segregation, and maize lncRNAs are less affected by cis- than by trans-genetic factors.

CONCLUSIONS

We integrate all available transcriptomic datasets to identify a comprehensive set of maize lncRNAs, provide a unique annotation resource of the maize genome and a genome-wide characterization of maize lncRNAs, and explore the genetic control of their expression using expression quantitative trait locus mapping.

De novo assembly and transcriptome analysis of contrasting sugarcane varieties

Sugarcane is an important crop and a major source of sugar and alcohol. In this study, we performed de novo assembly and transcriptome annotation for six sugarcane genotypes involved in bi-parental crosses. The de novo assembly of the sugarcane transcriptome was performed using short reads generated using the Illumina RNA-Seq platform. We produced more than 400 million reads, which were assembled into 72,269 unigenes. Based on a similarity search, the unigenes showed significant similarity to more than 28,788 sorghum proteins, including a set of 5,272 unigenes that are not present in the public sugarcane EST databases; many of these unigenes are likely putative undescribed sugarcane genes. From this collection of unigenes, a large number of molecular markers were identified, including 5,106 simple sequence repeats (SSRs) and 708,125 single-nucleotide polymorphisms (SNPs). This new dataset will be a useful resource for future genetic and genomic studies in this species.

Regulation of Mouse Retroelement MuERV-L/MERVL Expression by REX1 and Epigenetic Control of Stem Cell Potency

About half of the mammalian genome is occupied by DNA sequences that originate from transposable elements. Retrotransposons can modulate gene expression in different ways and, particularly retrotransposon-derived long terminal repeats, profoundly shape expression of both surrounding and distant genomic loci. This is especially important in pre-implantation development, during which extensive reprograming of the genome takes place and cells pass through totipotent and pluripotent states. At this stage, the main mechanism responsible for retrotransposon silencing, i.e., DNA methylation, is inoperative. A particular retrotransposon called muERV-L/MERVL is expressed during pre-implantation stages and contributes to the plasticity of mouse embryonic stem cells. This review will focus on the role of MERVL-derived sequences as controlling elements of gene expression specific for pre-implantation development, two-cell stage-specific gene expression, and stem cell pluripotency, the epigenetic mechanisms that control their expression, and the contributions of the pluripotency marker REX1 and the related Yin Yang 1 family of transcription factors to this regulation process.

A-to-I RNA editing occurs at over a hundred million genomic sites, located in a majority of human genes

RNA molecules transmit the information encoded in the genome and generally reflect its content. Adenosine-to-inosine (A-to-I) RNA editing by ADAR proteins converts a genomically encoded adenosine into inosine. It is known that most RNA editing in human takes place in the primate-specific Alu sequences, but the extent of this phenomenon and its effect on transcriptome diversity are not yet clear. Here, we analyzed large-scale RNA-seq data and detected ∼1.6 million editing sites. As detection sensitivity increases with sequencing coverage, we performed ultradeep sequencing of selected Alu sequences and showed that the scope of editing is much larger than anticipated. We found that virtually all adenosines within Alu repeats that form double-stranded RNA undergo A-to-I editing, although most sites exhibit editing at only low levels (<1%). Moreover, using high coverage sequencing, we observed editing of transcripts resulting from residual antisense expression, doubling the number of edited sites in the human genome. Based on bioinformatic analyses and deep targeted sequencing, we estimate that there are over 100 million human Alu RNA editing sites, located in the majority of human genes. These findings set the stage for exploring how this primate-specific massive diversification of the transcriptome is utilized.

Inhibition of LINE-1 retrotransposon-encoded reverse transcriptase modulates the expression of cell differentiation genes in breast cancer cells

Long Interspersed Elements (L1 elements) are biologically active retrotransposons that are capable of autonomous replication using their own reverse transcriptase (RT) enzyme. Expression of the normally repressed RT has been implicated in cancer cell growth. However, at present, little is known about the expression of L1-encoded RT activity or the molecular changes that are associated with RT activity in the development of breast cancer. Here, we report that RT activity is widespread in breast cancer cells. The expression of RT protein decreased markedly in breast cancer cells after treatment with the antiretroviral drug, efavirenz. While the majority of cells showed a significant reduction in proliferation, inhibition of RT was also accompanied by cell-specific differences in morphology. MCF7 cells displayed elongated microtubule extensions that adhered tightly to their substrate, while a large fraction of the T47D cells that we studied formed long filopodia projections. These morphological changes were reversible upon cessation of RT inhibition, confirming their dependence on RT activity. We also carried out gene expression profiling with microarrays and determined the genes that were differentially expressed during the process of cellular differentiation. Genes involved in proliferation, cell migration, and invasive activity were repressed in RT-inhibited cells. Concomitantly, genes involved in cell projection, formation of vacuolar membranes, and cell-to-cell junctions were significantly upregulated in RT-inhibited cells. qRT-PCR examination of the mRNA expression of these genes in additional cell lines yielded close correlation between their differential expression and the degree of cellular differentiation. Our study demonstrates that the inhibition of L1-encoded RT can reduce the rate of proliferation and promote differentiation of breast cancer cells. Together, these results provide a direct functional link between the expression of L1 retrotransposons and the development of breast cancer.

Long non-coding RNAs: new players in cell differentiation and development

Genomes of multicellular organisms are characterized by the pervasive expression of different types of non-coding RNAs (ncRNAs). Long ncRNAs (lncRNAs) belong to a novel heterogeneous class of ncRNAs that includes thousands of different species. lncRNAs have crucial roles in gene expression control during both developmental and differentiation processes, and the number of lncRNA species increases in genomes of developmentally complex organisms, which highlights the importance of RNA-based levels of control in the evolution of multicellular organisms. In this Review, we describe the function of lncRNAs in developmental processes, such as in dosage compensation, genomic imprinting, cell differentiation and organogenesis, with a particular emphasis on mammalian development.

Functional insights into the role of nuclear-retained long noncoding RNAs in gene expression control in mammalian cells

The mammalian genome harbors thousands of long noncoding RNA (lncRNA) genes. Recent studies have indicated the involvement of several of these lncRNAs in the regulation of gene expression. lncRNAs play crucial roles in various biological processes ranging from epigenetic gene regulation, transcriptional control,to post-transcriptional regulation. lncRNAs are localized in various subcellular compartments, and major proportion of these are retained in the cell nucleus and could be broadly classified as nuclear-retained lncRNAs (nrRNAs). Based on the identified functions,members of the nrRNAs execute diverse roles, including providing architectural support to the hierarchical subnuclear organization and influencing the recruitment of chromatin modifier factors to specific chromatin sites. In this review, we will summarize the recently described roles of mammalian nrRNAs in controlling gene expression by influencing chromatin organization, transcription,pre-mRNA processing, nuclear organization, and their involvement in disease.

Chromatin signatures at transcriptional start sites separate two equally populated yet distinct classes of intergenic long noncoding RNAs

BACKGROUND

Mammalian transcriptomes contain thousands of long noncoding RNAs (lncRNAs). Some lncRNAs originate from intragenic enhancers which, when active, behave as alternative promoters producing transcripts that are processed using the canonical signals of their host gene. We have followed up this observation by analyzing intergenic lncRNAs to determine the extent to which they might also originate from intergenic enhancers.

RESULTS

We integrated high-resolution maps of transcriptional initiation and transcription to annotate a conservative set of intergenic lncRNAs expressed in mouse erythroblasts. We subclassified intergenic lncRNAs according to chromatin status at transcriptional initiation regions, defined by relative levels of histone H3K4 mono- and trimethylation. These transcripts are almost evenly divided between those arising from enhancer-associated (elncRNA) or promoter-associated (plncRNA) elements. These two classes of 5' capped and polyadenylated RNA transcripts are indistinguishable with regard to their length, number of exons or transcriptional orientation relative to their closest neighboring gene. Nevertheless, elncRNAs are more tissue-restricted, less highly expressed and less well conserved during evolution. Of considerable interest, we found that expression of elncRNAs, but not plncRNAs, is associated with enhanced expression of neighboring protein-coding genes during erythropoiesis.

CONCLUSIONS

We have determined globally the sites of initiation of intergenic lncRNAs in erythroid cells, allowing us to distinguish two similarly abundant classes of transcripts. Different correlations between the levels of elncRNAs, plncRNAs and expression of neighboring genes suggest that functional lncRNAs from the two classes may play contrasting roles in regulating the transcript abundance of local or distal loci.

DUX4 binding to retroelements creates promoters that are active in FSHD muscle and testis

The human double-homeodomain retrogene DUX4 is expressed in the testis and epigenetically repressed in somatic tissues. Facioscapulohumeral muscular dystrophy (FSHD) is caused by mutations that decrease the epigenetic repression of DUX4 in somatic tissues and result in mis-expression of this transcription factor in skeletal muscle. DUX4 binds sites in the human genome that contain a double-homeobox sequence motif, including sites in unique regions of the genome as well as many sites in repetitive elements. Using ChIP-seq and RNA-seq on myoblasts transduced with DUX4 we show that DUX4 binds and activates transcription of mammalian apparent LTR-retrotransposons (MaLRs), endogenous retrovirus (ERVL and ERVK) elements, and pericentromeric satellite HSATII sequences. Some DUX4-activated MaLR and ERV elements create novel promoters for genes, long non-coding RNAs, and antisense transcripts. Many of these novel transcripts are expressed in FSHD muscle cells but not control cells, and thus might contribute to FSHD pathology. For example, HEY1, a repressor of myogenesis, is activated by DUX4 through a MaLR promoter. DUX4-bound motifs, including those in repetitive elements, show evolutionary conservation and some repeat-initiated transcripts are expressed in healthy testis, the normal expression site of DUX4, but more rarely in other somatic tissues. Testis expression patterns are known to have evolved rapidly in mammals, but the mechanisms behind this rapid change have not yet been identified: our results suggest that mobilization of MaLR and ERV elements during mammalian evolution altered germline gene expression patterns through transcriptional activation by DUX4. Our findings demonstrate a role for DUX4 and repetitive elements in mammalian germline evolution and in FSHD muscular dystrophy.

Transposable elements (TEs) are found in most genomes, and many TEs create extra copies of themselves in new genomic locations by a process called retrotransposition. TEs are often thought of as genomic parasites that must be suppressed, because retrotransposition can cause great harm to their host organism. However, during evolution, the functions encoded by TEs have sometimes been co-opted to the advantage of the host genome as novel genes or as gene regulatory regions. We studied a human transcription factor called DUX4 that is normally expressed in testis and repressed in muscle. Sometimes muscle repression fails, causing the disease facioscapulohumeral muscular dystrophy (FSHD). We find that DUX4 binds many TE types and can activate their transcription. Some activated TEs have been co-opted as novel promoters for human genes. DUX4's activation of these genes via TEs might be important in the biology of normal testis and may contribute to the FSHD disease process. Our findings raise the possibility that DUX4 and TEs co-evolved, as TEs may have hijacked DUX4 to aid their retrotransposition while DUX4 may have utilized TEs to modify its transcriptional network in the evolving germline.

ANRIL/CDKN2B-AS shows two-stage clade-specific evolution and becomes conserved after transposon insertions in simians

BACKGROUND

Many long non-coding RNA (lncRNA) genes identified in mammals have multiple exons and functional domains, allowing them to bind to polycomb proteins, DNA methyltransferases, and specific DNA sequences to regulate genome methylation. Little is known about the origin and evolution of lncRNAs. ANRIL/CDKN2B-AS consists of 19 exons on human chromosome 9p21 and regulates the expression of three cyclin-dependent kinase inhibitors (CDKN2A/ARF/CDKN2B).

RESULTS

ANRIL/CDKN2B-AS originated in placental mammals, obtained additional exons during mammalian evolution but gradually lost them during rodent evolution, and reached 19 exons only in simians. ANRIL lacks splicing signals in mammals. In simians, multiple transposons were inserted and transformed into exons of the ANRIL gene, after which ANRIL became highly conserved. A further survey reveals that multiple transposons exist in many lncRNAs.

CONCLUSIONS

ANRIL shows a two-stage, clade-specific evolutionary process and is fully developed only in simians. The domestication of multiple transposons indicates an impressive pattern of "evolutionary tinkering" and is likely to be important for ANRIL's structure and function. The evolution of lncRNAs and that of transposons may be highly co-opted in primates. Many lncRNAs may be functional only in simians.

Long non-coding RNA identification over mouse brain development by integrative modeling of chromatin and genomic features

In silico prediction of genomic long non-coding RNAs (lncRNAs) is prerequisite to the construction and elucidation of non-coding regulatory network. Chromatin modifications marked by chromatin regulators are important epigenetic features, which can be captured by prevailing high-throughput approaches such as ChIP sequencing. We demonstrate that the accuracy of lncRNA predictions can be greatly improved when incorporating high-throughput chromatin modifications over mouse embryonic stem differentiation toward adult Cerebellum by logistic regression with LASSO regularization. The discriminating features include H3K9me3, H3K27ac, H3K4me1, open reading frames and several repeat elements. Importantly, chromatin information is suggested to be complementary to genomic sequence information, highlighting the importance of an integrated model. Applying integrated model, we obtain a list of putative lncRNAs based on uncharacterized fragments from transcriptome assembly. We demonstrate that the putative lncRNAs have regulatory roles in vicinity of known gene loci by expression and Gene Ontology enrichment analysis. We also show that the lncRNA expression specificity can be efficiently modeled by the chromatin data with same developmental stage. The study not only supports the biological hypothesis that chromatin can regulate expression of tissue-specific or developmental stage-specific lncRNAs but also reveals the discriminating features between lncRNA and coding genes, which would guide further lncRNA identifications and characterizations.

Hypomethylation marks enhancers within transposable elements

Transposable elements (TEs) make up 50% of the human genome and are usually considered a mutational burden. A new study uses signatures of DNA hypomethylation to identify tissue-specific enhancers within TEs, providing fresh evidence that mobile DNA has a non-negligible role in genome regulation and evolution.