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.     

Polymorphism in structure of the retrotransposable element 412 in Drosophila simulans and D. melanogaster populations.

The structure of the 412 retrotransposable element was investigated in various natural populations of D. melanogaster and D. simulans by a restriction enzyme analysis. We show that although the canonical structure of the 412 element was the same in both species, a high structural polymorphism existed with various rearranged elements. A 412 family was thus composed of heterogeneous copies of different sizes, with a large proportion of full-size copies. D. simulans had more rearranged copies than D. melanogaster, with some specific copies, such as a 5.6-kb BsrBI fragment, present in all populations of D. simulans. Full-size and rearranged copies were detected in both the euchromatin and the heterochromatin, with many rearranged copies in D. simulans, suggesting a recent mobilization of the 412 element in this 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.     

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.     

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.      

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.      

Molecular variation at the In(2L)t proximal breakpoint site in natural populations of Drosophila melanogaster and D. simulans

A previous study of nucleotide polymorphism in a Costa Rican population of Drosophila melanogaster found evidence for a nonneutral deficiency in the number of haplotypes near the proximal breakpoint of In(2L)t, a common inversion polymorphism in this species. Another striking feature of the data was a window of unusually high nucleotide diversity spanning the breakpoint site. To distinguish between selective and neutral demographic explanations for the observed patterns in the data, we sample alleles from three additional populations of D. melanogaster and one population of D. simulans. We find that the strength of associations among sites found at the breakpoint varies between populations of D. melanogaster. In D. simulans, analysis of the homologous region reveals unusually elevated levels of nucleotide polymorphism spanning the breakpoint site. As with American populations of D. melanogaster, our D. simulans sample shows a marked reduction in the number of haplotypes but not in nucleotide diversity. Haplotype tests reveal a significant deficiency in the number of haplotypes relative to the neutral expectation in the D. simulans sample and some populations of D. melanogaster. At the breakpoint site, the level of divergence between haplotype classes is comparable to interspecific divergence. The observation of interspecific polymorphisms that differentiate major haplotype classes in both species suggests that haplotype classes at this locus are considerably old. When considered in the context of other studies on patterns of variation within and between populations of D. melanogaster and D. simulans, our data appear more consistent with the operation of selection than with simple demographic explanations.     

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.     

DNA variability and divergence at the notch locus in Drosophila melanogaster and D. simulans: a case of accelerated synonymous site divergence

DNA diversity in two segments of the Notch locus was surveyed in four populations of Drosophila melanogaster and two of D. simulans. In both species we observed evidence of non-steady-state evolution. In D. simulans we observed a significant excess of intermediate frequency variants in a non-African population. In D. melanogaster we observed a disparity between levels of sequence polymorphism and divergence between one of the Notch regions sequenced and other neutral X chromosome loci. The striking feature of the data is the high level of synonymous site divergence at Notch, which is the highest reported to date. To more thoroughly investigate the pattern of synonymous site evolution between these species, we developed a method for calibrating preferred, unpreferred, and equal synonymous substitutions by the effective (potential) number of such changes. In D. simulans, we find that preferred changes per "site" are evolving significantly faster than unpreferred changes at Notch. In contrast we observe a significantly faster per site substitution rate of unpreferred changes in D. melanogaster at this locus. These results suggest that positive selection, and not simply relaxation of constraint on codon bias, has contributed to the higher levels of unpreferred divergence along the D. melanogaster lineage at Notch.     

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.     

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.     

Wolbachia Infections in Drosophila Melanogaster and D. Simulans: Polymorphism and Levels of Cytoplasmic Incompatibility

