Dynamics Of Transposable Elements Research Articles

Comparative Analysis of Transposable Elements Dynamics in Fish with Different Sex Chromosome Systems.

Transposable elements (TEs) are widespread genomic components with substantial roles in genome evolution and sex chromosome differentiation. In this study, we compared the TE composition of three closely related fish with different sex chromosome systems: Megaleporinus elongatus (Z1Z1Z2Z2/Z1W1Z2W2), Megaleporinus macrocephalus (ZZ/ZW) (both with highly differentiated W sex chromosomes), and Leporinus friderici (without heteromorphic sex chromosomes). We created custom TE libraries for each species using clustering methods and manual annotation and prediction, and we predicted TE temporal dynamics through divergence-based analysis. The TE abundance ranged from 16 to 21% in the three mobilomes, with L. friderici having the lowest overall. Despite the recent amplification of TEs in all three species, we observed differing expansion activities, particularly between the two genera. Both Megaleporinus recently experienced high retrotransposon activity, with a reduction in DNA TEs, which could have implications in sex chromosome composition. In contrast, L. friderici showed the opposite pattern. Therefore, despite having similar TE compositions, Megaleporinus and Leporinus exhibit distinct TE histories that likely evolved after their separation, highlighting a rapid TE expansion over short evolutionary periods.

CAM evolution is associated with gene family expansion in an explosive bromeliad radiation.

The subgenus Tillandsia (Bromeliaceae) belongs to one of the fastest radiating clades in the plant kingdom and is characterised by the repeated evolution of Crassulacean acid metabolism (CAM). Despite its complex genetic basis, this water-conserving trait has evolved independently across many plant families and is regarded as a key innovation trait and driver of ecological diversification in Bromeliaceae. By producing high-quality genome assemblies of a Tillandsia species pair displaying divergent photosynthetic phenotypes, and combining genome-wide investigations of synteny, transposable element (TE) dynamics, sequence evolution, gene family evolution and temporal differential expression, we were able to pinpoint the genomic drivers of CAM evolution in Tillandsia. Several large-scale rearrangements associated with karyotype changes between the two genomes and a highly dynamic TE landscape shaped the genomes of Tillandsia. However, our analyses show that rewiring of photosynthetic metabolism is mainly obtained through regulatory evolution rather than coding sequence evolution, as CAM-related genes are differentially expressed across a 24-hour cycle between the two species but are not candidates of positive selection. Gene orthology analyses reveal that CAM-related gene families manifesting differential expression underwent accelerated gene family expansion in the constitutive CAM species, further supporting the view of gene family evolution as a driver of CAM evolution.

Blessing or curse: how the epigenetic resolution of host-transposable element conflicts shapes their evolutionary dynamics.

Transposable elements (TEs) are selfish genetic elements whose antagonistic interactions with hosts represent a common genetic conflict in eukaryotes. To resolve this conflict, hosts have widely adopted epigenetic silencing that deposits repressive marks at TEs. However, this mechanism is imperfect and fails to fully halt TE replication. Furthermore, TE epigenetic silencing can inadvertently spread repressive marks to adjacent functional sequences, a phenomenon considered a 'curse' of this conflict resolution. Here, we used forward simulations to explore how TE epigenetic silencing and its harmful side effects shape the evolutionary dynamics of TEs and their hosts. Our findings reveal that epigenetic silencing allows TEs and their hosts to stably coexist under a wide range of conditions, because the underlying molecular mechanisms give rise to copy-number dependency of the strength of TE silencing. Interestingly, contrary to intuitive expectations that TE epigenetic silencing should evolve to be as strong as possible, we found a selective benefit for modifier alleles that weaken TE silencing under biologically feasible conditions. These results reveal that the dual nature of TE epigenetic silencing, with both positive and negative effects, complicates its evolutionary trajectory and makes it challenging to determine whether TE epigenetic silencing is a 'blessing' or a 'curse'.

Quantitative pathogenicity and host adaptation in a fungal plant pathogen revealed by whole-genome sequencing

