Maintenance of a male-killing Wolbachia in Drosophila innubila by male-killing dependent and male-killing independent mechanisms

Many maternally inherited endosymbionts manipulate their host's reproduction in various ways to enhance their own fitness. One such mechanism is male killing (MK), in which sons of infected mothers are killed by the endosymbiont during development. Several hypotheses have been proposed to explain the advantages of MK, including resource reallocation from sons to daughters of infected females, avoidance of inbreeding by infected females, and, if transmission is not purely maternal, the facilitation of horizontal transmission to uninfected females. We tested these hypotheses in Drosophila innubila, a mycophagous species infected with MK Wolbachia. There was no evidence of horizontal transmission in the wild and no evidence Wolbachia reduced levels of inbreeding. Resource reallocation does appear to be operative, as Wolbachia-infected females are slightly larger, on average, than uninfected females, although the selective advantage of larger size is insufficient to account for the frequency of infection in natural populations. Wolbachia-infected females from the wild-although not those from the laboratory-were more fecund than uninfected females. Experimental studies revealed that Wolbachia can boost the fecundity of nutrient-deprived flies and reduce the adverse effect of RNA virus infection. Thus, this MK endosymbiont can provide direct, MK-independent fitness benefits to infected female hosts in addition to possible benefits mediated via MK.     

Expression and modulation of embryonic male-killing in Drosophila innubila: opportunities for multilevel selection

Organisms and the symbionts they harbor may experience opposing forces of selection. In particular, the contrasting inheritance patterns of maternally transmitted symbionts and their host's nuclear genes can engender conflict among organizational levels over the optimal host offspring sex ratio. This study uses a male-killing Wolbachia endosymbiont and its host Drosophila innubila to experimentally address the potential for multilevel selection in a host-symbiont system. We show that bacterial density can vary among infected females, and that females with a higher density have a more female-biased offspring sex ratio. Furthermore, bacterial density is an epigenetic and heritable trait: females with a low bacterial load have daughters with a lower-than-average bacterial density, whose offspring then experience less severe male-killing. For infected sons, the probability of embryonic mortality increases with the bacterial density in their mothers. The frequency distribution of Wolbachia density among individual D. innubila females, and therefore the dynamics of infection within populations of these flies, results both from processes affecting the growth and regulation of bacterial populations within cytoplasmic lineages and from selection among cytoplasmic lineages that vary in bacterial density. Estimates of effective population size of Wolbachia within cytoplasmic lineages and of D. innubila at the host population level suggest that selection among cytoplasmic lineages is likely to overwhelm the results of selection within lineages.     

The effective founder effect in a spatially expanding population

The gradual loss of diversity and the establishment of clines in allele frequencies associated with range expansions are patterns observed in many species, including humans. These patterns can result from a series of founder events occurring as populations colonize previously unoccupied areas. We develop a model of an expanding population and, using a branching process approximation, show that spatial gradients reflect different amounts of genetic drift experienced by different subpopulations. We then use this model to measure the net average strength of the founder effect, and we demonstrate that the predictions from the branching process model fit simulation results well. We further show that estimates of the effective founder size are robust to potential confounding factors such as migration between subpopulations. We apply our method to data from Arabidopsis thaliana. We find that the average founder effect is approximately three times larger in the Americas than in Europe, possibly indicating that a more recent, rapid expansion occurred.     

Population structure leads to male‐biased population sex ratios under environmental sex determination

Spatial structure has been shown to favor female‐biased sex allocation, but current theory fails to explain male biases seen in many taxa, particularly those with environmental sex determination (ESD). We present a theory and accompanying individual‐based simulation model that demonstrates how population structure leads to male‐biased population sex ratios under ESD. Our simulations agree with earlier work showing that the high productivity of female‐producing habitats creates a net influx of sex‐determining alleles into male‐producing habitats, causing larger sex ratio biases, and lower productivity in male‐producing environments (Harts et al. 2014). In contrast to previous findings, we show that male‐biasing habitats disproportionately impact the global sex ratio, resulting in stable male‐biased population sex ratios under ESD. The failure to detect a male bias in earlier work can be attributed to small subpopulation sizes leading to local mate competition, a condition unlikely to be met in most ESD systems. Simulations revealed that consistent male biases are expected over a wide range of population structures, environmental conditions, and genetic architectures of sex determination, with male excesses as large as 30 percent under some conditions. Given the ubiquity of genetic structure in natural populations, we predict that modest, enduring male biased allocation should be common in ESD species, a pattern consistent with reviews of ESD sex ratios.     

