Estimates of gene flow in Drosophila pseudoobscura determined from nucleotide sequence analysis of the alcohol dehydrogenase region

A history of studies of reproductive isolation between Drosophila pseudoobscura and D. persimilis

Drosophila pseudoobscura and D. persimilis are a sister species pair that have been used as a model for studies of reproductive isolation and speciation for almost 100 years owing to their close evolutionary history, well characterized genetic differences, and overlapping geographic distribution. There are extensive analyses of both pre- and post-zygotic isolation, including studies of courtship divergence, conspecific sperm precedence (CSP) and how reinforcement by natural selection may or may not act to strengthen isolation in sympatry. Post-zygotic analyses explore the underlying mechanics of reproductive isolation; how inversions may give rise to initial speciation events and misexpression of key genes typically found within inversion regions render hybrid offspring unfit or inviable. We aim here to present a history of studies of reproductive isolation between this species pair, looking at how the field has developed over the last century and identifying the open questions and gaps within the literature.

How chromosomal inversions reorient the evolutionary process

Inversions are structural mutations that reverse the sequence of a chromosome segment and reduce the effective rate of recombination in the heterozygous state. They play a major role in adaptation, as well as in other evolutionary processes such as speciation. Although inversions have been studied since the 1920s, they remain difficult to investigate because the reduced recombination conferred by them strengthens the effects of drift and hitchhiking, which in turn can obscure signatures of selection. Nonetheless, numerous inversions have been found to be under selection. Given recent advances in population genetic theory and empirical study, here we review how different mechanisms of selection affect the evolution of inversions. A key difference between inversions and other mutations, such as single nucleotide variants, is that the fitness of an inversion may be affected by a larger number of frequently interacting processes. This considerably complicates the analysis of the causes underlying the evolution of inversions. We discuss the extent to which these mechanisms can be disentangled, and by which approach.

Ancient gene drives: an evolutionary paradox

Selfish genetic elements such as selfish chromosomes increase their transmission rate relative to the rest of the genome and can generate substantial cost to the organisms that carry them. Such segregation distorters are predicted to either reach fixation (potentially causing population extinction) or, more commonly, promote the evolution of genetic suppression to restore transmission to equality. Many populations show rapid spread of segregation distorters, followed by the rapid evolution of suppression. However, not all drivers display such flux, some instead persisting at stable frequencies in natural populations for decades, perhaps hundreds of thousands of years, with no sign of suppression evolving or the driver spreading to fixation. This represents a major evolutionary paradox. How can drivers be maintained in the long term at stable frequencies? And why has suppression not evolved as in many other gene drive systems? Here, we explore potential factors that may explain the persistence of drive systems, focusing on the ancient sex-ratio driver in the fly Drosophila pseudoobscura. We discuss potential solutions to the evolutionary mystery of why suppression does not appear to have evolved in this system, and address how long-term stable frequencies of gene drive can be maintained. Finally, we speculate whether ancient drivers may be functionally and evolutionarily distinct to young drive systems.

How chromosomal rearrangements shape adaptation and speciation: Case studies in Drosophila pseudoobscura and its sibling species Drosophila persimilis

The gene arrangements of Drosophila have played a prominent role in the history of evolutionary biology from the original quantification of genetic diversity to current studies of the mechanisms for the origin and establishment of new inversion mutations within populations and their subsequent fixation between species supporting reproductive barriers. This review examines the genetic causes and consequences of inversions as recombination suppressors and the role that recombination suppression plays in establishing inversions in populations as they are involved in adaptation within heterogeneous environments. This often results in the formation of clines of gene arrangement frequencies among populations. Recombination suppression leads to the differentiation of the gene arrangements which may accelerate the accumulation of fixed genetic differences among populations. If these fixed mutations cause incompatibilities, then inversions pose important reproductive barriers between species. This review uses the evolution of inversions in Drosophila pseudoobscura and D. persimilis as a case study for how inversions originate, establish and contribute to the evolution of reproductive isolation.

The evolution of recombination rates in finite populations during ecological speciation

Recombination can impede ecological speciation with gene flow by mixing locally adapted genotypes with maladapted migrant genotypes from a divergent population. In such a scenario, suppression of recombination can be selectively favoured. However, in finite populations evolving under the influence of random genetic drift, recombination can also facilitate adaptation by reducing Hill-Robertson interference between loci under selection. In this case, increased recombination rates can be favoured. Although these two major effects on recombination have been studied individually, their joint effect on ecological speciation with gene flow remains unexplored. Using a mathematical model, we investigated the evolution of recombination rates in two finite populations that exchange migrants while adapting to contrasting environments. Our results indicate a two-step dynamic where increased recombination is first favoured (in response to the Hill-Robertson effect), and then disfavoured, as the cost of recombining locally with maladapted migrant genotypes increases over time (the maladaptive gene flow effect). In larger populations, a stronger initial benefit for recombination was observed, whereas high migration rates intensify the long-term cost of recombination. These dynamics may have important implications for our understanding of the conditions that facilitate incipient speciation with gene flow and the evolution of recombination in finite populations.

Genomics of Natural Populations: How Differentially Expressed Genes Shape the Evolution of Chromosomal Inversions in Drosophila pseudoobscura

