The genome sequence of the grape phylloxera provides insights into the evolution, adaptation, and invasion routes of an iconic pest.

BACKGROUND: Although native to North America, the invasion of the aphid-like grape phylloxera Daktulosphaira vitifoliae across the globe altered the course of grape cultivation. For the past 150 years, viticulture relied on grafting-resistant North American Vitis species as rootstocks, thereby limiting genetic stocks tolerant to other stressors such as pathogens and climate change. Limited understanding of the insect genetics resulted in successive outbreaks across the globe when rootstocks failed. Here we report the 294-Mb genome of D. vitifoliae as a basic tool to understand host plant manipulation, nutritional endosymbiosis, and enhance global viticulture. RESULTS: Using a combination of genome, RNA, and population resequencing, we found grape phylloxera showed high duplication rates since its common ancestor with aphids, but similarity in most metabolic genes, despite lacking obligate nutritional symbioses and feeding from parenchyma. Similarly, no enrichment occurred in development genes in relation to viviparity. However, phylloxera evolved > 2700 unique genes that resemble putative effectors and are active during feeding. Population sequencing revealed the global invasion began from the upper Mississippi River in North America, spread to Europe and from there to the rest of the world. CONCLUSIONS: The grape phylloxera genome reveals genetic architecture relative to the evolution of nutritional endosymbiosis, viviparity, and herbivory. The extraordinary expansion in effector genes also suggests novel adaptations to plant feeding and how insects induce complex plant phenotypes, for instance galls. Finally, our understanding of the origin of this invasive species and its genome provide genetics resources to alleviate rootstock bottlenecks restricting the advancement of viticulture.

Correction to: The genome sequence of the grape phylloxera provides insights into the evolution, adaptation, and invasion routes of an iconic pest.

An amendment to this paper has been published and can be accessed via the original article.

Experimental evolution reveals hyperparasitic interactions among transposable elements.

Transposable elements (TEs) are repeated DNA sequences that can constitute a substantial part of genomes. Studying TEs' activity, interactions, and accumulation dynamics is thus of major interest to understand genome evolution. Here, we describe the transposition dynamics of cut-and-paste mariner elements during experimental (short- and longer-term) evolution in Drosophila melanogaster Flies with autonomous and nonautonomous mariner copies were introduced in populations containing no active mariner, and TE accumulation was tracked by quantitative PCR for up to 100 generations. Our results demonstrate that (i) active mariner elements are highly invasive and characterized by an elevated transposition rate, confirming their capacity to spread in populations, as predicted by the "selfish-DNA" mechanism; (ii) nonautonomous copies act as parasites of autonomous mariner elements by hijacking the transposition machinery produced by active mariner, which can be considered as a case of hyperparasitism; (iii) this behavior resulted in a failure of active copies to amplify which systematically drove the whole family to extinction in less than 100 generations. This study nicely illustrates how the presence of transposition-competitive variants can deeply impair TE dynamics and gives clues to the extraordinary diversity of TE evolutionary histories observed in genomes.

Diversity and evolution of mariner-like elements in aphid genomes.

BACKGROUND: Although transposons have been identified in almost all organisms, genome-wide information on mariner elements in Aphididae remains unknown. Genomes of Acyrthosiphon pisum, Diuraphis noxia and Myzus persicae belonging to the Macrosiphini tribe, actually available in databases, have been investigated. RESULTS: A total of 22 lineages were identified. Classification and phylogenetic analysis indicated that they were subdivided into three monophyletic groups, each of them containing at least one putative complete sequence, and several non-autonomous sublineages corresponding to Miniature Inverted-Repeat Transposable Elements (MITE), probably generated by internal deletions. A high proportion of truncated and dead copies was also detected. The three clusters can be defined from their catalytic site: (i) mariner DD34D, including three subgroups of the irritans subfamily (Macrosiphinimar, Batmar-like elements and Dnomar-like elements); (ii) rosa DD41D, found in A. pisum and D. noxia; (iii) a new clade which differs from rosa through long TIRs and thus designated LTIR-like elements. Based on its catalytic domain, this new clade is subdivided into DD40D and DD41D subgroups. Compared to other Tc1/mariner superfamily sequences, rosa DD41D and LTIR DD40-41D seem more related to maT DD37D family. CONCLUSION: Overall, our results reveal three clades belonging to the irritans subfamily, rosa and new LTIR-like elements. Data on structure and specific distribution of these transposable elements in the Macrosiphini tribe contribute to the understanding of their evolutionary history and to that of their hosts.

