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.

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.

The somatic mobilization of transposable element mariner-Mos1 during the Drosophila lifespan and its biological consequences

Transposable elements (TEs) are mobile DNA sequences on genomes. Some elements are able to transpose in somatic cells, a process known as somatic transposition (ST), which has been associated with detrimental biological effects. The mariner-Mos1 element of Drosophila promotes transposition in somatic and germline cells and is an excellent model for studies related to the biological consequence of somatic excision (SE). In this work, we used temperature stress to induce increasing transposition of mariner-Mos1 during different stages of the development of D. simulans, aiming to quantify SE during lifespan. Furthermore, strains of D. melanogaster exhibiting differential expression of mariner-Mos1 were employed for estimating some biological consequences of mariner mobilization. It is shown that SE of mariner-Mos1 was not constant during development; the larval phase had the highest rates while the pupal stage exhibited lower rates, and in the embryonic stage, no difference was detected. SE can be detrimental, as suggested by correlation in SE level and reduction in behavioral activities and embryonic viability. This study showed that mariner-Mos1 SE accumulates during the Drosophila life cycle, and can be involved in detrimental effects.

Molecular Evolution of piggyBac Superfamily: From Selfishness to Domestication

The piggyBac transposable element was originally isolated from the cabbage looper moth, Trichoplusia ni, in the 1980s. Despite its early discovery and specificity compared to the other Class II elements, the diversity and evolution of this superfamily have been only partially analyzed. Two main types of elements can be distinguished: the piggyBac-like elements (PBLE) with terminal inverted repeats, untranslated region, and an open reading frame encoding a transposase, and the piggyBac-derived sequences (PGBD), containing a sequence derived from a piggyBac transposase, and which correspond to domesticated elements. To define the distribution, their structural diversity and phylogenetic relationships, analyses were conducted using known PBLE and PGBD sequences to scan databases. From this data mining, numerous new sequences were characterized (50 for PBLE and 396 for PGBD). Structural analyses suggest that four groups of PBLE can be defined according to the presence/absence of sub-terminal repeats. The transposase is characterized by highly variable catalytic domain and C-terminal region. There is no relationship between the structural groups and the phylogeny of these PBLE elements. The PGBD are clearly structured into nine main groups. A new group of domesticated elements is suspected in Neopterygii and the remaining eight previously described elements have been investigated in more detail. In all cases, these sequences are no longer transposable elements, the catalytic domain of the ancestral transposase is not always conserved, but they are under strong purifying selection. The phylogeny of both PBLE and PGBD suggests multiple and independent domestication events of PGBD from different PBLE ancestors.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3856-6) contains supplementary material, which is available to authorized users.

ALCOHOL TOLERANCE, ADH ACTIVITY, AND ECOLOGICAL NICHE OF DROSOPHILA SPECIES

In vitro alcohol dehydrogenase (ADH) activity was measured in adults of species belonging to Drosophila and to the related genus Zaprionus. Data were analyzed according to the known breeding sites and the level of ethanol tolerance of these species. Alcohol dehydrogenase activity was assayed with both ethanol (E) and isopropanol (I). Our results show a very broad range of activities among the 71 species investigated, the ratio of the highest value observed (D. melanogaster) to the lowest (D. pruinosa) being 65:1. A general positive correlation was found between the level of ADH activity and the capacity to detoxify ethanol. Nevertheless, many species show exceptions to this rule. Contrary to a logical expectation, adaptation to high alcoholic resources, which has been a recurrent evolutionary event, was not mediated by a more efficient use of ethanol, that is, an increase of the E/I ratio. This ratio seems to be quite variable according to the phylogeny and is especially low in the subgenus Sophophora as well as in Zaprionus. Alcohol tolerance clearly is related to the larval habitat of the species and shows that adaptation to alcoholic resources has been a major evolutionary challenge in drosophilids. This adaptation is not related to phylogeny, having occurred independently several times during the evolution of the group. Finally, it should be borne in mind that, besides metabolization and detoxification, other physiological processes such as nervous‐system tolerance or ethanol excretion may be involved in ethanol tolerance, and such functions also should be investigated. Environmental ethanol, which is certainly a major ecological parameter for many drosophilids, has selected a diversity of physiological adaptations, all related to the Adh locus, but presumably much more complicated than was previously believed.

