Distribution of the retrotransposable element 412 in Drosophila species

Copy numbers of sequences homologous to the Drosophila melanogaster retrotransposable element 412, their distribution between the chromosome arms and the chromocenter, and whether they contain full- size copies were analyzed for 55 species of the Drosophila genus. Element 412 insertion sites were detected on the chromosome arms of D. melanogaster, Drosophila simulans, and a few species of the obscura group, but the chromocenter was labeled in almost all species. The presence of element 412 sequences in the majority of species shows that this element has a long evolutionary history in Drosophilidae, although it may have recently invaded the chromosomes in some species, such as D. simulans. Differences in copy number between species may be due to population size or specific endogenous or environmental factors and may follow the worldwide invasion of the species. Putative full-length copies were detected in the chromocenters of some species with no copies on the chromosome arms, suggesting that the chromocenter may be a shelter for such copies and not only for deleted ones.      

Maintenance of transposable element copy number in natural populations of Drosophila melanogaster and D. simulans

To investigate the main forces controlling the containment of transposable elements (TE) in natural populations, we analyzed the copia, mdg1, and 412 elements in various populations of Drosophila melanogaster and D. simulans. A lower proportion of insertion sites on the X chromosome in comparison with the autosomes suggests that selection against the detrimental effects of TE insertions is the major force containing TE copies in populations of Drosophila. This selection effect hypothesis is strengthened by the absence of the negative correlation between recombination rate and TE copy number along the chromosomes, which was expected under the alternative ectopic exchange model (selection against the deleterious rearrangements promoted by recombination between TE insertions). A cline in 412 copy number in relation to latitude was observed among the natural populations of D. simulans, with very high numbers existing in some local populations (around 60 copies in a sample from Canberra, Australia). An apparent absence of selection effects in this Canberra sample and a value of transposition rate equal to 1–2 × 10^-3 whatever the population and its copy number agree with the idea of recent but temporarily drastic TE movements in local populations. The high values of transposition rate in D. simulans clearly disfavor the hypothesis that the low amount of transposable elements in this species could result from a low transposition rate. This revised version was published online in August 2006 with corrections to the Cover Date.     

Transposition rate of the 412 retrotransposable element is independent of copy number in natural populations of Drosophila simulans

The transposition and excision rates of the 412 retrotransposable element were estimated in five populations of Drosophila simulans differing in their average 412 copy numbers, which ranged from 2 to 54. The transposition rate was found to equal 1 x 10(-3) to 2 x 10(-3), independently of copy number. No excision was detected. These values eliminate autoregulation as a force maintaining copy number of the 412 element in natural populations of D. simulans.      

Ancestral polymorphism and recent invasion of transposable elements in Drosophila species

ABSTRACT: BACKGROUND: During the evolutionary history of transposable elements, some processes, such as ancestral polymorphisms and horizontal transfer of sequences between species, can produce incongruences in phylogenies. We investigated the evolutionary history of the transposable elements Bari and 412 in the sequenced genomes of the Drosophila melanogaster group and in the sibling species D. melanogaster and D. simulans using traditional phylogenetic and network approaches. RESULTS: The maximum likelihood (ML) phylogenetic analyses revealed incongruences and unresolved relationships for both the Bari and 412 elements. The DNA transposon Bari within the D. ananassae genome is more closely related to the element of the melanogaster complex than to the sequence in D. erecta, which is inconsistent with the species phylogeny. Divergence analysis and the comparison of the rate of synonymous substitutions per synonymous site of the Bari and host gene sequences explain the incongruence as an ancestral polymorphism inherited stochastically by the derived species. Unresolved relationships were observed in the ML phylogeny of both elements involving D. melanogaster, D. simulans and D. sechellia. A network approach was used to attempt to resolve these relationships. The resulting tree suggests recent transfers of both elements between D. melanogaster and D. simulans. The divergence values of the elements between these species support this conclusion. CONCLUSIONS: We showed that an ancestral polymorphism and recent invasion of genomes due to introgression or horizontal transfer between species occurred during the evolutionary history of the Bari and 412 elements in the melanogaster group. These invasions likely occurred in Africa during the Pleistocene, before the worldwide expansion of D. melanogaster and D. simulans.     

Wake up of transposable elements following Drosophila simulans worldwide colonization.