Chromosomal rearrangements can shape the structure of genetic variation in the genome directly through alteration of genes at breakpoints or indirectly by holding combinations of genetic variants together due to reduced recombination. The third chromosome of Drosophila pseudoobscura is a model system to test hypotheses about how rearrangements are established in populations because its third chromosome is polymorphic for >30 gene arrangements that were generated by a series of overlapping inversion mutations. Circumstantial evidence has suggested that these gene arrangements are selected. Despite the expected homogenizing effects of extensive gene flow, the frequencies of arrangements form gradients or clines in nature, which have been stable since the system was first described >80 years ago. Furthermore, multiple arrangements exist at appreciable frequencies across several ecological niches providing the opportunity for heterokaryotypes to form. In this study, we tested whether genes are differentially expressed among chromosome arrangements in first instar larvae, adult females and males. In addition, we asked whether transcriptional patterns in heterokaryotypes are dominant, semidominant, overdominant, or underdominant. We find evidence for a significant abundance of differentially expressed genes across the inverted regions of the third chromosome, including an enrichment of genes involved in sensory perception for males. We find the majority of loci show additivity in heterokaryotypes. Our results suggest that multiple genes have expression differences among arrangements that were either captured by the original inversion mutation or accumulated after it reached polymorphic frequencies, providing a potential source of genetic variation for selection to act upon. These data suggest that the inversions are favored because of their indirect effect of recombination suppression that has held different combinations of differentially expressed genes together in the various gene arrangement backgrounds.

Temperature can shape a cline in polyandry, but only genetic variation can sustain it over time

Multiple mating by females (polyandry) is a widespread behavior occurring in diverse taxa, species, and populations. Polyandry can also vary widely within species, and individual populations, so that both monandrous and polyandrous females occur together. Genetic differences can explain some of this intraspecific variation in polyandry, but environmental factors are also likely to play a role. One environmental factor that influences many fundamental biological processes is temperature. Higher temperatures have been shown to directly increase remating in laboratory studies of insects. In the longer term, high temperature could also help to drive the evolution of larger-scale patterns of behavior by changing the context-dependent balance of costs and benefits of polyandry across environments. We examined the relative influence of rearing and mating temperatures on female remating in populations of Drosophila pseudoobscura that show a latitudinal cline in polyandry in nature, using a range of ecologically relevant temperatures. We found that females of all genotypes remated more at cooler temperatures, which fits with the observation of higher average frequencies of polyandry at higher latitudes in this species. However, the impact of temperature was outweighed by the strong genetic control of remating in females in this species. It is likely that genetic factors provide the primary explanation for the latitudinal cline in polyandry in this species.

Recombining without Hotspots: A Comprehensive Evolutionary Portrait of Recombination in Two Closely Related Species of Drosophila

Meiotic recombination rate varies across the genome within and between individuals, populations, and species in virtually all taxa studied. In almost every species, this variation takes the form of discrete recombination hotspots, determined in some mammals by a protein called PRDM9. Hotspots and their determinants have a profound effect on the genomic landscape, and share certain features that extend across the tree of life. Drosophila, in contrast, are anomalous in their absence of hotspots, PRDM9, and other species-specific differences in the determination of recombination. To better understand the evolution of meiosis and general patterns of recombination across diverse taxa, we present a truly comprehensive portrait of recombination across time, combining recently published cross-based contemporary recombination estimates from each of two sister species with newly obtained linkage-disequilibrium-based historic estimates of recombination from both of these species. Using Drosophila pseudoobscura and Drosophila miranda as a model system, we compare recombination rate between species at multiple scales, and we suggest that Drosophila replicate the pattern seen in human-chimpanzee in which recombination rate is conserved at broad scales. We also find evidence of a species-wide recombination modifier(s), resulting in both a present and historic genome-wide elevation of recombination rates in D. miranda, and identify broad scale effects on recombination from the presence of an inversion. Finally, we reveal an unprecedented view of the distribution of recombination in D. pseudoobscura, illustrating patterns of linked selection and where recombination is taking place. Overall, by combining these estimation approaches, we highlight key similarities and differences in recombination between Drosophila and other organisms.

Evidence for stabilizing selection on codon usage in chromosomal rearrangements of Drosophila pseudoobscura

There has been a renewed interest in investigating the role of stabilizing selection acting on genome-wide traits such as codon usage bias. Codon bias, when synonymous codons are used at unequal frequencies, occurs in a wide variety of taxa. Standard evolutionary models explain the maintenance of codon bias through a balance of genetic drift, mutation and weak purifying selection. The efficacy of selection is expected to be reduced in regions of suppressed recombination. Contrary to observations in Drosophila melanogaster, some recent studies have failed to detect a relationship between the recombination rate, intensity of selection acting at synonymous sites, and the magnitude of codon bias as predicted under these standard models. Here, we examined codon bias in 2798 protein coding loci on the third chromosome of D. pseudoobscura using whole-genome sequences of 47 individuals, representing five common third chromosome gene arrangements. Fine-scale recombination maps were constructed using more than 1 million segregating sites. As expected, recombination was demonstrated to be significantly suppressed between chromosome arrangements, allowing for a direct examination of the relationship between recombination, selection, and codon bias. As with other Drosophila species, we observe a strong mutational bias away from the most frequently used codons. We find the rate of synonymous and nonsynonymous polymorphism is variable between different amino acids. However, we do not observe a reduction in codon bias or the strength of selection in regions of suppressed recombination as expected. Instead, we find that the interaction between weak stabilizing selection and mutational bias likely plays a role in shaping the composition of synonymous codons across the third chromosome in D. pseudoobscura.

Evidence that natural selection on codon usage in Drosophila pseudoobscura varies across codons

