Local recombination and mutation effects on molecular evolution in Drosophila.

I studied the cause of the significant difference in the synonymous-substitution pattern found in the achaete-scute complex genes in two Drosophila lineages, higher codon bias in Drosophila yakuba, and lower bias in D. melanogaster. Besides these genes, the functionally unrelated yellow gene showed the same substitution pattern, suggesting a region-dependent phenomenon in the X-chromosome telomere. Because the numbers of A/T --> G/C substitutions were not significantly different from those of G/C --> A/T in the yellow noncoding regions of these species, a AT/GC mutational bias could not completely account for the synonymous-substitution biases. In contrast, we did find an approximately 14-fold difference in recombination rates in the X-chromosome telomere regions between the two species, suggesting that the reduction of recombination rates in this region resulted in the reduction of the efficacy of selection in D. melanogaster. In addition, the D. orena yellow showed a 5% increase in the G + C content at silent sites in the coding and noncoding regions since the divergence from D. erecta. This pattern was significantly different from those at the orena Adh and Amy loci. These results suggest that local changes in recombination rates and mutational pressures are contributing to the irregular synonymous-substitution patterns in Drosophila.     

Reproductive Isolation and Morphogenetic Evolution in Drosophila Analyzed by Breakage of Ethological Barriers

This article reports the breaking of ethological barriers through the constitution of soma-germ line chimeras between species of the melanogaster subgroup of Drosophila, which are ethologically isolated. Female Drosophila yakuba and D. teissieri germ cells in a D. melanogaster ovary produced functional oocytes that, when fertilized by D. melanogaster sperm, gave rise to sterile yakuba-melanogaster and teissieri-melanogaster male and female hybrids. However, the erecta-melanogaster and orena-melanogaster hybrids were lethal, since female D. erecta and D. orena germ cells in a D. melanogaster ovary failed to form oocytes with the capacity to develop normally. This failure appears to be caused by an altered interaction between the melanogaster soma and the erecta and orena germ lines. Germ cells of D. teissieri and D. orena in a D. melanogaster testis produced motile sperm that was not stored in D. melanogaster females. This might be due to incompatibility between the teissieri and orena sperm and the melanogaster seminal fluid. A morphological analysis of the terminalia of yakuba-melanogaster and teissieri-melanogaster hybrids was performed. The effect on the terminalia of teissieri-melanogaster hybrids of a mutation in doublesex, a regulatory gene that controls the development of the terminalia, was also investigated.     

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.     

Intron length evolution in Drosophila

I present data on the evolution of intron lengths among 3 closely related Drosophila species, D. melanogaster, Drosophila simulans, and Drosophila yakuba. Using D. yakuba as an outgroup, I mapped insertion and deletion mutations in 148 introns (spanning approximately 30 kb) to the D. melanogaster and D. simulans lineages. Intron length evolution in the 2 sister species has been different: in D. melanogaster, X-linked introns have increased slightly in size, whereas autosomal ones have decreased slightly in size; in D. simulans, both X-linked and autosomal introns have decreased in size. To understand the possible evolutionary causes of these lineage- and chromosome-specific patterns of intron evolution, I studied insertion-deletion (indel) polymorphism and divergence in D. melanogaster. Small insertion mutations segregate at elevated frequencies and enjoy elevated probabilities of fixation, particularly on the X chromosome. In contrast, there is no detectable X chromosome effect on fixations in D. simulans. These findings suggest X chromosome-specific selection or biased gene conversion-gap repair favoring insertions in D. melanogaster but not in D. simulans. These chromosome- and lineage-specific patterns of indel substitution are not easily explained by existing general population genetic models of intron length evolution. Genomic data from D. melanogaster further suggest that the forces described here affect introns and intergenic regions similarly.     

Divergence and heterogeneity of the histone gene repeating units in the Drosophila melanogaster species subgroup