VHICA, a New Method to Discriminate between Vertical and Horizontal Transposon Transfer: Application to the Mariner Family within Drosophila.

Transposable elements (TEs) are genomic repeated sequences that display complex evolutionary patterns. They are usually inherited vertically, but can occasionally be transmitted between sexually independent species, through so-called horizontal transposon transfers (HTTs). Recurrent HTTs are supposed to be essential in life cycle of TEs, which are otherwise destined for eventual decay. HTTs also impact the host genome evolution. However, the extent of HTTs in eukaryotes is largely unknown, due to the lack of efficient, statistically supported methods that can be applied to multiple species sequence data sets. Here, we developed a new automated method available as a R package "vhica" that discriminates whether a given TE family was vertically or horizontally transferred, and potentially infers donor and receptor species. The method is well suited for TE sequences extracted from complete genomes, and applicable to multiple TEs and species at the same time. We first validated our method using Drosophila TE families with well-known evolutionary histories, displaying both HTTs and vertical transmission. We then tested 26 different lineages of mariner elements recently characterized in 20 Drosophila genomes, and found HTTs in 24 of them. Furthermore, several independent HTT events could often be detected within the same mariner lineage. The VHICA (Vertical and Horizontal Inheritance Consistence Analysis) method thus appears as a valuable tool to analyze the evolutionary history of TEs across a large range of species.

Transcriptional polymorphism of piRNA regulatory genes underlies the mariner activity in Drosophila simulans testes.

During colonization of new areas, natural populations have to deal with changing environments, and transposable elements (TEs) can be useful "tools" in the adaptation process as they are major contributor to the structural and functional evolution of genomes. In this general context, the activity (copy number, transcriptional and excision rate) of the mariner mos1 element was estimated in 19 natural populations of D. simulans. It is shown (i) that mos1 expression is always higher and more variable in testes than in ovaries; (ii) that mos1 activity is higher in colonizing populations compared to the sub-Saharan African ones (ancestral populations); (iii) that mos1 variations in transcript levels and copy number are negatively correlated with transcriptional variations of piRNA genes, aubergine and argonaute3. Furthermore, mos1 levels of expression in testes highly contrast with the low expression patterns of ago3. These results strongly suggest that the expression polymorphism of piRNA genes could be responsible for the mos1 variations, first between male and female germlines and second, according to the status of natural populations (colonizing or not). These results provide new perspectives about TEs and piRNA genes co-evolution in Drosophila germlines.

Characterization of irritans mariner-like elements in the olive fruit fly Bactrocera oleae (Diptera: Tephritidae): evolutionary implications.

Genomic variation among species is commonly driven by transposable element (TE) invasion; thus, the pattern of TEs in a genome allows drawing an evolutionary history of the studied species. This paper reports in vitro and in silico detection and characterization of irritans mariner-like elements (MLEs) in the genome and transcriptome of Bactrocera oleae (Rossi) (Diptera: Tephritidae). Eleven irritans MLE sequences have been isolated in vitro using terminal inverted repeats (TIRs) as primers, and 215 have been extracted in silico from the sequenced genome of B. oleae. Additionally, the sequenced genomes of Bactrocera tryoni (Froggatt) and Bactrocera cucurbitae (Diptera: Tephritidae) have been explored to identify irritans MLEs. A total of 129 sequences from B. tryoni have been extracted, while the genome of B. cucurbitae appears probably devoid of irritans MLEs. All detected irritans MLEs are defective due to several mutations and are clustered together in a monophyletic group suggesting a common ancestor. The evolutionary history and dynamics of these TEs are discussed in relation with the phylogenetic distribution of their hosts. The knowledge on the structure, distribution, dynamic, and evolution of irritans MLEs in Bactrocera species contributes to the understanding of both their evolutionary history and the invasion history of their hosts. This could also be the basis for genetic control strategies using transposable elements.

S-palmitoylation and s-oleoylation of rabbit and pig sarcolipin.