ADAPTATION TO FERMENTING RESOURCES IN DROSOPHILA MELANOGASTER: ETHANOL AND ACETIC ACID TOLERANCES SHARE A COMMON GENETIC BASIS

Ethanol and acetic acid tolerances were compared in a French, highly tolerant population, and in a Congolese, very sensitive population. For both tolerances, chromosome substitutions demonstrated a major effect on chromosome 3, a lesser effect on chromosome 2, and no effect on chromosome 1, except in interactions. Directional selection experiments led to significant increases of tolerance to both toxics. Of greater interest, a strong correlated response was observed in each line: increased ethanol tolerance was accompanied by higher acetic acid tolerance and vice versa. A high genetic correlation (average value r = 0.77) was found between the two traits. These data suggest that alcohol dehydrogenase (ADH) activity does not play a major role in explaining the physiological differences known between Afrotropical and European populations. The metabolic flux permitting the detoxification of ethanol and acetic acid seems to be mainly controlled by acetyl-coA synthetase (ACS) at least in adult flies. Acetic acid adaptation could be as important as ethanol adaptation in the ecology of Drosophila melanogaster and other Drosophila species.

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.

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.

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.

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.

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.

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.

How does selfing affect the dynamics of selfish transposable elements?

Background

Many theoretical models predicting the dynamics of transposable elements (TEs) in genomes, populations, and species have already been proposed. However, most of them only focus on populations of sexual diploid individuals, and TE dynamics in populations partly composed by autogamous individuals remains poorly investigated. To estimate the impact of selfing on TE dynamics, the short- and long-term evolution of TEs was simulated in outcrossing populations with various proportions of selfing individuals.

Results

Selfing has a deep impact on TE dynamics: the higher the selfing rate, the lower the probability of invasion. Already known non-equilibrium dynamics (complete loss, domestication, cyclical invasion of TEs) can all be described whatever the mating system. However, their pattern and their respective frequencies greatly depend on the selfing rate. For instance, in cyclical dynamics resulting from interactions between autonomous and non-autonomous copies, cycles are faster when the selfing rate increases. Interestingly, an abrupt change in the mating system from sexuality to complete asexuality leads to the loss of all the elements over a few hundred generations. In general, for intermediate selfing rates, the transposition activity remains maintained.

Conclusions

Our theoretical results evidence that a clear and systematic contrast in TE content according to the mating system is expected, with a smooth transition for intermediate selfing rates. Several parameters impact the TE copy number, and all dynamics described in allogamous populations can be also observed in partly autogamous species. This study thus provides new insights to understand the complex signal from empirical comparison of closely related species with different mating systems.

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.

The struggle for life of the genome's selfish architects

Transposable elements (TEs) were first discovered more than 50 years ago, but were totally ignored for a long time. Over the last few decades they have gradually attracted increasing interest from research scientists. Initially they were viewed as totally marginal and anecdotic, but TEs have been revealed as potentially harmful parasitic entities, ubiquitous in genomes, and finally as unavoidable actors in the diversity, structure, and evolution of the genome. Since Darwin's theory of evolution, and the progress of molecular biology, transposable elements may be the discovery that has most influenced our vision of (genome) evolution. In this review, we provide a synopsis of what is known about the complex interactions that exist between transposable elements and the host genome. Numerous examples of these interactions are provided, first from the standpoint of the genome, and then from that of the transposable elements. We also explore the evolutionary aspects of TEs in the light of post-Darwinian theories of evolution.

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.

Applying mobile genetic elements for genome analysis and evolution

Transposable elements (TEs) are ubiquitous components of all living organisms, and in the course of their coexistence with their respective host genomes, these parasitc DNAs have played important roles in the evolution of complex genetic networks. The interaction between mobile DNAs and their host genomes are quite diverse, ranging from modifications of gene structure and regulation to alterations in general genome architecture. Thus during evolutionary time these elements can be regarded as natural molecular tools in shaping the organization, structure, and function of eukaryotic genes and genomes. Based on their intrinsic properties and features, mobile DNAs are widely applied at present as a technical "toolbox," essential for studying a diverse spectrum of biological questions. In this review, we aim to summarize both the evolutionary impact of TEs on genome evolution and their valuable and diverse methodological applications as molecular tools.