The repeating units of the histone gene cluster containing the H1, H2A, H2B and H4 genes were amplified by PCR from the Drosophila melanogaster species subgroup, i.e., D. yakuba, D. erecta, D. sechellia, D. mauritiana, D. teissieri and D. orena. The PCR products were cloned and their nucleotide sequences of about 4.6-4.8kbp were determined to elucidate the mechanism of molecular evolution of the histone gene family. The heterogeneity among the histone gene repeating units was 0.6% and 0.7% for D. yakuba and D. sechellia, respectively, indicating the same level of heterogeneity as in the H3 gene region of D. melanogaster. Divergence of the genes among species even in the most closely related ones was much greater than the heterogeneity among family members, indicating a concerted mode of evolution for the histone gene repeating units. Among the species in the D. melanogaster species subgroup, the histone gene regions as well as 3rd codon position of the coding region showed nearly the same GC contents. These results suggested that the previous conclusion on analysis of the H3 gene regions, the gene family evolution in a concerted fashion, holds true for the whole histone gene repeating unit.     

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.      

Evolutionary relationships to parasitism by seven species of the Drosophila melanogaster subgroup

Parasitic wasps are an important component of the niche of Drosophila species. The susceptibility to the Cynipid Leptopilina boulardi was estimated in the seven sibling species of Drosophila belonging to the melanogaster subgroup. Three categories of Hies can be distinguished, according to the level of cellular immune reaction and success of parasitism. Drosophila melanogaster and D. mauritiana belong to the category 1, specified by no encapsulative reaction and a high rate of successful parasitism. Category 2, characterized by a moderate encapsulation rate and a high mortality include D. simulans.5, D. erecta and D. orena . Category 3, with D. yakuba and D. tcissien , is specified by a very low rate or absence of successful parasitism due to a highly efficient immune cellular reaction. This classification parallels the phylogenic relationship based upon polytene chromosome banding sequences. Such specific ditferences in susceptibility to parasites may plan an important role in the competition between these species in Africa.     

Changes in the recombinational environment affect divergence in the yellow gene of Drosophila

The complete coding region of the yellow (y) gene was sequenced in different Drosophila species. In the species of the melanogaster subgroup (D. melanogaster, D. simulans, D. mauritiana, D. yakuba, and D. erecta), this gene is located at the tip of the X chromosome in a region with a strong reduction in recombination rate. In contrast, in D. ananassae (included in the ananassae subgroup of the melanogaster group) and in the obscura group species (D. subobscura, D. madeirensis, D. guanche, and D. pseudoobscura), the y gene is located in regions with normal recombination rates. As predicted by the hitchhiking and background selection models, this change in the recombinational environment affected synonymous divergence in the y-gene-coding region. Estimates of the number of synonymous substitutions per site were much lower between the obscura group species and D. ananassae than between the species of the obscura group and the melanogaster subgroup. In fact, a highly significant increase in the rate of synonymous substitution was detected in all lineages leading to the species of the melanogaster subgroup relative to the D. ananassae lineage. This increase can be explained by a higher fixation rate of mutations from preferred to unpreferred codons (slightly deleterious mutations). The lower codon bias detected in all species of the melanogaster subgroup relative to D. ananassae (or to the obscura group species) would be consistent with this proposal. Therefore, at least in Drosophila, changes in the recombination rate in different lineages might cause deviations of the molecular-clock hypothesis and contribute to the overdispersion of the rate of synonymous substitution. In contrast, the change in the recombinational environment of the y gene has no detectable effect on the rate of amino acid replacement in the Yellow protein.     

Microevolutionary divergence pattern of the segmentation gene hunchback in Drosophila

To study the microevolutionary processes shaping the evolution of the segmentation gene hunchback (hb) from Drosophila melanogaster, we cloned and sequenced the gene from 12 isofemale lines representing wild-type populations of D. melanogaster, as well as from the closely related species Drosophila sechellia, Drosophila orena, and Drosophila yakuba. We find a relatively low degree of sequence variation in D. melanogaster (theta = 0.0017), which is, however, consistent with its chromosomal location in a region of low recombination. Tests of neutrality do not reject a neutral-evolution model for the whole region. However, pairwise tests with different subregions indicate that there is a relative excess of polymorphic sites in the leader and the intron. Codon usage pattern analysis shows a particularly biased codon usage in the highly conserved regions, which is in line with the hypothesis that selection on translational accuracy is the driving force behind such a bias. A comparison of the expression pattern of hb in different sibling species of D. melanogaster reveals some regulatory changes in D. yakuba, which could be interpreted as changes in the timing of secondary expression domains.     

Paleo-demography of the Drosophila melanogaster subgroup: application of the maximum likelihood method

