Skip to main content

Main menu

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home
  • Log in
  • My Cart

Advanced Search

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
Research Article

Dimensionality of mate choice, sexual isolation, and speciation

Paul A. Hohenlohe and Stevan J. Arnold
  1. Department of Zoology, Oregon State University, Corvallis, OR 97331-2914

See allHide authors and affiliations

PNAS September 21, 2010 107 (38) 16583-16588; https://doi.org/10.1073/pnas.1003537107
Paul A. Hohenlohe
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: hohenlop@science.oregonstate.edu
Stevan J. Arnold
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  1. Edited by Mark W. Blows, University of Queensland, Brisbane, Queensland, Australia, and accepted by the Editorial Board August 12, 2010 (received for review March 17, 2010)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Abstract

Multiple cues, across multiple sensory modalities, are involved in mate choice in a wide range of animal taxa. This multiplicity leads to the prediction that, in adaptive radiations, sexual isolation results from divergence in multiple dimensions. However, difficulties in directly measuring preferences and detecting multiple effects limit our ability to empirically assess the number of independent traits contributing to mate choice and sexual isolation. We present an approach to estimate the dimensionality of sexual isolation using mating trials across groups of related populations. We analyze nine radiations: seven in fruit flies (Drosophila) and one each in salamanders (Desmognathus) and cichlid fishes (Pseudotropheus). We find strong evidence that multiple latent traits—linear combinations of phenotypic traits and preferences—are responsible for the patterns of sexual isolation in all nine radiations but that dimensionality has a strong upper limit. Just two latent traits are implicated in the majority of cases. Mapping along latent trait axes tests predictions of sexual-selection models and allows correlation with specific phenotypic traits and functional components of mate choice. We find support for the role of stabilizing natural selection on the sexually selected (male) traits. In the cichlids, latent-trait axes incorporate male-coloration patterns and exhibit convergence as well as divergence among populations. In the salamanders, temporal patterning in sensory modalities and male vs. female preferences are reflected in different latent-trait axes.

  • behavioral isolation
  • female preference
  • multidimensional scaling
  • prezygotic isolation
  • sexual selection

Sexual isolation arising from disruptions in sexual communication is central to the process of speciation (1, 2), and such disruptions may arise from trait divergence driven by sexual selection (3, 4). Mate choice can involve a wide range of sensory modalities (chemical, auditory, visual, tactile, etc.), and multitrait, multimodal sexual communication is apparent across many animal taxa (5–8). Even within a single sensory modality (e.g., acoustic signaling in anurans), temporal patterning of signals is often complex and multivariate (7, 9). A few direct measurements of sexual selection have shown multivariate mate choice (10–16).

Unfortunately, most studies of mate choice focus on just one or two traits because of difficulties in detecting multivariate effects and directly measuring mate preference (7, 17). Furthermore, we lack a way to estimate the total number of functionally independent traits involved—the dimensionality of mate choice. This issue is analogous to the dimensionality of genetic variation for phenotypic traits; whereas genetic variation may be observed in many single traits, the prevalence of genetic correlation constrains the response to selection and phenotypic divergence (14, 18–21). We do not know to what extent the dimensionality of mate choice constrains the evolution of sexual isolation and speciation by sexual selection.

Here, we describe a method for determining the dimensionality of mate choice and sexual isolation. We use data from surveys of sexual isolation rather than direct examination of candidate traits. Sexual (behavioral) isolation is the degree to which the incidence of mating between populations is less than that within populations, which is assessed in mating trials (22). Ideally, multiple populations are surveyed so that sexual isolation is measured in all or most of the possible pair-wise combinations.

Our analysis builds on quantitative genetic models that describe the joint evolution of male traits and female preference, typically assuming both to be polygenic and normally distributed (23). Previous work extended such models to assess the effect of a single linear combination of traits on sexual isolation among populations (24, 25). We refer to such a combination as a latent trait. Here, we allow any number of independent latent traits (Fig. 1). Mating probability between males and females is a function of the distance between male traits and female preferences along all of the latent-trait axes (or equivalently, female traits and male preferences, some combination of both, or any set of quantitative traits for which similarity between males and females determines probability of mating). We use information-criterion statistics to determine the number of latent traits that best explains the observed pattern of sexual isolation and mate choice within and among populations in nine radiations across three taxa (24, 26–31).

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

Schematic diagram of the mate-choice model for two populations (A and B) in two latent-trait dimensions, with population means (solid circles) and within-population variation (dashed circles) for male traits (blue; z) and female preferences (red; y). The probability of mating between, for example, a female chosen randomly from population B (star 1) and a male chosen from either the same population (star 2) or population A (star 3) is a function of the distance between them (black dotted lines) in this phenotypic space.

Results

Dimensionality of Sexual Isolation.

We estimated the dimensionality of sexual isolation in nine radiations—seven from Drosophila and one each from species complexes in the salamander genus Desmognathus and the cichlid fish genus Pseudotropheus. All nine radiations exhibit divergence along multiple latent-trait dimensions (Table 1). However, there is strong evidence for an upper limit to dimensionality. In all cases, the best-fit number of dimensions [i.e., with the lowest Akaike Information Criterion (AIC) score] is much less than the maximum; all but two of the datasets require just two latent-trait dimensions. The remaining datasets require three and four latent-trait dimensions, respectively, and in the latter case (D. prosaltans), the AIC scores are nearly indistinguishable between three and four dimensions. In addition, both the log-likelihood values and the effective number of dimensions nD (19) reach a plateau at a dimensionality d well below the maximum (Table 1). Thus, as the points are placed in higher-dimensional spaces with d approaching the number of populations, the best-fit arrangement still occupies a subspace of dimensionality less than d. In fact, nD does not exceed three for any radiation and shows no correlation with the number of populations sampled (Spearman ρ = 0.085, P = 0.83) (32).