Knowledge of genetic determinism and evolutionary dynamics mediating host-pathogen interactions is essential to manage fungal plant diseases. Studies on the genetic architecture of fungal pathogenicity often focus on large-effect effector genes triggering strong, qualitative resistance. It is not clear how this translates to predominately quantitative interactions. Here, we use the Zymoseptoria tritici-wheat model to elucidate the genetic architecture of quantitative pathogenicity and mechanisms mediating host adaptation. With a multi-host genome-wide association study, we identify 19 high-confidence candidate genes associated with quantitative pathogenicity. Analysis of genetic diversity reveals that sequence polymorphism is the main evolutionary process mediating differences in quantitative pathogenicity, a process that is likely facilitated by genetic recombination and transposable element dynamics. Finally, we use functional approaches to confirm the role of an effector-like gene and a methyltransferase in phenotypic variation. This study highlights the complex genetic architecture of quantitative pathogenicity, extensive diversifying selection and plausible mechanisms facilitating pathogen adaptation.

Modeling transposable elements dynamics during polyploidization in plants

In this work we study the proliferation of transposable elements (TEs) and the epigenetic response of plants during the process of polyploidization. Through a deterministic model, expanding on our previous work on TE proliferation under epigenetic regulation, we study the long-term TE distribution and TE stability in the subgenomes of both autopolyploids and allopolyploids. We also explore different small-interfering RNA (siRNA) action modes on the subgenomes, including a model where siRNAs are not directed to specific genomes and one where siRNAs are directed – i.e. more active – in subgenomes with higher TE loads. In the autopolyploid case, we find long-term stable equilbria that tend to equilibrate the number of active TEs between subgenomes. In the allopolyploid case, directed siRNA action is fundamental to avoid a ”winner takes all” outcome of the competition between the TEs in the different subgenomes. We also show that decaying oscillations in the number of TEs occur naturally in all cases, perhaps explaining some of the observed features of ’genomic shock’ after hybridization events, and that the balance in the dynamics of the different types of siRNA is determinant for the synchronization of these oscillations.

Intraspecific Variation of Transposable Elements Reveals Differences in the Evolutionary History of Fungal Phytopathogen Pathotypes.

Transposable elements (TEs) contribute to intraspecific variation and play important roles in the evolution of fungal genomes. However, our understanding of the processes that shape TE landscapes is limited, as is our understanding of the relationship between TE content, population structure, and evolutionary history of fungal species. Fungal plant pathogens, which often have host-specific populations, are useful systems in which to study intraspecific TE content diversity. Here, we describe TE dynamics in five lineages of Magnaporthe oryzae, the fungus that causes blast disease of rice, wheat, and many other grasses. We identified differences in TE content across these lineages and showed that recent lineage-specific expansions of certain TEs have contributed to overall greater TE content in rice-infecting and Setaria-infecting lineages. We reconstructed the evolutionary histories of long terminal repeat-retrotransposon expansions and found that in some cases they were caused by complex proliferation dynamics of one element and in others by multiple elements from an older population of TEs multiplying in parallel. Additionally, we found evidence suggesting the recent transfer of a DNA transposon between rice- and wheat-infecting M. oryzae lineages and a region showing evidence of homologous recombination between those lineages, which could have facilitated such a transfer. By investigating intraspecific TE content variation, we uncovered key differences in the proliferation dynamics of TEs in various pathotypes of a fungal plant pathogen, giving us a better understanding of the evolutionary history of the pathogen itself.

The impact of paramutations on the invasion dynamics of transposable elements.

