Chromosomal rearrangements directly cause underdominant F1 pollen sterility in Mimulus lewisii–Mimulus cardinalis hybrids

Chromosomal rearrangements can contribute to the evolution of postzygotic reproductive isolation directly, by disrupting meiosis in F₁ hybrids, or indirectly, by suppressing recombination among genic incompatibilities. Because direct effects of rearrangements on fertility imply fitness costs during their spread, understanding the mechanism of F₁ hybrid sterility is integral to reconstructing the role(s) of rearrangements in speciation. In hybrids between monkeyflowers Mimulus cardinalis and Mimulus lewisii, rearrangements contain all quantitative trait loci (QTLs) for both premating barriers and pollen sterility, suggesting that they may have facilitated speciation in this model system. We used artificial chromosome doubling and comparative mapping to test whether heterozygous rearrangements directly cause underdominant male sterility in M. lewisii–M. cardinalis hybrids. Consistent with a direct chromosomal basis for hybrid sterility, synthetic tetraploid F₁s showed highly restored fertility (83.4% pollen fertility) relative to diploids F₁s (36.0%). Additional mapping with Mimulus parishii–M. cardinalis and M. parishii–M. lewisii hybrids demonstrated that underdominant male sterility is caused by one M. lewisii specific and one M. cardinalis specific reciprocal translocation, but that inversions had no direct effects on fertility. We discuss the importance of translocations as causes of reproductive isolation, and consider models for how underdominant rearrangements spread and fix despite intrinsic fitness costs.     

Floral volatile alleles can contribute to pollinator‐mediated reproductive isolation in monkeyflowers (Mimulus)

Pollinator‐mediated reproductive isolation is a major factor in driving the diversification of flowering plants. Studies of floral traits involved in reproductive isolation have focused nearly exclusively on visual signals, such as flower color. The role of less obvious signals, such as floral scent, has been studied only recently. In particular, the genetics of floral volatiles involved in mediating differential pollinator visitation remains unknown. The bumblebee‐pollinated Mimulus lewisii and hummingbird‐pollinated Mimulus cardinalis are a model system for studying reproductive isolation via pollinator preference. We have shown that these two species differ in three floral terpenoid volatiles – d ‐limonene, β‐myrcene, and E‐β‐ocimene – that are attractive to bumblebee pollinators. By genetic mapping and in vitro analysis of enzyme activity we demonstrate that these interspecific differences are consistent with allelic variation at two loci, LIMONENE‐MYRCENE SYNTHASE (LMS) and OCIMENE SYNTHASE (OS). Mimulus lewisii LMS (MlLMS) and OS (MlOS) are expressed most strongly in floral tissue in the last stages of floral development. Mimulus cardinalis LMS (McLMS) is weakly expressed and has a nonsense mutation in exon 3. Mimulus cardinalis OS (McOS) is expressed similarly to MlOS, but the encoded McOS enzyme produces no E‐β‐ocimene. Recapitulating the M. cardinalis phenotype by reducing the expression of MlLMS by RNA interference in transgenic M. lewisii produces no behavioral difference in pollinating bumblebees; however, reducing MlOS expression produces a 6% decrease in visitation. Allelic variation at the OCIMENE SYNTHASE locus is likely to contribute to differential pollinator visitation, and thus promote reproductive isolation between M. lewisii and M. cardinalis. OCIMENE SYNTHASE joins a growing list of ‘speciation genes’ (‘barrier genes’) in flowering plants.     

RFLP Mapping in Soybean: Association between Marker Loci and Variation in Quantitative Traits

Quantitative trait loci (QTLs) have been mapped to small intervals along the chromosomes of tomato (Lycopersicon esculentum), by a method we call substitution mapping. The size of the interval to which a QTL can be mapped is determined primarily by the number and spacing of previously mapped genetic markers in the region surrounding the QTL. We demonstrate the method using tomato genotypes carrying chromosomal segments from Lycopersicon chmielewskii, a wild relative of tomato with high soluble solids concentration but small fruit and low yield. Different L. chmielewskii chromosomal segments carrying a common restriction fragment length polymorphism were identified, and their regions of overlap determined using all available genetic markers. The effect of these chromosomal segments on soluble solids concentration, fruit mass, yield, and pH, was determined in the field. Many overlapping chromosomal segments had very different phenotypic effects, indicating QTLs affecting the phenotype(s) to lie in intervals of as little as 3 cM by which the segments differed. Some associations between different traits were attributed to close linkage between two or more QTLs, rather than pleiotropic effects of a single QTL: in such cases, recombination should separate desirable QTLs from genes with undesirable effects. The prominence of such trait associations in wide crosses appears partly due to infrequent reciprocal recombination between heterozygous chromosomal segments flanked by homozygous regions. Substitution mapping is particularly applicable to gene introgression from wild to domestic species, and generally useful in narrowing the gap between linkage mapping and physical mapping of QTLs.     

