Regulatory activities of transposable elements: from conflicts to benefits

Integrated analysis sheds light on evolutionary trajectories of young transcription start sites in the human genome

Understanding the molecular mechanisms and evolution of the gene regulatory system remains a major challenge in biology. Transcription start sites (TSSs) are especially interesting because they are central to initiating gene expression. Previous studies revealed widespread transcription initiation and fast turnover of TSSs in mammalian genomes. Yet, how new TSSs originate and how they evolve over time remain poorly understood. To address these questions, we analyzed ∼200,000 human TSSs by integrating evolutionary (inter- and intra-species) and functional genomic data, particularly focusing on evolutionarily young TSSs that emerged in the primate lineage. TSSs were grouped according to their evolutionary age using sequence alignment information as a proxy. Comparisons of young and old TSSs revealed that (1) new TSSs emerge through a combination of intrinsic factors, like the sequence properties of transposable elements and tandem repeats, and extrinsic factors such as their proximity to existing regulatory modules; (2) new TSSs undergo rapid evolution that reduces the inherent instability of repeat sequences associated with a high propensity of TSS emergence; and (3) once established, the transcriptional competence of surviving TSSs is gradually enhanced, with evolutionary changes subject to temporal (fewer regulatory changes in younger TSSs) and spatial constraints (fewer regulatory changes in more isolated TSSs). These findings advance our understanding of how regulatory innovations arise in the genome throughout evolution and highlight the genomic robustness and evolvability in these processes.

Proliferation of Regulatory DNA Elements Derived from Transposable Elements in the Maize Genome

Genomic regions free of nucleosomes, which are hypersensitive to DNase I digestion, are known as DNase I hypersensitive sites (DHSs) and frequently contain cis-regulatory DNA elements. To investigate their prevalence and characteristics in maize (Zea mays), we developed high-resolution genome-wide DHS maps using a modified DNase-seq technique. Maize DHSs exhibit depletion of nucleosomes and low levels of DNA methylation and are enriched with conserved noncoding sequences (CNSs). We developed a protoplast-based transient transformation assay to assess the potential gene expression enhancer and/or promoter functions associated with DHSs, which showed that more than 80% of DHSs overlapping with CNSs showed an enhancer function. Strikingly, nearly 25% of maize DHSs were derived from transposable elements (TEs), including both class I and class II transposons. Interestingly, TE-derived DHSs (teDHSs) homologous to retrotransposons were enriched with sequences related to the intrinsic cis-regulatory elements within the long terminal repeats of retrotransposons. We demonstrate that more than 80% of teDHSs can drive transcription of a reporter gene in protoplast assays. These results reveal the widespread occurrence of TE-derived cis-regulatory sequences and suggest that teDHSs play a major role in transcriptional regulation in maize.

Drosophila relics hobo and hobo-MITEs transposons as raw material for new regulatory networks

Hypermutable strains of Drosophila simulans have been studied for 20 years. Several mutants were isolated and characterized, some of which had phenotypes associated with alteration in development; for example, showing ectopic legs with eyes being expressed in place of antennae. The causal agent of this hypermutability is a non-autonomous hobo-related sequence (hoboVA). Around 100 mobilizable copies of this element are present in the D. simulans genome, and these are likely mobilized by the autonomous and canonical hobo element. We have shown that hoboVA has transcription factor binding sites for the developmental genes, hunchback and even-skipped, and that this transposon is expressed in embryos, following the patterns of these genes. We suggest that hobo and hobo-related elements can be material for the emergence of new regulatory networks.

Integrative analysis of large scale transcriptome data draws a comprehensive landscape of Phaeodactylum tricornutum genome and evolutionary origin of diatoms

Diatoms are one of the most successful and ecologically important groups of eukaryotic phytoplankton in the modern ocean. Deciphering their genomes is a key step towards better understanding of their biological innovations, evolutionary origins, and ecological underpinnings. Here, we have used 90 RNA-Seq datasets from different growth conditions combined with published expressed sequence tags and protein sequences from multiple taxa to explore the genome of the model diatom Phaeodactylum tricornutum, and introduce 1,489 novel genes. The new annotation additionally permitted the discovery of extensive alternative splicing in diatoms, including intron retention and exon skipping, which increase the diversity of transcripts generated in changing environments. In addition, we have used up-to-date reference sequence libraries to dissect the taxonomic origins of diatom genes. We show that the P. tricornutum genome is enriched in lineage-specific genes, with up to 47% of the gene models present only possessing orthologues in other stramenopile groups. Finally, we have performed a comprehensive de novo annotation of repetitive elements showing novel classes of transposable elements such as SINE, MITE and TRIM/LARD. This work provides a solid foundation for future studies of diatom gene function, evolution and ecology.

NuRD and CAF-1-mediated silencing of the D4Z4 array is modulated by DUX4-induced MBD3L proteins

The DUX4 transcription factor is encoded by a retrogene embedded in each unit of the D4Z4 macrosatellite repeat. DUX4 is normally expressed in the cleavage-stage embryo, whereas chromatin repression prevents DUX4 expression in most somatic tissues. Failure of this repression causes facioscapulohumeral muscular dystrophy (FSHD) due to mis-expression of DUX4 in skeletal muscle. In this study, we used CRISPR/Cas9 engineered chromatin immunoprecipitation (enChIP) locus-specific proteomics to characterize D4Z4-associated proteins. These and other approaches identified the Nucleosome Remodeling Deacetylase (NuRD) and Chromatin Assembly Factor 1 (CAF-1) complexes as necessary for DUX4 repression in human skeletal muscle cells and induced pluripotent stem (iPS) cells. Furthermore, DUX4-induced expression of MBD3L proteins partly relieved this repression in FSHD muscle cells. Together, these findings identify NuRD and CAF-1 as mediators of DUX4 chromatin repression and suggest a mechanism for the amplification of DUX4 expression in FSHD muscle cells.

Exaptation of transposable element coding sequences

Transposable elements (TEs) are mobile genetic elements that were once perceived as merely selfish, but are now recognized as potent agents of adaptation. One way TEs contribute to genome evolution is through TE exaptation, a process whereby TEs, which usually persist by replicating in the genome, transform into novel host genes, which thereafter persist by conferring phenotypic benefits. Exapted TEs are known to contribute diverse and vital functions, and may facilitate punctuated equilibrium, yet we have little understanding about the process of TE exaptation. In order to facilitate our understanding of how TE coding sequences may become exapted, here we incorporate the findings of recent publications into a framework and six-step model.

RNAi is a critical determinant of centromere evolution in closely related fungi

