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.     

Horizontal transfer of hobo transposable elements within the Drosophila melanogaster species complex: evidence from DNA sequencing.

The hobo family of transposable elements, one of three transposable-element families that cause hybrid dysgenesis in Drosophila melanogaster, appears to be present in all members of the D. melanogaster species complex: D. melanogaster, D. simulans, D. mauritiana, and D. sechellia. Some hobo-hybridizing sequences are also found in the other members of the melanogaster subgroup and in many members of the related montium subgroup. Surveys of older isofemale lines of D. melanogaster suggest that complete hobo elements were absent prior to 50 years ago and that hobo has recently been introduced into the species by horizontal transfer. To test the horizontal transfer hypothesis, the 2.6-kb XhoI fragments of hobo elements from D. melanogaster, D. simulans, and D. mauritiana were cloned and sequenced. The DNA sequences reveal an extremely low level of divergence and support the conclusion that the active hobo element has been horizontally transferred into or among these species in the recent past.     

Distribution of hobo transposable elements in the genus Drosophila

This study describes the distribution of hobo-hybridizing sequences in the genus Drosophila. Southern blot analysis of 134 species revealed that hobo sequences are limited to the melanogaster and montium subgroups of the melanogaster-species group. Of the hobo-bearing species, only D. melanogaster and two of its sibling species, D. simulans and D. mauritiana, were found to contain potentially complete hobo elements. The distribution of hobo sequences is one of the narrowest distributions thus far described for any Drosophila transposable element.     

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.     

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.     

Horizontal transmission, vertical inactivation, and stochastic loss of mariner-like transposable elements

Horizontal transmission has been well documented as a major mechanism for the dissemination of mariner-like elements (MLEs) among species. Less well understood are mechanisms that limit vertical transmission of MLEs resulting in the "spotty" or discontinuous distribution observed in closely related species. In this article we present evidence that the genome of the common ancestor of the melanogaster species subgroup of Drosophila contained an MLE related to the mellifera (honey bee) subfamily. Horizontal transmission, approximately 3-10 MYA, is strongly suggested by the observation that the sequence of the MLE in Drosophila erecta is 97% identical in nucleotide sequence with that of an MLE in the cat flea, Ctenocephalides felis. The D. erecta MLE has a spotty distribution among species in the melanogaster subgroup. The element has a high copy number in D. erecta and D. orena, a moderate copy number in D. teissieri and D. yakuba, and was apparently lost ("stochastic loss") in the lineage leading to D. melanogaster, D. simulans, D. mauritiana, and D. sechellia. In D. erecta, most copies are concentrated in the heterochromatin. Two copies from D. erecta, denoted De12 and De19, were cloned and sequenced, and they appear to be nonfunctional ("vertical inactivation"). It therefore appears that the predominant mode of MLE evolution is vertical inactivation and stochastic loss balanced against occasional reinvasion of lineages by horizontal transmission.      

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.     

Analysis of copia sequence variation within and between Drosophila species

The sequences of the 5' long-terminal repeat (LTR) and adjacent leader regions of 27 full-length copia elements isolated from natural populations of Drosophila melanogaster, D. simulans, and D. mauritiana are presented. Phylogenetic analyses indicate that although D. melanogaster copia elements are distinct from those of D. simulans and D. mauritiana, the elements of these latter two species are not distinguishable from one another. LTRs and adjacent 5' leader regions of elements isolated from D. simulans and D. mauritiana are structurally similar to one another and carry substantial deletional variation mapping to regions previously identified as being of potential importance for copia expression.      

Population genetics and phylogenetics of DNA sequence variation at multiple loci within the Drosophila melanogaster species complex

Two regions of the genome, a 1-kbp portion of the zeste locus and a 1.1- kbp portion of the yolk protein 2 locus, were sequenced in six individuals from each of four species: Drosophila melanogaster, D. simulans, D. mauritiana, and D. sechellia. The species and strains were the same as those of a previous study of a 1.9-kbp region of the period locus. No evidence was found for recent balancing or directional selection or for the accumulation of selected differences between species. Yolk protein 2 has a high level of amino acid replacement variation and a low level of synonymous variation, while zeste has the opposite pattern. This contrast is consistent with information on gene function and patterns of codon bias. Polymorphism levels are consistent with a ranking of effective population sizes, from low to high, in the following order: D. sechellia, D. melanogaster, D.mauritiana, and D. simulans. The apparent species relationships are very similar to those suggested by the period locus study. In particular, D. simulans appears to be a large population that is still segregating variation that arose before the separation of D. mauritiana and D. sechellia. It is estimated that the separation of ancestral D. melanogaster from the other species occurred 2.5-3.4 Mya. The separations of D. sechellia and D. mauritiana from ancestral D. simulans appear to have occurred 0.58- 0.86 Mya, with D. mauritiana having diverged from ancestral D. simulans 0.1 Myr more recently than D. sechellia.      