EVIDENCE FOR DOBZHANSKY-MULLER INCOMPATIBILITES CONTRIBUTING TO THE STERILITY OF HYBRIDS BETWEEN MIMULUS GUTTATUS AND M. NASUTUS

Abstract Both chromosomal rearrangements and negative interactions among loci (Dobzhansky‐Muller incompatibilities) have been advanced as the genetic mechanism underlying the sterility of interspecific hybrids. These alternatives invoke very different evolutionary histories during speciation and also predict different patterns of sterility in artificial hybrids. Chromosomal rearrangements require drift, inbreeding, or other special conditions for initial fixation and, because heterozygosity per se generates any problems with gamete formation, F1 hybrids will be most infertile. In contrast, Dobzhansky‐Muller incompatibilities may arise as byproducts of adaptive evolution and often affect the segregating F2 generation most severely. To distinguish the effects of these two mechanisms early in divergence, we investigated the quantitative genetics of hybrid sterility in a line cross between two members of the Mimulus guttatus species complex (M. guttatus and M. nasutus). Hybrids showed partial male and female sterility, and the patterns of infertility were not consistent with the action of chromosomal rearrangements alone. F2 and F1 hybrids exhibited equal decreases in pollen viability (> 40%) relative to the highly fertile parental lines. A large excess of completely pollen‐sterile F2 genotypes also pointed to the segregation of Dobzhansky‐Muller incompatibility factors affecting male fertility. Female fertility showed a pattern similarly consistent with epistatic interactions: F2 hybrids produced far fewer seeds per flower than F1 hybrids (88.0 ± 2.8 vs. 162.9 ± 8.5 SE, respectively) and either parental line, and many F2 genotypes were completely female sterile. Dobzhansky‐Muller interactions also resulted in the breakdown of several nonreproductive characters and appear to contribute to correlations between male and female fertility in the F2 generation. These results parallel and contrast with the genetics of postzygotic isolation in model animal systems and are a first step toward understanding the process of speciation in this well‐studied group of flowering plants.     

Determinants of leaf temperature in California Mimulus species at different altitudes

Summary Leaf energy balance and gas-exchange characteristics were studied in Mimulus cardinalis at 400 m and Mimulus lewisii at 2,700 m in the Sierra Nevada of central California. In contrast to previous observations, leaf temperatures were not near 30° C at air temperatures from 20 to 40° C but were coupled quite closely to air temperature. Stomatal conductance in both species decreased in response to increases in the water vapor concentration gradient, a response opposite that required to establish 30°C leaf temperatures over a wide range of air temperatures. The temperature optima for photosynthesis were broad in both species but 5° C higher for M. cardinalis than for M. lewisii. The direct or indirect effects of altitude did not contribute significantly to the maintenance of constant leaf temperatures. For both species, maintaining constant leaf temperatures appears to be less important than avoiding inhibitory water stress or diffusion limitation of photosynthesis.     

Growth and leaf physiology of monkeyflowers with different altitude ranges

