Reconstructing the Invasion Route of the P-Element in Drosophila melanogaster Using Extant Population Samples

On the origin of the P-element invasion in Drosophila simulans

The horizontal transfer (HT) of the P-element is one of the best documented cases of the HT of a transposable element. The P-element invaded natural D. melanogaster populations between 1950 and 1980 following its HT from Drosophila willistoni, a species endemic to South and Central America. Subsequently, it spread in D. simulans populations between 2006 and 2014, following a HT from D. melanogaster. The geographic region where the spread into D. simulans occurred is unclear, as both involved species are cosmopolitan. The P-element differs between these two species by a single base substitution at site 2040, where D. melanogaster carries a 'G' and D. simulans carries an 'A'. It has been hypothesized that this base substitution was a necessary adaptation that enabled the spread of the P-element in D. simulans, potentially explaining the 30-50-year lag between the invasions of D. melanogaster and D. simulans. To test this hypothesis, we monitored the invasion dynamics of P-elements with both alleles in experimental populations of D. melanogaster and D. simulans. Our results indicate that the allele at site 2040 has a minimal impact on the invasion dynamics of the P-element and, therefore, was not necessary for the invasion of D. simulans. However, we found that the host species significantly influenced the invasion dynamics, with higher P-element copy numbers accumulating in D. melanogaster than in D. simulans. Finally, based on SNPs segregating in natural D. melanogaster populations, we suggest that the horizontal transfer of the P-element from D. melanogaster to D. simulans likely occurred around Tasmania.

The Landscape of the DNA Transposons in the Genome of the Horezu_LaPeri Strain of Drosophila melanogaster

Natural transposons (NTs) represent mobile DNA sequences found in both prokaryotic and eukaryotic genomes. Drosophila melanogaster (the fruit fly) is a eukaryotic model organism with NTs standing for about 20% of its genome and has contributed significantly to the understanding of various aspects of transposon biology. Our study describes an accurate approach designed to map class II transposons (DNA transposons) in the genome of the Horezu_LaPeri fruit fly strain, consecutive to Oxford Nanopore Technology sequencing. A whole genome bioinformatics analysis was conducted using Genome ARTIST_v2, LoRTE and RepeatMasker tools to identify DNA transposons insertions. Then, a gene ontology enrichment analysis was performed in order to evaluate the potential adaptive role of some DNA transposons insertions. Herein, we describe DNA transposon insertions specific for the Horezu_LaPeri genome and a predictive functional analysis of some insertional alleles. The PCR validation of P-element insertions specific for this fruit fly strain, along with a putative consensus sequence for the KP element, is also reported. Overall, the genome of the Horezu_LaPeri strain contains several insertions of DNA transposons associated with genes known to be involved in adaptive processes. For some of these genes, insertional alleles obtained via mobilization of the artificial transposons were previously reported. This is a very alluring aspect, as it suggests that insertional mutagenesis experiments conducting adaptive predictions for laboratory strains may be confirmed by mirroring insertions which are expected to be found at least in some natural fruit fly strains.

The transposition rate has little influence on the plateauing level of the P-element

The popular trap model assumes that the invasions of transposable elements (TEs) in mammals and invertebrates are stopped by piRNAs that emerge after insertion of the TE into a piRNA cluster. It remains, however, still unclear which factors influence the dynamics of TE invasions. The activity of the TE (i.e., transposition rate) is one frequently discussed key factor. Here we take advantage of the temperature-dependent activity of the P-element, a widely studied eukaryotic TE, to test how TE activity affects the dynamics of a TE invasion. We monitored P-element invasion dynamics in experimental Drosophila simulans populations at hot and cold culture conditions. Despite marked differences in transposition rates, the P-element reached very similar copy numbers at both temperatures. The reduction of the insertion rate upon approaching the copy number plateau was accompanied by similar amounts of piRNAs against the P-element at both temperatures. Nevertheless, we also observed fewer P-element insertions in piRNA clusters than expected, which is not compatible with a simple trap model. The ping-pong cycle, which degrades TE transcripts, becomes typically active after the copy number plateaued. We generated a model, with few parameters, that largely captures the observed invasion dynamics. We conclude that the transposition rate has at the most only a minor influence on TE abundance, but other factors, such as paramutations or selection against TE insertions are shaping the TE composition.