View this table:
  • View inline
  • View popup
Table 1.

Dimensionality of sexual isolation for nine datasets

Models of Sexual Selection.

Patterns of divergence along latent-trait axes (Fig. 2) test some predictions of sexual-selection models. Important factors are natural selection on the male trait(s), typically modeled as a Gaussian fitness function with width ω2, and strength of female mating preference, a Gaussian function with width v2 (23, 33). In seven of nine datasets, female latent-trait means exhibit greater divergence across populations than male means [measured as the trace of the within-sex, among-population (co)variance matrix] (Fig. 2 C and G). The model by Lande (33) predicts greater divergence among female means (equation 20 in ref. 33), and this effect is more pronounced if stabilizing natural selection on the male trait is strong relative to female preference (i.e., ω2 is small relative to v2). Thus, our results suggest that natural selection generally constrains the male traits.

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Best-fit male (blue) and female (red) population means for latent traits underlying sexual isolation in nine radiations. Black lines link male and female means for each population. For all but A and G, these two axes represent the best-fit number of dimensions; for A and G, these are the first two of three or four axes, respectively. Canonical correlation-redundancy indices (bottom right of each plot) show the proportion of female variance explained by males of the same population (red) and male variance explained by females (blue). (A) Desmognathus (24); dashed lines enclose populations in three currently recognized species. (B) D. paulistorum, Amazonian race (26). (C) D. paulistorum, Andean–South Brazilian race (26). (D) D. paulistorum, multiple races (26). (E) D. sturtevanti (27). (F) D. willistoni (28). (G) D. prosaltans (29). (H) D. auraria complex (31). (I) Pseudotropheus (30); labels indicate populations.

The datasets differ in relative amounts of among-population versus within-population (between-sex) divergence; redundancy indices report the proportion of variance in each sex that is explained by the other sex (Fig. 2). The model by Lande (33) predicts (equation 10 in ref. 33) that, at equilibrium, the expected distance between the male and female mean should be roughly proportional to v2/ω2 and the deviation of the male mean from its optimum. This result suggests that, in some taxa (Fig. 2H), in which the means are close together, female preferences are strong (small v2), and/or stabilizing natural selection on the male traits is weak (large ω2). In other groups (Fig. 2C), female preferences are apparently weak (large v2), and/or natural selection on the male traits is strong (small ω2).

