Exaptation of transposable elements into novel cis-regulatory elements: is the evidence always strong?

On causal roles and selected effects: our genome is mostly junk

The idea that much of our genome is irrelevant to fitness-is not the product of positive natural selection at the organismal level-remains viable. Claims to the contrary, and specifically that the notion of "junk DNA" should be abandoned, are based on conflating meanings of the word "function". Recent estimates suggest that perhaps 90% of our DNA, though biochemically active, does not contribute to fitness in any sequence-dependent way, and possibly in no way at all. Comparisons to vertebrates with much larger and smaller genomes (the lungfish and the pufferfish) strongly align with such a conclusion, as they have done for the last half-century.

Regulatory remodeling in the allo-tetraploid frog Xenopus laevis

BACKGROUND

Genome duplication has played a pivotal role in the evolution of many eukaryotic lineages, including the vertebrates. A relatively recent vertebrate genome duplication is that in Xenopus laevis, which resulted from the hybridization of two closely related species about 17 million years ago. However, little is known about the consequences of this duplication at the level of the genome, the epigenome, and gene expression.

RESULTS

The X. laevis genome consists of two subgenomes, referred to as L (long chromosomes) and S (short chromosomes), that originated from distinct diploid progenitors. Of the parental subgenomes, S chromosomes have degraded faster than L chromosomes from the point of genome duplication until the present day. Deletions appear to have the largest effect on pseudogene formation and loss of regulatory regions. Deleted regions are enriched for long DNA repeats and the flanking regions have high alignment scores, suggesting that non-allelic homologous recombination has played a significant role in the loss of DNA. To assess innovations in the X. laevis subgenomes we examined p300-bound enhancer peaks that are unique to one subgenome and absent from X. tropicalis. A large majority of new enhancers comprise transposable elements. Finally, to dissect early and late events following interspecific hybridization, we examined the epigenome and the enhancer landscape in X. tropicalis × X. laevis hybrid embryos. Strikingly, young X. tropicalis DNA transposons are derepressed and recruit p300 in hybrid embryos.

CONCLUSIONS

The results show that erosion of X. laevis genes and functional regulatory elements is associated with repeats and non-allelic homologous recombination and furthermore that young repeats have also contributed to the p300-bound regulatory landscape following hybridization and whole-genome duplication.

Insertion Polymorphisms of Mobile Genetic Elements in Sexual and Asexual Populations of Daphnia pulex

Transposable elements (TEs) constitute a substantial portion of many eukaryotic genomes, and can in principle contribute to evolutionary innovation as well as genomic deterioration. Daphnia pulex serves as a useful model for studying TE dynamics as a potential cause and/or consequence of asexuality. We analyzed insertion polymorphisms of TEs in 20 sexual and 20 asexual isolates of D. pulex across North American from their available whole-genome sequencing data. Our results show that the total fraction of the derived sequences of TEs is not substantially different between asexual and sexual D. pulex isolates. However, in general, sexual clones contain fewer fixed TE insertions but more total insertion polymorphisms than asexual clones, supporting the hypothesis that sexual reproduction facilitates the spread and elimination of TEs. We identified nine asexual-specific fixed TE insertions, eight long terminal repeat retrotransposons, and one DNA transposon. By comparison, no sexual-specific fixed TE insertions were observed in our analysis. Furthermore, except one TE insertion located on a contig from chromosome 7, the other eight asexual-specific insertion sites are located on contigs from chromosome 9 that is known to be associated with obligate asexuality in D. pulex. We found that all nine asexual-specific fixed TE insertions can also be detected in some Daphnia pulicaria isolates, indicating that a substantial number of TE insertions in asexual D. pulex have been directly inherited from D. pulicaria during the origin of obligate asexuals.

Human evolution: the non-coding revolution

What made us human? Gene expression changes clearly played a significant part in human evolution, but pinpointing the causal regulatory mutations is hard. Comparative genomics enabled the identification of human accelerated regions (HARs) and other human-specific genome sequences. The major challenge in the past decade has been to link diverged sequences to uniquely human biology. This review discusses approaches to this problem, progress made at the molecular level, and prospects for moving towards genetic causes for uniquely human biology.

Repetitive DNA in the Catfish Genome: rDNA, Microsatellites, and Tc1-Mariner Transposon Sequences in Imparfinis Species (Siluriformes, Heptapteridae)

Physical mapping of repetitive DNA families in the karyotypes of fish is important to understand the organization and evolution of different orders, families, genera, or species. Fish in the genus Imparfinis show diverse karyotypes with various diploid numbers and ribosomal DNA (rDNA) locations. Here we isolated and characterized Tc1-mariner nucleotide sequences from Imparfinis schubarti, and mapped their locations together with 18S rDNA, 5S rDNA, and microsatellite probes in Imparfinis borodini and I. schubarti chromosomes. The physical mapping of Tc1/Mariner on chromosomes revealed dispersed signals in heterochromatin blocks with small accumulations in the terminal and interstitial regions of I. borodini and I. schubarti. Tc1/Mariner was coincident with rDNA chromosomes sites in both species, suggesting that this transposable element may have participated in the dispersion and evolution of these sequences in the fish genome. Our analysis suggests that different transposons and microsatellites have accumulated in the I. borodini and I. schubarti genomes and that the distribution patterns of these elements may be related to karyotype evolution within Imparfinis.

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 terminal repeats power evolution of genes and gene expression programs in mammalian oocytes and zygotes

Retrotransposons are "copy-and-paste" insertional mutagens that substantially contribute to mammalian genome content. Retrotransposons often carry long terminal repeats (LTRs) for retrovirus-like reverse transcription and integration into the genome. We report an extraordinary impact of a group of LTRs from the mammalian endogenous retrovirus-related ERVL retrotransposon class on gene expression in the germline and beyond. In mouse, we identified more than 800 LTRs from ORR1, MT, MT2, and MLT families, which resemble mobile gene-remodeling platforms that supply promoters and first exons. The LTR-mediated gene remodeling also extends to hamster, human, and bovine oocytes. The LTRs function in a stage-specific manner during the oocyte-to-embryo transition by activating transcription, altering protein-coding sequences, producing noncoding RNAs, and even supporting evolution of new protein-coding genes. These functions result, for example, in recycling processed pseudogenes into mRNAs or lncRNAs with regulatory roles. The functional potential of the studied LTRs is even higher, because we show that dormant LTR promoter activity can rescue loss of an essential upstream promoter. We also report a novel protein-coding gene evolution-D6Ertd527e-in which an MT LTR provided a promoter and the 5' exon with a functional start codon while the bulk of the protein-coding sequence evolved through a CAG repeat expansion. Altogether, ERVL LTRs provide molecular mechanisms for stochastically scanning, rewiring, and recycling genetic information on an extraordinary scale. ERVL LTRs thus offer means for a comprehensive survey of the genome's expression potential, tightly intertwining with gene expression and evolution in the germline.