The “centromere paradox” refers to rapidly evolving and highly diverse centromere DNA sequences even in closely related eukaryotes. However, factors contributing to this rapid divergence are largely unknown. Here, we identified large regional, LTR retrotransposon-rich centromeres in a group of human fungal pathogens belonging to the Cryptococcus species complex. We provide evidence that loss-of-functional RNAi machinery and possibly cytosine DNA methylation trigger instability of the genome by activation of centromeric retrotransposons presumably suppressed by RNAi. We propose that RNAi, together with cytosine DNA methylation, serves as a critical determinant that maintains repetitive transposon-rich centromere structures. This study explores the direct link between RNAi and centromere structure evolution.

The centromere DNA locus on a eukaryotic chromosome facilitates faithful chromosome segregation. Despite performing such a conserved function, centromere DNA sequence as well as the organization of sequence elements is rapidly evolving in all forms of eukaryotes. The driving force that facilitates centromere evolution remains an enigma. Here, we studied the evolution of centromeres in closely related species in the fungal phylum of Basidiomycota. Using ChIP-seq analysis of conserved inner kinetochore proteins, we identified centromeres in three closely related Cryptococcus species: two of which are RNAi-proficient, while the other lost functional RNAi. We find that the centromeres in the RNAi-deficient species are significantly shorter than those of the two RNAi-proficient species. While centromeres are LTR retrotransposon-rich in all cases, the RNAi-deficient species lost all full-length retroelements from its centromeres. In addition, centromeres in RNAi-proficient species are associated with a significantly higher level of cytosine DNA modifications compared with those of RNAi-deficient species. Furthermore, when an RNAi-proficient Cryptococcus species and its RNAi-deficient mutants were passaged under similar conditions, the centromere length was found to be occasionally shortened in RNAi mutants. In silico analysis of predicted centromeres in a group of closely related Ustilago species, also belonging to the Basidiomycota, were found to have undergone a similar transition in the centromere length in an RNAi-dependent fashion. Based on the correlation found in two independent basidiomycetous species complexes, we present evidence suggesting that the loss of RNAi and cytosine DNA methylation triggered transposon attrition, which resulted in shortening of centromere length during evolution.

Transposable element activity, genome regulation and human health

A convergence of novel genome analysis technologies is enabling population genomic studies of human transposable elements (TEs). Population surveys of human genome sequences have uncovered thousands of individual TE insertions that segregate as common genetic variants, i.e. TE polymorphisms. These recent TE insertions provide an important source of naturally occurring human genetic variation. Investigators are beginning to leverage population genomic data sets to execute genome-scale association studies for assessing the phenotypic impact of human TE polymorphisms. For example, the expression quantitative trait loci (eQTL) analytical paradigm has recently been used to uncover hundreds of associations between human TE insertion variants and gene expression levels. These include population-specific gene regulatory effects as well as coordinated changes to gene regulatory networks. In addition, analyses of linkage disequilibrium patterns with previously characterized genome-wide association study (GWAS) trait variants have uncovered TE insertion polymorphisms that are likely causal variants for a variety of common complex diseases. Gene regulatory mechanisms that underlie specific disease phenotypes have been proposed for a number of these trait associated TE polymorphisms. These new population genomic approaches hold great promise for understanding how ongoing TE activity contributes to functionally relevant genetic variation within and between human populations.

Individual retrotransposon integrants are differentially controlled by KZFP/KAP1-dependent histone methylation, DNA methylation and TET-mediated hydroxymethylation in naïve embryonic stem cells

Background

The KZFP/KAP1 (KRAB zinc finger proteins/KRAB-associated protein 1) system plays a central role in repressing transposable elements (TEs) and maintaining parent-of-origin DNA methylation at imprinting control regions (ICRs) during the wave of genome-wide reprogramming that precedes implantation. In naïve murine embryonic stem cells (mESCs), the genome is maintained highly hypomethylated by a combination of TET-mediated active demethylation and lack of de novo methylation, yet KAP1 is tethered by sequence-specific KZFPs to ICRs and TEs where it recruits histone and DNA methyltransferases to impose heterochromatin formation and DNA methylation.

Results

Here, upon removing either KAP1 or the cognate KZFP, we observed rapid TET2-dependent accumulation of 5hmC at both ICRs and TEs. In the absence of the KZFP/KAP1 complex, ICRs lost heterochromatic histone marks and underwent both active and passive DNA demethylation. For KAP1-bound TEs, 5mC hydroxylation correlated with transcriptional reactivation. Using RNA-seq, we further compared the expression profiles of TEs upon Kap1 removal in wild-type, Dnmt and Tet triple knockout mESCs. While we found that KAP1 represents the main effector of TEs repression in all three settings, we could additionally identify specific groups of TEs further controlled by DNA methylation. Furthermore, we observed that in the absence of TET proteins, activation upon Kap1 depletion was blunted for some TE integrants and increased for others.

Conclusions

Our results indicate that the KZFP/KAP1 complex maintains heterochromatin and DNA methylation at ICRs and TEs in naïve embryonic stem cells partly by protecting these loci from TET-mediated demethylation. Our study further unveils an unsuspected level of complexity in the transcriptional control of the endovirome by demonstrating often integrant-specific differential influences of histone-based heterochromatin modifications, DNA methylation and 5mC oxidation in regulating TEs expression.

Electronic supplementary material

The online version of this article (10.1186/s13072-018-0177-1) contains supplementary material, which is available to authorized users.

The epigenetic alterations of endogenous retroelements in aging

Endogenous retroelements, transposons that mobilize through RNA intermediates, include some of the most abundant repetitive sequences of the human genome, such as Alu and LINE-1 sequences, and human endogenous retroviruses. Recent discoveries demonstrate that these mobile genetic elements not only act as intragenomic parasites, but also exert regulatory roles in living cells. The risk of genomic instability represented by endogenous retroelements is normally counteracted by a series of epigenetic control mechanisms which include, among the most important, CpG DNA methylation. Indeed, most of the genomic CpG sites subjected to DNA methylation in the nuclear DNA are carried by these repetitive elements. As other parts of the genome, endogenous retroelements and other transposable elements are subjected to deep epigenetic alterations during aging, repeatedly observed in the context of organismal and cellular senescence, in human and other species. This review summarizes the current status of knowledge about the epigenetic alterations occurring in this large, non-genic portion of the genome in aging and age-related conditions, with a focus on the causes and the possible functional consequences of these alterations.

Transposable elements shape the human proteome landscape via formation of cis-acting upstream open reading frames