Multiple invasions of Gypsy and Micropia retroelements in genus Zaprionus and melanogaster subgroup of the genus Drosophila

Background

The Zaprionus genus shares evolutionary features with the melanogaster subgroup, such as space and time of origin. Although little information about the transposable element content in the Zaprionus genus had been accumulated, some of their elements appear to be more closely related with those of the melanogaster subgroup, indicating that these two groups of species were involved in horizontal transfer events during their evolution. Among these elements, the Gypsy and the Micropia retroelements were chosen for screening in seven species of the two Zaprionus subgenera, Anaprionus and Zaprionus.

Results

Screening allowed the identification of diverse Gypsy and Micropia retroelements only in species of the Zaprionus subgenus, showing that they are transcriptionally active in the sampled species. The sequences of each retroelement were closely related to those of the melanogaster species subgroup, and the most parsimonious hypothesis would be that 15 horizontal transfer events shaped their evolution. The Gypsy retroelement of the melanogaster subgroup probably invaded the Zaprionus genomes about 11 MYA. In contrast, the Micropia retroelement may have been introduced into the Zaprionus subgenus and the melanogaster subgroup from an unknown donor more recently (~3 MYA).

Conclusion

Gypsy and Micropia of Zaprionus and melanogaster species share similar evolutionary patterns. The sharing of evolutionary, ecological and ethological features probably allowed these species to pass through a permissive period of transposable element invasion, explaining the proposed waves of horizontal transfers.

Evolutionary genetics of Zaprionus. II. Mitochondrial DNA and chromosomal variation of the invasive drosophilid Zaprionus indianus in Egypt

Zaprionus indianus is an Afrotropical drosophilid species that has expanded its geographical range in the Palearctic region and the Americas during the second half of the last century. It has invaded Egypt within the past two decades from East Africa or Asia and became a dominant species in the drosophilid fauna therein, but the exact date of introduction and source of the propagule remain unknown. Here, we investigate the genetic structure of eight geographical populations within and around the Nile Delta using mitochondrial DNA sequences of the cox2 gene and chromosomal inversion polymorphism. A very low level of genetic variability was detected for both markers, mainly attributed to the introduction bottleneck. Nonetheless, both indicate a significant population structure, with a southeastern-northwestern cline. Demographic history analysis suggested northwestern populations to be younger (expanding in ca. 1992) than southeastern ones (expanding in ca. 1985). The In(II)A polymorphism was only observed in the northwestern population, but one-year interval analysis of the Alexandria population revealed the lack of seasonal fluctuation and a trend toward the loss of the polymorphism. Based on these data and faunistic records, we propose a multiple introduction scenario for Z. indianus in Egypt-according to which a first wave in the early 1980s from Sudan through normal northward range expansion or fruit trade, and a second wave in the early 1990s from Asia via fruit trade. We also suggest, from ecological observations, fruit trade data and known adaptive versatility of Z. indianus, date palm, the dominant fruit in Egypt and in the oases where Z. indianus predominates, to play a major role in the spread of the species in the Middle East.

Analysis of the DDE motif in the Mutator superfamily

The eukaryotic Mutator family of transposable elements is widespread in plants. Active or potentially active copies are also found in fungi and protozoans, and sequences related to this family have been detected in metazoans as well. Members of this family are called Mutator-like elements (MULEs). They encode transposases, which contain a region conserved with transposases of the IS256 prokaryotic family, known to harbor a DDE catalytic domain. Different DDE or D34E motifs have been proposed in some groups of eukaryotic MULEs based on primary sequence conservation. On a large number of protein sequences related to, and representative of, all MULE families, we analyzed global conservation, the close environment of different acidic residues and the secondary structure. This allowed us to identify a potential DDE motif that is likely to be homologous to the one in IS256-like transposases. The characteristics of this motif are depicted in each known family of MULEs. Different hypotheses about the evolution of this triad are discussed.

The mariner transposable element in natural populations of Drosophila simulans