The species divergence times and demographic histories of Drosophila melanogaster and its three sibling species, D. mauritiana, D. simulans, and D. yakuba, were investigated using a maximum likelihood (ML) method. Thirty-nine orthologous loci for these four species were retrieved from DDBJ/EMBL/GenBank database. Both autosomal and X-linked loci were used in this study. A significant degree of rate heterogeneity across loci was observed for each pair of species. Most loci have the GC content greater than 50% at the third codon position. The codon usage bias in Drosophila loci is considered to result in the high GC content and the heterogenous rates across loci. The chi-square, G, and Fisher's exact tests indicated that data sets with 11, 23, and 9 pairs of DNA sequences for the comparison of D. melanogaster with D. mauritiana, D. simulans, and D. yakuba, respectively, retain homogeneous rates across loci. We applied the ML method to these data sets to estimate the DNA sequence divergences before and after speciation of each species pair along with their standard deviations. Using 1.6 x 10(-8) as the rate of nucleotide substitutions per silent site per year, our results indicate that the D. melanogaster lineage split from D. yakuba approximately 5.1 +/- 0.8 million years ago (mya), D. mauritiana 2.7 +/- 0.4 mya, and D. simulans 2.3 +/- 0.3 mya. It implies that D. melanogaster became distinct from D. mauritiana and D. simulans at approximately the same time and from D. yakuba no earlier than 10 mya. The effective ancestral population size of D. melanogaster appears to be stable over evolutionary time. Assuming 10 generations per year for Drosophila, the effective population size in the ancestral lineage immediately prior to the time of species divergence is approximately 3 x 10(6), which is close to that estimated for the extant D. melanogaster population. The D. melanogaster did not encounter any obvious bottleneck during the past 10 million years.     

Phylogenetic incongruence in the Drosophila melanogaster species group

Drosophila melanogaster and its close relatives are used extensively in comparative biology. Despite the importance of phylogenetic information for such studies, relationships between some melanogaster species group members are unclear due to conflicting phylogenetic signals at different loci. In this study, we use twelve nuclear loci (eleven coding and one non-coding) to assess the degree of phylogenetic incongruence in this model system. We focus on two nodes: (1) the node joining the Drosophila erecta-Drosophila orena, Drosophila melanogaster-Drosophila simulans, and Drosophila yakuba-Drosophila teissieri lineages, and (2) the node joining the lineages leading to the melanogaster, takahashii, and eugracilis subgroups. We find limited evidence for incongruence at the first node; our data, as well as those of several previous studies, strongly support monophyly of a clade consisting of D. erecta-D. orena and D. yakuba-D. teissieri. By contrast, using likelihood based tests of congruence, we find robust evidence for topological incongruence at the second node. Different loci support different relationships among the melanogaster, takahashii, and eugracilis subgroups, and the observed incongruence is not easily attributable to homoplasy, non-equilibrium base composition, or positive selection on a subset of loci. We argue that lineage sorting in the common ancestor of these three subgroups is the most plausible explanation for our observations. Such lineage sorting may lead to biased estimation of tree topology and evolutionary rates, and may confound inferences of positive selection.     

Molecular Evolution between Drosophila Melanogaster and D. Simulans: Reduced Codon Bias, Faster Rates of Amino Acid Substitution, and Larger Proteins in D. Melanogaster

Both natural selection and mutational biases contribute to variation in codon usage bias within Drosophila species. This study addresses the cause of codon bias differences between the sibling species, Drosophila melanogaster and D. simulans. Under a model of mutation-selection-drift, variation in mutational processes between species predicts greater base composition differences in neutrally evolving regions than in highly biased genes. Variation in selection intensity, however, predicts larger base composition differences in highly biased loci. Greater differences in the G+C content of 34 coding regions than 46 intron sequences between D. melanogaster and D. simulans suggest that D. melanogaster has undergone a reduction in selection intensity for codon bias. Computer simulations suggest at least a fivefold reduction in N(e)s at silent sites in this lineage. Other classes of molecular change show lineage effects between these species. Rates of amino acid substitution are higher in the D. melanogaster lineage than in D. simulans in 14 genes for which outgroup sequences are available. Surprisingly, protein sizes are larger in D. melanogaster than in D. simulans in the 34 genes compared between the two species. A substantial fraction of silent, replacement, and insertion/deletion mutations in coding regions may be weakly selected in Drosophila.     

