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

Diversity-enhancing selection acts on a female reproductive protease family in four subspecies of Drosophila mojavensis

Protein components of the Drosophila male ejaculate are critical modulators of reproductive success, several of which are known to evolve rapidly. Recent evidence of adaptive evolution in female reproductive tract proteins suggests this pattern may reflect sexual selection at the molecular level. Here we explore the evolutionary dynamics of a five-paralog gene family of female reproductive proteases within geographically isolated subspecies of Drosophila mojavensis. Remarkably, four of five paralogs show exceptionally low differentiation between subspecies and unusually structured haplotypes that suggest the retention of old polymorphisms. These gene genealogies are accompanied by deviations from neutrality consistent with diversifying selection. While diversifying selection has been observed among the reproductive molecules of mammals and marine invertebrates, our study provides the first evidence of this selective regime in any Drosophila reproductive protein, male or female.

Assigning sequences to species in the absence of large interspecific differences

Barcoding is an initiative to define a standard fragment of DNA to be used to assign unknown sequences to existing known species groups that have been pre-identified externally (by a taxonomist). Several methods have been described that attempt to place this assignment into a Bayesian statistical framework. Here we describe an algorithm that makes use of segregating sites and we examine how well these methods perform in the absence of an interspecific 'barcoding gap'. When a barcoding gap exists, that is when the data are clearly delimited, most methods perform well. Here we have used data from the Drosophila genus because this genus includes sibling species and the species relationships within this species while complex are, arguably, better understood than in any other group. The results show that the Bayesian methods perform well even in the absence of a barcoding gap. The sequences from Drosophila are correctly identified and only when the degree of incomplete lineage sorting is extreme in simulations or within the Drosophila species, do they fail in their identifications and even then, the "correct" species has a high posterior probability.

Comparative genomics of Drosophila mtDNA: Novel features of conservation and change across functional domains and lineages

To gain insight on mitochondrial DNA (mtDNA) evolution, we assembled and analyzed the mitochondrial genomes of Drosophila erecta, D. ananassae, D. persimilis, D. willistoni, D. mojavensis, D. virilis and D. grimshawi together with the sequenced mtDNAs of the melanogaster subgroup. Genomic comparisons across the well-defined Drosophila phylogeny impart power for detecting conserved mtDNA regions that maintain metabolic function and regions that evolve uniquely on lineages. Evolutionary rate varies across intergenic regions of the mtDNA. Rapidly evolving intergenic regions harbor the majority of mitochondrial indel divergence. In contrast, patterns of nearly perfect conservation within intergenic regions reveal a refined set of nucleotides underlying the binding of transcription termination factors. Sequencing of 5' cDNA ends indicates that cytochrome C oxidase I (CoI) has a novel (T/C)CG start codon and that perfectly conserved regions upstream of two NADH dehydrogenase (ND) genes are transcribed and likely extend these protein sequences. Substitutions at synonymous sites in the Drosophila mitochondrial proteomes reflect a mutation process that is biased toward A and T nucleotides and differs between mtDNA strands. Differences in codon usage bias across genes reveal that weak selection at silent sites may offset the mutation bias. The mutation-selection balance at synonymous sites has also diverged between the Drosophila and Sophophora lineages. Rates of evolution are highly heterogeneous across the mitochondrial proteome, with ND accumulating many more amino acid substitutions than CO. These oxidative phosphorylation complex-specific rates of evolution vary across lineages and may reflect physiological and ecological change across the Drosophila phylogeny.

Post-zygotic isolation in cactophilic Drosophila: larval viability and adult life-history traits of D. mojavensis/D. arizonae hybrids

Drosophila mojavensis and Drosophila arizonae are cactophilic flies that have been used extensively in speciation studies. Incomplete premating isolation, evidence of reinforcement, and a lack of recent introgression between these species point to a potentially important role for post-zygotic isolating barriers in this system. Other than hybrid male sterility, however, post-zygotic isolation between D. mojavensis and D. arizonae has received little attention. In this study, we examined viability and life-history traits of D. mojavensis/D. arizonae F(1) hybrids from sympatric crosses. Specifically, we reared hybrids and purebreds on the natural host cacti of each parental species and compared viability, development time, thorax length, and desiccation resistance between hybrids and purebreds. Interestingly, hybrid females from both crosses performed similarly or even better than purebred females. In contrast, hybrid sons of D. arizonae mothers, in addition to being sterile, had shorter average thorax length than males of both parental species, and hybrid males from both crosses had substantially lower desiccation resistance than D. mojavensis males. The probable cost to hybridization for D. mojavensis females resulting from reduced desiccation resistance of hybrid sons may have been an important selective factor in the history of reinforcement for crosses involving these females.

Multiple introductions of the Spiroplasma bacterial endosymbiont into Drosophila

