Genetics of hybrid male sterility among strains and species in the Drosophila pseudoobscura species group

Taxa in the early stages of speciation may bear intraspecific allelic variation at loci conferring barrier traits in hybrids such as hybrid sterility. Additionally, hybridization may spread alleles that confer barrier traits to other taxa. Historically, few studies examine within- and between-species variation at loci conferring reproductive isolation. Here, we test for allelic variation within Drosophila persimilis and within the Bogota subspecies of D. pseudoobscura at regions previously shown to contribute to hybrid male sterility. We also test whether D. persimilis and the USA subspecies of D. pseudoobscura share an allele conferring hybrid sterility in a D. pseudoobscura bogotana genetic background. All loci conferred similar hybrid sterility effects across all strains studied, although we detected some statistically significant quantitative effect variation among D. persimilis alleles of some hybrid incompatibility QTLs. We also detected allelism between D. persimilis and D. pseudoobscura USA at a second chromosome hybrid sterility QTL. We hypothesize that either the QTL is ancestral in D. persimilis and D. pseudoobscura USA and lost in D. pseudoobscura bogotana, or gene flow transferred the QTL from D. persimilis to D. pseudoobscura USA. We discuss our findings in the context of population features that may contribute to variation in hybrid incompatibilities.     

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.     

A test of the chromosomal rearrangement model of speciation in Drosophila pseudoobscura

Recent studies suggest that chromosomal rearrangements play a significant role in speciation by preventing recombination and maintaining species persistence despite interspecies gene flow. Factors conferring adaptation or reproductive isolation are maintained in rearranged regions in the face of hybridization, while such factors are eliminated from collinear regions. As a direct test of this rearrangement model, we evaluated the genetic basis of hybrid male sterility in a sympatric species pair, Drosophila pseudoobscura pseudoobscura and D. persimilis, and an allopatric species pair, D. pseudoobscura bogotana and D. persimilis. Our results are consistent with the proposed model: virtually all of the sterility factors in the former pair are associated with three inverted regions, whereas sterility factors are present in the collinear regions in the latter pair. These findings indicate recombination and selection may have eliminated sterility factors outside the inverted regions between D. p. pseudoobscura and D. persimilis, suggesting chromosomal rearrangements may facilitate species persistence despite hybridization.     

Variation in Activities of Amylase Allozymes Associated with Chromosome Inversions in DROSOPHILA PSEUDOOBSCURA, D. PERSIMILIS and D. MIRANDA

Different electrophoretic alleles of amylase show associations with particular chromosome 3 inversions in D. pseudoobscura and D. persimilis. Relative adult amylase activities were compared in 37, 37 and 10 strains of D. pseudoobscura, D. persimilis and D. miranda, respectively. Strains carrying the same electrophoretic allele were compared by crossing these lines individually to a reference strain carrying a different electrophoretic mobility allele. This procedure allows comparisons among species, inversions, electromorphs and strains for genetic variation in amylase activity. F2 analysis established that the activity variation co-segregates with the structural amylase locus. This type of variation could be due to either structural gene differences or differences in closely linked, cis-acting regulatory regions. Variation has been detected among and within electrophoretic mobility classes. Moreover, this variation is clearly nonrandom and reveals more of the genetic structure associated with the chromosomal inversion phylogeny of D. pseudoobscura and D. persimilis. ----Some of the findings are: (1) Similar electromorphs in D. pseudoobscura and D. persimilis usually show different activities. These species show nearly complete differentiation of amylase alleles, based on activities. (2) D. persimilis has the broadest range of variation in amylase activity, about four-fold between the highest and lowest alleles. D. pseudoobscura and D. miranda are also polymorphic for activity, but have more constrained ranges of variation. D. miranda alleles show on the average about four times the activity of D. pseudoobscura alleles. (3) Some association of electrophoretic mobility and activity has been found. Alleles 1.09 of D. persimilis, as well as 1.43 and 1.55 of D. miranda, have relatively high activity. It may be that these high activity alleles are part of an adaptation to cooler habitats. (4) Within electrophoretic classes, associations of activities with inversions have been found. These are especially strong in D. persimilis. The 1.00 alleles in the ST, KL, MD and WT inversions, the 0.92 allele in the ST and MD inversions and the 1.09 allele in the WT and KL inversions have levels of activities that depend upon the arrangement in which they are located. These results demonstrate that suppression of recombination in inversion heterokaryotypes can result in extensive genic divergence between inversions.     

