Recurrent evolution of vertebrate transcription factors by transposase capture

Whole-genome sequencing reveals the genetic mechanisms of domestication in classical inbred mice

Background

The laboratory mouse was domesticated from the wild house mouse. Understanding the genetics underlying domestication in laboratory mice, especially in the widely used classical inbred mice, is vital for studies using mouse models. However, the genetic mechanism of laboratory mouse domestication remains unknown due to lack of adequate genomic sequences of wild mice.

Results

We analyze the genetic relationships by whole-genome resequencing of 36 wild mice and 36 inbred strains. All classical inbred mice cluster together distinctly from wild and wild-derived inbred mice. Using nucleotide diversity analysis, Fst, and XP-CLR, we identify 339 positively selected genes that are closely associated with nervous system function. Approximately one third of these positively selected genes are highly expressed in brain tissues, and genetic mouse models of 125 genes in the positively selected genes exhibit abnormal behavioral or nervous system phenotypes. These positively selected genes show a higher ratio of differential expression between wild and classical inbred mice compared with all genes, especially in the hippocampus and frontal lobe. Using a mutant mouse model, we find that the SNP rs27900929 (T>C) in gene Astn2 significantly reduces the tameness of mice and modifies the ratio of the two Astn2 (a/b) isoforms.

Conclusion

Our study indicates that classical inbred mice experienced high selection pressure during domestication under laboratory conditions. The analysis shows the positively selected genes are closely associated with behavior and the nervous system in mice. Tameness may be related to the Astn2 mutation and regulated by the ratio of the two Astn2 (a/b) isoforms.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02772-1.

Functional Characterization of the N-Terminal Disordered Region of the piggyBac Transposase

The piggyBac DNA transposon is an active element initially isolated from the cabbage looper moth, but members of this superfamily are also present in most eukaryotic evolutionary lineages. The functionally important regions of the transposase are well described. There is an RNase H-like fold containing the DDD motif responsible for the catalytic DNA cleavage and joining reactions and a C-terminal cysteine-rich domain important for interaction with the transposon DNA. However, the protein also contains a ~100 amino acid long N-terminal disordered region (NTDR) whose function is currently unknown. Here we show that deletion of the NTDR significantly impairs piggyBac transposition, although the extent of decrease is strongly cell-type specific. Moreover, replacing the NTDR with scrambled but similarly disordered sequences did not rescue transposase activity, indicating the importance of sequence conservation. Cell-based transposon excision and integration assays reveal that the excision step is more severely affected by NTDR deletion. Finally, bioinformatic analyses indicated that the NTDR is specific for the piggyBac superfamily and is also present in domesticated, transposase-derived proteins incapable of catalyzing transposition. Our results indicate an essential role of the NTDR in the "fine-tuning" of transposition and its significance in the functions of piggyBac-originated co-opted genes.

Transcriptional dynamics of transposable elements in the type I IFN response in Myotis lucifugus cells

Background

Bats are a major reservoir of zoonotic viruses, and there has been growing interest in characterizing bat-specific features of innate immunity and inflammation. Recent studies have revealed bat-specific adaptations affecting interferon (IFN) signaling and IFN-stimulated genes (ISGs), but we still have a limited understanding of the genetic mechanisms that have shaped the evolution of bat immunity. Here we investigated the transcriptional and epigenetic dynamics of transposable elements (TEs) during the type I IFN response in little brown bat (Myotis lucifugus) primary embryonic fibroblast cells, using RNA-seq and CUT&RUN.

Results

We found multiple bat-specific TEs that undergo both locus-specific and family-level transcriptional induction in response to IFN. Our transcriptome reassembly identified multiple ISGs that have acquired novel exons from bat-specific TEs, including NLRC5, SLNF5 and a previously unannotated isoform of the IFITM2 gene. We also identified examples of TE-derived regulatory elements, but did not find strong evidence supporting genome-wide epigenetic activation of TEs in response to IFN.

Conclusion

Collectively, our study uncovers numerous TE-derived transcripts, proteins, and alternative isoforms that are induced by IFN in Myotis lucifugus cells, highlighting candidate loci that may contribute to bat-specific immune function.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-022-00277-z.

Mammalian genome innovation through transposon domestication

Since the discovery of transposons, their sheer abundance in host genomes has puzzled many. While historically viewed as largely harmless 'parasitic' DNAs during evolution, transposons are not a mere record of ancient genome invasion. Instead, nearly every element of transposon biology has been integrated into host biology. Here we review how host genome sequences introduced by transposon activities provide raw material for genome innovation and document the distinct evolutionary path of each species.

New insights into the evolutionary dynamic and lineage divergence of gasdermin E in metazoa

Gasdermin (GSDM) is a family of pore-forming proteins that induce pyroptosis. To date, the origin and evolution of GSDM in Metazoa remain elusive. Here, we found that GSDM emerged early in Placozoa but is absent in a large number of invertebrates. In the lower vertebrate, fish, three types of GSDME, i.e., GSDMEa, GSDMEb, and a previously unreported type (designated GSDMEc), were idenitied. Evolutionarily, the three GSDMEs are distinctly separated: GSDMEa is closely related to tetrapod GSDME; GSDMEb exists exclusively in fish; GSDMEc forms the lineage root of tetrapod GSDMA/B/C/D. GSDMEc shares conserved genomic features with and is probably the prototype of GSDMA, which we found existing in all tetrapod classes. GSDMEc displays fast evolutionary dynamics, likely as a result of genomic transposition. A cross-metazoan analysis of GSDME revealed that GSDMEa shares a conserved caspase recognition motif with the GSDME of tetrapods and cnidarians, whereas GSDMEb has a unique caspase recognition motif similar to that of mammalian GSDMD, and GSDMEc exhibits no apparent caspase recognition motif. Through functional test, four highly conserved residues in vertebrate GSDME proved to be essential to auto-inhibition. Together our results provide new insights into the origin, evolution, and function of metazoan GSDMs.