Transposons are major drivers of mammalian genome evolution. To obtain new insights into the contribution of transposons to the regulation of protein translation, we here examined how transposons affected the genesis and function of upstream open reading frames (uORFs), which serve as cis-acting elements to regulate translation from annotated ORFs (anORFs) located downstream of the uORFs in eukaryotic mRNAs. Among 39,786 human uORFs, 3,992 had ATG trinucleotides of a transposon origin, termed "transposon-derived upstream ATGs" or TuATGs. Luciferase reporter assays suggested that many TuATGs modulate translation from anORFs. Comparisons with transposon consensus sequences revealed that most TuATGs were generated by nucleotide substitutions in non-ATG trinucleotides of integrated transposons. Among these non-ATG trinucleotides, GTG and ACG were converted into TuATGs more frequently, indicating a CpG methylation-mediated process of TuATG formation. Interestingly, it is likely that this process accelerated human-specific upstream ATG formation within transposon sequences in 5' untranslated regions after divergence between human and nonhuman primates. Methylation-mediated TuATG formation seems to be ongoing in the modern human population and could alter the expression of disease-related proteins. This study shows that transposons have potentially been shaping the human proteome landscape via cis-acting uORF creation.

The role of transposable elements in functional evolution of amphioxus genome: the case of opsin gene family

Transposable elements (TEs) are able to jump to new locations (transposition) in the genome, usually after replication. They constitute the so-called selfish or junk DNA and take over large proportions of some genomes. Due to their ability to move around they can change the DNA landscape of genomes and are therefore a rich source of innovation in genes and gene regulation. Surge of sequence data in the past years has significantly facilitated large scale comparative studies. Cephalochordates have been regarded as a useful proxy to ancestral chordate condition partially due to the comparatively slow evolutionary rate at morphological and genomic level. In this study, we used opsin gene family from three Branchiostoma species as a window into cephalochordate genome evolution. We compared opsin complements in terms of family size, gene structure and sequence allowing us to identify gene duplication and gene loss events. Furthermore, analysis of the opsin containing genomic loci showed that they are populated by TEs. In summary, we provide evidence of the way transposable elements may have contributed to the evolution of opsin gene family and to the shaping of cephalochordate genomes in general.

Mammalian transposable elements and their impacts on genome evolution

Transposable elements (TEs) are genetic elements with the ability to mobilize and replicate themselves in a genome. Mammalian genomes are dominated by TEs, which can reach copy numbers in the hundreds of thousands. As a result, TEs have had significant impacts on mammalian evolution. Here we summarize the current understanding of TE content in mammal genomes and find that, with a few exceptions, most fall within a predictable range of observations. First, one third to one half of the genome is derived from TEs. Second, most mammalian genomes are dominated by LINE and SINE retrotransposons, more limited LTR retrotransposons, and minimal DNA transposon accumulation. Third, most mammal genome contains at least one family of actively accumulating retrotransposon. Finally, horizontal transfer of TEs among lineages is rare. TE exaptation events are being recognized with increasing frequency. Despite these beneficial aspects of TE content and activity, the majority of TE insertions are neutral or deleterious. To limit the deleterious effects of TE proliferation, the genome has evolved several defense mechanisms that act at the epigenetic, transcriptional, and post-transcriptional levels. The interaction between TEs and these defense mechanisms has led to an evolutionary arms race where TEs are suppressed, evolve to escape suppression, then are suppressed again as the defense mechanisms undergo compensatory change. The result is complex and constantly evolving interactions between TEs and host genomes.

Bioinformatic analysis of regulation of natural antisense transcripts by transposable elements in human mRNA

Non-coding RNA is no longer considered to be "junk" DNA, based on evidence uncovered in recent decades. In particular, the important role played by natural antisense transcripts (NATs) in regulating the expression of genes is receiving increasing attention. However, the regulatory mechanisms of NATs remain incompletely understood. It is well-known that the insertion of transposable elements (TEs) can affect gene transcription. Using a bioinformatics approach, we identified NATs using human mRNA sequences from the UCSC Genome Browser Database. Our in silico analysis identified 1079 NATs and 700 sense-antisense gene pairs. We identified 179 NATs that showed evidence of having been affected by TEs during cellular gene expression. These findings may provide an understanding of the complex regulation mechanisms of NATs. If our understanding of NATs as modulators of gene expression is further enhanced, we can develop ways to control gene expression.

Profiling of Human Molecular Pathways Affected by Retrotransposons at the Level of Regulation by Transcription Factor Proteins

Endogenous retroviruses and retrotransposons also termed retroelements (REs) are mobile genetic elements that were active until recently in human genome evolution. REs regulate gene expression by actively reshaping chromatin structure or by directly providing transcription factor binding sites (TFBSs). We aimed to identify molecular processes most deeply impacted by the REs in human cells at the level of TFBS regulation. By using ENCODE data, we identified ~2 million TFBS overlapping with putatively regulation-competent human REs located in 5-kb gene promoter neighborhood (~17% of all TFBS in promoter neighborhoods; ~9% of all RE-linked TFBS). Most of REs hosting TFBS were highly diverged repeats, and for the evolutionary young (0–8% diverged) elements we identified only ~7% of all RE-linked TFBS. The gene-specific distributions of RE-linked TFBS generally correlated with the distributions for all TFBS. However, several groups of molecular processes were highly enriched in the RE-linked TFBS regulation. They were strongly connected with the immunity and response to pathogens, with the negative regulation of gene transcription, ubiquitination, and protein degradation, extracellular matrix organization, regulation of STAT signaling, fatty acids metabolism, regulation of GTPase activity, protein targeting to Golgi, regulation of cell division and differentiation, development and functioning of perception organs and reproductive system. By contrast, the processes most weakly affected by the REs were linked with the conservative aspects of embryo development. We also identified differences in the regulation features by the younger and older fractions of the REs. The regulation by the older fraction of the REs was linked mainly with the immunity, cell adhesion, cAMP, IGF1R, Notch, Wnt, and integrin signaling, neuronal development, chondroitin sulfate and heparin metabolism, and endocytosis. The younger REs regulate other aspects of immunity, cell cycle progression and apoptosis, PDGF, TGF beta, EGFR, and p38 signaling, transcriptional repression, structure of nuclear lumen, catabolism of phospholipids, and heterocyclic molecules, insulin and AMPK signaling, retrograde Golgi-ER transport, and estrogen signaling. The immunity-linked pathways were highly represented in both categories, but their functional roles were different and did not overlap. Our results point to the most quickly evolving molecular pathways in the recent and ancient evolution of human genome.

Activation of Lineage Regulators and Transposable Elements across a Pluripotent Spectrum