In cosmopolitan species, geographical variations in copy number and/or level of transposition activity have been observed for several transposable elements (TEs). Environment, history and population structure can contribute to such variation in ways that are difficult to tease apart. For the mariner element, previous studies of the geographic variation of its somatic activity in natural populations of Drosophila simulans have shown contradictory results (latitudinal clines of divergent orientations or no apparent structure). To try and resolve these inconsistencies, we gathered all available data on the mariner somatic activity of worldwide natural populations. This includes previously published results by different groups and also new data. The correlations between the level of activity and several geoclimatic factors were tested. Although no general effect of temperature was found, a relationship with the invasion history was detected. It was also shown that recent invasive populations have a higher level of activity than the putative ancestral ones. Our results strongly suggest that variability of the mariner somatic activity among natural populations of D. simulans is mainly due to populational and historical factors probably related to the recent world colonization of this species. Indeed, this activity is correlated to the main route out of Africa (the Nile route) and the recent colonization of continents such as Australia and South America.

DNA barcode discovers two cryptic species and two geographical radiations in the invasive drosophilid Zaprionus indianus

Comparing introduced to ancestral populations within a phylogeographical context is crucial in any study aiming to understand the ecological genetics of an invasive species. Zaprionus indianus is a cosmopolitan drosophilid that has recently succeeded to expand its geographical range upon three continents (Africa, Asia and the Americas). We studied the distribution of mitochondrial DNA (mtDNA) haplotypes for two genes (CO-I and CO-II) among 23 geographical populations. mtDNA revealed the presence of two well-supported phylogenetic lineages (phylads), with bootstrap value of 100%. Phylad I included three African populations, reinforcing the African-origin hypothesis of the species. Within phylad II, a distinct phylogeographical pattern was discovered: Atlantic populations (from the Americas and Madeira) were closer to the ancestral African populations than to Eastern ones (from Madagascar, Middle East and India). This means that during its passage from endemism to cosmopolitanism, Z. indianus exhibited two independent radiations, the older (the Eastern) to the East, and the younger (the Atlantic) to the West. Discriminant function analysis using 13 morphometrical characters was also able to discriminate between the two molecular phylads (93.34 ± 1.67%), although detailed morphological analysis of male genitalia using scanning electron microscopy showed no significant differences. Finally, crossing experiments revealed the presence of reproductive barrier between populations from the two phylads, and further between populations within phylad I. Hence, a bona species status was assigned to two new, cryptic species: Zaprionus africanus and Zaprionus gabonicus, and both were encompassed along with Z. indianus and Zaprionus megalorchis into the indianus complex. The ecology of these two species reveals that they are forest dwellers, which explains their restricted endemic distribution, in contrast to their relative cosmopolitan Z. indianus, known to be a human-commensal. Our results reconfirm the great utility of mtDNA at both inter- and intraspecific analyses within the frame of an integrated taxonomical project.

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.

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.

Mesosternal bristle number in a cosmopolitan drosophilid: an X-linked variable trait independent of sternopleural bristles

Mesosternal (MS) bristles in Drosophila are a pair of machrochaetae found at the sternal end of the sternopleural (STP) microchaetae, and are thought to be invariable. In a closely related drosophilid genus, Zaprionus, their number is four and, in contrast to Drosophila, they show interspecific and intraspecific variability. The genetic basis of MS bristle number variability was studied in Z. indianus, the only cosmopolitan species of the genus. The trait responded rapidly to selection and two lines were obtained, one lacking any bristles (0-0) and the other bearing the normal phenotype (2-2). Other symmetrical phenotypes, (1-1) and (3-3), could also be selected for, but with lesser success. By contrast, STP bristle number did not vary significantly between the two lines (0-0) and (2-2), revealing its genetic independence from MS bristle number. Reciprocal crosses between these two lines showed that MS bristle number is mainly influenced by a major gene on the X chromosome (i.e. F(1) males always resembled their mothers) with codominant expression (i.e. heterozygous F(1) females harboured an average phenotype of 2 bristles). However, trait penetrance was incomplete and backcrosses revealed that this variability was partly due to genetic modifiers, most likely autosomal. The canalization of MS bristle number was investigated under different temperatures, and the increased appearance of abnormal phenotypes mainly occurred at extreme temperatures. There was a bias, however, towards bristle loss, as shown by a liability (developmental map) analysis. Finally, when ancestral and introduced populations were compared, the latter were far less stable, suggesting that genetic bottlenecks may perturb the MS bristle number canalization system. MS bristle number, thus, appears to be an excellent model for investigating developmental canalization at both the quantitative and the molecular level.