Analysis of repeat elements in the Pristionchus pacificus genome reveals an ancient invasion by horizontally transferred transposons

BACKGROUND

Repetitive sequences and mobile elements make up considerable fractions of individual genomes. While transposition events can be detrimental for organismal fitness, repetitive sequences form an enormous reservoir for molecular innovation. In this study, we aim to add repetitive elements to the annotation of the Pristionchus pacificus genome and assess their impact on novel gene formation.

RESULTS

Different computational approaches define up to 24% of the P. pacificus genome as repetitive sequences. While retroelements are more frequently found at the chromosome arms, DNA transposons are distributed more evenly. We found multiple DNA transposons, as well as LTR and LINE elements with abundant evidence of expression as single-exon transcripts. When testing whether transposons disproportionately contribute towards new gene formation, we found that roughly 10-20% of genes across all age classes overlap transposable elements with the strongest trend being an enrichment of low complexity regions among the oldest genes. Finally, we characterized a horizontal gene transfer of Zisupton elements into diplogastrid nematodes. These DNA transposons invaded nematodes from eukaryotic donor species and experienced a recent burst of activity in the P. pacificus lineage.

CONCLUSIONS

The comprehensive annotation of repetitive elements in the P. pacificus genome builds a resource for future functional genomic analyses as well as for more detailed investigations of molecular innovations.

Roles of transposable elements in the regulation of mammalian transcription

Transposable elements (TEs) comprise about half of the mammalian genome. TEs often contain sequences capable of recruiting the host transcription machinery, which they use to express their own products and promote transposition. However, the regulatory sequences carried by TEs may affect host transcription long after the TEs have lost the ability to transpose. Recent advances in genome analysis and engineering have facilitated systematic interrogation of the regulatory activities of TEs. In this Review, we discuss diverse mechanisms by which TEs contribute to transcription regulation. Notably, TEs can donate enhancer and promoter sequences that influence the expression of host genes, modify 3D chromatin architecture and give rise to novel regulatory genes, including non-coding RNAs and transcription factors. We discuss how TEs spur regulatory evolution and facilitate the emergence of genetic novelties in mammalian physiology and development. By virtue of their repetitive and interspersed nature, TEs offer unique opportunities to dissect the effects of mutation and genomic context on the function and evolution of cis-regulatory elements. We argue that TE-centric studies hold the key to unlocking general principles of transcription regulation and evolution.

Zebrafish transposable elements show extensive diversification in age, genomic distribution, and developmental expression

There is considerable interest in understanding the effect of transposable elements (TEs) on embryonic development. Studies in humans and mice are limited by the difficulty of working with mammalian embryos and by the relative scarcity of active TEs in these organisms. The zebrafish is an outstanding model for the study of vertebrate development, and over half of its genome consists of diverse TEs. However, zebrafish TEs remain poorly characterized. Here we describe the demography and genomic distribution of zebrafish TEs and their expression throughout embryogenesis using bulk and single-cell RNA sequencing data. These results reveal a highly dynamic genomic ecosystem comprising nearly 2000 distinct TE families, which vary in copy number by four orders of magnitude and span a wide range of ages. Longer retroelements tend to be retained in intergenic regions, whereas short interspersed nuclear elements (SINEs) and DNA transposons are more frequently found nearby or within genes. Locus-specific mapping of TE expression reveals extensive TE transcription during development. Although two-thirds of TE transcripts are likely driven by nearby gene promoters, we still observe stage- and tissue-specific expression patterns in self-regulated TEs. Long terminal repeat (LTR) retroelements are most transcriptionally active immediately following zygotic genome activation, whereas DNA transposons are enriched among transcripts expressed in later stages of development. Single-cell analysis reveals several endogenous retroviruses expressed in specific somatic cell lineages. Overall, our study provides a valuable resource for using zebrafish as a model to study the impact of TEs on vertebrate development.

Neofunctionalization of a non-coding portion of a DNA transposon in the coding region of the chimerical sex-determining gene dm-W in Xenopus frogs

Most vertebrate sex-determining genes (SDGs) emerge as neofunctionalized genes through duplication and/or mutation of ancestral genes that are involved with sexual differentiation. We previously demonstrated dm-W to be the SDG in the African clawed frog Xenopus laevis and found that a portion of this gene emerged from the masculinization gene dmrt1 after allotetraploidization by interspecific hybridization between two ancestral species around 17–18 Ma. dm-W has four exons consisting of a noncoding exon 1, dmrt1-derived exons 2 and 3, and an orphan exon 4 (Ex4) of unknown origin that includes coding sequence (CDS). In this study, we searched for the origin of Ex4 and investigated the function of the CDS of this exon. We found that the Ex4-CDS is derived from a noncoding portion of the hAT-10 family of DNA transposon. Evolutionary analysis of transposons and determination of the Ex4 sequences from three other species indicated that Ex4 was generated before the diversification of most or all extant allotetraploid species in subgenus Xenopus, during which time we hypothesize that transposase activity of this hAT superfamily was active. Using DNA–protein binding and transfection assays, we further demonstrate that the Ex4-encoded amino acid sequence increases the DNA-binding ability and transrepression activity of DM-W. These findings suggest that the conversion of the noncoding transposon sequence to the CDS of dm-W contributed to neofunctionalization of a new chimeric SDG in the ancestor of the allotetraploid Xenopus species, offering new insights into de novo origin and functional evolution of chimerical genes.