Enhancer evolution and the origins of morphological novelty

A central goal of evolutionary biology is to understand the genetic origin of morphological novelties-i.e. anatomical structures unique to a taxonomic group. Elaboration of morphology during development depends on networks of regulatory genes that activate patterned gene expression through transcriptional enhancer regions. We summarize recent case studies and genome-wide investigations that have uncovered diverse mechanisms though which new enhancers arise. We also discuss how these enhancer-originating mechanisms have clarified the history of genetic networks underlying diversification of genital structures in flies, limbs and neural crest in chordates, and plant leaves. These studies have identified enhancers that were pivotal for morphological divergence and highlighted how novel genetic networks shaping form emerged from pre-existing ones.

LTR-retrotransposons in plants: Engines of evolution

LTR retrotransposons are the most abundant group of transposable elements (TEs) in plants. These elements can fall inside or close to genes, and therefore influence their expression and evolution. This review aims to examine how LTR retrotransposons, especially Ty1-copia elements, mediate gene regulation and evolution. Various stimuli, including polyploidization and biotic and abiotic elicitors, result in the transcription and movement of these retrotransposons, and can facilitate adaptation. The presence of cis-regulatory motifs in the LTRs are central to their stress-mediated responses and are shared with host stress-responsive genes, showing a complex evolutionary history in which TEs provide new regulatory units to genes. The presence of retrotransposon remnants in genes that are necessary for normal gene function, demonstrates the importance of exaptation and co-option, and is also a consequence of the abundance of these elements in plant genomes. Furthermore, insertions of LTR retrotransposons in and around genes provide potential for alternative splicing, epigenetic control, transduction, duplication and recombination. These characteristics can become an active part of the evolution of gene families as in the case of resistance genes (R-genes). The character of TEs as exclusively selfish is now being re-evaluated. Since genome-wide reprogramming via TEs is a long evolutionary process, the changes we can examine are case-specific and their fitness advantage may not be evident until TE-derived motifs and domains have been completely co-opted and fixed. Nevertheless, the presence of LTR retrotransposons inside genes and as part of gene promoter regions is consistent with their roles as engines of plant genome evolution.

Independent evolution of genomic characters during major metazoan transitions

Metazoan evolution encompasses a vast evolutionary time scale spanning over 600 million years. Our ability to infer ancestral metazoan characters, both morphological and functional, is limited by our understanding of the nature and evolutionary dynamics of the underlying regulatory networks. Increasing coverage of metazoan genomes enables us to identify the evolutionary changes of the relevant genomic characters such as the loss or gain of coding sequences, gene duplications, micro- and macro-synteny, and non-coding element evolution in different lineages. In this review we describe recent advances in our understanding of ancestral metazoan coding and non-coding features, as deduced from genomic comparisons. Some genomic changes such as innovations in gene and linkage content occur at different rates across metazoan clades, suggesting some level of independence among genomic characters. While their contribution to biological innovation remains largely unclear, we review recent literature about certain genomic changes that do correlate with changes to specific developmental pathways and metazoan innovations. In particular, we discuss the origins of the recently described pharyngeal cluster which is conserved across deuterostome genomes, and highlight different genomic features that have contributed to the evolution of this group. We also assess our current capacity to infer ancestral metazoan states from gene models and comparative genomics tools and elaborate on the future directions of metazoan comparative genomics relevant to evo-devo studies.

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.

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.

CTRL+INSERT: retrotransposons and their contribution to regulation and innovation of the transcriptome

The human genome contains millions of fragments from retrotransposons-highly repetitive DNA sequences that were once able to "copy and paste" themselves to other regions in the genome. However, the majority of retrotransposons have lost this capacity through acquisition of mutations or through endogenous silencing mechanisms. Without this imminent threat of transposition, retrotransposons have the potential to act as a major source of genomic innovation. Indeed, large numbers of retrotransposons have been found to be active in specific contexts: as gene regulatory elements and promoters for protein-coding genes or long noncoding RNAs, among others. In this review, we summarise recent findings about retrotransposons, with implications in gene expression regulation, the expansion of gene isoform diversity and the generation of long noncoding RNAs. We highlight key examples that demonstrate their role in cellular identity and their versatility as markers of cell states, and we discuss how their dysregulation may contribute to the formation of and possibly therapeutic response in human cancers.

Evolutionary direction of processed pseudogenes

While some pseudogenes have been reported to play important roles in gene regulation, little is known about the possible relationship between pseudogene functions and evolutionary process of pseudogenes, or about the forces responsible for the pseudogene evolution. In this study, we characterized human processed pseudogenes in terms of evolutionary dynamics. Our results show that pseudogenes tend to evolve toward: lower GC content, strong dinucleotide bias, reduced abundance of transcription factor binding motifs and short palindromes, and decreased ability to form nucleosomes. We explored possible evolutionary forces that shaped the evolution pattern of pseudogenes, and concluded that mutations in pseudogenes are likely determined, at least partially, by neighbor-dependent mutational bias and recombination-associated selection.

Integrated small copy number variations and epigenome maps of disorders of sex development