Every species is limited both geographically and ecologically to a subset of available habitats, yet for many species the causes of distribution limits are unknown. Temperature is thought to be one of the primary determinants of species distributions along latitudinal and altitudinal gradients. This study examined leaf physiology and plant performance under contrasting temperature regimes of sister species of monkeyflower, Mimulus cardinalis and Mimulus lewisii (Phrymaceae), that differ in altitude distribution to test the hypothesis that temperature is the primary determinant of differences in fitness versus altitude. Each species attained greatest aboveground biomass, net photosynthetic rate, and effective quantum yield of photosystem II when grown under temperatures characteristic of the altitudinal range center. Although both species exhibited greater stem length, stomatal conductance, and intercellular CO₂ concentration in hot than in cold temperatures, these traits showed much greater reductions under cold temperature for M. cardinalis than for M. lewisii. Survival of M. lewisii was also sensitive to temperature, showing a striking decrease in hot temperatures. Within each temperature regime, the species native to that temperature displayed greatest growth and leaf physiological capacity. Populations from the altitude range center and range margin of each species were used to examine population differentiation, but central and marginal populations did not differ in most growth or leaf physiological responses to temperature. This study provides evidence that M. cardinalis and M. lewisii differ in survival, growth, and leaf physiology under temperature regimes characterizing their contrasting low and high altitude range centers, and suggests that the species’ altitude range limits may arise, in part, due to metabolic limitations on growth that ultimately decrease survival and limit reproduction.Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Ancient and recent introgression shape the evolutionary history of pollinator adaptation and speciation in a model monkeyflower radiation (Mimulus section Erythranthe)

Inferences about past processes of adaptation and speciation require a gene-scale and genome-wide understanding of the evolutionary history of diverging taxa. In this study, we use genome-wide capture of nuclear gene sequences, plus skimming of organellar sequences, to investigate the phylogenomics of monkeyflowers in Mimulus section Erythranthe (27 accessions from seven species). Taxa within Erythranthe, particularly the parapatric and putatively sister species M. lewisii (bee-pollinated) and M. cardinalis (hummingbird-pollinated), have been a model system for investigating the ecological genetics of speciation and adaptation for over five decades. Across >8000 nuclear loci, multiple methods resolve a predominant species tree in which M. cardinalis groups with other hummingbird-pollinated taxa (37% of gene trees), rather than being sister to M. lewisii (32% of gene trees). We independently corroborate a single evolution of hummingbird pollination syndrome in Erythranthe by demonstrating functional redundancy in genetic complementation tests of floral traits in hybrids; together, these analyses overturn a textbook case of pollination-syndrome convergence. Strong asymmetries in allele sharing (Patterson’s D-statistic and related tests) indicate that gene tree discordance reflects ancient and recent introgression rather than incomplete lineage sorting. Consistent with abundant introgression blurring the history of divergence, low-recombination and adaptation-associated regions support the new species tree, while high-recombination regions generate phylogenetic evidence for sister status for M. lewisii and M. cardinalis. Population-level sampling of core taxa also revealed two instances of chloroplast capture, with Sierran M. lewisii and Southern Californian M. parishii each carrying organelle genomes nested within respective sympatric M. cardinalis clades. A recent organellar transfer from M. cardinalis, an outcrosser where selfish cytonuclear dynamics are more likely, may account for the unexpected cytoplasmic male sterility effects of selfer M. parishii organelles in hybrids with M. lewisii. Overall, our phylogenomic results reveal extensive reticulation throughout the evolutionary history of a classic monkeyflower radiation, suggesting that natural selection (re-)assembles and maintains species-diagnostic traits and barriers in the face of gene flow. Our findings further underline the challenges, even in reproductively isolated species, in distinguishing re-use of adaptive alleles from true convergence and emphasize the value of a phylogenomic framework for reconstructing the evolutionary genetics of adaptation and speciation.     

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.     

UNDERSTANDING THE BASIS OF DIMINISHED GENE FLOW BETWEEN HYBRIDIZING CHROMOSOME RACES OF THE HOUSE MOUSE

Speciation may be promoted in hybrid zones if there is an interruption to gene flow between the hybridizing forms. For hybridizing chromosome races of the house mouse in Valtellina (Italy), distinguished by whole‐arm chromosomal rearrangements, previous studies have shown that there is greater interruption to gene flow at the centromeres of chromosomes that differ between the races than at distal regions of the same chromosome or at the centromeres of other chromosomes. Here, by increasing the number of markers along race‐specific chromosomes, we reveal a decay in between‐race genetic differentiation from the centromere to the distal telomere. For the first time, we use simulation models to investigate the possible role of recombination suppression and hybrid breakdown in generating this pattern. We also consider epistasis and selective sweeps as explanations for isolated chromosomal regions away from the centromere showing differentiation between the races. Hybrid breakdown alone is the simplest explanation for the decay in genetic differentiation with distance from the centromere. Robertsonian fusions/whole‐arm reciprocal translocations are common chromosomal rearrangements characterizing both closely related species and races within species, and this fine‐scale empirical analysis suggests that the unfitness associated with these rearrangements in the heterozygous state may contribute to the speciation process.     

