Adaptive evolution of recently duplicated accessory gland protein genes in desert Drosophila

The relationship between animal mating system variation and patterns of protein polymorphism and divergence is poorly understood. Drosophila provides an excellent system for addressing this issue, as there is abundant interspecific mating system variation. For example, compared to D. melanogaster subgroup species, repleta group species have higher remating rates, delayed sexual maturity, and several other interesting differences. We previously showed that accessory gland protein genes (Acp's) of Drosophila mojavensis and D. arizonae evolve more rapidly than Acp's in the D. melanogaster subgroup and that adaptive Acp protein evolution is likely more common in D. mojavensis/D. arizonae than in D. melanogaster/D. simulans. These findings are consistent with the idea that greater postcopulatory selection results in more adaptive evolution of seminal fluid proteins in the repleta group flies. Here we report another interesting evolutionary difference between the repleta group and the D. melanogaster subgroup Acp's. Acp gene duplications are present in D. melanogaster, but their high sequence divergence indicates that the fixation rate of duplicated Acp's has been low in this lineage. Here we report that D. mojavensis and D. arizonae genomes contain several very young duplicated Acp's and that these Acp's have experienced very rapid, adaptive protein divergence. We propose that rapid remating of female desert Drosophila generates selection for continuous diversification of the male Acp complement to improve male fertilization potential. Thus, mating system variation may be associated with adaptive protein divergence as well as with duplication of Acp's in Drosophila.     

Positive selection and the molecular evolution of a gene of male reproduction, Acp26Aa of Drosophila

The gene for a male ejaculatory protein, Acp26Aa, in four sibling species of the Drosophila melanogaster subgroup has previously been shown to have a nonsynonymous rate (Ka) of nucleotide substitution that is indistinguishable from the synonymous rate (Ks). By examining this gene in two other species of this subgroup, we found that Ka is generally large and can sometimes be more than twice as large as Ks. This suggests that positive selection may be operating at this locus of male reproduction.      

Multilocus test for introgression between the cactophilic species Drosophila mojavensis and Drosophila arizonae

Information obtained from laboratory studies regarding the efficacy of barriers to gene flow (reproductive isolation) between species is often incomplete or misleading, so detailed genetic analyses are needed to determine whether hybridization and introgression occur in nature. Previous laboratory studies of the cactophilic species Drosophila mojavensis and Drosophila arizonae suggest that reproductive isolation is incomplete and that gene flow may occur in sympatry. We sampled 18 nuclear and one mitochondrial loci from multiple populations of D. arizonae and D. mojavensis to test for the signature of recent or historic gene flow between these two species. We located chromosomal regions that were inverted between these species and analyzed those regions independently of others. Statistical tests for introgression using all loci or only collinear loci failed to reject expectations of an isolation model. Further tests using average nucleotide differences between species and phylogenetic analyses also failed to find support for introgression between D. mojavensis and D. arizonae. Additional ecological and behavioral studies of these species in their natural habitats are required to explain why the signature of gene flow was not detected at the DNA sequence level in populations when laboratory studies suggest such gene flow should be possible.     

Evolutionary relationships of Drosophila mojavensis geographic host races and their sister species Drosophila arizonae

The cactophilic Drosophila mojavensis species group living in the deserts and dry tropical forests of the southwestern United States and Mexico provides a valuable system for studies in diversification and speciation. Rigorous studies of the relationships between host races of D. mojavensis and the relationships among the members of the species group (D. mojavensis, Drosophila arizona, and Drosophila navojoa) are lacking. We used mitochondrial CO1 sequence data to address the phylogenetics and population genetics of this species group. In this study we have found that the sister species D. mojavensis and D. arizonae share no mitochondrial haplotypes and thus show no evidence for recent introgression. We estimate the divergence time between D. mojavensis and D. arizonae to be between 1.91 and 2.97 million years ago. D. arizonae shows little structure in our population genetic analyses but there is phylogenetic differentiation between southeastern and northern populations of D. arizonae. Drosophila mojavensis shows significant population and phylogenetic structure across the four geographic regions of its distribution. The mitochondrial data support an origin of D. mojavensis on the mainland with early differentiation into the populations now found in the Mojave Desert and the Mainland Sonoran Desert and later colonization of the Baja Peninsula, in contrast to previous models. Also, the sister clade to D. mojavensis/D. arizonae includes D. navojoa and Drosophila huaylasi. By defining the genetic relationships among these populations, we provide a foundation for more sophisticated hypothesis testing regarding the timing of early speciation events and host switches in this species group.     

Sequence variation of alcohol dehydrogenase (Adh) paralogs in cactophilic Drosophila

This study focuses on the population genetics of alcohol dehydrogenase (Adh) in cactophilic Drosophila. Drosophila mojavensis and D. arizonae utilize cactus hosts, and each host contains a characteristic mixture of alcohol compounds. In these Drosophila species there are two functional Adh loci, an adult form (Adh-2) and a larval and ovarian form (Adh-1). Overall, the greater level of variation segregating in D. arizonae than in D. mojavensis suggests a larger population size for D. arizonae. There are markedly different patterns of variation between the paralogs across both species. A 16-bp intron haplotype segregates in both species at Adh-2, apparently the product of an ancient gene conversion event between the paralogs, which suggests that there is selection for the maintenance of the intron structure possibly for the maintenance of pre-mRNA structure. We observe a pattern of variation consistent with adaptive protein evolution in the D. mojavensis lineage at Adh-1, suggesting that the cactus host shift that occurred in the divergence of D. mojavensis from D. arizonae had an effect on the evolution of the larval expressed paralog. Contrary to previous work we estimate a recent time for both the divergence of D. mojavensis and D. arizonae (2.4 +/- 0.7 MY) and the age of the gene duplication (3.95 +/- 0.45 MY).     

Sympatry, allopatry and sexual isolation between Drosophila mojavensis and D. arizonae

Populations of the North American cactophilic fruitfly Drosophila mojavensis and its sibling species D. arizonae exist both in sympatry and in allopatry. Females of D. arizonae, regardless of their population of origin, are effectively completely isolated behaviorally from D. mojavensis males. On the other hand, females of D. mojavensis from the sympatric populations in Sonora, Mexico exhibit significantly stronger premating isolation from D. arizonae males than do D. mojavensis females from allopatric populations from the Baja California peninsula. Earlier studies interpreted these limited observations as support for reinforcement. Since the time of those studies, additional allopatric populations of D. mojavensis have been collected from southern California and from Santa Catalina Island, off the coast of southern California. Here, we tested the prediction that if sympatry is in fact associated with increased isolation in D. mojavensis, these additional allopatric populations also should show, relative to the sympatric ones, less isolation from D. arizonae. Our results are consistent with this prediction and suggest that isolation is in fact stronger in sympatry.     

An Autosomal Factor from Drosophila Arizonae Restores Normal Spermatogenesis in Drosophila Mojavensis Males Carrying the D. Arizonae Y Chromosome

Males of Drosophila mojavensis whose Y chromosome is replaced by the Y chromosome of the sibling species Drosophila arizonae are sterile. It is shown that genetic material from the fourth chromosome of D. arizonae is necessary and sufficient, in single dose, to restore fertility in these males. In introgression and mapping experiments this material segregates as a single Mendelian factor (sperm motility factor, SMF). Light and electron microscopy studies of spermatogenesis in D. mojavensis males whose Y chromosome is replaced by introgression with the Y chromosome of D. arizonae (these males are symbolized as mojY(a)) revealed postmeiotic abnormalities all of which are restored when the SMF of D. arizonae is co-introgressed (these males are symbolized as mojY(a)SMF(a)). The number of mature sperm per bundle in mojY(a)SMF(a) is slightly less than in pure D. mojavensis and is even smaller in males whose fertility is rescued by introgression of the entire fourth chromosome of D. arizonae. These observations establish an interspecific incompatibility between the Y chromosome and an autosomal factor (or more than one tightly linked factors) that can be useful for the study of the evolution of male hybrid sterility in Drosophila and the genetic control of spermatogenesis.     

Activity variation in alcohol dehydrogenase paralogs is associated with adaptation to cactus host use in cactophilic Drosophila

Drosophila mojavensis and Drosophila arizonae are species of cactophilic flies that share a recent duplication of the alcohol dehydrogenase (Adh) locus. One paralog (Adh-2) is expressed in adult tissues and the other (Adh-1) in larvae and ovaries. Enzyme activity measurements of the ADH-2 amino acid polymorphism in D. mojavensis suggest that the Fast allozyme allele has a higher activity on 2-propanol than 1-propanol. The Fast allele was found at highest frequency in populations that utilize hosts with high proportions of 2-propanol, while the Slow allele is most frequent in populations that utilize hosts with high proportions of 1-propanol. This suggests that selection for ADH-2 allozyme alleles with higher activity on the most abundant alcohols is occurring in each D. mojavensis population. In the other paralog, ADH-1, significant differences between D. mojavensis and D. arizonae are associated with a previously shown pattern of adaptive protein evolution in D. mojavensis. Examination of protein sequences showed that a large number of amino acid fixations between the paralogs have occurred in catalytic residues. These changes are potentially responsible for the significant difference in substrate specificity between the paralogs. Both functional and sequence variation within and between paralogs suggests that Adh has played an important role in the adaptation of D. mojavensis and D. arizonae to their cactophilic life.     

Hierarchical structure in the Drosophila mojavensis cluster (Diptera: Drosophilidae)

Sequences from three gene regions from the nuclear and mitochondrial genomes were used to examine the extent and nature of hierarchical structure in the Drosophilamojavensis cluster (Drosophila arizonae, D. mojavensis and D. navojoa) of the D. repleta species group. To determine the genetic divergence of these three species, sequence data were analyzed using maximum parsimony and population aggregation analysis. Individual and combined gene genealogies indicate that D. arizonae and D. mojavensis are neither diagnosable nor monophyletic with respect to one another. Although D. navojoa has differentiated from D. arizonae and D. mojavensis, as diagnosed by nuclear gene sequences, it may have undergone a reticulation event with D. arizonae. Our results suggest that either these taxa are still undergoing differentiation at the molecular level or have experienced gene flow in the recent past.     

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.     

Molecular evolution of sex-biased genes in Drosophila

Studies of morphology, interspecific hybridization, protein/DNA sequences, and levels of gene expression have suggested that sex-related characters (particularly those involved in male reproduction) evolve rapidly relative to non-sex-related characters. Here we report a general comparison of evolutionary rates of sex-biased genes using data from cDNA microarray experiments and comparative genomic studies of Drosophila. Comparisons of nonsynonymous/synonymous substitution rates (d(N)/d(S)) between species of the D. melanogaster subgroup revealed that genes with male-biased expression had significantly faster rates of evolution than genes with female-biased or unbiased expression. The difference was caused primarily by a higher d(N) in the male-biased genes. The same pattern was observed for comparisons among more distantly related species. In comparisons between D. melanogaster and D. pseudoobscura, genes with highly biased male expression were significantly more divergent than genes with highly biased female expression. In many cases, orthologs of D. melanogaster male-biased genes could not be identified in D. pseudoobscura through a Blast search. In contrast to the male-biased genes, there was no clear evidence for accelerated rates of evolution in female-biased genes, and most comparisons indicated a reduced rate of evolution in female-biased genes relative to unbiased genes. Male-biased genes did not show an increased ratio of nonsynonymous/synonymous polymorphism within D. melanogaster, and comparisons of polymorphism/divergence ratios suggest that the rapid evolution of male-biased genes is caused by positive selection.     

Patterns of divergence among conifer ESTs and polymorphism in Pinus sylvestris identify putative selective sweeps

Finding genes that are under positive selection is a difficult task, especially in non-model organisms. Here, we have analyzed expressed sequence tag (EST) data from 4 species (Pinus pinaster, Pinus taeda, Picea glauca, and Pseudotsuga menziesii) to investigate selection patterns during their evolution and to identify genes likely to be under positive selection. To confirm selection, population samples of these genes have been sequenced in Pinus sylvestris, a species that was not included in the EST data set. The estimates of branch-specific Ka/Ks (nonsynonymous/synonymous substitution rates) across all genes in the EST data set were similar or smaller than estimates from other higher plant species. There was no evidence for the traditional indication of positive selection, Ka/Ks above 1. However, several lines of evidence based on polymorphism patterns suggest that genes with high Ka/Ks (0.20-0.52) in the EST data set are in fact more affected by positive selection in P. sylvestris than genes with low Ka/Ks (0.01-0.04). The high Ka/Ks genes have a lower level of polymorphism and more negative Tajima's D than the low Ka/Ks genes. Further, in the high Ka/Ks group, the Hudson-Kreitman-Aguade test is significant. This suggests that the EST data set is a good starting point for finding genes under positive selection in conifers and that even moderate Ka/Ks values could be indicative of selection. A group of 5 genes with high Ka/Ks collectively show evidence for positive selection within P. sylvestris.     

Localization and characterization of X chromosome inversion breakpoints separating Drosophila mojavensis and Drosophila arizonae

Ectopic exchange between transposable elements or other repetitive sequences along a chromosome can produce chromosomal inversions. As a result, genome sequence studies typically find sequence similarity between corresponding inversion breakpoint regions. Here, we identify and investigate the breakpoint regions of the X chromosome inversion distinguishing Drosophila mojavensis and Drosophila arizonae. We localize one inversion breakpoint to 13.7 kb and localize the other to a 1-Mb interval. Using this localization and assuming microsynteny between Drosophila melanogaster and D. arizonae, we pinpoint likely positions of the inversion breakpoints to windows of less than 3000 bp. These breakpoints define the size of the inversion to approximately 11 Mb. However, in contrast to many other studies, we fail to find significant sequence similarity between the 2 breakpoint regions. The localization of these inversion breakpoints will facilitate future genetic and molecular evolutionary studies in this species group, an emerging model system for ecological genetics.     

The hemoglobin genes of Drosophila

We recently reported the unprecedented occurrence of a hemoglobin gene (glob1) in the fruitfly Drosophila melanogaster. Here we investigate the structure and evolution of the glob1 gene in other Drosophila species. We cloned and sequenced glob1 genes and cDNA from D. pseudoobscura and D. virilis, and identified the glob1 gene sequences of D. simulans, D. yakuba, D. erecta, D. ananassae, D. mojavensis and D. grimshawi in the databases. Gene structure (introns in helix positions D7.0 and G7.0), gene synteny and sequence of glob1 are highly conserved, with high ds/dn ratios indicating strong purifying selection. The data suggest an important role of the glob1 protein in Drosophila, which may be the control of oxygen flow from the tracheal system. Furthermore, we identified two additional globin genes (glob2 and glob3) in the Drosophilidae. Although the sequences are highly derived, the amino acids required for heme- and oxygen-binding are conserved. In contrast to other known insect globin, the glob2 and glob3 genes harbour both globin-typical introns at positions B12.2 and G7.0. Both genes are conserved in various drosophilid species, but only expression of glob2 could be demonstrated by western blotting and RT-PCR. Phylogenetic analyses show that the clade leading to glob2 and glob3, which are sistergroups, diverged first in the evolution of dipteran globins. glob1 is closely related to the intracellular hemoglobin of the botfly Gasterophilus intestinalis, and the extracellular hemoglobins from the chironomid midges derive from this clade.