Tissue-specific infection dynamics of male-killing and nonmale-killing spiroplasmas in Drosophila melanogaster

The male-killing spiroplasma strain NSRO causes an extremely female-biased sex ratio of the host, Drosophila melanogaster, as a result of selective death of male offspring during embryogenesis. The spiroplasma strain NSRO-A, a variant of NSRO, does not cause such symptoms. In an attempt to gain insights into the mechanism underlying the symbiont-induced reproductive phenotype, infection densities of the spiroplasmas in different tissues were monitored during host aging using a quantitative PCR technique. The density dynamics in the hemolymph were reminiscent of those in the whole body, whereas the density dynamics in the fat body, intestine and ovary were not. These results suggest that the majority of the spiroplasmas colonize and proliferate in the hemolymph of the host. In the hemolymph and whole body, the infection densities of NSRO were generally higher than those of NSRO-A, which may be related to the different reproductive phenotypes caused by the spiroplasmas.     

两性具有不同出生率和死亡率的种群动态

以往的种群动力学模型均隐含假设性比为1：1,而实际上并非总如此,不同性别的出生和死亡是不完全相同的。文中就考虑了两性具有不同出生率和死亡率的种群动态问题。可以知道种群动态只是受雌性控制,与雄性无关。     

Population and evolutionary dynamics in spatially structured seasonally varying environments

Increasingly imperative objectives in ecology are to understand and forecast population dynamic and evolutionary responses to seasonal environmental variation and change. Such population and evolutionary dynamics result from immediate and lagged responses of all key life‐history traits, and resulting demographic rates that affect population growth rate, to seasonal environmental conditions and population density. However, existing population dynamic and eco‐evolutionary theory and models have not yet fully encompassed within‐individual and among‐individual variation, covariation, structure and heterogeneity, and ongoing evolution, in a critical life‐history trait that allows individuals to respond to seasonal environmental conditions: seasonal migration. Meanwhile, empirical studies aided by new animal‐tracking technologies are increasingly demonstrating substantial within‐population variation in the occurrence and form of migration versus year‐round residence, generating diverse forms of ‘partial migration’ spanning diverse species, habitats and spatial scales. Such partially migratory systems form a continuum between the extreme scenarios of full migration and full year‐round residence, and are commonplace in nature. Here, we first review basic scenarios of partial migration and associated models designed to identify conditions that facilitate the maintenance of migratory polymorphism. We highlight that such models have been fundamental to the development of partial migration theory, but are spatially and demographically simplistic compared to the rich bodies of population dynamic theory and models that consider spatially structured populations with dispersal but no migration, or consider populations experiencing strong seasonality and full obligate migration. Second, to provide an overarching conceptual framework for spatio‐temporal population dynamics, we define a ‘partially migratory meta‐population’ system as a spatially structured set of locations that can be occupied by different sets of resident and migrant individuals in different seasons, and where locations that can support reproduction can also be linked by dispersal. We outline key forms of within‐individual and among‐individual variation and structure in migration that could arise within such systems and interact with variation in individual survival, reproduction and dispersal to create complex population dynamics and evolutionary responses across locations, seasons, years and generations. Third, we review approaches by which population dynamic and eco‐evolutionary models could be developed to test hypotheses regarding the dynamics and persistence of partially migratory meta‐populations given diverse forms of seasonal environmental variation and change, and to forecast system‐specific dynamics. To demonstrate one such approach, we use an evolutionary individual‐based model to illustrate that multiple forms of partial migration can readily co‐exist in a simple spatially structured landscape. Finally, we summarise recent empirical studies that demonstrate key components of demographic structure in partial migration, and demonstrate diverse associations with reproduction and survival. We thereby identify key theoretical and empirical knowledge gaps that remain, and consider multiple complementary approaches by which these gaps can be filled in order to elucidate population dynamic and eco‐evolutionary responses to spatio‐temporal seasonal environmental variation and change.     