Genomic impacts of chromosomal inversions in parapatric Drosophila species

Chromosomal inversions impact genetic variation and facilitate speciation in part by reducing recombination in heterokaryotypes. We generated multiple whole-genome shotgun sequences of the parapatric species pair Drosophila pseudoobscura and Drosophila persimilis and their sympatric outgroup (Drosophila miranda) and compared the average pairwise differences for neutral sites within, just outside and far outside of the three large inversions. Divergence between D. pseudoobscura and D. persimilis is high inside the inversions and in the suppressed recombination regions extending 2.5 Mb outside of inversions, but significantly lower in collinear regions further from the inversions. We observe little evidence of decreased divergence predicted to exist in the centre of inversions, suggesting that gene flow through double crossovers or gene conversion is limited within the inversion, or selection is acting within the inversion to maintain divergence in the face of gene flow. In combination with past studies, we provide evidence that inversions in this system maintain areas of high divergence in the face of hybridization, and have done so for a substantial period of time. The left arm of the X chromosome and chromosome 2 inversions appear to have arisen in the lineage leading to D. persimilis approximately 2 Ma, near the time of the split of D. persimilis-D. pseudoobscura-D. miranda, but likely fixed within D. persimilis much more recently, as diversity within D. persimilis is substantially reduced inside and near these two inversions. We also hypothesize that the inversions in D. persimilis may provide an empirical example of the 'mixed geographical mode' theory of inversion origin and fixation, whereby allopatry and secondary contact both play a role.     

The genetics of hybrid male sterility between the allopatric species pair Drosophila persimilis and D. pseudoobscura bogotana: dominant sterility alleles in collinear autosomal regions

F(1) hybrid male sterility is thought to result from interactions between loci on the X chromosome and dominant-acting loci on the autosomes. While X-linked loci that contribute to hybrid male sterility have been precisely localized in many animal taxa, their dominant autosomal interactors have been more difficult to localize precisely and/or have been shown to be of relatively smaller effect. Here, we identified and mapped at least four dominant autosomal factors contributing to hybrid male sterility in the allopatric species pair Drosophila persimilis and D. pseudoobscura bogotana. Using these results, we tested predictions of reduced recombination models of speciation. Consistent with these models, three of the four QTL associated with hybrid male sterility occur in collinear (uninverted) regions of these genomes. Furthermore, these QTL do not contribute significantly to hybrid male sterility in crosses between the sympatric species D. persimilis and D. pseudoobscura pseudoobscura. The autosomal loci identified in this study provide the basis for introgression mapping and, ultimately, for molecular cloning of interacting genes that contribute to F(1) hybrid sterility.     

Molecular population genetics of sequence length diversity in the Adh region of Drosophila pseudoobscura

Positive and negative selection on indel variation may explain the correlation between intron length and recombination levels in natural populations of Drosophila. A nucleotide sequence analysis of the 3.5 kilobase sequence of the alcohol dehydrogenase (Adh) region from 139 Drosophila pseudoobscura strains and one D. miranda strain was used to determine whether positive or negative selection acts on indel variation in a gene that experiences high levels of recombination. A total of 30 deletion and 36 insertion polymorphisms were segregating within D. pseudoobscura populations and no indels were fixed between D. pseudoobscura and its two sibling species D. miranda and D. persimilis. The ratio of Tajima's D to its theoretical minimum value (D(min)) was proposed as a metric to assess the heterogeneity in D among D. pseudoobscura loci when the number of segregating sites differs among loci. The magnitude of the D/D(min) ratio was found to increase as the rate of population expansion increases, allowing one to assess which loci have an excess of rare variants due to population expansion versus purifying selection. D. pseudoobscura populations appear to have had modest increases in size accounting for some of the observed excess of rare variants. The D/D(min) ratio rejected a neutral model for deletion polymorphisms. Linkage disequilibrium among pairs of indels was greater than between pairs of segregating nucleotides. These results suggest that purifying selection removes deletion variation from intron sequences, but not insertion polymorphisms. Genome rearrangement and size-dependent intron evolution are proposed as mechanisms that limit runaway intron expansion.     

