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

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.     

Frequency Distribution of Esterase-5 Alleles in Two Populations of DROSOPHILA PSEUDOOBSCURA

Statistical tests comparing allele frequencies in natural populations with those predicted by various theories of genic variation depend critically on the accurate enumeration of alleles. This study used a series of five sequential electrophoretic conditions to characterize the allele frequency distributions of esterase-5 in two large population samples of Drosophila pseudoobscura from California. In Standard chromosome lines 12 electromorphs were discriminated using a single electrophoretic condition. When four additional criteria were used, the number of electromorphs increased to 41, 33 in one population and 22 in the other. Both populations had the same two alleles in high frequency, with other alleles present in frequencies of 6% or less. Although each population had a number of unique alleles, a χ² contingency test demonstrated no significant genetic divergence between them. A statistical comparison of allele frequencies in both populations with that predicted by neutral models suggests that the individual and combined distributions deviate from neutrality in the direction of purifying selection.—Sex-Ratio chromosomes differed markedly from Standard chromosomes in both allelic content and diversity. In 32 Sex-Ratio chromosomes from one population only three alleles were found, all of which were detected under the initial "standard" electrophoretic conditions. Moreover, none of these alleles was found in the Standard chromosome lines.     

Rates and patterns of chromosomal evolution in Drosophila pseudoobscura and D. miranda