A phylogenetic and proteomic reconstruction of eukaryotic chromatin evolution

Histones and associated chromatin proteins have essential functions in eukaryotic genome organization and regulation. Despite this fundamental role in eukaryotic cell biology, we lack a phylogenetically-comprehensive understanding of chromatin evolution. Here, we combine comparative proteomics and genomics analysis of chromatin in eukaryotes and archaea. Proteomics uncovers the existence of histone post-translational modifications in Archaea. However, archaeal histone modifications are scarce, in contrast with the highly conserved and abundant marks we identify across eukaryotes. Phylogenetic analysis reveals that chromatin-associated catalytic functions (e.g., methyltransferases) have pre-eukaryotic origins, whereas histone mark readers and chaperones are eukaryotic innovations. We show that further chromatin evolution is characterized by expansion of readers, including capture by transposable elements and viruses. Overall, our study infers detailed evolutionary history of eukaryotic chromatin: from its archaeal roots, through the emergence of nucleosome-based regulation in the eukaryotic ancestor, to the diversification of chromatin regulators and their hijacking by genomic parasites.

Structural insights into the evolution of the RAG recombinase

Adaptive immunity in jawed vertebrates relies on the assembly of antigen receptor genes by the recombination activating gene 1 (RAG1)-RAG2 (collectively RAG) recombinase in a reaction known as V(D)J recombination. Extensive biochemical and structural evidence indicates that RAG and V(D)J recombination evolved from the components of a RAG-like (RAGL) transposable element through a process known as transposon molecular domestication. This Review describes recent advances in our understanding of the functional and structural transitions that occurred during RAG evolution. We use the structures of RAG and RAGL enzymes to trace the evolutionary adaptations that yielded a RAG recombinase with exquisitely regulated cleavage activity and a multilayered array of mechanisms to suppress transposition. We describe how changes in modes of DNA binding, alterations in the dynamics of protein-DNA complexes, single amino acid mutations and a modular design likely enabled RAG family enzymes to survive and spread in the genomes of eukaryotes. These advances highlight the insight that can be gained from viewing evolution of vertebrate immunity through the lens of comparative genome analyses coupled with structural biology and biochemistry.

Intraspecific de novo gene birth revealed by presence-absence variant genes in Caenorhabditis elegans

Genes embed their evolutionary history in the form of various alleles. Presence-absence variants (PAVs) are extreme cases of such alleles, where a gene present in one haplotype does not exist in another. Because PAVs may result from either birth or death of a gene, PAV genes and their alternative alleles, if available, can represent a basis for rapid intraspecific gene evolution. Using long-read sequencing technologies, this study traced the possible evolution of PAV genes in the PD1074 and CB4856 C. elegans strains as well as their alternative alleles in 14 other wild strains. We updated the CB4856 genome by filling 18 gaps and identified 46 genes and 7,460 isoforms from both strains not annotated previously. We verified 328 PAV genes, out of which 46 were C. elegans-specific. Among these possible newly born genes, 12 had alternative alleles in other wild strains; in particular, the alternative alleles of three genes showed signatures of active transposons. Alternative alleles of three other genes showed another type of signature reflected in accumulation of small insertions or deletions. Research on gene evolution using both species-specific PAV genes and their alternative alleles may provide new insights into the process of gene evolution.

Diversity of Harbinger-like Transposons in Teleost Fish Genomes

Simple Summary

The study of transposable elements, which are repeated DNA sequences that can insert into new locations in genomes, is of particular interest to genome evolution, as they are sources of mutations but also of new regulatory and coding sequences. Teleost fish are a species-rich clade presenting a high diversity of transposable elements, both quantitatively and qualitatively, making them a very attractive group to investigate the evolution of mobile sequences. We studied Harbinger-like DNA transposons, which are widespread from plants to vertebrates but absent from mammalian genomes. These elements code for both a transposase and a Myb-like protein. We observed high variability in the genomic composition of Harbinger-like sequences in teleost fish. While Harbinger transposons might have been present in a common ancestor of all the fish species studied, ISL2EU elements were possibly gained by horizontal transfer at the base of teleost fish. Transposase and Myb-like protein phylogenies of Harbinger transposons indicated unique origins of the association between both genes and suggests recombination was rare between transposon sublineages. Finally, we report one case of Harbinger horizontal transfer between divergent fish species and the transcriptional activity of both Harbinger and ISL2EU transposons in teleost fish. There was male-biased expression in the gonads of the medaka fish.

Abstract

Harbinger elements are DNA transposons that are widespread from plants to vertebrates but absent from mammalian genomes. Among vertebrates, teleost fish are the clade presenting not only the largest number of species but also the highest diversity of transposable elements, both quantitatively and qualitatively, making them a very attractive group to investigate the evolution of mobile sequences. We studied Harbinger DNA transposons and the distantly related ISL2EU elements in fish, focusing on representative teleost species compared to the spotted gar, the coelacanth, the elephant shark and the amphioxus. We observed high variability in the genomic composition of Harbinger-like sequences in teleost fish, as they covered 0.002–0.14% of the genome, when present. While Harbinger transposons might have been present in a common ancestor of all the fish species studied here, with secondary loss in elephant shark, our results suggests that ISL2EU elements were gained by horizontal transfer at the base of teleost fish 200–300 million years ago, and that there was secondary loss in a common ancestor of pufferfishes and stickleback. Harbinger transposons code for a transposase and a Myb-like protein. We reconstructed and compared molecular phylogenies of both proteins to get insights into the evolution of Harbinger transposons in fish. Transposase and Myb-like protein phylogenies showed global congruent evolution, indicating unique origin of the association between both genes and suggesting rare recombination between transposon sublineages. Finally, we report one case of Harbinger horizontal transfer between divergent fish species and the transcriptional activity of both Harbinger and ISL2EU transposons in teleost fish. There was male-biased expression in the gonads of the medaka fish.