Sarcolipin (SLN) is a regulatory peptide present in sarcoplasmic reticulum (SR) from skeletal muscle of animals. We find that native rabbit SLN is modified by a fatty acid anchor on Cys-9 with a palmitic acid in about 60% and, surprisingly, an oleic acid in the remaining 40%. SLN used for co-crystallization with SERCA1a (Winther, A. M., Bublitz, M., Karlsen, J. L., Moller, J. V., Hansen, J. B., Nissen, P., and Buch-Pedersen, M. J. (2013) Nature 495, 265-2691; Ref. 1) is also palmitoylated/oleoylated, but is not visible in crystal structures, probably due to disorder. Treatment with 1 m hydroxylamine for 1 h removes the fatty acids from a majority of the SLN pool. This treatment did not modify the SERCA1a affinity for Ca(2+) but increased the Ca(2+)-dependent ATPase activity of SR membranes indicating that the S-acylation of SLN or of other proteins is required for this effect on SERCA1a. Pig SLN is also fully palmitoylated/oleoylated on its Cys-9 residue, but in a reverse ratio of about 40/60. An alignment of 67 SLN sequences from the protein databases shows that 19 of them contain a cysteine and the rest a phenylalanine at position 9. Based on a cladogram, we postulate that the mutation from phenylalanine to cysteine in some species is the result of an evolutionary convergence. We suggest that, besides phosphorylation, S-acylation/deacylation also regulates SLN activity.

Characterization of mariner-like transposons of the mauritiana Subfamily in seven tree aphid species.

Mariner-like elements (MLEs) are Class II transposons present in all eukaryotic genomes in which MLEs have been searched for. This article reports the detection of MLEs in seven of the main fruit tree aphid species out of eight species studied. Deleted MLE sequences of 916-919 bp were characterized, using the terminal-inverted repeats (TIRs) of mariner elements belonging to the mauritiana Subfamily as primers. All the sequences detected were deleted copies of full-length elements that included the 3'- and 5'-TIRs but displayed internal deletions affecting Mos1 activity. Networks based on the mtDNA cytochrome oxidase subunit-I (CO-I) and MLE sequences were incongruent, suggesting that mutations in transposon sequences had accumulated before speciation of tree aphid species occurred, and that they have been maintained in this species via vertical transmissions. This is the first evidence of the widespread occurrence of MLEs in aphids.

Genomic landscape and evolutionary dynamics of mariner transposable elements within the Drosophila genus.

BACKGROUND: The mariner family of transposable elements is one of the most widespread in the Metazoa. It is subdivided into several subfamilies that do not mirror the phylogeny of these species, suggesting an ancient diversification. Previous hybridization and PCR studies allowed a partial survey of mariner diversity in the Metazoa. In this work, we used a comparative genomics approach to access the genus-wide diversity and evolution of mariner transposable elements in twenty Drosophila sequenced genomes. RESULTS: We identified 36 different mariner lineages belonging to six distinct subfamilies, including a subfamily not described previously. Wide variation in lineage abundance and copy number were observed among species and among mariner lineages, suggesting continuous turn-over. Most mariner lineages are inactive and contain a high proportion of damaged copies. We showed that, in addition to substitutions that rapidly inactivate copies, internal deletion is a major mechanism contributing to element decay and the generation of non-autonomous sublineages. Hence, 23% of copies correspond to several Miniature Inverted-repeat Transposable Elements (MITE) sublineages, the first ever described in Drosophila for mariner. In the most successful MITEs, internal deletion is often associated with internal rearrangement, which sheds light on the process of MITE origin. The estimation of the transposition rates over time revealed that all lineages followed a similar progression consisting of a rapid amplification burst followed by a rapid decrease in transposition. We detected some instances of multiple or ongoing transposition bursts. Different amplification times were observed for mariner lineages shared by different species, a finding best explained by either horizontal transmission or a reactivation process. Different lineages within one species have also amplified at different times, corresponding to successive invasions. Finally, we detected a preference for insertion into short TA-rich regions, which appears to be specific to some subfamilies. CONCLUSIONS: This analysis is the first comprehensive survey of this family of transposable elements at a genus scale. It provides precise measures of the different evolutionary processes that were hypothesized previously for this family based on PCR data analysis. mariner lineages were observed at almost all "life cycle" stages: recent amplification, subsequent decay and potential (re)-invasion or invasion of genomes.

Fossil rhabdoviral sequences integrated into arthropod genomes: ontogeny, evolution, and potential functionality.