Bacterial endosymbionts are common in insects and can have dramatic effects on their host's evolution. So far, the only heritable symbionts found in Drosophila have been Wolbachia and Spiroplasma. While the incidence and effects of Wolbachia have been studied extensively, the prevalence and significance of Spiroplasma infections in Drosophila are less clear. These small, gram-positive, helical bacteria infect a diverse array of plant and arthropod hosts, conferring a variety of fitness effects. Male-killing Spiroplasma are known from certain Drosophila species; however, in others, Spiroplasma appear not to affect sex ratio. Previous studies have identified different Spiroplasma haplotypes in Drosophila populations, although no extensive surveys have yet been reported. We used a multilocus sequence analysis to reconstruct a robust Spiroplasma endosymbiont phylogeny, assess genetic diversity, and look for evidence of recombination. Six loci were sequenced from over 65 Spiroplasma-infected individuals from nine different Drosophila species. Analysis of these sequences reveals at least five separate introductions of four phylogenetically distinct Spiroplasma haplotypes, indicating that more extensive sampling will likely reveal an even greater Spiroplasma endosymbiont diversity. Patterns of variation in Drosophila mitochondrial haplotypes in Spiroplasma-infected and uninfected flies imply imperfect vertical transmission in host populations and possible horizontal transmission.

Duplication, selection and gene conversion in a Drosophila mojavensis female reproductive protein family

Protein components of the Drosophila male ejaculate, several of which evolve rapidly, are critical modulators of reproductive success. Recent studies of female reproductive tract proteins indicate they also are extremely divergent between species, suggesting that reproductive molecules may coevolve between the sexes. Our current understanding of intersexual coevolution, however, is severely limited by the paucity of genetic and evolutionary studies on the female molecules involved. Physiological evidence of ejaculate-female coadaptation, paired with a promiscuous mating system, makes Drosophila mojavensis an exciting model system in which to study the evolution of reproductive proteins. Here we explore the evolutionary dynamics of a five-paralog gene family of female reproductive proteases within populations of D. mojavensis and throughout the repleta species group. We show that the proteins have experienced ongoing gene duplication and adaptive evolution and further exhibit dynamic patterns of pseudogenation, copy number variation, gene conversion, and selection within geographically isolated populations of D. mojavensis. The integration of these patterns in a single gene family has never before been documented in a reproductive protein.

Genetic diversification and demographic history of the cactophilic pseudoscorpion Dinocheirus arizonensis from the Sonoran Desert

Sequence data from a segment of the mitochondrial cytochrome c oxidase subunit I (COI) gene were used to examine phylogenetic relationships, estimate gene flow and infer demographic history of the cactophilic chernetid pseudoscorpion, Dinocheirus arizonensis (Banks), from the Sonoran Desert. Phylogenetic trees resolved two clades of D. arizonensis, one from mainland Sonora, Mexico and southern Arizona (clade I) and the other from the Baja California peninsula and southern Arizona (clade II). The two clades were separated by a mean genetic distance (d) of approximately 2.6%. Hierarchical analysis of molecular variance indicated highly significant population structuring in D. arizonensis (overall Phi(ST)=0.860; P<0.0001), with 80% of the genetic variation distributed among the two clades. Most pairwise comparisons of Phi(ST) among populations within each clade, however, were not significant. The results suggest that phoretic dispersal on vagile cactophilic insects such as the neriid cactus fly Odontoloxozus longicornis (Coquillett) provides sufficient gene flow to offset the accumulation of unique haplotypes within each clade of the non-vagile pseudoscorpion. Preliminary results on dispersal capability of O. longicornis were consistent with this conclusion. Tests designed to reconstruct demographic history from sequence data indicated that both clades of D. arizonensis, as well as O. longicornis, have experienced historical population expansions. Potential barriers to gene flow that may have led to genetic isolation and diversification in clades I and II of D. arizonensis are discussed.

Genetic architecture of hybrid male sterility in Drosophila: analysis of intraspecies variation for interspecies isolation

BACKGROUND

The genetic basis of postzygotic isolation is a central puzzle in evolutionary biology. Evolutionary forces causing hybrid sterility or inviability act on the responsible genes while they still are polymorphic, thus we have to study these traits as they arise, before isolation is complete.

METHODOLOGY/PRINCIPAL FINDINGS

Isofemale strains of D. mojavensis vary significantly in their production of sterile F(1) sons when females are crossed to D. arizonae males. We took advantage of the intraspecific polymorphism, in a novel design, to perform quantitative trait locus (QTL) mapping analyses directly on F(1) hybrid male sterility itself. We found that the genetic architecture of the polymorphism for hybrid male sterility (HMS) in the F(1) is complex, involving multiple QTL, epistasis, and cytoplasmic effects.