A retrotransposon storm marks clinical phenoconversion to late-onset Alzheimer's disease

Recent reports have suggested that the reactivation of otherwise transcriptionally silent transposable elements (TEs) might induce brain degeneration, either by dysregulating the expression of genes and pathways implicated in cognitive decline and dementia or through the induction of immune-mediated neuroinflammation resulting in the elimination of neural and glial cells. In the work we present here, we test the hypothesis that differentially expressed TEs in blood could be used as biomarkers of cognitive decline and development of AD. To this aim, we used a sample of aging subjects (age > 70) that developed late-onset Alzheimer’s disease (LOAD) over a relatively short period of time (12–48 months), for which blood was available before and after their phenoconversion, and a group of cognitive stable subjects as controls. We applied our developed and validated customized pipeline that allows the identification, characterization, and quantification of the differentially expressed (DE) TEs before and after the onset of manifest LOAD, through analyses of RNA-Seq data. We compared the level of DE TEs within more than 600,000 TE-mapping RNA transcripts from 25 individuals, whose specimens we obtained before and after their phenotypic conversion (phenoconversion) to LOAD, and discovered that 1790 TE transcripts showed significant expression differences between these two timepoints (logFC ± 1.5, logCMP > 5.3, nominal p value < 0.01). These DE transcripts mapped both over- and under-expressed TE elements. Occurring before the clinical phenoconversion, this TE storm features significant increases in DE transcripts of LINEs, LTRs, and SVAs, while those for SINEs are significantly depleted. These dysregulations end with signs of manifest LOAD. This set of highly DE transcripts generates a TE transcriptional profile that accurately discriminates the before and after phenoconversion states of these subjects. Our findings suggest that a storm of DE TEs occurs before phenoconversion from normal cognition to manifest LOAD in risk individuals compared to controls, and may provide useful blood-based biomarkers for heralding such a clinical transition, also suggesting that TEs can indeed participate in the complex process of neurodegeneration.

The state of Medusozoa genomics: current evidence and future challenges

Medusozoa is a widely distributed ancient lineage that harbors one-third of Cnidaria diversity divided into 4 classes. This clade is characterized by the succession of stages and modes of reproduction during metagenic lifecycles, and includes some of the most plastic body plans and life cycles among animals. The characterization of traditional genomic features, such as chromosome numbers and genome sizes, was rather overlooked in Medusozoa and many evolutionary questions still remain unanswered. Modern genomic DNA sequencing in this group started in 2010 with the publication of the Hydra vulgaris genome and has experienced an exponential increase in the past 3 years. Therefore, an update of the state of Medusozoa genomics is warranted. We reviewed different sources of evidence, including cytogenetic records and high-throughput sequencing projects. We focused on 4 main topics that would be relevant for the broad Cnidaria research community: (i) taxonomic coverage of genomic information; (ii) continuity, quality, and completeness of high-throughput sequencing datasets; (iii) overview of the Medusozoa specific research questions approached with genomics; and (iv) the accessibility of data and metadata. We highlight a lack of standardization in genomic projects and their reports, and reinforce a series of recommendations to enhance future collaborative research.

Dynamics and Impacts of Transposable Element Proliferation in the Drosophila nasuta Species Group Radiation

Transposable element (TE) mobilization is a constant threat to genome integrity. Eukaryotic organisms have evolved robust defensive mechanisms to suppress their activity, yet TEs can escape suppression and proliferate, creating strong selective pressure for host defense to adapt. This genomic conflict fuels a never-ending arms race that drives the rapid evolution of TEs and recurrent positive selection of genes involved in host defense; the latter has been shown to contribute to postzygotic hybrid incompatibility. However, how TE proliferation impacts genome and regulatory divergence remains poorly understood. Here, we report the highly complete and contiguous (N50 = 33.8-38.0 Mb) genome assemblies of seven closely related Drosophila species that belong to the nasuta species group-a poorly studied group of flies that radiated in the last 2 My. We constructed a high-quality de novo TE library and gathered germline RNA-seq data, which allowed us to comprehensively annotate and compare TE insertion patterns between the species, and infer the evolutionary forces controlling their spread. We find a strong negative association between TE insertion frequency and expression of genes nearby; this likely reflects survivor bias from reduced fitness impact of TEs inserting near lowly expressed, nonessential genes, with limited TE-induced epigenetic silencing. Phylogenetic analyses of insertions of 147 TE families reveal that 53% of them show recent amplification in at least one species. The most highly amplified TE is a nonautonomous DNA element (Drosophila INterspersed Element; DINE) which has gone through multiple bouts of expansions with thousands of full-length copies littered throughout each genome. Across all TEs, we find that TEs expansions are significantly associated with high expression in the expanded species consistent with suppression escape. Thus, whereas horizontal transfer followed by the invasion of a naïve genome has been highlighted to explain the long-term survival of TEs, our analysis suggests that evasion of host suppression of resident TEs is a major strategy to persist over evolutionary times. Altogether, our results shed light on the heterogenous and context-dependent nature in which TEs affect gene regulation and the dynamics of rampant TE proliferation amidst a recently radiated species group.

icaR and icaT Are Ancient Chromosome Genes Encoding Substrates of the Type III Secretion Apparatus in Shigella flexneri