Embryonic stem cells (ESCs) are characterized by the pluripotent capacity to generate all embryonic lineages. Here, we show that ESCs can occupy a spectrum of distinct transcriptional and epigenetic states in response to varied extrinsic conditions. This spectrum broadly corresponds to a developmental continuum of pluripotency and is coupled with a gradient of increasing global DNA methylation. Each pluripotent state is linked with activation of distinct classes of transposable elements (TEs), which in turn influence ESCs through generating chimeric transcripts. Moreover, varied ESC culture parameters differentially license heterogeneous activation of master lineage regulators, including Sox1, Gata4, or Blimp1, and influence differentiation. Activation of Blimp1 is prevalent in 2i (without LIF) conditions, and marks a dynamic primordial germ cell (PGC)-like sub-state that is directly repressed by Klf4 downstream of LIF/STAT3 signaling. Thus, extrinsic cues establish a spectrum of pluripotent states, in part by modulating sub-populations, as well as directing the transcriptome, epigenome, and TE.

•

Diverse culture parameters establish a spectrum of ESC pluripotency

•

Activation of distinct transposable elements in each pluripotent state

•

ESC conditions influence DNA methylation and response to differentiation cues

•

Distinct heterogeneities including a PGC-like state regulated by LIF and KLF4

De Novo characterization of transcriptomes from two North American Papaipema stem-borers (Lepidoptera: Noctuidae)

Stem-borers in the genus Papaipema (Lepidoptera: Noctuidae) range from highly polyphagous agricultural pests to specialists on more than 20 families of flowering plants, many of them highly toxic. Papaipema is the largest genus of noctuids endemic to North America and provides an excellent study system for the evolution of noctuid host plant use. To improve the availability of genomic resources for such investigations, we performed de novo transcriptome sequencing and assembly for two specialist Papaipema with unusual larval hosts: P. speciosissima, which is associated with ferns, and the undescribed P. “sp. 4,” which is associated with bamboo. The resulting transcriptomes were similar in terms of completeness, gene count, and gene identity, but we identified some 8,000 genes (~17% of each transcriptome) not shared between the two species. While some of these have identifiable orthologs in other Lepidoptera, ~5% of each transcriptome consists of species-specific genes. We examine the function of these genes and find that almost half have retrotransposon-related functional domains. The potential role of species-specific genes is discussed, and the expansion of certain retrotransposon families in Papaipema is examined.

Identifying co-opted transposable elements using comparative epigenomics

The human genome gives rise to different epigenomic landscapes that define each cell type and can be deregulated in disease. Recent efforts by ENCODE, the NIH Roadmap and the International Human Epigenome Consortium (IHEC) have made significant advances towards assembling reference epigenomic maps of various tissues. Notably, these projects have found that approximately 80% of human DNA was biochemically active in at least one epigenomic assay while only approximately 10% of the sequence displayed signs of purifying selection. Given that transposable elements (TEs) make up at least 50% of the human genome and can be actively transcribed or act as regulatory elements either for their own purposes or be co-opted for the benefit of their host; we are interested in exploring their overall contribution to the "functional" genome. Traditional methods used to identify functional DNA have relied on comparative genomics, conservation analysis and low throughput validation assays. To discover co-opted TEs, and distinguish them from noisy genomic elements, we argue that comparative epigenomic methods will also be important.

Analysis of Transposable Elements in Coccidioides Species

Coccidioides immitis and C. posadasii are primary pathogenic fungi that cause disease in immunologically-normal animals and people. The organism is found exclusively in arid regions of the Southwestern United States, Mexico, and South America, but not in other parts of the world. This study is a detailed analysis of the transposable elements (TE) in Coccidioides spp. As is common in most fungi, Class I and Class II transposons were identified and the LTR Gypsy superfamily is the most common. The minority of Coccidioides Gypsy transposons contained regions highly homologous to polyprotein domains. Phylogenetic analysis of the integrase and reverse transcriptase sequences revealed that many, but not all, of the Gypsy reverse transcriptase and integrase domains clustered by species suggesting extensive transposition after speciation of the two Coccidiodies spp. The TEs were clustered and the distribution is enriched for the ends on contigs. Analysis of gene expression data from C. immitis found that protein-coding genes within 1 kB of hAT or Gypsy TEs were poorly expressed. The expression of C. posadasii genes within 1 kB of Gypsy TEs was also significantly lower compared to all genes but the difference in expression was smaller than C. immitis. C. posadasii orthologs of C. immitis Gyspsy-associated genes were also likely to be TE-associated. In both C. immitis and C. posadasii the TEs were preferentially associated with genes annotated with protein kinase gene ontology terms. These observations suggest that TE may play a role in influencing gene expression in Coccidioides spp. Our hope is that these bioinformatic studies of the potential TE influence on expression and evolution of Coccidioides will prompt the development of testable hypotheses to better understand the role of TEs in the biology and gene regulation of Coccidioides spp.

Evidence for convergent evolution of SINE-directed Staufen-mediated mRNA decay

Primate-specific Alu short interspersed elements (SINEs) as well as rodent-specific B and ID (B/ID) SINEs can promote Staufen-mediated decay (SMD) when present in mRNA 3'-untranslated regions (3'-UTRs). The transposable nature of SINEs, their presence in long noncoding RNAs, their interactions with Staufen, and their rapid divergence in different evolutionary lineages suggest they could have generated substantial modification of posttranscriptional gene-control networks during mammalian evolution. Some of the variation in SMD regulation produced by SINE insertion might have had a similar regulatory effect in separate mammalian lineages, leading to parallel evolution of the Staufen network by independent expansion of lineage-specific SINEs. To explore this possibility, we searched for orthologous gene pairs, each carrying a species-specific 3'-UTR SINE and each regulated by SMD, by measuring changes in mRNA abundance after individual depletion of two SMD factors, Staufen1 (STAU1) and UPF1, in both human and mouse myoblasts. We identified and confirmed orthologous gene pairs with 3'-UTR SINEs that independently function in SMD control of myoblast metabolism. Expanding to other species, we demonstrated that SINE-directed SMD likely emerged in both primate and rodent lineages >20-25 million years ago. Our work reveals a mechanism for the convergent evolution of posttranscriptional gene regulatory networks in mammals by species-specific SINE transposition and SMD.

Transposable elements: genome innovation, chromosome diversity, and centromere conflict

Although it was nearly 70 years ago when transposable elements (TEs) were first discovered “jumping” from one genomic location to another, TEs are now recognized as contributors to genomic innovations as well as genome instability across a wide variety of species. In this review, we illustrate the ways in which active TEs, specifically retroelements, can create novel chromosome rearrangements and impact gene expression, leading to disease in some cases and species-specific diversity in others. We explore the ways in which eukaryotic genomes have evolved defense mechanisms to temper TE activity and the ways in which TEs continue to influence genome structure despite being rendered transpositionally inactive. Finally, we focus on the role of TEs in the establishment, maintenance, and stabilization of critical, yet rapidly evolving, chromosome features: eukaryotic centromeres. Across centromeres, specific types of TEs participate in genomic conflict, a balancing act wherein they are actively inserting into centromeric domains yet are harnessed for the recruitment of centromeric histones and potentially new centromere formation.