Evidence for de novo evolution of testis-expressed genes in the Drosophila yakuba/Drosophila erecta clade

The mutational origin and subsequent evolution of de novo genes, which are hypothesized to be genes of recent origin that are not obviously related to ancestral coding sequence, are poorly understood. However, accumulating evidence suggests that such genes may often function in male reproduction. Here we use testis-derived expressed sequence tags (ESTs) from Drosophila yakuba to identify genes that have likely arisen either in D. yakuba or in the D. yakuba/D. erecta ancestor. We found several such genes, which show testis-biased expression and are often X-linked. Comparative data indicate that three of these genes have very short open reading frames, which suggests the possibility that a significant number of testis-biased de novo genes in the D. yakuba/D. erecta clade may be noncoding RNA genes. These data, along with previously published data from D. melanogaster, support the idea that many de novo Drosophila genes function in male reproduction and that a small region of the X chromosome in the melanogaster subgroup may be a hotspot for the evolution of novel testis-biased genes.     

Mitochondrial DNA evolution in themelanogaster species subgroup ofDrosophila

Detailed restriction maps (40 cleavage sites on average) of mitochondrial DNAs (mtDNAs) from the eight species of the melanogaster species subgroup of Drosophila were established. Comparison of the cleavage sites allowed us to build a phylogenetic tree based on the matrix of nucleotide distances and to select the most parsimonious network. The two methods led to similar results, which were compared with those in the literature obtained from nuclear characters. The three chromosomally homosequential species D. simulans, D. mauritiana, and D. sechellia are mitochondrially very related, but exhibit complex phylogenetic relationships. D. melanogaster is their closest relative, and the four species form a monophyletic group (the D. melanogaster complex), which is confirmed by the shared unusual length of their mt genomes (18-19 kb). The other four species of the subgroup (D. yakuba, D. teissieri, D. erecta, and D. orena) are characterized by a much shorter mt genome (16-16.5 kb). The monophyletic character of the D. yakuba complex, however, is questionable. Two species of this complex, D. yakuba and D. teissieri, are mitochondrially indistinguishable (at the level of our investigation) in spite of their noticeable allozymic and chromosomal divergence. Finally, mtDNA distances were compared with the nuclear-DNA distances thus far established. These sequences seem to evolve at rather similar rates, the mtDNA rate being barely double that of nuclear DNA.     

Nucleotide variation in the triosephosphate isomerase (Tpi) locus of Drosophila melanogaster and Drosophila simulans

DNA sequence variation in a 1.1-kb region including the coding portion of the Tpi locus was examined in 25 homozygous third-chromosome lines of Drosophila melanogaster, nine lines of Drosophila simulans, and one line of Drosophila yakuba. Our data show that the widespread allozyme polymorphism observed in cosmopolitan D. melanogaster is due to a glutamic acid substitution occurring in a phylogenetically conserved lysine that has been identified as part of the "hinged-lid" active site of the enzyme. This observation suggests that the replacement polymorphism may have important functional consequences. One replacement polymorphism was also observed in D. simulans, although its functional relevance is more difficult to assess, since it affects a site that is not strongly conserved. This amino acid change in D. simulans is associated with a single lineage possessing seven unique silent substitutions, which may be indicative of balancing selection or population subdivision. The absence of fixed amino acid differences between D. melanogaster and D. simulans and only a single difference with D. yakuba suggests that triose phosphate isomerase is under strong functional constraint. Silent variation is slightly higher for D. melanogaster than for D. simulans. Finally, we outline the general lack of evidence for old balanced polymorphisms at allozyme loci in D. melanogaster.      

Synthesis patterns of a set of follicle cell proteins in Drosophila melanogaster sibling species

The protein synthesis pattern of a set of stage and tissue specific proteins has previously been described in Drosophila melanogaster. The analysis of this set of follicle cell proteins (Fc proteins) is here extended to cover several sibling species of Drosophila melanogaster, namely D. simulans, D. mauritiana, D. erecta and D. yakuba. Even though a similar set of proteins were synthesized in these species, minor differences in size of the proteins were found between the species. Some of the species exhibited variation within species.