Small copy number variations (CNVs) have typically not been analyzed or reported in clinical settings and hence have remained underrepresented in databases and the literature. Here, we focused our investigations on these small CNVs using chromosome microarray analysis (CMA) data previously obtained from patients with atypical characteristics or disorders of sex development (DSD). Using our customized CMA track targeting 334 genes involved in the development of urogenital and reproductive structures and a less stringent analysis filter, we uncovered small genes with recurrent and overlapping CNVs as small as 1 kb, and small regions of homozygosity (ROHs), imprinting and position effects. Detailed analysis of these high-resolution data revealed CNVs and ROHs involving structural and functional domains, repeat elements, active transcription sites and regulatory regions. Integration of these genomic data with DNA methylation, histone modification and predicted RNA expression profiles in normal testes and ovaries suggested spatiotemporal and tissue-specific gene regulation. This study emphasized a DSD-specific and gene-targeted CMA approach that uncovered previously unanalyzed or unreported small genes and CNVs, contributing to the growing resources on small CNVs and facilitating the narrowing of the genomic gap for identifying candidate genes or regions. This high-resolution analysis tool could improve the diagnostic utility of CMA, not only in patients with DSD but also in other clinical populations. These integrated data provided a better genomic-epigenomic landscape of DSD and greater opportunities for downstream research.

Long Terminal Repeats: From Parasitic Elements to Building Blocks of the Transcriptional Regulatory Repertoire

The life cycle of endogenous retroviruses (ERVs), also called long terminal repeat (LTR) retrotransposons, begins with transcription by RNA polymerase II followed by reverse transcription and re-integration into the host genome. While most ERVs are relics of ancient integration events, "young" proviruses competent for retrotransposition-found in many mammals, but not humans-represent an ongoing threat to host fitness. As a consequence, several restriction pathways have evolved to suppress their activity at both transcriptional and post-transcriptional stages of the viral life cycle. Nevertheless, accumulating evidence has revealed that LTR sequences derived from distantly related ERVs have been exapted as regulatory sequences for many host genes in a wide range of cell types throughout mammalian evolution. Here, we focus on emerging themes from recent studies cataloging the diversity of ERV LTRs acting as important transcriptional regulatory elements in mammals and explore the molecular features that likely account for LTR exaptation in developmental and tissue-specific gene regulation.

Enhancer activity of Helitron in sericin-1 gene promoter from Bombyx mori

Sericin is a kind of water-soluble protein expressed specifically in the middle silk gland of Bombyx mori. When the sericin-1 gene promoter was cloned and a transgenic vector was constructed to express a foreign protein, a specific Helitron, Bmhel-8, was identified in the sericin-1 gene promoter sequence in some genotypes of Bombyx mori and Bombyx mandarina. Given that the Bmhel-8 Helitron transposon was present only in some genotypes, it could be the source of allelic variation in the sericin-1 promoter. The length of the sericin-1 promoter sequence is approximately 1063 or 643 bp. The larger size of the sequence or allele is ascribed to the presence of Bmhel-8. Silkworm genotypes can be homozygous for either the shorter or larger promoter sequence or heterozygous, containing both alleles. Bmhel-8 in the sericin-1 promoter exhibits enhancer activity, as demonstrated by a dual-luciferase reporter system in BmE cell lines. Furthermore, Bmhel-8 displays enhancer activity in a sericin-1 promoter-driven gene expression system but does not regulate the tissue-specific expression of sericin-1.

Massive contribution of transposable elements to mammalian regulatory sequences

Barbara McClintock discovered the existence of transposable elements (TEs) in the late 1940s and initially proposed that they contributed to the gene regulatory program of higher organisms. This controversial idea gained acceptance only much later in the 1990s, when the first examples of TE-derived promoter sequences were uncovered. It is now known that half of the human genome is recognizably derived from TEs. It is thus important to understand the scope and nature of their contribution to gene regulation. Here, we provide a timeline of major discoveries in this area and discuss how transposons have revolutionized our understanding of mammalian genomes, with a special emphasis on the massive contribution of TEs to primate evolution. Our analysis of primate-specific functional elements supports a simple model for the rate at which new functional elements arise in unique and TE-derived DNA. Finally, we discuss some of the challenges and unresolved questions in the field, which need to be addressed in order to fully characterize the impact of TEs on gene regulation, evolution and disease processes.

Evolution of Epigenetic Regulation in Vertebrate Genomes

Empirical models of sequence evolution have spurred progress in the field of evolutionary genetics for decades. We are now realizing the importance and complexity of the eukaryotic epigenome. While epigenome analysis has been applied to genomes from single-cell eukaryotes to human, comparative analyses are still relatively few and computational algorithms to quantify epigenome evolution remain scarce. Accordingly, a quantitative model of epigenome evolution remains to be established. We review here the comparative epigenomics literature and synthesize its overarching themes. We also suggest one mechanism, transcription factor binding site (TFBS) turnover, which relates sequence evolution to epigenetic conservation or divergence. Lastly, we propose a framework for how the field can move forward to build a coherent quantitative model of epigenome evolution.

Mariner Transposons Contain a Silencer: Possible Role of the Polycomb Repressive Complex 2

Transposable elements are driving forces for establishing genetic innovations such as transcriptional regulatory networks in eukaryotic genomes. Here, we describe a silencer situated in the last 300 bp of the Mos1 transposase open reading frame (ORF) which functions in vertebrate and arthropod cells. Functional silencers are also found at similar locations within three other animal mariner elements, i.e. IS630-Tc1-mariner (ITm) DD34D elements, Himar1, Hsmar1 and Mcmar1. These silencers are able to impact eukaryotic promoters monitoring strong, moderate or low expression as well as those of mariner elements located upstream of the transposase ORF. We report that the silencing involves at least two transcription factors (TFs) that are conserved within animal species, NFAT-5 and Alx1. These cooperatively act with YY1 to trigger the silencing activity. Four other housekeeping transcription factors (TFs), neuron restrictive silencer factor (NRSF), GAGA factor (GAF) and GTGT factor (GTF), were also found to have binding sites within mariner silencers but their impact in modulating the silencer activity remains to be further specified. Interestingly, an NRSF binding site was found to overlap a 30 bp motif coding a highly conserved PHxxYSPDLAPxD peptide in mariner transposases. We also present experimental evidence that silencing is mainly achieved by co-opting the host Polycomb Repressive Complex 2 pathway. However, we observe that when PRC2 is impaired another host silencing pathway potentially takes over to maintain weak silencer activity. Mariner silencers harbour features of Polycomb Response Elements, which are probably a way for mariner elements to self-repress their transcription and mobility in somatic and germinal cells when the required TFs are expressed. At the evolutionary scale, mariner elements, through their exaptation, might have been a source of silencers playing a role in the chromatin configuration in eukaryotic genomes.