Retroelements represent a considerable fraction of many eukaryotic genomes and are considered major drives for adaptive genetic innovations. Recent discoveries showed that despite not normally using DNA intermediates like retroviruses do, Mononegaviruses (i.e., viruses with nonsegmented, negative-sense RNA genomes) can integrate gene fragments into the genomes of their hosts. This was shown for Bornaviridae and Filoviridae, the sequences of which have been found integrated into the germ line cells of many vertebrate hosts. Here, we show that Rhabdoviridae sequences, the major Mononegavirales family, have integrated only into the genomes of arthropod species. We identified 185 integrated rhabdoviral elements (IREs) coding for nucleoproteins, glycoproteins, or RNA-dependent RNA polymerases; they were mostly found in the genomes of the mosquito Aedes aegypti and the blacklegged tick Ixodes scapularis. Phylogenetic analyses showed that most IREs in A. aegypti derived from multiple independent integration events. Since RNA viruses are submitted to much higher substitution rates as compared with their hosts, IREs thus represent fossil traces of the diversity of extinct Rhabdoviruses. Furthermore, analyses of orthologous IREs in A. aegypti field mosquitoes sampled worldwide identified an integrated polymerase IRE fragment that appeared under purifying selection within several million years, which supports a functional role in the host's biology. These results show that A. aegypti was subjected to repeated Rhabdovirus infectious episodes during its evolution history, which led to the accumulation of many integrated sequences. They also suggest that like retroviruses, integrated rhabdoviral sequences may participate actively in the evolution of their hosts.

Genomic parasites or symbionts? Modeling the effects of environmental pressure on transposition activity in asexual populations.

Transposable elements are DNA segments capable of persisting in host genomes by self-replication in spite of deleterious mutagenic effects. The theoretical dynamics of these elements within genomes has been studied extensively, and population genetic models predict that they can invade and maintain as a result of both intra-genomic and inter-individual selection in sexual species. In asexuals, the success of selfish DNA is more difficult to explain. However, most theoretical work assumes constant environment. Here, we analyze the impact of environmental change on the dynamics of transposition activity when horizontal DNA exchange is absent, based on a stochastic computational model of transposable element proliferation. We argue that repeated changes in the phenotypic optimum in a multidimensional fitness landscape may induce explosive bursts of transposition activity associated with faster adaptation. However, long-term maintenance of transposition activity is unlikely. This could contribute to the significant variation in the transposable element copy number among closely related species.

Internal deletions of transposable elements: the case of Lemi elements.

Mobile elements using a "cut and paste" mechanism of transposition (Class II) are frequently prone to internal deletions and the question of the origin of these copies remains elusive. In this study, we looked for copies belonging to the Lemi Family (Tc1-mariner-IS630 SuperFamily) in the plant genomes, and copies within internal deletions were analyzed in detail. Lemi elements are found exclusively in Eudicots, and more than half of the copies have been deleted. All deletions occur between microhomologies (direct repeats from 2 to 13 bp). Copies less than 500 bp long, similar to MITEs, are frequent. These copies seem to result from large deletions occurring between microhomologies present within a region of 300 bp at both extremities of the element. These regions are particularly A/T rich, compared to the internal part of the element, which increases the probability of observing short direct repeats. Most of the molecular mechanisms responsible for double strand break repair are able to induce deletions between microhomologies during the repair process. This could be a quick way to reduce the population of active copies within a genome and, more generally, to reduce the overall activity of the element after it has entered a naive genome.

A unified classification system for eukaryotic transposable elements.

Our knowledge of the structure and composition of genomes is rapidly progressing in pace with their sequencing. The emerging data show that a significant portion of eukaryotic genomes is composed of transposable elements (TEs). Given the abundance and diversity of TEs and the speed at which large quantities of sequence data are emerging, identification and annotation of TEs presents a significant challenge. Here we propose the first unified hierarchical classification system, designed on the basis of the transposition mechanism, sequence similarities and structural relationships, that can be easily applied by non-experts. The system and nomenclature is kept up to date at the WikiPoson web site.

Copia retrotransposon in the Zaprionus genus: another case of transposable element sharing with the Drosophila melanogaster subgroup.

Copia is a retrotransposon that appears to be distributed widely among the Drosophilidae subfamily. Evolutionary analyses of regulatory regions have indicated that the Copia retrotransposon evolved through both positive and purifying selection, and that horizontal transfer (HT) could also explain its patchy distribution of the among the subfamilies of the melanogaster subgroup. Additionally, Copia elements could also have transferred between melanogaster subgroup and other species of Drosophilidae-D. willistoni and Z. tuberculatus. In this study, we surveyed seven species of the Zaprionus genus by sequencing the LTR-ULR and reverse transcriptase regions, and by using RT-PCR in order to understand the distribution and evolutionary history of Copia in the Zaprionus genus. The Copia element was detected, and was transcriptionally active, in all species investigated. Structural and selection analysis revealed Zaprionus elements to be closely related to the most ancient subfamily of the melanogaster subgroup, and they seem to be evolving mainly under relaxed purifying selection. Taken together, these results allowed us to classify the Zaprionus sequences as a new subfamily-ZapCopia, a member of the Copia retrotransposon family of the melanogaster subgroup. These findings indicate that the Copia retrotransposon is an ancient component of the genomes of the Zaprionus species and broaden our understanding of the diversity of retrotransposons in the Zaprionus genus.