Transposable elements (TEs) make up around 10%-15% of the Drosophila melanogaster genome, but its sibling species Drosophila simulans carries only one third as many such repeat sequences. We do not, however, have an overall view of copy numbers of the various classes of TEs (long terminal repeat [LTR] retrotransposons, non-LTR retrotransposons, and transposons) in genomes of natural populations of both species. We analyzed 34 elements in individuals from various natural populations of these species. We show that D. melanogaster has higher average chromosomal insertion site numbers per genome than D. simulans for all TEs except five. The LTR retrotransposons gypsy, ZAM, and 1731 and the transposon bari-1 present similar low copy numbers in both species. The transposon hobo has a large number of insertion sites, with significantly more sites in D. simulans. High variation between populations in number of insertion sites of some elements of D. simulans suggests that these elements can invade the genome of the entire species starting from a local population. We propose that TEs in the D. simulans genome are being awakened and amplified as they had been a long time ago in D. melanogaster.     

Comparative analysis of transposable elements in the melanogaster subgroup sequenced genomes

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

A Reanalysis of Protein Polymorphism in Drosophila Melanogaster, D. Simulans, D. Sechellia and D. Mauritiana: Effects of Population Size and Selection

Comparison of synonymous and nonsynonymous variation/substitution within and between species at individual genes has become a widely used general approach to detect the effect of selection versus drift. The sibling species group comprised of two cosmopolitan (Drosophila melanogaster and Drosophila simulans) and two island (Drosophila mauritiana and Drosophila sechellia) species has become a model system for such studies. In the present study we reanalyzed the pattern of protein variation in these species, and the results were compared against the patterns of nucleotide variation obtained from the literature, mostly available for melanogaster and simulans. We have mainly focused on the contrasting patterns of variation between the cosmopolitan pair. The results can be summarized as follows: (1) As expected the island species D. mauritiana and D. sechellia showed much less variation than the cosmopolitan species D. melanogaster and D. simulans. (2) The chromosome 2 showed significantly less variation than chromosome 3 and X in all four species which may indicate effects of past selective sweeps. (3) In contrast to its overall low variation, D. mauritiana showed highest variation for X-linked loci which may indicate introgression from its sibling, D. simulans. (4) An average population of D. simulans was as heterozygous as that of D. melanogaster (14.4% v.s. 13.9%) but the difference was large and significant when considering only polymorphic loci (37.2% v.s. 26.1%). (5) The species-wise pooled populations of these two species showed similar results (all loci = 18.3% v.s. 20.0%, polymorphic loci = 47.2% v.s. 37.6%). (6) An average population of D. simulans had more low-frequency alleles than D. melanogaster, and the D. simulans alleles were found widely distributed in all populations whereas the D. melanogaster alleles were limited to local populations. As a results of this, pooled populations of D. melanogaster showed more polymorphic loci than those of D. simulans (48.0% v.s. 32.0%) but the difference was reduced when the comparison was made on the basis of an average population (29.1% v.s. 21.4%). (7) While the allele frequency distributions within populations were nonsignificant in both D. melanogaster and D. simulans, melanogaster had fewer than simulans, but more than expected from the neutral theory, low frequency alleles. (8) Diallelic loci with the second allele with a frequency less than 20% had similar frequencies in all four species but those with the second allele with a frequency higher than 20% were limited to only melanogaster the latter group of loci have clinal (latitudinal) patterns of variation indicative of balancing selection. (9) The comparison of D. simulans/D. melanogaster protein variation gave a ratio of 1.04 for all loci and 1.42 for polymorphic loci, against a ratio of approximately 2-fold difference for silent nucleotide sites. This suggests that the species ratios of protein and silent nucleotide polymorphism are too close to call for selective difference between silent and allozyme variation in D. simulans. In conclusion, the contrasting levels of allozyme polymorphism, distribution of rare alleles, number of diallelic loci and the patterns of geographic differentiation between the two species suggest the role of natural selection in D. melanogaster, and of possibly ancient population structure and recent worldwide migration in D. simulans. Population size differences alone are insufficient as an explanation for the patterns of variation between these two species.     

Geographical variation in insertion site number of retrotransposon 412 inDrosophila simulans