Sequences homologous to the hobo transposable element in E strains of Drosophila melanogaster

Hobo is one of the three Drosophila melanogaster transposable elements, together with the P and I elements, that seem to have recently invaded the genome of this species. Surveys of the presence of hobo in strains from different geographical and temporal origins have shown that recently collected strains contain complete and deleted elements with high sequence similarity (H strains), but old strains lack hobo elements (E strains). Besides the canonical hobo sequences, both H and E strains show other poorly known hobo-related sequences. In the present work, we analyze the presence, cytogenetic location, and structure of some of these sequences in E strains of D. melanogaster. By in situ hybridization, we found that euchromatic hobo-related sequences were in fixed positions in all six E strains analyzed: 38C in the 2L arm; 42B and 55A in the 2R arm; 79E and 80B in the 3L arm; and 82C, 84C, and 84D in the 3R arm. Sequence comparison shows that some of the hobo-related sequences from Oregon-R and iso-1 strains are similar to the canonical hobo element, but their analysis reveals that they are substantially diverged and rearranged and cannot code for a functional transposase. Our results suggest that these ubiquitous hobo-homologous sequences are immobile and are distantly related to the modern hobo elements from D. melanogaster.     

DNA Sequence Variation at the Period Locus within and among Species of the Drosophila Melanogaster Complex

A 1.9-kilobase region of the period locus was sequenced in six individuals of Drosophila melanogaster and from six individuals of each of three sibling species: Drosophila simulans, Drosophila sechellia and Drosophila mauritiana. Extensive genealogical analysis of 174 polymorphic sites reveals a complex history. It appears that D. simulans, as a large population still segregating very old lineages, gave rise to the island species D. mauritiana and D. sechellia. Rather than considering these speciation events as having produced ``sister' taxa, it seems more appropriate to consider D. simulans a parent species to D. sechellia and D. mauritiana. The order, in time, of these two phylogenetic events remains unclear. D. mauritiana supports a large number of polymorphisms, many of which are shared with D. simulans, and so appears to have begun and persisted as a large population. In contrast, D. sechellia has very little variation and seems to have experienced a severe population bottleneck. Alternatively, the low variation in D. sechellia could be due to recent directional selection and genetic hitchhiking at or near the per locus.     

Intra- and interspecies variation among Bari-1 elements of the melanogaster species group.

We have investigated the distribution of sequences homologous to Bari-1, a Tc1-like transposable element first identified in Drosophila melanogaster, in 87 species of the Drosophila genus. We have also isolated and sequenced Bari-1 homologues from D. simulans, D. mauritiana, and D. sechellia, the species constituting with D. melanogaster the melanogaster complex, and from D. diplacantha and D. erecta, two phylogenetically more distant species of the melanogaster group. Within the melanogaster complex the Bari-1 elements are extremely similar to each other, showing nucleotide identity values of at least 99.3%. In contrast, Bari-1-like elements from D. diplacantha and D. erecta are on average only 70% similar to D. melanogaster Bari-1 and are usually defective due to nucleotide deletions and/or insertions in the ORFs encoding their transposases. In D. erecta the defective copies are all located in the chromocenter and on chromosome 4. Surprisingly, while D. melanogaster Bari-1 elements possess 26-bp inverted terminal repeats, their D. diplacantha and D. erecta homologues possess long inverted terminal repeats similar to the terminal structures observed in the S elements of D. melanogaster and in several other Tc1-like elements of different organisms. This finding, together with the nucleotide and amino acid identity level between D. diplacantha and D. erecta elements and Bari-1 of D. melanogaster, suggests a common evolutionary origin and a rapid diversification of the termini of these Drosophila Tc1-like elements.     

Patterns of Microsatellite Variability in the Drosophila Melanogaster Complex

Forty-seven microsatellite loci were amplified in Drosophila melanogaster, Drosophila simulans, Drosophila mauritiana and Drosophila sechellia. The two cosmopolitan species D. melanogaster and D. simulans were found to be the most variable ones, followed by D. mauritiana and D. sechellia. A model based clustering algorithm was applied to the population samples of D. melanogaster, D. simulans and D. sechellia. No evidence for population substructure was detected within species--most likely due to insufficient power. A Markov chain Monte Carlo method developed for demographic inference based on microsatellites provided unambiguous evidence for population contraction in D. melanogaster, D. simulans and D. sechellia, despite that the D. melanogaster and D. simulans population samples were of non-African origin and represented recently expanded populations.     

Differences in Crossover Frequency and Distribution among Three Sibling Species of Drosophila