The Neuronal Gene Arc Encodes a Repurposed Retrotransposon Gag Protein that Mediates Intercellular RNA Transfer

The neuronal gene Arc is essential for long-lasting information storage in the mammalian brain, mediates various forms of synaptic plasticity, and has been implicated in neurodevelopmental disorders. However, little is known about Arc's molecular function and evolutionary origins. Here, we show that Arc self-assembles into virus-like capsids that encapsulate RNA. Endogenous Arc protein is released from neurons in extracellular vesicles that mediate the transfer of Arc mRNA into new target cells, where it can undergo activity-dependent translation. Purified Arc capsids are endocytosed and are able to transfer Arc mRNA into the cytoplasm of neurons. These results show that Arc exhibits similar molecular properties to retroviral Gag proteins. Evolutionary analysis indicates that Arc is derived from a vertebrate lineage of Ty3/gypsy retrotransposons, which are also ancestors to retroviruses. These findings suggest that Gag retroelements have been repurposed during evolution to mediate intercellular communication in the nervous system.

Transposable Element Mediated Innovation in Gene Regulatory Landscapes of Cells: Re-Visiting the "Gene-Battery" Model

Transposable elements (TEs) are no longer considered to be "junk" DNA. Here, we review how TEs can impact gene regulation systematically. TEs encode various regulatory elements that enables them to regulate gene expression. RJ Britten and EH Davidson hypothesized that TEs can integrate the function of various transcriptional regulators into gene regulatory networks. Uniquely TEs can deposit regulatory sites across the genome when they transpose, and thereby bring multiple genes under control of the same regulatory logic. Several studies together have robustly established that TEs participate in embryonic development and oncogenesis. We discuss the regulatory characteristics of TEs in context of evolution to understand the extent of their impact on gene networks. Understanding these features of TEs is central to future investigations of TEs in cellular processes and phenotypic presentations, which are applicable to development and disease studies. We re-visit the Britten-Davidson "gene-battery" model and understand the genetic and transcriptional impact of TEs in innovating gene regulatory networks.

PlaNC-TE: a comprehensive knowledgebase of non-coding RNAs and transposable elements in plants