Some models of sexual selection predict an equilibrium deviation between male and female values within populations (24, 25, 33, 34). Across related species lying along an equilibrium line or plane, consistent directional natural selection on the male trait(s) would produce a consistent direction of deviation between males and females, maintained by a tension between sexual selection and natural selection on the male trait. We found scant evidence for this prediction; the mean within-population male to female vector differed significantly from zero in only one of nine radiations (Fig. 2H) (Hotelling's t2 = 113.8, P = 0.03).

Testing Other Hypotheses.

Divergence along latent-trait axes mirrors taxonomic assignment of populations in some cases. The largest redundancy indices occur in a dataset with multiple species (Fig. 2H), illustrating among-species divergence. In the Desmognathus group, which was split into three species (35) after mating trials were conducted, the latent-trait values cluster the populations concordant with taxonomic revision (Fig. 2A). In contrast, in Pseudotropheus, our results support previous inferences of parallel evolution and a lack of phylogenetic signal in traits involved in mate choice; phylogenetic analysis places the Chisumulu (CH) population outside a clade containing the other four populations (36), contrary to the clustering along latent-trait axes (Fig. 2I). The datasets differ widely in the distribution of divergence among latent-trait axes, such that the primary axis accounts for anywhere from 41.6% to 81.6% of the total variance across traits (Fig. S1). They also differ substantially in the degree to which males and females diverge along separate axes (Fig. 2H versus 2C and Fig. S1).

Linking Latent Traits to Phenotypes.

Latent-trait mapping in Pseudotropheus indicates some correlates with male nuptial coloration (Fig. 2I). Our results support a hypothesis (30) that two populations, Nkhata Bay (NB) and CH, show a low level of sexual isolation because of similar coloration in the males, evidenced by the males’ similar latent-trait values, despite phylogenetic distance. Two discrete male coloration traits vary conspicuously among these five populations (36). Yellow chest occurs only in Mara Rocks (MR), and this difference explains 98.8% of the variance among males along latent-trait axis 1. In contrast, the orange dorsal phenotype occurs only in Mphanga Rocks emmiltos (EM), accounting for minimal proportions of variance in either latent trait.

In Desmognathus, courtship is a multiphase process, and termination at different phases may reflect the contribution of different sensory modalities or functional components of mate choice to overall sexual isolation. We assessed the contribution of these factors to the three latent traits using correlation analysis on male–female distances along each axis and four behavioral measures of sexual isolation (37): the proportions of mating trials ending with no sexual activity, ending after the pursuit phase, and ending after the persuasion phase, and the time taken to reach the persuasion phase in pairs that reached it (Table 2). One correlation remained highly significant after correction for false discovery rate (FDR) (38)—that between latent trait 1 and termination after the pursuit phase (Spearman ρ = 0.569, P = 0.0008). The correlation between latent trait 2 and time to persuasion was elevated but not significant after FDR correction because of a smaller sample size (ρ = 0.496, P = 0.028).

View this table:
  • View inline
  • View popup
Table 2.

Correlation of latent traits with phases of courtship in Desmognathus

Discussion

Dimensionality of Mate Choice and Sexual Isolation.

We found consistent support for the view that multiple cues are used in mate choice, but we also found that only a small number of independent latent traits are actually responsible for sexual isolation among populations in a radiation. Each latent trait may be a linear combination of many actual phenotypic traits, which are strongly correlated with each other across populations in their contribution to sexual isolation.

A central assumption in this analysis is that the linear combinations of traits and preferences that make up the latent traits remain constant across the radiation. A few lines of evidence support this assumption and the conclusion that divergence may proceed repeatedly along a limited number of latent-trait axes. Radiation among similar environments may lead to parallel evolution in mate-choice traits, either as a byproduct of divergent selection on other traits or because of direct environmental correlates with mate-choice factors (39–42). Male nuptial coloration patterns have evolved repeatedly in Lake Malawi cichlids, including the populations analyzed here (36). Hypotheses for constraints on divergence in nuptial coloration include the physical aspects of color vision in the lake and feeding ecology (36). Correlation among mate-choice traits may also arise, because multiple cues are indicators of a smaller number of underlying axes, such as body condition of prospective mates (43–47). Finally, strong genetic correlations among traits may constrain divergence by both sexual selection and genetic drift by limiting the dimensionality of the G matrix of additive genetic variation. In Drosophila species, such strong correlations have been observed among both cuticular hydrocarbon (CHC) traits in males and female preferences for those traits (16, 48).

Implicit in the analysis is that the phenotypic covariance structure among traits and preferences within males and females, as well as the covariance structure of the multivariate preference functions, is constant across populations within a radiation (Methods). Comparative empirical data are limited, although similarity in the genetic covariance structure of CHCs has been observed across populations in Drosophila (49, 50). We conducted numerical simulations to test the effect of this assumption of constant covariance structure on estimates of dimensionality. Variable covariance structure creates variation in, but does not systematically bias, the estimated dimensionality of sexual isolation. Thus, the upper limit to dimensionality of sexual isolation among populations is robust to relaxation of this assumption (SI Text and Fig. S2).

Finally, our analysis assumed unimodal preference functions (Methods). The functional form of individual preferences is difficult to estimate (7), but a few studies suggest that unimodal preference functions are common in nature (9, 23, 51–53). Also observed are open-ended preference functions (7, 54). Further analysis is needed to test whether open-ended preferences result in substantially different estimates of dimensionality. Better still, a comprehensive analysis would allow multiple functional forms of preference to co-occur for different traits (13, 55).

From Latent Traits to Specific Phenotypes.

Our method does not identify traits responsible for sexual isolation; however, it suggests testable hypotheses about specific phenotypes. After populations are mapped onto latent-trait space, correlation with traits, preferences, and functional components of the mate-choice process can be assessed. In many organisms, observational and experimental studies provide lists of candidate traits. Because of the difficulties in assessing preferences directly (7), this correlational exercise may be most fruitful using more easily measured (usually male) phenotypic traits.

In plethodontid salamanders such as the Desmognathus species, courtship consists of several distinct phases: male identification and pursuit of the female, a persuasion phase during which the male delivers pheromones to the female, and spermatophore deposition and sperm transfer (34, 56). Different sensory modalities are used: vision during the pursuit phase, chemical communication during pursuit and persuasion, and tactual communication during sperm transfer (57). Breakdown in sexual communication at different phases implicates different traits and preferences as well as different sensory modalities in sexual isolation. For instance, breakdown of mating after the pursuit phase suggests incompatibility between male pheromones and their perception by females. This aspect of isolation correlated strongly with latent trait 1. Thus, latent trait 1 may primarily reflect the role of chemical communication in sexual isolation. We also found a weak relationship between divergence along latent trait 2 and time to reach the persuasion phase. During this period, the male uses visual and chemical cues to evaluate a female as a potential mate. Latent trait 2 may incorporate female traits assessed during this period (e.g., size and pheromones) and male perception of them. These relationships suggest that mate choice and sexual isolation may be somewhat modular—that evolution of different functional clusters of traits and preferences, grouped both by sensory modality and temporal phase during courtship, may occur independently.

Male coloration plays an important role in courtship and sexual isolation in the Lake Malawi cichlid fishes (36, 58). Our analysis is consistent with the role of male coloration patterns. We also confirmed that populations can diverge as well as converge along trait axes (36). Despite the myriad of visual, behavioral, auditory, and olfactory traits apparently involved in sexual selection in these taxa (59–61), divergence among populations is limited to just a few combinations of traits. In our analysis, this limit is manifest as a relatively small dimensionality; in a phylogenetic approach (36), it is apparent that replicate populations converge on similar phenotypes, and the number of male coloration patterns is surprisingly limited.

Limits to Dimensionality of Sexual Isolation.

In all nine radiations, we observed an apparent upper bound to the dimensionality of sexual isolation. These populations likely diverge across a large number of phenotypic traits or preferences, and these traits together form a high-dimensional space. Within this space, however, the populations occupy a much smaller-dimensional subspace. This result may find an explanation in sexual-selection models. A general prediction is that the dimensionality of the equilibrium equals the number of genetically independent traits minus the number of selective constraints (62, 63). For example, when two traits (e.g., a male trait and female preference for that trait) evolve but only one is under selection, the expected equilibrium is a one-dimensional line (33). At any point on the equilibrium, the forces of natural and sexual selection balance, and the population drifts freely (25, 33). The dispersion of populations (Fig. 2) may reflect stochastic differentiation on a two-, three-, or four-dimensional equilibrium. Under this interpretation, the latent traits that we identified represent combinations of multiple traits and preferences, and the correlations among these traits are maintained by the forces of natural and sexual selection leading to equilibrium. Alternatively, if all traits and preferences are subject to direct natural selection, the equilibrium reduces to a point (64). In this case, divergence among populations could be the result of movement of the natural-selective optima for the traits. The pattern of movement of these optima would determine the dimensionality of sexual isolation. Comparative analysis of natural selection on key traits may be able to distinguish these scenarios.

Our finding that multiple traits are involved in sexual isolation but that the number of functionally independent traits is limited carries several implications for understanding speciation driven by sexual selection. First, it reinforces the empirical strategy of testing multiple traits in studies of mate choice, sexual selection, and speciation. Second, because quantitative genetic models of sexual selection have almost always focused on a single male trait and a single female preference (23, 65, 66), these models have likely underestimated the rates at which ornamental traits and sexual isolation can evolve. We need theoretical models that explore the consequences of multivariate sexual isolation to test this supposition. Third, claiming that particular traits are important to speciation simply because trait differentiation coincides with species boundaries leaves much to be desired. Surveys of sexual isolation, combined with a multivariate study of candidate traits, could determine which traits actually confer sexual isolation and facilitate speciation.

Finally, a longstanding question in evolution is how genetic variation is maintained in the presence of strong selection. One resolution to this issue lies in a multivariate perspective: genetic correlation and correlational selection mean that genetic variance can persist in many single traits, even while genetic variation is limited and evolution is constrained in the multivariate direction in which selection acts (18, 21). This view led Walsh and Blows (21) to conclude that “trait-by-trait explanations of the natural world are doomed to fail” (21). Analogously, we found that, while mate choice may involve a large number of traits, sexual isolation is constrained to just a few dimensions. The clarion call of Walsh and Blows (21) for a multivariate view of genetic variation applies equally to understanding mate choice and speciation by sexual selection.

Methods

Mating Trial Datasets.

We used published data on mating trials in three taxa: salamanders, fruit flies, and cichlid fishes. We modeled mating probability as a function of distance between male and female trait values in a multidimensional space (Fig. 1). Fitting points to male and female population means relies on a matrix of pair-wise isolation estimates. Our method is most robust when this matrix is complete or nearly so (i.e., when females of each population are tested against males of each population, including their own, and vice versa). In some cases, we analyzed only a subset of published data to provide a dense matrix (extracted datasets in Tables S1–S9). Experimental design differs among mating studies, described below, and is implicit in the mating-probability functions. Conclusions about sexual isolation may be sensitive to experimental design (no-choice, male-choice, etc.); however, we found no indication of bias among the three experimental designs represented here, and mating trial results seem broadly comparable across designs (22). Differences in experimental design did necessitate alternative statistical models: binomial and Poisson.

Arnold et al. (24) used a no-choice design in the D. ochrophaeus complex in the southeastern United States. A single male and single female were paired for a standard period, and mating (insemination) was scored as present or absent. We used data on nine populations, now grouped as three species: D. carolinensis, D. ocoee, and D. orestes (35).

Several researchers have assessed sexual isolation in Drosophila species using a male-choice design, pairing males of one population with females from their own and another population. Mating is assessed as presence or absence of insemination in dissected females. From a massive dataset (26) on populations of D. paulistorum, we extracted three subsets: four populations in the Amazonian race, eight populations in the Andean–South Brazilian race, and seven populations from across six races. We used data on three other New World species: five populations of D. sturtevanti (27), five populations of D. willistoni (28), and seven populations of D. prosaltans (29). Finally, we analyzed data on four species in the Asian D. auraria complex (31).

Knight and Turner (30) used a multiple-choice experimental design among five populations of cichlid fish in the P. zebra species complex (taxonomy of these populations remains in flux) (30, 60). Males and females of all five populations were placed together in replicate enclosures, and the total number of matings of each pair-wise combination was recorded over a standardized period. Population labels are from ref. 30: CH, MR, NB, EM, and RW (Ruarwe).

Model.

We begin with a univariate model of sexual isolation (24). Let z be a latent male trait (linear combination of traits with invariant coefficients) and y be a latent female preference. Assume that y and z are normally distributed within each population with means Embedded Image and Embedded Image and variances Embedded Image and Embedded Image, respectively. Following an absolute-preference model (or equivalently, a relative-preference model where the female trait is Embedded Image) (33), probability of mating between a female y and a male z is proportional to a Gaussian function with deviation parameter (z – y) and width parameter Embedded Image. Thus, the probability of mating between a female randomly chosen from population i and a male from population j is Embedded Image (24).

We extend this model to d traits so that πij is now a multivariate Gaussian function of the Euclidean distance between trait values. We assume that the within-population covariance among male traits and among female preferences, as well as the covariance structure of the multivariate preference function, is constant across populations (mathematical details and a test of this assumption in SI Text). This assumption allows us to scale and rotate the axes of trait space by the combined (co)variance structure so that z and y are column vectors of independent traits and preferences, respectively. We assume that cij, the probability of mating when zj = yi, takes a constant value c. Then, the probability of mating between a female and a male randomly chosen from populations i and j, respectively, is Embedded Image, where T represents the vector transpose. For the salamander and fruit fly datasets, we assume that mating trials are independent and use a binomial model. The probability of observing kij matings in nij trials is a binomial function of πij. The log likelihood of the complete set of population mean-trait values Embedded Image and Embedded Image is (Eq. 1):

Embedded Image

For the cichlids, we use a Poisson model for the total number of matings kij over the experimental period. The probability distribution of kij is Poisson with parameter Embedded Image, where p is the expected number of matings when zj = yi, again assumed constant across population pairs. The log likelihood of the complete set of trait values is (Eq. 2)

Embedded Image

Determining Dimensionality.

Given a number of dimensions d and mating-trial results, finding the maximum-likelihood values of Z and Y is analogous to metric multidimensional unfolding (67) but with a cost function given by Eqs. 1 or 2. Because the likelihood surfaces are relatively rugged (more so at lower dimensions), we alternated a multivariate Newton–Raphson algorithm and gradient hill climbing with systematic swapping of points to find the global best fit. For each dataset at each dimension, we also estimated the maximum-likelihood value of c (by iterative hill climbing) or p (with an explicit solution) for the binomial or Poisson models, respectively. The best-fit points at dimension d provided the seed for the dimension d + 1 (computer code to conduct this analysis is available from P.A.H.).

We determined the optimal d for each dataset using the corrected AIC (cAIC) (68). For d dimensions and m populations, Z and Y represent a total of 2dm coordinates. However, the likelihood of a set of points depends only on the pair-wise distances between them (Eqs. 1 and 2), and therefore, there are Embedded Image free parameters, plus one for either c or p. This value is given in Table 1. For m populations, in principle, a maximum of d = 2m − 1 dimensions may be required to fit all 2m points in Z and Y. However, the lowest AIC score was always achieved at fewer than m dimensions, and above this number, the best-fit configuration of points did not change substantially. Thus, we present results only up to dimension m − 1 for each dataset. At each d, we calculated the effective number of dimensions Embedded Image, where Embedded Image is the ith eigenvalue of the (co)variance matrix of all male and female population means (19).

For each best-fit set of points Z and Y, we centered and rotated the axes analogous to a principal components analysis so that the first latent-trait axis represents the maximum variance among all trait values. Plots and all further analyses were done in this rotated space. We used canonical correlation-redundancy analysis (32) between Z and Y to assess within- vs. among-population divergence. To test for consistent deviation between males and females of each population, we calculated Hotelling's t2 statistic for the m deviation vectors (yi – zi) in each dataset. Significance was assessed by evaluating (m − d)t2/d(m – 1) against an F(d, m – d) distribution.

Acknowledgments

We thank members of the Houck–Arnold laboratory and two anonymous reviewers for very helpful and insightful comments. This work was supported by National Science Foundation Grant IOS-0843392 (to P.A.H. and S.J.A.).

Footnotes

  • 1To whom correspondence should be addressed. E-mail: hohenlop{at}science.oregonstate.edu.
  • Author contributions: P.A.H. and S.J.A. designed research; P.A.H. performed research; P.A.H. contributed new reagents/analytic tools; P.A.H. analyzed data; and P.A.H. and S.J.A. wrote the paper.

  • The authors declare no conflict of interest.

  • This article is a PNAS Direct Submission. M.W.B. is a guest editor invited by the Editorial Board.

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1003537107/-/DCSupplemental.

    References

    1. ↵
      1. Coyne JA,
      2. Orr HA
      (2004) Speciation (Sinauer, Sunderland, MA).
    2. ↵
      1. Mayr E
      (1942) Systematics and the Origin of Species, from the Viewpoint of a Zoologist (Columbia University Press, New York).
    3. ↵
      1. Losos JB
      (2009) Ecology and Adaptive Radiation of Anoles (University of California Press, Berkeley, CA).
    4. ↵
      1. Panhuis TM,
      2. Butlin RK,
      3. Zuk M,
      4. Tregenza T
      (2001) Sexual selection and speciation. Trends Ecol Evol 16:364–371.
      OpenUrlCrossRefPubMed
    5. ↵
      1. Andersson M
      (1994) Sexual Selection (Princeton University Press, Princeton).
    6. ↵
      1. Candolin U
      (2003) The use of multiple cues in mate choice. Biol Rev Camb Philos Soc 78:575–595.
      OpenUrlPubMed
    7. ↵
      1. Chenoweth SF,
      2. Blows MW
      (2006) Dissecting the complex genetic basis of mate choice. Nat Rev Genet 7:681–692.
      OpenUrlCrossRefPubMed
    8. ↵
      1. Jennions MD,
      2. Petrie M
      (1997) Variation in mate choice and mating preferences: A review of causes and consequences. Biol Rev Camb Philos Soc 72:283–327.
      OpenUrlPubMed
    9. ↵
      1. Gerhardt HC,
      2. Huber F
      (2002) Acoustic Communication in Insects and Anurans: Common Problems and Diverse Solutions (University of Chicago Press, Chicago).
    10. ↵
      1. Bentsen CL,
      2. Hunt J,
      3. Jennions MD,
      4. Brooks R
      (2006) Complex multivariate sexual selection on male acoustic signaling in a wild population of Teleogryllus commodus. Am Nat 167:E102–E116.
      OpenUrlCrossRefPubMed
    11. ↵
      1. Blows MW,
      2. Brooks R,
      3. Kraft PG
      (2003) Exploring complex fitness surfaces: Multiple ornamentation and polymorphism in male guppies. Evolution 57:1622–1630.
      OpenUrlCrossRefPubMed
    12. ↵
      1. Blows MW,
      2. Chenoweth SF,
      3. Hine E
      (2004) Orientation of the genetic variance-covariance matrix and the fitness surface for multiple male sexually selected traits. Am Nat 163:329–340.
      OpenUrlCrossRefPubMed
    13. ↵
      1. Brooks R,
      2. et al.
      (2005) Experimental evidence for multivariate stabilizing sexual selection. Evolution 59:871–880.
      OpenUrlCrossRefPubMed
    14. ↵
      1. Chenoweth SF,
      2. Blows MW
      (2008) Q(St) meets the G matrix: The dimensionality of adaptive divergence in multiple correlated quantitative traits. Evolution 62:1437–1449.
      OpenUrlCrossRefPubMed
    15. ↵
      1. Gerhardt HC,
      2. Brooks R
      (2009) Experimental analysis of multivariate female choice in gray treefrogs (Hyla versicolor): Evidence for directional and stabilizing selection. Evolution 63:2504–2512.
      OpenUrlCrossRefPubMed
    16. ↵
      1. McGuigan K,
      2. Van Homrigh A,
      3. Blows MW
      (2008) Genetic analysis of female preference functions as function-valued traits. Am Nat 172:194–202.
      OpenUrlCrossRefPubMed
    17. ↵
      1. Kingsolver JG,
      2. et al.
      (2001) The strength of phenotypic selection in natural populations. Am Nat 157:245–261.
      OpenUrlCrossRefPubMed
    18. ↵
      1. Hansen TF,
      2. Houle D
      (2008) Measuring and comparing evolvability and constraint in multivariate characters. J Evol Biol 21:1201–1219.
      OpenUrlCrossRefPubMed
    19. ↵
      1. Kirkpatrick M
      (2009) Patterns of quantitative genetic variation in multiple dimensions. Genetica 136:271–284.
      OpenUrlCrossRefPubMed
    20. ↵
      1. Mezey JG,
      2. Houle D
      (2005) The dimensionality of genetic variation for wing shape in Drosophila melanogaster. Evolution 59:1027–1038.
      OpenUrlCrossRefPubMed
    21. ↵
      1. Walsh B,
      2. Blows MW
      (2009) Abundant genetic variation + strong selection = multivariate genetic constraints: A geometric view of adaptation. Annu Rev Ecol Evol Syst 40:41–59.
      OpenUrlCrossRef
    22. ↵
      1. Coyne JA,
      2. Orr HA
      (1989) Patterns of speciation in Drosophila. Evolution 43:362–381.
      OpenUrlCrossRef
    23. ↵
      1. Mead LS,
      2. Arnold SJ
      (2004) Quantitative genetic models of sexual selection. Trends Ecol Evol 19:264–271.
      OpenUrlCrossRefPubMed
    24. ↵
      1. Arnold SJ,
      2. Verrell PA,
      3. Tilley SG
      (1996) The evolution of asymmetry in sexual isolation: A model and a test case. Evolution 50:1024–1033.
      OpenUrlCrossRef
    25. ↵
      1. Uyeda JC,
      2. Arnold SJ,
      3. Hohenlohe PA,
      4. Mead LS
      (2009) Drift promotes speciation by sexual selection. Evolution 63:583–594.
      OpenUrlCrossRefPubMed
    26. ↵
      1. Carmody G,
      2. et al.
      (1962) Mating preferences and sexual isolation within and between the incipient species of Drosophila paulistorum. Am Midl Nat 68:67–82.
      OpenUrlCrossRef
    27. ↵
      1. Dobzhansky T
      (1944) Experiments on sexual isolation in Drosophila: III. Geographic strains of Drosophila sturtevanti. Proc Natl Acad Sci USA 30:335–339.
      OpenUrlFREE Full Text
    28. ↵
      1. Dobzhansky T,
      2. Mayr E
      (1944) Experiments on sexual isolation in Drosophila: I. Geographic strains of Drosophila willistoni. Proc Natl Acad Sci USA 30:238–244.
      OpenUrlFREE Full Text
    29. ↵
      1. Dobzhansky T,
      2. Streisinger G
      (1944) Experiments on sexual isolation in Drosophila: II. Geographic strains of Drosophila prosaltans. Proc Natl Acad Sci USA 30:340–345.
      OpenUrlFREE Full Text
    30. ↵
      1. Knight ME,
      2. Turner GF
      (2004) Laboratory mating trials indicate incipient speciation by sexual selection among populations of the cichlid fish Pseudotropheus zebra from Lake Malawi. Proc Biol Sci 271:675–680.
      OpenUrlAbstract/FREE Full Text
    31. ↵
      1. Kurokawa H,
      2. Oguma Y,
      3. Tachibana N
      (1982) Sexual isolation among four species of D. auraria complex. Drosoph Inf Serv 58:98–99.
      OpenUrl
    32. ↵
      1. Team RDC
      (2007) R: A Language and Environment for Statistical Computing, Version 2.6.0 (Foundation for Statistical Computing, Vienna).
    33. ↵
      1. Lande R
      (1981) Models of speciation by sexual selection on polygenic traits. Proc Natl Acad Sci USA 78:3721–3725.
      OpenUrlAbstract/FREE Full Text
    34. ↵
      1. Mead LS,
      2. Verrell PA
      (2002) Evolution of courtship behaviour patterns and reproductive isolation in the Desmognathus ochrophaeus complex. Ethology 108:403–427.
      OpenUrlCrossRef
    35. ↵
      1. Tilley SG,
      2. Mahoney MJ
      (1996) Patterns of genetic differentiation in salamanders of the Desmognathus ochrophaeus complex (Amphibia: Plethodontidae) Herpet Monogr 10:1–42.
      OpenUrlCrossRef
    36. ↵
      1. Allender CJ,
      2. Seehausen O,
      3. Knight ME,
      4. Turner GF,
      5. Maclean N
      (2003) Divergent selection during speciation of Lake Malawi cichlid fishes inferred from parallel radiations in nuptial coloration. Proc Natl Acad Sci USA 100:14074–14079.
      OpenUrlAbstract/FREE Full Text
    37. ↵
      1. Verrell PA,
      2. Arnold SJ
      (1989) Behavioral observations of sexual isolation among allopatric populations of the mountain dusky salamander, Desmognathus ochrophaeus. Evolution 43:745–755.
      OpenUrlCrossRef
    38. ↵
      1. Benjamini Y,
      2. Hochberg Y
      (1995) Controlling the false discovery rate: A practical and powerful approach to multiple testing. J R Statist Soc B 57:289–300.
      OpenUrl
    39. ↵
      1. McKinnon JS,
      2. et al.
      (2004) Evidence for ecology's role in speciation. Nature 429:294–298.
      OpenUrlCrossRefPubMed
    40. ↵
      1. Nosil P,
      2. Crespi BJ,
      3. Sandoval CP
      (2002) Host-plant adaptation drives the parallel evolution of reproductive isolation. Nature 417:440–443.
      OpenUrlCrossRefPubMed
    41. ↵
      1. Rundle HD,
      2. Chenoweth SF,
      3. Doughty P,
      4. Blows MW
      (2005) Divergent selection and the evolution of signal traits and mating preferences. PLoS Biol 3:e368.
      OpenUrlCrossRefPubMed
    42. ↵
      1. Rundle HD,
      2. Nagel L,
      3. Wenrick Boughman J,
      4. Schluter D
      (2000) Natural selection and parallel speciation in sympatric sticklebacks. Science 287:306–308.
      OpenUrlAbstract/FREE Full Text
    43. ↵
      1. David P,
      2. Bjorksten T,
      3. Fowler K,
      4. Pomiankowski A
      (2000) Condition-dependent signalling of genetic variation in stalk-eyed flies. Nature 406:186–188.
      OpenUrlCrossRefPubMed
    44. ↵
      1. Hine E,
      2. Chenoweth SF,
      3. Blows MW
      (2004) Multivariate quantitative genetics and the lek paradox: Genetic variance in male sexually selected traits of Drosophila serrata under field conditions. Evolution 58:2754–2762.
      OpenUrlCrossRefPubMed
    45. ↵
      1. Kotiaho JS,
      2. Simmons LW,
      3. Tomkins JL
      (2001) Towards a resolution of the lek paradox. Nature 410:684–686.
      OpenUrlCrossRefPubMed
    46. ↵
      1. Petfield D,
      2. Chenoweth SF,
      3. Rundle HD,
      4. Blows MW
      (2005) Genetic variance in female condition predicts indirect genetic variance in male sexual display traits. Proc Natl Acad Sci USA 102:6045–6050.
      OpenUrlAbstract/FREE Full Text
    47. ↵
      1. Proulx SR,
      2. Servedio MR
      (2009) Dissecting selection on female mating preferences during secondary contact. Evolution 63:2031–2046.
      OpenUrlCrossRefPubMed
    48. ↵
      1. Hine E,
      2. Blows MW
      (2006) Determining the effective dimensionality of the genetic variance-covariance matrix. Genetics 173:1135–1144.
      OpenUrlCrossRefPubMed
    49. ↵
      1. Blows MW,
      2. Higgie M
      (2003) Genetic constraints on the evolution of mate recognition under natural selection. Am Nat 161:240–253.
      OpenUrlCrossRefPubMed
    50. ↵
      1. Higgie M,
      2. Blows MW
      (2007) Are traits that experience reinforcement also under sexual selection? Am Nat 170:409–420.
      OpenUrlCrossRefPubMed
    51. ↵
      1. Ritchie MG
      (1996) The shape of female mating preferences. Proc Natl Acad Sci USA 93:14628–14631.
      OpenUrlAbstract/FREE Full Text
    52. ↵
      1. Shaw KL
      (2000) Interspecific genetics of mate recognition: Inheritance of female acoustic preference in Hawaiian crickets. Evolution 54:1303–1312.
      OpenUrlCrossRefPubMed
    53. ↵
      1. Sinervo B,
      2. Zamudio KR
      (2001) The evolution of alternative reproductive strategies: Fitness differential, heritability, and genetic correlation between the sexes. J Hered 92:198–205.
      OpenUrlAbstract/FREE Full Text
    54. ↵
      1. Boake CRB,
      2. DeAngelis MP,
      3. Andreadis DK
      (1997) Is sexual selection and species recognition a continuum? Mating behavior of the stalk-eyed fly Drosophila heteroneura. Proc Natl Acad Sci USA 94:12442–12445.
      OpenUrlAbstract/FREE Full Text
    55. ↵
      1. Ritchie MG,
      2. Saarikettu M,
      3. Livingstone S,
      4. Hoikkala A
      (2001) Characterization of female preference functions for Drosophila montana courtship song and a test of the temperature coupling hypothesis. Evolution 55:721–727.
      OpenUrlCrossRefPubMed
    56. ↵
      1. Sever DM
      1. Houck LD,
      2. Arnold SJ
      (2003) in Phylogeny and Reproductive Biology of Urodela (Amphibia) ed Sever DM (Science Publishers, Enfield, NH), pp 383–424.
    57. ↵
      1. Houck LD,
      2. Verrell PA
      (1993) Studies of courtship behavior in plethodontid salamanders: A review. Herpetologica 49:175–184.
      OpenUrl
    58. ↵
      1. Genner MJ,
      2. Turner GF
      (2005) The mbuna cichlids of Lake Malawi: A model for rapid speciation and adaptive radiation. Fish Fish 6:1–34.
      OpenUrl
    59. ↵
      1. Amorim MCP,
      2. Simões JM,
      3. Fonseca PJ,
      4. Turner GF
      (2008) Species differences in courtship acoustic signals among five Lake Malawi cichlid species (Pseudotropheus spp.) J Fish Biol 72:1355–1368.
      OpenUrlCrossRef
    60. ↵
      1. Blais J,
      2. et al.
      (2009) Assortative mating among Lake Malawi cichlid fish populations is not simply predictable from male nuptial colour. BMC Evol Biol 9:53.
      OpenUrlCrossRefPubMed
    61. ↵
      1. Plenderleith M,
      2. van Oosterhout C,
      3. Robinson RL,
      4. Turner GF
      (2005) Female preference for conspecific males based on olfactory cues in a Lake Malawi cichlid fish. Biol Lett 1:411–414.
      OpenUrlAbstract/FREE Full Text
    62. ↵
      1. Kirkpatrick M
      (1986) The handicap mechanism of sexual selection does not work. Am Nat 127:222–240.
      OpenUrlCrossRef
    63. ↵
      1. Lande R,
      2. Arnold SJ
      (1985) Evolution of mating preference and sexual dimorphism. J Theor Biol 117:651–664.
      OpenUrlCrossRefPubMed
    64. ↵
      1. Pomiankowski A,
      2. Iwasa Y,
      3. Nee S
      (1991) The evolution of costly mate preferences I. Fisher and biased mutation. Evolution 45:1422–1430.
      OpenUrlCrossRef
    65. ↵
      1. Iwasa Y,
      2. Pomiankowski A
      (1994) The evolution of mate preferences for multiple sexual ornaments. Evolution 48:853–867.
      OpenUrlCrossRef
    66. ↵
      1. Pomiankowski A,
      2. Iwasa Y
      (1993) Evolution of multiple sexual preferences by Fisher's runaway process of sexual selection. Proc R Soc London, Ser B 253:173–181.
      OpenUrlAbstract/FREE Full Text
    67. ↵
      1. Cox TF,
      2. Cox MAA
      (2001) Multidimensional Scaling (Chapman & Hall/CRC, Boca Raton, FL).
    68. ↵
      1. Burnham KP,
      2. Anderson DR
      (1998) Model Selection and Inference: A Practical Information-Theoretic Approach (Springer, New York).
    PreviousNext
    Back to top
    Article Alerts
    Email Article

    Thank you for your interest in spreading the word on PNAS.

    NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

    Enter multiple addresses on separate lines or separate them with commas.
    Dimensionality of mate choice, sexual isolation, and speciation
    (Your Name) has sent you a message from PNAS
    (Your Name) thought you would like to see the PNAS web site.
    CAPTCHA
    This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
    Citation Tools
    Dimensionality of mate choice, sexual isolation, and speciation
    Paul A. Hohenlohe, Stevan J. Arnold
    Proceedings of the National Academy of Sciences Sep 2010, 107 (38) 16583-16588; DOI: 10.1073/pnas.1003537107

    Citation Manager Formats

    • BibTeX
    • Bookends
    • EasyBib
    • EndNote (tagged)
    • EndNote 8 (xml)
    • Medlars
    • Mendeley
    • Papers
    • RefWorks Tagged
    • Ref Manager
    • RIS
    • Zotero
    Request Permissions
    Share
    Dimensionality of mate choice, sexual isolation, and speciation
    Paul A. Hohenlohe, Stevan J. Arnold
    Proceedings of the National Academy of Sciences Sep 2010, 107 (38) 16583-16588; DOI: 10.1073/pnas.1003537107
    del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
    • Tweet Widget
    • Facebook Like
    • Mendeley logo Mendeley

    Article Classifications

    • Biological Sciences
    • Evolution
    Proceedings of the National Academy of Sciences: 107 (38)
    Table of Contents

    Submit

    Sign up for Article Alerts

    Jump to section

    • Article
      • Abstract
      • Results
      • Discussion
      • Methods
      • Acknowledgments
      • Footnotes
      • References
    • Figures & SI
    • Info & Metrics
    • PDF

    You May Also be Interested in

    Water from a faucet fills a glass.
    News Feature: How “forever chemicals” might impair the immune system
    Researchers are exploring whether these ubiquitous fluorinated molecules might worsen infections or hamper vaccine effectiveness.
    Image credit: Shutterstock/Dmitry Naumov.
    Reflection of clouds in the still waters of Mono Lake in California.
    Inner Workings: Making headway with the mysteries of life’s origins
    Recent experiments and simulations are starting to answer some fundamental questions about how life came to be.
    Image credit: Shutterstock/Radoslaw Lecyk.
    Cave in coastal Kenya with tree growing in the middle.
    Journal Club: Small, sharp blades mark shift from Middle to Later Stone Age in coastal Kenya
    Archaeologists have long tried to define the transition between the two time periods.
    Image credit: Ceri Shipton.
    Illustration of groups of people chatting
    Exploring the length of human conversations
    Adam Mastroianni and Daniel Gilbert explore why conversations almost never end when people want them to.
    Listen
    Past PodcastsSubscribe
    Panda bear hanging in a tree
    How horse manure helps giant pandas tolerate cold
    A study finds that giant pandas roll in horse manure to increase their cold tolerance.
    Image credit: Fuwen Wei.

    Similar Articles

    Site Logo
    Powered by HighWire
    • Submit Manuscript
    • Twitter
    • Facebook
    • RSS Feeds
    • Email Alerts

    Articles

    • Current Issue
    • Special Feature Articles – Most Recent
    • List of Issues

    PNAS Portals

    • Anthropology
    • Chemistry
    • Classics
    • Front Matter
    • Physics
    • Sustainability Science
    • Teaching Resources

    Information

    • Authors
    • Editorial Board
    • Reviewers
    • Subscribers
    • Librarians
    • Press
    • Cozzarelli Prize
    • Site Map
    • PNAS Updates
    • FAQs
    • Accessibility Statement
    • Rights & Permissions
    • About
    • Contact

    Feedback    Privacy/Legal

    Copyright © 2021 National Academy of Sciences. Online ISSN 1091-6490