Transposons are mobile DNA sequences that have long co-evolved with the genome of their hosts. Consequently, they are involved in the generation of mutations, as well as the creation of genes and regulatory networks. Controlling the transposon activity, and consequently its negative effects on both the host soma and germ line, is a challenge for the survival of both the host and the transposon. To silence transposons, hosts often use defence mechanisms involving DNA methylation and RNA interference pathways. Here we show that mariner transposons can self-regulate their activity by using a silencer element located in their DNA sequence. The silencer element interferes with host housekeeping protein transcription factors involved in the polycomb silencing pathways. As the regulation of chromatin configuration by polycomb is an important regulator of animal development, our findings open the possibility that mariner silencers might have been exapted during animal evolution to participate in certain regulation pathways of their hosts. Since some of the TFs involved in mariner silencer activity play a role at different stages of nervous system development and neuron differentiation, it might be possible that mariner transposons can be active during some steps of cell differentiation. Interestingly, mariner transposons (i.e. IS630-Tc1-mariner (ITm) DD34D transposons) have so far only been found in genomes of animals having a nervous system.

The domain structure and distribution of Alu elements in long noncoding RNAs and mRNAs

Approximately 75% of the human genome is transcribed and many of these spliced transcripts contain primate-specific Alu elements, the most abundant mobile element in the human genome. The majority of exonized Alu elements are located in long noncoding RNAs (lncRNAs) and the untranslated regions of mRNA, with some performing molecular functions. To further assess the potential for Alu elements to be repurposed as functional RNA domains, we investigated the distribution and evolution of Alu elements in spliced transcripts. Our analysis revealed that Alu elements are underrepresented in mRNAs and lncRNAs, suggesting that most exonized Alu elements arising in the population are rare or deleterious to RNA function. When mRNAs and lncRNAs retain exonized Alu elements, they have a clear preference for Alu dimers, left monomers, and right monomers. mRNAs often acquire Alu elements when their genes are duplicated within Alu-rich regions. In lncRNAs, reverse-oriented Alu elements are significantly enriched and are not restricted to the 3' and 5' ends. Both lncRNAs and mRNAs primarily contain the Alu J and S subfamilies that were amplified relatively early in primate evolution. Alu J subfamilies are typically overrepresented in lncRNAs, whereas the Alu S dimer is overrepresented in mRNAs. The sequences of Alu dimers tend to be constrained in both lncRNAs and mRNAs, whereas the left and right monomers are constrained within particular Alu subfamilies and classes of RNA. Collectively, these findings suggest that Alu-containing RNAs are capable of forming stable structures and that some of these Alu domains might have novel biological functions.

The 3D folding of metazoan genomes correlates with the association of similar repetitive elements

The potential roles of the numerous repetitive elements found in the genomes of multi-cellular organisms remain speculative. Several studies have suggested a role in stabilizing specific 3D genomic contacts. To test this hypothesis, we exploited inter-chromosomal contacts frequencies obtained from Hi-C experiments and show that the folding of the human, mouse and Drosophila genomes is associated with a significant co-localization of several specific repetitive elements, notably many elements of the SINE family. These repeats tend to be the oldest ones and are enriched in transcription factor binding sites. We propose that the co-localization of these repetitive elements may explain the global conservation of genome folding observed between homologous regions of the human and mouse genome. Taken together, these results support a contribution of specific repetitive elements in maintaining and/or reshaping genome architecture over evolutionary times.

Evolution of the EKA family of powdery mildew avirulence-effector genes from the ORF 1 of a LINE retrotransposon

Background

The Avrk1 and Avra10 avirulence (AVR) genes encode effectors that increase the pathogenicity of the fungus Blumeria graminis f.sp. hordei (Bgh), the powdery mildew pathogen, in susceptible barley plants. In resistant barley, MLK1 and MLA10 resistance proteins recognize the presence of AVRK1 and AVRA10, eliciting the hypersensitive response typical of gene for gene interactions. Avrk1 and Avra10 have more than 1350 homologues in Bgh genome, forming the EKA (Effectors homologous to Avrk1 and Avra10) gene family.

Results

We tested the hypothesis that the EKA family originated from degenerate copies of Class I LINE retrotransposons by analysing the EKA family in the genome of Bgh isolate DH14 with bioinformatic tools specially developed for the analysis of Transposable Elements (TE) in genomes. The Class I LINE retrotransposon copies homologous to Avrk1 and Avra10 represent 6.5 % of the Bgh annotated genome and, among them, we identified 293 AVR/effector candidate genes. We also experimentally identified peptides that indicated the translation of several predicted proteins from EKA family members, which had higher relative abundance in haustoria than in hyphae.

Conclusions

Our analyses indicate that Avrk1 and Avra10 have evolved from part of the ORF1 gene of Class I LINE retrotransposons. The co-option of Avra10 and Avrk1 as effectors from truncated copies of retrotransposons explains the huge number of homologues in Bgh genome that could act as dynamic reservoirs from which new effector genes may evolve. These data provide further evidence for recruitment of retrotransposons in the evolution of new biological functions.

Electronic supplementary material

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

Evolution of DNA-Binding Sites of a Floral Master Regulatory Transcription Factor

Flower development is controlled by the action of key regulatory transcription factors of the MADS-domain family. The function of these factors appears to be highly conserved among species based on mutant phenotypes. However, the conservation of their downstream processes is much less well understood, mostly because the evolutionary turnover and variation of their DNA-binding sites (BSs) among plant species have not yet been experimentally determined. Here, we performed comparative ChIP (chromatin immunoprecipitation)-seq experiments of the MADS-domain transcription factor SEPALLATA3 (SEP3) in two closely related Arabidopsis species: Arabidopsis thaliana and A. lyrata which have very similar floral organ morphology. We found that BS conservation is associated with DNA sequence conservation, the presence of the CArG-box BS motif and on the relative position of the BS to its potential target gene. Differences in genome size and structure can explain that SEP3 BSs in A. lyrata can be located more distantly to their potential target genes than their counterparts in A. thaliana. In A. lyrata, we identified transposition as a mechanism to generate novel SEP3 binding locations in the genome. Comparative gene expression analysis shows that the loss/gain of BSs is associated with a change in gene expression. In summary, this study investigates the evolutionary dynamics of DNA BSs of a floral key-regulatory transcription factor and explores factors affecting this phenomenon.