Highly similar piggyBac transposase-like sequences from various Bactrocera (Diptera, Tephritidae) species

The piggyBac transposable element is currently the vector of choice for transgenesis, enhancer trapping, gene discovery and gene function determination in both insects and mammals. However, the recent discovery of sequences with similarity to piggyBac in a wide diversity of organisms suggests that piggyBac may be horizontally transferred to distantly related species. This has raised concern on the wide-range application of piggyBac-based transformation vectors and their stability. In this paper, the presence of sequences homologous to the piggyBac transposase was investigated in 17 species belonging to six genera within the Tephritidae family, including many pest species for which transformation has already been achieved. piggyBac-like sequences, with a high degree of similarity to the original Trichoplusia ni transposase sequence were identified only in six species of the Bactrocera genus.

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.

Genome ecosystem and transposable elements species

Transposable elements are known to be "selfish DNA" sequences able to spread and be maintained in all genomes analyzed so far. Their evolution depends on the interaction they have with the other components of the genome, including genes and other transposable elements. These relationships are complex and have often been compared to those of species living and competing in an ecosystem. The aim of this current work is a proposition to fill the conceptual gap existing between genome biology and ecology, assuming that genomic components, such as transposable elements families, can be compared to species interacting in an ecosystem. Using this framework, some of the main models defined in the population genetics of transposable elements can then been reformulated, and some new kinds of realistic relationships, such as symbiosis between different genomic components, can then be modelled and explored.

Population genetics models of competition between transposable element subfamilies

Transposable elements are one of the major components of genomes. Some copies are fully efficient; i.e., they are able to produce the proteins needed for their own transposition, and they can move and duplicate into the genome. Other copies are mutated. They may have lost their moving ability, their coding capacity, or both, thus becoming pseudogenes slowly eliminated from the genome through deletions and natural selection. Little is known about the dynamics of such mutant elements, particularly concerning their interactions with autonomous copies. To get a better understanding of the transposable elements' evolution after their initial invasion, we have designed a population genetics model of transposable elements dynamics including mutants or nonfunctional sequences. We have particularly focused on the case where these sequences are nonautonomous elements, known to be able to use the transposition machinery produced by the autonomous ones. The results show that such copies generally prevent the system from achieving a stable transposition-selection equilibrium and that nonautonomous elements can invade the system at the expense of autonomous ones. The resulting dynamics are mainly cyclic, which highlights the similarities existing between genomic selfish DNA sequences and host-parasite systems.

Sexual dimorphism of body size and sternopleural bristle number: a comparison of geographic populations of an invasive cosmopolitan drosophilid

Zaprionus indianus is a cosmopolitan drosophilid, of Afrotropical origin, which has recently colonized South America. The sexual dimorphism (SD) of body size is low, males being almost as big as females. We investigated 10 natural populations, 5 from America and 5 from Africa, using the isofemale line technique. Three traits were measured on each fly: wing and thorax length and sternopleural bristle number. Two indices of SD were compared, and found to be highly correlated (r > 0.99). For the sake of simplicity, only the female/male (F/M) ratio was further considered. A significant genetic variability of SD was found in all cases, although with a low heritability (intra-class correlation of 0.13), about half the value found for the traits themselves. For size SD, we did not find any variation among continents or any latitudinal trend, and average values were 1.02 for wing length and 1.01 for thorax length. Bristle number SD was much greater (1.07). Among mass laboratory strains, SD was genetically much more variable than in recently collected populations, a likely consequence of laboratory drift. Altogether, SD, although genetically variable and prone to laboratory drift, is independent of size variations and presumably submitted to a stabilizing selection in nature.

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.

Pre-reproductive isolation as a consequence of allopatric differentiation between populations of Drosophila melanogaster

While pre-reproductive isolations are more and more frequently described between closely related species or within species, very little is known about their conditions of emergence. In Brazzaville, two populations (Kronenbourg and Loua) of Drosophila melanogaster show a premating isolation. Two hypotheses were proposed to explain such a situation: a local sympatric differentiation or an allopatric divergence followed by a secondary contact. A microsatellite analysis, using markers on all chromosomes, strongly suggests that the Kronenbourg population has a European origin. Therefore, the allopatric divergence between Kronenbourg and Loua populations is probably responsible for the sexual isolation observed today in sympatry, after a recent introduction of a European propagule in Brazzaville.