The insertion site number of the retrotransposable element 412 was analyzed in natural populations ofDrosophila simulans of worldwide origin by in situ hybridization. We observe a gradient in copy number ranging from as high as 23 in Europe to 1–10 in South Africa, while populations in Madagascar and the Indian Islands, which are the cradle ofD. simulans, have only 3–7 copies. We find very different copy numbers in some local populations of Australia and the Pacific Islands (with around 60 copies in 1 sample and only 5 in another), suggesting spontaneous transposition bursts in local populations. Such bursts occurring now and then in local natural populations followed by fly migration could lead to the progressive invasion of the entire species by the transposable element mobilized, explaining the gradient in 412 copy number between northern and southern hemispheres.     

A temperature cline in copy number for 412 but not roo/B104 retrotransposons in populations of Drosophila simulans.

The copy number of the retrotransposable element 412 of Drosophila simulans from populations collected worldwide shows a negative correlation with minimum temperature. No association was detected for the roo/B104 element. The possibility that selective pressures might regulate the 412 copy number in these natural populations is supported by detection of selection against the detrimental effects of 412 insertions (estimated by the proportion of insertions on the X chromosome in comparison with the autosomes) but not roo/B104. These data reveal different spatial patterns for two element families, and strongly suggest that some factors in the environment, such as temperature, may interfere with the control of retrotransposition, thus affecting important aspects of genomic evolution.     

Sequence Analysis of Active Mariner Elements in Natural Populations of Drosophila Simulans

Active and inactive mariner elements from natural and laboratory populations of Drosophila simulans were isolated and sequenced in order to assess their nucleotide variability and to compare them with previously isolated mariner elements from the sibling species Drosophila mauritiana and Drosophila sechellia. The active elements of D. simulans are very similar among themselves (average 99.7% nucleotide identity), suggesting that the level of mariner expression in different natural populations is largely determined by position effects, dosage effects and perhaps other factors. Furthermore, the D. simulans elements exhibit nucleotide identities of 98% or greater when compared with mariner elements from the sibling species. Parsimony analysis of mariner elements places active elements from the three species into separate groups and suggests that D. simulans is the species from which mariner elements in D. mauritiana and D. sechellia are most likely derived. This result strongly suggests that the ancestral form of mariner among these species was an active element. The two inactive mariner elements sequenced from D. simulans are very similar to the inactive peach element from D. mauritiana. The similarity may result from introgression between D. simulans and D. mauritiana or from selective constraints imposed by regulatory effects of inactive elements.     

Comparative analysis of the localization and mobility of retrotransposons in sibling species Drosophila simulans and Drosophila melanogaster]

The distribution of four retrotransposon families (MDG1, MDG3, MDG4 and copia) on polytene chromosomes of different (from 9 to 15) Drosophila simulans strains is studied. The mean number of MDG1 and copia euchromatic hybridization sites (3 sites for each element) is drastically decreased in D. simulans in comparison with D. melanogaster (24 and 18 sites respectively). The mean number of MDG3 sites of hybridization is 5 in D. simulans against 12 in D. melanogaster. As for MDG4 both species have on the average about 2-3 euchromatic sites. The majority of MDG1 and copia and about a half of MDG3 euchromatic copies are localized in restricted number of sites (hot spots) on D. simulans polytene chromosomes. In D. melanogaster these elements are scattered along the chromosomes though there are some hot spots too. It appears that euchromatic copies of MDG1 and copia are considerably less mobile in D. simulans in contrast to D. melanogaster. Some common hot spots of retrotransposon localization in D. simulans and D. melanogaster were earlier described as intercalary heterochromatin regions in D. melanogaster. The level of interstrain variability of MDG4 hybridization sites is comparable in both species. Comparative blot-analysis of adult and larval salivary gland DNA shows that MDG1 and copia are situated mainly in euchromatic regions of D. melanogaster chromosomes. In D. simulans genome they are located mainly in heterochromatic regions underreplicated in salivary gland polytene chromosomes. There are interspecies differences in the distribution of retrotransposons in beta-heterochromatic chromosome regions.     

African Drosophila melanogaster and D. simulans populations have similar levels of sequence variability, suggesting comparable effective population sizes