No evidence for Z‐chromosome rearrangements between the pied flycatcher and the collared flycatcher as judged by gene‐based comparative genetic maps

Revealing the genetic basis of reproductive isolation is fundamental for understanding the speciation process. Chromosome speciation models propose a role for chromosomal rearrangements in promoting the build up of reproductive isolation between diverging populations and empirical data from several animal and plant taxa support these models. The pied flycatcher and the collared flycatcher are two closely related species that probably evolved reproductive isolation during geographical separation in Pleistocene glaciation refugia. Despite the short divergence time and current hybridization, these two species demonstrate a high degree of intrinsic post‐zygotic isolation and previous studies have shown that traits involved in mate choice and hybrid viability map to the Z‐chromosome. Could rearrangements of the Z‐chromosome between the species explain their reproductive isolation? We developed high coverage Z‐chromosome linkage maps for both species, using gene‐based markers and large‐scale SNP genotyping. Best order maps contained 57–62 gene markers with an estimated average density of one every 1–1.5 Mb. We estimated the recombination rates in flycatcher Z‐chromosomes to 1.1–1.3 cM/Mb. A comparison of the maps of the two species revealed extensive co‐linearity with no strong evidence for chromosomal rearrangements. This study does therefore not provide support the idea that sex chromosome rearrangements have caused the relatively strong post‐zygotic reproductive isolation between these two Ficedula species.     

Dissecting the role of a large chromosomal inversion in life history divergence throughout the Mimulus guttatus species complex

Chromosomal inversions can play an important role in adaptation, but the mechanism of their action in many natural populations remains unclear. An inversion could suppress recombination between locally beneficial alleles, thereby preventing maladaptive reshuffling with less‐fit, migrant alleles. The recombination suppression hypothesis has gained much theoretical support but empirical tests are lacking. Here, we evaluated the evolutionary history and phenotypic effects of a chromosomal inversion which differentiates annual and perennial forms of Mimulus guttatus. We found that perennials likely possess the derived orientation of the inversion. In addition, this perennial orientation occurs in a second perennial species, M. decorus, where it is strongly associated with life history differences between co‐occurring M. decorus and annual M. guttatus. One prediction of the recombination suppression hypothesis is that loci contributing to local adaptation will predate the inversion. To test whether the loci influencing perenniality pre‐date this inversion, we mapped QTLs for life history traits that differ between annual M. guttatus and a more distantly related, collinear perennial species, M. tilingii. Consistent with the recombination suppression hypothesis, we found that this region is associated with life history in the absence of the inversion, and this association can be broken into at least two QTLs. However, the absolute phenotypic effect of the LG8 inversion region on life history is weaker in M. tilingii than in perennials which possess the inversion. Thus, while we find support for the recombination suppression hypothesis, the contribution of this inversion to life history divergence in this group is likely complex.     

Temperature and photoperiodic control of developmental responses in climatic races of Mimulus

Two strongly differentiated climatic races of the Mimulus cardinalis-lewisiicomplex were grown at a variety of temperatures (3–27°C)and photoperiods (8 and 16 hr) under controlled environmentalconditions. M. cardinalis (the lowland race, 400 m) and M. lewisii(the sub-alpine race, 3200 m) were found to differ in theirphysiological responses to the varied environments in severalsignificant ways: 1) At 27°C (16 hrphotoperiod), M. lewisiisustained 100% mortality in contrast to the substantial growthand flowering of M. cardinalis under these conditions; 2) In8 hr photoperiods at all temperatures, there was little growthand no flowering in M. lewisii whereas there was considerablegrowth at all temperatures, and flowering at 23 and 27°Cin M. cardinalis; 3) At low temperatures (7–15°C),16 hr photoperiods, flowering occurred a week or two earlierin M. lewisii than in M. cardinalis. The lowland race has asignificantly wider temperature and photoperiodic tolerancethan has the sub-alpine race. Applications of gibberellic acidto rosette Mimulus plants under non-inductive conditions (15°C,8 hr photoperiod) promoted vigorous stem elongation withoutflowering. The application of steroids, other hormones and metaboliteshad no observable effects. ¹Present address: Faculty of Botany, Ohio State University,Columbus, Ohio, U.S.A. (Received March 16, 1970; )     