CONCLUSIONS/SIGNIFICANCE

The role of extensive intraspecific polymorphism, multiple QTL, and epistatic interactions in HMS in this young species pair shows that HMS is arising as a complex trait in this system. Directional selection alone would be unlikely to maintain polymorphism at multiple loci, thus we hypothesize that directional selection is unlikely to be the only evolutionary force influencing postzygotic isolation.

Sexual isolation and mating propensity among allopatric Drosophila mettleri populations

Drosophila mettleri is found in deserts of North America breeding in soil soaked by the juices of necrotic cacti. Saguaro (Carnegiea gigantea) and cardón (Pachycereus pringlei) are the usual host cacti in Mexico and Arizona, while prickly pear (Opuntia spp.) is used by an isolated population on Santa Catalina Island off the southern California Coast. Populations of D. mettleri show significant local genetic differentiation, especially when geographical isolation is coupled with host shifts. We tested for evidence of sexual isolation among allopatric populations of D. mettleri using a variety of choice and no-choice tests. Populations exhibited significant differences in mating propensity, which translated into significant deviations from random mating. While in some cases these deviations were consistent with sexual isolation, in others, negative assortative mating was observed. No relationship between degree of genetic differentiation and the appearance of sexual isolation was detected.

Molecular evolution and population genetics of two Drosophila mettleri cytochrome P450 genes involved in host plant utilization

Understanding the genetic basis of adaptation is one of the primary goals of evolutionary biology. The evolution of xenobiotic resistance in insects has proven to be an especially suitable arena for studying the genetics of adaptation, and resistant phenotypes are known to result from both coding and regulatory changes. In this study, we examine the evolutionary history and population genetics of two Drosophila mettleri cytochrome P450 genes that are putatively involved in the detoxification of alkaloids present in two of its cactus hosts: saguaro (Carnegiea gigantea) and senita (Lophocereus schottii). Previous studies demonstrated that Cyp28A1 was highly up-regulated following exposure to rotting senita tissue while Cyp4D10 was highly up-regulated following exposure to rotting saguaro tissue. Here, we show that a subset of sites in Cyp28A1 experienced adaptive evolution specifically in the D. mettleri lineage. Moreover, neutrality tests in several populations were also consistent with a history of selection on Cyp28A1. In contrast, we did not find evidence for positive selection on Cyp4D10, although this certainly does not preclude its involvement in host plant use. A surprising result that emerged from our population genetic analyses was the presence of significant genetic differentiation between flies collected from different host plant species (saguaro and senita) at Organ Pipe National Monument, Arizona, USA. This preliminary evidence suggests that D. mettleri may have evolved into distinctive host races that specialize on different hosts, a possibility that warrants further investigation.

Species-specific shifts in centromere sequence composition are coincident with breakpoint reuse in karyotypically divergent lineages

The evolution of three classes of centromere sequences across nine species of macropodine marsupials were compared with that of other genes, showing that each species has experienced differential expansion and contraction of individual classes.

Background

It has been hypothesized that rapid divergence in centromere sequences accompanies rapid karyotypic change during speciation. However, the reuse of breakpoints coincident with centromeres in the evolution of divergent karyotypes poses a potential paradox. In distantly related species where the same centromere breakpoints are used in the independent derivation of karyotypes, centromere-specific sequences may undergo convergent evolution rather than rapid sequence divergence. To determine whether centromere sequence composition follows the phylogenetic history of species evolution or patterns of convergent breakpoint reuse through chromosome evolution, we examined the phylogenetic trajectory of centromere sequences within a group of karyotypically diverse mammals, macropodine marsupials (wallabies, wallaroos and kangaroos).

Results

The evolution of three classes of centromere sequences across nine species within the genus Macropus (including Wallabia) were compared with the phylogenetic history of a mitochondrial gene, Cytochrome b (Cyt b), a nuclear gene, selenocysteine tRNA (TRSP), and the chromosomal histories of the syntenic blocks that define the different karyotype arrangements. Convergent contraction or expansion of predominant satellites is found to accompany specific karyotype rearrangements. The phylogenetic history of these centromere sequences includes the convergence of centromere composition in divergent species through convergent breakpoint reuse between syntenic blocks.

Conclusion

These data support the 'library hypothesis' of centromere evolution within this genus as each species possesses all three satellites yet each species has experienced differential expansion and contraction of individual classes. Thus, we have identified a correlation between the evolution of centromere satellite sequences, the reuse of syntenic breakpoints, and karyotype convergence in the context of a gene-based phylogeny.

The molecular basis of host adaptation in cactophilic Drosophila: molecular evolution of a glutathione S-transferase gene (GstD1) in Drosophila mojavensis