Like other species of Drosophila, Drosophila pseudoobscura has a distinct bias toward the usage of C- and G-ending codons. Previous studies have indicated that this bias is due, at least in part, to natural selection. Codon bias clearly differs among amino acids (and other codon classes) in Drosophila, which may reflect differences in the intensity of selection on codon usage. Ongoing natural selection on synonymous codon usage should be reflected in the shapes of the site frequency spectra of derived states at polymorphic positions. Specifically, regardless of other demographic effects on the spectrum, it should be shifted toward higher values for changes from less-preferred to more-preferred codons, and toward lower values for the converse. If the intensity of natural selection is increased, shifts in the site frequency spectra should be more pronounced. A total of 33,729 synonymous polymorphic sites on Chromosome 2 in D. pseudoobscura were analyzed. Shifts in the site frequency spectra are consistent with differential intensity of natural selection on codon usage, with stronger shifts associated with higher codon bias. The shifts, in general, are greater for polymorphic synonymous sites than for polymorphic intron sites, also consistent with natural selection. However, unlike observations in D. melanogaster, codon bias is not reduced in areas of low recombination in D. pseudoobscura; the site frequency spectrum signal for selection on codon usage remains strong in these regions. However, diversity is reduced, as expected. It is possible that estimates of low recombination reflect a recent change in recombination rate.

Does polyandry control population sex ratio via regulation of a selfish gene?

The extent of female multiple mating (polyandry) can strongly impact on the intensity of sexual selection, sexual conflict, and the evolution of cooperation and sociality. More subtly, polyandry may protect populations against intragenomic conflicts that result from the invasion of deleterious selfish genetic elements (SGEs). SGEs commonly impair sperm production, and so are likely to be unsuccessful in sperm competition, potentially reducing their transmission in polyandrous populations. Here, we test this prediction in nature. We demonstrate a heritable latitudinal cline in the degree of polyandry in the fruitfly Drosophila pseudoobscura across the USA, with northern population females remating more frequently in both the field and the laboratory. High remating was associated with low frequency of a sex-ratio-distorting meiotic driver in natural populations. In the laboratory, polyandry directly controls the frequency of the driver by undermining its transmission. Hence we suggest that the cline in polyandry represents an important contributor to the cline in sex ratio in nature. Furthermore, as the meiotic driver causes sex ratio bias, variation in polyandry may ultimately determine population sex ratio across the USA, a dramatic impact of female mating decisions. As SGEs are ubiquitous it is likely that the reduction of intragenomic conflict by polyandry is widespread.

Recent selection on the Y-to-dot translocation in Drosophila pseudoobscura

The Drosophila pseudoobscura dot chromosome acquired genes from the ancestral Drosophila Y chromosome in a Y-to-dot translocation event that occurred between 12.7 and 20.8 Ma. The formerly Y-linked genes mostly retained their testis-specific expression but shrank drastically in size, mostly through intron reduction, since becoming part of the dot chromosome in this species. We investigated the impact of this translocation on the evolution of the both the Y-to-dot translocated region and the original segments of the dot chromosome in D. pseudoobscura. Our survey of polymorphism and divergence across the chromosome reveals a reduction in variation, a deletion polymorphism segregating at high frequency, and a shift in the frequency spectra, all consistent with a history of recent selective sweeps in the Y-to-dot translocated region but not on the rest of the dot chromosome. We do find evidence for recombination primarily as gene conversion on the dot chromosome; however, predicted recombination events are restricted to the part of the dot chromosome outside the translocation. It therefore appears that recombination has resulted in a degree of decoupling between the ancestral Y region and the conserved region of the dot chromosome.

Restricted gene flow within and between rapidly diverging Neotropical plant species

Speciation involves the evolution of traits and genetic differences that contribute to reproductive isolation and the cessation of gene flow, and studying closely related species and divergent populations gives insight into how these phenomena proceed. Here, we document patterns of gene flow within and between two members of a rapid Neotropical species radiation, Costus pulverulentus and Costus scaber (Costaceae). These species co-occur in the tropical rainforest and share pollinators, but are reproductively isolated by a series of prezygotic barriers, some of which show evidence of reinforcement at sympatric sites. Here, we genotype microsatellite markers in plants from eight sites that span the geographical range of both species, including four sympatric sites. We also genotype putative hybrids found at two sympatric sites. We find high levels of genetic isolation among populations within each species and low but detectable levels of introgression between species at sympatric sites. Putative hybrids identified by morphology are consistent with F1 or more advanced hybrids. Our results highlight the effectiveness of prezygotic isolating mechanisms at maintaining species boundaries in young radiations and provide empirical data on levels of gene flow consistent with reinforcement.

Within-niche constant selection of Drosophila inversions?

Extensive and fruitful work is being devoted for more than 70 years to elucidate the fine points of the maintenance of inversion polymorphisms of the genus Drosophila. Recent studies have resumed selection in heterogeneous environments (or niches) as a major underlying mechanism for these balanced polymorphisms. In those studies, constant selection within niches is assumed throughout although this assumption is since long known not to hold. In the present communication it is sustained that the results in those studies are robust in the face of this fact. To that end, this communication deals with a particular long-lasting question within this topic-whether the minimal model of constant viability selection (MCV, assuming frequency-, sex-, and stage-independent adaptive values) suffices to reproduce the trajectories of frequencies of Drosophila chromosomal arrangements observed in experimental populations along generations under homogeneous environments. Fitness estimates are here obtained from published trajectories of frequencies using a maximum likelihood approach, and relevant literature is revised in the light of these new analyses, pointing to an affirmative answer to that question.