According to the prevailing view, the trap model, the activity of invading transposable elements (TEs) is greatly reduced when a TE copy jumps into a piRNA cluster, which triggers the emergence of piRNAs that silence the TE. One crucial component in the host defence are paramutations. Mediated by maternally deposited piRNAs, paramutations convert TE insertions into piRNA producing loci, thereby transforming selfish TEs into agents of the host defence. Despite this significant effect, the impact of paramutations on the dynamics of TE invasions remains unknown. To address this issue, we performed extensive forward simulations of TE invasions with piRNA clusters and paramutations. We found that paramutations significantly affect TE dynamics, by accelerating the silencing of TE invasions, reducing the number of insertions accumulating during the invasions and mitigating the fitness cost of TEs. We also demonstrate that piRNA production induced by paramutations, an epigenetically inherited trait, may be positively selected. Finally, we show that paramutations may account for three important open problems with the trap model. Firstly, paramutated TE insertions may compensate for the insufficient number of insertions in piRNA clusters observed in previous studies. Secondly, paramutations may explain the discrepancy between the observed and the expected abundance of different TE families in Drosophila melanogaster. Thirdly, piRNA clusters may be crucial to trigger the host defence, but paramutations render the clusters dispensable once the defence has been established. This could account for the lack of TE activation when three major piRNA clusters were deleted in a previous study.

Selection signatures and population dynamics of transposable elements in lima bean

The domestication process in lima bean (Phaseolus lunatus L.) involves two independent events, within the Mesoamerican and Andean gene pools. This makes lima bean an excellent model to understand convergent evolution. The mechanisms of adaptation followed by Mesoamerican and Andean landraces are largely unknown. Genes related to these adaptations can be selected by identification of selective sweeps within gene pools. Previous genetic analyses in lima bean have relied on Single Nucleotide Polymorphism (SNP) loci, and have ignored transposable elements (TEs). Here we show the analysis of whole-genome sequencing data from 61 lima bean accessions to characterize a genomic variation database including TEs and SNPs, to associate selective sweeps with variable TEs and to predict candidate domestication genes. A small percentage of genes under selection are shared among gene pools, suggesting that domestication followed different genetic avenues in both gene pools. About 75% of TEs are located close to genes, which shows their potential to affect gene functions. The genetic structure inferred from variable TEs is consistent with that obtained from SNP markers, suggesting that TE dynamics can be related to the demographic history of wild and domesticated lima bean and its adaptive processes, in particular selection processes during domestication.

Telomere-to-telomere assemblies of 142 strains characterize the genome structural landscape in Saccharomyces cerevisiae.

Pangenomes provide access to an accurate representation of the genetic diversity of species, both in terms of sequence polymorphisms and structural variants (SVs). Here we generated the Saccharomyces cerevisiae Reference Assembly Panel (ScRAP) comprising reference-quality genomes for 142 strains representing the species' phylogenetic and ecological diversity. The ScRAP includes phased haplotype assemblies for several heterozygous diploid and polyploid isolates. We identified circa (ca.) 4,800 nonredundant SVs that provide a broad view of the genomic diversity, including the dynamics of telomere length and transposable elements. We uncovered frequent cases of complex aneuploidies where large chromosomes underwent large deletions and translocations. We found that SVs can impact gene expression near the breakpoints and substantially contribute to gene repertoire evolution. We also discovered that horizontally acquired regions insert at chromosome ends and can generate new telomeres. Overall, the ScRAP demonstrates the benefit of a pangenome in understanding genome evolution at population scale.

The composition dynamics of transposable elements in human blastocysts.

Transposable elements (TEs) are mobile DNA sequences that can replicate themselves and play significant roles in embryo development and chromosomal structure remodeling. In this study, we investigated the variation of TEs in blastocysts with different parental genetic backgrounds. We analyzed the proportions of 1137 TEs subfamilies from six classes at the DNA level using Bowtie2 and PopoolationTE2 in 196 blastocysts with abnormal parental chromosomal diseases. Our findings revealed that the parental karyotype was the dominant factor influencing TEs frequencies. Out of the 1116 subfamilies, different frequencies were observed in blastocysts with varying parental karyotypes. The development stage of blastocysts was the second most crucial factor influencing TEs proportions. A total of 614 subfamilies exhibited different proportions at distinct blastocyst stages. Notably, subfamily members belonging to the Alu family showed a high proportion at stage 6, while those from the LINE class exhibited a high proportion at stage 3 and a low proportion at stage 6. Moreover, the proportions of some TEs subfamilies also varied depending on blastocyst karyotype, inner cell mass status, and outer trophectoderm status. We found that 48 subfamilies displayed different proportions between balanced and unbalanced blastocysts. Additionally, 19 subfamilies demonstrated varying proportions among different inner cell mass scores, and 43 subfamilies exhibited different proportions among outer trophectoderm scores. This study suggests that the composition of TEs subfamilies may be influenced by various factors and undergoes dynamic modulation during embryo development.