Silence of the killers: discovery of male-killing suppression in a rearing strain of the small brown planthopper,Laodelphax striatellus

According to evolutionary theory, sex ratio distortions caused by reproductive parasites such as Wolbachia and Spiroplasma are predicted to be rapidly normalized by the emergence of host nuclear suppressors. However, such processes in the evolutionary arms race are difficult to observe because sex ratio biases will be promptly hidden and become superficially unrecognizable. The evolution of genetic suppressors has been reported in just two insect species so far. In the small brown planthopper, Laodelphax striatellus, female-biases caused by Spiroplasma, which is a ‘late’ male-killer, have been found in some populations. During the continuous rearing of L. striatellus, we noted that a rearing strain had a 1 : 1 sex ratio even though it harboured Spiroplasma. Through introgression crossing experiments with a strain lacking suppressors, we revealed that the L. striatellus strain had the zygotic male-killing suppressor acting as a dominant trait. The male-killing phenotype was hidden by the suppressor even though Spiroplasma retained its male-killing ability. This is the first study to demonstrate the existence of a late male-killing suppressor and its mode of inheritance. Our results, together with those of previous studies, suggest that the inheritance modes of male-killing suppressors are similar regardless of insect order or early or late male killing.     

Natural variation of the Y chromosome suppresses sex ratio distortion and modulates testis-specific gene expression in Drosophila simulans

X-linked sex-ratio distorters that disrupt spermatogenesis can cause a deficiency in functional Y-bearing sperm and a female-biased sex ratio. Y-linked modifiers that restore a normal sex ratio might be abundant and favored when a X-linked distorter is present. Here we investigated natural variation of Y-linked suppressors of sex-ratio in the Winters systems and the ability of these chromosomes to modulate gene expression in Drosophila simulans. Seventy-eight Y chromosomes of worldwide origin were assayed for their resistance to the X-linked sex-ratio distorter gene Dox. Y chromosome diversity caused males to sire ∼63% to ∼98% female progeny. Genome-wide gene expression analysis revealed hundreds of genes differentially expressed between isogenic males with sensitive (high sex ratio) and resistant (low sex ratio) Y chromosomes from the same population. Although the expression of about 75% of all testis-specific genes remained unchanged across Y chromosomes, a subset of post-meiotic genes was upregulated by resistant Y chromosomes. Conversely, a set of accessory gland-specific genes and mitochondrial genes were downregulated in males with resistant Y chromosomes. The D. simulans Y chromosome also modulated gene expression in XXY females in which the Y-linked protein-coding genes are not transcribed. The data suggest that the Y chromosome might exert its regulatory functions through epigenetic mechanisms that do not require the expression of protein-coding genes. The gene network that modulates sex ratio distortion by the Y chromosome is poorly understood, other than that it might include interactions with mitochondria and enriched for genes expressed in post-meiotic stages of spermatogenesis.     

Evolutionary dynamics and population biology of a polymorphic insect

Conspicuous heritable polymorphisms are useful to address the question if morph frequencies are stable or whether they fluctuate between generations. Ecological geneticists have studied colour polymorphisms in the past, but there are few long-term studies of genetic dynamics across multiple generations. We studied morph-frequency dynamics and female fecundity in the trimorphic blue-tailed damselfly (Ischnura elegans). The morphs include a male-coloured (androchrome) type of female, which is thought to be maintained by frequency-dependent sexual conflict. Morph frequencies changed significantly between years across all populations. There was evidence for directional frequency change since androchrome females increased in 9 of 10 populations across a 4-year period. There was heterogeneity between populations in their evolutionary trajectories, partly caused by population age: androchrome frequencies were initially high in young populations but gradually decreased and approached the level of old populations. We discuss the possible causes of morph-frequency fluctuations, and the role of morph-specific fecundity, dispersal and other forces influencing evolutionary dynamics in this system.     

Differential impacts of habitat heterogeneity on male and female connectivity in a spatially structured pest system