Quantitative trait analysis and geographic variability of natural populations of Zaprionus indianus, a recent invader in Brazil

Five natural samples of a recent South America invader, the drosophilid Zaprionus indianus, were investigated with the isofemale line technique. These samples were compared to five African mainland populations, investigated with the same method. The results were also compared to data obtained on mass cultures of other populations from Africa and India. Three quantitative traits were measured on both sexes, wing and thorax length and sternopleural bristle number. We did not find any latitudinal trend among the American samples, while a significant increase in body size with latitude was observed in the Indian and, to a lesser degree, in the African populations. American populations were also characterized by their bigger size. Genetic variability, estimated by the intraclass correlation among isofemale lines, was similar in American and African populations. The intraline, nongenetic variability was significantly less in the American samples, suggesting a better developmental stability, the origin of which is unclear. A positive relationship was evident between intraline variability of size traits and the wing-thorax length correlation. Altogether, our data suggest that the colonizing propagule introduced to Brazil had a fairly large size, preventing any bottleneck effect being detected. The big body size of American flies suggests that they came from a high-latitude African country. The lack of a latitudinal dine in America seems to be related to the short time elapsed since introduction. The very rapid spread of Z. indianus all over South America suggests that it might rapidly invade North America.

Abundance, distribution and dynamics of retrotransposable elements and transposons: similarities and differences

Retrotransposable elements and transposons are generally both found in most eukaryotes. These two classes of elements are usually distinguished on the basis of their differing mechanisms of transposition. However, their respective frequencies, their intragenomic dynamics and distributions, and the frequencies of their horizontal transfer from one species to another can also differ. The main objective of this review is to compare these two types of elements from a new perspective, using data provided by genome sequencing projects and relating this to the theoretical and observed dynamics. It is shown that the traditional division into two classes, based on the transposition mechanisms, becomes less obvious when other factors are taken into consideration. A great diversity in distribution and dynamics within each class is observed. In contrast, the impact on and the interactions with the genome can show striking similarities between families of the two classes.

Classification and nomenclature of retrotransposable elements

The classification and nomenclature of retrotransposable elements is reviewed. A comparison is made between the initial classification summarized in Capy et al. (1997b), and the more recent proposal based on the classification of the viruses (Hull, 2001). Several problems, mainly relating to the position of elements belonging to the DIRS-like or Bel-like groups, are discussed. The first classification is now out of date, and must be revisited to take account of the discovery of new elements, however the second cannot be extended to the DNA elements. There is therefore, clear evidence of the need to adopt a general and a common classification.

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.

Isofemale lines in Drosophila: an empirical approach to quantitative trait analysis in natural populations

Founding isofemale lines from wild collected females is a basic tool for investigating the genetic architecture of Drosophila natural populations. The method permits the analysis of quantitative traits under laboratory conditions, with a much broader scope than the mere evidence of a significant genetic heterogeneity among lines. Genetic variability is generally demonstrated by a significant coefficient of intraclass correlation, but several experimental precautions are needed and explained here. The relationship between classical (additive) heritability and intraclass correlation is not straightforward, presumably because the genetic bottlenecks due to the initiation of the lines unravel a significant, nonadditive genetic variance due to dominance and epistatic effects. It is thus suggested to consider intraclass correlation as a specific genetic parameter that enables comparisons between different traits, different populations or different environments. The use of isofemale lines is, however, not restricted to the calculation of an intraclass correlation. It can be used to estimate genetic correlations among traits or environments. The method is also convenient for the analysis of phenotypic plasticity in relation to an environmental gradient. A precise description of the response curves (the reaction norms) is possible, distinguishing trait parameters and plasticity parameters. A fairly general conclusion is that, for a given trait, mean value and plasticity are genetically independent. It is also possible to analyze traits, which, like sexual dimorphism, must be measured on different individuals, and even to demonstrate their genetic variability. In many cases, further empirical and theoretical analyses are possible and needed. It is argued that, in the future, isofemale lines will have an increasing significance among the various techniques appropriate to the analysis of quantitative evolutionary genetics in a diversity of species.

Cchobo, a hobo-related sequence in Ceratitis capitata

