A comprehensive study of genic variation in natural populations of Drosophila melanogaster. VII. Varying rates of genic divergence as revealed by two-dimensional electrophoresis.

Four sibling species from the melanogaster subgroup (Drosophila melanogaster, D. simulans, D. sechellia, and D. mauritiana) were studied for genetic divergence, by high-resolution two-dimensional protein electrophoresis (2DE) coupled with ultrasensitive silver staining. A total of eight tissues from larval and adult developmental stages representing both gonadal (germ-line) and nongonadal (somatic) tissues were analyzed for protein divergence between species. Close to 400 polypeptides (protein spots) were scored from each tissue and species, and protein divergence was measured on the basis of qualitative differences (presence/absence) of protein spots in pairwise species comparisons. The observed levels of genic divergence varied among tissues and among species. When larval hemolymph proteins (which are known to be highly polymorphic) were excluded, there was no evidence to suggest that either the larval or adult-stage proteins, as a whole, are more diverged than the other; variation between different tissues rather than between developmental stages appears to be the most significant factor affecting genetic divergence between species. The reproductive tissue (testis and accessory gland) showed more divergence than did the nonreproductive tissue; D. melanogaster testis (from both larvae and adult males) showed the highest level of divergence. In view of the previous observation that D. simulans, D. mauritiana, and D. sechellia show similar but significantly less reproductive isolation from each other than from D. melanogaster, the present results suggest a correlation between the levels of reproductive-tract-protein divergence and the degree of reproductive isolation in these species.     

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.     

Effect of Divergence Time and Recombination Rate on Molecular Evolution of <Emphasis Type="Italic">Drosophila</Emphasis> INE-1 Transposable Elements and Other Candidates for Neutrally Evolving Sites

Interspecies divergence of orthologous transposable element remnants is often assumed to be simply due to genetic drift of neutral mutations that occurred after the divergence of the species. However, divergence may also be affected by other factors, such as variation in the mutation rate, ancestral polymorphisms, or selection. Here we attempt to determine the impact of these forces on divergence of three classes of sites that are often assumed to be selectively unconstrained (INE-1 TE remnants, sites within short introns, and fourfold degenerate sites) in two different pairwise comparisons of Drosophila (D. melanogaster vs. D. simulans and D. simulans vs. D. sechellia). We find that divergence of these three classes of sites is strongly influenced by the recombination environment in which they are located, and this is especially true for the closer D. simulans vs. D. sechellia comparison. We suggest that this is mainly a result of the contribution of ancestral polymorphisms in different recombination regions. We also find that intergenic INE-1 elements are significantly more diverged than intronic INE-1 in both pairwise comparisons, implying the presence of either negative selection or lower mutation rates in introns. Furthermore, we show that substitution rates in INE-1 elements are not associated with the length of the noncoding sequence in which they are located, suggesting that reduced divergence in long noncoding sequences is not due to reduced mutation rates in these regions. Finally, we show that GC content for each site within INE-1 sequences has evolved toward an equilibrium value (approximately 33%) since insertion.     

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.     

High level of divergence of male-reproductive-tract proteins, between Drosophila melanogaster and its sibling species, D. simulans

We compared male-reproductive-tract polypeptides of Drosophila melanogaster and D. simulans by using two-dimensional gel electrophoresis. Approximately 64% of male-reproductive-tract polypeptides were identical between two randomly chosen isofemale lines from these two species, compared with 83% identity for third-instar imaginal wing-disc polypeptides. Qualitatively similar differences were found between reproductive tracts and imaginal discs when D. sechellia was compared with D. melanogaster and with D. simulans. When genic polymorphism was taken into account, approximately 10% of male- reproductive-tract polypeptides were apparently fixed for different alleles between D. melanogaster and D. simulans; this proportion is the same as that found for soluble enzymes by one-dimensional gel electrophoresis. Strikingly, approximately 20% of male-reproductive- tract polypeptides of either D. melanogaster or D. simulans had no detectable homologue in the other species. We propose that proteins of the Drosophila male reproductive tract may have diverged more extensively between species than have other types of proteins and that much of this divergence may involve large changes in levels of polypeptide expression.      

An interspecific QTL study of Drosophila wing size and shape variation to investigate the genetic basis of morphological differences

The Drosophila wing has been used as a model in studies of morphogenesis and evolution; the use of such models can contribute to our understanding of mechanisms that promote morphological divergence among populations and species. We mapped quantitative trait loci (QTL) affecting wing size and shape traits using highly inbred introgression lines between D. simulans and D. sechellia, two sibling species of the melanogaster subgroup. Eighteen QTL peaks that are associated with 12 wing traits were identified, including two principal components. The wings of D. simulans and D. sechellia significantly diverged in size; two of the QTL peaks could account for part of this interspecific divergence. Both of these putative QTLs were mapped at the same cytological regions as other QTLs for intraspecific wing size variation identified in D. melanogaster studies. In these regions, one or more loci could account for intra- and interspecific variation in the size of Drosophila wings. Three other QTL peaks were related to a pattern of interspecific variation in wing size and shape traits that is summarized by one principal component. In addition, we observed that female wings are significantly larger and longer than male wings and the second, fourth and fifth longitudinal veins are closer together at the distal wing area. This pattern was summarized by another principal component, for which one QTL was mapped.     