Transposable Elements Shape the Genome Diversity and the Evolution of Noctuidae Species.

Noctuidae is known to have high species diversity, although the genomic diversity of Noctuidae species has yet to be studied extensively. Investigation of transposable elements (TEs) in this family can improve our understanding of the genomic diversity of Noctuidae. In this study, we annotated and characterized genome-wide TEs in ten noctuid species belonging to seven genera. With multiple annotation pipelines, we constructed a consensus sequence library containing 1038-2826 TE consensus. The genome content of TEs showed high variation in the ten Noctuidae genomes, ranging from 11.3% to 45.0%. The relatedness analysis indicated that the TE content, especially the content of LINEs and DNA transposons, is positively correlated with the genome size (r = 0.86, p-value = 0.001). We identified SINE/B2 as a lineage-specific subfamily in Trichoplusia ni, a species-specific expansion of the LTR/Gypsy subfamily in Spodoptera exigua, and a recent expansion of SINE/5S subfamily in Busseola fusca. We further revealed that of the four TE classes, only LINEs showed phylogenetic signals with high confidence. We also examined how the expansion of TEs contributed to the evolution of noctuid genomes. Moreover, we identified 56 horizontal transfer TE (HTT) events among the ten noctuid species and at least three HTT events between the nine Noctuidae species and 11 non-noctuid arthropods. One of the HTT events of a Gypsy transposon might have caused the recent expansion of the Gypsy subfamily in the S. exigua genome. By determining the TE content, dynamics, and HTT events in the Noctuidae genomes, our study emphasized that TE activities and HTT events substantially impacted the Noctuidae genome evolution.

All families of transposable elements were active in the recent wheat genome evolution and polyploidy had no impact on their activity.

Bread wheat (Triticum aestivum L.) is a major crop and its genome is one of the largest ever assembled at reference-quality level. It is 15Gb, hexaploid, with 85% of transposable elements (TEs). Wheat genetic diversity was mainly focused on genes and little is known about the extent of genomic variability affecting TEs, transposition rate, and the impact of polyploidy. Multiple chromosome-scale assemblies are now available for bread wheat and for its tetraploid and diploid wild relatives. In this study, we computed base pair-resolved, gene-anchored, whole genome alignments of A, B, and D lineages at different ploidy levels in order to estimate the variability that affects the TE space. We used assembled genomes of 13 T.aestivum cultivars (6x=AABBDD) and a single genome for Triticum durum (4x=AABB), Triticum dicoccoides (4x=AABB), Triticum urartu (2x=AA), and Aegilops tauschii (2x=DD). We show that 5%-34% of the TE fraction is variable, depending on the species divergence. Between 400 and 13,000 novel TE insertions per subgenome were detected. We found lineage-specific insertions for nearly all TE families in di-, tetra-, and hexaploids. No burst of transposition was observed and polyploidization did not trigger any boost of transposition. This study challenges the prevailing idea of wheat TE dynamics and is more in agreement with an equilibrium model of evolution.

Phenotype stability and dynamics of transposable elements in a strain of the microalga Tisochrysis lutea with improved lipid traits.

Microalgal domestication is an expanding research field that aims to multiply and accelerate the potential of microalgae for various biotechnological purposes. We investigated the stability of improved lipid traits and genetic changes of a domesticated strain of the haptophyte Tisochrysis lutea, TisoS2M2, previously obtained by a mutation-selection improvement program. After 7 years of maintenance, TisoS2M2 still displayed improved lipid traits compared with the native strain, demonstrating that a mutation-selection improvement program is suitable for obtaining a domesticated strain with stable, improved phenotype over time. We identified specific genetic variations between the native and domesticated strains and focused on the dynamics of transposable elements (TEs). DNA transposons mainly caused specific TE indels of the domesticated strain TisoS2M2, and some specific TE indels may have impacted genes associated to the neutral lipid pathway. We revealed transposition events for TEs in T. lutea and discussed on the potential role of the improvement program on their activity.