Geographic selection in the small heat shock gene complex differentiating populations of Drosophila pseudoobscura

Environmental temperature plays a crucial role in determining a species distribution and abundance by affecting individual physiological processes, metabolic activities, and developmental rates. Many studies have identified clinal variation in phenotypes associated with response to environmental stresses, but variation in traits associated with climatic adaptation directly attributed to sequence variation within candidate gene regions has been difficult to identify. Insect heat shock genes are possible agents of thermal tolerance because of their involvement in protein folding, traffic, protection, and renaturation at the cellular level in response to temperature stress. Previously, members of the Drosophila small heat shock protein (sHSP) complex (Hsp23, Hsp26, Hsp27, Hsp67Ba) have been implicated as candidate climatic adaptation genes; therefore, this research examines sequence variation at these genes in 2 distant populations of Drosophila pseudoobscura. Flies from Tempe, AZ (n = 30) and Cheney, WA (n = 17) were used in the study. We identify high differentiation in the heat-shock complex (F(ST) : 0.219**, 0.262*, 0.279***, 0.166 not significant) as compared with neighboring genes and Tajima's D values indicative of balancing selection (Mann-Whitney U = 38, n(1) = 10 n(2) = 4, P < 0.05 two-tailed), both of which are suggestive of such climatic adaptation.

Adaptive evolution of genes duplicated from the Drosophila pseudoobscura neo-X chromosome

Drosophila X chromosomes are disproportionate sources of duplicated genes, and these duplications are usually the result of retrotransposition of X-linked genes to the autosomes. The excess duplication is thought to be driven by natural selection for two reasons: X chromosomes are inactivated during spermatogenesis, and the derived copies of retroposed duplications tend to be testis expressed. Therefore, autosomal derived copies of retroposed genes provide a mechanism for their X-linked paralogs to "escape" X inactivation. Once these duplications have fixed, they may then be selected for male-specific functions. Throughout the evolution of the Drosophila genus, autosomes have fused with X chromosomes along multiple lineages giving rise to neo-X chromosomes. There has also been excess duplication from the two independent neo-X chromosomes that have been examined--one that occurred prior to the common ancestor of the willistoni species group and another that occurred along the lineage leading to Drosophila pseudoobscura. To determine what role natural selection plays in the evolution of genes duplicated from the D. pseudoobscura neo-X chromosome, we analyzed DNA sequence divergence between paralogs, polymorphism within each copy, and the expression profiles of these duplicated genes. We found that the derived copies of all duplicated genes have elevated nonsynonymous polymorphism, suggesting that they are under relaxed selective constraints. The derived copies also tend to have testis- or male-biased expression profiles regardless of their chromosome of origin. Genes duplicated from the neo-X chromosome appear to be under less constraints than those duplicated from other chromosome arms. We also find more evidence for historical adaptive evolution in genes duplicated from the neo-X chromosome, suggesting that they are under a unique selection regime in which elevated nonsynonymous polymorphism provides a large reservoir of functional variants, some of which are fixed by natural selection.

Polyandry prevents extinction

Females of most animal species are polyandrous, with individual females usually mating with more than one male. However, the ubiquity of polyandry remains enigmatic because of the potentially high costs to females of multiple mating. Current theory to account for the high prevalence of polyandry largely focuses on its benefits to individual females. There are also higher-level explanations for the high incidence of polyandry-polyandrous clades may speciate more rapidly. Here we test the hypothesis that polyandry may also reduce population extinction risk. We demonstrate that mating with multiple males protects populations of the fruit fly Drosophila pseudoobscura against extinction caused by a "selfish" sex-ratio-distorting element. Thus, the frequency of female multiple mating in nature may be associated not only with individual benefits to females of this behavior but also with increased persistence over time of polyandrous species and populations. Furthermore, we show that female remating behavior can determine the frequency of sex-ratio distorters in populations. This may also be true for many other selfish genetic elements in natural populations.

Fine-scale mapping of recombination rate in Drosophila refines its correlation to diversity and divergence

Regional rates of recombination often correlate with levels of nucleotide diversity, and either selective or neutral hypotheses can explain this relationship. Regional recombination rates also correlate with nucleotide differences between human and chimpanzee, consistent with models where recombination is mutagenic; however, a lack of correlation is observed in the Drosophila melanogaster group, consistent with models invoking natural selection. Here, we revisit the relationship among recombination, diversity, and interspecies difference by generating empirical estimates of these parameters in Drosophila pseudoobscura. To measure recombination rate, we genotyped 1,294 backcross hybrids at 50 markers across the largest assembled linkage group in this species. Genome-wide diversity was estimated by sequencing a second isolate of D. pseudoobscura at shallow coverage. Alignment to the sequenced genome of the closely related species, Drosophila persimilis, provided nucleotide site orthology. Our findings demonstrate that scale is critical in determining correlates to recombination rate: fine-scale cross-over rate estimates are far stronger predictors of both diversity and interspecies difference than broad-scale estimates. The correlation of fine-scale recombination rate to diversity and interspecies difference appears to be genome-wide, evidenced by examination of an X-linked region in greater detail. Because we observe a strong correlation of cross-over rate with interspecies difference, even after correcting for segregating ancestral variation, we suggest that both mutagenic and selective forces generate these correlations, the latter in regions of low crossing over. We propose that it is not cross-overs per se that are mutagenic, but rather repair of DNA double-strand break precursors via crossing over and gene conversion.