Wolbachia are endosymbiotic bacteria, widespread in terrestrial Arthropods. They are mainly transmitted vertically, from mothers to offspring and induce various alterations of their hosts' sexuality and reproduction, the most commonly reported phenomenon being Cytoplasmic Incompatibility (CI), observed in Drosophila melanogaster and D. simulans. Basically, CI results in a more or less intense embryonic mortality, occurring in crosses between males infected by Wolbachia and uninfected females. In D. simulans, Wolbachia and CI were observed in 1986. Since then, this host species has become a model system for investigating the polymorphism of Wolbachia infections and CI. In this review we describe the different Wolbachia infections currently known to occur in D. melanogaster and D. simulans. The two species are highly contrasting with regard to symbiotic diversity: while five Wolbachia variants have been described in D. simulans natural populations, D. melanogaster seems to harbor one Wolbachia variant only. Another marked difference between these two Drosophila species is their permissiveness with regard to CI, which seems to be fully expressed in D. simulans but partially or totally repressed in D. melanogaster, demonstrating the involvement of host factors in the control of CI levels. The potential of the two host species regarding the understanding of CI and its evolution is also discussed.     

Historicity and the Population Genetics of Drosophila Melanogaster and D. Simulans

We summarize data showing that there is population structure in African populations of Drosophila from the melanogaster-simulans complex. In D. melanogaster, population structuring is found at individual loci, but is obscured by population structuring for large inversions that simultaneously affect several loci. In D. simulans, molecular polymorphism at the X-linked vermilion locus suggests that different groups of populations have been geographically isolated for some time. Invading populations are probably derived from different areas in Africa. European populations originate from an east African population that was probably not at a demographic equilibrium. The origin of the Antilles population is apparently different and is as yet unknown. In south-western France, populations from these two species undergo different population structuring at the scale of a few kilometres: D. melanogaster makes up a large panmictic population, whereas D. simulans forms a metapopulation that is divided into smaller demes.     

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.     

A Snapshot of Histone Modifications within Transposable Elements in Drosophila Wild Type Strains

Transposable elements (TEs) are a major source of genetic variability in genomes, creating genetic novelty and driving genome evolution. Analysis of sequenced genomes has revealed considerable diversity in TE families, copy number, and localization between different, closely related species. For instance, although the twin species Drosophila melanogaster and D. simulans share the same TE families, they display different amounts of TEs. Furthermore, previous analyses of wild type derived strains of D. simulans have revealed high polymorphism regarding TE copy number within this species. Several factors may influence the diversity and abundance of TEs in a genome, including molecular mechanisms such as epigenetic factors, which could be a source of variation in TE success. In this paper, we present the first analysis of the epigenetic status of four TE families (roo, tirant, 412 and F) in seven wild type strains of D. melanogaster and D. simulans. Our data shows intra- and inter-specific variations in the histone marks that adorn TE copies. Our results demonstrate that the chromatin state of common TEs varies among TE families, between closely related species and also between wild type strains.     

Ecological and Genetic Interactions in Drosophila–parasitoids Communities: A Case Study with D. Melanogaster, D. Simulans and their Common Leptopilina Parasitoids in Southe-astern France

Drosophila species are attacked by a number of parasitoid wasps, which constitute an important factor of population regulation. Since Drosophila melanogaster and Drosophila simulans share common parasitoid species, their ecology and evolution can hardly be understood without considering parasitoids. After a short review of data available on Drosophila-parasitoid interactions involving D. melanogaster and D. simulans as hosts, we report field and laboratory experiments investigating the ecological role of Leptopilina parasitoids in Drosophila communities of southern France. Seasonal survey of species abundance shows that strong interspecific interactions occur at both tropic levels. D. simulans progressively replaces D. melanogaster in southern areas suggesting competitive displacement. Parasitoids are responsible for very high Drosophila mortality (up to 90% in some fruits). Field data emphasize the importance of selective pressure that parasitoids exert on Drosophila communities. The two Leptopilina parasites (L. heterotoma and L boulardi) have different local abundances, which vary in time, and they also compete for hosts. We show that parasitoids can mediate the coexistence of D. melanogaster and D. simulans in the laboratory, and thus may contribute to their puzzling coexistence in the field. Conversely, hosts exert selective pressures on parasitoids, and development on either D. melanogaster or D. simulans strongly affects fitness of adult wasps in a temperature-dependent fashion. Local variation in host species abundance and diversity could thus account for the genetic differentiation we observed in one parasitoid species. Despite laboratory studies cannot fully explain complex field situations, it is clear that the ecology and evolution of Drosophila populations and communities, especially D. melanogaster and D. simulans, are strongly constrained by parasitoids, which should receive more attention.