Worldwide distribution of transposable element copy number in natural populations of Drosophila simulans

Catalogue of the Diptera (Insecta) of Morocco— an annotated checklist, with distributions and a bibliography

The faunistic knowledge of the Diptera of Morocco recorded from 1787 to 2021 is summarized and updated in this first catalogue of Moroccan Diptera species. A total of 3057 species, classified into 948 genera and 93 families (21 Nematocera and 72 Brachycera), are listed. Taxa (superfamily, family, genus and species) have been updated according to current interpretations, based on reviews in the literature, the expertise of authors and contributors, and recently conducted fieldwork. Data to compile this catalogue were primarily gathered from the literature. In total, 1225 references were consulted and some information was also obtained from online databases. Each family was reviewed and the checklist updated by the respective taxon expert(s), including the number of species that can be expected for that family in Morocco. For each valid species, synonyms known to have been used for published records from Morocco are listed under the currently accepted name. Where available, distribution within Morocco is also included. One new combination is proposed: Assuaniamelanoleuca (Séguy, 1941), comb. nov. (Chloropidae).

A Transposon Story: From TE Content to TE Dynamic Invasion of Drosophila Genomes Using the Single-Molecule Sequencing Technology from Oxford Nanopore

Transposable elements (TEs) are the main components of genomes. However, due to their repetitive nature, they are very difficult to study using data obtained with short-read sequencing technologies. Here, we describe an efficient pipeline to accurately recover TE insertion (TEI) sites and sequences from long reads obtained by Oxford Nanopore Technology (ONT) sequencing. With this pipeline, we could precisely describe the landscapes of the most recent TEIs in wild-type strains of Drosophila melanogaster and Drosophila simulans. Their comparison suggests that this subset of TE sequences is more similar than previously thought in these two species. The chromosome assemblies obtained using this pipeline also allowed recovering piRNA cluster sequences, which was impossible using short-read sequencing. Finally, we used our pipeline to analyze ONT sequencing data from a D. melanogaster unstable line in which LTR transposition was derepressed for 73 successive generations. We could rely on single reads to identify new insertions with intact target site duplications. Moreover, the detailed analysis of TEIs in the wild-type strains and the unstable line did not support the trap model claiming that piRNA clusters are hotspots of TE insertions.

Transposable elements in individual genotypes of Drosophila simulans

Transposable elements are abundant, dynamic components of the genome that affect organismal phenotypes and fitness. In Drosophila melanogaster, they have increased in abundance as the species spread out of Africa, and different populations differ in their transposable element content. However, very little is currently known about how transposable elements differ between individual genotypes, and how that relates to the population dynamics of transposable elements overall. The sister species of D. melanogaster, D. simulans, has also recently become cosmopolitan, and panels of inbred genotypes exist from cosmopolitan and African flies. Therefore, we can determine whether the differences in colonizing populations are repeated in D. simulans, what the dynamics of transposable elements are in individual genotypes, and how that compares to wild flies. After estimating copy number in cosmopolitan and African D. simulans, I find that transposable element load is higher in flies from cosmopolitan populations. In addition, transposable element load varies considerably between populations, between genotypes, but not overall between wild and inbred lines. Certain genotypes either contain active transposable elements or are more permissive of transposition and accumulate copies of particular transposable elements. Overall, it is important to quantify genotype-specific transposable element dynamics as well as population averages to understand the dynamics of transposable element accumulation over time.

SINE Retrotransposon variation drives Ecotypic disparity in natural populations of Coilia nasus

Background

SINEs are a type of nonautonomous retrotransposon that can transpose from one site to be integrated elsewhere in an organism genome. SINE insertion can give rise to genetic variants and regulate gene expression, allowing organisms to acquire new adaptive capacity. Studies on this subject have focused on the impacts of SINEs on genes. However, ecological disparities in fish have not yet been explained by SINEs.

Results