Causes and consequences of linkage disequilibrium among transposable elements within eukaryotic genomes.

Sex and recombination can affect the dynamics of transposable elements (TEs) in various ways: while sex is expected to help TEs to spread within populations, the deleterious effect of ectopic recombination among transposons represents a possible source of purifying selection limiting their number. Furthermore, recombination may also increase the efficiency of selection against TEs by reducing selective interference among loci. In order to better understand the effects of recombination and reproductive systems on TE dynamics, this article provides analytical expressions for the linkage disequilibrium among TEs in a classical model in which TE number is stabilized by synergistic purifying selection. The results show that positive linkage disequilibrium is predicted in infinite populations despite negative epistasis, due to the effect of the transposition process. Positive linkage disequilibrium may substantially inflate the variance in the number of elements per genome in the case of partially selfing or partially clonal populations. Finite population size tends to generate negative linkage disequilibrium (Hill-Robertson effect), the relative importance of this effect increasing with the degree of linkage among loci. The model is then extended in order to explore how TEs may affect selection for recombination. While positive linkage disequilibrium generated by transposition generally disfavors recombination, the Hill-Robertson effect may represent a non-negligible source of indirect selection for recombination when TEs are abundant. However, the direct fitness cost imposed by ectopic recombination among elements generally drives the population towards low-recombination regimes, at which TEs cannot be maintained at a stable equilibrium.

TrEMOLO: accurate transposable element allele frequency estimation using long-read sequencing data combining assembly and mapping-based approaches

Transposable Element MOnitoring with LOng-reads (TrEMOLO) is a new software that combines assembly- and mapping-based approaches to robustly detect genetic elements called transposable elements (TEs). Using high- or low-quality genome assemblies, TrEMOLO can detect most TE insertions and deletions and estimate their allele frequency in populations. Benchmarking with simulated data revealed that TrEMOLO outperforms other state-of-the-art computational tools. TE detection and frequency estimation by TrEMOLO were validated using simulated and experimental datasets. Therefore, TrEMOLO is a comprehensive and suitable tool to accurately study TE dynamics. TrEMOLO is available under GNU GPL3.0 at https://github.com/DrosophilaGenomeEvolution/TrEMOLO.

Genome Size Changes by Duplication, Divergence, and Insertion in Caenorhabditis Worms.

Genome size has been measurable since the 1940s but we still do not understand genome size variation. Caenorhabditis nematodes show strong conservation of chromosome number but vary in genome size between closely related species. Androdioecy, where populations are composed of males and self-fertile hermaphrodites, evolved from outcrossing, female-male dioecy, three times in this group. In Caenorhabditis, androdioecious genomes are 10-30% smaller than dioecious species, but in the nematode Pristionchus, androdioecy evolved six times and does not correlate with genome size. Previous hypotheses include genome size evolution through: 1) Deletions and "genome shrinkage" in androdioecious species; 2) Transposable element (TE) expansion and DNA loss through large deletions (the "accordion model"); and 3) Differing TE dynamics in androdioecious and dioecious species. We analyzed nematode genomes and found no evidence for these hypotheses. Instead, nematode genome sizes had strong phylogenetic inertia with increases in a few dioecious species, contradicting the "genome shrinkage" hypothesis. TEs did not explain genome size variation with the exception of the DNA transposon Mutator which was twice as abundant in dioecious genomes. Across short and long evolutionary distances Caenorhabditis genomes evolved through small structural mutations including gene-associated duplications and insertions. Seventy-one protein families had significant, parallel decreases across androdioecious Caenorhabditis including genes involved in the sensory system, regulatory proteins and membrane-associated immune responses. Our results suggest that within a dynamic landscape of frequent small rearrangements in Caenorhabditis, reproductive mode mediates genome evolution by altering the precise fates of individual genes, proteins, and the phenotypes they underlie.