Codon usage in twelve species of Drosophila

Background

Codon usage bias (CUB), the uneven use of synonymous codons, is a ubiquitous observation in virtually all organisms examined. The pattern of codon usage is generally similar among closely related species, but differs significantly among distantly related organisms, e.g., bacteria, yeast, and Drosophila. Several explanations for CUB have been offered and some have been supported by observations and experiments, although a thorough understanding of the evolutionary forces (random drift, mutation bias, and selection) and their relative importance remains to be determined. The recently available complete genome DNA sequences of twelve phylogenetically defined species of Drosophila offer a hitherto unprecedented opportunity to examine these problems. We report here the patterns of codon usage in the twelve species and offer insights on possible evolutionary forces involved.

Results

(1) Codon usage is quite stable across 11/12 of the species: G- and especially C-ending codons are used most frequently, thus defining the preferred codons. (2) The only amino acid that changes in preferred codon is Serine with six species of the melanogaster group favoring TCC while the other species, particularly subgenus Drosophila species, favor AGC. (3) D. willistoni is an exception to these generalizations in having a shifted codon usage for seven amino acids toward A/T in the wobble position. (4) Amino acids differ in their contribution to overall CUB, Leu having the greatest and Asp the least. (5) Among two-fold degenerate amino acids, A/G ending amino acids have more selection on codon usage than T/C ending amino acids. (6) Among the different chromosome arms or elements, genes on the non-recombining element F (dot chromosome) have the least CUB, while genes on the element A (X chromosome) have the most. (7) Introns indicate that mutation bias in all species is approximately 2:1, AT:GC, the opposite of codon usage bias. (8) There is also evidence for some overall regional bias in base composition that may influence codon usage.

Conclusion

Overall, these results suggest that natural selection has acted on codon usage in the genus Drosophila, at least often enough to leave a footprint of selection in modern genomes. However, there is evidence in the data that random forces (drift and mutation) have also left patterns in the data, especially in genes under weak selection for codon usage for example genes in regions of low recombination. The documentation of codon usage patterns in each of these twelve genomes also aids in ongoing annotation efforts.

Meiotic transmission of Drosophila pseudoobscura chromosomal arrangements

Drosophila pseudoobscura harbors a rich gene arrangement polymorphism on the third chromosome generated by a series of overlapping paracentric inversions. The arrangements suppress recombination in heterokaryotypic individuals, which allows for the selective maintenance of coadapted gene complexes. Previous mapping experiments used to determine the degree to which recombination is suppressed in gene arrangement heterozygotes produced non-recombinant progeny in non-Mendelian ratios. The deviations from Mendelian expectations could be the result of viability differences between wild and mutant chromosomes, meiotic drive because of achiasmate pairing of homologues in heterokaryotypic females during meiosis, or a combination of both mechanisms. The possibility that the frequencies of the chromosomal arrangements in natural populations are affected by mechanisms other than adaptive selection led us to consider these hypotheses. We performed reciprocal crosses involving both heterozygous males and females to determine if the frequency of the non-recombinant progeny deviates significantly from Mendelian expectations and if the frequencies deviate between reciprocal crosses. We failed to observe non-Mendelian ratios in multiple crosses, and the frequency of the non-recombinant classes differed in only one of five pairs of reciprocal crosses despite sufficient power to detect these differences in all crosses. Our results indicate that deviations from Mendelian expectations in recombination experiments involving the D. pseudoobscura inversion system are most likely due to fitness differences of gene arrangement karyotypes in different environments.

Evaluation of the genomic extent of effects of fixed inversion differences on intraspecific variation and interspecific gene flow in Drosophila pseudoobscura and D. persimilis

There is increasing evidence that chromosomal inversions may facilitate the formation or persistence of new species by allowing genetic factors conferring species-specific adaptations or reproductive isolation to be inherited together and by reducing or eliminating introgression. However, the genomic domain of influence of the inverted regions on introgression has not been carefully studied. Here, we present a detailed study on the consequences that distance from inversion breakpoints has had on the inferred level of gene flow and divergence between Drosophila pseudoobscura and D. persimilis. We identified the locations of the inversion breakpoints distinguishing D. pseudoobscura and D. persimilis in chromosomes 2, XR, and XL. Population genetic data were collected at specific distances from the inversion breakpoints of the second chromosome and at two loci inside the XR and XL inverted regions. For loci outside the inverted regions, we found that distance from the nearest inversion breakpoint had a significant effect on several measures of divergence and gene flow between D. pseudoobscura and D. persimilis. The data fitted a logarithmic relationship, showing that the suppression of crossovers in inversion heterozygotes also extends to loci located outside the inversion but close to it (within 1-2 Mb). Further, we detected a significant reduction in nucleotide variation inside the inverted second chromosome region of D. persimilis and near one breakpoint, consistent with a scenario in which this inversion arose and was fixed in this species by natural selection.

Estimating the contribution of mutation, recombination and gene conversion in the generation of haplotypic diversity