Functions of noncoding sequences in mammalian genomes

Most of the mammalian genome consists of nucleotide sequences not coding for proteins. Exons of genes make up only 3% of the human genome, while the significance of most other sequences remains unknown. Recent genome studies with high-throughput methods demonstrate that the so-called noncoding part of the genome may perform important functions. This hypothesis is supported by three groups of experimental data: 1) approximately 10% of the sequences, most of which are located in noncoding parts of the genome, is evolutionarily conserved and thus can be of functional importance; 2) up to 99% of the mammalian genome is being transcribed forming short and long noncoding RNAs in addition to common mRNA; and 3) mutations in noncoding parts of the genome can be accompanied by progression of pathological states of the organism. In the light of these data, in the review we consider the functional role of numerous known sequences of noncoding parts of the genome including introns, DNA methylation regions, enhancers and locus control regions, insulators, S/MAR sequences, pseudogenes, and genes of noncoding RNAs, as well as transposons and simple repeats of centromeric and telomeric regions of chromosomes. The assumption is made that the intergenic noncoding sequences without definite/clear functions can be involved in spatial organization of genetic loci in interphase nuclei.

Genomic approaches to studying human-specific developmental traits

Changes in developmental regulatory programs drive both disease and phenotypic differences among species. Linking human-specific traits to alterations in development is challenging, because we have lacked the tools to assay and manipulate regulatory networks in human and primate embryonic cells. This field was transformed by the sequencing of hundreds of genomes – human and non-human – that can be compared to discover the regulatory machinery of genes involved in human development. This approach has identified thousands of human-specific genome alterations in developmental genes and their regulatory regions. With recent advances in stem cell techniques, genome engineering, and genomics, we can now test these sequences for effects on developmental gene regulation and downstream phenotypes in human cells and tissues.

Acytota - associated kingdom of neglected life

There is a huge variety of RNA- and DNA-containing entities that multiply within and propagate between cells across all kingdoms of life, having no cells of their own. Apart from cellular organisms, these entities (viroids, plasmids, mobile elements and viruses among others) are the only ones with distinct genetic identities but which are not included in any traditional tree of life. We suggest to introduce or, rather, revive the distinct category of acellular organisms, Acytota, as an additional, undeservedly ignored full-fledged kingdom of life. Acytota are indispensable players in cellular life and its evolution. The six traditional kingdoms (Cytota) and Acytota together complete the classification of the biological world (Biota), leaving nothing beyond.

Retrotransposition creates sloping shores: a graded influence of hypomethylated CpG islands on flanking CpG sites

Long interspersed elements (LINEs), through both self-mobilization and trans-mobilization of short interspersed elements and processed pseudogenes, have made an indelible impact on the structure and function of the human genome. One consequence is the creation of new CpG islands (CGIs). In fact, more than half of all CGIs in the genome are associated with repetitive DNA, three-quarters of which are derived from retrotransposons. However, little is known about the epigenetic impact of newly inserted CGIs. We utilized a transgenic LINE-1 mouse model and tracked DNA methylation dynamics of individual germline insertions during mouse development. The retrotransposed GFP marker sequence, a strong CGI, is hypomethylated in male germ cells but hypermethylated in somatic tissues, regardless of genomic location. The GFP marker is similarly methylated when delivered into the genome via the Sleeping Beauty DNA transposon, suggesting that the observed methylation pattern may be independent of the mode of insertion. Comparative analyses between insertion- and non-insertion-containing alleles further reveal a graded influence of the retrotransposed CGI on flanking CpG sites, a phenomenon that we described as "sloping shores." Computational analyses of human and mouse methylomic data at single-base resolution confirm that sloping shores are universal for hypomethylated CGIs in sperm and somatic tissues. Additionally, the slope of a hypomethylated CGI can be affected by closely positioned CGI neighbors. Finally, by tracing sloping shore dynamics through embryonic and germ cell reprogramming, we found evidence of bookmarking, a mechanism that likely determines which CGIs will be eventually hyper- or hypomethylated.

[ENCODE apophenia or a panglossian analysis of the human genome]

In September 2012, a batch of more than 30 articles presenting the results of the ENCODE (Encyclopaedia of DNA Elements) project was released. Many of these articles appeared in Nature and Science, the two most prestigious interdisciplinary scientific journals. Since that time, hundreds of other articles dedicated to the further analyses of the Encode data have been published. The time of hundreds of scientists and hundreds of millions of dollars were not invested in vain since this project had led to an apparent paradigm shift: contrary to the classical view, 80% of the human genome is not junk DNA, but is functional. This hypothesis has been criticized by evolutionary biologists, sometimes eagerly, and detailed refutations have been published in specialized journals with impact factors far below those that published the main contribution of the Encode project to our understanding of genome architecture. In 2014, the Encode consortium released a new batch of articles that neither suggested that 80% of the genome is functional nor commented on the disappearance of their 2012 scientific breakthrough. Unfortunately, by that time many biologists had accepted the idea that 80% of the genome is functional, or at least, that this idea is a valid alternative to the long held evolutionary genetic view that it is not. In order to understand the dynamics of the genome, it is necessary to re-examine the basics of evolutionary genetics because, not only are they well established, they also will allow us to avoid the pitfall of a panglossian interpretation of Encode. Actually, the architecture of the genome and its dynamics are the product of trade-offs between various evolutionary forces, and many structural features are not related to functional properties. In other words, evolution does not produce the best of all worlds, not even the best of all possible worlds, but only one possible world.

Putative regulation mechanism for the MSTN gene by a CpG island generated by the SINE marker Ins227bp

Background

A single nucleotide polymorphism (SNP) in the first intron of the myostatin gene (MSTN) is associated with aptness of elite Thoroughbreds to race over sprint, middle or long distances. This intronic marker (g.66493737 T ≻ C), a short interspersed nuclear element (SINE) of 227 bp (Ins227bp) insertion polymorphism in the MSTN promoter, and the adjacent SNP BIEC2-417495 have not been studied for their association with racing aptness of the average Thoroughbreds raced in countries with lower status of the racing industry. This study investigated these markers regarding their prevalence and association with performance in common race horses. Markers were genotyped by amplification refractory mutation system-quantitative PCR (ARMS-qPCR) or amplicon melting. Furthermore, we asked whether the Ins227bp marker might theoretically regulate the expression of myostatin by generating a novel target for DNA methylation or by changing binding sites for transcription factors. Putative sites for DNA methylation or binding of transcription factors were predicted by MethPrimer and by the softwares JASPAR, MatInspector and UniPROBE, respectively.