New SINEs were isolated from Coilia nasus, which has two ecotypes—migratory and resident—that differ in their spawning and migration behaviors. The SINEs possess two structures that resemble a tRNA gene and a LINE retrotransposon tail. Comparison of olfactory tissue transcriptomes, intact SINE transcript copies were detected in only the migratory fish at the initial retrotransposition stage. The SINE DNA copy numbers were higher in the resident type than in the migratory type, while the frequency of SINE insertion was higher in the migratory type than in the resident type. Furthermore, SINE insertions can lead to new repeats of short DNA fragments in the genome, along with target site duplications. SINEs in the resident type have undergone excision via a mechanism in which predicted cleavage sites are formed by mutations, resulting in gaps that are then filled by microsatellites via microhomology-induced replication.

Conclusions

Notably, SINEs in the resident type have undergone strong natural selection, causing genomic heteroplasmy and driving ecological diversity of C. nasus. Our results reveal possible evolutionary mechanisms underlying the ecological diversity at the interface between SINE mobilization and organism defense.

TEtools facilitates big data expression analysis of transposable elements and reveals an antagonism between their activity and that of piRNA genes

Over recent decades, substantial efforts have been made to understand the interactions between host genomes and transposable elements (TEs). The impact of TEs on the regulation of host genes is well known, with TEs acting as platforms of regulatory sequences. Nevertheless, due to their repetitive nature it is considerably hard to integrate TE analysis into genome-wide studies. Here, we developed a specific tool for the analysis of TE expression: TEtools. This tool takes into account the TE sequence diversity of the genome, it can be applied to unannotated or unassembled genomes and is freely available under the GPL3 (https://github.com/l-modolo/TEtools). TEtools performs the mapping of RNA-seq data obtained from classical mRNAs or small RNAs onto a list of TE sequences and performs differential expression analyses with statistical relevance. Using this tool, we analyzed TE expression from five Drosophila wild-type strains. Our data show for the first time that the activity of TEs is strictly linked to the activity of the genes implicated in the piwi-interacting RNA biogenesis and therefore fits an arms race scenario between TE sequences and host control genes.

Tempo and Mode of Transposable Element Activity in Drosophila

The evolutionary dynamics of transposable element (TE) insertions have been of continued interest since TE activity has important implications for genome evolution and adaptation. Here, we infer the transposition dynamics of TEs by comparing their abundance in natural D. melanogaster and D. simulans populations. Sequencing pools of more than 550 South African flies to at least 320-fold coverage, we determined the genome wide TE insertion frequencies in both species. We suggest that the predominance of low frequency insertions in the two species (>80% of the insertions have a frequency <0.2) is probably due to a high activity of more than 58 families in both species. We provide evidence for 50% of the TE families having temporally heterogenous transposition rates with different TE families being affected in the two species. While in D. melanogaster retrotransposons were more active, DNA transposons showed higher activity levels in D. simulans. Moreover, we suggest that LTR insertions are mostly of recent origin in both species, while DNA and non-LTR insertions are older and more frequently vertically transmitted since the split of D. melanogaster and D. simulans. We propose that the high TE activity is of recent origin in both species and a consequence of the demographic history, with habitat expansion triggering a period of rapid evolution.

Variable expression levels detected in the Drosophila effectors of piRNA biogenesis

piRNAs (piwi-interacting RNAs) are a class of small interfering RNAs that play a major role in the regulation of transposable elements (TEs) in Drosophila and are considered of fundamental importance in gonadal development. Genes encoding the effectors of the piRNA machinery are thus often thought to be highly constrained. On the contrary, as actors of genetic immunity, these genes have also been shown to evolve rapidly and display a high level of sequence variability. In order to assess the support for these competing models, we analyzed seven genes of the piRNA pathway using a collection of wild-type strains of Drosophila simulans, which are known to display significant variability in their TE content between strains. We showed that these genes exhibited wide variation in transcript levels, and we discuss some evolutionary considerations regarding the observed variability in TE copy numbers.

tirant, a newly discovered active endogenous retrovirus in Drosophila simulans

Endogenous retroviruses have the ability to become permanently integrated into the genomes of their host, and they are generally transmitted vertically from parent to progeny. With the exception of gypsy, few endogenous retroviruses have been identified in insects. In this study, we describe the tirant endogenous retrovirus in a subset of Drosophila simulans natural populations. By focusing on the envelope gene, we show that the entire retroviral cycle (transcription, translation, and retrotransposition) can be completed for tirant within one population of this species.

Comparative analysis of transposable elements in the melanogaster subgroup sequenced genomes

Transposable elements (TEs) are indwelling components of genomes, and their dynamics have been a driving force in genome evolution. Although we now have more information concerning their amounts and characteristics in various organisms, we still have little data from overall comparisons of their sequences in very closely-related species. While the Drosophila melanogaster genome has been extensively studied, we have only limited knowledge regarding the precise TE sequences in the genomes of the related species Drosophila simulans, Drosophila sechellia and Drosophila yakuba. In this study we analyzed the number and structure of TE copies in the sequenced genomes of these four species. Our findings show that, unexpectedly, the number of TE insertions in D. simulans is greater than that in D. melanogaster, but that most of the copies in D. simulans are degraded and in small fragments, as in D. sechellia and D. yakuba. This suggests that all three species were invaded by numerous TEs a long time ago, but have since regulated their activity, as the present TE copies are degraded, with very few full-length elements. In contrast, in D. melanogaster, a recent activation of TEs has resulted in a large number of almost-identical TE copies. We have detected variants of some TEs in D. simulans and D. sechellia, that are almost identical to the reference TE sequences in D. melanogaster, suggesting that D. melanogaster has recently been invaded by active TE variants from the other species. Our results indicate that the three species D. simulans, D. sechellia, and D. yakuba seem to be at a different stage of their TE life cycle when compared to D. melanogaster. Moreover, we show that D. melanogaster has been invaded by active TE variants for several TE families likely to come from D. simulans or the ancestor of D. simulans and D. sechellia. The numerous horizontal transfer events implied to explain these results could indicate introgression events between these species.

Genomic distribution of retrotransposons 297, 1731, copia, mdg1 and roo in the Drosophila melanogaster species subgroup

The intragenomic distribution of five retrotransposon families (297, 1731, copia, mdg1 and roo) in the species of the melanogaster complex was studied by comparing results of the Southern blotting technique in males and females with those of in situ hybridization. The degree of structural polymorphism of each family in the different species was also investigated by restriction enzyme analysis. It was found that genomic distribution is a trait that depends on the family and species. The distribution of roo is mainly euchromatic in the four species and 1731 is heterochromatic, but the distribution of families 297, copia and mdg1 is markedly different in the melanogaster and simulans clades. These families were mainly euchromatic in D. melanogaster but heterochromatic in its sibling species. In the simulans clade most copia and mdg1 elements are located on chromosome Y. Differences in genomic distribution are unrelated with structural conservation. The relation of intragenomic distribution to phylogeny, transpositional activity and the role of the host genome are discussed.

The protist Trichomonas vaginalis harbors multiple lineages of transcriptionally active Mutator-like elements

Background

For three decades the Mutator system was thought to be exclusive of plants, until the first homolog representatives were characterized in fungi and in early-diverging amoebas earlier in this decade.

Results

Here, we describe and characterize four families of Mutator-like elements in a new eukaryotic group, the Parabasalids. These Trichomonas vaginalis Mutator- like elements, or TvMULEs, are active in T. vaginalis and patchily distributed among 12 trichomonad species and isolates. Despite their relatively distinctive amino acid composition, the inclusion of the repeats TvMULE1, TvMULE2, TvMULE3 and TvMULE4 into the Mutator superfamily is justified by sequence, structural and phylogenetic analyses. In addition, we identified three new TvMULE-related sequences in the genome sequence of Candida albicans. While TvMULE1 is a member of the MuDR clade, predominantly from plants, the other three TvMULEs, together with the C. albicans elements, represent a new and quite distinct Mutator lineage, which we named TvCaMULEs. The finding of TvMULE1 sequence inserted into other putative repeat suggests the occurrence a novel TE family not yet described.

Conclusion

These findings expand the taxonomic distribution and the range of functional motif of MULEs among eukaryotes. The characterization of the dynamics of TvMULEs and other transposons in this organism is of particular interest because it is atypical for an asexual species to have such an extreme level of TE activity; this genetic landscape makes an interesting case study for causes and consequences of such activity. Finally, the extreme repetitiveness of the T. vaginalis genome and the remarkable degree of sequence identity within its repeat families highlights this species as an ideal system to characterize new transposable elements.

Genomic environment influences the dynamics of the tirant LTR retrotransposon in Drosophila

Combining genome sequence analysis and functional analysis, we show that some full-length copies of tirant are present in heterochromatic regions in Drosophila simulans and that when tested in vitro, these copies have a functional promoter. However, when inserted in heterochromatic regions, tirant copies are inactive in vivo, and only transcription of euchromatic copies can be detected. Thus, our data indicate that the localization of the element is a hallmark of its activity in vivo and raise the question of genomic invasions by transposable elements and the importance of their genomic integration sites.

An old bilbo-like non-LTR retroelement insertion provides insight into the relationship of species of the virilis group

In Drosophila, at the population and species level, fixation of a TE insertion is an unlikely fate. Of the few reported fixations at the species level most involve non-LTR retroelements. In this work we report the fixation of a non-LTR retroelement in five species (Drosophila littoralis, Drosophila virilis, Drosophila lummei, Drosophila americana and Drosophila novamexicana) of the virilis group of Drosophila. In most species, this TE insertion is being lost through the accumulation of small deletions, but there is also evidence for the accumulation of large deletions. In the americana lineage an insertion of about 900 bp of the non-LTR retroelement is a marker for the Xc inversion. This insertion is, at most, 80 kb away from the basal Xc inversion breakpoint. The presence of a bilbo-like element in D. littoralis but not in D. kanekoi, suggests that D. littoralis is more closely related to species of the virilis phylad than to species of the montana phylad, which is in contrast with the traditional view. Nevertheless, the phylogenetic analyses here performed using a 7 gene dataset suggest that D. littoralis is indeed more closely related to species of the virilis phylad than to species of the montana phylad. The re-evaluation of the phylogenetic relationship of the species of the virilis group, under the assumption of a relaxed molecular clock, results in an estimated age of the bilbo-like element insertion of at least 7.5 Mya.

Losing helena: the extinction of a drosophila line-like element

Background

Transposable elements (TEs) are major players in evolution. We know that they play an essential role in genome size determination, but we still have an incomplete understanding of the processes involved in their amplification and elimination from genomes and populations. Taking advantage of differences in the amount and distribution of the Long Interspersed Nuclear Element (LINE), helena in Drosophila melanogaster and D. simulans, we analyzed the DNA sequences of copies of this element in samples of various natural populations of these two species.

Results

In situ hybridization experiments revealed that helena is absent from the chromosome arms of D. melanogaster, while it is present in the chromosome arms of D. simulans, which is an unusual feature for a TE in these species. Molecular analyses showed that the helena sequences detected in D. melanogaster were all deleted copies, which diverged from the canonical element. Natural populations of D. simulans have several copies, a few of them full-length, but most of them internally deleted.

Conclusion

Overall, our data suggest that a mechanism that induces internal deletions in the helena sequences is active in the D. simulans genome.

Amplification of the 1731 LTR retrotransposon in Drosophila melanogaster cultured cells: Origin of neocopies and impact on the genome

Transposable elements (TEs), represent a large fraction of the eukaryotic genome. In Drosophila melanogaster, about 20% of the genome corresponds to such middle repetitive DNA dispersed sequences. A fraction of TEs is composed of elements showing a retrovirus-like structure, the LTR-retrotransposons, the first TEs to be described in the Drosophila genome. Interestingly, in D. melanogaster embryonic immortal cell culture genomes the copy number of these LTR-retrotransposons was revealed to be higher than the copy number in the Drosophila genome, presumably as the result of transposition of some copies to new genomic locations [Potter, S.S., Brorein Jr., W.J., Dunsmuir, P., Rubin, G.M., 1979. Transposition of elements of the 412, copia and 297 dispersed repeated gene families in Drosophila. Cell 17, 415-427; Junakovic, N., Di Franco, C., Best-Belpomme, M., Echalier, G., 1988. On the transposition of copia-like nomadic elements in cultured Drosophila cells. Chromosoma 97, 212-218]. This suggests that so many transpositions modified the genome organisation and consequently the expression of targeted genes. To understand what has directed the transposition of TEs in Drosophila cell culture genomes, a search to identify the newly transposed copies was undertaken using 1731, a LTR-retrotransposon. A comparison between 1731 full-length elements found in the fly sequenced genome (y(1); cn(1)bw(1), sp(1) stock) and 1731 full-length elements amplified by PCR in the two cell line was done. The resulting data provide evidence that all 1731 neocopies were derived from a single copy slightly active in the Drosophila genome and subsequently strongly activated in cultured cells; and that this active copy is related to a newly evolved genomic variant (Kalmykova, A.I., et al., 2004. Selective expansion of the newly evolved genomic variants of retrotransposon 1731 in the Drosophila genomes. Mol. Biol. Evol. 21, 2281-2289). Moreover, neocopies are shown to be inserted in different sets of genes in the two cell lines suggesting they might be involved in the biological and physiological differences observed between Kc and S2 cell lines.

Evolutionary pathways of the tirant LTR retrotransposon in the Drosophila melanogaster subgroup of species

Tirant, a LTR retrotransposon with copies scattered over the chromosome arms of Drosophila melanogaster, is in the process of being lost from the chromosome arms of most natural populations of the sister species D. simulans. In an attempt to clarify the dynamics and evolution of tirant, we have studied the regulatory and reverse transcriptase regions in copies of the nine closely related species of the D. melanogaster subgroup. We show that tirant is mainly vertically transmitted in these species, with the exception of a horizontal transfer event from an ancestor of D. melanogaster to D. teissieri. We propose that, in four of the species (D. melanogaster, D. simulans, D. sechellia, and D. mauritiana), the observed patterns of evolution of the regulatory region vary with genome constraints and with the history and biogeography of the species.

The evolutionary history of Drosophila buzzatii. XXXVI. Molecular structural analysis of Osvaldo retrotransposon insertions in colonizing populations unveils drift effects in founder events

Previous work on transposable element distribution in colonizing populations of Drosophila buzzatii revealed a high frequency of occupancy in several chromosomal sites. Two explanatory hypotheses were advanced: the founder hypothesis, by which founder genetic drift was responsible, and the unstable hypothesis that assigns this unusual distribution to bursts of transposition toward some chromosomal sites. Here, we study the molecular structure of three euchromatic Osvaldo clones isolated from sites occupied at high (A4 and B9) and low frequency (B4) in colonizing populations, to test these hypotheses. Large insertions, duplications, and indels in the Osvaldo coding region and LTR were detected in the A4 clone and a truncated Osvaldo with many substitutions was found in the B9 clone. These altered sequences indicate that the two copies of this retroelement are precolonization insertions. Interestingly, the LTR of the A4 clone and the reverse transcriptase region of B9 show identical sequences in all colonizing populations indicating, most probably, that they are identical by descent. Moreover, Osvaldo is inserted at the same nucleotide site in all colonizing populations. On the other hand an almost identical LTR sequence, except by 1 base deletion, was found in the B4 clone compared to the canonical active Osvaldo element. These results suggest that Osvaldo copies in highly occupied sites are, most probably, identical by descent and strongly favor the founder hypothesis. On the other hand, low-insertion-frequency sites could represent recent transposition events. This work emphasizes the importance of molecular population studies to disentangle the effects of genetic drift and transposition in colonization.

Ongoing loss of the tirant transposable element in natural populations of Drosophila simulans

Tirant is a long terminal repeat (LTR) retrotransposon with an average of 11 insertion sites on the chromosome arms of Drosophila melanogaster flies collected from natural populations worldwide. In the sibling species Drosophila simulans, tirant is found only in African populations, which harbor a few insertion sites (1 to 5) on the chromosome arms, although some tirant sequences are present in the heterochromatin of most populations. This distribution in D. simulans reflects either the recent genomic invasion of African populations by a new variant of tirant, or a loss of tirant from the entire species apart from some sequence relics still present in Africa. In an attempt to clarify the situation, we focused on the LTR-UTR region of tirant copies from various populations of both D. melanogaster and D. simulans. We found two distinct types of regulatory region: one type was present in both D. melanogaster and D. simulans, and the other was present only in D. simulans. Copies of this latter type of tirant were transcriptionally inactive in gonads. Here we propose that the present day distribution of tirant in D. simulans populations reflects an ancient invasion of D. simulans by tirant copies followed by the loss of active copies from most populations, apart from the African ones, suggesting that this loss is still ongoing in this species.

Reversible introduction of transgenes in natural populations of insects

The most serious challenge concerning genetically modified insects remains their invasion ability. Indeed, transgenic insects often show lower fitness than wild individuals, and the transgene does not seem able to spread through a natural population without a driving system. The use of remobilizable vectors, based on the invading properties of transposable elements, has been frequently suggested. Simulations show that this strategy can be efficient. Moreover, if the transgene is designed to use transposition machinery already present in the genome, the transgene invasion appears to be potentially reversible after a few hundred generations, leading to new experimental perspectives.

The first steps of transposable elements invasion: parasitic strategy vs. genetic drift

Transposable elements are often considered as selfish DNA sequences able to invade the genome of their host species. Their evolutive dynamics are complex, due to the interaction between their intrinsic amplification capacity, selection at the host level, transposition regulation, and genetic drift. Here, we propose modeling the first steps of TE invasion, i.e., just after a horizontal transfer, when a single copy is present in the genome of one individual. If the element has a constant transposition rate, it will disappear in most cases: the elements with low-transposition rate are frequently lost through genetic drift, while those with high-transposition rate may amplify, leading to the sterility of their host. Elements whose transposition rate is regulated are able to successfully invade the populations, thanks to an initial transposition burst followed by a strong limitation of their activity. Self-regulation or hybrid dysgenesis may thus represent some genome-invasion parasitic strategies.

Insertional polymorphism of a non-LTR mobile element (NLRCth1) in European populations of Chironomus riparius (Diptera, Chironomidae) as detected by transposon insertion display

The midge Chironomus riparius is distributed all over the Palearctic region and is well characterized both at the morphological and cytogenetic levels. Here we describe a population study based on the insertional polymorphism of the retroposon NLRCth1, by means of a S-SAP (sequence-specific amplification polymorphism) derived technique (transposon insertion display; TID). While a previous study of allozyme polymorphism in Russian samples showed little variability, all the amplicons we identified are polymorphic. Genetic distances between 6 natural populations were calculated according to Nei and did not show a positive correlation with geographic distances. The genetic diversity detected among individuals of a given population was one order of magnitude higher than that among populations. However, the value of phi(ST) was significant (p < 0.001) and indicates that natural populations are more genetically differentiated than random samples of individuals.

Is genome size influenced by colonization of new environments in dipteran species?

Genome size differences are usually attributed to the amplification and deletion of various repeated DNA sequences, including transposable elements (TEs). Because environmental changes may promote modifications in the amount of these repeated sequences, it has been postulated that when a species colonizes new environments this could be followed by an increase in its genome size. We tested this hypothesis by estimating the genome size of geographically distinct populations of Drosophila ananassae, Drosophila malerkotliana, Drosophila melanogaster, Drosophila simulans, Drosophila subobscura, and Zaprionus indianus, all of which have known colonization capacities. There was no strong statistical differences between continents for most species. However, we found that populations of D. melanogaster from east Africa have smaller genomes than more recent populations. For species in which colonization is a recent event, the differences between genome sizes do not thus seem to be related to colonization history. These findings suggest either that genome size is seldom modified in a significant way during colonization or that it takes time for genome size of invading species to change significantly.

New regulatory regions of Drosophila 412 retrotransposable element generated by recombination

There are no doubts that transposable elements (TEs) have greatly influenced genomes evolution. They have, however, evolved in different ways throughout mammals, plants, and invertebrates. In mammals they have been shown to be widely present but with low transposition activity; in plants they are responsible for large increases in genome size. In Drosophila, despite their low amount, transposition seems to be higher. Therefore, to understand how these elements have evolved in different genomes and how host genomes have proposed to go around them, are major questions on genome evolution. We analyzed sequences of the retrotransposable elements 412 in natural populations of the Drosophila simulans and D. melanogaster species that greatly differ in their amount of TEs. We identified new subfamilies of this element that were the result of mutation or insertion-deletion process, but also of interfamily recombinations. These new elements were well conserved in the D. simulans natural populations. The new regulatory regions produced by recombination could give rise to new elements able to overcome host control of transposition and, thus, become potential genome invaders.