Recombination occurs through both homologous crossing over and homologous gene conversion during meiosis. The contribution of recombination relative to mutation is expected to be dramatically reduced in inbreeding organisms. We report coalescent-based estimates of the recombination parameter (rho) relative to estimates of the mutation parameter (theta) for 18 genes from the highly self-fertilizing grass, wild barley, Hordeum vulgare ssp. spontaneum. Estimates of rho/theta are much greater than expected, with a mean rho/theta approximately 1.5, similar to estimates from outcrossing species. We also estimate rho with and without the contribution of gene conversion. Genotyping errors can mimic the effect of gene conversion, upwardly biasing estimates of the role of conversion. Thus we report a novel method for identifying genotyping errors in nucleotide sequence data sets. We show that there is evidence for gene conversion in many large nucleotide sequence data sets including our data that have been purged of all detectable sequencing errors and in data sets from Drosophila melanogaster, D. simulans, and Zea mays. In total, 13 of 27 loci show evidence of gene conversion. For these loci, gene conversion is estimated to contribute an average of twice as much as crossing over to total recombination.

Mechanisms of genetic exchange within the chromosomal inversions of Drosophila pseudoobscura

We have used the inversion system of Drosophila pseudoobscura to investigate how genetic flux occurs among the gene arrangements. The patterns of nucleotide polymorphism at seven loci were used to infer gene conversion events between pairs of different gene arrangements. We estimate that the average gene conversion tract length is 205 bp and that the average conversion rate is 3.4 x 10(-6), which is 2 orders of magnitude greater than the mutation rate. We did not detect gene conversion events between all combinations of gene arrangements even though there was sufficient nucleotide variation for detection and sufficient opportunity for exchanges to occur. Genetic flux across the inverted chromosome resulted in higher levels of differentiation within 0.1 Mb of inversion breakpoints, but a slightly lower level of differentiation in central inverted regions. No gene conversion events were detected within 17 kb of an inversion breakpoint suggesting that the formation of double-strand breaks is reduced near rearrangement breakpoints in heterozygotes. At least one case where selection rather than proximity to an inversion breakpoint is responsible for reduction in polymorphism was identified.

Sexual Isolation Between Sympatric and Allopatric Populations of Drosophila pseudoobscura and D. persimilis*

According to reinforcement theory, sexual isolation between species in sympatry is strengthened by natural selection against maladaptive hybrids. Drosophila pseudoobscura and D. persimilis from four locations where these species are sympatric, and from three locations where only D. pseudoobscura has been found, were utilized in studies of sexual isolation. Multiple-choice observation chambers were used to record matings between sympatric and allopatric strains of the two species. There was a wide variation in sexual isolation between the two species in the four localities we studied. The average isolation index for sympatric strains of the species was not significantly different from the average index for allopatric strains. There were no meaningful differences between the isolation indices in sympatric and allopatric strains of the species. The failure to find a relationship is likely the result of gene flow among populations within the two species.

The genetics of speciation by reinforcement

Reinforcement occurs when natural selection strengthens behavioral discrimination to prevent costly interspecies matings, such as when matings produce sterile hybrids. This evolutionary process can complete speciation, thereby providing a direct link between Darwin's theory of natural selection and the origin of new species. Here, by examining a case of speciation by reinforcement in Drosophila, we present the first high-resolution genetic study of variation within species for female mating discrimination that is enhanced by natural selection. We show that reinforced mating discrimination is inherited as a dominant trait, exhibits variability within species, and may be influenced by a known set of candidate genes involved in olfaction. Our results show that the genetics of reinforced mating discrimination is different from the genetics of mating discrimination between species, suggesting that overall mating discrimination might be a composite phenomenon, which in Drosophila could involve both auditory and olfactory cues. Examining the genetics of reinforcement provides a unique opportunity for both understanding the origin of new species in the face of gene flow and identifying the genetic basis of adaptive female species preferences, two major gaps in our understanding of speciation.

Mating discrimination between two species of Drosophila is more pronounced where they hybridize and genes related to odor appear responsible for this "reinforcement," thus providing insights into the genetics of speciation

Evolutionary genomics of inversions in Drosophila pseudoobscura: evidence for epistasis

Drosophila pseudoobscura harbors a rich polymorphism for paracentric inversions on the third chromosome, and the clines in the inversion frequencies across the southwestern United States indicate that strong natural selection operates on them. Isogenic inversion strains were made from isofemale lines collected from four localities, and eight molecular markers were mapped on the third chromosome. Nucleotide diversity was measured for these loci and formed the basis of an evolutionary genomic analysis. The loci were differentiated among inversions. The inversions did not show significant differences among populations, however, likely the result of extensive gene flow among populations. Some loci had significant reductions in nucleotide diversity within inversions compared with interspecies divergence, suggesting that these loci are near inversion breakpoints or are near targets of directional selection. Linkage disequilibrium (LD) levels tended to decrease with distance between loci, indicating that some genetic exchange occurs among gene arrangements despite the presence of inversions. In some cases, however, adjacent genes had low levels of interlocus LD and loosely linked genes had high levels of interlocus LD, suggesting strong epistatic selection. Our results support the hypothesis that the inversions of D. pseudoobscura have emerged as suppressors of recombination to maintain positive epistatic relationships among loci within gene arrangements that developed as the species adapted to a heterogeneous environment.

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).

Inferring the history of speciation from multilocus DNA sequence data: the case of Drosophila pseudoobscura and close relatives