Results

Pairwise linkage disequilibrium between g.66493737 T ≻ C and Ins227bp was high (r² = 0.93). A lower linkage was determined for g.66493737 T ≻ C and BIEC2-417495 (r² = 0.69) as well as for BIEC2-417495 and Ins227bp (r² = 0.76). The estimated frequencies for the presence of Ins227bp (I) indel and the C alleles at g.66493737 T ≻ C and BIEC2-417495 were 0.46, 0.47 and 0.43, respectively. Heterozygotes represented the most abundant genotype at each locus. The best racing distance (BRD) was significantly different between the homozygotes of each SNP (p = 0.01 to 0.03). C allele homozygotes at BIEC2-417495 or g.66493737 T ≻ C, as well as Ins227bp homozygotes earned most money on a mean distance ranging from 1211 to 1230 m. Heterozygotes earned most money on races over 1690 to 1709 m. The BRD for the T/T carriers at both SNP loci and for the SINE-free genotype was 1812 to 1854 m. Other performance parameters were not significantly different between the genotypes, except of the relative success score (RSS). The RSS was significantly slightly better on a distance of ≤1300 m for all carriers of the C allele and the Ins227bp compared to homozygous T genotypes and SINE-negative horses (p = 0.037 to 0.046). For distances of more than 1300 m the RSS was not significantly different between genotypes.

In silico assessment indicated that the Ins227bp promoter insertion might have generated a CpG island and a few novel putative binding sites for transcription factors.

Conclusions

All three target polymorphisms (Ins227bp, g.66493737 T ≻ C, BIEC2-417495) are suitable markers to assess the ability of non-elite Thoroughbreds to race at short or longer distances. The CpG island generated by Ins227bp may cause training-induced silencing of MSTN expression.

Electronic supplementary material

The online version of this article (doi:10.1186/s12917-015-0428-3) contains supplementary material, which is available to authorized users.

A transposable element within the Non-canonical telomerase RNA of Arabidopsis thaliana modulates telomerase in response to DNA damage [corrected]

Long noncoding RNAs (lncRNAs) have emerged as critical factors in many biological processes, but little is known about how their regulatory functions evolved. One of the best-studied lncRNAs is TER, the essential RNA template for telomerase reverse transcriptase. We previously showed that Arabidopsis thaliana harbors three TER isoforms: TER1, TER2 and TER2S. TER1 serves as a canonical telomere template, while TER2 is a novel negative regulator of telomerase activity, induced in response to double-strand breaks (DSBs). TER2 contains a 529 nt intervening sequence that is removed along with 36 nt at the RNA 3’ terminus to generate TER2S, an RNA of unknown function. Here we investigate how A. thaliana TER2 acquired its regulatory function. Using data from the 1,001 Arabidopsis genomes project, we report that the intervening sequence within TER2 is derived from a transposable element termed DSB responsive element (DRE). DRE is found in the TER2 loci of most but not all A. thaliana accessions. By analyzing accessions with (TER2) and without DRE (TER2Δ) we demonstrate that this element is responsible for many of the unique properties of TER2, including its enhanced binding to TERT and telomerase inhibitory function. We show that DRE destabilizes TER2, and further that TER2 induction by DNA damage reflects increased RNA stability and not increased transcription. DRE-mediated changes in TER2 stability thus provide a rapid and sensitive switch to fine-tune telomerase enzyme activity. Altogether, our data shows that invasion of the TER2 locus by a small transposon converted this lncRNA into a DNA damage sensor that modulates telomerase enzyme activity in response to genome assault.

Telomerase is a highly regulated enzyme whose activity is essential for long-term cellular proliferation. In the presence of DNA double-strand breaks (DSBs), telomerase activity must be curtailed to promote faithful DNA repair. We previously showed that the flowering plant Arabidopsis thaliana rapidly down-regulates telomerase in response to DSBs, and further that this mode of regulation is dependent on TER2, a non-canonical telomerase RNA subunit. Here we demonstrate that the unique regulatory properties of TER2 are conveyed by a transposable element (TE) embedded in the TER2 gene. A comparison of A. thaliana accessions with and without the TE revealed that the element increases the binding affinity of TER2 for the telomerase catalytic subunit TERT relative to the canonical telomerase RNA subunit. The TE also increases TER2 turnover. In response to DSBs, TER2 is induced and accumulates in TERT containing complexes in vivo. Thus, invasion of a TE endows TER2 with a DNA damage sensor to rapidly and reversibly modulate enzyme activity in response to genotoxic stress. These findings provide an example of how exaptation of a TE altered the function of a long noncoding RNA. In this case, a duplicated gene (TER2) was used as the platform, and the TE as the tool to engineer a novel mode of telomerase regulation.

Unraveling the Tangled Skein: The Evolution of Transcriptional Regulatory Networks in Development

The molecular and genetic basis for the evolution of anatomical diversity is a major question that has inspired evolutionary and developmental biologists for decades. Because morphology takes form during development, a true comprehension of how anatomical structures evolve requires an understanding of the evolutionary events that alter developmental genetic programs. Vast gene regulatory networks (GRNs) that connect transcription factors to their target regulatory sequences control gene expression in time and space and therefore determine the tissue-specific genetic programs that shape morphological structures. In recent years, many new examples have greatly advanced our understanding of the genetic alterations that modify GRNs to generate newly evolved morphologies. Here, we review several aspects of GRN evolution, including their deep preservation, their mechanisms of alteration, and how they originate to generate novel developmental programs.

What does it take to evolve an enhancer? A simulation-based study of factors influencing the emergence of combinatorial regulation