Shigella is an Escherichia coli pathovar that colonizes the cytosol of mucosal cells in the human large intestine. To do this, Shigella uses a Type III Secretion Apparatus (T3SA) to translocate several proteins into host cells. The T3SA and its substrates are encoded by genes of the virulence plasmid pINV or by chromosomal genes derived thereof. We recently discovered two chromosomal genes, which seem unrelated to pINV, although they are activated by MxiE and IpgC similarly to some of the canonical substrates of the T3SA. Here, we showed that the production of the corresponding proteins depended on the conservation of a MxiE box in their cognate promoters. Furthermore, both proteins were secreted by the T3SA in a chaperone-independent manner through the recognition of their respective amino-terminal secretion signal. Based on these observations, we named these new genes icaR and icaT, which stand for invasion chromosome antigen with homology for a transcriptional regulator and a transposase, respectively. icaR and icaT have orthologs in commensal and pathogenic E. coli strains belonging mainly to phylogroups A, B1, D and E. Finally, we demonstrated that icaR and icaT orthologs could be activated by the coproduction of IpgC and MxiE in strains MG1655 K-12 (phylogroup A) and O157:H7 ATCC 43888 (phylogroup E). In contrast, the coproduction of EivF and YgeG, which are homologs of MxiE and IpgC in the E. coli T3SS 2 (ETT2), failed to activate icaR and icaT.

IMPORTANCEicaR and icaT are the latest members of the MxiE regulon discovered in the chromosome. The proteins IcaR and IcaT, albeit produced in small amounts, are nonetheless secreted by the T3SA comparably to canonical substrates. The high occurrence of icaR and icaT in phylogroups A, B1, D, and E coupled with their widespread absence in their B2 counterparts agree with the consensus E. coli phylogeny. The widespread conservation of the MxiE box among icaR and icaT orthologs supports the notion that both genes had already undergone coevolution with transcriptional activators ipgC and mxiE- harbored in pINV or a relative- in the last common ancestor of Shigella and of E. coli from phylogroups A, B1, D, and E. The possibility that icaR and icaT may contribute to Shigella pathogenesis cannot be excluded, although some of their characteristics suggest they are fossil genes.

Transposable Elements: Major Players in Shaping Genomic and Evolutionary Patterns

Transposable elements (TEs) are ubiquitous genetic elements, able to jump from one location of the genome to another, in all organisms. For this reason, on the one hand, TEs can induce deleterious mutations, causing dysfunction, disease and even lethality in individuals. On the other hand, TEs can increase genetic variability, making populations better equipped to respond adaptively to environmental change. To counteract the deleterious effects of TEs, organisms have evolved strategies to avoid their activation. However, their mobilization does occur. Usually, TEs are maintained silent through several mechanisms, but they can be reactivated during certain developmental windows. Moreover, TEs can become de-repressed because of drastic changes in the external environment. Here, we describe the ‘double life’ of TEs, being both ‘parasites’ and ‘symbionts’ of the genome. We also argue that the transposition of TEs contributes to two important evolutionary processes: the temporal dynamic of evolution and the induction of genetic variability. Finally, we discuss how the interplay between two TE-dependent phenomena, insertional mutagenesis and epigenetic plasticity, plays a role in the process of evolution.

The evolution and role of eukaryotic-like domains in environmental intracellular bacteria: the battle with a eukaryotic cell

Intracellular pathogens that are able to thrive in different environments, such as Legionella spp. which preferentially live in protozoa in aquatic environments or environmental Chlamydiae which replicate either within protozoa or a range of animals, possess a plethora of cellular biology tools to influence their eukaryotic host. The host manipulation tools that evolved in the interaction with protozoa, confer these bacteria the capacity to also infect phylogenetically distinct eukaryotic cells, such as macrophages and thus they can also be human pathogens. To manipulate the host cell, bacteria use protein secretion systems and molecular effectors. Although these molecular effectors are encoded in bacteria, they are expressed and function in a eukaryotic context often mimicking or inhibiting eukaryotic proteins. Indeed, many of these effectors have eukaryotic-like domains. In this review we propose that the main pathways environmental intracellular bacteria need to subvert in order to establish the host eukaryotic cell as a replication niche are chromatin remodelling, ubiquitination signalling, and modulation of protein-protein interactions via tandem repeat domains. We then provide mechanistic insight into how these proteins might have evolved as molecular weapons. Finally, we highlight that in environmental intracellular bacteria the number of eukaryotic-like domains and proteins is considerably higher than in intracellular bacteria specialised to an isolated niche, such as obligate intracellular human pathogens. As mimics of eukaryotic proteins are critical components of host pathogen interactions, this distribution of eukaryotic-like domains suggests that the environment has selected them.

Taming, Domestication and Exaptation: Trajectories of Transposable Elements in Genomes

During evolution, several types of sequences pass through genomes. Along with mutations and internal genetic tinkering, they are a useful source of genetic variability for adaptation and evolution. Most of these sequences are acquired by horizontal transfers (HT), but some of them may come from the genomes themselves. If they are not lost or eliminated quickly, they can be tamed, domesticated, or even exapted. Each of these processes results from a series of events, depending on the interactions between these sequences and the host genomes, but also on environmental constraints, through their impact on individuals or population fitness. After a brief reminder of the characteristics of each of these states (taming, domestication, exaptation), the evolutionary trajectories of these new or acquired sequences will be presented and discussed, emphasizing that they are not totally independent insofar as the first can constitute a step towards the second, and the second is another step towards the third.

Genome Stability Is in the Eye of the Beholder: CR1 Retrotransposon Activity Varies Significantly across Avian Diversity