Comparisons of gene orders between species permit estimation of the rate of chromosomal evolution since their divergence from a common ancestor. We have compared gene orders on three chromosomes of Drosophila pseudoobscura with its close relative, D. miranda, and the distant outgroup species, D. melanogaster, by using the public genome sequences of D. pseudoobscura and D. melanogaster and approximately 50 in situ hybridizations of gene probes in D. miranda. We find no evidence for extensive transfer of genes among chromosomes in D. miranda. The rates of chromosomal rearrangements between D. miranda and D. pseudoobscura are far higher than those found before in Drosophila and approach those for nematodes, the fastest rates among higher eukaryotes. In addition, we find that the D. pseudoobscura chromosome with the highest level of inversion polymorphism (Muller's element C) does not show an unusually fast rate of evolution with respect to chromosome structure, suggesting that this classic case of inversion polymorphism reflects selection rather than mutational processes. On the basis of our results, we propose possible ancestral arrangements for the D. pseudoobscura C chromosome, which are different from those in the current literature. We also describe a new method for correcting for rearrangements that are not detected with a limited set of markers.     

Patterns of evolution of genes disrupted in expression in Drosophila species hybrids

Divergence between species in regulatory pathways may contribute to hybrid incompatibilities such as sterility. Consistent with this idea, genes involved in male fertility often evolve faster than most other genes both in amino acid sequence and in expression. Previously, we identified a panel of male-specific genes under-expressed in sterile male hybrids of Drosophila simulans and D. mauritiana relative to pure species, and we showed that this under-expression is associated with infertility. In a preliminary effort to assess the generalities in the patterns of evolution of these genes, I examined patterns of mRNA expression in three of these genes in sterile F 1 hybrid males of D. pseudoobscura and D. persimilis . F 1 hybrid males bearing D. persimilis X chromosomes under-expressed all these genes relative to the parental species, while hybrids bearing D. pseudoobscura X chromosomes under-expressed two of these three genes. Interestingly, the third gene, CG5762 , has undergone extensive amino acid evolution within the D. pseudoobscura species group, possibly driven by positive natural selection. We conclude that some of the same genes exhibit disruptions in expression within each of the two species groups, which could suggest commonalities in the regulatory architecture of sterility in these groups. Alternative explanations are also considered.     

The genetics of reproductive isolation and the potential for gene exchange between Drosophila pseudoobscura and D. persimilis via backcross hybrid males

Hybrid male sterility, hybrid inviability, sexual isolation, and a hybrid male courtship dysfunction reproductively isolate Drosophila pseudoobscura and D. persimilis. Previous studies of the genetic bases of these isolating mechanisms have yielded only limited information about how much and what areas of the genome are susceptible to interspecies introgression. We have examined the genetic basis of these barriers to gene exchange in several thousand backcross hybrid male progeny of these species using 14 codominant molecular genetic markers spanning the five chromosomes of these species, focusing particularly on the autosomes. Hybrid male sterility, hybrid inviability, and the hybrid male courtship dysfunction were all associated with X-autosome interactions involving primarily the inverted regions on the left arm of the X-chromosome and the center of the second chromosome. Sexual isolation from D. pseudoobscura females was primarily associated with the left arm of the X-chromosome, although both the right arm and the center of the second chromosome also contributed to it. Sexual isolation from D. persimilis females was primarily associated with the second chromosome. The absence of isolating mechanisms being associated with many autosomal regions, including some large inverted regions that separate the strains, suggests that these phenotypes may not be caused by genes spread throughout the genome. We suggest that gene flow between these species via hybrid males may be possible at loci spread across much of the autosomes.     

Evidence for Extensive Genetic Differentiation between the Sex-Ratio and the Standard Arrangement of DROSOPHILA PSEUDOOBSCURA and D. PERSIMILIS and Identification of Hybrid Sterility Factors

This study deals with sex-ratio genes tightly linked within the Sex-Ratio inversion. By taking advantage of the fact that the Sex-Ratio chromosome of Drosophila persimilis [SR(B)] is homosequential to the Standard chromosome of D. pseudoobscura [ST(A)], we carried out two reciprocal introgression experiments. Individual segments of SR(B) or ST(A) were introgressed into the genome of D. pseudoobscura or D. persimilis, respectively. Males possessing a hybrid SR(B)-ST(A) X chromosome and a genetic background derived from either of the two species were tested for fertility and sex-ratio expression.—It was found that, in terms of the meiotic drive genes, the Sex-Ratio chromosome differs extensively from the Standard chromosome. Because recombinations of these genes result in a complete loss of sex-ratio expression, this finding lends strong support to the hypothesis of gene coadaptation. Coadaptation, in this context, is the advantage of being transmitted preferentially. In light of this finding, the evolution of the sex-ratio system in these two sibling species is discussed.—Introgression experiments also yielded information about hybrid sterility. With reciprocal introgression, sterility interactions were found to be "asymmetric." The asymmetry is fully expected from the viewpoint of evolution of postmating reproductive isolation.     

Molecular Organization of the X Chromosome in Different Species of the Obscura Group of Drosophila

Nine single copy regions located on the X chromosome have been mapped by in situ hybridization in six species of the obscura group of Drosophila. Three Palearctic species, D. subobscura, D. madeirensis and D. guanche, and three Nearctic species, D. pseudoobscura, D. persimilis and D. miranda, have been studied. Eight of the regions include known genes from D. melanogaster (Pgd, zeste, white, cut, vermilion, RNA polymerase II 215, forked and suppressor of forked) and the ninth region (lambda DsubF6) has not yet been characterized. In all six species, as in D. melanogaster, all probes hybridize to a single site. Established chromosomal arm homologies of Muller's element A are only partly supported by present results since two of the probes (Pgd and zeste) hybridize at the proximal end of the XR chromosomal arm in the three Nearctic species. In addition to the centric fusion of Muller's A (= XL) and D (= XR) elements, the metacentric X chromosome of the Nearctic species requires a pericentric inversion to account for this result. Previously proposed homologies of particular chromosomal regions of the A (= X) chromosome in the three species of the D. subobscura cluster and of the XL chromosomal arm in the three species of the D. pseudoobscura cluster are discussed in light of the present results. Location of the studied markers has changed drastically not only since the divergence between the melanogaster and obscura groups but also since the Palearctic and Nearctic species of the obscura group diverged.     

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.     

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.     

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.     

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.     

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.     

A re-evaluation of 12S ribosomal RNA variability in Drosophila pseudoobscura

Two recent studies have presented conflicting views on variation present within the 294 base third domain of the 12S rRNA gene in the genus Drosophila, and in D. pseudoobscura in particular. One study suggested that this gene is highly invariant across the genus, while another recovered 22 distinct haplotypes from 22 strains of D. pseudoobscura. We have sequenced this gene in numerous lines of D. pseudoobscura and its relatives, noting only two haplotypes in the third domain, and we failed to confirm any of the published sequences. Second, we note that the published sequence divergence between strains of D. pseudoobscura was as great as that documented between distantly related Drosophila species. Third, we show that the published polymorphisms of this region within D. pseudoobscura would disrupt the secondary structure of the resulting molecule. We conclude that the published 12S rRNA sequences of D. pseudoobscura do not accurately reflect variability of the functional gene, and that this gene is relatively invariant in D. pseudoobscura and D. persimilis.     

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.     

Chromosomal and behavioral studies of Mexican Drosophila. III. Inversion polymorphism of D. pseudoobscura

Four new gene arrangements of chromosome 3 of Drosophila pseudoobscura are reported, as well as an updated phylogenetic sequence of the inversion polymorphism in this species. Evidence that the Tree Line gene arrangement of chromosome 3 was the original form in D. pseudoobscura, and that Central Mexico is the center of distribution of this species is discussed.     

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.     

Localization of ecs, dor and swl genes in eight Drosophila species]

The DNA sequences from Drosophila melanogaster early ecdysterone-inducible puff 2B have been located in 8 Drosophila species by in situ hybridization. The location site of the ecs, dor and swi genes in D. funebris, D. virilis, D. hydei, D. repleta, D. mercatorum, D. paranaensis is a puff on the telomeric and of X chromosome; in D. kanekoi it is the puff in distal part of X chromosome; and in D. pseudoobscura it is the puff in proximal portion of X chromosome. So, conservative organization of DNA sequences located in D. melanogaster 2B puff could be suggested. Dispersed distribution of some DNA segments from the region studied in D. hydei chromosomes was revealed.