Transposable elements (TEs) play an essential role in the genetic variability of eukaryotic species. In plants, they may comprise up to 90% of the total genome. Non-coding RNAs (ncRNAs) are known to control gene expression and regulation. Although the relationship between ncRNAs and TEs is known, obtaining the organized data for sequenced genomes is not straightforward. In this study, we describe the PlaNC-TE (http://planc-te.cp.utfpr.edu.br), a user-friendly portal harboring a knowledgebase created by integrating and analysing plant ncRNA-TE data. We identified a total of 14 350 overlaps between ncRNAs and TEs in 40 plant genomes. The database allows users to browse, search and download all ncRNA and TE data analysed. Overall, PlaNC-TE not only organizes data and provides insights about the relationship between ncRNA and TEs in plants but also helps improve genome annotation strategies. Moreover, this is the first database to provide resources to broadly investigate functions and mechanisms involving TEs and ncRNAs in plants.

Recombination rate plasticity: revealing mechanisms by design

For over a century, scientists have known that meiotic recombination rates can vary considerably among individuals, and that environmental conditions can modify recombination rates relative to the background. A variety of external and intrinsic factors such as temperature, age, sex and starvation can elicit 'plastic' responses in recombination rate. The influence of recombination rate plasticity on genetic diversity of the next generation has interesting and important implications for how populations evolve. Further, many questions remain regarding the mechanisms and molecular processes that contribute to recombination rate plasticity. Here, we review 100 years of experimental work on recombination rate plasticity conducted in Drosophila melanogaster We categorize this work into four major classes of experimental designs, which we describe via classic studies in D. melanogaster Based on these studies, we highlight molecular mechanisms that are supported by experimental results and relate these findings to studies in other systems. We synthesize lessons learned from this model system into experimental guidelines for using recent advances in genotyping technologies, to study recombination rate plasticity in non-model organisms. Specifically, we recommend (1) using fine-scale genome-wide markers, (2) collecting time-course data, (3) including crossover distribution measurements, and (4) using mixed effects models to analyse results. To illustrate this approach, we present an application adhering to these guidelines from empirical work we conducted in Drosophila pseudoobscuraThis article is part of the themed issue 'Evolutionary causes and consequences of recombination rate variation in sexual organisms'.

Living Organisms Author Their Read-Write Genomes in Evolution

Evolutionary variations generating phenotypic adaptations and novel taxa resulted from complex cellular activities altering genome content and expression: (i) Symbiogenetic cell mergers producing the mitochondrion-bearing ancestor of eukaryotes and chloroplast-bearing ancestors of photosynthetic eukaryotes; (ii) interspecific hybridizations and genome doublings generating new species and adaptive radiations of higher plants and animals; and, (iii) interspecific horizontal DNA transfer encoding virtually all of the cellular functions between organisms and their viruses in all domains of life. Consequently, assuming that evolutionary processes occur in isolated genomes of individual species has become an unrealistic abstraction. Adaptive variations also involved natural genetic engineering of mobile DNA elements to rewire regulatory networks. In the most highly evolved organisms, biological complexity scales with "non-coding" DNA content more closely than with protein-coding capacity. Coincidentally, we have learned how so-called "non-coding" RNAs that are rich in repetitive mobile DNA sequences are key regulators of complex phenotypes. Both biotic and abiotic ecological challenges serve as triggers for episodes of elevated genome change. The intersections of cell activities, biosphere interactions, horizontal DNA transfers, and non-random Read-Write genome modifications by natural genetic engineering provide a rich molecular and biological foundation for understanding how ecological disruptions can stimulate productive, often abrupt, evolutionary transformations.

Selective silencing of euchromatic L1s revealed by genome-wide screens for L1 regulators

Transposable elements, also known as transposons, are now recognized not only as parasitic DNA, the spread of which in the genome must be controlled by the host, but also as major players in genome evolution and regulation. Long interspersed element-1 (LINE-1, also known as L1), the only currently autonomous mobile transposon in humans, occupies 17% of the genome and generates inter- and intra-individual genetic variation, in some cases resulting in disease. However, how L1 activity is controlled and the function of L1s in host gene regulation are not completely understood. Here we use CRISPR-Cas9 screening strategies in two distinct human cell lines to provide a genome-wide survey of genes involved in the control of L1 retrotransposition. We identify functionally diverse genes that either promote or restrict L1 retrotransposition. These genes, which are often associated with human diseases, control the L1 life cycle at the transcriptional or the post-transcriptional level in a manner that can depend on the endogenous L1 nucleotide sequence, underscoring the complexity of L1 regulation. We further investigate the restriction of L1 by the protein MORC2 and by the human silencing hub (HUSH) complex subunits MPP8 and TASOR. HUSH and MORC2 can selectively bind evolutionarily young, full-length L1s located within transcriptionally permissive euchromatic environments, and promote deposition of histone H3 Lys9 trimethylation (H3K9me3) for transcriptional silencing. Notably, these silencing events often occur within introns of transcriptionally active genes, and lead to the downregulation of host gene expression in a HUSH-, MORC2-, and L1-dependent manner. Together, these results provide a rich resource for studies of L1 retrotransposition, elucidate a novel L1 restriction pathway and illustrate how epigenetic silencing of transposable elements rewires host gene expression programs.

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.

A Horizontally Transferred Autonomous Helitron Became a Full Polydnavirus Segment in Cotesia vestalis

Bracoviruses associate symbiotically with thousands of parasitoid wasp species in the family Braconidae, working as virulence gene vectors, and allowing the development of wasp larvae within hosts. These viruses are composed of multiple DNA circles that are packaged into infective particles, and injected together with wasp's eggs during parasitization. One of the viral segments of Cotesia vestalis bracovirus contains a gene that has been previously described as a helicase of unknown origin. Here, we demonstrate that this gene is a Rep/Helicase from an intact Helitron transposable element that covers the viral segment almost entirely. We also provide evidence that this element underwent at least two horizontal transfers, which appear to have occurred consecutively: first from a Drosophila host ancestor to the genome of the parasitoid wasp C. vestalis and its bracovirus, and then from C. vestalis to a lepidopteran host (Bombyx mori). Our results reinforce the idea of parasitoid wasps as frequent agents of horizontal transfers in eukaryotes. Additionally, this Helitron-bracovirus segment is the first example of a transposable element that effectively became a whole viral circle.

Advances towards Controlling Meiotic Recombination for Plant Breeding

Meiotic homologous recombination generates new combinations of preexisting genetic variation and is a crucial process in plant breeding. Within the last decade, our understanding of plant meiotic recombination and genome diversity has advanced considerably. Innovation in DNA sequencing technology has led to the exploration of high-resolution genetic and epigenetic information in plant genomes, which has helped to accelerate plant breeding practices via high-throughput genotyping, and linkage and association mapping. In addition, great advances toward understanding the genetic and epigenetic control mechanisms of meiotic recombination have enabled the expansion of breeding programs and the unlocking of genetic diversity that can be used for crop improvement. This review highlights the recent literature on plant meiotic recombination and discusses the translation of this knowledge to the manipulation of meiotic recombination frequency and location with regards to crop plant breeding.

The Intra-Dependence of Viruses and the Holobiont

Animals live in symbiosis with the microorganisms surrounding them. This symbiosis is necessary for animal health, as a symbiotic breakdown can lead to a disease state. The functional symbiosis between the host, and associated prokaryotes, eukaryotes, and viruses in the context of an environment is the holobiont. Deciphering these holobiont associations has proven to be both difficult and controversial. In particular, holobiont association with viruses has been of debate even though these interactions have been occurring since cellular life began. The controversy stems from the idea that all viruses are parasitic, yet their associations can also be beneficial. To determine viral involvement within the holobiont, it is necessary to identify and elucidate the function of viral populations in symbiosis with the host. Viral metagenome analyses identify the communities of eukaryotic and prokaryotic viruses that functionally associate within a holobiont. Similarly, analyses of the host in response to viral presence determine how these interactions are maintained. Combined analyses reveal how viruses interact within the holobiont and how viral symbiotic cooperation occurs. To understand how the holobiont serves as a functional unit, one must consider viruses as an integral part of disease, development, and evolution.

Transposable Element Exaptation into Regulatory Regions Is Rare, Influenced by Evolutionary Age, and Subject to Pleiotropic Constraints

Transposable element (TE)-derived sequences make up approximately half of most mammalian genomes, and many TEs have been co-opted into gene regulatory elements. However, we lack a comprehensive tissue- and genome-wide understanding of how and when TEs gain regulatory activity in their hosts. We evaluated the prevalence of TE-derived DNA in enhancers and promoters across hundreds of human and mouse cell lines and primary tissues. Promoters are significantly depleted of TEs in all tissues compared with their overall prevalence in the genome (P < 0.001); enhancers are also depleted of TEs, though not as strongly as promoters. The degree of enhancer depletion also varies across contexts (1.5-3×), with reproductive and immune cells showing the highest levels of TE regulatory activity in humans. Overall, in spite of the regulatory potential of many TE sequences, they are significantly less active in gene regulation than expected from their prevalence. TE age is predictive of the likelihood of enhancer activity; TEs originating before the divergence of amniotes are 9.2 times more likely to have enhancer activity than TEs that integrated in great apes. Context-specific enhancers are more likely to be TE-derived than enhancers active in multiple tissues, and young TEs are more likely to overlap context-specific enhancers than old TEs (86% vs. 47%). Once TEs obtain enhancer activity in the host, they have similar functional dynamics to one another and non-TE-derived enhancers, likely driven by pleiotropic constraints. However, a few TE families, most notably endogenous retroviruses, have greater regulatory potential. Our observations suggest a model of regulatory co-option in which TE-derived sequences are initially repressed, after which a small fraction obtains context-specific enhancer activity, with further gains subject to pleiotropic constraints.

Two modes of targeting transposable elements by piRNA pathway in human testis

PIWI proteins and their partner small RNAs, termed piRNAs, are known to control transposable elements (TEs) in the germline. Here, we provide evidence that in humans this control is exerted in two different modes. On the one hand, production of piRNAs specifically targeting evolutionarily youngest TEs (L1HS, L1PA2-L1PA6, LTR12C, SVA) is present both at prenatal and postnatal stages of spermatogenesis and is performed without involvement of piRNA clusters. On the other hand, at postnatal stages, piRNAs deriving from pachytene clusters target "older" TEs and thus complement cluster-independent piRNA production to achieve relevant targeting of virtually all TEs expressed in postnatal testis. We also find that converging transcription of antisense-oriented genes contributes to the origin of genic postnatal prepachytene clusters. Finally, while a fraction of pachytene piRNAs was previously shown to arise from long intergenic noncoding RNAs (lincRNAs, i.e., pachytene piRNA cluster primary transcripts), we ascertain that these are a specific set of lincRNAs that both possess distinguishing epigenetic features and are expressed exclusively in testis.

Plasticity of the MFS1 Promoter Leads to Multidrug Resistance in the Wheat Pathogen Zymoseptoria tritici

The ascomycete Zymoseptoria tritici is the causal agent of Septoria leaf blotch on wheat. Disease control relies mainly on resistant wheat cultivars and on fungicide applications. The fungus displays a high potential to circumvent both methods. Resistance against all unisite fungicides has been observed over decades. A different type of resistance has emerged among wild populations with multidrug-resistant (MDR) strains. Active fungicide efflux through overexpression of the major facilitator gene MFS1 explains this emerging resistance mechanism. Applying a bulk-progeny sequencing approach, we identified in this study a 519-bp long terminal repeat (LTR) insert in the MFS1 promoter, a relic of a retrotransposon cosegregating with the MDR phenotype. Through gene replacement, we show the insert as a mutation responsible for MFS1 overexpression and the MDR phenotype. Besides this type I insert, we found two different types of promoter inserts in more recent MDR strains. Type I and type II inserts harbor potential transcription factor binding sites, but not the type III insert. Interestingly, all three inserts correspond to repeated elements present at different genomic locations in either IPO323 or other Z. tritici strains. These results underline the plasticity of repeated elements leading to fungicide resistance in Z. tritici and which contribute to its adaptive potential. IMPORTANCE Disease control through fungicides remains an important means to protect crops from fungal diseases and to secure the harvest. Plant-pathogenic fungi, especially Zymoseptoria tritici, have developed resistance against most currently used active ingredients, reducing or abolishing their efficacy. While target site modification is the most common resistance mechanism against single modes of action, active efflux of multiple drugs is an emerging phenomenon in fungal populations reducing additionally fungicides' efficacy in multidrug-resistant strains. We have investigated the mutations responsible for increased drug efflux in Z. tritici field strains. Our study reveals that three different insertions of repeated elements in the same promoter lead to multidrug resistance in Z. tritici. The target gene encodes the membrane transporter MFS1 responsible for drug efflux, with the promoter inserts inducing its overexpression. These results underline the plasticity of repeated elements leading to fungicide resistance in Z. tritici.

Chromosome Evolution in Connection with Repetitive Sequences and Epigenetics in Plants

Chromosome evolution is a fundamental aspect of evolutionary biology. The evolution of chromosome size, structure and shape, number, and the change in DNA composition suggest the high plasticity of nuclear genomes at the chromosomal level. Repetitive DNA sequences, which represent a conspicuous fraction of every eukaryotic genome, particularly in plants, are found to be tightly linked with plant chromosome evolution. Different classes of repetitive sequences have distinct distribution patterns on the chromosomes. Mounting evidence shows that repetitive sequences may play multiple generative roles in shaping the chromosome karyotypes in plants. Furthermore, recent development in our understanding of the repetitive sequences and plant chromosome evolution has elucidated the involvement of a spectrum of epigenetic modification. In this review, we focused on the recent evidence relating to the distribution pattern of repetitive sequences in plant chromosomes and highlighted their potential relevance to chromosome evolution in plants. We also discussed the possible connections between evolution and epigenetic alterations in chromosome structure and repatterning, such as heterochromatin formation, centromere function, and epigenetic-associated transposable element inactivation.

KRAB zinc finger proteins

Krüppel-associated box domain zinc finger proteins (KRAB-ZFPs) are the largest family of transcriptional regulators in higher vertebrates. Characterized by an N-terminal KRAB domain and a C-terminal array of DNA-binding zinc fingers, they participate, together with their co-factor KAP1 (also known as TRIM28), in repression of sequences derived from transposable elements (TEs). Until recently, KRAB-ZFP/KAP1-mediated repression of TEs was thought to lead to irreversible silencing, and the evolutionary selection of KRAB-ZFPs was considered to be just the host component of an arms race against TEs. However, recent advances indicate that KRAB-ZFPs and their TE targets also partner up to establish species-specific regulatory networks. Here, we provide an overview of the KRAB-ZFP gene family, highlighting how its evolutionary history is linked to that of TEs, and how KRAB-ZFPs influence multiple aspects of development and physiology.

Integrating networks and comparative genomics reveals retroelement proliferation dynamics in hominid genomes

Retroelements (REs) are mobile DNA sequences that multiply and spread throughout genomes by a copy-and-paste mechanism. These parasitic elements are active in diverse genomes, from yeast to humans, where they promote diversity, cause disease, and accelerate evolution. Because of their high copy number and sequence similarity, studying their activity and tracking their proliferation dynamics is a challenge. It is particularly difficult to pinpoint the few REs in a genome that are still active in the haystack of degenerate and suppressed elements. We develop a computational framework based on network theory that tracks the path of RE proliferation throughout evolution. We analyze SVA (SINE-VNTR-Alu), the youngest RE family in human genomes, to understand RE dynamics across hominids. Integrating comparative genomics and network tools enables us to track the course of SVA proliferation, identify yet unknown active communities, and detect tentative "master REs" that played key roles in SVA propagation, providing strong support for the fundamental "master gene" model of RE proliferation. The method is generic and thus can be applied to REs of any of the thousands of available genomes to identify active RE communities and master REs that were pivotal in the evolution of their host genomes.

Rapid Increase in Genome Size as a Consequence of Transposable Element Hyperactivity in Wood-White (Leptidea) Butterflies

Characterizing and quantifying genome size variation among organisms and understanding if genome size evolves as a consequence of adaptive or stochastic processes have been long-standing goals in evolutionary biology. Here, we investigate genome size variation and association with transposable elements (TEs) across lepidopteran lineages using a novel genome assembly of the common wood-white (Leptidea sinapis) and population re-sequencing data from both L. sinapis and the closely related L. reali and L. juvernica together with 12 previously available lepidopteran genome assemblies. A phylogenetic analysis confirms established relationships among species, but identifies previously unknown intraspecific structure within Leptidea lineages. The genome assembly of L. sinapis is one of the largest of any lepidopteran taxon so far (643 Mb) and genome size is correlated with abundance of TEs, both in Lepidoptera in general and within Leptidea where L. juvernica from Kazakhstan has considerably larger genome size than any other Leptidea population. Specific TE subclasses have been active in different Lepidoptera lineages with a pronounced expansion of predominantly LINEs, DNA elements, and unclassified TEs in the Leptidea lineage after the split from other Pieridae. The rate of genome expansion in Leptidea in general has been in the range of four Mb/Million year (My), with an increase in a particular L. juvernica population to 72 Mb/My. The considerable differences in accumulation rates of specific TE classes in different lineages indicate that TE activity plays a major role in genome size evolution in butterflies and moths.

Arc - An endogenous neuronal retrovirus?

The neuronal gene Arc is essential for long-lasting information storage in the mammalian brain and has been implicated in various neurological disorders. However, little is known about Arc's evolutionary origins. Recent studies suggest that mammalian Arc originated from a vertebrate lineage of Ty3/gypsy retrotransposons, which are also ancestral to retroviruses. In particular, Arc contains homology to the Gag polyprotein that forms the viral capsid and is essential for viral infectivity. This surprising connection raises the intriguing possibility that Arc may share molecular characteristics of retroviruses.

Analysis of the spotted gar genome suggests absence of causative link between ancestral genome duplication and transposable element diversification in teleost fish

Teleost fish have been shown to contain many superfamilies of transposable elements (TEs) that are absent from most tetrapod genomes. Since theories predict an increase in TE activity following polyploidization, such diversity might be linked to the 3R whole-genome duplication that occurred approximately 300 million years ago before the teleost radiation. To test this hypothesis, we have analyzed the genome of the spotted gar Lepisosteus oculatus, which diverged from the teleost lineage before the 3R duplication. Our results indicate that TE diversity and copy numbers are similar in gar and teleost genomes, suggesting that TE diversity was ancestral and not linked to the 3R whole-genome duplication. We propose that about 25 distinct superfamilies of TEs were present in the last ancestor of gars and teleost fish about 300 million years ago in the ray-finned fish lineage.

RNase H As Gene Modifier, Driver of Evolution and Antiviral Defense

Retroviral infections are 'mini-symbiotic' events supplying recipient cells with sequences for viral replication, including the reverse transcriptase (RT) and ribonuclease H (RNase H). These proteins and other viral or cellular sequences can provide novel cellular functions including immune defense mechanisms. Their high error rate renders RT-RNases H drivers of evolutionary innovation. Integrated retroviruses and the related transposable elements (TEs) have existed for at least 150 million years, constitute up to 80% of eukaryotic genomes and are also present in prokaryotes. Endogenous retroviruses regulate host genes, have provided novel genes including the syncytins that mediate maternal-fetal immune tolerance and can be experimentally rendered infectious again. The RT and the RNase H are among the most ancient and abundant protein folds. RNases H may have evolved from ribozymes, related to viroids, early in the RNA world, forming ribosomes, RNA replicases and polymerases. Basic RNA-binding peptides enhance ribozyme catalysis. RT and ribozymes or RNases H are present today in bacterial group II introns, the precedents of TEs. Thousands of unique RTs and RNases H are present in eukaryotes, bacteria, and viruses. These enzymes mediate viral and cellular replication and antiviral defense in eukaryotes and prokaryotes, splicing, R-loop resolvation, DNA repair. RNase H-like activities are also required for the activity of small regulatory RNAs. The retroviral replication components share striking similarities with the RNA-induced silencing complex (RISC), the prokaryotic CRISPR-Cas machinery, eukaryotic V(D)J recombination and interferon systems. Viruses supply antiviral defense tools to cellular organisms. TEs are the evolutionary origin of siRNA and miRNA genes that, through RISC, counteract detrimental activities of TEs and chromosomal instability. Moreover, piRNAs, implicated in transgenerational inheritance, suppress TEs in germ cells. Thus, virtually all known immune defense mechanisms against viruses, phages, TEs, and extracellular pathogens require RNase H-like enzymes. Analogous to the prokaryotic CRISPR-Cas anti-phage defense possibly originating from TEs termed casposons, endogenized retroviruses ERVs and amplified TEs can be regarded as related forms of inheritable immunity in eukaryotes. This survey suggests that RNase H-like activities of retroviruses, TEs, and phages, have built up innate and adaptive immune systems throughout all domains of life.

Ancient human miRNAs are more likely to have broad functions and disease associations than young miRNAs

Background

microRNAs (miRNAs) are essential to the regulation of gene expression in eukaryotes, and improper expression of miRNAs contributes to hundreds of diseases. Despite the essential functions of miRNAs, the evolutionary dynamics of how they are integrated into existing gene regulatory and functional networks is not well understood. Knowledge of the origin and evolutionary history a gene has proven informative about its functions and disease associations; we hypothesize that incorporating the evolutionary origins of miRNAs into analyses will help resolve differences in their functional dynamics and how they influence disease.

Results

We computed the phylogenetic age of miRNAs across 146 species and quantified the relationship between human miRNA age and several functional attributes. Older miRNAs are significantly more likely to be associated with disease than younger miRNAs, and the number of associated diseases increases with age. As has been observed for genes, the miRNAs associated with different diseases have different age profiles. For example, human miRNAs implicated in cancer are enriched for origins near the dawn of animal multicellularity. Consistent with the increasing contribution of miRNAs to disease with age, older miRNAs target more genes than younger miRNAs, and older miRNAs are expressed in significantly more tissues. Furthermore, miRNAs of all ages exhibit a strong preference to target older genes; 93% of validated miRNA gene targets were in existence at the origin of the targeting miRNA. Finally, we find that human miRNAs in evolutionarily related families are more similar in their targets and expression profiles than unrelated miRNAs.

Conclusions

Considering the evolutionary origin and history of a miRNA provides useful context for the analysis of its function. Consistent with recent work in Drosophila, our results support a model in which miRNAs increase their expression and functional regulatory interactions over evolutionary time, and thus older miRNAs have increased potential to cause disease. We anticipate that these patterns hold across mammalian species; however, comprehensively evaluating them will require refining miRNA annotations across species and collecting functional data in non-human systems.

Electronic supplementary material

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

Gapless genome assembly of Colletotrichum higginsianum reveals chromosome structure and association of transposable elements with secondary metabolite gene clusters

BACKGROUND

The ascomycete fungus Colletotrichum higginsianum causes anthracnose disease of brassica crops and the model plant Arabidopsis thaliana. Previous versions of the genome sequence were highly fragmented, causing errors in the prediction of protein-coding genes and preventing the analysis of repetitive sequences and genome architecture.

RESULTS

Here, we re-sequenced the genome using single-molecule real-time (SMRT) sequencing technology and, in combination with optical map data, this provided a gapless assembly of all twelve chromosomes except for the ribosomal DNA repeat cluster on chromosome 7. The more accurate gene annotation made possible by this new assembly revealed a large repertoire of secondary metabolism (SM) key genes (89) and putative biosynthetic pathways (77 SM gene clusters). The two mini-chromosomes differed from the ten core chromosomes in being repeat- and AT-rich and gene-poor but were significantly enriched with genes encoding putative secreted effector proteins. Transposable elements (TEs) were found to occupy 7% of the genome by length. Certain TE families showed a statistically significant association with effector genes and SM cluster genes and were transcriptionally active at particular stages of fungal development. All 24 subtelomeres were found to contain one of three highly-conserved repeat elements which, by providing sites for homologous recombination, were probably instrumental in four segmental duplications.

CONCLUSION

The gapless genome of C. higginsianum provides access to repeat-rich regions that were previously poorly assembled, notably the mini-chromosomes and subtelomeres, and allowed prediction of the complete SM gene repertoire. It also provides insights into the potential role of TEs in gene and genome evolution and host adaptation in this asexual pathogen.