In a previous study, a model of landscape heterogeneity was developed and applied to a spatially structured wild rabbit (Oryctolagus cuniculus) population. That study showed clearly the influence of resource heterogeneity on connectivity levels. The simulation study was based on female movements and used population genetic validation data appropriate for a female study. Most models assume that males and females will exhibit similar patterns, although this has rarely been tested. In the current study we extend the analysis to consider differences between female and male connectivity in the same spatially structured pest system. Amplified fragment length polymorphism (AFLP) markers were screened on the same samples used previously for mtDNA analysis. The mtDNA data were used to validate female results, and AFLP data were used to validate combined male and female results. Connectivity patterns from the two simulations (female, and combined male and female) connectivity patterns showed no association. However, each was concordant with appropriate validation data, showing highly significant associations between pairwise population connectivity and the genetic data. A relative connectivity metric for the combined simulation was regressed against the mean of pairwise Φ_ST values, with almost 70% of the variation explained by a linear model. Demonstrating differential effects of habitat heterogeneity on male and female connectivity provides further evidence that spatial resource heterogeneity impacts on connectivity. Understanding differences in population connectivity will allow improved predictions of disease spread, local extinctions and recolonizations. Furthermore, modelling such differences in pest systems will allow management plans to be better targeted, for example by strategically introducing diseases for control purposes into populations which exhibit high male connectivity to aid spread, but low female connectivity to inhibit recolonization potential after control.     

Mitochondrial DNA diversity, population structure, and gender association in the gynodioecious plant Silene vulgaris

A highly variable mitochondrial DNA (mtDNA) restriction fragment length polymorphism (RFLP) locus is used to assess the population structure of mitochondrial genomes in the gynodioecious plant Silene vulgaris at two spatial scales. Thirteen mtDNA haplotypes were identified within 250 individuals from 18 populations in a 20-km diameter region of western Virginia. The population structure of these mtDNA haplotypes was estimated as thetaST = 0.574 (+/- 0.066 SE) and, surprisingly, genetic differentiation among populations was negatively correlated with geographic distance (Mantel r = -0.246, P < 0.002). Additionally, mtDNA haplotypes were spatially clumped at the scale of meters within one population. Gender in S. vulgaris is determined by an interaction between autosomal male fertility restorers and cytoplasmic male sterility (CMS) factors, and seed fitness is affected by an interaction between gender and population sex ratio; thus, selection acting on gender could influence the distribution of mtDNA RFLP haplotypes. The sex ratio (females:hermaphrodites) varied among mtDNA haplotypes across the entire metapopulation, possibly because the haplotypes were in linkage disequilibrium with different CMS factors. The gender associated with some of the most common haplotypes varied among populations, suggesting that there is also population structure in male fertility restorer genes. In comparison with reports of mtDNA variation from other published studies, we found that S. vulgaris exhibits a large number of mtDNA haplotypes relative to that observed in other species.     

Low genetic diversity and strong but shallow population differentiation suggests genetic homogenization by metapopulation dynamics in a social spider

Mating systems and population dynamics influence genetic diversity and structure. Species that experience inbreeding and limited gene flow are expected to evolve isolated, divergent genetic lineages. Metapopulation dynamics with frequent extinctions and colonizations may, on the other hand, deplete and homogenize genetic variation, if extinction rate is sufficiently high compared to the effect of drift in local demes. We investigated these theoretical predictions empirically in social spiders that are highly inbred. Social spiders show intranest mating, female‐biased sex ratio, and frequent extinction and colonization events, factors that deplete genetic diversity within nests and populations and limit gene flow. We characterized population genetic structure in Stegodyphus sarasinorum, a social spider distributed across the Indian subcontinent. Species‐wide genetic diversity was estimated over approximately 2800 km from Sri Lanka to Himalayas, by sequencing 16 protein‐coding nuclear loci. We found 13 SNPs in 6592 bp (π = 0.00045) indicating low species‐wide nucleotide diversity. Three genetic lineages were strongly differentiated; however, only one fixed difference among them suggests recent divergence. This is consistent with a scenario of metapopulation dynamics that homogenizes genetic diversity across the species' range. Ultimately, low standing genetic variation may hamper a species' ability to track environmental change and render social inbreeding spiders ‘evolutionary dead‐ends’.     