Putative chromosomal rearrangements are associated primarily with ecotype divergence rather than geographic separation in an intertidal,poorly dispersing snail

Littorina saxatilis is becoming a model system for understanding the genomic basis of ecological speciation. The parallel formation of crab‐adapted ecotypes that exhibit partial reproductive isolation from wave‐adapted ecotypes has enabled genomic investigation of conspicuous shell traits. Recent genomic studies suggest that chromosomal rearrangements may enable ecotype divergence by reducing gene flow. However, the genomic architecture of traits that are divergent between ecotypes remains poorly understood. Here, we use 11,504 single nucleotide polymorphism (SNP) markers called using the recently released L. saxatilis genome to genotype 462 crab ecotype, wave ecotype and phenotypically intermediate Spanish L. saxatilis individuals with scored phenotypes. We used redundancy analysis to study the genetic architecture of loci associated with shell shape, shape corrected for size, shell size and shell ornamentation, and to compare levels of co‐association among different traits. We discovered 341 SNPs associated with shell traits. Loci associated with trait divergence between ecotypes were often located inside putative chromosomal rearrangements recently characterized in Swedish L. saxatilis. In contrast, we found that shell shape corrected for size varied primarily by geographic site rather than by ecotype and showed little association with these putative rearrangements. We conclude that genomic regions of elevated divergence inside putative rearrangements were associated with divergence of L. saxatilis ecotypes along steep environmental axes—consistent with models of adaptation with gene flow—but were not associated with divergence among the three geographical sites. Our findings support predictions from models indicating the importance of genomic regions of reduced recombination allowing co‐association of loci during ecological speciation with ongoing gene flow.     

A Widespread Chromosomal Inversion Polymorphism Contributes to a Major Life-History Transition,Local Adaptation,and Reproductive Isolation

The role of chromosomal inversions in adaptation and speciation is controversial. Historically, inversions were thought to contribute to these processes either by directly causing hybrid sterility or by facilitating the maintenance of co-adapted gene complexes. Because inversions suppress recombination when heterozygous, a recently proposed local adaptation mechanism predicts that they will spread if they capture alleles at multiple loci involved in divergent adaptation to contrasting environments. Many empirical studies have found inversion polymorphisms linked to putatively adaptive phenotypes or distributed along environmental clines. However, direct involvement of an inversion in local adaptation and consequent ecological reproductive isolation has not to our knowledge been demonstrated in nature. In this study, we discovered that a chromosomal inversion polymorphism is geographically widespread, and we test the extent to which it contributes to adaptation and reproductive isolation under natural field conditions. Replicated crosses between the prezygotically reproductively isolated annual and perennial ecotypes of the yellow monkeyflower, Mimulus guttatus, revealed that alternative chromosomal inversion arrangements are associated with life-history divergence over thousands of kilometers across North America. The inversion polymorphism affected adaptive flowering time divergence and other morphological traits in all replicated crosses between four pairs of annual and perennial populations. To determine if the inversion contributes to adaptation and reproductive isolation in natural populations, we conducted a novel reciprocal transplant experiment involving outbred lines, where alternative arrangements of the inversion were reciprocally introgressed into the genetic backgrounds of each ecotype. Our results demonstrate for the first time in nature the contribution of an inversion to adaptation, an annual/perennial life-history shift, and multiple reproductive isolating barriers. These results are consistent with the local adaptation mechanism being responsible for the distribution of the two inversion arrangements across the geographic range of M. guttatus and that locally adaptive inversion effects contribute directly to reproductive isolation. Such a mechanism may be partially responsible for the observation that closely related species often differ by multiple chromosomal rearrangements.     

A microsatellite linkage map for Drosophila montana shows large variation in recombination rates,and a courtship song trait maps to an area of low recombination