Since the sequencing of the zebra finch genome it has become clear that avian genomes, while largely stable in terms of chromosome number and gene synteny, are more dynamic at an intrachromosomal level. A multitude of intrachromosomal rearrangements and significant variation in transposable element (TE) content have been noted across the avian tree. TEs are a source of genome plasticity, because their high similarity enables chromosomal rearrangements through nonallelic homologous recombination, and they have potential for exaptation as regulatory and coding sequences. Previous studies have investigated the activity of the dominant TE in birds, chicken repeat 1 (CR1) retrotransposons, either focusing on their expansion within single orders, or comparing passerines with nonpasserines. Here, we comprehensively investigate and compare the activity of CR1 expansion across orders of birds, finding levels of CR1 activity vary significantly both between and within orders. We describe high levels of TE expansion in genera which have speciated in the last 10 Myr including kiwis, geese, and Amazon parrots; low levels of TE expansion in songbirds across their diversification, and near inactivity of TEs in the cassowary and emu for millions of years. CR1s have remained active over long periods of time across most orders of neognaths, with activity at any one time dominated by one or two families of CR1s. Our findings of higher TE activity in species-rich clades and dominant families of TEs within lineages mirror past findings in mammals and indicate that genome evolution in amniotes relies on universal TE-driven processes.

Characterization of DNA methylation and promoter activity of long terminal repeat elements of feline endogenous retrovirus RDRS C2a

Endogenous retroviruses (ERVs) are genomic elements derived from retroviral infections in ancestral germ lines. Most ERVs are inactivated by genetic or epigenetic mechanisms, such as DNA methylation. RD-114-virus-related sequence (RDRS) C2a is a feline endogenous retrovirus present in all domestic cats; however, its expression and function are not clearly known. DNA methylation at CpG dinucleotides is a hallmark of silenced ERVs. This study aimed to investigate whether long terminal repeats (LTRs) of RDRS C2a function as a gene regulatory region. The DNA methylation status of RDRS C2a was examined by bisulfite sequencing, and CpG sites in 5' LTR of RDRS C2a were found hypomethylated, whereas those in 3' LTR were hypermethylated in feline cells. Several transcription factor-binding sites were identified in LTRs of RDRS C2a. Luciferase assay suggested that 5' LTR of RDRS C2a exhibited strong transcriptional activity, which was suppressed by in vitro DNA methylation. The study indicates that 5' LTR of RDRS C2a possibly functions as a promoter for itself or neighboring genes.

Transposable elements shape the evolution of mammalian development

Transposable elements (TEs) promote genetic innovation but also threaten genome stability. Despite multiple layers of host defence, TEs actively shape mammalian-specific developmental processes, particularly during pre-implantation and extra-embryonic development and at the maternal-fetal interface. Here, we review how TEs influence mammalian genomes both directly by providing the raw material for genetic change and indirectly via co-evolving TE-binding Krüppel-associated box zinc finger proteins (KRAB-ZFPs). Throughout mammalian evolution, individual activities of ancient TEs were co-opted to enable invasive placentation that characterizes live-born mammals. By contrast, the widespread activity of evolutionarily young TEs may reflect an ongoing co-evolution that continues to impact mammalian development.

The Evolutionary Volte-Face of Transposable Elements: From Harmful Jumping Genes to Major Drivers of Genetic Innovation

Transposable elements (TEs) are self-replicating DNA elements that constitute major fractions of eukaryote genomes. Their ability to transpose can modify the genome structure with potentially deleterious effects. To repress TE activity, host cells have developed numerous strategies, including epigenetic pathways, such as DNA methylation or histone modifications. Although TE neo-insertions are mostly deleterious or neutral, they can become advantageous for the host under specific circumstances. The phenomenon leading to the appropriation of TE-derived sequences by the host is known as TE exaptation or co-option. TE exaptation can be of different natures, through the production of coding or non-coding DNA sequences with ultimately an adaptive benefit for the host. In this review, we first give new insights into the silencing pathways controlling TE activity. We then discuss a model to explain how, under specific environmental conditions, TEs are unleashed, leading to a TE burst and neo-insertions, with potential benefits for the host. Finally, we review our current knowledge of coding and non-coding TE exaptation by providing several examples in various organisms and describing a method to identify TE co-option events.

Reconstruction of human genome evolution in yeast: an educational primer for use with "systematic humanization of the yeast cytoskeleton discerns functionally replaceable from divergent human genes"

The evolution of eukaryotic organisms starting with the last eukaryotic common ancestor was accompanied by lineage-specific expansion of gene families. A paper by Garge et al. provides an excellent opportunity to have students explore how expansion of gene families via gene duplication results in protein specialization, in this case in the context of eukaryotic cytoskeletal organization . The authors tested hypotheses about conserved protein function by systematic "humanization" of the yeast cytoskeletal components while employing a wide variety of methodological approaches. We outline several exercises to promote students' ability to explore the genomic databases, perform bioinformatic analyses, design experiments for functional analysis of human genes in yeast and critically interpret results to address both specific and general questions.

The Landscapes of Full-Length Transcripts and Splice Isoforms as Well as Transposons Exonization in the Lepidopteran Model System, Bombyx mori