Gene Flow and Natural Selection in the Origin of Drosophila Pseudoobscura and Close Relatives

The divergence of Drosophila pseudoobscura and close relatives D. persimilis and D. pseudoobscura bogotana has been studied using comparative DNA sequence data from multiple nuclear loci. New data from the Hsp82 and Adh regions, in conjunction with existing data from Adh and the Period locus, are examined in the light of various models of speciation. The principal finding is that the three loci present very different histories, with Adh indicating large amounts of recent gene flow among the taxa, while little or no gene flow is apparent in the data from the other loci. The data were compared with predictions from several isolation models of divergence. These models include no gene flow, and they were found to be incompatible with the data. Instead the DNA data, taken together with other evidence, seem consistent with divergence models in which natural selection acts against gene flow at some loci more than at others. This family of models includes some sympatric and parapatric speciation models, as well as models of secondary contact and subsequent reinforcement of sexual isolation.     

Mitochondrial DNA variation in Drosophila pseudoobscura and related species in Pacific northwest populations

We have analysed mitochondrial DNA (mtDNA) from Pacific Northwest populations of Drosophila pseudoobscura, D. persimilis, and D. miranda using six restriction enzymes. We find that HpaII restriction sites are hypervariable compared to the other enzymes used. This hypervariability allows construction of a maximum parsimony map linking each mtDNA genotype. Small insertions, possibly tandem duplications, appear to have arisen concomitantly with, or subsequent to, speciation events, perhaps within the A + T rich region. Convergence of mtDNA genotypes is also evident. Unlike findings for other populations of these species, we find little evidence of mitochondrial introgression between D. pseudoobscura and D. persimilis, despite their ability to produce fertile hybrid females.     

A survey of chromosomal and nucleotide sequence variation in Drosophila miranda

There have recently been several studies of the evolution of Y chromosome degeneration and dosage compensation using the neo-sex chromosomes of Drosophila miranda as a model system. To understand these evolutionary processes more fully, it is necessary to document the general pattern of genetic variation in this species. Here we report a survey of chromosomal variation, as well as polymorphism and divergence data, for 12 nuclear genes of D. miranda. These genes exhibit varying levels of DNA sequence polymorphism. Compared to its well-studied sibling species D. pseudoobscura, D. miranda has much less nucleotide sequence variation, and the effective population size of this species is inferred to be several-fold lower. Nevertheless, it harbors a few inversion polymorphisms, one of which involves the neo-X chromosome. There is no convincing evidence for a recent population expansion in D. miranda, in contrast to D. pseudoobscura. The pattern of population subdivision previously observed for the X-linked gene period is not seen for the other loci, suggesting that there is no general population subdivision in D. miranda. However, data on an additional region of period confirm population subdivision for this gene, suggesting that local selection is operating at or near period to promote differentiation between populations.     

Unusual pattern of single nucleotide polymorphism at the exuperantia2 locus of Drosophila pseudoobscura

We have investigated the pattern of DNA sequence variation at the exuperantia2 locus in Drosophila pseudoobscura. This adds to the increasing dataset of genetic variation in D. pseudoobscura, a useful model species for evolutionary genetic studies. The level of silent site nucleotide diversity and the divergence from an outgroup Drosophila miranda are comparable with those for other X-linked loci. One peculiar pattern at the exu2 locus of D. pseudoobscura is a complete linkage disequilibrium between two SNPs, one of which is a replacement site. As a result, there are two distinct haplotype groups in our dataset. Based upon the comparisons with the outgroup sequences from D. miranda and Drosophila persimilis, we show that the newly derived haplotype group has lower diversity than the ancestral haplotype group. The pattern of protein evolution at exu2 shows some deviation from the neutral model. Together, these and other characteristics of the exu2 locus suggest the action of selection on the pattern of SNP variation, consistent with a partial selective sweep associated with the newly derived haplotype.     

Microsatellite variation in populations of Drosophila pseudoobscura and Drosophila persimilis