Hidden suppression of sex ratio distortion suggests Red queen dynamics between Wolbachia and its dwarf spider host

Genetic conflict theory predicts strong selection for host nuclear factors suppressing endosymbiont effects on reproduction; however, evidence of these suppressors is currently scarce. This can either be caused by a low suppressor evolution rate, or if suppressors originate frequently, by rapid spread and concurrent masking of their activity by silencing the endosymbiont effect. To explore this, we use two populations of a dwarf spider with a similar female bias, caused by a Wolbachia infection. Using inter‐ and intrapopulation crosses, we determine that one of these populations demonstrates a higher suppressing capability towards Wolbachia despite having a similar population sex ratio. This suggests that spider and endosymbiont are locked in so‐called red queen dynamics where, despite continuous coevolution, average fitness remains the same, hence hiding the presence of the suppressor. Finding different suppressor activity in populations that even lack phenotypic differentiation (i.e. similar sex ratio) further supports the hypothesis that suppressors originate often, but are often hidden by their own mode of action by countering endosymbiont effects.     

Genetic diversity of a successful colonizer: isolated populations of Metrioptera roeselii regain variation at an unusually rapid rate

Newly founded isolated populations need to overcome detrimental effects of low genetic diversity. The establishment success of a population may therefore depend on various mechanisms such as assortative mating, purging of deleterious alleles, creation of new mutations and/or repeated inflow of new genotypes to reduce the effects of inbreeding and further loss of genetic variation. We compared the level of genetic variation in introduced populations of an insect species (Metrioptera roeselii) far beyond its natural distribution with levels found in their respective founder populations and coupled the data with timing since establishment. This allowed us to analyze if the introduced populations showed signs of temporal changes in genetic variation and have made it possible to evaluate underlying mechanisms. For this, we used neutral genetic markers, seven microsatellite loci and a 676–bp‐long sequence of the mtDNA COI gene. All tested indices (allelic richness, unbiased expected heterozygosity, effective size, haplotype diversity, and nucleotide diversity) except inbreeding coefficient had significantly higher values in populations within the founding populations inside the continuous area of the species distribution compared with the introduced populations. A logarithmic model showed a significant correlation of both allelic richness and unbiased expected heterozygosity with age of the isolated populations. Considering the species' inferred colonization history and likely introduction pathways, we suggest that multiple introductions are the main mechanism behind the temporal pattern observed. However, we argue that influences of assortative mating, directional selection, and effects of an exceptional high intrapopulation mutation rate may have impacts. The ability to regain genetic diversity at this level may be one of the main reasons why M. roeselii successfully continue to colonize northern Europe.     

Extremely low genetic variability and highly structured local populations of Arabidopsis thaliana at higher latitudes

The genetic diversity and population structure of Arabidopsis thaliana populations from Norway were studied and compared to a worldwide sample of A. thaliana to investigate the demographic history and elucidate possible colonization routes of populations at the northernmost species limit. We genotyped 282 individuals from 31 local populations using 149 single nucleotide polymorphism markers. A high level of population subdivision (F(ST) = 0.85 ± 0.007) was found indicating that A. thaliana is highly structured at the regional level. Significant relationships between genetic and geographical distances were found, suggesting an isolation by distance mode of evolution. Genetic diversity was much lower, and the level of linkage disequilibrium was higher in populations from the north (65-68°N) compared to populations from the south (59-62°N); this is consistent with a northward expansion pattern. A neighbour-joining tree showed that populations from northern Norway form a separate cluster, while the remaining populations are distributed over a few minor clusters. Minimal gene flow seems to have occurred between populations in different regions, especially between the geographically distant northern and southern populations. Our data suggest that northern populations represent a homogenous group that may have been established from a few founders during northward expansions, while populations in the central part of Norway constitute an admixed group established by founders of different origins, most probably as a result of human-mediated gene flow. Moreover, Norwegian populations appeared to be homogenous and isolated compared to a worldwide sample of A. thaliana, but they are still grouped with Swedish populations, which may indicate common colonization histories.     