Rapid evolution of two odorant-binding protein genes, Obp57d and Obp57e, in the Drosophila melanogaster species group

Genes encoding odorant-binding protein (OBP) form a large family in an insect genome. Two OBP genes, Obp57d and Obp57e, were previously identified to be involved in host-plant recognition in Drosophila sechellia. Here, by comparing the genomic sequences at the Obp57d/e locus from 27 Drosophila species, we found large differences in gene number between species. Phylogenetic analysis revealed that Obp57d and Obp57e in the D. melanogaster species group arose by gene duplication of an ancestral OBP gene that remains single in the obscura species group. Further gain and loss of OBP genes were observed in several lineages in the melanogaster group. Site-specific analysis of evolutionary rate suggests that Obp57d and Obp57e have functionally diverged from each other. Thus, there are two classes of gene number differences in the Obp57d/e region: the difference of the genes that have functionally diverged from each other and the difference of the genes that appear to be functionally identical. Our analyses demonstrate that these two classes of differences can be distinguished by comparisons of many genomic sequences from closely related species.     

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.      

Hsp70 duplication in the Drosophila melanogaster species group: how and when did two become five?

To determine how the modern copy number (5) of hsp70 genes in Drosophila melanogaster evolved, we localized the duplication events that created the genes in the phylogeny of the melanogaster group, examined D. melanogaster genomic sequence to investigate the mechanisms of duplication, and analyzed the hsp70 gene sequences of Drosophila orena and Drosophila mauritiana. The initial two-to-four hsp70 duplication occurred 10--15 MYA, according to fixed in situ hybridization to polytene chromosomes, before the origin and divergence of the melanogaster and five other species subgroups of the melanogaster group. Analysis of more than 30 kb of flanking sequence surrounding the hsp70 gene clusters suggested that this duplication was likely a retrotransposition. For the melanogaster subgroup, Southern hybridization and an hsp70 restriction map confirmed the conserved number (4) and arrangement of hsp70 genes in the seven species other than D. melanogaster. Drosophila melanogaster is unique; tandem duplication and gene conversion at the derived cluster yielded a fifth hsp70 gene. The four D. orena hsp70 genes are highly similar and concertedly evolving. In contrast, the D. mauritiana hsp70 genes are divergent, and many alleles are nonfunctional. The proliferation, concerted evolution, and maintenance of functionality in the D. melanogaster hsp70 genes is consistent with the action of natural selection in this species.     

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.     

BAC library construction and BAC end sequencing of five Drosophila species: the comparative map with the D. melanogaster genome

We constructed and characterized arrayed bacterial artificial chromosome (BAC) libraries of five Drosophila species (D. melanogaster, D. simulans, D. sechellia, D. auraria, and D. ananassae), which are genetically well characterized in the studies of meiosis, evolution, population genetics, and developmental biology. The BAC libraries comprise 8,000 to 12,500 clones for each species, estimated to cover the most of the genomes. We sequenced both ends of most of these BAC clones with a success rate of 91%. Of these, 53,701 clones consisting of non-repetitive BAC end sequences (BESs) were mapped with reference of the public D. melanogaster genome sequences. The BES mapping estimated that the BAC libraries of D. auraria and D. ananassae covered 47% and 57% of the D. melanogaster genome, respectively, and those of D. melanogaster, D. sechellia, and D. simulans covered 94-97%. The low coverage by BESs of D. auraria and D. ananassae may be due to the high sequence divergence with D. melanogaster. From the comparative BES mapping, 111 possible breakpoints of chromosomal rearrangements were identified in these four species. The breakpoints of the major chromosome rearrangement between D. simulans and D. melanogaster on the third chromosome were determined within 20 kb in 84E and 30 kb in 93E/F. Corresponding breakpoints were also identified in D. sechellia. The BAC clones described here will be an important addition to the Drosophila genomic resources.     

Analysis of nucleotide substitutions of mitochondrial DNAs in Drosophila melanogaster and its sibling species

To study the rate and pattern of nucleotide substitution in mitochondrial DNA (mtDNA), we cloned and sequenced a 975-bp segment of mtDNA from Drosophila melanogaster, D. simulans, and D. mauritiana containing the genes for three transfer RNAs and parts of two protein- coding genes, ND2 and COI. Statistical analysis of synonymous substitutions revealed a predominance of transitions over transversions among the three species, a finding differing from previous results obtained from a comparison of D. melanogaster and D. yakuba. The number of transitions observed was nearly the same for each species comparison, including D. yakuba, despite the differences in divergence times. However, transversions seemed to increase steadily with increasing divergence time. By contrast, nonsynonymous substitutions in the ND2 gene showed a predominance of transversions over transitions. Most transversions were between A and T and seemed to be due to some kind of mutational bias to which the A + T-rich mtDNA of Drosophila species may be subject. The overall rate of nucleotide substitution in Drosophila mtDNA appears to be slightly faster (approximately 1.4 times) than that of the Adh gene. This contrasts with the result obtained for mammals, in which the mtDNA evolves approximately 10 times faster than single-copy nuclear DNA. We have also shown that the start codon of the COI gene is GTGA in D. simulans and GTAA in D. mauritiana. These codons are different from that of D. melanogaster (ATAA).      