Drosophila melanogaster and D. simulans are two closely related species with a similar distribution range. Many studies suggested that D. melanogaster has a smaller effective population size than D. simulans. As most evidence was derived from non-African populations, we readdressed this question by sequencing 10 X-linked loci in five African D. simulans and six African D. melanogaster populations. Contrary to previous results, we found no evidence for higher variability, and thus larger effective population size, in D. simulans. Our observation of similar levels of variability of both species will have important implications for the interpretation of patterns of molecular evolution.     

Molecular identification of the active ninja retrotransposon and the inactive aurora element in Drosophila simulans and D. melanogaster.

How transposable elements evolve is a key facet in understanding of spontaneous mutation and genomic rearrangements in various organisms. One of the best ways to approach this question is to study a newly evolved transposable element whose presence is restricted to a specific population or strain. The retrotransposons ninja and aurora may provide insights into the process of their evolution, because of their contrasting characteristics, even though they show high sequence identity. The ninja retrotransposon was found in a Drosophila simulans strain in high copy number and is potent in transposition. On the other hand, aurora elements are distributed widely among the species belonging to the Drosophila melanogaster species complex, but are immobile at least in D. melanogaster. In order to distinguish the two closely resembled retrotransposons by molecular means, we determined and compared DNA sequence of the elements, and identified characteristic internal deletions and nucleotide substitutions in 5'-long terminal repeats (LTR). Analyses of the structure of ninja homologs and LTR sequences amplified from both genomic and cloned DNA revealed that the actively transposable ninja elements were present only in D. simulans strains, but inactive aurora elements exist in both D. melanogaster and D. simulans.     

hobo transposable elements in Drosophila melanogaster and D. simulans.

Genomic patterns of occurrence of the transposable element hobo are polymorphic in the sibling species Drosophila melanogaster and D. simulans. Most tested strains of both species have apparently complete (3.0 kb) and smaller hobo elements (H lines), but in both species some strains completely lack such canonical hobo elements (E lines). The occurrence of H and E lines in D. simulans as well as in D. melanogaster implies that an hypothesis of recent introduction in the latter species is inadequate to explain the phylogenetic occurrence of hobo. Particular internally deleted elements, the approximately 1.5 kb Th1 and Th2 elements, are abundant in many lines of D. melanogaster, and an analogous 1.1 kb internally deleted element, h del sim, is abundant in most lines of D. simulans. Besides the canonical hobo sequences, both species (and their sibling species D. sechellia and D. mauritiana) have many hobo-hybridizing sequences per genome that do not appear to be closely related to the canonical hobo sequence.     

Genetic variation and differentiation at microsatellite loci in Drosophila simulans. Evidence for founder effects in new world populations.

Drosophila simulans isofemale lines from Africa, South America, and two locations in North America were surveyed for variation at 16 microsatellite loci on the X, second, and third chromosomes, and 18 microsatellites, which are unmapped. D. simulans is thought to have colonized New World habitats only relatively recently (within the last few hundred years). Consistent with a founder effect occurring as colonizers moved into these New World habitats, we find less microsatellite variability in North and South American D. simulans populations than for an African population. Population subdivision as measured at microsatellites is moderate when averaged across all loci (FST = 0.136), but contrasts sharply with previous studies of allozyme variation, which have showed significantly less differentiation in D. simulans than in D. melanogaster. There are substantially fewer private alleles observed in New World populations of D. simulans than seen in a similar survey of D. melanogaster. In addition to possible differences in population size during their evolutionary histories, varying colonization histories or other demographic events may be necessary to explain discrepancies in the patterns of variation observed at various genetic markers between these closely related species.     

Population structure among African and derived populations of Drosophila simulans: evidence for ancient subdivision and recent admixture.

Previous studies based on allozyme variation have found little evidence for genetic differentiation in Drosophila simulans. On the basis of DNA sequence variation at two nuclear loci in four African populations of D. simulans, we show that there is significant structure to D. simulans populations within Africa. Variation at one of the loci, vermilion, appears to be neutral and supports an eastern African origin for European and American populations. Samples from the West Indies, Europe, and North America had a nucleotide diversity lower than that of African populations at vermilion and show nonequilibrium haplotype distributions at both vermilion and G6pd, consistent with a hypothesis of recent bottleneck and possibly also admixture in the history of these populations. Directional selection, previously documented at G6pd, appears to have occurred within the coalescence time of the species, obscuring deep population history.