There is widespread interest today in understanding enhancers, which are regulatory elements typically harboring several transcription factor binding sites and mediating the combinatorial effect of transcription factors on gene expression. The evolution of enhancers poses interesting unanswered questions, for example, the evolutionary time taken for a typical enhancer to emerge or the factors shaping its evolution. Existing approaches to cis-regulatory evolution have often ignored the combinatorial nature and varied biochemical mechanisms of gene regulation encoded in enhancers. We report on our investigation of enhancer evolution through the use of PEBCRES, a framework for evolutionary simulation of enhancers that employs a mechanistic and well-supported sequence-to-expression model to assign fitness to the evolving enhancer genotype. We estimated the time necessary to evolve, from genomic background, enhancers capable of driving complex gene expression patterns similar to those involved in early development in Drosophila. We found the time-to-evolve to range between 0.5 and 10 Myr, and to vary greatly with the target expression pattern, complexity of the real enhancer known to encode that pattern, and the strength of input from specific transcription factors. To our knowledge, this is the first estimate of waiting times for realistic enhancers to evolve. The in silico evolved enhancers had, with a few interesting exceptions, site compositions similar to those seen in real enhancers for the same patterns. Our simulations also revealed that certain features of an enhancer might evolve not due to their biological function but as aids to the evolutionary process itself.

Long-range regulatory interactions at the 4q25 atrial fibrillation risk locus involve PITX2c and ENPEP

BACKGROUND

Recent genome-wide association studies have uncovered genomic loci that underlie an increased risk for atrial fibrillation, the major cardiac arrhythmia in humans. The most significant locus is located in a gene desert at 4q25, approximately 170 kilobases upstream of PITX2, which codes for a transcription factor involved in embryonic left-right asymmetry and cardiac development. However, how this genomic region functionally and structurally relates to PITX2 and atrial fibrillation is unknown.

RESULTS

To characterise its function, we tested genomic fragments from 4q25 for transcriptional activity in a mouse atrial cardiomyocyte cell line and in transgenic mouse embryos, identifying a non-tissue-specific potentiator regulatory element. Chromosome conformation capture revealed that this region physically interacts with the promoter of the cardiac specific isoform of Pitx2. Surprisingly, this regulatory region also interacts with the promoter of the next neighbouring gene, Enpep, which we show to be expressed in regions of the developing mouse heart essential for cardiac electrical activity.

CONCLUSIONS

Our data suggest that de-regulation of both PITX2 and ENPEP could contribute to an increased risk of atrial fibrillation in carriers of disease-associated variants, and show the challenges that we face in the functional analysis of genome-wide disease associations.

Differential regulation of genes by retrotransposons in rice promoters

Rice genome harbors genes and promoters with retrotransposon insertions. There is very little information about their function. The effect of retrotransposon insertions in four rice promoter regions on gene regulation, was investigated using promoter-reporter gene constructs with and without retrotransposons. Differences in expression levels of gus and egfp reporter genes in forward orientation and rfp in reverse orientation were evaluated in rice plants with transient expression employing quantitative RT-PCR analysis, histochemical GUS staining, and eGFP and RFP fluorescent microscopy. The presence of SINE in the promoter 1 (P1) resulted in higher expression levels of the reporter genes, whereas the presence of LINE in P2 or gypsy LTR retrotransposon in P3 reduced expression of the reporter genes. Furthermore, the SINE in P1 acts as an enhancer in contrast with the LINE in P2 and the gypsy LTR retrotransposon in P3 which act as silencers. CTAA and CGG motifs in these retrotransposons are the likely candidates for the downregulation compared to TCTT motif (SINE) which is a candidate for the upregulation of gene expression. The effect of retrotransposons on gene regulation correlated with the earlier investigation of conservation patterns of these four retrotransposon insertions in several rice accessions implying their evolutionary significance.

Non-coding RNA: what is functional and what is junk?

The genomes of large multicellular eukaryotes are mostly comprised of non-protein coding DNA. Although there has been much agreement that a small fraction of these genomes has important biological functions, there has been much debate as to whether the rest contributes to development and/or homeostasis. Much of the speculation has centered on the genomic regions that are transcribed into RNA at some low level. Unfortunately these RNAs have been arbitrarily assigned various names, such as “intergenic RNA,” “long non-coding RNAs” etc., which have led to some confusion in the field. Many researchers believe that these transcripts represent a vast, unchartered world of functional non-coding RNAs (ncRNAs), simply because they exist. However, there are reasons to question this Panglossian view because it ignores our current understanding of how evolution shapes eukaryotic genomes and how the gene expression machinery works in eukaryotic cells. Although there are undoubtedly many more functional ncRNAs yet to be discovered and characterized, it is also likely that many of these transcripts are simply junk. Here, we discuss how to determine whether any given ncRNA has a function. Importantly, we advocate that in the absence of any such data, the appropriate null hypothesis is that the RNA in question is junk.

Transposable element-driven transcript diversification and its relevance to genetic disorders

The human genome project and subsequent gene annotation projects have shown that the human genome contains 22,000-25,000 functional genes. Therefore, it is believed that the diversity of protein repertoire is achieved by the alternative splicing (AS) mechanism. Transposable elements (TEs) are mobile in nature and can therefore alter their position in the genome. The insertion of TEs into a new gene region can result in AS of a particular transcript through various mechanisms, including intron retention, and alternative donor or acceptor splice sites. TE-derived AS is thought to have played a part in primate evolution and in hominid radiation. However, TE-derived AS or genetic instability may sometimes result in genetic disorders. For the past two decades, numerous studies have been performed on TEs and their role in genomes. Accumulating evidence shows that the term 'junk DNA', previously used for TEs is a misnomer. Recent research has indicated that TEs may have clinical potential. However, to explore the feasibility of using TEs in clinical practice, additional studies are required. This review summarizes the available literature on TE-derived AS, alternative promoter, and alternative polyadenylation. The review covers the effects of TEs on coding genes and their clinical implications, and provides our perspectives and directions for future research.

Transposable elements contribute to activation of maize genes in response to abiotic stress