Local changes in GC/AT substitution biases and in crossover frequencies on Drosophila chromosomes

I present here evidence of remarkable local changes in GC/AT substitution biases and in crossover frequencies on Drosophila chromosomes. The substitution pattern at 10 loci in the telomeric region of the X chromosome was studied for four species of the Drosophila melanogaster species subgroup. Drosophila orena and Drosophila erecta are clearly the most closely related species pair (the erecta complex) among the four species studied; however, the overall data at the 10 loci revealed a clear dichotomy in the silent substitution patterns between the AT-biased- substitution melanogaster and erecta lineages and the GC-biased-substitution yakuba and orena lineages, suggesting two or more independent changes in GC/AT substitution biases. More importantly, the results indicated a between- loci heterogeneity in GC/AT substitution bias in this small region independently in the yakuba and orena lineages. Indeed, silent substitutions in the orena lineage were significantly biased toward G and C at the consecutive yellow, lethal of scute, and asense loci, but they were significantly biased toward A and T at sta. The substitution bias toward G and C was centered in different areas in yakuba (significantly biased at EG:165H7.3, EG:171D11.2, and suppressor of sable). The similar silent substitution patterns in coding and noncoding regions, furthermore, suggested mutational biases as a cause of the substitution biases. On the other hand, previous study reveals that Drosophila yakuba has about 20-fold higher crossover frequencies in the telomeric region of the X chromosome than does D. melanogaster; this study revealed that the total genetic map length of the yakuba X chromosome was only about 1.5 times as large as that of melanogaster and that the map length of the X-telomeric y-sta region did not differ between Drosophila yakuba and D. erecta. Taken together, the data strongly suggested that an approximately 20- fold reduction in the X-telomeric crossover frequencies occurred in the ancestral population of D. melanogaster after the melanogaster-yakuba divergence but before the melanogaster-simulans divergence.     

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.     

Distribution and conservation of sequences homologous to the 1731 retrotransposon in Drosophila

The distribution of 1731 retrotransposon-hybridizing sequences in the family Drosophilidae has been studied using a 1731 probe from Drosophila melanogaster. Squash blot and Southern blot analyses of 42 species reveal that the 1731 sequences are widespread within both the Sophophora and Drosophila subgenera and are also present in the genera Scaptomyza and Zaprionus. Hence the 1731 retrotransposon family appears to have a long evolutionary history in the Drosophilidae genome. Differences of hybridization signal intensity suggested that the 1731 sequence is well conserved only in the three species most closely related to D. melanogaster (D. simulans, D. mauritiana, and D. sechellia). A survey of insertion sites in numerous different populations of the previous four species by in situ hybridization to polytene chromosomes has shown in all cases both chromocentric hybridizations and a low number of sites (0-5) on the chromosomal arms. This number of sites is among the lowest observed in D. melanogaster and D. simulans when 1731 is compared with other retrotransposon families. In addition, we have observed species-specific patterns of the chromocentric hybridization signal, suggesting rapid modifications of the beta-heterochromatin components since the radiation of the melanogaster subgroup.      

Nup96-dependent hybrid lethality occurs in a subset of species from the simulans clade of Drosophila

The cross of Drosophila melanogaster females to D. simulans males typically produces lethal F(1) hybrid males. F(1) male lethality is suppressed when the D. simulans Lhr(1) hybrid rescue strain is used. Viability of these F(1) males carrying Lhr(1) is in turn substantially reduced when the hybrids are heterozygous for some mutant alleles of the D. melanogaster Nup96 gene. I show here that similar patterns of Nup96-dependent lethality occur when other hybrid rescue mutations are used to create F(1) males, demonstrating that Nup96 does not reduce hybrid viability by suppressing the Lhr(1) rescue effect. The penetrance of this Nup96-dependent lethality does not correlate with the penetrance of the F(1) hybrid rescue, arguing that these two phenomena reflect genetically independent processes. D. simulans, together with two additional sister species, forms a clade that speciated after the divergence of their common ancestor from D. melanogaster. I report here that Nup96(-) reduces F(1) viability in D. melanogaster hybrids with one of these sister species, D. sechellia, but not with the other, D. mauritiana. These results suggest that Nup96-dependent lethality evolved after the speciation of D. melanogaster from the common ancestor of the simulans clade and is caused by an interaction among Nup96, unknown gene(s) on the D. melanogaster X chromosome, and unknown autosomal gene(s), at least some of which have diverged in D. simulans and D. sechellia but not in D. mauritiana. The genetic properties of Nup96 are also discussed relative to other hybrid lethal genes.