Tests of the role of sexual selection in genitalic divergence with mutiple hybrid clines.
Determining whether sexually selected characters affect reproductive isolation is key to understanding the signifcance of such selectively driven divergence in speciation. Divergence in the shape of male genitalia is ubiquitous in many insect taxa and can have important fitness consequences for males within a species. Yet, it is unclear whether this selectively divergent character affects gene flow among species. We test this hypothesis using explicit predictions about the clinal transition in genitalic shape across a hybrid zone between the grasshopper subspecies Barytettix humphreysii humphreysii and B. h. cochisei. Slight morphological differences in genitalic shape across their ranges raised the possibility that local processes might give misleading interpretations about sexual selection's potential contribution to reproductive isolation. Therefore, we examined multiple transects to explore the effectiveness of genitalic divergence as a barrier to gene flow on a background of potentially different extrinsic (e.g., environmental conditions) or intrinsic (e.g., genetic makeup) factors. If sexual selection is the predominant factor maintaining the hybrid zone, the shape of the clines between transects is expected to be concordant, thereby overriding potential dispersal differences associated with varying environmental and ecological conditions or variation in population densities that would cause cline shape to vary among localities. Furthermore, if selection against hybrids with intermediate genitalic morphologies does indeed represent a strong barrier to gene flow, the morphological transition between the adjacent subspecies is predicted to be quite steep. These predictions were supported by the geometric morphometric analyses, suggesting that genitalic divergence plays a direct role in reproductive isolation and that sexual selection plays a predominant role in the maintenance of the subspecies differences. Thus, the results of our study provide compelling evidence that postmating-prezygotic interactions can indeed limit gene flow, and consequently may play an important role in speciation. Using this framework, we suggest future studies that can address a number of remaining questions about the nature of selection, and ultimately, how sexual selection is operating in these taxa.
cryptic female choice, geometric morphometrics, introgression, postmating-prezygotic interactions, reproductive isolation
The shape of character clines, as well as the degree of concor-dance and coincidence of different characters across hybrid zones, can provide valuable insights into the nature of barriers to gene exchange (Barton & Hewitt 1985, 1989, Harrison 1990). Yet, sexually-selected traits are surprisingly rarely studied in hybrid zones (Butlin & Neems 1994), despite the potential of this context to reveal whether evolutionary divergence of such characters has important consequences for reproductive isolation, as predicted by theoretical work (e.g., Lande 1981, 1982, Schluter & Price 1993, Gavrilets 2003).
In a hybrid zone the steepness of character transitions between species can inform about the strength of selective forces separating the taxa (Barton & Gale 1993). The frequently observed concordance of multiple characters across hybrid zones has emphasized the role of selective pressures arising from genetic associations between alleles, or breakdown of coadapted-gene complexes, where the movement of alleles is impeded by unfavorable, or inviable, genetic backgrounds (Barton & Hewitt 1985, 1989, Harrison 1990). With such intrinsic barriers to gene flow, character transitions between the adjacent species arise as a result of the strong linkage disequilibrium generated by dispersal into the zone--that is, a steep cline reflects the strong net selective pressure of the whole genome. However, like the effects of disruptive-ecological selection, an effective barrier to gene flow can be generated if hybrids have a strong mating disadvantage (e.g., Stratton & Uetz 1986, Vamosi & Schluter 1999, Cruz et al. 2001, Naisbit et al. 2001). Such factors may play a much more important role than the traditional criteria for reproductive isolation--hybrid inviability and sterility (Coyne & Orr 1989)--which are frequently absent in recently diverged taxa suggesting that such postzygotic effects are the by-products of genetic divergence rather than a driving force in speciation.
In this study we examine two clines in a morphological character--the shape of the male genitalia--and we compare cline shape between transects to examine the role sexual selection might play in the evolution of this character and in reproductively isolating the grasshopper subspecies, Barytettix humphreysii humphreysii and B. h. cochisei (Orthoptera: Acrididae: Melanoplinae). While differences in genitalic morphology in insects is widely believed to reflect divergence driven by sexual selection (Eberhard 1996, Eberhard & Cordero 2003, House & Simmons 2003, Hosken & Stockley 2004), what, if any, consequence such divergence has for speciation is not clear (e.g., Oneal & Knowles 2012). By studying divergence in the genitalic character across a hybrid zone, we gain insight into how strong selection might be (i.e., would genitalic divergence pose a significant barrier via a mechanism of strong selection against hybrid individuals with intermediate genitalic phenotypes) and whether it predominates (as measured by similar cline shapes across independent transects).
The hybrid zone between B. h. humphreysii and B. h. cochisei extends from southern Arizona to northern Sonora, Mexico (Cohn & Cantrall 1974). These grasshoppers are flightless and do not sing or exhibit obvious courtship displays or ecological differences. They are morphologically quite similar--the primary character distinguishing them is differences in the male genitalia, specifically the aedeagus, or intromittent organ. Slight morphological differences in the genitalia across the range of each subspecies (Cohn & Cantrall 1974) raises the possibility that interpretations about general hybrid zone maintenance may be misled by local processes (Butlin 1994, Shapiro 1998), including different extrinsic (e.g., environmental conditions) or intrinsic (e.g., genetic makeup) factors (Hairston et al. 1992, Buerkle & Rieseberg 2001, Morgan-Richards & Wallis 2003, Smadja et al. 2004). However, by studying multiple transects across the hybrid zone, we can explore whether selection is the primary factor maintaining the zone, and consequently the effectiveness of genitalic divergence as a barrier to gene flow. If sexual selection is the predominant factor maintaining the hybrid zone, we predict that the shape of the clines between transects to be concordant, thereby overriding potential dispersal differences associated with varying environmental and ecological conditions or variation in population densities that would cause cline shape to vary among localities (e.g., Szymura & Barton 1991, Butlin et al. 1991). Furthermore, if selection against hybrids with intermediate genitalic morphologies does indeed represent a strong barrier to gene flow, then we predict the morphological transition between the adjacent subspecies to be quite steep. Results from morphometric analyses are indeed consistent with these predictions, suggesting that genitalic divergence plays a direct role in reproductive isolation and that sexual selection plays a predominant role in the maintenance of the subspecies differences, although our work does not identify the mechanism by which sexual selection might act (see Eberhard 1996). However, there are a number of remaining questions about the nature of this selection (i.e., we cannot rule out selection on other correlated characters, or that the steepness of the cline is unique to genitalic characters, for example). Nevertheless, using these findings as a framework, which is informative without invoking a specific mechanism about how sexual selection might operate, we discuss additional studies that could potentially reveal how sexual selection is operating in these taxa.
Material and methods
Samples.--Adult male specimens were collected from two transects across the hybrid zone between grasshopper populations of the subspecies Barytettix h. humphreysii and B. h. cochisei (Orthoptera: Acrididae: Melanoplinae) in southern Arizona (Appendix 1). 106 and 230 museum specimens (Museum of Zoology, Insect Division, University of Michigan) were examined from two transects of 18.2 km and 22.6 km in length, which we refer to as the "Bisbee Junction" and "Don Luis" transects, respectively; specimens are deposited in the collection with genitalia extruded for digitization. The transects are approximately 4.2 km apart and oriented in an east-westerly direction (see Fig. 4 for orientation of specimens across hybrid zone). Additional collections of non-hybrid individuals were also made from two sites about 50 km from the transects (Appendix 1). Patterns of morphological variation were quantified using morphometric analysis of outlines of the genitalic structures (e.g., Arnqvist 1998). In the context of clinal variation, morphological shape can be analyzed in a manner similar to the method of plotting allele frequency as a function of position along a hybrid transect to characterize cline shape across a hybrid zone. Note that while we focus on a specific trait, this does not preclude that other characters, including those of the genitalic complex, might be subject to selection (see Marquez & Knowles 2007). Whether such characters do or do not show strong clinal variation is reserved for future investigations and we highlight what such findings would imply for our results (see Discussion).
Collection of morphometric outlines.--Photographs of the extruded male genitalia taken from a dorsal perspective were digitized (16-bit grayscale). Before the shape of specimen could be quantified, each image was processed according to a standardized procedure to produce silhouettes that tracked the margin of the right ventral valve of the aedeagus (Fig. 1). To produce a closed outline so that the shape of the genitalia could be described using an elliptical Fourier analysis (Ferson et al. 1985), the outlines were closed with a half circle at the junction between the aedeagal sheath and the aedeagus itself--a clearly identifable and consistent landmark (see also Marquez & Knowles 2007). The program SHAPE (Iwata & Ukai 2002) was used to conduct the elliptic Fourier analysis. Fourier analyses were made invariant of size, position and rotation (i.e., the first three harmonics were not included in the statistical analyses, Rohlf & Archie 1984) using the "Chc2Nef" component of SHAPE. For each specimen, 11 harmonics were included in analyses, where each harmonic yields 4 variables or Fourier coefficients that were analyzed in a principal-components analysis (Liu et al. 1996, Rohlf & Archie 1984) using the "PrinComp" module of SHAPE; ordination of the specimens indicated that the first five principal components, each of which accounted for 1% or more of the total variation of each transect, captured 96% of the variance among the specimens (scores: PC1 = 82%, PC2 = 7%, PC3 = 4%, PC4 = 2%, and PC5 = 1%). Further analyses were therefore based on the first five principal components of shape variation of each sample. We note that elliptical Fourier analysis is effective for characterizing shape variation (see Marquez & Knowles 2007), especially when landmarks are difficult to identify (although see Oneal & Knowles 2012 for an alternative based on semi-landmarks).
Statistical Analysis.--To quantify shape differences among individuals collected along a transect, a measure of morphological distance between each specimen and a parental species was used and calculated separately for each transect. Specifically, the distance D was computed as the Euclidean distance between each specimen and the parental (i.e., non-hybrid) species at the eastern end of its respective transect (i.e., B. h. cochisei that is geographically adjacent to the study transect), based on principal component scores, according to the equation:
D= [[(PC[1.sub.parental] - PC[1.sub.hybrid]).sup.2] + [(PC[2.sub.parental]--PC[2.sub.hybrid]).sup.2] + ... + [[(PC[n.sub.parental]--PC[n.sub.hybrid]).sup.2 ].sup. 1\2],
for the n selected principal components. Finally, D was normalized so that all distances were in the interval [0,1], thus providing a scale-free representation of each specimen in a common multivariate-morphometric space, allowing the direct comparison of the two transects in terms of morphological distance.
Regression of D against corresponding geographical distances between each specimen and the western end of the transect, normalized to the [0, 1] interval, was used to characterize the geographic distribution of morphological change along the transects. Prior to this regression, a logit transformation was applied to morphological distances (Barton & Gale 1993). To compare the clines of the two transects, slopes from the linear portion of the regression of normalized D against normalized geographical distances (e.g., Szymura & Barton 1986), were first calculated. Then each transect was bootstrapped (N = 10,000) to obtain standard errors and 95% confidence intervals for each of the transect slopes. To test for significant differences between slopes, data from the two transects were pooled and 10,000 bootstrap subsamples with the same size as the original samples were obtained from the pooled dataset and their slopes calculated as described above. The absolute difference between slopes was computed for each pair of bootstrap subsamples to produce an empirical distribution of the expected difference between slopes under the null hypothesis of clinal homogeneity across transects. Observed absolute slope difference was then compared to this distribution. An observed value larger than 95% of the bootstrapped differences would reject the hypothesis that both transects belong to clines with equal slopes.
Additionally, principal component scores were plotted against each other for each transect to identifly trends in the relationships of the principal components among localities. Analyses were carried out in Matlab v. 6.5 (Mathworks 2003). Compiled versions of the software employed are freely available upon request.
Principal-components analysis indicates that not only is the variation among specimens captured by the elliptical-Fourier analysis, but also that the character transition between B. h. humphreysii and B. h. cochisei can be precisely described in this quantitative framework. The shape of the male genitalia of the parental subspecies was similar between transects; average non-normalized morphological distance (D) was 0.241 ([+ or -] 0.019) and 0.210 ([+ or -] 0.007) for B. h. humphreysii, and 0.057 ([+ or -] 0.005) and 0.042 ([+ or -] 0.004) for B. h. cochisei, from the Don Luis and Bisbee transects, respectively. Scores of hybrid individuals from the two transects widely overlap in the first two principal components of the pooled data, which account for 84.5% and 7.0% of total variance, respectively (Fig. 2), indicating that specimens from both transects share a large portion of their shape variation. Specific differences between the male genitalia of B. h. humphreysii and B. h. cochisei, as well as the shape of the hybrid individuals, quantified by the elliptical Fourier analysis (Fig. 3) matched the qualitative descriptions of Cohn & Cantrall (1974).
Analysis of the shape of the male genitalia across the hybrid zones revealed that both transects were characterized by sigmoid curves with a steep clinal transition between the two adjacent taxa (Fig. 4). The slope at the Bisbee Junction transect was 9.03 (95% confidence interval: 8.4635-9.4522) and did not differ significantly from the slope of the Don Luis transect, which was 9.37 (95% confidence interval: 8.8874-9.8979); the observed difference of 0.35 between the slopes is not large enough to reject the hypothesis that the transects' slopes were equal (P = 0.75) based on 10,000 bootstrap subsamples. Although we have a fairly large number of observations to calculate the slope (with more than 100 specimens per transect), it is possible that additional samples might provide power for detecting a difference in slope that was not detected here.
The steep character transition between B. h. humphreysii and B. h. cochisei is indicative of strong selection acting on the hybrids. Moreover, the smooth cline shape and concordance between the transects confirms that this selection dominates the dynamics of the hybrid zone, overriding the effects of other potential factors such as environmental and ecological differences or varying population densities between transects (Barton & Hewitt 1985). Because the steep cline was observed in the male genitalia, a character widely postulated to be under sexual selection (Eberhard 1996), the results have a number of intriguing implications.
Even without specifying the mechanisms of selection, not only are the results consistent with the hypothesized direct role of male-genitalic shape differences in reproductive isolation, but the study also provides evidence for a predominant role of sexual selection in the maintenance of species differences (because we are studying a character involved in reproductive interactions, sexual rather than natural selection is the focus of our conclusions; see Eberhard 1994). As such, the work complements other work that aims to understand the drivers of genitalic divergence and the role of genitalic divergence in speciation. However, a number of unanswered questions about the nature of selection remain. The framework provided by these museum specimens highlights that future investigations could elucidate longstanding issues regarding the genetic basis of reproductive isolation, and ultimately how sexual selection may be involved in the isolation of species (Jiggins et al. 1996).
Role of genitalic divergence in reproductive isolation.--Cline widths vary in response to the levels of selection a particular character experiences (Barton & Gale 1993). If the shape of the male genitalia was not related to hybrid fitness, then a broad cline is expected. Moreover, lack of significant difference in slopes of character transition across transects implies that selection against hybrids is independent of environmental or ecological conditions (Buerkle & Rieseberg 2001). Otherwise, introgression patterns would differ between transects (e.g., Butlin et al. 1991, Szymura & Barton 1991).
The implications of this steep clinal transition in the shape of the male genitalia across the hybrid zone are particularly intriguing. Unlike other sexually selected traits, genitalic characters are not subject to other selective influences so they can provide a clear window into sexual selection without being confounded by other factors. This contrasts with cases where sexual selection affecting divergence and reproductive isolation is mediated by ecological divergence--that is, the selection on the mating characters directly affecting reproductive isolation as a pleiotropic effect of ecological change (e.g., in Heliconius butterflies, Naisbit et al. 2001; in Littorina snails, Cruz et al. 2001; in wolf spiders, Stratton & Uetz 1986). For example, divergent habitat preferences in sticklebacks bias encounter rates between hybrids and parental species producing the conditions for sexual selection to contribute to reproductive isolation because mating success is tied to habitat use (Hatfield & Schluter 1996, Vamosi & Schluter 1999).
Evidence of sexual selection acting on hybrids of B. h. humphreysii and B. h. cochisei is also uniquely significant; sexual selection in this case is operating on a postmating-prezygotic character. In studies that have shown reproductive isolation is directly affected by sexual selection, gene flow between the hybridizing taxa is impeded by premating characters (e.g., Brumfield et al. 2001, McDonald et al. 2 0 01 , Bailey et al. 2003, Bronson et al. 2003, Smadja et al. 2004). Attention to postmating-prezygotic characters, including male genitalia (e.g., Arnqvist 1998, Eberhard & Cordero 2003), and their potential role in species divergence has increased recently, with the primary focus on antagonistic interactions between the sexes (e.g., Wolfner 2002, Civetta & Clark 2000, Knowles & Markow 2001, Swanson et al. 2001). Yet, whether selection acting on postmating-prezygotic traits is strong enough to produce reproductive isolation remains largely unknown (Knowles et al. 2004a, b). The results of our study provide compelling evidence that postmating-prezygotic interactions can indeed limit gene flow, and consequently may play an important role in speciation.
The nature of selection.--How selection is acting on the male genitalia is an interesting question that requires further study. Despite the conspicuous and fairly ubiquitous pattern of genitalic divergence among insect taxa (Eberhard 1996), surprisingly little is known about how selection might operate (e.g., Oneal & Knowles 2015). Similar to sexual selection on premating characters, reduced fitness of hybrids could reflect female choice (e.g., Saetre et al. 1997, Babik et al. 2003, Smadja et al. 2004), as well as male-male interactions (e.g., Brumfield et al. 2001, McDonald et al. 2001, Bronson et al. 2003). For example, the assortative fertilization associated with differences in the shape of the male genitalia could result from the reduced success of sperm storage of the hybrid males, possibly reflecting lowered stimulation of females to store sperm or misplacements of spermatophores by the males in the female genital tract (Eberhard 1996). Understanding how selection operates against hybrids could also provide insights into the driving forces underlying the divergence of male-genitalic characters.
While our results are consistent with models of speciation that posit a direct role for divergent sexual selection (e.g., Lande 1981, 1982, Schluter & Price 1993, Gavrilets 2003), important aspects about the nature of this selection still need to be addressed. A steep cline may not necessarily reflect the strength of selection acting on any single character alone (Gavrilets & Cruzan 1998). Selection on other characters that are in linkage disequilibrium with the character under study can contribute to the shape of the cline (Szymura & Barton 1991, Barton & Gale 1993, Barton & Baird 1996, Bridle et al. 2001). Consequently, to determine sexual selection's specific contribution to reproductive isolation (i.e., the strength of selection operating against hybrids based on differences in male-genitalic morphology), the potential contribution of the net selective pressure on the genome needs to be considered (e.g., Dasmahapatra et al. 2002). This study provides a framework in which the clinal transitions of other characters, and specifically those not under sexual selection, can be studied to distinguish these alternative interpretations about the nature of selection. If strong selection against hybrids of intermediate genitalic morphology constitutes the primary factor affecting reproductive isolation, then coincident clines of multiple characters are not expected since the fitness associated with each character across the hybrid zone is likely to differ (Barton & Hewitt 1985, 1989, Harrison 1990, Butlin et al. 1991). A steeper cline is also predicted for the male genitalia compared to non-sexually selected traits reflecting the different intensities of selection impeding the movement of alleles between the adjacent parental taxa (e.g., Szymura & Barton 1991, Harrison & Bogdanowicz 1997, Brumfield et al. 2001, Babik et al. 2003). Further study of this hybrid zone will not only reveal the extent to which genitalic differences, and therefore sexual selection, contribute to reproductive isolation, but will also contribute interesting insights into longstanding questions regarding the evolution of reproductive isolation, and specifically the role of selectively driven divergence (Wu 2001, Gray & Cade 2000, Vines et al. 2003).
We would like to express our utmost gratitude to Theodore J. Cohn who provided the initial specimens and images that formed the basis for this study, as well as thoughtful contributions regarding the hybrid zone. Thanks to members of the Knowles lab for their input. The research was part of TMC's undergraduate honors project and was funded by a grant from the University of Michigan to LLK.
Arnqvist G. 1998. Comparative evidence for the evolution of genitalia by sexual selection. Nature 393: 784-786.
Babik W., Szymura J.M., Rafnski J. 2003. Nuclear markers, mitochondrial DNA and male secondary sexual traits variation in a newt hybrid zone (Triturus vulgaris' T. montandoni). Molecular Ecology 12:1913-1930.
Bailey R.I., Thomas C.D., Butlin R.K. 2003. Premating barriers to gene exchange and their implications for the structure of a mosaic hybrid zone between Chorthippus brunneus and C. jacobsi (Orthoptera: Acrididae). Journal of Evolutionary Biology 17:108-119.
Barton N.H., Hewitt G.M. 1985. Analysis of hybrid zones. Annual Review of Ecology and Systematics 16:113-148.
Barton N.H., Hewitt G.M. 1989. Adaptation, speciation and hybrid zones. Nature 341: 497-503.
Barton N.H., Gale K.S. 1993. Genetic Analysis of Hybrid Zones, pp. 13-14. In: Harrison R.G., (Ed.) Hybrid Zones and the Evolutionary Process. Oxford University Press, New York.
Barton N.H., Baird S.J.E. 1996. Development of software for the analysis of geographic variation and hybrid zones. University of Edinburgh, Edinburgh, U.K. Available via http:/helios.bto.ed.ac.uk/evolgen/.
Bridle J.R., Baird S.J.E., Butlin R.K. 2001. Spatial structure and habitat variation in a grasshopper hybrid zone. Evolution 55: 1832-1843.
Bronson C.L., Grubb T.C. Jr., Sattler G.D., Braun M.J. 2003. Mate preference: a possible causal mechanism for a moving hybrid zone. Animal Behaviour 65: 489-500.
Brumfield R.T., Jernigan R.W., McDonald D.B., Braun M.J. 2001. Evolutionary implications of divergent clines in an avian (Manacus: Aves) hybrid zone. Evolution 55: 2070-2087.
Buerkle C.A., Rieseberg L.H. 2001. Low intraspecific variation for genomic isolation between hybridizing sunflower species. Evolution 55: 84-691.
Butlin R.K. 1994. Genetic variation in mating signals and responses. Pp. 327-366. In: Lambert D.M., Spencer H.G. (Eds) Speciation and the recognition concept. Johns Hopkins University Press, Baltimore, London.
Butlin R.K., Neems R.M. 1994. Hybrid zones and sexual selection: Science 265: 122.
Butlin R.K., Richie M.G., Hewitt G.M. 1991. Comparisons among morphological characters and between localitites in the Chorthippus parallelus hybrid zone (Orthoptera: Acrididae). Philosophical Transactions of the Royal Society of London, B 334: 297-308.
Civetta A., Clark A.G. 2000. Correlated effects of sperm competition and postmating female mortality. Proceedings of the National Academy of Sciences, USA 97: 13162-13165.
Cruz R., Rolan-Alvarez E., Garcia C. 2001. Sexual selection on phenotypic traits in a hybrid zone of Littorina saxatilis (Olivi). Journal of Evolutionary Biology 14: 773-785.
Cohn T.J., Cantrall I.J. 1974. Variation and Speciation in the Grasshoppers of the Conalcaeini (Orthoptera: Acrididae: Melanoplinae): The Lowland Forms of Western Mexico, the Genus Barytettix. San Diego Society of Natural History Memoirs 6:1-131.
Coyne, J. A. & Orr, H. A. 1989. Patterns of speciation in Drosophila. Evolution 43: 362-381.
Dasmahapatra K.K., Blum M.J., Aiello A., Hackwell S., Davies N., Bermingham E.P., Mallet J. 2002. Inferences from a rapidly moving hybrid zone. Evolution 56: 741-753.
Eberhard W.G. 1996. Female Control: Sexual Selection by Cryptic Female Choice. Princeton Univ. Press, New Jersey.
Eberhard W.G., Cordero C. 2003. Sexual conflict and female choice. Trends in Ecology and Evolution 18: 438-439.
Eberhard W.G. 2004. Male-female conflict and genitalia: failure to confirm predictions ininsects and spiders. Biological Reviews 79: 121-186.
Ferson S., Rohlf F.J., Koehn R.K. 1985. Measuring shape variation of two-dimensional outlines. Systematic Zoology 34: 59-68.
Gavrilets S. 2003. Models of speciation: what have we learned in 40 years? Evolution 57: 2197-2215.
Gavrilets S., Cruzan M.B. 1998. Neutral gene flow across single locus clines. Evolution 52: 1277-1284.
Gray D.A., Cade W.H. 2000. Sexual selection and speciation in field crickets. Proceedings of the National Academy of Sciences, USA 97: 14449-14454.
Hairston Sr. N.G., Wiley R.H., Smith C.K. 1992. The dynamics of two hybrid zones in Appalachian salamanders of the genus Plethodon. Evolution 46: 930-938.
Harrison R.G. 1990. Hybrid zones: windows on evolutionary process, pp. 69-128. In: Futuyma D., Antonovics J. (Eds) Oxford Surveys in Evolutionary Biology, vol. 7. Oxford University Press, New York.
Harrison R.G., Bogdanowicz S.M. 1997. Patterns of variation and linkage disequilibrium in a field cricket hybrid zone. Evolution 51: 493-505.
Hatfield T., Schluter D. 1996. A test for sexual selection on hybrids of two sympatric sticklebacks. Evolution 50: 2429-2434.
Hosken D.J., Stockley P. 2004. Sexual selection and genital evolution. Trends in Ecology and Evolution 19: 87-93.
House C.M., Simmons L.W. 2003. Genital morphology and fertilization success in the dung beetle Onthaphoagus taurus: an example of sexually selected male genitalia. Proceedings of the Royal Society of London, B 270: 447-455.
Iwata H., Ukai Y. 2002. SHAPE: A computer program package for quantitative evaluation of biological shapes based on elliptic Fourier descriptors. Journal of Heredity 93: 384-485.
Jiggins C.D., Mcmillan W.O., Neukirchen W., Mallet J. 1996. What can hybrid zones tell us about speciation? The case of Heliconius erato and H. himera (Lepidoptera: Nymphalidae). Biological Journal of the Linnean Society 59:221-242.
Knowles L.L., Markow T.A. 2001 Sexually antagonistic coevolution of a postmating-prezygotic reproductive character in desert Drosophila. Proceedings of the National Academy of Sciences, USA 98: 8692-8696.
Knowles L.L., Hernandez B.B., Markow T.A. 2004a. Exploring the consequences of postmating-prezygotic interactions between the sexes. Proceedings of the Royal Society of London, B 271:S357-S359.
Knowles L.L., Hernandez B.B., Markow T.A. 2004b. Non-antagonistic interactions between the sexes revealed by the ecological consequences of reproductive traits. Journal of Evolutionary Biology 18:156-161.
Lande R. 1981. Models of speciation by sexual selection on polygenic traits. Proceedings of the National Academy of Sciences, USA 78:3721-3725.
Lande R. 1982. Rapid origin of sexual isolation and character divergence in a cline. Evolution 36: 213-223.
Liu J., Mercer J.M., Stam L.F., Gibson G.C., Zeng Z.-B., Laurie CC. 1996. Genetic analysis of a morphological shape difference in the male genitalia of Drosophila simulans and D. mauritiana. Genetics 142: 1129-1145.
Marquez E.J., Knowles L.L. 2007. Correlated Evolution in multivariate traits: detecting codivergence across multiple dimensions. Journal of Evolutionary Biology 20:2334-2348. doi:10.1111/j.1420-9101.2007.01415x
Mathworks. 2003. Matlab version 6.5, release 13. The Mathworks, Inc., Natick, MA.
McDonald D.B., Clay R.P., Brumfield R.T., Braun M.J. 2001. Sexual selection on plumage and behavior in an avian hybrid zone: experimental tests of male-male interactions. Evolution 55: 1443-1451.
Morgan-Richards M., Wallis G.P. 2003. A comparison of five hybrid zones of the weta Hemideina thoracica (Orthoptera: Anostostomatidae): degree of cytogenetic differentiation fails to predict zone width. Evolution 57: 849-861.
Naisbit R.E., Jiggins C.D., Mallet J. 2001. Disruptive sexual selection against hybrids contributes to speciation between Heliconius cydno and Heliconius melpomene. Proceedings of the Royal Society of London, B 268: 1-6.
Oneal E., Knowles L.L. 2012. Ecological selection as the cause and sexual differentiation as the consequence of species divergence? Proceedings of the Royal Society of London, B 280: 20122236.
Oneal E., Knowles L.L. 2015. Paternity analyses in wild-caught and lab-reared Caribbean cricket females reveal the influence of mating environment on post-copulatory sexual selection. Journal of Evolutionary Biology 28: 2300-2307.
Rohlf F.J., Archie J.W. 1984. A comparison of Fourier methods for the description of wing shape in mosquitos (Diptera, Culcidae). Systematic Zoology 33: 302-317.
Saetre G-P., Kral M., Bures S. 1997. Differential species recognition abilities of males and females in a flycatcher hybrid zone. Journal of Avian Biology 28: 259-263.
Schluter D., Price T. 1993. Honesty, perception and population divergence in sexually selected traits. Proceedings of the Royal Society of London, B 253: 117-122.
Shapiro L.H. 1998. Hybridization and geographic variation in two meadow katydid contact zones. Evolution 52: 784-796.
Smadja, C., Catalan J., Ganem G. 2004. Strong premating divergence in a unimodal hybrid zone between two subspecies of the house mouse. Journal of Evolutionary Biology 17: 165-176.
Stratton G.E., Uetz G.W. 1986. The inheritance of courtship behavior and its role as a reproductive isolating mechanism in two species of Schizocosa wolf spiders (Araneae: Lycosidea). Evolution 40:129-141.
Swanson W.J., Clark A.G., Waldrip-Dali H.M., Wolfner M.F., Aquadro C.F. 2001. Evolutionary EST analysis identifes rapidly evolving male reproductive proteins in Drosophila. Proceedings of the National Academy of Sciences, USA 98: 7375-7379.
Szymura J.M., Barton N.H. 1991. The genetic structure of the hybrid zone between the fre-bellied toads Bombina bombina and B. variegata: comparisons between transects and between loci. Evolution 45: 237-261.
Vamosi S.M., Schluter D. 1999. Sexual selection against hybrids between sympatric stickleback species: evidence from a field experiment. Evolution 53: 874-879.
Vines T.H., Kohler S.C., Thiel M., Ghira I., Sands T. R., MacCallum C.J., Barton N.H., Nurnberger B. 2003. The maintenance of reproductive isolation in a mosaic hybrid zone between the fre-bellied toads Bombina bombina and B. variegata. Evolution 57: 1876-1888.
Wu C.-I. 2001. The genic view of the process of speciation. Journal of Evolutionary Biology 14:851-865.
Wolfner M.F. 2002. The gifts that keep on giving: physiological functions and evolutionary dynamics of male seminal proteins in Drosophila. Heredity. 88:85?93.
L. LACEY KNOWLES, THOMAS M. CHAPPELL, ELADIO J. MARQUEZ, THEODORE J. COHN [dagger]
Department of Ecology and Evolution, Museum of Zoology, University of Michigan, Ann Arbor, MI 48109-1079. Email: knowlesl@umich. edu
[dagger] deceased November 2012
Appendix 1. Specific locality information for specimens (Cochise Co., AZ).
specimen locality number label locality description 95-99 Bisbee5 1.5 mi. E. Bisbee Jct.. 100-118 Don1 12 mi. W. Don Luis (E. side of river) 0.6 mi. E. Palominas 119-147 Don2 11.7 mi. W. Don Luis (W. side of river), 0.9 E. Palominas 148-170 Don3 8.7 mi. W. Don Luis, 3.9 mi. E. Palominas 171-193 Don4 1.9 mi. W. Don Luis, 10.7 mi. E. Palominas 194-218 Don5 0.9 mi. W. Don Luis, 11.7 mi. E. Palominas 219-245 Don6 1.9 mi. NE. Don Luis 246-268 Don7 0.7 mi. NE. Don Luis 304-324 Don8 4.3 mi. SW. Don Luis 269-302 Don0 Don Luis (just N. Naco Rd. Jct.) 524-532 Bisbee1 6.3 mi. E. Bisbee (High Lonesome Rd.) 515-516 Bisbee2 4.8 mi. E. Bisbee Jct. 517-523 Bisbee3 2.6 mi. E. Bisbee Jct. 513-514 Bisbee4 2.2 mi. E. Bisbee Jct. 349-366 Bisbee5 1.5 mi. E. Bisbee Jct. 367-383 Bisbee6 2.0 mi. W. Bisbee Jct. 384-398 Bisbee7 3.5 mi. W. Bisbee Jct. 471-491 Bisbee8 3.3 mi. N. Bisbee Jct. 492-512 Bisbee0 Bisbee Jct. 533-537 Doug1 4.4 mi. W. Douglas 543-546 Doug2 27.2 mi. E. Douglas specimen number latitude longitude 95-99 31.35 -109.86 100-118 31.40 -110.11 119-147 31.40 -110.11 148-170 31.40 -110.06 171-193 31.40 -109.94 194-218 31.40 -109.93 219-245 31.42 -109.89 246-268 31.40 -109.90 304-324 31.40 -109.96 269-302 31.40 -109.91 524-532 31.34 -109.77 515-516 31.34 -109.80 517-523 31.34 -109.84 513-514 31.35 -109.84 349-366 31.35 -109.86 367-383 31.35 -109.91 384-398 31.35 -109.94 471-491 31.39 -109.88 492-512 31.35 -109.88 533-537 31.34 -109.61 543-546 31.34 -109.08