A hobo-related sequence, Cchobo, with high similarity to the Drosophila melanogaster HFL1 and hobo108 elements was isolated from the medfly. Thirteen PCR-derived clones, which share 97.9-100% DNA identity, were sequenced, seven of which do not show frame-shift or stop codon mutations in their conceptual translations. The consensus sequence has 99.7% DNA identity with the D. melanogaster hobo element HFLI. In a phylogenetic analysis with other hobo-related elements, Cchobo clusters with the HFL1 and hobo108 elements from D. melanogaster and hobo-related elements from D. simulans, D. mauritiana and Mamestra brassicae. These elements may have undergone horizontal transfer in the recent past. The genomic distribution of Cchobo was studied by FISH to mitotic and polytene chromosomes, which revealed that Cchobo is distributed within both the heterochromatin and euchromatin. Intra- and interstrain polymorphisms were detected both at euchromatic and heterochromatic sites. These findings suggest that active copies of the element may be present in the medfly genome.

Comparative life histories and ecophysiology of Drosophila melanogaster and D. simulans

Numerous laboratory investigations have compared Drosophila melanogaster and D. simulans for various life history traits and fitness related ecophysiological parameters. From presently available information, it is however difficult to get a general comparative pattern describing the divergence of their ecological niches and understanding their demographic success. Two environmental factors seem however to have played a major role: temperature and alcoholic resources. From an ecophysiological approach, D. simulans may be described as generally more sensitive to stresses; other results point to this species as more cold adapted than its sibling; in some cases, however, D. simulans may appear as better adapted to a warm environment. When investigated, ecophysiological traits show a lesser geographic variability in D. simulans than in D. melanogaster. Presently available information does not explain the ecological prevalence of D. simulans in many places with a mild temperate or subtropical climate. This is presumably due to the fact that most comparisons have been done at a single, standard temperature of 25 degrees C. Comparative studies should be undertaken, spanning the thermal ranges of the two species, and the phenotypic plasticity of ecophysiological traits should now be considered.

Comparative analysis of morphological traits among Drosophila melanogaster and D. simulans: genetic variability, clines and phenotypic plasticity

The two sibling cosmopolitan species, Drosophila melanogaster and D. simulans, are able to proliferate under very different climatic conditions. This has resulted in local adaptations, which are often arranged in latitudinal clines. Such clines are documented for body weight, wing and thorax length, sternopleural and abdominal bristle number, ovariole number and thoracic pigmentation. The overall magnitude of geographical differentiation is, however, much less in D. simulans than in D. melanogaster, and latitudinal clines are less pronounced. The fact that natural populations live under different climates raises the problem of interaction between temperature and phenotype. The reaction norms of morphometrical traits have been investigated as a function of growth temperature. The shapes of the response curves vary according to the investigated trait. They are generally curvilinear and can be described by calculating characteristic values after polynomial adjustments. For a given trait, the reaction norms of the two species are similar in their shape, although some significant differences may be observed. Within each species, significant differences are also observed between geographic populations: reaction norms are not parallel and the divergence is better marked when more distant populations (e.g., temperate and tropical) are compared. It thus appears that besides mean trait value, phenotypic plasticity is also a target of natural selection. A specific analysis of wing shape variation according to growth temperature was also undertaken. Reaction norms with different shapes may be observed in various parts of the wing: the major effect is found between the basis and the tip of the wing, but in a similar way in the two species. By contrast, some ratios, called wing indices by taxonomists, may exhibit completely different reaction norms in the two species. For a single developmental temperature (25 degrees C) the phenotypic variability of morphometrical traits is generally similar in the two species, and also the genetic variability, estimated by the intraclass correlation. A difference exists, however, for the ovariole number which is less variable in D. simulans. Variance parameters may vary according to growth temperature, and a detailed analysis was made on wing dimensions. An increase of environmental variability at extreme, heat or cold temperatures, has been found in both species. Opposite trends were, however, observed for the genetic variability: a maximum heritability in D. simulans at middle temperatures, corresponding to a minimum heritability in D. melanogaster. Whether such a difference exists for other traits and in other populations deserves further investigations. In conclusion, morphometrical analyses reveal a large amount of significant differences which may be related to speciation and to the divergence of ecological niches. Within each species, numerous geographic variations are also observed which, in most cases, reflect some kinds of climatic adaptation.