Evolutionary divergence and possible incipient speciation in post-glacial populations of a cosmopolitan aquatic plant

Habitat configuration is expected to have a major influence on genetic exchange and evolutionary divergence among populations. Aquatic organisms occur in two fundamentally different habitat types, the sea and freshwater lakes, making them excellent models to study the contrasting effects of continuity vs. isolation on genetic divergence. We compared the divergence in post-glacial populations of a cosmopolitan aquatic plant, the pondweed Potamogeton pectinatus that simultaneously occurs in freshwater lakes and coastal marine sites. Relative levels of gene flow were inferred in 12 lake and 14 Baltic Sea populations in northern Germany using nine highly polymorphic microsatellite markers developed for P. pectinatus. We found highly significant isolation-by-distance in both habitat types (P < 0.001). Genetic differentiation increased approximately 2.5-times faster among freshwater populations compared with those from the Baltic Sea. As different levels of genetic drift or population history cannot explain these differences, higher population connectivity in the sea relative to freshwater populations is the most likely source of contrasting evolutionary divergence. These findings are consistent with the notion that freshwater angiosperms are more conducive to allopatric speciation than their life-history counterparts in the sea, the relative species poor seagrasses. Surprisingly, population pairs from different habitat types revealed almost maximal genetic divergence expected for complete reproductive isolation, regardless of their respective geographical distance. Hence, the barrier to gene flow between lake and sea habitat types cannot be due to dispersal limitation. We may thus have identified a case of rapid incipient speciation in post-glacial populations of a widespread aquatic plant.     

Between migration load and evolutionary rescue: dispersal,adaptation and the response of spatially structured populations to environmental change

The evolutionary potential of populations is mainly determined by population size and available genetic variance. However, the adaptability of spatially structured populations may also be affected by dispersal: positively by spreading beneficial mutations across sub-populations, but negatively by moving locally adapted alleles between demes. We develop an individual-based, two-patch, allelic model to investigate the balance between these opposing effects on a population''s evolutionary response to rapid climate change. Individual fitness is controlled by two polygenic traits coding for local adaptation either to the environment or to climate. Under conditions of selection that favour the evolution of a generalist phenotype (i.e. weak divergent selection between patches) dispersal has an overall positive effect on the persistence of the population. However, when selection favours locally adapted specialists, the beneficial effects of dispersal outweigh the associated increase in maladaptation for a narrow range of parameter space only (intermediate selection strength and low linkage among loci), where the spread of beneficial climate alleles is not strongly hampered by selection against non-specialists. Given that local selection across heterogeneous and fragmented landscapes is common, the complex effect of dispersal that we describe will play an important role in determining the evolutionary dynamics of many species under rapidly changing climate.     

Adaptation to a seasonally varying environment: a strong latitudinal cline in reproductive diapause combined with high gene flow in Drosophila montana

Adaptation to seasonal changes in the northern hemisphere includes an ability to predict the forthcoming cold season from gradual changes in environmental cues early enough to prepare for the harsh winter conditions. The magnitude and speed of changes in these cues vary between the latitudes, which induces strong selection pressures for local adaptation.We studied adaptation to seasonal changes in Drosophila montana, a northern maltfly, by defining the photoperiodic conditions leading to adult reproductive diapause along a latitudinal cline in Finland and by measuring genetic differentiation and the amount of gene flow between the sampling sites with microsatellites. Our data revealed a clear correlation between the latitude and the critical day length (CDL), in which half of the females of different cline populations enter photoperiodic reproductive diapause. There was no sign of limited gene flow between the cline populations, even though these populations showed isolation by distance. Our results show that local adaptation may occur even in the presence of high gene flow, when selection for locally adaptive life-history traits is strong. A wide range of variation in the CDLs of the fly strains within and between the cline populations may be partly due to gene flow and partly due to the opposing selection pressures for fly reproduction and overwinter survival. This variation in the timing of diapause will enhance populations' survival over the years that differ in the severity of the winter and in the length of the warm period and may also help them respond to long-term changes in environmental conditions.     

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.