The divergence of Drosophila pseudoobscura from its close relatives, D. persimilis and D. pseudoobscura bogotana, was examined using the pattern of DNA sequence variation in a common set of 50 inbred lines at 11 loci from diverse locations in the genome. Drosophila pseudoobscura and D. persimilis show a marked excess of low-frequency variation across loci, consistent with a model of recent population expansion in both species. The different loci vary considerably, both in polymorphism levels and in the levels of polymorphisms that are shared by different species pairs. A major question we address is whether these patterns of shared variation are best explained by gene flow or by persistence since common ancestry. A new test of gene flow, based on patterns of linkage disequilibrium, is developed. The results from these, and other tests, support a model in which D. pseudoobscura and D. persimilis have exchanged genes at some loci. However, the pattern of variation suggests that most gene flow, although occurring after speciation began, was not recent. There is less evidence of gene flow between D. pseudoobscura and D. p. bogotana. The results are compared with recent work on the genomic locations of genes that contribute to reproductive isolation between D. pseudoobscura and D. persimilis. We show that there is a good correspondence between the genomic regions associated with reproductive isolation and the regions that show little or no evidence of gene flow.

Protein variation in Adh and Adh-related in Drosophila pseudoobscura. Linkage disequilibrium between single nucleotide polymorphisms and protein alleles

A 3.5-kb segment of the alcohol dehydrogenase (Adh) region that includes the Adh and Adh-related genes was sequenced in 139 Drosophila pseudoobscura strains collected from 13 populations. The Adh gene encodes four protein alleles and rejects a neutral model of protein evolution with the McDonald-Kreitman test, although the number of segregating synonymous sites is too high to conclude that adaptive selection has operated. The Adh-related gene encodes 18 protein haplotypes and fails to reject an equilibrium neutral model. The populations fail to show significant geographic differentiation of the Adh-related haplotypes. Eight of 404 single nucleotide polymorphisms (SNPs) in the Adh region were in significant linkage disequilibrium with three ADHR protein alleles. Coalescent simulations with and without recombination were used to derive the expected levels of significant linkage disequilibrium between SNPs and 18 protein haplotypes. Maximum levels of linkage disequilibrium are expected for protein alleles at moderate frequencies. In coalescent models without recombination, linkage disequilibrium decays between SNPs and high frequency haplotypes because common alleles mutate to haplotypes that are rare or that reach moderate frequency. The implication of this study is that linkage disequilibrium mapping has the highest probability of success with disease-causing alleles at frequencies of 10%.

Adult female Drosophila pseudoobscura survive and carry fertile sperm through long periods in the cold: populations are unlikely to suffer substantial bottlenecks in overwintering

To assess whether, while overwintering, natural populations of Drosophila pseudoobscura are likely to experience substantial bottlenecks in their numbers and genotypes, laboratory tests of the cold sensitivities of each stage of the life history and reproduction were undertaken. Three genetically distinctive lineages established from flies caught at high elevation were used for testing in temperatures likely to persist in protected pockets of fermenting deciduous leaf fall in overwintering sites. Sensitivities to cold of each stage in development were measured as frequencies of survival to adulthood following a period in 5 degrees C in a particular stage. The cold sensitivity of adults was measured as the survival in and following cold stays in adulthood. It was found that cold sensitivity decreases as development progresses, but that only adults (females more than males) are able to withstand long periods in the cold. The cold sensitivity of reproductive capacity of males was scored as their success in mating following a two-month cold stay, and of females as the numbers laying fertile eggs following periods of months in the cold. Both males and females maintain reproductive capacity. Of particular significance, however, is that even after six months in the cold females are able to restart production of eggs and these eggs may be fertilized by the sperm of matings prior to their cold stay. Thus, a substantial proportion of overwintering genomes must be those of adult females and those of the sperm carried by females from matings in the previous summer. This simple finding strongly suggests that populations are not likely to suffer substantial bottlenecks while overwintering. Further, it indicates how arrays of genetic variation may be maintained through winters and largely avoid winter selective pressures. Frequent migration between populations is therefore not required to maintain the variation commonly found in populations throughout the species range.

Nucleotide polymorphism and natural selection at the pantophysin (Pan I) locus in the Atlantic cod, Gadus morhua (L.)

Molecular studies of nucleotide sequence variation have rarely attempted to test hypotheses related to geographically varying patterns of natural selection. The present study tested the role of spatially varying selection in producing significant linkage disequilibrium and large differences in the frequencies of two common alleles at the pantophysin (Pan I) locus among five populations of the Atlantic cod, Gadus morhua. Nucleotide sequences of 124 Pan I alleles showed strong evidence for an unusual mix of balancing and directional selection but no evidence of stable geographically varying selection. The alleles were highly divergent at both the nucleotide level (differing on average by 19 mutations) and at amino acid level (each having experienced three amino acid substitutions since diverging from a common ancestral allele). All six amino acid substitutions occurred in a 56-residue intravesicular loop (IV1 domain) of the vesicle protein and each involved a radical change. An analysis of molecular variation revealed significant heterogeneity in the frequencies of recently derived mutations segregating within both allelic classes, suggesting that two selective sweeps may be presently occurring among populations. The dynamic nature of the Pan I polymorphism in G. morhua and clear departure from equilibrium conditions invalidate a simple model of spatially varying selection.

Molecular population genetics of X-linked genes in Drosophila pseudoobscura