We have isolated, characterized and mapped 33 dinucleotide, three trinucleotide and one tetranucleotide repeat loci from the four major chromosomes of Drosophila pseudoobscura. Average inferred repeat unit length of the dinucleotide repeats is 12 repeat units, similar to D. melanogaster. Assays of D. pseudoobscura and populations of its sibling species, D. persimilis, using 10 of these loci show extremely high levels of variation compared with similar studies of dinucleotide repeat variation in D. melanogaster populations. The high levels of variation are consistent with an average mutation rate of approximately 10(-6) per locus per generation and an effective population size of D. pseudoobscura approximately four times larger than that of D. melanogaster. Consistent with allozymes and nucleotide sequence polymorphism, the dinucleotide repeat loci reveal minimal structure across four populations of D. pseudoobscura. Finally, our preliminary recombinational mapping of 24 of these microsatellites suggests that the total recombinational genome size may be larger than previously inferred using morphological mutant markers.     

Examination of Allelic Variation at the Hexokinase Loci of DROSOPHILA PSEUDOOBSCURA and D. PERSIMILIS by Different Methods

Recently a number of electrophoretic techniques have been applied to reveal the presence of additional genetic variation among the electrophoretic mobility classes of the highly polymorphic xanthine dehydrogenase (XDH ) and esterase-5 (est-5) loci. We examined the hexokinase loci of Drosophila pseudoobscura and D. persimilis using a variety of techniques to determine whether further allelic variation could be revealed for these much less polymorphic loci and to analyze the nature of the known variation at the hexokinase-1 (hex-1) locus. The following studies were conducted: 135 strains of the two species from six localities were examined with buffer pH ranging from 5.5 to 10.0; 40 strains of D. pseudoobscura and 9 strains of D. persimilis from Mather were studied using starch gel concentrations ranging from 8.5 to 15.5% and were examined for differences in heat stability and reactivity to the thiol reagent pCMSA; strains were also tested for susceptibility to urea denaturation and differences in relative activities. Major findings of the work are: (1) No additional allelic variation could be detected at any of the hexokinase loci by applying these techniques. The finding of abundant hidden genetic variation in XDH and est-5 does not extend to all enzyme loci. (2) Evidence from studies using pCMSA indicates that the hex-1 alleles 0.6, 0.8, 1.0 and 1.2 of the two species form a series of unit charge steps. Since the 0.94 allele of D. persimilis has mobility intermediate between 0.8 and 1.0, it is argued that routine electrophoretic techniques are sensitive to at least some conservative amino acid substitutions. (3) Strong correlations were found at the hex-1 locus between low allelic frequency, reduced relative activity and reduced stability to heat and urea denaturation. Since the three sibling species, D. pseudoobscura, D. persimilis and D. miranda, all appear to share a common high frequency allele (1.0) at that locus, these findings are taken as evidence that the observed allelic frequencies are a result of directional selection and mutation, rather than any form of balancing selection.     

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.     

Molecular evolution of the Sex-Ratio inversion complex in Drosophila pseudoobscura: analysis of the Esterase-5 gene region

The Sex-Ratio chromosome in Drosophila pseudoobscura is subject to meiotic drive. It is associated with a series of three nonoverlapping paracentric inversions on the right arm of the X chromosome. The esterase-5 gene region has been localized to section 23 within the subbasal inversion of the Sex-Ratio inversion complex, making esterase- 5 a convenient locus for molecular evolutionary analyses of the Sex- Ratio inversion complex and the associated drive system. A 504-bp fragment of noncoding, intergenic DNA from the esterase-5 gene region was amplified and sequenced from 14 Sex-Ratio and 14 Standard X chromosomes of D. pseudoobscura, and from 9 X chromosomes of its two sibling species, Drosophila persimilis and Drosophila miranda. There is extensive sequence differentiation between the Sex-Ratio and Standard chromosomal types. The common Standard chromosome is highly polymorphic, while, as expected from either the neutral mutation theory or the selective sweep hypothesis, the rarer Sex-Ratio chromosome has much less within-chromosome nucleotide polymorphism. We estimate that the Standard and Sex-Ratio chromosomes in D. pseudoobscura diverged between 700,000 and 1.3 Mya, or at least 2 million generations ago. The clustering of D. pseudoobscura Sex-Ratio chromosomes in a neighbor- joining phylogeny indicates a fairly old, monophyletic origin in this species. It appears from these data that Sex-Ratio genes were present prior to the divergence of D. pseudoobscura and D. persimilis and that both the Standard and Sex-Ratio chromosomes of D. persimilis were derived from the Standard chromosome of D. pseudoobscura after the inversion events that isolated the D. pseudoobscura Sex-Ratio chromosome.      