The domesticated silkworm, Bombyx mori, is an important model system for the order Lepidoptera. Currently, based on third-generation sequencing, the chromosome-level genome of Bombyx mori has been released. However, its transcripts were mainly assembled by using short reads of second-generation sequencing and expressed sequence tags which cannot explain the transcript profile accurately. Here, we used PacBio Iso-Seq technology to investigate the transcripts from 45 developmental stages of Bombyx mori. We obtained 25,970 non-redundant high-quality consensus isoforms capturing ∼60% of previous reported RNAs, 15,431 (∼47%) novel transcripts, and identified 7,253 long non-coding RNA (lncRNA) with a large proportion of novel lncRNA (∼56%). In addition, we found that transposable elements (TEs) exonization account for 11,671 (∼45%) transcripts including 5,980 protein-coding transcripts (∼32%) and 5,691 lncRNAs (∼79%). Overall, our results expand the silkworm transcripts and have general implications to understand the interaction between TEs and their host genes. These transcripts resource will promote functional studies of genes and lncRNAs as well as TEs in the silkworm.

ExOrthist: a tool to infer exon orthologies at any evolutionary distance

Several bioinformatic tools have been developed for genome-wide identification of orthologous and paralogous genes. However, no corresponding tool allows the detection of exon homology relationships. Here, we present ExOrthist, a fully reproducible Nextflow-based software enabling inference of exon homologs and orthogroups, visualization of evolution of exon-intron structures, and assessment of conservation of alternative splicing patterns. ExOrthist evaluates exon sequence conservation and considers the surrounding exon-intron context to derive genome-wide multi-species exon homologies at any evolutionary distance. We demonstrate its use in different evolutionary scenarios: whole genome duplication in frogs and convergence of Nova-regulated splicing networks (https://github.com/biocorecrg/ExOrthist).

A domesticated Harbinger transposase forms a complex with HDA6 and promotes histone H3 deacetylation at genes but not TEs in Arabidopsis

In eukaryotes, histone acetylation is a major modification on histone N-terminal tails that is tightly connected to transcriptional activation. HDA6 is a histone deacetylase involved in the transcriptional regulation of genes and transposable elements (TEs) in Arabidopsis thaliana. HDA6 has been shown to participate in several complexes in plants, including a conserved SIN3 complex. Here, we uncover a novel protein complex containing HDA6, several Harbinger transposon-derived proteins (HHP1, SANT1, SANT2, SANT3, and SANT4), and MBD domain-containing proteins (MBD1, MBD2, and MBD4). We show that mutations of all four SANT genes in the sant-null mutant cause increased expression of the flowering repressors FLC, MAF4, and MAF5, resulting in a late flowering phenotype. Transcriptome deep sequencing reveals that while the SANT proteins and HDA6 regulate the expression of largely overlapping sets of genes, TE silencing is unaffected in sant-null mutants. Our global histone H3 acetylation profiling shows that SANT proteins and HDA6 modulate gene expression through deacetylation. Collectively, our findings suggest that Harbinger transposon-derived SANT domain-containing proteins are required for histone deacetylation and flowering time control in plants.

Genome-wide mapping of binding sites of the transposase-derived SETMAR protein in the human genome

Throughout evolution, DNA transposons provide a recurrent supply of genetic information to give rise to novel gene functions by fusion of their transposase domain to various domains of host-encoded proteins. One of these "domesticated", transposase-derived factors is SETMAR/Metnase which is a naturally occurring fusion protein that consists of a histone-lysine methyltransferase domain and an HsMar1 transposase. To elucidate the biological role of SETMAR, it is crucial to identify genomic targets to which SETMAR specifically binds and link these sites to the regulation of gene expression. Herein, we mapped the genomic landscape of SETMAR binding in a near-haploid human leukemia cell line (HAP1) in order to identify on-target and off-target binding sites at high resolution and to elucidate their role in terms of gene expression. Our analysis revealed a perfect correlation between SETMAR and inverted terminal repeats (ITRs) of HsMar1 transposon remnants, which are considered as natural target sites for SETMAR binding. However, we did not detect any untargeted events at non-ITR sequences, calling into question previously proposed off-target binding sites. We identified sequence fidelity of the ITR motif as a key factor for determining the binding affinity of SETMAR for chromosomes, as higher conservation of ITR sequences resulted in increased affinity for chromatin and stronger repression of SETMAR-bound gene loci. These associations highlight how SETMAR's chromatin binding fine-tune gene regulatory networks in human tumour cells.

New Genes Born-In or Invading Vertebrate Genomes

Which is the origin of genes is a fundamental question in Biology, indeed a question older than the discovery of genes itself. For more than a century, it was uneven to think in origins other than duplication and divergence from a previous gene. In recent years, however, the intersection of genetics, embryonic development, and bioinformatics, has brought to light that de novo generation from non-genic DNA, horizontal gene transfer and, noticeably, virus and transposon invasions, have shaped current genomes, by integrating those newcomers into old gene networks, helping to shape morphological and physiological innovations. We here summarized some of the recent research in the field, mostly in the vertebrate lineage with a focus on protein-coding novelties, showing that the placenta, the adaptative immune system, or the highly developed neocortex, among other innovations, are linked to de novo gene creation or domestication of virus and transposons. We provocatively suggest that the high tolerance to virus infections by bats may also be related to previous virus and transposon invasions in the bat lineage.

Transposable Elements Contribute to Genome Dynamics and Gene Expression Variation in the Fungal Plant Pathogen Verticillium dahliae

Transposable elements (TEs) are a major source of genetic and regulatory variation in their host genome and are consequently thought to play important roles in evolution. Many fungal and oomycete plant pathogens have evolved dynamic and TE-rich genomic regions containing genes that are implicated in host colonization and adaptation. TEs embedded in these regions have typically been thought to accelerate the evolution of these genomic compartments, but little is known about their dynamics in strains that harbor them. Here, we used whole-genome sequencing data of 42 strains of the fungal plant pathogen Verticillium dahliae to systematically identify polymorphic TEs that may be implicated in genomic as well as in gene expression variation. We identified 2,523 TE polymorphisms and characterize a subset of 8% of the TEs as polymorphic elements that are evolutionary younger, less methylated, and more highly expressed when compared with the remaining 92% of the total TE complement. As expected, the polyrmorphic TEs are enriched in the adaptive genomic regions. Besides, we observed an association of polymorphic TEs with pathogenicity-related genes that localize nearby and that display high expression levels. Collectively, our analyses demonstrate that TE dynamics in V. dahliae contributes to genomic variation, correlates with expression of pathogenicity-related genes, and potentially impacts the evolution of adaptive genomic regions.

The Structural, Functional and Evolutionary Impact of Transposable Elements in Eukaryotes

Transposable elements (TEs) are nearly ubiquitous in eukaryotes. The increase in genomic data, as well as progress in genome annotation and molecular biology techniques, have revealed the vast number of ways mobile elements have impacted the evolution of eukaryotes. In addition to being the main cause of difference in haploid genome size, TEs have affected the overall organization of genomes by accumulating preferentially in some genomic regions, by causing structural rearrangements or by modifying the recombination rate. Although the vast majority of insertions is neutral or deleterious, TEs have been an important source of evolutionary novelties and have played a determinant role in the evolution of fundamental biological processes. TEs have been recruited in the regulation of host genes and are implicated in the evolution of regulatory networks. They have also served as a source of protein-coding sequences or even entire genes. The impact of TEs on eukaryotic evolution is only now being fully appreciated and the role they may play in a number of biological processes, such as speciation and adaptation, remains to be deciphered.

The IDIP framework for assessing protein function and its application to the prion protein

The quest to determine the function of a protein can represent a profound challenge. Although this task is the mandate of countless research groups, a general framework for how it can be approached is conspicuously lacking. Moreover, even expectations for when the function of a protein can be considered to be 'known' are not well defined. In this review, we begin by introducing concepts pertinent to the challenge of protein function assignments. We then propose a framework for inferring a protein's function from four data categories: 'inheritance', 'distribution', 'interactions' and 'phenotypes' (IDIP). We document that the functions of proteins emerge at the intersection of inferences drawn from these data categories and emphasise the benefit of considering them in an evolutionary context. We then apply this approach to the cellular prion protein (PrP^C ), well known for its central role in prion diseases, whose function continues to be considered elusive by many investigators. We document that available data converge on the conclusion that the function of the prion protein is to control a critical post-translational modification of the neural cell adhesion molecule in the context of epithelial-to-mesenchymal transition and related plasticity programmes. Finally, we argue that this proposed function of PrP^C has already passed the test of time and is concordant with the IDIP framework in a way that other functions considered for this protein fail to achieve. We anticipate that the IDIP framework and the concepts analysed herein will aid the investigation of other proteins whose primary functional assignments have thus far been intractable.

Distribution of Merlin in eukaryotes and first report of DNA transposons in kinetoplastid protists

DNA transposons are defined as repeated DNA sequences that can move within the host genome through the action of transposases. The transposon superfamily Merlin was originally found mainly in animal genomes. Here, we describe a global distribution of the Merlin in animals, fungi, plants and protists, reporting for the first time their presence in Rhodophyceae, Metamonada, Discoba and Alveolata. We identified a great variety of potentially active Merlin families, some containing highly imperfect terminal inverted repeats and internal tandem repeats. Merlin-related sequences with no evidence of mobilization capacity were also observed and may be products of domestication. The evolutionary trees support that Merlin is likely an ancient superfamily, with early events of diversification and secondary losses, although repeated re-invasions probably occurred in some groups, which would explain its diversity and discontinuous distribution. We cannot rule out the possibility that the Merlin superfamily is the product of multiple horizontal transfers of related prokaryotic insertion sequences. Moreover, this is the first account of a DNA transposon in kinetoplastid flagellates, with conserved Merlin transposase identified in Bodo saltans and Perkinsela sp., whereas it is absent in trypanosomatids. Based on the level of conservation of the transposase and overlaps of putative open reading frames with Merlin, we propose that in protists it may serve as a raw material for gene emergence.

Transposons Increase Transcriptional Complexity: The Good Parasite?

A recent study by Cosby et al. sheds light on the role of transposons in the adaptive evolution of their hosts. These genetic elements were thought to be largely deleterious. However, when coupled with alternative splicing, there appears to be an exponential increase in the diversity of proteins encoded, which display novel functions and are conserved by natural selection.

Insertion of Badnaviral DNA in the Late Blight Resistance Gene (R1a) of Brinjal Eggplant (Solanum melongena)

Endogenous viral elements (EVEs) are widespread in plant genomes. They result from the random integration of viral sequences into host plant genomes by horizontal DNA transfer and have the potential to alter host gene expression. We performed a large-scale search for co-transcripts including caulimovirid and plant sequences in 1,678 plant and 230 algal species and characterized 50 co-transcripts in 45 distinct plant species belonging to lycophytes, ferns, gymnosperms and angiosperms. We found that insertion of badnavirus EVEs along with Ty-1 copia mobile elements occurred into a late blight resistance gene (R1) of brinjal eggplant (Solanum melongena) and wild relatives in genus Solanum and disrupted R1 orthologs. EVEs of two previously unreported badnaviruses were identified in the genome of S. melongena, whereas EVEs from an additional novel badnavirus were identified in the genome of S. aethiopicum, the cultivated scarlet eggplant. Insertion of these viruses in the ancestral lineages of the direct wild relatives of the eggplant would have occurred during the last 3 Myr, further supporting the distinctiveness of the group of the eggplant within the giant genus Solanum.