Competition between Mitochondrial Haplotypes in Distinct Nuclear Genetic Environments: Drosophila Pseudoobscura Vs. D. Persimilis

A test for coadaptation of nuclear and mitochondrial genomes was performed using the sibling species, Drosophila pseudoobscura and D. persimilis. Two lines of flies with ``disrupted' cytonuclear genotypes were constructed by repeated backcrossing of males from one species to females carrying mitochondrial DNA (mtDNA) from the other species. Each ``disrupted' strain was competed in population cages with the original stock of each species from which the recurrent males were obtained during the backcrossing. As such, the two species' mitochondrial types were competed reciprocally in the nuclear genetic environments of each species. The trajectories of mtDNA haplotypes were followed in discrete-generation population cages using a PCR-four-cutter approach. A significant increase in the frequency of D. pseudoobscura mtDNA was observed in each of four replicate cages with a D. pseudoobscura nuclear background. In the D. persimilis nuclear background, one cage actually showed an increase in frequency of D. pseudoobscura mtDNA, although together the four replicate cages show little change in frequency. These results were repeated after frequency perturbations and reinitiation of each cage. An analysis of fitness components revealed that fertility selection greatly outweighed viability selection in these cytonuclear competition experiments. The asymmetry of the fitnesses of the mtDNA haplotypes on the two genetic backgrounds is consistent in direction with the previously reported asymmetry of female fertility in backcrosses between these two species. While our experiments do not allow us to identify mtDNA as the sole source of fitness variation, at a minimum the data indicate a fitness association between nuclear fertility factors and the D. pseudoobscura mtDNA on its own genetic background.     

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.     

Postmating-prezygotic isolation is not an important source of selection for reinforcement within and between species in Drosophila pseudoobscura and D. persimilis

Abstract Most work on adaptive speciation to date has focused on the role of low hybrid fitness as the force driving reinforcement (the evolution of premating isolation after secondary contact that reduces the likelihood of matings between populations). However, recent theoretical work has shown that postmating, prezygotic incompatibilities may also be important in driving premating isolation. We quantified premating, postmating-prezygotic, and early postzygotic fitness effects in crosses among three populations: Drosophila persimilis, D. pseudoobscura USA (sympatric to D. persimilis ), and D. pseudoobscura Bogotá (allopatric to D. persimilis ). Interspecific matings were more likely to fail when they involved the sympatric populations than when they involved the allopatric populations, consistent with reinforcement. We also found that failure rate in sympatric mating trials depended on whether D. persimilis females were paired with D. pseudoobscura males or the reverse. This asymmetry most likely indicates differences in discrimination against heterospecific males by females. By measuring egg laying rate, fertilization success and hatching success, we also compared components of postmating-prezygotic and early postzygotic isolation. Postmating-prezygotic fitness costs were small and not distinguishable between hetero- and conspecific crosses. Early postzygotic fitness effects due to hatching success differences were also small in between-population crosses. There was, however, a postzygotic fitness effect that may have resulted from an X-linked allele found in one of the two strains of D. pseudoobscura USA. We conclude that the postmating-prezygotic fitness costs we measured probably did not drive premating isolation in these species. Premating isolation is most likely driven in sympatric populations by previously known hybrid male sterility.     

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.     

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.     

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.     

上海及邻近地区果蝇（Drosophila albomicans）的迁入及其线粒体DNA多态性研究

本文报道１９９１年秋果蝇（Ｄ．ａｌｂｏｍｉｃａｎｓ）突然迁入上海的事实。我们用９种限制性核酸内切酶对上海、嘉定、青浦和南汇４个群体的２９个单雌系的线粒体ＤＮＡ进行了酶切片段长度多态性研究,发现每个单雌系各有其特征的酶切图谱。在与广东、昆明、台湾珊瑚潭、日本洲本及马来西亚吉隆坡群体比较后,我们认为上海及邻近地区突然发现的Ｄ．ａｌｂｏｍｉｃａｎｓ来自中国大陆南部的多个群体。这一结论支持该种目前正处于向北迁移的分化阶段的观点。我们还对线粒体ＤＮＡ的种内多态性进行了讨论     

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.     

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.     

Estimates of Gene Flow in Drosophila Pseudoobscura Determined from Nucleotide Sequence Analysis of the Alcohol Dehydrogenase Region

The genetic structure of Drosophila pseudoobscura populations was inferred from a nucleotide sequence analysis of a 3.4-kb segment of the alcohol dehydrogenase (Adh) region. A total of 99 isochromosomal strains collected from 13 populations in North and South America were used to determine if any population departed from a neutral model and to estimate levels of gene flow between populations. This study also included the nucleotide sequences from two sibling species, D. persimilis and D. miranda. We estimated the neutral mutation parameter, 4N mu, in synonymous and noncoding sites for 17 subregions of Adh in each of nine populations with sample sizes greater than three. The nucleotide diversity data in the nine populations was tested for departures from an equilibrium neutral model with two statistical tests. The Tajima and the Hudson, Kreitman, Aguade tests showed that each population fails to reject a neutral model. Tests for genetic differentiation between populations fail to show any population substructure among the North American populations of D. pseudoobscura. The nucleotide diversity data is consistent with direct and indirect measures of gene flow that show extensive dispersal between populations of D. pseudoobscura.     

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.     

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.     

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.      

Cytoplasmic Incompatibility and Mating Preference in Colombian Drosophila Pseudoobscura

MACRAE and ANDERSON observed a large frequency change of mitochondrial DNA (mtDNA) haplotypes in a population initiated with two allopatric strains of Drosophila pseudoobscura, BogER from Colombia and AH162 from California. They concluded that mtDNA haplotypes in D. pseudoobscura are not always selectively neutral. NIGRO and PROUT suggested, however, that a maternally transmitted incompatibility system, similar to the one they observed in two strains of D. simulans from Italy, could account for the observed mtDNA frequency changes. SINGH and HALE postulated that a mating preference between the strains BogER and AH162 in MACRAE and ANDERSON's experiment, in the form of negative assortative mating, could also account for the mtDNA frequency changes. We report two experiments designed to test the hypotheses: that a maternally transmitted cytoplasmic incompatibility system exists between D. pseudoobscura strains BogER and AH162; and, that BogER females mate preferentially with AH162 males. Our results do not support either hypothesis.     

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.     

Is There Selection on RFLP Differences in Mitochondrial Dna?

Experimental populations of Drosophila simulans were established for the purpose of detecting the presence or absence of selection on a restriction fragment length polymorphism in mitochondrial DNA (mtDNA). It was then discovered that the founding strains differed with respect to the Rickettsia-mediated incompatibility system in this species, which is maternally transmitted together with the mtDNA differences. A population model was constructed using the known fitness effects of the incompatibility system, with the result that the population trajectories can be completely explained by the effects of the microorganism with no need to invoke selection on mtDNA. The strong conclusion is that in this case we can rule out the strong selection proposed by MacRae and Anderson to explain the "dramatic mtDNA changes" in their Drosophila pseudoobscura populations. The population theory used for the experiments is discussed in the context of natural populations. Estimated parameters include the possibility that with two populations, one with the organism and one without it, there may be no bias as to which will invade the other, which in turn suggests no global tendency for the infection to spread or decline.     

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.     

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.