<title>Introduction</title> <p>During the modern synthesis, Dobzhansky and Mayr proposed the Biological Species Concept that defines species by their ability to exchange genetic material within their group while being prevented from exchanging genetic material between groups ##UREF##0##[1]##, ##UREF##1##[2]##. Reproductive isolating mechanisms subsequently have been the focus of extensive research such that a great deal now is understood about mechanisms of both pre and postzygotic isolation. Capturing the process of speciation early enough to determine the initial genetic causes of reproductive isolation, however, is challenging. Consequently a gap still remains in our understanding of the mechanisms underlying the very first stages of postzygotic isolation.</p> <p>Postzygotic isolation exists when a hybrid, or its progeny, experience a reduction in fitness relative to the parental types. Whether the genetic variation needed for postzygotic isolation is segregating within species in the form of epistatic variation ##REF##10937210##[3]##, ##REF##28565332##[4]## or instead arises as de novo mutations in allopatric populations ##REF##10937209##[5]## remains controversial. Thus, there is need for further empirical study to determine the roles of segregating and de novo polymorphism in speciation, and the underlying architecture of such variation ##REF##10937210##[3]##.</p> <p>Three general patterns of postzygotic isolation are well supported. First, Haldane's rule ##UREF##2##[6]## states that decreases in hybrid fitness are more common in the heterogametic sex; second, hybrid sterility often appears earlier in species divergence than hybrid inviability ##REF##28569002##[7]##, ##REF##28568795##[8]##; and third, heterogametic isolation generally appears earlier in speciation than homogametic isolation ##REF##28568795##[8]##, ##REF##28568554##[9]##. Taken together, these patterns suggest that the earliest manifestation of postzygotic isolation likely will be hybrid sterility in the heterogametic sex. These patterns hold across diverse taxa in which the male is the heterogametic sex, such as mammals and <italic>Drosophila</italic> and interestingly, also for systems where the female is the heterogametic sex, such as birds and <italic>Lepidoptera</italic>
##UREF##3##[10]##. Various theories exist about the basis for Haldane's rule ##REF##7635302##[11]##–##UREF##4##[13]##. The explanations with the most consistent empirical support revolve around the X (Z)-chromosome ##REF##17850182##[14]##, leading to the expectation that the X-chromosome plays a substantial genetic role in the occurrence of postzygotic isolation in any given system.</p> <p>The problem of how a trait that decreases the fitness of hybrids could evolve between two populations sharing a common ancestral genome without either population passing through a state of reduced fitness is solved by the Dobzhansky-Muller (DM) model of reproductive isolation ##UREF##0##[1]##, ##UREF##5##[15]##–##UREF##7##[17]##. For example, if derived alleles of two loci arise independently in two populations then, upon meeting in a hybrid, they may be incompatible, leading to decreased fitness ##REF##7789779##[18]##. Several known interacting genomic regions have been implicated in hybrid incompatibilities ##REF##8765678##[19]##–##REF##9624532##[23]##. In addition, the incompatible alleles of the <italic>Lhr</italic> and <italic>Hmr</italic> genes in the <italic>Drosophila melanogaster – D. simulans</italic> pair have been demonstrated functionally to cause hybrid male lethality ##REF##17124320##[24]##.</p> <p>Expanding the Dobzhansky-Muller model, Orr ##REF##7789779##[18]## proposed that the number of two locus incompatibilities between two populations will increase as a function of the square of their divergence times. If more complex genetic scenarios are considered, the accumulation of incompatibilities is even faster. Those genetic incompatibilities that initially defined the two species therefore, can rapidly become obscured by subsequent genetic differentiation. Efforts to elucidate the origin of isolation can be facilitated by identifying the architecture of postzygotic isolating mechanisms in incipient species, where the number of postzygotic incompatibilities are still few and, ideally, not yet fixed within the populations</p> <p>A major lingering question in speciation genetics concerns the actual genetic architecture of postzygotic isolation as it first emerges. Does it have a simple or a complex genetic basis at the earliest stages of the speciation process? A burdensome disadvantage to genetic studies of postzygotic isolation is that the phenotype of interest, sterility or inviability, precludes direct genetic crossing schemes. Many successful studies of species pairs, in which the factors underlying isolation are effectively fixed, utilized the common occurrence of asymmetry ##REF##17435235##[25]## and/or Haldane's Rule. While the genetic dissections of incompatibility phenotypes by backcrosses and introgression ##UREF##3##[10]##, ##REF##16415357##[21]##, ##REF##17124320##[24]##, ##REF##12586709##[26]##, ##REF##12663535##[27]## have contributed significantly to our current understanding of the genetics of postzygotic isolation, it is not without limitations. Since the crossing schemes employed potentially disrupted the very co-adapted gene complexes that might be critical to the F<sub>1</sub> phenotype we cannot be certain that the genetic architecture for the incompatibility observed in the F<sub>1</sub> hybrid has been fully characterized. Evidence of substantial within species epistasis contributing to postzygotic isolation has been demonstrated in <italic>Tribolium</italic> by joint-scaling analysis ##REF##17348915##[28]##, ##REF##17348932##[29]## showing that co-adapted gene complexes may indeed be important in the architecture of hybrid incompatibilities.</p> <p>In the present study, we circumvented the limitations of the hybrid phenotype by exploiting within-species polymorphism for between-species postzygotic isolation to determine the genetic architecture of segregating variation for F<sub>1</sub> hybrid male sterility directly. We used <italic>Drosophila mojavensis</italic> and <italic>D. arizonae</italic>, a recently diverged species pair (0.66–1.2 my ##REF##14660699##[30]##, ##REF##17305857##[31]##) and a model system for the study of speciation ##UREF##10##[32]##. Crosses between <italic>D. mojavensis</italic> females and <italic>D. arizonae</italic> males exhibit widespread between-population ##REF##15184657##[33]##, ##REF##2332612##[34]## and within-population variation for the HMS phenotype suggesting polymorphism at multiple loci within <italic>D. mojavensis</italic>
##REF##15184657##[33]##. The polymorphism for HMS shows that this trait is not yet fixed in <italic>D. mojavensis</italic>. Since intraspecific evolution is the process by which gene frequencies change, only a trait exhibiting genetic variation is capable of evolving. Thus HMS in this system is still capable of experiencing the forces of evolution (selection and/or drift) and can give insight into how HMS evolves. In a novel design, we have utilized the within-population variation for production of F<sub>1</sub> HMS to determine the genomic regions and epistatic interactions that shape the architecture of postzygotic isolation within the F<sub>1</sub> hybrid males. Below we address the following questions: What is the nature and number of polymorphic loci within <italic>D. mojavensis</italic> contributing to the HMS phenotype? And what can that tell us, in conjunction with other information about this species pair, about the evolution of HMS?</p>
<title>Materials and Methods</title> <title>Lines used for mapping</title> <p>Genetic lines were derived from lines used in the HMS polymorphism study by Reed and Markow ##REF##15184657##[33]##. Isofemale <italic>D. mojavensis</italic> lines were collected from Santa Catalina Island (CI) in April 2001. The Low (CI-10) line exhibiting low levels of HMS was inbred through full-sib mating for six generations and the High (CI-12) line exhibiting high levels of HMS was inbred for eight generations. The lines had minimum inbreeding coefficients of 0.734 and 0.826 respectively ##UREF##11##[35]##. They showed significant differences from each other for HMS (presence of motile sperm) observed when females from those lines were crossed to <italic>D. arizonae</italic> males from an inbred line (minimum inbreeding coefficient of 0.859) collected at Peralta Canyon Trail Head, east of Phoenix, Arizona (AZ), in April 1997. All lines show normal male fertility when the males are genetically pure <italic>D. mojavensis</italic> or <italic>D. arizonae</italic>
##REF##15184657##[33]##. The sterility phenotype is only manifested in the species hybrid state. Homozygosity of the High line, used in the <italic>D. mojavensis</italic> genome sequencing project, had been confirmed by sequencing five marker loci in 10 females (L. Matzkin, <italic>unpublished</italic>) and an additional seven loci in 10 individuals by Agencourt Biosciences, and none showed sequence variation. The two <italic>D. mojavensis</italic> QTL lines were homokaryotypic, lacking inversion polymorphism. The lines used showed extreme contrasting HMS (sperm motility) phenotypes in heterospecific crosses (##TAB##0##Table 1##) as well as strong viability.</p> <title>Identifying Markers and Genotyping</title> <p>We designed fluorescently-tagged PCR primers for flanking sequences of microsatellite loci identified by Ross <italic>et al</italic>. ##REF##12777536##[36]##. Additional candidate microsatellite sequences were drawn from Staten <italic>et al.</italic>
##UREF##12##[37]##. Potential loci were genotyped in a panel of four individuals from each <italic>D. mojavensis</italic> line to confirm that the lines contained different alleles for the locus. We genotyped 132 potential loci but only found 25 that distinguished the High and Low lines (primer sequences not already described in Staten <italic>et al</italic>. ##UREF##12##[37]## shown in ##SUPPL##0##Table S1##). Polymorphic marker discovery was inhibited by the general deficit of neutral microsatellite variation in the wild population from which the mapping lines were derived ##REF##16910998##[38]##. We had markers on the X-chromosome and the four major autosomes but failed to find a marker on the 6<sup>th</sup> dot chromosome. Markers were assigned to chromosome using one or more of the following methods. 1) BLASTing ##REF##8162065##[39]## the flanking sequence against the <italic>D. melanogaster</italic> genome sequence to identify which Muller element and thereby to what chromosome in <italic>D. mojavensis</italic> it belonged; 2) linkage mapping to allozyme markers assigned to chromosome by Zouros ##UREF##13##[40]##, ##UREF##14##[41]##; 3) looking for close linkage between assigned markers and unassigned makers. Informative loci were multiplexed and PCR amplified using HotStarTaq DNA Polymerase (Qiagen Inc.). Genotyping was performed by the Genomic Analysis and Technology Core (GATC) at the University of Arizona using an ABI 3730 DNA Analyzer. Genotypes were read using Genotyper ver. 1.1 from Applied Biosystems.</p> <title>Crossing Design</title> <p>The crossing design used was an innovative modification of a standard F<sub>2</sub> line cross (##FIG##0##Figure 1##). Five virgin F<sub>2</sub> females, at six days post-eclosion, were paired with five sexually mature virgin male <italic>D. arizonae</italic> (AZ) for 48 hours. After 48 hours, the females were placed individually into vials, allowed to oviposit for four days, then removed and frozen for later genotyping. Virgin hybrid males were collected from those vials producing viable adults. Hybrid males were allowed to reach sexual maturity at 9–12 days post-eclosion before phenotyping. Use of the F<sub>2</sub> mothers' genotypes allowed for a largest possible sample of recombinant genotypes (given limited resources for genotyping) while the use of the sons' phenotypes allowed for a good estimate of their mother's average phenotype, maximizing the quality and quantity of the data. A reciprocal design was also performed; instead of using Low females crossed to High males at the first step (F<sub>2</sub> cross) High females were crossed to Low males (RF<sub>2</sub> cross). A large number of attempted matings between F<sub>2</sub> and RF<sub>2</sub> females and <italic>D. arizonae</italic> males failed to produce offspring as was expected considering the significant prezygotic isolation between the species. Seven percent of the F<sub>2</sub> (224 of 3214) and the 12% of the RF<sub>2</sub> (278 of 2339) mothers produced larvae after being mated to <italic>D. arizonae</italic> males.</p> <p>Mating and/or oviposition took place on standard opuntia-banana food (<ext-link ext-link-type="uri" xlink:href="http://flyfood.arl.arizona.edu/opuntia.php3#s2.1">http://flyfood.arl.arizona.edu/opuntia.php3#s2.1</ext-link>) while virgins were maintained on standard cornmeal food (<ext-link ext-link-type="uri" xlink:href="http://flyfood.arl.arizona.edu/cornmeal.php3">http://flyfood.arl.arizona.edu/cornmeal.php3</ext-link>). Cultures were maintained at 24+/−1 degree Celsius and under a strict 12h: 12h light: dark cycle.</p> <title>Phenotyping</title> <p>Male <italic>Drosophila</italic> store mature sperm in a seminal vesicle, a structure at the base of each testis. To assess sperm motility, males were anesthetized with ether just prior to dissection. Testes and seminal vesicles were dissected from the male in sperm buffer (0.05 M Tris, 1.1% NaCL, 0.1% Glucose, 0.01% L-arginine, 0.01% L-lysine, pH 8.7). Reproductive organs were transferred to a fresh slide with an 11 microliter droplet of sperm buffer, stretched with forceps to be linear, and then placed under a coverslip. Pressing the coverslip down on the sample forced all sperm out of the seminal vesicle, allowing motility of the sperm to be observed under dark-field microscopy. The entire area where potentially motile sperm could be found was visually scanned and each seminal vesicle assigned an integer score from 0 to 6 on a sperm motility scale.</p> <p>Unlike mammalian systems that have short sperm, <italic>Drosophila</italic> sperm are extremely long, forming a tangled wiggling knot, in which the movement of any individual sperm cannot be tracked. The zero-to-six scale is a reasonable means to approximate the amount of motility exhibited by <italic>Drosophila</italic> sperm, and is a refinement of the presence/absence of motility scheme used in many previous studies of sperm motility in <italic>Drosophila</italic>
##REF##15184657##[33]##, ##REF##11454758##[42]##, ##REF##16175455##[43]##. A quantitative measurement of sperm motility more suitable for quantitative genetic analysis results from the zero-to-six scale. A score of “0” meant there was no motility observed, “1” was one or two motile sperm, “2” being several motile sperm, “3” being several motile sperm in two or more areas of the field, “4” is three to four areas of moderate motility, “5” being several large areas of high motility, and “6” being wild-type where most of the field showed high levels of sperm motility (as is observed in non-hybrid males of these species). A male's phenotype was the average score of his two seminal vesicles. We phenotyped up to ten sons for each mother (mean family size of 3.4) and found no trends of decreasing motility scores with time (first verses second testis scored). The genotype of the male being scored was not known during phenotyping. A subset of samples were also scored and rescored in a randomized order while blind to sample identity. Replicate scores were highly correlated (R<sup>2</sup> = 0.79, p&lt;0.0001), validating this method for estimating <italic>Drosophila</italic> sperm motility suggesting its utility for future studies. Sperm motility is one phenotypic component of sterility. We have shown in previous work that motile sperm are required though not always sufficient for males to be able to reproduce in this system due to additional genetic polymorphism for other phenotypic components of sterility ##REF##15184657##[33]##. Studies of other phenotypic components of sterility would be interesting though not required to assess the importance of natural genetic variation for HMS as estimated by sperm motility in this system.</p> <title>Linkage Mapping</title> <p>F<sub>2</sub> or RF<sub>2</sub> females that successfully produced hybrid male offspring were genotyped for the 25 microsatellite markers described above. Maximum likelihood marker (##SUPPL##1##Figure S1##) order was calculated for each chromosome and the Haldane mapping function was used to assign linkage distances between markers in centiMorgans (cM) using Mapmaker 3.0 ##UREF##15##[44]##.</p> <p>Autosomal genotypes were pooled from the F<sub>2</sub> and RF<sub>2</sub> crosses for linkage map calculation while they are separated for the X-chromosome due to differences in X-chromosome inheritance patterns in the two crosses. Marker density was an average of 27.6 cM.</p> <title>QTL Mapping</title> <p>QTL mapping was conducted using Windows QTL Cartographer 2.5 ##UREF##16##[45]##, ##UREF##17##[46]##. We conducted composite interval mapping ##REF##8248199##[47]##, ##REF##8013918##[48]## on the F<sub>2</sub>, RF<sub>2</sub>, and combined data sets using a common linkage maps for the autosomes and their respective linkage maps for the X-chromosome. Composite interval mapping (CIM) uses cofactors from other regions to control for genome wide influences when regression analysis at any given genomic location is being tested ##REF##8248199##[47]##, ##REF##8013918##[48]##. For our study a multiple regression analysis at each 2 cM segment of each chromosome was performed using five randomly selected markers from elsewhere in the genome as cofactors to control for genetic effects in other genomic regions (Model 6 in QTLCartographer). A 10 cM region around the test location was excluded for selection of the five control markers. A likelihood ratio (LR) test was performed for each position comparing the null hypothesis that no QTL is present at that position to the alternative that there is a QTL present. To control for any violations of normality, multiple testing, or the possible disruption of the asymptotic distribution of the chi-squared test statistic by idiosyncrasies of the study (due to factors such as multiple QTL per chromosome or distorted segregation ratios), QTL peak significance was assessed by 1000 permutations ##REF##7851788##[49]##.</p> <p>Two datasets were used in mapping. In the one called <italic>sons</italic>, sons of each mother were analyzed with their mother's genotype and their individual motility phenotype to account for the within mother variation in her sons' motility scores. Analyses were also conducted on the second dataset using the mother's Best Linear Unbiased Prediction (BLUP) ##UREF##18##[50]##, ##UREF##19##[51]##, which considers both the mean and variance in offspring scores, as a summary of her “breeding quality” for hybrid male sperm motility calculated in the SAS/STAT Mixed Procedure (SAS 9.1.3, SAS Institute Inc. Cary, NC 2000–2004). The BLUP data is well suited to determine the basis of additive (heritable) variation while the <italic>sons</italic> data set is better for detecting genetic loci experiencing non-additive effects.</p> <title>Epistatic Interactions</title> <p>We assessed epistatic interactions using methods modified from Moehring and Mackay ##REF##15280239##[52]## and Morgan and Mackay ##REF##16404413##[53]##. We tested for interaction effects between all pairs of markers using a repeated measures model of the form:where F is the effect of mother, C is the effect of the direction of the cross (F<sub>2</sub> vs. RF<sub>2</sub>), M<sub>1</sub> is the effect of genotype at marker 1 and M<sub>2</sub> is the effect of genotype at marker 2. All 300 possible two-way interactions were tested and significance of the interaction term (M<sub>1</sub>*M<sub>2</sub>) was assessed by Bonferroni correction and False Discovery Rate (FDR). All calculations were implemented using the SAS Mixed procedure (SAS 9.1.3, SAS Institute Inc., Cary, NC 2000–2004).</p> <title>Maternal Effects</title> <p>Maternal effects were identified as a significant cross by main effect interaction. Hybrid male phenotype (n = 304) was regressed against the interaction between marker genotypes as the four main effect loci and the direction of the cross (F<sub>2</sub> vs. RF<sub>2</sub>) in the multiple regression model of the form:Where C in the direction of the cross and M<sub>1</sub> is the genotype at the marker associated with one of the main effect QTL.</p>
<title>Results</title> <p>Inbred <italic>D. mojavensis</italic> lines showing High and Low phenotypic levels of HMS when crossed to <italic>D. arizonae</italic> were used in a standard F<sub>2</sub> (and reciprocal) intercross QTL mapping design (##FIG##0##Figure 1##). Sperm motility in the parental, F<sub>1</sub>, RF<sub>1</sub>, F<sub>2</sub>, and RF<sub>2</sub> lines are presented in ##TAB##0##Table 1##. High and Low lines showed significant differences in the proportion of inter-species hybrid sons with motile sperm (91% and 46% respectively), and motility clearly was dominant in all hybrids (##TAB##0##Table 1##). Sperm motility is necessary but not always sufficient for fertility ##REF##15184657##[33]##, thus motility is used as a correlate to absolute fertility. Sperm motility scores were normally distributed (Shapiro-Wilks W = 0.99, p = 0.38), with a mean of 2.35+/−0.08.</p> <p>Female F<sub>2</sub>s (149 F<sub>2</sub> and 155 RF<sub>2</sub>), of the species <italic>D. mojavensis</italic>, were crossed to male <italic>D. arizonae</italic> for QTL analysis. Sperm motility was scored, from 0 to 6, in the interspecific hybrid males. 304 F<sub>2</sub> mothers were genotyped at 25 microsatellite markers across the genome and used to generate the linkage map. Phenotypes were mapped relative to genotype in two ways. One was to use each son's motility score, hereafter referred to as <italic>sons</italic>, as a sample of the mother's phenotype and the other was to map a mother's Best Linear Unbiased Prediction (BLUP) for her sons' motility scores. <italic>Sons</italic> better captures non-additive effects, while BLUP highlights the additive effects in the architecture.</p> <title>Main effects</title> <p>There were two striking outcomes of the main effects analysis; first that the F<sub>2</sub> and RF<sub>2</sub> mapping populations did not match more closely, and second, no main effects were detected on the X-chromosome. Main effect QTL were detected, however, on the 2<sup>nd</sup>, 3<sup>rd</sup>, and 5<sup>th</sup> chromosomes (##TAB##1##Table 2##, ##FIG##1##Figure 2##). Note that estimated effect sizes (percent of phenotypic variation explained) of all detected QTL are likely to be overestimates due to the Beavis effect ##UREF##20##[54]##, ##REF##12618414##[55]##, thus, there are likely to be additional genetic factors contributing to observed differences between lines that were not detected in this study. Also the actual genotype of the hybrid males at a given marker can only fall into two possible genotypic classes, High or Low for the <italic>D. mojavensis</italic> half of their genome: they will be monomorphic with respect to the <italic>D. arizonae</italic> half of their genome. Thus, the variance due to non-additive effects identified in the QTL mapping corresponds to intra-genome epistasis rather then intra-locus dominance. Finally, remember that these two lines are only a sample from a larger population. There may be additional polymorphic factors on the X-chromosome or autosomes segregating in the population as a whole that were not captured by this biparental cross. One QTL with estimated effect size of 10.9–22.1% occurred on the 3<sup>rd</sup> chromosome at about 60 cM and was observed in both datasets (<italic>sons</italic> and BLUP) and across all population samples (F<sub>2</sub>, RF<sub>2</sub>, and combined). The RF<sub>2</sub> population sample did not reveal any other significant QTL.</p> <p>QTL were also identified on the 2<sup>nd</sup> chromosome. Analysis of the <italic>sons</italic> dataset in the F<sub>2</sub> and combined populations found two and one QTL on the 2<sup>nd</sup> chromosome respectively (##TAB##1##Table 2##, ##FIG##1##Figure 2##). The combined population <italic>sons</italic> QTL at 58 cM has a large non-additive effect of 2.23 motility units and estimated effect size is 55.3%. The two QTL found in the F<sub>2</sub> population <italic>sons</italic> dataset (at 32 and 66 cM) also had large (and opposite) non-additive effects of −2.17 and 2.16 motility units, respectively. The RF<sub>2</sub>
<italic>sons</italic> dataset had an elevated likelihood ratio on the 2<sup>nd</sup> chromosome but it did not reach significance. The BLUP datasets exposed a significant QTL on 2nd chromosome for the combined population (60 cM) and the RF<sub>2</sub> population (50 cM) with a non-additive effect of 0.92 to 1.10 motility units, with estimated effect size of 44.3 to 59.3%. The lack of a significant 2<sup>nd</sup> chromosome QTL in the F<sub>2</sub> BLUP dataset can be explained by the apparent equal and opposite non-additive effects of the two QTL found in the <italic>sons</italic> dataset.</p> <p>Chromosome five also contained a QTL, but only in the F<sub>2</sub> population <italic>sons</italic> dataset. It was at 54.5 cM with a relatively large additive effect of −0.80 and moderate non-additive effect of −0.43 motility units with an estimated effect size of only 9%.</p> <p>We tested for biases due to variation in family size by rerunning all analyses of both the BLUP and <italic>sons</italic> datasets on families limited to three randomly selected sons. All QTLs replicated with controlled family size (data not shown) so we are confident that the mapping results are not due to biases in family size. We also looked for any biases introduced by segregation bias due to prezygotic isolation between the species or hybrid male over dominance in viability. We found that 15 of the 25 markers were not in the expected 1:2:1 (autosomes) or 1:1 (x-chromosome) ratios having excess heterozygotes, but there was no association between violation of these ratios and the location of the QTL (X<sup>2</sup> = 0.45, p = 0.50). Thus, there is evidence of genetic variation for prezygotic isolation with a possible heterozygote advantage that may be worth mapping in future studies, but that variation did not introduce bias into the mapping for HMS. The linkage map, though, is likely only accurate for the population that was genotyped and cannot be applied to the F<sub>2</sub>-females that failed to mother hybrid sons, because two or more segregation distortion loci within a given marker interval can bias estimates of recombination rate ##UREF##21##[56]##–##REF##24173787##[58]##.</p> <title>Architectural Complexity</title> <p>The lack of a segregating main effect on the X-chromosome was surprising, as this chromosome has been generally shown to play a major role in the evolution of hybrid male sterility ##REF##17850182##[14]##. There may still be factors on the X-chromosome (or autosomes) segregating in the ancestral population as a whole whose effects were not captured by the two lines used in this study, and there also may be fixed factors on the X-chromosome that we cannot detect with this design. Additionally, there were several major discrepancies between the F<sub>2</sub> and RF<sub>2</sub> populations. Specifically, four main effect QTL were identified in the F<sub>2</sub> population; two had largely additive effects on the 3<sup>rd</sup> and 5<sup>th</sup> chromosomes and two had largely non-additive effects on the 2<sup>nd</sup> chromosome. In the RF<sub>2</sub> population, however, only a single, largely additive QTL was found on the 3<sup>rd</sup> chromosome. The combined population data set shows one 2<sup>nd</sup> chromosome QTL and the 3<sup>rd</sup> chromosome QTL. These unexpected results motivated further analyses of the underlying architectural complexity of HMS.</p> <p>There are several possible explanations for the discrepancy between F<sub>2</sub> and RF<sub>2</sub> populations. First, we may have failed to detect QTL due to simple sampling variance between the F<sub>2</sub> and RF<sub>2</sub> populations. Second, the different compositions of the X-chromosomes in the reciprocal populations could have lead to an X-by-autosome interaction. Finally, an interaction between the nature of the cross (<italic>e.g.</italic> cytoplasm) and the QTLs could account for (or underlie) the observed differences between reciprocal crosses. We address each of these below.</p> <p>
<italic>Sampling Variance</italic>: To test whether the inevitable force of sampling variance was driving the differences observed between the F<sub>2</sub> and RF<sub>2</sub> populations, we subsampled from the datasets to test for consistency of the signal. From each population 100 datasets were created by randomly selecting 80% of the data points from the original dataset. Complete mapping analyses, including permutations for significance thresholds, were then performed on the subsampled datasets. Despite the reduced power of the smaller datasets, the QTLs found on the full dataset remained detectable in between 27% and 70% of the subsampled datasets (##SUPPL##2##Figure S2##), while spurious QTL were only detected in a maximum of 24% (<italic>e.g</italic>. see left portion of chromosome 5 in the F<sub>2</sub> dataset). Subsampling indicates that QTL identified in the complete dataset remain detectable in a large portion of the subsets despite variation in sampling and spurious QTL (those not detected in the complete dataset but in the subsets) rarely arise due to sampling variance. The consistency between our subsampled datasets and our complete dataset suggests that our findings are fairly robust to variation in sampling and that the differences between the F<sub>2</sub> and RF<sub>2</sub> are likely biologically real, though some of the differences could still be due in part to sampling variance.</p> <p>
<italic>Inheritance of the X-chromosome</italic>: If an important X-linked genetic difference between the parental lines interacts with the autosomes, the cross direction will affect the frequency of that factor in the mapping population. This, in turn, would influence the detection and nature of autosomal QTL that interact with the X. The X-linked factor could primarily influence the penetrance or expression of the genetic variation at the autosomal factor and thus have a significant epistatic effect without having a main effect. To test for X-by-autosome interactions, as well as for other epistasis, we looked for significant pairwise interactions between all markers across the combined dataset (n = 304). We found seven interactions with a FDR significance of 0.05 or less, two of which were between the X-chromosome and the autosomes (##TAB##2##Table 3##, ##FIG##1##Figure 2##). Epistatic interactions between the X-chromosome and the autosomes therefore could be contributing to some of the differences observed between the F<sub>2</sub> and RF<sub>2</sub> populations. We also found 5 inter-autosomal epistatic interactions. Notice that markers on the 4<sup>th</sup> chromosome were involved in five of the seven interactions, thus the 4<sup>th</sup> chromosome appears to be playing an important epistatic role in the HMS in this system worthy of additional study. Only one of nine markers involved in these pairwise interactions (m4_9_1 on the 5<sup>th</sup> chromosome) resides in or near a main effect QTL. Two additional markers (a1_2_1 on chromosome 3 and m1_10_11 on chromosome 4) did not appear in the list of significant marker by marker interactions, but exhibited significant deviations from the additive expectation after a Bonferroni correction, indicating that there are at least two other unidentified epistatic interactions. Taken together, not only does the X-chromosome appear to play an epistatic role in shaping the genetic architecture of HMS in this system, but so does within-species epistasis throughout the rest of the genome.</p> <p>
<italic>Maternal effects:</italic> Cytoplasmic factors may also underlie interactions between the direction of the cross and the QTLs. We tested for marker by cross interactions for the four markers showing main effects and found a significant interaction effect of marker dmoj2210 on the 2nd chromosome with cross (p = 0.0142), after controlling for the main effect of the 3<sup>rd</sup> chromosome (##FIG##2##Figure 3##). This marker did not show evidence of interacting with the X-chromosome in the test for epistasis above. Chromosome 2 is where the large and opposing non-additive QTL were found in the F<sub>2</sub> population but not in the RF<sub>2</sub> populations. The significant cross by marker interaction on the 2<sup>nd</sup> chromosome is a likely contributor to the differences seen in the mapping for the two crosses. The cross effect is not necessarily independent of the epistatic effect of the X –chromosome found above. But since dmoj2210 did not show any evidence of interacting with the X-chromosome, it is likely that the cytoplasm plays a role in HMS. Thus, both the X-chromosome and the cytoplasmic environment should be further explored in future studies of HMS.</p>
<title>Discussion</title> <p>The first question we proposed to answer in this study was: what is the basic architecture of HMS in the early stages of its evolution? Is HMS caused by only a few factors of large effect or does it have a more complex basis, involving multiple loci and epistatic interactions? An earlier study ##REF##15184657##[33]## showed evidence of multiple polymorphic factors for HMS within <italic>D. mojavensis</italic> and additional data show genetic polymorphisms for HMS segregating in the sister species <italic>D. arizonae</italic> (Reed, <italic>unpublished</italic>). Our present QTL study reveals that HMS, even at the earliest stages of speciation, can be surprisingly complex. We found within-population variation not only for several QTL of main effect, but also for significant within-species epistatic interactions. Similar complexity has been seen in other systems, such as with the <italic>Lhr-Hmr</italic> gene pair, where the incompatibility requires an epistatic interaction with the hybrid genetic background to be expressed ##UREF##22##[59]##. There is likely to be more segregating genetic variation in this species that was not captured by these two particular mapping lines, and fixed genetic factors that could not be detected by this design; potentially both of which further contribute to the complexity of this trait. While previous investigations on more diverged species pairs ##UREF##3##[10]## have revealed a multi-factored genetic basis for post-zygotic isolation, these studies were unable to determine whether the polygenicity had arisen early or was a byproduct of long divergence time. We show that HMS has a polygenic and epistatic basis even early in the evolutionary process, before the ability to produce HMS is fixed within a species. These findings are consistent with the complex genetic basis of early postzygotic isolation found in <italic>Tribolium</italic>
##REF##17348915##[28]##, ##REF##17348932##[29]##
</p> <p>In addition to the unexpected complexity of HMS in this incipient species pair, we were surprised to discover no main-effect QTL segregating on the X-chromosome, though we did find an epistatic effect. While we may have lacked the statistical power to detect X-linked main effects, main effects were detected on other chromosomes indicating that sample size alone is not the full explanation. There may still be fixed X-linked effects that could not be detected by this design, and fixed effects on the X are more likely due its smaller effective population size. Demuth and Wade did demonstrate the role of epistatic variation on the X-chromosome on Haldane's Rule both theoretically ##UREF##4##[13]## and empirically ##REF##17348915##[28]##, ##REF##17348932##[29]##. We also found such a role for the X-chromosome in our epistasis tests.</p> <p>Finally, our data show a striking difference between reciprocal crosses, suggesting a potential role for maternal effects and genome-wide epistatic variation on HMS. Cytoplasmic and epistatic effects have been demonstrated in other systems as well ##REF##17124320##[24]##, ##REF##17348915##[28]##, ##REF##17348932##[29]##, ##UREF##22##[59]##. Zouros and colleagues have mapped the genetic basis of HMS resulting from the reciprocal cross (female <italic>D. arizonae</italic> mated to male <italic>D. mojavensis</italic>) to the third, fourth and Y-chromosomes ##UREF##13##[40]##, ##UREF##14##[41]##, ##REF##28563497##[60]##–##REF##28581068##[62]##. Whether or not the third chromosome effect in the cross using <italic>D. arizonae</italic> mothers is a function of the same factor or factors as revealed in our study remains unknown.</p> <p>The second question we proposed to answer was: what can the genetic architecture of HMS tell us, in conjunction with other information, about its evolution? Substantial within-species natural polymorphism for between species postzygotic isolation has been characterized in many independent taxa thus far, including at least nine <italic>Drosophila</italic> species groups, ##REF##15184657##[33]##, ##REF##2332612##[34]##, ##REF##16175455##[43]##, ##REF##17246918##[63]##–##UREF##26##[70]##, <italic>Tribolium</italic> beetles ##REF##8150655##[71]##, ##REF##9383066##[72]##, <italic>Chorthippus parallelus</italic> grasshoppers ##UREF##27##[73]##, <italic>Mus musculus</italic>
##REF##18005156##[74]##, <italic>Xiphophorus</italic> fish ##REF##8917541##[75]##, <italic>Crepis</italic> hawksberg weed ##REF##17246594##[76]##, <italic>Gossypium</italic> cotton ##REF##24538920##[77]## and <italic>Mimulus</italic> monkey flowers ##REF##17300433##[78]##. Thus, intraspecific polymorphism for postzygotic isolation appears to be a common phenomenon. Shuker <italic>et al</italic>. ##UREF##27##[73]## argue that witnessing substantial within-species polymorphism for hybrid incompatibilities means that loci contributing to the hybrid phenotype are likely to be neutral or nearly-neutral in the parental species' genetic background. Finding that a trait is polymorphic at multiple loci, supports the role of drift and/or balancing selection in the genetics of the trait. If directional selection was occurring on the genes related to the trait, no matter what the within- species function of the genes, they would be polymorphic for only the briefest of times ##REF##17248440##[79]##. The probability of catching multiple polymorphic loci under directional selection would be minute. Considering that there is evidence of balancing selection on male reproductive proteins in several systems ##REF##17033959##[80]##–##REF##16614223##[83]## it is not unreasonable to hypothesize that it might also be playing a role in HMS.</p> <p>A handful of genes causing hybrid incompatibilities have been identified in other species and all but one (a gene transposition ##REF##16960009##[84]##) show evidence of being under positive selection (<italic>OdsH</italic> in <italic>D. mauritiana</italic>, ##REF##9822383##[85]##; <italic>Xmrk-2</italic> in <italic>Xiphophorus</italic>, ##REF##10545466##[86]##; <italic>Nup96</italic>, <italic>Hmr</italic>, <italic>Lhr</italic> in the <italic>D. melanogaster/D. simulans</italic> species pair ##REF##17124320##[24]##, ##REF##12695567##[87]##, ##REF##12802326##[88]##), leading to the impression that most postzygotic speciation genes experience directional selection. Given our findings with <italic>D. mojavensis</italic> and in other systems discussed above, however, we believe it is premature to discount the role of drift and balancing selection in speciation. Empirically, we need to characterize more completely the genetic architecture of incompatibility phenotypes during early stages to find reproductive isolation genes before they are fixed and test those loci for evidence of selection. In the present study we have begun the empirical aspect of this process.</p>
<p>Conceived and designed the experiments: LKR TAM. Performed the experiments: LKR BAL. Analyzed the data: LKR. Wrote the paper: LKR TAM.</p> <title>Background</title> <p>The genetic basis of postzygotic isolation is a central puzzle in evolutionary biology. Evolutionary forces causing hybrid sterility or inviability act on the responsible genes while they still are polymorphic, thus we have to study these traits as they arise, before isolation is complete.</p> <title>Methodology/Principal Findings</title> <p>Isofemale strains of <italic>D. mojavensis</italic> vary significantly in their production of sterile F<sub>1</sub> sons when females are crossed to <italic>D. arizonae</italic> males. We took advantage of the intraspecific polymorphism, in a novel design, to perform quantitative trait locus (QTL) mapping analyses directly on F<sub>1</sub> hybrid male sterility itself. We found that the genetic architecture of the polymorphism for hybrid male sterility (HMS) in the F<sub>1</sub> is complex, involving multiple QTL, epistasis, and cytoplasmic effects.</p> <title>Conclusions/Significance</title> <p>The role of extensive intraspecific polymorphism, multiple QTL, and epistatic interactions in HMS in this young species pair shows that HMS is arising as a complex trait in this system. Directional selection alone would be unlikely to maintain polymorphism at multiple loci, thus we hypothesize that directional selection is unlikely to be the only evolutionary force influencing postzygotic isolation.</p>
<title>Supporting Information</title>
<p>We greatly appreciated the assistance from M. Wasserman in confirming homosequentiality of the two <italic>D. mojavensis</italic> QTL lines by inspection of polytene chromosomes. We thank C. Ross for assistance in collecting specimens and M. Lowrance for assistance with inbreeding. We also appreciate advice on this study and manuscript from I. Dworkin, G. Gibson, E. Kelleher, C. Machado, T. Mackay, L. Matzkin, T. Morgan, D. Nielsen, D. Papaj, B. Payseur, and B. Walsh. We also are greatly indebted to the Santa Catalina Island Conservancy for collecting permits and logistical support in collecting specimens.</p>
<fig id="pone-0003076-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003076.g001</object-id><label>Figure 1</label><caption><title>Crossing scheme.</title><p>Design used to produce the recombinant <italic>D. mojavensis</italic> females to be mated to inbred <italic>D. arizonae</italic> males.</p></caption></fig> <fig id="pone-0003076-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003076.g002</object-id><label>Figure 2</label><caption><title>Composite interval mapping results for intraspecific variation for interspecific hybrid male sterility.</title><p>Likelihood ratio values for the F<sub>2</sub> (dashed red), RF<sub>2</sub> (dashed blue), and combined (solid black) data sets are given for the X through 5<sup>th</sup> chromosome. The X-chromosome scale is normalized to “combined” map length for the F<sub>2</sub> and RF<sub>2</sub> maps. The threshold for significance is the horizontal line (red- F<sub>2</sub>, blue- RF<sub>2</sub>, and black- combined). (a) Mapping results considering each son as measure of the mother's phenotype. Black asterisks indicate the location of markers with significant pairwise epistatic effects and yellow asterisks indicate two additional markers showing a significant deviation from the additive expectation. There is clear evidence for a QTL of major effect on the 3<sup>rd</sup> chromosome for all data sets, on the 2<sup>nd</sup> chromosome for the F<sub>2</sub> and combined data sets, and on the 5<sup>th</sup> for the F<sub>2</sub> dataset. (b) Mapping results of the BLUP values for each mother. There is evidence for a QTL on major effect on the 3<sup>rd</sup> chromosome in all three datasets and on the 2<sup>nd</sup> chromosome in the combined dataset.</p></caption></fig> <fig id="pone-0003076-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003076.g003</object-id><label>Figure 3</label><caption><title>Interaction effect between marker genotype and cross-type on motility score.</title><p>The average motility score of the homozygous marker genotype on the second chromosome (dmoj2210-High or Low parental line) relative to which cross it was derived (F<sub>2</sub> or RF<sub>2</sub>). There is a significant interaction between marker and cross type.</p></caption></fig>
<table-wrap id="pone-0003076-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003076.t001</object-id><label>Table 1</label><caption><title>Hybrid Sperm Motility in sons of mothers from the High and Low <italic>D. mojavensis</italic> lines and in the intercross populations of the High and Low lines, when females crossed to <italic>D. arizonae</italic> fathers.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1"/><td colspan="2" align="left" rowspan="1">Proportion (N)</td><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">Cross<xref ref-type="table-fn" rid="nt101">1</xref>
</td><td align="left" rowspan="1" colspan="1">No Motility</td><td align="left" rowspan="1" colspan="1">Motility<xref ref-type="table-fn" rid="nt102">2</xref>
</td><td align="left" rowspan="1" colspan="1">Standard Error</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Low</td><td align="left" rowspan="1" colspan="1">0.54 (35)</td><td align="left" rowspan="1" colspan="1">0.46 (30)</td><td align="left" rowspan="1" colspan="1">0.062</td></tr><tr><td align="left" rowspan="1" colspan="1">F<sub>1</sub>
</td><td align="left" rowspan="1" colspan="1">0.07 (11)</td><td align="left" rowspan="1" colspan="1">0.93 (154)</td><td align="left" rowspan="1" colspan="1">0.019</td></tr><tr><td align="left" rowspan="1" colspan="1">F<sub>2</sub>
</td><td align="left" rowspan="1" colspan="1">0.15 (62)</td><td align="left" rowspan="1" colspan="1">0.85 (353)</td><td align="left" rowspan="1" colspan="1">0.017</td></tr><tr><td align="left" rowspan="1" colspan="1">RF<sub>1</sub>
</td><td align="left" rowspan="1" colspan="1">0.07 (13)</td><td align="left" rowspan="1" colspan="1">0.93 (177)</td><td align="left" rowspan="1" colspan="1">0.018</td></tr><tr><td align="left" rowspan="1" colspan="1">RF<sub>2</sub>
</td><td align="left" rowspan="1" colspan="1">0.13 (75)</td><td align="left" rowspan="1" colspan="1">0.87 (483)</td><td align="left" rowspan="1" colspan="1">0.014</td></tr><tr><td align="left" rowspan="1" colspan="1">High</td><td align="left" rowspan="1" colspan="1">0.09 (13)</td><td align="left" rowspan="1" colspan="1">0.91 (131)</td><td align="left" rowspan="1" colspan="1">0.024</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0003076-t002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003076.t002</object-id><label>Table 2</label><caption><title>Maximum Likelihood Position and Effect Size for individual <italic>son</italic> and BLUP QTL found in Composite Interval Mapping.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Data Type<xref ref-type="table-fn" rid="nt103">1</xref>
</td><td align="left" rowspan="1" colspan="1">Cross Type</td><td align="left" rowspan="1" colspan="1">Chromosome</td><td align="left" rowspan="1" colspan="1">Position (cM<xref ref-type="table-fn" rid="nt104">2</xref>)</td><td align="left" rowspan="1" colspan="1">LR<xref ref-type="table-fn" rid="nt105">3</xref>
</td><td align="left" rowspan="1" colspan="1">Additive Effect</td><td align="left" rowspan="1" colspan="1">Non-Additive Effect</td><td align="left" rowspan="1" colspan="1">R<xref ref-type="table-fn" rid="nt104">2</xref>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">IS</td><td align="left" rowspan="1" colspan="1">F<sub>2</sub>
</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">32</td><td align="left" rowspan="1" colspan="1">23.68</td><td align="left" rowspan="1" colspan="1">0.24</td><td align="left" rowspan="1" colspan="1">−2.165</td><td align="left" rowspan="1" colspan="1">0.488</td></tr><tr><td align="left" rowspan="1" colspan="1">IS</td><td align="left" rowspan="1" colspan="1">F<sub>2</sub>
</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">66</td><td align="left" rowspan="1" colspan="1">27.12</td><td align="left" rowspan="1" colspan="1">−0.381</td><td align="left" rowspan="1" colspan="1">2.164</td><td align="left" rowspan="1" colspan="1">0.506</td></tr><tr><td align="left" rowspan="1" colspan="1">IS</td><td align="left" rowspan="1" colspan="1">F<sub>2</sub>
</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">61.2</td><td align="left" rowspan="1" colspan="1">41.64</td><td align="left" rowspan="1" colspan="1">−0.902</td><td align="left" rowspan="1" colspan="1">−0.137</td><td align="left" rowspan="1" colspan="1">0.157</td></tr><tr><td align="left" rowspan="1" colspan="1">IS</td><td align="left" rowspan="1" colspan="1">F<sub>2</sub>
</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">54.5</td><td align="left" rowspan="1" colspan="1">28.39</td><td align="left" rowspan="1" colspan="1">−0.802</td><td align="left" rowspan="1" colspan="1">−0.427</td><td align="left" rowspan="1" colspan="1">0.096</td></tr><tr><td align="left" rowspan="1" colspan="1">IS</td><td align="left" rowspan="1" colspan="1">RF<sub>2</sub>
</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">59.2</td><td align="left" rowspan="1" colspan="1">46.24</td><td align="left" rowspan="1" colspan="1">−0.79</td><td align="left" rowspan="1" colspan="1">0.357</td><td align="left" rowspan="1" colspan="1">0.127</td></tr><tr><td align="left" rowspan="1" colspan="1">IS</td><td align="left" rowspan="1" colspan="1">Combined</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">58</td><td align="left" rowspan="1" colspan="1">56.82</td><td align="left" rowspan="1" colspan="1">−0.392</td><td align="left" rowspan="1" colspan="1">2.225</td><td align="left" rowspan="1" colspan="1">0.553</td></tr><tr><td align="left" rowspan="1" colspan="1">IS</td><td align="left" rowspan="1" colspan="1">Combined</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">61.2</td><td align="left" rowspan="1" colspan="1">73.24</td><td align="left" rowspan="1" colspan="1">−0.739</td><td align="left" rowspan="1" colspan="1">0.132</td><td align="left" rowspan="1" colspan="1">0.109</td></tr><tr><td align="left" rowspan="1" colspan="1">BLUP</td><td align="left" rowspan="1" colspan="1">F<sub>2</sub>
</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">69.2</td><td align="left" rowspan="1" colspan="1">18.98</td><td align="left" rowspan="1" colspan="1">−0.344</td><td align="left" rowspan="1" colspan="1">−0.046</td><td align="left" rowspan="1" colspan="1">0.122</td></tr><tr><td align="left" rowspan="1" colspan="1">BLUP</td><td align="left" rowspan="1" colspan="1">RF<sub>2</sub>
</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">50</td><td align="left" rowspan="1" colspan="1">12.07</td><td align="left" rowspan="1" colspan="1">−0.255</td><td align="left" rowspan="1" colspan="1">1.101</td><td align="left" rowspan="1" colspan="1">0.593</td></tr><tr><td align="left" rowspan="1" colspan="1">BLUP</td><td align="left" rowspan="1" colspan="1">RF<sub>2</sub>
</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">59.2</td><td align="left" rowspan="1" colspan="1">23.64</td><td align="left" rowspan="1" colspan="1">−0.476</td><td align="left" rowspan="1" colspan="1">0.294</td><td align="left" rowspan="1" colspan="1">0.221</td></tr><tr><td align="left" rowspan="1" colspan="1">BLUP</td><td align="left" rowspan="1" colspan="1">Combined</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">60</td><td align="left" rowspan="1" colspan="1">13.8</td><td align="left" rowspan="1" colspan="1">−0.197</td><td align="left" rowspan="1" colspan="1">0.921</td><td align="left" rowspan="1" colspan="1">0.443</td></tr><tr><td align="left" rowspan="1" colspan="1">BLUP</td><td align="left" rowspan="1" colspan="1">Combined</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">65.2</td><td align="left" rowspan="1" colspan="1">37.95</td><td align="left" rowspan="1" colspan="1">−0.406</td><td align="left" rowspan="1" colspan="1">0.138</td><td align="left" rowspan="1" colspan="1">0.153</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0003076-t003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003076.t003</object-id><label>Table 3</label><caption><title>Significant Marker-Marker Interactions</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Marker A</td><td align="left" rowspan="1" colspan="1">Marker B</td><td align="left" rowspan="1" colspan="1">Chromosome A</td><td align="left" rowspan="1" colspan="1">Chromosome B</td><td align="left" rowspan="1" colspan="1">p-value</td><td align="left" rowspan="1" colspan="1">FDR</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">
<italic>dmojx030</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>a3_11_1</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>X</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>2</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>0.00112</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>0.0496</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>m2_18_2</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>m2_19_2</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>X</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>4</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>0.00072</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>0.0496</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>a3_11_1</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>m2_19_2</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>2</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>4</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>0.0006</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>0.0496</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>m2_17_15</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>m2_19_2</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>3</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>4</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>0.00116</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>0.0496</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>dmoj4050</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>dmoj4060</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>4</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>4</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>0.00116</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>0.0496</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold>m2_19_2</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>m4_9_1</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>4</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>5</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>0.00013</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>0.0393</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>a3_12_6</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>m4_9_1</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>5</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>5</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>0.0003</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>0.0453</italic>
</td></tr></tbody></table></alternatives></table-wrap>
<disp-formula></disp-formula> <disp-formula></disp-formula>
<supplementary-material content-type="local-data" id="pone.0003076.s001"><label>Table S1</label><caption><p>Primers for Microsatellites.</p><p>(0.06 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0003076.s002"><label>Figure S1</label><caption><p>Linkage map for D. mojavensis. Distance between markers given in centiMorgans.</p><p>(0.30 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0003076.s003"><label>Figure S2</label><caption><p>Subsampling effect on QTL detection. Bar graphs of count (right x-axis, out of 100) of 80%-subsampled datasets at each mapping interval showing a likelihood ratio greater than the significance threshold determined by permutation for each subset. Solid lines indicate likelihood ratio (left x-axis) calculated from the complete dataset. A. F2-sons dataset. B. RF2-sons dataset. Subsampling indicates that QTL identified in the complete dataset remain detectable in a large portion of the subsets and spurious QTL (those not detected in the completely dataset) rarely arise due to sampling variance. Thus, patterns observed in the complete dataset appear robust to sampling variance.</p><p>(0.16 MB DOC)</p></caption></supplementary-material>
<table-wrap-foot><fn id="nt101"><label>1</label><p>Cross is Genotype of the <italic>D. mojavensis</italic> mother.</p></fn><fn id="nt102"><label>2</label><p>Motility is defined as a male having one or more motile sperm.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt103"><label>1</label><p>Individual Son (IS) or the Best Linear Unbiased Prediction (BLUP) values.</p></fn><fn id="nt104"><label>2</label><p>cM is centiMorgans.</p></fn><fn id="nt105"><label>3</label><p>LR is likelihood ratio.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt106"><label/><p>Markers and their corresponding chromosomes with interaction effects with a False Discovery Rate (FDR) of 0.05 or less. Interactions significant after a Bonferroni correction are in bold and interactions significant at an FDR of 0.05 are italicized.</p></fn></table-wrap-foot> <fn-group><fn fn-type="COI-statement"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p><bold>Funding: </bold>National Science Foundation support for this work was provided by Interdisciplinary Research Training Group on Plant-Insect Interactions (DBI-9602249) to the University of Arizona and the Integrative Graduate Education and Research Training Grant in Evolutionary, Computational, and Molecular Approaches to Genome Structure and Function (DGE-0114420) to the University of Arizona, Doctoral Dissertation Improvement Grant (DEB 0414184) to LKR, and grants (DEB 95-10645 and DEB 0075312) to TAM. In addition, the Undergraduate Biological Research Program (UBRP) at the University of Arizona supported BAL. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</p></fn></fn-group>
<graphic id="pone-0003076-t001-1" xlink:href="pone.0003076.t001"/> <graphic xlink:href="pone.0003076.g001"/> <graphic xlink:href="pone.0003076.e001.jpg" mimetype="image" position="float"/> <graphic xlink:href="pone.0003076.e002.jpg" mimetype="image" position="float"/> <graphic xlink:href="pone.0003076.g002"/> <graphic id="pone-0003076-t002-2" xlink:href="pone.0003076.t002"/> <graphic id="pone-0003076-t003-3" xlink:href="pone.0003076.t003"/> <graphic xlink:href="pone.0003076.g003"/>
<media xlink:href="pone.0003076.s001.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0003076.s002.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0003076.s003.doc"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["1"], "element-citation": ["\n"], "surname": ["Dobzhansky"], "given-names": ["T"], "year": ["1937"], "article-title": ["Genetics and the origin of species."], "publisher-loc": ["New York"], "publisher-name": ["Columbia University"]}, {"label": ["2"], "element-citation": ["\n"], "surname": ["Mayr"], "given-names": ["E"], "year": ["1942"], "article-title": ["Systematics and the Origin of Species."], "publisher-loc": ["New York"], "publisher-name": ["Columbia University Press"]}, {"label": ["6"], "element-citation": ["\n"], "surname": ["Haldane"], "given-names": ["JBS"], "year": ["1922"], "article-title": ["Sex ratio and unisexual sterility in animal hybrids."], "source": ["Journal of Genetics"], "volume": ["12"], "fpage": ["101"], "lpage": ["109"]}, {"label": ["10"], "element-citation": ["\n"], "surname": ["Coyne", "Orr"], "given-names": ["J", "HA"], "year": ["2004"], "article-title": ["Speciation."], "publisher-loc": ["Sunderland, MA"], "publisher-name": ["Sinauer Associates Inc"]}, {"label": ["13"], "element-citation": ["\n"], "surname": ["Demuth", "Wade"], "given-names": ["JP", "MJ"], "year": ["2005"], "article-title": ["On the theoretical and empirical framework for studying genetic interactions within and among species."], "source": ["American Naturalist"], "volume": ["165"], "fpage": ["524"], "lpage": ["536"]}, {"label": ["15"], "element-citation": ["\n"], "surname": ["Muller"], "given-names": ["HJ"], "year": ["1942"], "article-title": ["Isolating mechanisms, evolution and temperature."], "source": ["Biological Symposia"], "volume": ["6"], "fpage": ["71"], "lpage": ["125"]}, {"label": ["16"], "element-citation": ["\n"], "surname": ["Muller", "Huxley"], "given-names": ["HJ", "J"], "year": ["1940"], "article-title": ["Bearing of the Drosophila work on systematics."], "source": ["The New Systematics"], "publisher-loc": ["Oxford"], "publisher-name": ["Clarendon"], "fpage": ["185"], "lpage": ["268"]}, {"label": ["17"], "element-citation": ["\n"], "surname": ["Bateson", "Seward"], "given-names": ["W", "AC"], "year": ["1909"], "article-title": ["Heredity and variation in modern lights."], "source": ["Darwin and Modern Science"], "publisher-loc": ["Cambridge, UK"], "publisher-name": ["Cambridge University Press"], "fpage": ["85"], "lpage": ["101"]}, {"label": ["20"], "element-citation": ["\n"], "surname": ["Christie", "Macnair"], "given-names": ["P", "MR"], "year": ["1984"], "article-title": ["Complementary Lethal Factors in 2 North-American Populations of the Yellow Monkey Flower."], "source": ["Journal of Heredity"], "volume": ["75"], "fpage": ["510"], "lpage": ["511"]}, {"label": ["22"], "element-citation": ["\n"], "surname": ["Hutchinson"], "given-names": ["JB"], "year": ["1932"], "article-title": ["The genetics of cotton. Part VII. \u201cCrumpled\u201d: A new dominant in Asiatic cottons produced by complementary factors."], "source": ["Journal of Genetics"], "volume": ["25"], "fpage": ["281"], "lpage": ["291"]}, {"label": ["32"], "element-citation": ["\n"], "surname": ["Markow", "Hocutt", "J.", "S.H."], "given-names": ["TA", "GD", "HD", "B"], "year": ["1998"], "article-title": ["Reproductive isolation in Sonoran Desert Drosophila: testing the limits of the rules."], "source": ["Endless Forms: species and speciation"], "publisher-loc": ["Oxford"], "publisher-name": ["Oxford University Press"], "fpage": ["234"], "lpage": ["244"]}, {"label": ["35"], "element-citation": ["\n"], "surname": ["Falconer", "Mackay"], "given-names": ["DS", "TFC"], "year": ["1996"], "article-title": ["Introduction to quantitative genetics."], "source": ["Harlow"], "publisher-loc": ["England"], "publisher-name": ["Prentice Hall"]}, {"label": ["37"], "element-citation": ["\n"], "surname": ["Staten", "Schully", "Noor"], "given-names": ["R", "SD", "MAF"], "year": ["2004"], "article-title": ["A microsatellite linkage map of Drosophila mojavensis."], "source": ["Bmc Genetics"], "volume": ["5"]}, {"label": ["40"], "element-citation": ["\n"], "surname": ["Zouros"], "given-names": ["E"], "year": ["1981"], "article-title": ["An autosome-Y chromosome combination that causes sterility in D. mojavensis and D. arizonensis hybrids."], "source": ["Drosophila Information Service"], "volume": ["56"], "fpage": ["167"], "lpage": ["168"]}, {"label": ["41"], "element-citation": ["\n"], "surname": ["Zouros", "Young"], "given-names": ["E", "JPW"], "year": ["1991"], "article-title": ["Searching for speciaton genes in the species pair Drosophila mojavnesis and D. arizonae."], "source": ["Molecular Techniques in Taxonomy"], "publisher-loc": ["Berlin"], "publisher-name": ["Springer-Verlag"], "fpage": ["233"], "lpage": ["247"]}, {"label": ["44"], "element-citation": ["\n"], "surname": ["Lander", "Abrahamson", "Barlow", "Daly", "Lincoln"], "given-names": ["E", "J", "A", "M", "S"], "year": ["1987"], "article-title": ["Mapmaker a Computer Package for Constructing Genetic-Linkage Maps."], "source": ["Cytogenetics and Cell Genetics"], "volume": ["46"], "fpage": ["642"], "lpage": ["642"]}, {"label": ["45"], "element-citation": ["\n"], "surname": ["Basten", "Weir", "Zeng", "Smith", "Gavora", "Benkel", "Chesnais", "Fairfull"], "given-names": ["CJ", "BS", "ZB", "C", "JS", "B", "J", "W"], "year": ["1994"], "article-title": ["Zmap-a QTL cartographer."], "source": ["Proceedings of the 5th world congress on genetics applied to livestock production:Computing strategies and software"], "publisher-loc": ["Guelph, Ontario, Canada"], "publisher-name": ["Organizing committee, 5th world congress on genetics applied to livestock production"], "fpage": ["65"], "lpage": ["66"]}, {"label": ["46"], "element-citation": ["\n"], "surname": ["Basten", "Weir", "Zeng"], "given-names": ["CJ", "BS", "ZB"], "year": ["2002"], "article-title": ["QTL Cartographer. 1.16 ed."], "publisher-loc": ["Raleigh, NC"], "publisher-name": ["Department of Statistics, North Carolina State Univeristy"]}, {"label": ["50"], "element-citation": ["\n"], "surname": ["Mclean", "Sanders", "Stroup"], "given-names": ["RA", "WL", "WW"], "year": ["1991"], "article-title": ["A Unified Approach to Mixed Linear-Models."], "source": ["American Statistician"], "volume": ["45"], "fpage": ["54"], "lpage": ["64"]}, {"label": ["51"], "element-citation": ["\n"], "surname": ["Lynch", "Walsh"], "given-names": ["M", "B"], "year": ["1998"], "article-title": ["Genetics and analysis of quantitative traits."], "publisher-loc": ["Sunderland, MA"], "publisher-name": ["Sinauer Associates Inc"]}, {"label": ["54"], "element-citation": ["\n"], "surname": ["Beavis"], "given-names": ["W"], "year": ["1994"], "article-title": ["The power and deceit of QTL experiments: lessons from comparative QTL studies."], "source": ["Report of the 49th annual corn and sorghum reseach conference"], "publisher-loc": ["Washington, DC"], "publisher-name": ["American Seed Trade Association"], "fpage": ["252"], "lpage": ["268"]}, {"label": ["56"], "element-citation": ["\n"], "surname": ["Fisher", "Bailey"], "given-names": ["RA", "NTJ"], "year": ["1949"], "article-title": ["The Estimation of Linkage with Differential Viability."], "source": ["Heredity"], "volume": ["3"], "fpage": ["215"], "lpage": ["228"]}, {"label": ["59"], "element-citation": ["\n"], "surname": ["Pontecorvo"], "given-names": ["G"], "year": ["1943"], "article-title": ["Viability interactions between chromosomes of Drosophila melanogaster and Drosophila simulans."], "source": ["Journal of Genetics"], "volume": ["45"], "fpage": ["51"], "lpage": ["66"]}, {"label": ["65"], "element-citation": ["\n"], "surname": ["Patterson", "Stone"], "given-names": ["JT", "WS"], "year": ["1952"], "article-title": ["Evolution in the Genus Drosophila."], "publisher-loc": ["New York"], "publisher-name": ["Macmillian"]}, {"label": ["66"], "element-citation": ["\n"], "surname": ["Kimura"], "given-names": ["M"], "year": ["1987"], "article-title": ["Habitat differentiation and speciaiton in the "], "italic": ["Drosophila auraria"], "source": ["Kontyu"], "volume": ["55"], "fpage": ["429"], "lpage": ["436"]}, {"label": ["68"], "element-citation": ["\n"], "surname": ["Mainland"], "given-names": ["GB"], "year": ["1942"], "article-title": ["Genetic relationships in the "], "italic": ["Drosophila funebris"], "source": ["University of Texas Publications"], "volume": ["4228"], "fpage": ["74"], "lpage": ["112"]}, {"label": ["70"], "element-citation": ["\n"], "surname": ["Crow"], "given-names": ["JF"], "year": ["1942"], "article-title": ["Cross fertility and isolating mechanisms in the "], "italic": ["Drosophila mulleri"], "source": ["University of Texas Publications"], "volume": ["4228"], "fpage": ["53"], "lpage": ["67"]}, {"label": ["73"], "element-citation": ["\n"], "surname": ["Shuker", "Underwood", "King", "Butlin"], "given-names": ["DM", "K", "TM", "RK"], "year": ["2005"], "article-title": ["Patterns of male sterility in a grasshopper hybrid zone imply accumulation of hybrid incompatibilities without selection."], "source": ["Proceedings of the Royal Society B-Biological Sciences"], "volume": ["272"], "fpage": ["2491"], "lpage": ["2497"]}]
PLoS One. 2008 Aug 27; 3(8):e3076