Drosophila mojavensis is a cactophilic fly endemic to the northwestern deserts of North America. This species includes four genetically isolated cactus host races each individually specializing on the necrotic tissues of a different cactus species. The necrosis of each cactus species provides the resident D. mojavensis populations with a distinct chemical environment. A previous investigation of the role of transcriptional variation in the adaptation of D. mojavensis to its hosts produced a set of candidate loci that are differentially expressed in response to host shifts, and among them was glutathione S-transferase D1 (GstD1). In both D. melanogaster and Anopheles gambiae, GstD1 has been implicated in the resistance of these species to the insecticide dichloro-diphenyl-trichloroethane (DDT). The pattern of sequence variation of the GstD1 locus from all four D. mojavensis populations, D. arizonae (sister species), and D. navojoa (outgroup) has been examined. The data suggest that in two populations of D. mojavensis GstD1 has gone through a period of adaptive amino acid evolution. Further analyses indicate that of the seven amino acid fixations that occurred in the D. mojavensis lineage, two of them occur in the active site pocket, potentially having a significant effect on substrate specificity and in the adaptation to alternative cactus hosts.

Multilocus nuclear sequences reveal intra- and interspecific relationships among chromosomally polymorphic species of cactophilic Drosophila

Drosophila mojavensis and Drosophila arizonae, a pair of sibling species endemic to North America, constitute an important model system to study ecological genetics and the evolution of reproductive isolation. This species pair can produce fertile hybrids in some crosses and are sympatric in a large part of their ranges. Despite the potential for hybridization in nature, however, evidence of introgression has not been rigorously sought. Further, the evolutionary relationships within and among the geographically distant populations of the two species have not been characterized in detail using high-resolution molecular studies. Both species have six chromosomes: five large acrocentrics and one 'dot' chromosome. Fixed inversion differences between the species exist in three chromosomes (X, 2 and 3) while three are colinear (4, 5 and 6), suggesting that were introgression to occur, it would be most likely in the colinear chromosomes. We utilized nucleotide sequence variation at multiple loci on five chromosomes to test for evidence of introgression, and to test various scenarios for the evolutionary relationships of these two species and their populations. While we do not find evidence of recent introgression, loci in the colinear chromosomes appear to have participated in exchange in the past. We also found considerable population structure within both species. The level of differentiation discovered among D. arizonae populations was unexpectedly high and suggests that its populations, as well as those of D. mojavensis, may be themselves undergoing incipient speciation and merit further attention.

Genetic differentiation and demographic history in Drosophila pachea from the Sonoran Desert

Genetic variation at six microsatellite DNA loci and a segment of the mitochondrial cytochrome oxidase subunit I (COI) locus was used to estimate gene flow, population structure, and demographic history in the cactophilic Drosophila pachea from the Sonoran Desert of North America, a species that shows a strict association with its senita host cactus (genus Lophocereus). For microsatellite analyses, thirteen populations of D. pachea were sampled, five in mainland Mexico and the southwestern USA, and eight on the Baja California (Baja) peninsula, covering essentially the entire range of the species. Analysis of molecular variance (AMOVA) of microsatellite data revealed that populations from both the mainland and the Baja peninsula generally showed little structure, although there were a few exceptions, suggesting some local differentiation and restriction of gene flow within both regions. Pairwise comparisons of F(ST) among each of the mainland and Baja populations showed evidence of both panmixia and population subdivision. AMOVA performed on grouped populations from both the mainland and Baja, however, revealed significant partitioning of genetic variation among the two regions, but no partitioning among localities within each region. Bayesian skyline analyses of the COI data set, consisting of four mainland and seven peninsular populations, revealed population expansions dating to the Pleistocene or late Pliocene in D. pachea from both regions, although regional differences were seen in the estimated timing of the expansions and in changes in effective population size over time.

Functional genomics of cactus host shifts in Drosophila mojavensis

Understanding the genetic basis of adaptation to novel environments remains one of the major challenges confronting evolutionary biologists. While newly developed genomic approaches hold considerable promise for addressing this overall question, the relevant tools have not often been available in the most ecologically interesting organisms. Our study organism, Drosophila mojavensis, is a cactophilic Sonoran Desert endemic utilizing four different cactus hosts across its geographical range. Its well-known ecology makes it an attractive system in which to study the evolution of gene expression during adaptation. As a cactophile, D. mojavensis oviposits in the necrotic tissues of cacti, therefore exposing larvae and even adults to the varied and toxic compounds of rotting cacti. We have developed a cDNA microarray of D. mojavensis to examine gene expression associated with cactus host use. Using a population from the Baja California population we examined gene expression differences of third instar larvae when reared in two chemically distinct cactus hosts, agria (Stenocereus gummosus, native host) vs. organpipe (Stenocereus thurberi, alternative host). We have observed differential gene expression associated with cactus host use in genes involved in metabolism and detoxification.