This article presents a nucleotide sequence analysis of 500 bp determined in each of five X-linked genes, runt, sisterlessA, period, esterase 5, and Heat-shock protein 83, in 40 Drosophila pseudoobscura strains collected from two populations. Estimates of the neutral migration parameter for the five loci show that gene flow among D. pseudoobscura populations is sufficient to homogenize inversion frequencies across the range of the species. Nucleotide diversity at each locus fails to reject a neutral model of molecular evolution. The sample of 40 chromosomes included six Sex-ratio inversions, a series of three nonoverlapping inversions that are associated with a strong meiotic drive phenotype. The selection driven by the Sex-ratio meiotic drive element has not fixed variation across the X chromosome of D. pseudoobscura because, while significant linkage disequilibrium was observed within the sisterlessA, period, and esterase 5 genes, we did not find evidence for nonrandom association among loci. The Sex-ratio chromosome was estimated to be 25,000 years old based on the decomposition of linkage disequilibrium between esterase 5 and Heat-shock protein 83 or 1 million years old based on the net divergence of esterase 5 between Standard and Sex-ratio chromosomes. Genetic diversity was depressed within esterase 5 within Sex-ratio chromosomes, while the four other genes failed to show a reduction in heterozygosity in the Sex-ratio background. The reduced heterogeneity in esterase 5 is due either to its location near one of the Sex-ratio inversion breakpoints or that it is closely linked to a gene or genes responsible for the Sex-ratio meiotic drive system.

Reinforcement and other consequences of sympatry

The pattern of greater species mating discrimination between sympatric taxa than between allopatric taxa has been attributed to the strengthening of mate discrimination to avoid maladaptive hybridization. This process, termed reinforcement, has been highly contentious, particularly with regard to its role in speciation. Here, I review some recent studies of reinforcement, discuss alternative explanations for the pattern of greater species discrimination in sympatric taxa, and point to some new directions that may help to clarify the evolutionary forces involved. In particular, we need more ecological work on putative cases of reinforcement, more theoretical models that give diagnostic predictions of reinforcement relative to other modes of divergence, and empirical studies to evaluate these diagnostic predictions.

DNA sequence variation and the recombinational landscape in Drosophila pseudoobscura: a study of the second chromosome

The relationship between rates of recombination and DNA sequence polymorphism was analyzed for the second chromosome of Drosophila pseudoobscura. We constructed integrated genetic and physical maps of this chromosome using molecular markers at 10 loci spanning most of its physical length. The total length of the map was 128.2 cM, almost twice that of the homologous chromosome arm (3R) in D. melanogaster. There appears to be very little centromeric suppression of recombination, and rates of recombination are quite uniform across most of the chromosome. Levels of sequence variation (theta(W), based on the number of segregating sites) at seven loci (tropomyosin 1, Rhodopsin 3, Rhodopsin 1, bicoid, Xanthine dehydrogenase, Myosin light chain 1, and ribosomal protein 49) varied from 0.0036 to 0.0167. Generally consistent with earlier studies, the average estimate of theta(W) at total sites is 1.5-fold higher than that in D. melanogaster, while average theta(W) at silent sites is almost 3-fold higher. These estimates of variation were analyzed in the context of a background selection model under the same parameters of mutation rate and selection as have been proposed for D. melanogaster. It is likely that a significant fraction of the higher level of sequence variation in D. pseudoobscura can be explained by differences in regional rates of recombination rather than a larger species-level effective population size. However, the distribution of variation among synonymous, nonsynonymous, and noncoding sites appears to be quite different between the species, making direct comparisons of neutral variation, and hence inferences about effective population size, difficult. Tajima's D statistics for 6 out of the 7 loci surveyed are negative, suggesting that D. pseudoobscura may have experienced a rapid population expansion in the recent past or, alternatively, that slightly deleterious mutations constitute an important component of standing variation in this species.

Natural selection and the frequency distributions of "silent" DNA polymorphism in Drosophila

In Escherichia coli, Saccharomyces cerevisiae, and Drosophila melanogaster, codon bias may be maintained by a balance among mutation pressure, genetic drift, and natural selection favoring translationally superior codons. Under such an evolutionary model, silent mutations fall into two fitness categories: preferred mutations that increase codon bias and unpreferred changes in the opposite direction. This prediction can be tested by comparing the frequency spectra of synonymous changes segregating within populations; natural selection will elevate the frequencies of advantageous mutations relative to that of deleterious changes. The frequency distributions of preferred and unpreferred mutations differ in the predicted direction among 99 alleles of two D. pseudoobscura genes and five alleles of eight D. simulans genes. This result confirms the existence of fitness classes of silent mutations. Maximum likelihood estimates suggest that selection intensity at silent sites is, on average, very weak in both D. pseudoobscura and D. simulans (magnitude of NS approximately 1). Inference of evolutionary processes from within-species sequence variation is often hindered by the assumption of a stationary frequency distribution. This assumption can be avoided when identifying the action of selection and tested when estimating selection intensity.

The speciation history of Drosophila pseudoobscura and close relatives: inferences from DNA sequence variation at the period locus

Thirty-five period locus sequences from Drosophila pseudoobscura and its siblings species, D. p. bogotana, D. persimilis, and D. miranda, were studied. A large amount of variation was found within D. pseudoobscura and D. persimilis, consistent with histories of large effective population sizes. D. p. bogotana, however, has a severe reduction in diversity. Combined analysis of per with two other loci, in both D. p. bogotana and D. pseudoobscura, strongly suggest this reduction is due to recent directional selection at or near per within D. p. bogotana. Since D. p. bogotana is highly variable and shares variation with D. pseudoobscura at other loci, the low level of variation at per within D. p. bogotana can not be explained by a small effective population size or by speciation via founder event. Both D. pseudoobscura and D. persimilis have considerable intraspecific gene flow. A large portion of one D. persimilis sequence appears to have arisen via introgression from D. pseudoobscura. The time of this event appears to be well after the initial separation of these two species. The estimated times since speciation are one mya for D. pseudoobscura and D. persimilis and 2 mya since the formation of D. miranda.