The causes of phylogenetic conflict in a classic Drosophila species group

Bifurcating phylogenies are frequently used to describe the evolutionary history of groups of related species. However, simple bifurcating models may poorly represent the evolutionary history of species that have been exchanging genes. Here, we show that the history of three well-known closely related species, Drosophila pseudoobscura, D. persimilis and D. p. bogotana, is not well represented by a bifurcating phylogenetic tree. The phylogenetic relationships among these species vary widely between different genomic regions. Much of this phylogenetic variation can be explained by the potential of different genomic regions to introgress between species, as measured in laboratory studies. We argue that the utility of multiple markers in species-level phylogenetic studies can be greatly enhanced by knowledge of genomic location and, in the case of hybridizing species, by knowledge of the functional or linkage relationships among the markers and regions of the genome that reduce hybrid fitness.     

A recombinational portrait of the Drosophila pseudoobscura genome

Drosophila pseudoobscura has been intensively studied by evolutionary biologists for over 70 years. The recent publication of the genome sequence not only permits studies of comparative genomics with other dipterans but also opens the door to identifying genes associated with adaptive traits or speciation or testing for the signature of natural selection across the genome. Information on regional rates of recombination, localization of inversion breakpoints distinguishing it from its sibling species D. persimilis, and known polymorphic markers may be useful in such studies. Here, we present a molecular linkage map of four of the five major chromosome arms of D. pseudoobscura. In doing so, we order and orient several sequence contig groups, localize the inversion breakpoints on chromosome 2 to intervals of 200 kilobases, and identify one error in the published sequence assembly. Our results show that regional recombination rates in D. pseudoobscura are much higher than in D. melanogaster and significantly higher even than in D. persimilis. Furthermore, we detect a non-significant positive correlation between recombination rate and published DNA sequence variation. Finally, the online Appendix presents 200 primer sequence pairs for molecular markers that can be used for mapping of quantitative trait loci, of which 125 are known to be polymorphic within or between species.     

EPISTASIS MODIFIES THE DOMINANCE OF LOCI CAUSING HYBRID MALE STERILITY IN THE DROSOPHILA PSEUDOOBSCURA SPECIES GROUP

Speciation, the evolution of reproductive isolation between populations, serves as the driving force for generating biodiversity. Postzygotic barriers to gene flow, such as F ₁ hybrid sterility and inviability, play important roles in the establishment and maintenance of biological species. F ₁ hybrid incompatibilities in taxa that obey Haldane's rule, the observation that the heterogametic sex suffers greater hybrid fitness problems than the homogametic sex, are thought to often result from interactions between recessive-acting X-linked loci and dominant-acting autosomal loci. Because they play such prominent roles in producing hybrid incompatibilities, we examine the dominance and nature of epistasis between alleles derived from Drosophila persimilis that confer hybrid male sterility in the genetic background of its sister species, D. pseudoobscura bogotana . We show that epistasis elevates the apparent dominance of individually recessive-acting QTL such that they can contribute to F ₁ hybrid sterility. These results have important implications for assumptions underlying theoretical models of hybrid incompatibilities and may offer a possible explanation for why, to date, identification of dominant-acting autosomal "speciation genes" has been challenging.     

The role of gene conversion in determining sequence variation and divergence in the Est-5 gene family in Drosophila pseudoobscura.

Nucleotide sequences of eight Est-5A and Est-5C genes corresponding to previously sequenced Est-5B genes in Drosophila pseudoobscura were determined to compare patterns of polymorphism and divergence among members of this small gene family. The three esterase genes were also sequenced from D. persimilis and D. miranda for interspecific comparisons. The data provide evidence that gene conversion between loci contributes to polymorphism and to the homogenization of the Est5 genes. For Est-5B, which encodes one of the most highly polymorphic proteins in Drosophila, 12% of the segregating amino acid variants appear to have been introduced via gene conversion from other members of the gene family. Interlocus gene conversion can also explain high sequence similarity, especially at synonymous sites, between Est-5B and Est-5A. Tests of neutrality using interspecific comparisons show that levels of polymorphism conform to neutral expectations at each Est-5 locus. However, McDonald-Kreitman tests based on intraspecific gene comparisons indicate that positive selection on amino acids has accompanied Est-5 gene duplication and divergence in D. pseudoobscura.