The evolutionary history of mariner-like elements in Neotropical drosophilids.

The evolutionary history of mariner-like elements (MLEs) in 49 mainly Neotropical drosophilid species is described. So far, the investigations about the distribution of MLEs were performed mainly using hybridization assays with the Mos1 element (the first mariner active element described) in a widely range of drosophilid species and these sequences were found principally in species that arose in Afrotropical and Sino-Indian regions. Our analysis in mainly Neotropical drosophilid species shows that twenty-three species presented MLEs from three different subfamilies in their genomes: eighteen species had MLEs from subfamily mellifera, fifteen from subfamily mauritiana and three from subfamily irritans. Eleven of these species exhibited elements from more than one subfamily in their genome. In two subfamilies, the analyzed coding region was uninterrupted and contained conserved catalytic motifs. This suggests that these sequences were probably derived from active elements. The species with these putative active elements are Drosophila mediopunctata and D. busckii for the mauritiana subfamily, and D. paramediostriata for the mellifera subfamily. The phylogenetic analysis of MLE, shows a complex evolutionary pattern, exhibiting vertical transfer, stochastic loss and putative events of horizontal transmission occurring between different Drosophilidae species, and even those belonging to more distantly related taxa such as Bactrocera tryoni (Tephritidae family), Sphyracephala europaea (Diopsoidea superfamily) and Buenoa sp. (Hemiptera order). Moreover, our data show that the distribution of MLEs is not restricted to Afrotropical and Sino-Indian species. Conversely, these TEs are also widely distributed in drosophilid species arisen in the Neotropical region.

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.

Long-term evolution of transposable elements.

Transposable elements are often considered parasitic DNA sequences, able to invade the genome of their host thanks to their self-replicating ability. This colonization process has been extensively studied, both theoretically and experimentally, but their long-term coevolution with the genomes is still poorly understood. In this work, we aim to challenge previous population genetics models by considering features of transposable elements as quantitative, rather than discrete, variables. We also describe more realistic transposable element dynamics by accounting for the variability of the insertion effect, from deleterious to adaptive, as well as mutations leading to a loss of transposition activity and to nonautonomous copies. Individual-based simulations of the behavior of a transposable-element family over several thousand generations show different ways in which active or inactive copies can be maintained for a very long time. Results reveal an unexpected impact of genetic drift on the "junk DNA" content of the genome and strongly question the likelihood of the sustainable long-term stable transposition-selection equilibrium on which numerous previous works were based.

piRNA and Transposon Dynamics in Drosophila: A Female Story.

The germlines of metazoans contain transposable elements (TEs) causing genetic instability and affecting fitness. To protect the germline from TE activity, gonads of metazoans produce TE-derived PIWI-interacting RNAs (piRNAs) that silence TE expression. In Drosophila, our understanding of piRNA biogenesis is mainly based on studies of the Drosophila melanogaster female germline. However, it is not known whether piRNA functions are also important in the male germline or whether and how piRNAs are affected by the global genomic context. To address these questions, we compared genome sequences, transcriptomes, and small RNA libraries extracted from entire testes and ovaries of two sister species: D. melanogaster and Drosophila simulans. We found that most TE-derived piRNAs were produced in ovaries and that piRNA pathway genes were strongly overexpressed in ovaries compared with testes, indicating that the silencing of TEs by the piRNA pathway mainly took place in the female germline. To study the relationship between host piRNAs and TE landscape, we analyzed TE genomic features and how they correlate with piRNA production in the two species. In D. melanogaster, we found that TE-derived piRNAs target recently active TEs. In contrast, although Drosophila simulans TEs do not display any features of recent activity, the host still intensively produced silencing piRNAs targeting old TE relics. Together, our results show that the piRNA silencing response mainly takes place in Drosophila ovaries and indicate that the host piRNA response is implemented following a burst of TE activity and could persist long after the extinction of active TE families.