BACKGROUND
the zaprionus genus  is composed by two subgenera  which seem to be originated in asia about  <dig> mya  <cit> .
the anaprionus subgenus has first diversified in the oriental biogeographic region, being relatively less studied.
however, the zaprionus subgenus ancestor originated during the middle to early miocene in oriental regions before diversifying in tropical africa, sharing space and time with the best studied drosophilid species, the melanogaster subgroup  <cit> .
due to its evolutionary history, ecological and morphological diversity, the zaprionus genus seems to be a good model for comparative studies with the melanogaster subgroup.
although the phylogenetic relationships within the zaprionus genus had been recently proposed  <cit> , its taxonomic positioning in the drosophilidae family remains a matter for discussion.
currently, most reports agree that zaprionus belongs to the drosophila genus, and that it would be more closely related to drosophila than to sophophora subgenera  <cit> , to which the species of the melanogaster subgroup belong .
hence, zaprionus and melanogaster species would share a last common ancestor at least as old as the divergence of the sophophora and drosophila subgenera, i.e., between  <dig> and  <dig> mya  <cit> .
while the melanogaster species subgroup - mainly the two sibling species d. melanogaster and d. simulans - presents the most well known dipteran mobilomes, knowledge about the transposable elements  content of the zaprionus species is scant.
only seven te families have been identified in all zaprionus genus.
the elements  <dig> of z. tuberculatus,  <dig> of z. ornatus, bari- <dig> of z. tuberculatus and, mariner of  <dig> zaprionus species could only be detected by southern blot hybridization using probes derived from d. melanogaster  <cit> .
furthermore, partial sequences were produced for gypsy of z. indianus, copia of z. tuberculatus, mariner of z. tuberculatus and z. verruca, and rover of z. indianus  <cit> .
interestingly, evolutionary analyses demonstrated that four of these seven tes of zaprionus are highly related to those the melanogaster subgroup, sharing their common ancestor more recently than the species ones.
gypsy, copia, mariner and rover were therefore proposed to have been horizontally transferred between species of the zaprionus subgenus and the melanogaster subgroup of the genus drosophila  <cit> .
horizontal transfer  of tes has been inferred when three criteria are met:  high sequence conservation between te sequences of distantly related species;  incongruence between te and host phylogenies and,  discontinuous distribution of a te across a species group.
additionally, requirements like geographic, temporal and ecological overlapping and, exclusion of alternative hypothesis as ancestral polymorphism promoting differential subfamily fixation and high selective constraints over the sequences need to be considered  <cit> .
in a previous study, we have performed a search of d. melanogaster tes in the american continent invasive species z. indianus by dot blot methodology, aiming to broaden our understanding of te occurrence in species of the zaprionus genus  <cit> .
all the  <dig> tes analyzed were present in the z. indianus genome, but only five had strong hybridization signals with the probes that were used: roo/b <dig>  doc, mdg- <dig>  micropia and gypsy.
among them, we selected gypsy and micropia retroelements to expand our evolutionary analyses in the zaprionus genus.
this choice was due to the evidence of the close relationship of these z. indianus elements with others already described in the literature, inferred by the strong hybridization signals using heterologous probes, the easy amplification with the primers available in the literature and the proposition of a gypsy ht of between z. indianus and d. simulans  <cit> .
gypsy and micropia are classified as ty3/gypsy elements because of their gene arrangement in the pol coding region - reverse transcriptase/rnaseh/integrase  <cit> .
ty3/gypsy retroelements, which are similar to vertebrate retroviruses in both sequence and genomic structure, are broadly distributed among eukaryotes  <cit> .
they differ from vertebrate retroviruses mainly because they do not have an env gene necessary to complete the infectious cycle  <cit> .
gypsy is an exception to the typical ty3/gypsy structure because an active env is present, potentially capable of maintaining infectivity  <cit> .
hence, gypsy was classified as an errantivirus by the international committee on taxonomy of viruses   <cit> .
the micropia retroelement had already been studied in d. melanogaster, and several species of the repleta and the cardini groups of the subgenus drosophila  <cit> .
evolutionary analyses showed the existence of at least two micropia families that have ~30% divergence in dna sequence  <cit> .
the first family, present in the repleta group, comprises two subfamilies than were differentiated by nucleotide divergence and antisense rna expression.
the repleta species group seems to have transmitted one of its elements to the d. cardini group through ht  <cit> .
the second family is represented by the d. melanogaster element that also transcribes antisense rnas responsible for control of its expression  <cit> .
the data suggest that the micropia retroelement is an ancient component of drosophila genomes  <cit> .
in contrast to micropia, the gypsy retroelement has been analyzed in the  <dig> drosophila genomes  <cit> , and individual sequences from several drosophilidae species are available, such as the virilis and repleta groups of subgenus drosophila; the willistoni and melanogaster groups of the subgenus sophophora; and the scaptodrosophila latifasciaeformis and the z. indianus species  <cit> .
two gypsy families have been identified, one exclusive to d. willistoni and the other widely distributed among drosophilidae species  <cit> .
the latter is divided into  <dig> subfamilies, seven of which have been implicated in ht events, demonstrating that inheritance of the gypsy retroelement is complex  <cit> .
ht events involving gypsy could have been facilitated by the presence of env, assumed to confer autonomous infectivity on the retroelement, thereby allowing it to cross species barriers  <cit> .
specifically in the melanogaster species subgroup, two gypsy subfamilies have been identified, and seven hts have been proposed among melanogaster species subgroups, or between one species of that subgroup and s. latifasciaeformis and z. indianus  <cit> .
these inter-genera hts are assumed to have occurred because synonymous sites of gypsy env gene are three to five times more highly conserved than those of host genes in species that diverged longer than  <dig> mya  <cit> .
we have now surveyed the gypsy and micropia retroelements in seven species of the zaprionus genus, comprising one species of the anaprionus subgenus and six of the zaprionus subgenus.
our objective was to investigate the distribution and evolutionary history of these two elements in this scarcely studied genus by comparing their sequences with those of the  <dig> drosophila genomes and other sequences available in the nucleotide database.
the results show that the zaprionus gypsy and micropia retroelements are present and transcriptionally active only in the subgenus zaprionus, to which they were introduced by more than one horizontal transfer event from different donor species.
RESULTS
distribution of gypsy and micropia in the genus zaprionus
in an attempt to search for the presence of gypsy and micropia retroelements in seven species of the zaprionus genus , pcr reactions using primers that amplify gypsy env and micropia rnaseh sequences were carried out .
this revealed the presence of these elements in the six species of the zaprionus subgenus, but not in z. multistriatus .
three clones of each gypsy and micropia pcr fragments were sequenced and used in the evolutionary analyses.
additionally, rt-pcr reactions and southern blot hybridization showed that all zaprionus subgenus species analyzed harbor a few transcriptionally active insertions of gypsy and micropia elements, varying from one  to seven  insertions for gypsy, and from two  to seven  for micropia .
nd: not determined.
evolutionary analyses of the gypsy retroelement
to infer the evolutionary history of gypsy retroelement in the zaprionus subgenus, we carried out genomic searches for zaprionus homologous sequences in the  <dig> drosophila genomes and the genbank database.
the in silico search retrieved genomic sequences only from the melanogaster subgroup, indicating that the other seven drosophila genomes  do not harbor gypsy sequences closely related to those of the subgenus zaprionus.
eight sequences were identified in d. melanogaster, nine in d. simulans,  <dig> in d. sechellia, six in d. yakuba, and  <dig> in d. erecta .
a blast-based analysis of the structure sequence indicated that only d. melanogaster genome presents putatively full-length insertions .
inferring phylogenetic relationships and divergence
the phylogenetic relationships inferred using the maximum parsimony , maximum likelihood  and neighbor-joining  methods of phylogenetic reconstruction produced similar patterns with identical main branches.
the analyses were performed with three sequences of each zaprionus species and  <dig> sequences of species from the melanogaster subgroup.
the s. latifasciaeformis and d. willistoni gypsy sequences available in the genbank database were included in these reconstructions because of their respectively closer and more distant evolutionary relationships to those of the melanogaster subgroup  <cit> .
figure  <dig> shows the tree reconstructed by ml analysis.
the sequences of the genus zaprionus, the melanogaster species subgroup and s. latifasciaeformis were grouped together in three clades that were not congruent with the species phylogeny.
clade  <dig>  contains sequences of d. melanogaster, d. erecta, d. simulans and only one sequence of the genus zaprionus, z. tuberculatus <dig>  clade  <dig>  includes z. tuberculatus, z. camerounensis, s. latifasciaeformis, d. melanogaster, d. erecta and all sequences of d. yakuba, except d. yakubaa <dig>  which branches outside the  <dig> clades.
the reported clustering of d. erecta and s. latifasciaeformis has been explained as an ht event  <cit> .
finally, clade  <dig> comprises the gypsy sequences of z. davidi, z. indianus and the two cryptic species, z. gabonicus and z. africanus, as well as sequences of d. melanogaster, d. simulans and all sequences of d. sechellia, with no bootstrap support.
additionally, a tree was reconstructed using the zaprionus sequences, the genomic sequences of melanogaster subgroup, and gypsy sequences of  <dig> other drosophila species .
a similar topology of the figure  <dig> was obtained, with all zaprionus sequences clustering together with melanogaster elements in the same clade distribution, which reinforces the close relationship between zaprionus and melanogaster gypsy sequences .
distance analysis shows that the sequences within the three clades have relatively low divergence levels, with mean values of  <dig> ,  <dig>  and  <dig>  within clades  <dig>   <dig> and  <dig>  respectively .
the distances of the zaprionus vs melanogaster species within each group are similar to the mean values just cited.
the mean distances between the clades show that clade  <dig> is closer to clade  <dig> than to clade  <dig>  with mean values of  <dig>  ,  <dig>   and  <dig>  .
the drosophila gypsy sequences, which vary up to 20% in nucleotide composition, were previously classified into  <dig> subfamilies  <cit> .
among them, the sequences of z. indianus and the melanogaster species subgroup were included in two different subfamilies.
the addition of  <dig> sequences belonging to six zaprionus species allowed us to confirm this proposition.
the species distribution and divergence of the clade  <dig> from  <dig> and  <dig>  indicate that it corresponds to one of those subfamilies previously proposed, hereafter called subfamily  <dig> .
on the other hand, a divergence of  <dig>  between clades  <dig> and  <dig> led us to classify them as variants  <dig> and  <dig>  of another subfamily  that may correspond with the second subfamily already described due to its species composition.
note: the mcl pairwise distance values between the zaprionus subgenus and melanogaster species subgroup were summarized according the gypsy subfamilies and variants and, species groups for micropia.
s <dig> - subfamily 1; s <dig> - subfamily 2; v <dig> - variant 1; v <dig> - variant  <dig>  n: number of pairwise comparisons.
testing the ht hypothesis
three hypotheses could explain the high levels of conservation among gypsy sequences within each clade  and the phylogenetic inconsistencies compared to the species tree:  high selective constraints conserving the sequences;  ancestral polymorphism promoting differential subfamily fixation in the species; and  the occurrence of ht.
to test whether  can best explain the estimated conservation, the ratio of non-synonymous divergence per synonymous divergence  within the clades were compared to those of the gpdh nuclear gene, offering a comparative measurement of the selective constraints.
although active elements might suffer selective pressure to maintain their coding sequences, gpdh would be expected to be even more restricted because it plays an essential role in glycerophospholipid metabolism in drosophila.
if it is assumed that synonymous substitutions are under almost strictly neutral evolution, dn:ds =  <dig>  dn:ds < <dig> and dn:ds > <dig> will represent neutral evolution, purifying selection and positive darwinian selection, respectively.
the analyses indicate that the mean gypsy dn:ds ratios ranged from  <dig>  to  <dig> , while the mean gpdh values did not exceed  <dig>  .
these ratios suggest that, although dn:ds < <dig>  purifying selection over the gypsy env sequences is relaxed  compared to the host gene.
also, a z-test for the gypsy env region indicates neutrality for most of the pairwise comparisons  within the clades .
hence, high selective constraints do not explain the incongruities observed intra-clades.
on the other hand, the neutrality hypothesis was refuted for all pairwise comparisons between the clades, indicating that purifying selection could plays a role in the conservation of the env coding sequence among the gypsy subfamilies.
in that case, natural selection could act at the genomic level, conserving the te genes, since elements that transpose most efficiently, or at the highest rate are most likely to survive and propagate  <cit> .
note: pairwise ds and dn distances in additional files  <dig> and  <dig>
in the absence of intra-clade selective constraints, the comparisons between the ds values of the te and a host gene permit evaluation of the hypothesis of ancestral polymorphism within subfamilies.
if the gypsy retroelement is evolving vertically and the incongruence in the phylogeny are due to differential subfamily distribution in the analyzed species, the ds pairwise comparisons between the te and the host gene would be expected to be equivalent.
but, if the ds of te is lower than the ds of host gene under similar or higher levels of selective constraints and in the absence of strong codon usage bias, ht events can be inferred.
using this approach, the ds of gypsy and gpdh gene, as well as the codon bias index  were compared for one clone representing each species of zaprionus and the best blastn subject of the genomic sequences in each clade.
codon usage bias is a factor that may be responsible for low ds values for te sequences, since there is a negative correlation between codon usage bias and ds values  <cit> .
the codon bias index  indicates that the gypsy and gpdh sequences did not suffer high codon bias, although the mean value of  <dig>  for gpdh suggests that the small difference in magnitude of ds between gpdh and gypsy ought to be taken cautiously in inferring ht, since the variation in the gpdh ds distance could reflect codon usage.
hence, comparisons between gypsy and gpdh ds were used to infer ht only if the ds values of gypsy were at least twice as low as those of gpdh.
the gypsy ds values between zaprionus and melanogaster species ranged from  <dig>  to  <dig> times  lower than those of gpdh .
on the other hand, comparisons of ds between the zaprionus species within each clade gave gypsy ds values approximately equal or higher than gpdh ds, excepting for z. davidi <dig> vs z. indianus <dig> and z. africanus <dig> which presented ds of gypsy slightly lower than the host gene ones.
these results may indicate that the gypsy retroelement evolved mainly by vertical transmission within the zaprionus subgenus, but there may have been ht events between the zaprionus and the melanogaster species what could reject the hypothesis of ancestral polymorphism within the clades.
phylogenetic relationships of the micropia retroelement
in order to evaluate the evolutionary history of micropia, the same approaches of phylogeny and divergence estimations used for gypsy were applied.
sequences homologous to the zaprionus species were also identified only in species of the melanogaster subgroup, except for d. erecta, the unique species of the melanogaster subgroup for which no significant matches were obtained.
d. sechellia had the most hits , followed by d. simulans , d. yakuba  and d. melanogaster  .
structure analysis indicates that the d. melanogaster and d. simulans genomes harbor putatively full-length insertions .
inferring phylogenetic relationships and divergence
the mp, nj and ml reconstructions of the micropia sequences from the zaprionus subgenus, the genomic sequences from the melanogaster subgroup and the genbank sequences of the repleta group had similar topologies.
as for gypsy, only ml reconstruction is given .
the tree exhibits several incongruences according to species phylogeny, such as grouping species in the zaprionus genus together with those of the melanogaster subgroup.
the sequences of the micropia family of the repleta group were clustered outside a well-supported branch that harbors all sequences of the zaprionus subgenus and the melanogaster species subgroup, constituted by two clades.
clade  <dig>  with low support , grouped the d. sechelliaa <dig> sequence with three internal clades, corresponding to the z. camerounensis sequences , d. yakuba  and the sequences of z. davidi, z. indianus, z. africanus and z. gabonicus .
this last clade did not give a clear internal species-specific clustering.
clade  <dig> grouped sequences of the melanogaster subgroup  together with z. tuberculatus in 76% of the replications.
despite the clustering of all melanogaster sequences in clade  <dig>  the support was not high .
the mean divergence within clade  <dig> was much lower  than in clade  <dig>  or between them  .
analyzing the species groups separately, the mean divergence of the zaprionus sequences was  <dig> ; it was  <dig>  within the melanogaster species subgroup, and  <dig>  between the zaprionus and the melanogaster sequences .
the phylogenetic tree and distance estimates <20% indicate that micropia in the zaprionus and melanogaster species belong to the same subfamily.
testing the ht hypothesis
evaluation of the selective constraints showed that the mean dn:ds values of the micropia sequences were 14- <dig> times higher than those of gpdh .
between the zaprionus and melanogaster sequences, for example, the values were  <dig>  and  <dig>  for micropia and gpdh, respectively.
in addition, the z-test failed to reject the neutrality hypothesis in ~80% of the pairwise comparisons .
these results refute the hypothesis that the high similarity between the micropia sequences of zaprionus and melanogaster is due to high selective constraints conserving their nucleotide sequences.
to evaluate the ancestral polymorphism hypothesis for micropia, the ds values of micropia rnaseh and gpdh sequences were compared between the sequences of zaprionus and melanogaster species in each clade .
all the ds distances for the micropia were lower than those for the gpdh.
for instance, the ds between z. tuberculatus and the melanogaster subgroup sequences in clade  <dig> were from  <dig>  to  <dig>  times lower than that those for the gpdh.
those for the d. yakubaa <dig> sequence were from  <dig>  to  <dig>  times lower.
in addition, ds was zero in some comparisons, e.g.
those among the z. camerounensis, z. davidi, z. africanus and z. gabonicus sequences of clade  <dig>  the high conservation in the synonymous sites of the rnaseh sequences favors the hypothesis that their low divergence is due to recent acquisitions of micropia in the zaprionus and melanogaster genomes, and that ht has shaped its evolution, as already proposed for the gypsy retroelement.
network trees corroborate the inference of ht
network reconstructions were used as an additional approach to test the hypothesis that ht has shaped the evolution of gypsy and micropia in zaprionus and melanogaster species.
since te sequences can be considered as populations of sequences that share a common ancestor, this type of phylogeny lets one propose possible donors and receptors species, and to reconstruct the te dispersion routes by the proposition of the ancestral nodes.
the gypsy network indicated three clusters of sequences , which correspond to the s <dig> subfamily and the two variants of the s <dig> subfamily .
moreover, the s <dig> subfamily was closer to s2v <dig> than to s2v <dig>  which indicates that s2v <dig> diversified from s2v <dig>  the presence of median vectors - that correspond to a hypothetical  or unsampled sequences - connecting the zaprionus and melanogaster sequences suggests that the ht events occurred in the ancestors of the sequences sampled, indicating that the lateral transmissions were not recent.
the network tree allowed us to see a possible ht not clearly showed by the traditional phylogeny, but evident from the ds comparisons.
the z. davidi <dig> sequence is the ancestral node of the z. africanus <dig> sequence , suggesting a single ht event from z. davidi to z. africanus.
the network reconstruction also let the d. yakubaa <dig> sequence be included as a s2v <dig> member.
network reconstruction for micropia clustered the sequences of d. yakuba, d. melanogaster, z. camerounensis, z. davidi, z. indianus and z. africanus with central median vectors, suggesting that high sequence conservation was not due to hts involving the species analyzed, but probably one or more donor species not sampled .
the d. melanogastera <dig> sequence connects the median vectors with a large cluster that harbors d. sechellia, d. simulans and d. melanogaster sequences, which in turn links a median branch to the z. tuberculatus sequences.
this arrangement indicates at least two more recent transfers, one from d. melanogaster to d. simulans and/or d. sechellia or even to the ancestor of those sister species , and a second from any of these species to z. tuberculatus.
in addition, the branch of the d. sechelliaa <dig> sequence with a median vector that links the z. camerounensis, z. davidi, z. indianus and z. africanus sequences suggest an additional ht involving the ancestor of d. sechellia.
however, the absence of geographic sharing between d. sechellia and zaprionus species means that this proposal should be viewed cautiously.
an alternative explanation could be a high rate of evolution of the d. sechelliaa <dig> sequence, followed by random convergent evolution.
ht inference and estimation of the divergence time between te sequences
to reinforce the ht hypothesis, the divergence times of both retroelements were calculated using the molecular clock .
the divergence estimates for gypsy sequences showed that the hts may have started around  <dig>  mya .
regarding the s <dig> subfamily, two hts can be proposed.
the first may have occurred from a d. melanogaster/d.
simulans ancestor to z. tuberculatus about  <dig> - <dig>  mya, prior to the time of d. melanogaster/d.
simulans splitting .
this event was followed by a transfer between d. melanogaster and d. erecta , species that share a common ancestor at 8- <dig> mya  <cit> .
although an alternative hypothesis of ht between d. simulans and d. erecta can be formulated on the basis of the phylogeny branching, the greater difference in ds comparisons between the pair d. erecta/d.
melanogaster than the pair d. erecta/d.
simulans  reinforces the first hypothesis.
symbols: *: species ancestor.
#: time of divergence calculated using ds pairwise comparisons between sequences connected by the network phylogeny.
the gypsy s2v <dig> sequences seem to have been involved at least in four ht events.
we suggest that d. erecta was the donor of s2v <dig> in a more ancient transfer to z. tuberculatus , followed by introduction into the z. camerounensis/z.
davidi ancestor , and lately into d. yakuba  and d. melanogaster .
the transfer to z. camerounensis/z.
davidi probably occurred in their ancestor, because these species diverged around  <dig>  mya  <cit> , but the hypothesis of two independent hts cannot be disregarded.
the donor status of d. erecta was based on their basal positioning in the network, the repeated clustering to all the other species of this variant in the phylogeny, the divergence times between d. erecta and d. melanogaster/d.
yakuba .
the s2v <dig> variant is probably involved in one ht event that happened around  <dig> - <dig>  mya.
it could have occurred between the ancestor of the modern sequences of d. melanogaster/d.
simulans/d.
sechellia and the indianus/davidi complexes.
alternatively, the ht could have happened in each complex separately, since the davidi complex diverged about  <dig>  mya  <cit> .
the ht between the d. melanogaster/d.
simulans/d.
sechellia ancestor and the indianus complex ancestor is a broader scenario than the ht proposed between z. indianus and d. simulans  <cit> , exemplifying the importance of analyzing as many species as possible before inferring horizontal transfers  <cit> .
finally, the time of divergence also corroborates the ht from z. davidi to z. africanus s2v <dig> sequences, as seen in the network tree, which probably occurred ~ <dig>  mya, subsequent to the divergence of these species .
comparing the divergence times between the micropia retroelement of the zaprionus and melanogaster sequences with those of gypsy, the hts of micropia might have started more recently than those of gypsy .
the invariance of the synonymous sites between several pairwise comparisons  makes it difficult to date these events, but the highest ds levels indicate that, of the species we have analyzed, d. yakuba and the indianus complex species were the first to be invaded by micropia.
so, the evolutionary analyses lead to the inference that micropia diverged after the splitting of their species hosts, suggesting that at least an unknown species might have transferred micropia sequences independently in five hts to d. melanogaster, d. yakuba , z. camerounensis, z. davidi and the indianus complex .
species of neither the repleta nor the cardini species group are good candidates, since they harbor micropia retrotransposons of a different family  <cit> .
hence, the probable donor might have been a non-studied species.
based mainly on the network tree, two other ht events can be proposed.
d. melanogaster may have been the donor of micropia to the d. simulans and d. sechellia, or even to their ancestor.
however, introgression cannot be excluded in this case.
in turn, one of these three species might have been involved in the transmission of micropia to z. tuberculatus, about  <dig> - <dig>  mya.
discussion
this survey of the gypsy and micropia in seven species of zaprionus indicates that they are widely distributed and transcriptionally active only in species of subgenus zaprionus.
additionally, the evolutionary analyses demonstrate that these two retroelements are closely related to those of the melanogaster species subgroup, and that their histories might have been repeatedly marked by events of inter- and intra-subgenus hts.
several pieces of evidence have shown that the phenomenon of ht is frequent in eukaryotes  <cit> .
it is important to note that the hts proposed here were assumed only when the species involved in the event shared geographic, temporal and ecological environments, and when three separate pieces of evidence suggested their occurrence:  lower ds values between the te sequences than the gpdh sequences,  incongruence between host and te phylogenies, and  last common ancestor of te more recent than that of the species.
the studies of micropia and gypsy in subgenus zaprionus species highlight the ht events as an alternative evolutionary mechanism of the retroelements evolution.
among the hts inferred for gypsy, four events seemed to involve the zaprionus and melanogaster species, three involved only the melanogaster group species, and only one the zaprionus species.
for the micropia retrotransposon, seven hts were suggested: five between unknown donor and the zaprionus or the melanogaster species, one among the melanogaster species, and one between the zaprionus and the melanogaster species.
the  <dig> hts we have proposed probably occurred during two different waves of invasions.
first, the gypsy retroelement of the melanogaster subgroup had been introduced in the zaprionus subgenus ~ <dig> mya.
after that , a second ht wave involving the micropia retroelement may have introduced it into both the zaprionus subgenus and melanogaster subgroup from an unknown donor.
after the introduction, the tes could have been prone to transposition and re-introduction in other related species.
for example, z. davidi has donated their gypsy sequence to z. africanus after had received it from an ancestor of d. melanogaster, d. simulans and d. sechellia.
according to the te evolutionary cycle, the te history could start from ht events, followed by the initial transpositional burst in the new host, and then the accumulation of defective copies along with host-directed epigenetic silencing  <cit> .
the last step could lead to the loss of mobility and, finally, to molecular erosion by random mutations.
the high frequency of hts, combined with the transcriptional activity and low insertion numbers of gypsy and micropia retroelements in the zaprionus species subgenus indicates that these elements remain at the first or second steps of the evolutionary life history of the tes, i.e.
the invasion or the genomic spread period.
a question that arises is whether these elements have been transferred in similar ways.
hts of gypsy are easier to explain, since the active gypsy env can enable its dispersion and expansion of the genomic territory, but this mechanism does not to apply to micropia, which lacks env.
nevertheless, retroelements could have been transferred concomitantly with retrovirus infections, even gypsy; or parasitic infestations, as for example by intracellular wolbachia bacteria or mites  <cit> .
moreover, our data allowed us to propose there were ht of gypsy and micropia among sister species such as d. melanogaster, d. simulans and d. sechellia, although introgression could not be ruled out.
another concern is the evolutionary process of fixation of the tes after the horizontal introduction.
these elements could share similar mechanisms of regulation, which would permit them escaping from selective forces and allow fixation after invasion.
control of drosophila retroelement mobilization depends on the transcription rate, which is directly related to the presence of specific regulatory proteins  <cit>  and the rate of rna degradation, mediated by rnai systems  <cit> .
however, genetic information about the species of the genus zaprionus is scarce.
further, studies of the zaprionus genes involved in controlling mobilization could lead to hypotheses about the mechanisms that allow highly related sequences to be shared in the zaprionus and melanogaster species genomes.
CONCLUSIONS
our survey suggests that zaprionus species have experienced waves of retroelement invasions, particularly during the last  <dig> million years.
species of the melanogaster subgroup might have donated their gypsy sequences to species of the zaprionus subgenus, and these two species groups could have received the micropia retroelement from one or more unknown donor species.
after the initial introduction, these species could share their elements.
since the zaprionus genus and the melanogaster subgroup seem to share the same age of origin and diversification in tropical africa, as well as ecological features, our data suggest that they passed through a permissive period of transposable element invasion during the diversification period.

 the zaprionus genus is composed by two subgenera which seem to be originated in asia about mya. the anaprionus subgenus has first diversified in the oriental biogeographic region, being relatively less studied. however, the zaprionus subgenus ancestor originated during the middle to early miocene in oriental regions before diversifying in tropical africa, sharing space and time with the best studied drosophilid species, the melanogaster subgroup. due to its evolutionary history, ecological and morphological diversity, the zaprionus genus seems to be a good model for comparative studies with the melanogaster subgroup. although the phylogenetic relationships within the zaprionus genus had been recently proposed, its taxonomic positioning in the drosophilidae family remains a matter for discussion. currently, most reports agree that zaprionus belongs to the drosophila genus, and that it would be more closely related to drosophila than to sophophora subgenera, to which the species of the melanogaster subgroup belong. hence, zaprionus and melanogaster species would share a last common ancestor at least as old as the divergence of the sophophora and drosophila subgenera, i.e., between and mya. while the melanogaster species subgroup - mainly the two sibling species d. melanogaster and d. simulans - presents the most well known dipteran mobilomes, knowledge about the transposable elements content of the zaprionus species is scant. only seven te families have been identified in all zaprionus genus. the elements of z. tuberculatus, of z. ornatus, bari- of z. tuberculatus and, mariner of zaprionus species could only be detected by southern blot hybridization using probes derived from d. melanogaster. furthermore, partial sequences were produced for gypsy of z. indianus, copia of z. tuberculatus, mariner of z. tuberculatus and z. verruca, and rover of z. indianus. interestingly, evolutionary analyses demonstrated that four of these seven tes of zaprionus are highly related to those the melanogaster subgroup, sharing their common ancestor more recently than the species ones. gypsy, copia, mariner and rover were therefore proposed to have been horizontally transferred between species of the zaprionus subgenus and the melanogaster subgroup of the genus drosophila. horizontal transfer of tes has been inferred when three criteria are met: high sequence conservation between te sequences of distantly related species; incongruence between te and host phylogenies and, discontinuous distribution of a te across a species group. additionally, requirements like geographic, temporal and ecological overlapping and, exclusion of alternative hypothesis as ancestral polymorphism promoting differential subfamily fixation and high selective constraints over the sequences need to be considered. in a previous study, we have performed a search of d. melanogaster tes in the american continent invasive species z. indianus by dot blot methodology, aiming to broaden our understanding of te occurrence in species of the zaprionus genus. all the tes analyzed were present in the z. indianus genome, but only five had strong hybridization signals with the probes that were used: roo/b doc, mdg- micropia and gypsy. among them, we selected gypsy and micropia retroelements to expand our evolutionary analyses in the zaprionus genus. this choice was due to the evidence of the close relationship of these z. indianus elements with others already described in the literature, inferred by the strong hybridization signals using heterologous probes, the easy amplification with the primers available in the literature and the proposition of a gypsy ht of between z. indianus and d. simulans. gypsy and micropia are classified as ty3/gypsy elements because of their gene arrangement in the pol coding region - reverse transcriptase/rnaseh/integrase. ty3/gypsy retroelements, which are similar to vertebrate retroviruses in both sequence and genomic structure, are broadly distributed among eukaryotes. they differ from vertebrate retroviruses mainly because they do not have an env gene necessary to complete the infectious cycle. gypsy is an exception to the typical ty3/gypsy structure because an active env is present, potentially capable of maintaining infectivity. hence, gypsy was classified as an errantivirus by the international committee on taxonomy of viruses. the micropia retroelement had already been studied in d. melanogaster, and several species of the repleta and the cardini groups of the subgenus drosophila. evolutionary analyses showed the existence of at least two micropia families that have ~30% divergence in dna sequence. the first family, present in the repleta group, comprises two subfamilies than were differentiated by nucleotide divergence and antisense rna expression. the repleta species group seems to have transmitted one of its elements to the d. cardini group through ht. the second family is represented by the d. melanogaster element that also transcribes antisense rnas responsible for control of its expression. the data suggest that the micropia retroelement is an ancient component of drosophila genomes. in contrast to micropia, the gypsy retroelement has been analyzed in the drosophila genomes, and individual sequences from several drosophilidae species are available, such as the virilis and repleta groups of subgenus drosophila; the willistoni and melanogaster groups of the subgenus sophophora; and the scaptodrosophila latifasciaeformis and the z. indianus species. two gypsy families have been identified, one exclusive to d. willistoni and the other widely distributed among drosophilidae species. the latter is divided into subfamilies, seven of which have been implicated in ht events, demonstrating that inheritance of the gypsy retroelement is complex. ht events involving gypsy could have been facilitated by the presence of env, assumed to confer autonomous infectivity on the retroelement, thereby allowing it to cross species barriers. specifically in the melanogaster species subgroup, two gypsy subfamilies have been identified, and seven hts have been proposed among melanogaster species subgroups, or between one species of that subgroup and s. latifasciaeformis and z. indianus. these inter-genera hts are assumed to have occurred because synonymous sites of gypsy env gene are three to five times more highly conserved than those of host genes in species that diverged longer than mya. we have now surveyed the gypsy and micropia retroelements in seven species of the zaprionus genus, comprising one species of the anaprionus subgenus and six of the zaprionus subgenus. our objective was to investigate the distribution and evolutionary history of these two elements in this scarcely studied genus by comparing their sequences with those of the drosophila genomes and other sequences available in the nucleotide database. the results show that the zaprionus gypsy and micropia retroelements are present and transcriptionally active only in the subgenus zaprionus, to which they were introduced by more than one horizontal transfer event from different donor species.  distribution of gypsy and micropia in the genus zaprionus in an attempt to search for the presence of gypsy and micropia retroelements in seven species of the zaprionus genus, pcr reactions using primers that amplify gypsy env and micropia rnaseh sequences were carried out. this revealed the presence of these elements in the six species of the zaprionus subgenus, but not in z. multistriatus. three clones of each gypsy and micropia pcr fragments were sequenced and used in the evolutionary analyses. additionally, rt-pcr reactions and southern blot hybridization showed that all zaprionus subgenus species analyzed harbor a few transcriptionally active insertions of gypsy and micropia elements, varying from one to seven insertions for gypsy, and from two to seven for micropia. nd: not determined. evolutionary analyses of the gypsy retroelement to infer the evolutionary history of gypsy retroelement in the zaprionus subgenus, we carried out genomic searches for zaprionus homologous sequences in the drosophila genomes and the genbank database. the in silico search retrieved genomic sequences only from the melanogaster subgroup, indicating that the other seven drosophila genomes do not harbor gypsy sequences closely related to those of the subgenus zaprionus. eight sequences were identified in d. melanogaster, nine in d. simulans, in d. sechellia, six in d. yakuba, and in d. erecta. a blast-based analysis of the structure sequence indicated that only d. melanogaster genome presents putatively full-length insertions. inferring phylogenetic relationships and divergence the phylogenetic relationships inferred using the maximum parsimony, maximum likelihood and neighbor-joining methods of phylogenetic reconstruction produced similar patterns with identical main branches. the analyses were performed with three sequences of each zaprionus species and sequences of species from the melanogaster subgroup. the s. latifasciaeformis and d. willistoni gypsy sequences available in the genbank database were included in these reconstructions because of their respectively closer and more distant evolutionary relationships to those of the melanogaster subgroup. figure shows the tree reconstructed by ml analysis. the sequences of the genus zaprionus, the melanogaster species subgroup and s. latifasciaeformis were grouped together in three clades that were not congruent with the species phylogeny. clade contains sequences of d. melanogaster, d. erecta, d. simulans and only one sequence of the genus zaprionus, z. tuberculatus clade includes z. tuberculatus, z. camerounensis, s. latifasciaeformis, d. melanogaster, d. erecta and all sequences of d. yakuba, except d. yakubaa which branches outside the clades. the reported clustering of d. erecta and s. latifasciaeformis has been explained as an ht event. finally, clade comprises the gypsy sequences of z. davidi, z. indianus and the two cryptic species, z. gabonicus and z. africanus, as well as sequences of d. melanogaster, d. simulans and all sequences of d. sechellia, with no bootstrap support. additionally, a tree was reconstructed using the zaprionus sequences, the genomic sequences of melanogaster subgroup, and gypsy sequences of other drosophila species. a similar topology of the figure was obtained, with all zaprionus sequences clustering together with melanogaster elements in the same clade distribution, which reinforces the close relationship between zaprionus and melanogaster gypsy sequences. distance analysis shows that the sequences within the three clades have relatively low divergence levels, with mean values of, and within clades and respectively. the distances of the zaprionus vs melanogaster species within each group are similar to the mean values just cited. the mean distances between the clades show that clade is closer to clade than to clade with mean values of, and. the drosophila gypsy sequences, which vary up to 20% in nucleotide composition, were previously classified into subfamilies. among them, the sequences of z. indianus and the melanogaster species subgroup were included in two different subfamilies. the addition of sequences belonging to six zaprionus species allowed us to confirm this proposition. the species distribution and divergence of the clade from and indicate that it corresponds to one of those subfamilies previously proposed, hereafter called subfamily. on the other hand, a divergence of between clades and led us to classify them as variants and of another subfamily that may correspond with the second subfamily already described due to its species composition. note: the mcl pairwise distance values between the zaprionus subgenus and melanogaster species subgroup were summarized according the gypsy subfamilies and variants and, species groups for micropia. s - subfamily 1; s - subfamily 2; v - variant 1; v - variant n: number of pairwise comparisons. testing the ht hypothesis three hypotheses could explain the high levels of conservation among gypsy sequences within each clade and the phylogenetic inconsistencies compared to the species tree: high selective constraints conserving the sequences; ancestral polymorphism promoting differential subfamily fixation in the species; and the occurrence of ht. to test whether can best explain the estimated conservation, the ratio of non-synonymous divergence per synonymous divergence within the clades were compared to those of the gpdh nuclear gene, offering a comparative measurement of the selective constraints. although active elements might suffer selective pressure to maintain their coding sequences, gpdh would be expected to be even more restricted because it plays an essential role in glycerophospholipid metabolism in drosophila. if it is assumed that synonymous substitutions are under almost strictly neutral evolution, dn:ds = dn:ds < and dn:ds > will represent neutral evolution, purifying selection and positive darwinian selection, respectively. the analyses indicate that the mean gypsy dn:ds ratios ranged from to, while the mean gpdh values did not exceed. these ratios suggest that, although dn:ds < purifying selection over the gypsy env sequences is relaxed compared to the host gene. also, a z-test for the gypsy env region indicates neutrality for most of the pairwise comparisons within the clades. hence, high selective constraints do not explain the incongruities observed intra-clades. on the other hand, the neutrality hypothesis was refuted for all pairwise comparisons between the clades, indicating that purifying selection could plays a role in the conservation of the env coding sequence among the gypsy subfamilies. in that case, natural selection could act at the genomic level, conserving the te genes, since elements that transpose most efficiently, or at the highest rate are most likely to survive and propagate. note: pairwise ds and dn distances in additional files and in the absence of intra-clade selective constraints, the comparisons between the ds values of the te and a host gene permit evaluation of the hypothesis of ancestral polymorphism within subfamilies. if the gypsy retroelement is evolving vertically and the incongruence in the phylogeny are due to differential subfamily distribution in the analyzed species, the ds pairwise comparisons between the te and the host gene would be expected to be equivalent. but, if the ds of te is lower than the ds of host gene under similar or higher levels of selective constraints and in the absence of strong codon usage bias, ht events can be inferred. using this approach, the ds of gypsy and gpdh gene, as well as the codon bias index were compared for one clone representing each species of zaprionus and the best blastn subject of the genomic sequences in each clade. codon usage bias is a factor that may be responsible for low ds values for te sequences, since there is a negative correlation between codon usage bias and ds values. the codon bias index indicates that the gypsy and gpdh sequences did not suffer high codon bias, although the mean value of for gpdh suggests that the small difference in magnitude of ds between gpdh and gypsy ought to be taken cautiously in inferring ht, since the variation in the gpdh ds distance could reflect codon usage. hence, comparisons between gypsy and gpdh ds were used to infer ht only if the ds values of gypsy were at least twice as low as those of gpdh. the gypsy ds values between zaprionus and melanogaster species ranged from to times lower than those of gpdh. on the other hand, comparisons of ds between the zaprionus species within each clade gave gypsy ds values approximately equal or higher than gpdh ds, excepting for z. davidi vs z. indianus and z. africanus which presented ds of gypsy slightly lower than the host gene ones. these results may indicate that the gypsy retroelement evolved mainly by vertical transmission within the zaprionus subgenus, but there may have been ht events between the zaprionus and the melanogaster species what could reject the hypothesis of ancestral polymorphism within the clades. phylogenetic relationships of the micropia retroelement in order to evaluate the evolutionary history of micropia, the same approaches of phylogeny and divergence estimations used for gypsy were applied. sequences homologous to the zaprionus species were also identified only in species of the melanogaster subgroup, except for d. erecta, the unique species of the melanogaster subgroup for which no significant matches were obtained. d. sechellia had the most hits, followed by d. simulans, d. yakuba and d. melanogaster. structure analysis indicates that the d. melanogaster and d. simulans genomes harbor putatively full-length insertions. inferring phylogenetic relationships and divergence the mp, nj and ml reconstructions of the micropia sequences from the zaprionus subgenus, the genomic sequences from the melanogaster subgroup and the genbank sequences of the repleta group had similar topologies. as for gypsy, only ml reconstruction is given. the tree exhibits several incongruences according to species phylogeny, such as grouping species in the zaprionus genus together with those of the melanogaster subgroup. the sequences of the micropia family of the repleta group were clustered outside a well-supported branch that harbors all sequences of the zaprionus subgenus and the melanogaster species subgroup, constituted by two clades. clade with low support, grouped the d. sechelliaa sequence with three internal clades, corresponding to the z. camerounensis sequences, d. yakuba and the sequences of z. davidi, z. indianus, z. africanus and z. gabonicus. this last clade did not give a clear internal species-specific clustering. clade grouped sequences of the melanogaster subgroup together with z. tuberculatus in 76% of the replications. despite the clustering of all melanogaster sequences in clade the support was not high. the mean divergence within clade was much lower than in clade or between them. analyzing the species groups separately, the mean divergence of the zaprionus sequences was; it was within the melanogaster species subgroup, and between the zaprionus and the melanogaster sequences. the phylogenetic tree and distance estimates <20% indicate that micropia in the zaprionus and melanogaster species belong to the same subfamily. testing the ht hypothesis evaluation of the selective constraints showed that the mean dn:ds values of the micropia sequences were 14- times higher than those of gpdh. between the zaprionus and melanogaster sequences, for example, the values were and for micropia and gpdh, respectively. in addition, the z-test failed to reject the neutrality hypothesis in ~80% of the pairwise comparisons. these results refute the hypothesis that the high similarity between the micropia sequences of zaprionus and melanogaster is due to high selective constraints conserving their nucleotide sequences. to evaluate the ancestral polymorphism hypothesis for micropia, the ds values of micropia rnaseh and gpdh sequences were compared between the sequences of zaprionus and melanogaster species in each clade. all the ds distances for the micropia were lower than those for the gpdh. for instance, the ds between z. tuberculatus and the melanogaster subgroup sequences in clade were from to times lower than that those for the gpdh. those for the d. yakubaa sequence were from to times lower. in addition, ds was zero in some comparisons, e.g. those among the z. camerounensis, z. davidi, z. africanus and z. gabonicus sequences of clade the high conservation in the synonymous sites of the rnaseh sequences favors the hypothesis that their low divergence is due to recent acquisitions of micropia in the zaprionus and melanogaster genomes, and that ht has shaped its evolution, as already proposed for the gypsy retroelement. network trees corroborate the inference of ht network reconstructions were used as an additional approach to test the hypothesis that ht has shaped the evolution of gypsy and micropia in zaprionus and melanogaster species. since te sequences can be considered as populations of sequences that share a common ancestor, this type of phylogeny lets one propose possible donors and receptors species, and to reconstruct the te dispersion routes by the proposition of the ancestral nodes. the gypsy network indicated three clusters of sequences, which correspond to the s subfamily and the two variants of the s subfamily. moreover, the s subfamily was closer to s2v than to s2v which indicates that s2v diversified from s2v the presence of median vectors - that correspond to a hypothetical or unsampled sequences - connecting the zaprionus and melanogaster sequences suggests that the ht events occurred in the ancestors of the sequences sampled, indicating that the lateral transmissions were not recent. the network tree allowed us to see a possible ht not clearly showed by the traditional phylogeny, but evident from the ds comparisons. the z. davidi sequence is the ancestral node of the z. africanus sequence, suggesting a single ht event from z. davidi to z. africanus. the network reconstruction also let the d. yakubaa sequence be included as a s2v member. network reconstruction for micropia clustered the sequences of d. yakuba, d. melanogaster, z. camerounensis, z. davidi, z. indianus and z. africanus with central median vectors, suggesting that high sequence conservation was not due to hts involving the species analyzed, but probably one or more donor species not sampled. the d. melanogastera sequence connects the median vectors with a large cluster that harbors d. sechellia, d. simulans and d. melanogaster sequences, which in turn links a median branch to the z. tuberculatus sequences. this arrangement indicates at least two more recent transfers, one from d. melanogaster to d. simulans and/or d. sechellia or even to the ancestor of those sister species, and a second from any of these species to z. tuberculatus. in addition, the branch of the d. sechelliaa sequence with a median vector that links the z. camerounensis, z. davidi, z. indianus and z. africanus sequences suggest an additional ht involving the ancestor of d. sechellia. however, the absence of geographic sharing between d. sechellia and zaprionus species means that this proposal should be viewed cautiously. an alternative explanation could be a high rate of evolution of the d. sechelliaa sequence, followed by random convergent evolution. ht inference and estimation of the divergence time between te sequences to reinforce the ht hypothesis, the divergence times of both retroelements were calculated using the molecular clock. the divergence estimates for gypsy sequences showed that the hts may have started around mya. regarding the s subfamily, two hts can be proposed. the first may have occurred from a d. melanogaster/d. simulans ancestor to z. tuberculatus about - mya, prior to the time of d. melanogaster/d. simulans splitting. this event was followed by a transfer between d. melanogaster and d. erecta, species that share a common ancestor at 8- mya. although an alternative hypothesis of ht between d. simulans and d. erecta can be formulated on the basis of the phylogeny branching, the greater difference in ds comparisons between the pair d. erecta/d. melanogaster than the pair d. erecta/d. simulans reinforces the first hypothesis. symbols: *: species ancestor. #: time of divergence calculated using ds pairwise comparisons between sequences connected by the network phylogeny. the gypsy s2v sequences seem to have been involved at least in four ht events. we suggest that d. erecta was the donor of s2v in a more ancient transfer to z. tuberculatus, followed by introduction into the z. camerounensis/z. davidi ancestor, and lately into d. yakuba and d. melanogaster. the transfer to z. camerounensis/z. davidi probably occurred in their ancestor, because these species diverged around mya, but the hypothesis of two independent hts cannot be disregarded. the donor status of d. erecta was based on their basal positioning in the network, the repeated clustering to all the other species of this variant in the phylogeny, the divergence times between d. erecta and d. melanogaster/d. yakuba. the s2v variant is probably involved in one ht event that happened around - mya. it could have occurred between the ancestor of the modern sequences of d. melanogaster/d. simulans/d. sechellia and the indianus/davidi complexes. alternatively, the ht could have happened in each complex separately, since the davidi complex diverged about mya. the ht between the d. melanogaster/d. simulans/d. sechellia ancestor and the indianus complex ancestor is a broader scenario than the ht proposed between z. indianus and d. simulans, exemplifying the importance of analyzing as many species as possible before inferring horizontal transfers. finally, the time of divergence also corroborates the ht from z. davidi to z. africanus s2v sequences, as seen in the network tree, which probably occurred ~ mya, subsequent to the divergence of these species. comparing the divergence times between the micropia retroelement of the zaprionus and melanogaster sequences with those of gypsy, the hts of micropia might have started more recently than those of gypsy. the invariance of the synonymous sites between several pairwise comparisons makes it difficult to date these events, but the highest ds levels indicate that, of the species we have analyzed, d. yakuba and the indianus complex species were the first to be invaded by micropia. so, the evolutionary analyses lead to the inference that micropia diverged after the splitting of their species hosts, suggesting that at least an unknown species might have transferred micropia sequences independently in five hts to d. melanogaster, d. yakuba, z. camerounensis, z. davidi and the indianus complex. species of neither the repleta nor the cardini species group are good candidates, since they harbor micropia retrotransposons of a different family. hence, the probable donor might have been a non-studied species. based mainly on the network tree, two other ht events can be proposed. d. melanogaster may have been the donor of micropia to the d. simulans and d. sechellia, or even to their ancestor. however, introgression cannot be excluded in this case. in turn, one of these three species might have been involved in the transmission of micropia to z. tuberculatus, about - mya. discussion this survey of the gypsy and micropia in seven species of zaprionus indicates that they are widely distributed and transcriptionally active only in species of subgenus zaprionus. additionally, the evolutionary analyses demonstrate that these two retroelements are closely related to those of the melanogaster species subgroup, and that their histories might have been repeatedly marked by events of inter- and intra-subgenus hts. several pieces of evidence have shown that the phenomenon of ht is frequent in eukaryotes. it is important to note that the hts proposed here were assumed only when the species involved in the event shared geographic, temporal and ecological environments, and when three separate pieces of evidence suggested their occurrence: lower ds values between the te sequences than the gpdh sequences, incongruence between host and te phylogenies, and last common ancestor of te more recent than that of the species. the studies of micropia and gypsy in subgenus zaprionus species highlight the ht events as an alternative evolutionary mechanism of the retroelements evolution. among the hts inferred for gypsy, four events seemed to involve the zaprionus and melanogaster species, three involved only the melanogaster group species, and only one the zaprionus species. for the micropia retrotransposon, seven hts were suggested: five between unknown donor and the zaprionus or the melanogaster species, one among the melanogaster species, and one between the zaprionus and the melanogaster species. the hts we have proposed probably occurred during two different waves of invasions. first, the gypsy retroelement of the melanogaster subgroup had been introduced in the zaprionus subgenus ~ mya. after that, a second ht wave involving the micropia retroelement may have introduced it into both the zaprionus subgenus and melanogaster subgroup from an unknown donor. after the introduction, the tes could have been prone to transposition and re-introduction in other related species. for example, z. davidi has donated their gypsy sequence to z. africanus after had received it from an ancestor of d. melanogaster, d. simulans and d. sechellia. according to the te evolutionary cycle, the te history could start from ht events, followed by the initial transpositional burst in the new host, and then the accumulation of defective copies along with host-directed epigenetic silencing. the last step could lead to the loss of mobility and, finally, to molecular erosion by random mutations. the high frequency of hts, combined with the transcriptional activity and low insertion numbers of gypsy and micropia retroelements in the zaprionus species subgenus indicates that these elements remain at the first or second steps of the evolutionary life history of the tes, i.e. the invasion or the genomic spread period. a question that arises is whether these elements have been transferred in similar ways. hts of gypsy are easier to explain, since the active gypsy env can enable its dispersion and expansion of the genomic territory, but this mechanism does not to apply to micropia, which lacks env. nevertheless, retroelements could have been transferred concomitantly with retrovirus infections, even gypsy; or parasitic infestations, as for example by intracellular wolbachia bacteria or mites. moreover, our data allowed us to propose there were ht of gypsy and micropia among sister species such as d. melanogaster, d. simulans and d. sechellia, although introgression could not be ruled out. another concern is the evolutionary process of fixation of the tes after the horizontal introduction. these elements could share similar mechanisms of regulation, which would permit them escaping from selective forces and allow fixation after invasion. control of drosophila retroelement mobilization depends on the transcription rate, which is directly related to the presence of specific regulatory proteins and the rate of rna degradation, mediated by rnai systems. however, genetic information about the species of the genus zaprionus is scarce. further, studies of the zaprionus genes involved in controlling mobilization could lead to hypotheses about the mechanisms that allow highly related sequences to be shared in the zaprionus and melanogaster species genomes.  our survey suggests that zaprionus species have experienced waves of retroelement invasions, particularly during the last million years. species of the melanogaster subgroup might have donated their gypsy sequences to species of the zaprionus subgenus, and these two species groups could have received the micropia retroelement from one or more unknown donor species. after the initial introduction, these species could share their elements. since the zaprionus genus and the melanogaster subgroup seem to share the same age of origin and diversification in tropical africa, as well as ecological features, our data suggest that they passed through a permissive period of transposable element invasion during the diversification period. 
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. the gypsy retroelement of the melanogaster subgroup probably invaded the zaprionus genomes about mya. in contrast, the micropia retroelement may have been introduced into the zaprionus subgenus and the melanogaster subgroup from an unknown donor more recently. the sequences of each retroelement were closely related to those of the melanogaster species subgroup, and the most parsimonious hypothesis would be that horizontal transfer events shaped their evolution. gypsy and micropia of zaprionus and melanogaster species share similar evolutionary patterns. 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 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.
 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.  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 horizontal transfer events shaped their evolution. the gypsy retroelement of the melanogaster subgroup probably invaded the zaprionus genomes about mya. in contrast, the micropia retroelement may have been introduced into the zaprionus subgenus and the melanogaster subgroup from an unknown donor more recently.  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. 