Comparisons of the genetic and cytogenetic maps of three sibling species of Drosophila reveal marked differences in the frequency and cumulative distribution of crossovers during meiosis. The maps for two of these species, Drosophila melanogaster and D. simulans, have previously been described, while this report presents new map data for D. mauritiana, obtained using a set of P element markers. A genetic map covering nearly the entire genome was constructed by estimating the recombination fraction for each pair of adjacent inserts. The P-based genetic map of mauritiana is ~1.8 times longer than the standard melanogaster map. It appears that mauritiana has higher recombination along the entire length of each chromosome, but the difference is greatest in centromere-proximal regions of the autosomes. The mauritiana autosomes show little or no centromeric recombinational suppression, a characteristic that is prominent in melanogaster. D. simulans appears to be intermediate both in terms of total map length and intensity of the autosomal centromeric effect. These interspecific differences in recombination have important evolutionary implications for DNA sequence organization and variability. In particular, mauritiana is expected to differ from melanogaster in patterns and amounts of sequence variation and transposon insertions.     

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.     

Evolution of the Transposable Element Mariner in Drosophila Species

The distribution of the transposable element mariner was examined in the genus Drosophila. Among the eight species comprising the melanogaster species subgroup, the element is present in D. mauritiana, D. simulans, D. sechellia, D. yakuba and D. teissieri, but it is absent in D. melanogaster, D. erecta and D. orena. Multiple copies of mariner were sequenced from each species in which the element occurs. The inferred phylogeny of the elements and the pattern of divergence were examined in order to evaluate whether horizontal transfer among species or stochastic loss could better account for the discontinuous distribution of the element among the species. The data suggest that the element was present in the ancestral species before the melanogaster subgroup diverged and was lost in the lineage leading to D. melanogaster and the lineage leading to D. erecta and D. orena. This inference is consistent with the finding that mariner also occurs in members of several other species subgroups within the overall melanogaster species group. Within the melanogaster species subgroup, the average divergence of mariner copies between species was lower than the coding region of the alcohol dehydrogenase (Adh) gene. However, the divergence of mariner elements within species was as great as that observed for Adh. We conclude that the relative sequence homogeneity of mariner elements within species is more likely a result of rapid amplification of a few ancestral elements than of concerted evolution. The mariner element may also have had unequal mutation rates in different lineages.     

Full-size P and KP elements predominate in wild Drosophila melanogaster

We analyzed the genomic P elements of 57 wild-derived Drosophila melanogaster isofemale lines from Africa, Australia and Asia. All carried many P sequences per genome, and the full-size P and the internally deleted KP elements were very common or predominant in the populations. The genomic content of full-size P and KP elements does not correlate well with the P transposition-inducing and -repressing abilities of a line. Our results show that a large majority of type I repressor elements are full-size P elements, and almost all type II repressor elements are KP elements in the natural populations of D. melanogaster from these parts of the world.     

Distribution of transposable elements in neotropical species of Drosophila

The phylogenetic distribution of transposable families, P, gypsy, hobo, I, and mariner has been analyzed in 33 species of 11 groups of neotropical Drosophila and a Drosophilidae species Zygotrica vittimaculosa, using squash blot and dot blot. Genomic DNA of almost all neotropical species tested hybridized with gypsy probe and some species showed a particularly strong hybridization signal, as D. gaucha, D. virilis, and species of flavopilosa group. The hobo element was restricted to melanogaster group and some strains of D. willistoni. Only D. simulans DNA showed hybridization to mariner probe in all species tested and D. simulans and D. melanogaster showed hybridization with I element probe. P element homologous sequence was present in D. melanogaster and all species and strains of the willistoni and saltans groups tested. The presence of at least one P-homologous sequence was detected in Drosophila mediopunctata. This one was the only P-bearing species of all six tested from the tripunctata group. Four different pairs of primers homologous to segments of the canonical sequence of D. melanogaster's P were used to amplify specific sequences from D. mediopunctata DNA, showing the occurrence of seemingly well-conserved P-homologous sequences. This revised version was published online in July 2006 with corrections to the Cover Date.     

Molecular population genetics of male accessory gland proteins in the Drosophila simulans complex

Accessory gland proteins are a major component of Drosophila seminal fluid. These proteins have a variety of functions and may be subject to sexual selection and/or antagonistic evolution between the sexes. Most population genetic data from these proteins are from D. melanogaster and D. simulans. Here, we extend the population genetic analysis of Acp genes to the other simulans complex species, D. mauritiana and D. sechellia. We sequenced population samples of seven Acp's from D. mauritiana, D. sechellia, and D. simulans. We investigated the population genetics of these genes on individual simulans complex lineages and compared Acp polymorphism and divergence to polymorphism and divergence from a set of non-Acp loci in the same species. Polymorphism and divergence data from the simulans complex revealed little evidence for adaptive protein evolution at individual loci. However, we observed a dramatically inflated index of dispersion for amino acid substitutions in the simulans complex at Acp genes, but not at non-Acp genes. This pattern of episodic bursts of protein evolution in Acp's provides the strongest evidence to date that the population genetic mechanisms driving Acp divergence are different from the mechanisms driving evolution at most Drosophila genes.