Transposable elements (TEs) account for a large portion of the genome in many eukaryotic species. Despite their reputation as “junk” DNA or genomic parasites deleterious for the host, TEs have complex interactions with host genes and the potential to contribute to regulatory variation in gene expression. It has been hypothesized that TEs and genes they insert near may be transcriptionally activated in response to stress conditions. The maize genome, with many different types of TEs interspersed with genes, provides an ideal system to study the genome-wide influence of TEs on gene regulation. To analyze the magnitude of the TE effect on gene expression response to environmental changes, we profiled gene and TE transcript levels in maize seedlings exposed to a number of abiotic stresses. Many genes exhibit up- or down-regulation in response to these stress conditions. The analysis of TE families inserted within upstream regions of up-regulated genes revealed that between four and nine different TE families are associated with up-regulated gene expression in each of these stress conditions, affecting up to 20% of the genes up-regulated in response to abiotic stress, and as many as 33% of genes that are only expressed in response to stress. Expression of many of these same TE families also responds to the same stress conditions. The analysis of the stress-induced transcripts and proximity of the transposon to the gene suggests that these TEs may provide local enhancer activities that stimulate stress-responsive gene expression. Our data on allelic variation for insertions of several of these TEs show strong correlation between the presence of TE insertions and stress-responsive up-regulation of gene expression. Our findings suggest that TEs provide an important source of allelic regulatory variation in gene response to abiotic stress in maize.

Transposable elements are mobile DNA elements that are a prevalent component of many eukaryotic genomes. While transposable elements can often have deleterious effects through insertions into protein-coding genes they may also contribute to regulatory variation of gene expression. There are a handful of examples in which specific transposon insertions contribute to regulatory variation of nearby genes, particularly in response to environmental stress. We sought to understand the genome-wide influence of transposable elements on gene expression responses to abiotic stress in maize, a plant with many families of transposable elements located in between genes. Our analysis suggests that a small number of maize transposable element families may contribute to the response of nearby genes to abiotic stress by providing stress-responsive enhancer-like functions. The specific insertions of transposable elements are often polymorphic within a species. Our data demonstrate that allelic variation for insertions of the transposable elements associated with stress-responsive expression can contribute to variation in the regulation of nearby genes. Thus novel insertions of transposable elements provide a potential mechanism for genes to acquire cis-regulatory influences that could contribute to heritable variation for stress response.

Functional annotation of the Ophiostoma novo-ulmi genome: insights into the phytopathogenicity of the fungal agent of Dutch elm disease

The ascomycete fungus Ophiostoma novo-ulmi is responsible for the pandemic of Dutch elm disease that has been ravaging Europe and North America for 50 years. We proceeded to annotate the genome of the O. novo-ulmi strain H327 that was sequenced in 2012. The 31.784-Mb nuclear genome (50.1% GC) is organized into 8 chromosomes containing a total of 8,640 protein-coding genes that we validated with RNA sequencing analysis. Approximately 53% of these genes have their closest match to Grosmannia clavigera kw1407, followed by 36% in other close Sordariomycetes, 5% in other Pezizomycotina, and surprisingly few (5%) orphans. A relatively small portion (∼3.4%) of the genome is occupied by repeat sequences; however, the mechanism of repeat-induced point mutation appears active in this genome. Approximately 76% of the proteins could be assigned functions using Gene Ontology analysis; we identified 311 carbohydrate-active enzymes, 48 cytochrome P450s, and 1,731 proteins potentially involved in pathogen-host interaction, along with 7 clusters of fungal secondary metabolites. Complementary mating-type locus sequencing, mating tests, and culturing in the presence of elm terpenes were conducted. Our analysis identified a specific genetic arsenal impacting the sexual and vegetative growth, phytopathogenicity, and signaling/plant-defense-degradation relationship between O. novo-ulmi and its elm host and insect vectors.

Widespread contribution of transposable elements to the innovation of gene regulatory networks

Transposable elements (TEs) have been shown to contain functional binding sites for certain transcription factors (TFs). However, the extent to which TEs contribute to the evolution of TF binding sites is not well known. We comprehensively mapped binding sites for 26 pairs of orthologous TFs in two pairs of human and mouse cell lines (representing two cell lineages), along with epigenomic profiles, including DNA methylation and six histone modifications. Overall, we found that 20% of binding sites were embedded within TEs. This number varied across different TFs, ranging from 2% to 40%. We further identified 710 TF–TE relationships in which genomic copies of a TE subfamily contributed a significant number of binding peaks for a TF, and we found that LTR elements dominated these relationships in human. Importantly, TE-derived binding peaks were strongly associated with open and active chromatin signatures, including reduced DNA methylation and increased enhancer-associated histone marks. On average, 66% of TE-derived binding events were cell type-specific with a cell type-specific epigenetic landscape. Most of the binding sites contributed by TEs were species-specific, but we also identified binding sites conserved between human and mouse, the functional relevance of which was supported by a signature of purifying selection on DNA sequences of these TEs. Interestingly, several TFs had significantly expanded binding site landscapes only in one species, which were linked to species-specific gene functions, suggesting that TEs are an important driving force for regulatory innovation. Taken together, our data suggest that TEs have significantly and continuously shaped gene regulatory networks during mammalian evolution.

Volatile evolution of long noncoding RNA repertoires: mechanisms and biological implications

Thousands of genes encoding long noncoding RNAs (lncRNAs) have been identified in all vertebrate genomes thus far examined. The list of lncRNAs partaking in arguably important biochemical, cellular, and developmental activities is steadily growing. However, it is increasingly clear that lncRNA repertoires are subject to weak functional constraint and rapid turnover during vertebrate evolution. We discuss here some of the factors that may explain this apparent paradox, including relaxed constraint on sequence to maintain lncRNA structure/function, extensive redundancy in the regulatory circuits in which lncRNAs act, as well as adaptive and non-adaptive forces such as genetic drift. We explore the molecular mechanisms promoting the birth and rapid evolution of lncRNA genes, with an emphasis on the influence of bidirectional transcription and transposable elements, two pervasive features of vertebrate genomes. Together these properties reveal a remarkably dynamic and malleable noncoding transcriptome which may represent an important source of robustness and evolvability.

Population genomics of transposable elements in Drosophila

Studies of the population dynamics of transposable elements (TEs) in Drosophila melanogaster indicate that consistent forces are affecting TEs independently of their modes of transposition and regulation. New sequencing technologies enable biologists to sample genomes at an unprecedented scale in order to quantify genome-wide polymorphism for annotated and novel TE insertions. In this review, we first present new insights gleaned from high-throughput data for population genomics studies of D. melanogaster. We then consider the latest population genomics models for TE evolution and present examples of functional evidence revealed by genome-wide studies of TE population dynamics in D. melanogaster. Although most of the TE insertions are deleterious or neutral, some TE insertions increase the fitness of the individual that carries them and play a role in genome adaptation.

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.