Current advances in genetic analysis are opening up our knowledge of the genetics of species differences, but challenges remain, particularly for out‐bred natural populations. We constructed a microsatellite‐based linkage map for two out‐bred lines of Drosophila montana derived from divergent populations by taking advantage of the Drosophila virilis genome and available cytological maps of both species. Although the placement of markers was quite consistent with cytological predictions, the map indicated large heterogeneity in recombination rates along chromosomes. We also performed a quantitative trait locus (QTL) analysis on a courtship song character (carrier frequency), which differs between populations and is subject to strong sexual selection. Linkage mapping yielded two significant QTLs, which explained 3% and 14% of the variation in carrier frequency, respectively. Interestingly, as in other recent studies of traits which can influence speciation, the strongest QTL mapped to a genomic region partly covered by an inversion polymorphism.     

Chromosome mapping of Miller-Diecker,Smith-Magenis and RARA loci in non-human primates: implications in the evolution of human chromosome 17

Molecular cytogenetics allows to verify chromosomal homologies previously hypothesised on the base of banding pattern comparison in different species. So far only the chromosome painting technique has been extensively used in studies of chromosomal evolution. This technique allows to detect only interchromosomal rearrangements. Human and Great Apes chromosomes basically differ by intrachromosomal rearrangements, in particular inversions; with chromosome painting it has just been possible to confirm the origin by fusion of human chromosome 2 and a reciprocal translocation in Gorilla, involving the homologous of chromosome 5 and 17. In order to verify intrachromosomal rearrangements in human chromosomal evolution, chromosome mapping of human loci in non-human primates is a useful approach. We mapped Miller-Diecker, Smith-Magenis and RARA loci localised on human chromosome 17, in Gorilla gorilla, Pongo pygmaeus, Macaca fascicularis and Cercopithecus aethiops. On the base of the obtained results it was possible to verify chromosomal rearrangements previously identified by banding, to achieve new informations about the controversial evolution of human chromosome 17, and to detect the occurrence of a paracentric inversion in the homologous in Cercopithecus aethiops.     

Mapping QTLs for defective female gametophyte development in an inter-subspecific cross in Oryza sativa L.

The embryo-sac is an essential structure for angiosperm reproduction. The cytological and genetic characterization of embryo-sac sterility was examined in a cross between Oryza sativa ssp. indica cv. ZYQ8 and ssp. japonica cultivar, JX17. The arrest of embryo-sac development was manifested following meiosis in the F₁ hybrid. When the megaspore carried the lethal genotype, the nucleus either failed to divide or divided only once, and the immature embryo-sac degenerated. Abortion of the embryo-sac in the indica-japonica hybrid background was not observed in their original parents, and an effect of cytoplasmic gene(s) on embryo-sac sterility in the reciprocal F₁ hybrids was not detected. Using a rice molecular linkage map based on a doubled haploid (DH) population from the cross of ZYQ8 /JX17, we mapped quantitative trait loci (QTLs) for the defective development of the female gametophyte in backcross progenies from the DH lines. The result demonstrated that a polygenic system is involved in both megagametogenesis and postzygotic isolation in inter-subspecific hybrid rice. Received: 4 May 2000 / Accepted: 20 September 2000     

INTROGRESSION IN THE DROSOPHILA SUBOBSCURA—D. MADEIRENSIS SISTER SPECIES: EVIDENCE OF GENE FLOW IN NUCLEAR GENES DESPITE MITOCHONDRIAL DIFFERENTIATION

Species hybridization, and thus the potential for gene flow, was once viewed as reproductive mistake. However, recent analysis based on large datasets and newly developed models suggest that gene exchange is not as rare as originally suspected. To investigate the history and speciation of the closely related species Drosophila subobscura, D. madeirensis, and D. guanche, we obtained polymorphism and divergence data for 26 regions throughout the genome, including the Y chromosome and mitochondrial DNA. We found that the D. subobscura X/autosome ratio of silent nucleotide diversity is significantly smaller than the 0.75 expected under neutrality. This pattern, if held genomewide, may reflect a faster accumulation of beneficial mutations on the X chromosome than on autosomes. We also detected evidence of gene flow in autosomal regions, while sex chromosomes remain distinct. This is consistent with the large X effect on hybrid male sterility seen in this system and the presence of two X chromosome inversions fixed between species. Overall, our data conform to chromosomal speciation models in which rearrangements are proposed to serve as gene flow barriers. Contrary to other observations in Drosophila, the mitochondrial genome appears resilient to gene flow in the presence of nuclear exchange.