Genome sizes and repeatome evolution in zoantharians (Cnidaria: Hexacorallia: Zoantharia).

Across eukaryotes, large variations of genome sizes have been observed even between closely related species. Transposable elements as part of the repeated DNA have been proposed and confirmed as one of the most important contributors to genome size variation. However, the evolutionary implications of genome size variation and transposable element dynamics are not well understood. Together with phenotypic traits, they are commonly referred to as the "C-value enigma". The order Zoantharia are benthic cnidarians found from intertidal zones to the deep sea, and some species are particularly abundant in coral reefs. Despite their high ecological relevance, zoantharians have yet to be largely studied from the genomic point of view. This study aims at investigating the role of the repeatome (total content of repeated elements) in genome size variations across the order Zoantharia. To this end, whole-genomes of 32 zoantharian species representing five families were sequenced. Genome sizes were estimated and the abundances of different repeat classes were assessed. In addition, the repeat overlap between species was assessed by a sequence clustering method. The genome sizes in the dataset varied up to 2.4 fold magnitude. Significant correlations between genome size, repeated DNA content and transposable elements, respectively (Pearson's correlation test R2=0.47, p=0.0016; R2=0.22, p=0.05) were found, suggesting their involvement in the dynamics of genome expansion and reduction. In all species, long interspersed nuclear elements and DNA transposons were the most abundant identified elements. These transposable elements also appeared to have had a recent expansion event. This was in contrast to the comparative clustering analysis which revealed species-specific patterns of satellite elements' amplification. In summary, the genome sizes of zoantharians likely result from the complex dynamics of repeated elements. Finally, the majority of repeated elements (up to 70%) could not be annotated to a known repeat class, highlighting the need to further investigate non-model cnidarian genomes. More research is needed to understand how repeated DNA dynamics relate to zoantharian evolution and their biology.

A review of strategies used to identify transposition events in plant genomes.

Transposable elements (TEs) were initially considered redundant and dubbed 'junk DNA'. However, more recently they were recognized as an essential element of genome plasticity. In nature, they frequently become active upon exposition of the host to stress conditions. Even though most transposition events are neutral or even deleterious, occasionally they may happen to be beneficial, resulting in genetic novelty providing better fitness to the host. Hence, TE mobilization may promote adaptability and, in the long run, act as a significant evolutionary force. There are many examples of TE insertions resulting in increased tolerance to stresses or in novel features of crops which are appealing to the consumer. Possibly, TE-driven de novo variability could be utilized for crop improvement. However, in order to systematically study the mechanisms of TE/host interactions, it is necessary to have suitable tools to globally monitor any ongoing TE mobilization. With the development of novel potent technologies, new high-throughput strategies for studying TE dynamics are emerging. Here, we present currently available methods applied to monitor the activity of TEs in plants. We divide them on the basis of their operational principles, the position of target molecules in the process of transposition and their ability to capture real cases of actively transposing elements. Their possible theoretical and practical drawbacks are also discussed. Finally, conceivable strategies and combinations of methods resulting in an improved performance are proposed.

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

BackgroundBats 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.ResultsWe 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.ConclusionCollectively, 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.

Population genetics of transposable element load: A mechanistic account of observed overdispersion

In an empirical analysis of transposable element (TE) abundance within natural populations of Mimulus guttatus and Drosophila melanogaster, we found a surprisingly high variance of TE count (e.g., variance-to-mean ratio on the order of 10 to 300). To obtain insight regarding the evolutionary genetic mechanisms that underlie the overdispersed population distributions of TE abundance, we developed a mathematical model of TE population genetics that includes the dynamics of element proliferation and purifying selection on TE load. The modeling approach begins with a master equation for a birth-death process and extends the predictions of the classical theory of TE dynamics in several ways. In particular, moment-based analyses of population distributions of TE load reveal that overdispersion is likely to arise via copy-and-paste proliferation dynamics, especially when the elementary processes of proliferation and excision are approximately balanced. Parameter studies and analytic work confirm this result and further suggest that overdispersed population distributions of TE abundance are probably not a consequence of purifying selection on total element load.