Reinforcement shapes clines in female mate discrimination in Drosophila subquinaria

Divergence in genetic (co)variances and the alignment of gmax with phenotypic divergence

To better understand the sources of biological diversity in nature, we need information on the mechanisms underlying population divergence. Biological systems with patterns of naturally occurring adaptive variation among populations can provide insight into the genetic architecture of diverging traits and the influence of genetic constraints on responses to selection. Using a system of reproductive character displacement in the North American mushroom-feeding fly Drosophila subquinaria, we assessed patterns of genetic (co)variance among a suite of chemical signaling traits and divergence in this pattern among populations. D. subquinaria exhibits stronger reproductive isolation against the closely related Drosophila recens in sympatry, where both female mating preferences and male chemical signaling traits have diverged from the ancestral allopatric populations. We collected 3 wild populations from each region and, in the lab, characterized the phenotypic divergence in these traits, as well as the additive genetic (co)variance structure (G-matrix), via replicate breeding designs. We found divergence between allopatric and sympatric D. subquinaria in the shape and size of the G-matrix, and that the leading axis of genetic variance (gmax) had changed in sympatry to come into alignment with the primary axis of phenotypic divergence between the sympatric and allopatric regions.

Postmating prezygotic isolation occurs at two levels of divergence in Drosophila recens and D. subquinaria

Identifying the presence and strength of reproductive isolating barriers is necessary to understand how species form and then remain distinct in the face of ongoing gene flow. Here, we study reproductive isolation at two stages of the speciation process in the closely related mushroom-feeding species Drosophila recens and Drosophila subquinaria. We assess 3 isolating barriers that occur after mating, including the number of eggs laid, the proportion of eggs laid that hatched, and the number of adult offspring from a single mating. First, all 3 reproductive barriers are present between D. recens females and D. subquinaria males, which are at the late stages of speciation but still produce fertile daughters through which gene flow can occur. There is no evidence for geographic variation in any of these traits, concurrent with patterns of behavioural isolation. Second, all 3 of these reproductive barriers are strong between geographically distant conspecific populations of D. subquinaria, which are in the early stages of speciation and show genetic differentiation and asymmetric behavioural discrimination. The reduction in the number of eggs laid is asymmetric, consistent with patterns in behavioural isolation, and suggests the evolution of postmating prezygotic isolation due to cascade reinforcement against mating with D. recens. In summary, not only may postmating prezygotic reproductive barriers help maintain isolation between D. recens and D. subquinaria, but they may also drive the earliest stages of isolation within D. subquinaria.

Sexual signal evolution and patterns of assortative mating across an intraspecific contact zone

Contact zones provide important insights into the evolutionary processes that underlie lineage divergence and speciation. Here, we use a contact zone to ascertain speciation potential in the red-eyed treefrog (Agalychnis callidryas), a brightly coloured and polymorphic frog that exhibits unusually high levels of intraspecific variation. Populations of A. callidryas differ in a number of traits, several of which are known sexual signals that mediate premating reproductive isolation in allopatric populations. Along the Caribbean coast of Costa Rica, a ~100 km contact zone, situated between two phenotypically and genetically divergent parent populations, contains multiple colour pattern phenotypes and late-generation hybrids. This contact zone provides the opportunity to examine processes that are important in the earliest stages of lineage divergence. We performed analyses of colour pattern variation in five contact zone sites and six parental sites and found complex, continuous colour variation along the contact zone. We found discordance between the geographic distribution of colour pattern and previously described genomic population structure. We then used a parental site and contact zone site to measure assortative mating and directional selection from naturally-occurring amplectant mating pairs. We found assortative mating in a parental population, but no assortative mating in the contact zone. Furthermore, we uncovered evidence of directional preference towards the adjacent parental phenotype in the contact zone population, but no directional preference in the parent population. Combined, these data provide insights into potential dynamics at the contact zone borders and indicate that incipient speciation between parent populations will be slowed.

The population genetics of speciation by cascade reinforcement

Species interactions drive diverse evolutionary outcomes. Speciation by cascade reinforcement represents one example of how species interactions can contribute to the proliferation of species. This process occurs when the divergence of mating traits in response to selection against interspecific hybridization incidentally leads to reproductive isolation among populations of the same species. Here, we investigated the population genetic outcomes of cascade reinforcement in North American chorus frogs (Hylidae: Pseudacris). Specifically, we estimated the frequency of hybridization among three taxa, assessed genetic structure within the focal species, P. feriarum, and ascertained the directionality of gene flow within P. feriarum across replicated contact zones via coalescent modeling. Through field observations and preliminary experimental crosses, we assessed whether hybridization is possible under natural and laboratory conditions. We found that hybridization occurs among P. feriarum and two conspecifics at a low rate in multiple contact zones, and that gene flow within the former species is unidirectional from allopatry into sympatry with these other species in three of four contact zones studied. We found evidence of substantial genetic structuring within P. feriarum including a divergent western allopatric cluster, a behaviorally-distinct sympatric South Carolina cluster, and several genetically-overlapping clusters from the remainder of the distribution. Furthermore, we found sub-structuring between reinforced and nonreinforced populations in the two most intensely-sampled contact zones. Our literature review indicated that P. feriarum hybridizes with at least five heterospecifics at the periphery of its range providing a mechanism for further intraspecific diversification. This work strengthens the evidence for cascade reinforcement in this clade, revealing the geographic and genetic landscape upon which this process can contribute to the proliferation of species.

Secondary contact rather than coexistence-Erebia butterflies in the Alps

Secondary contact zones are ideal systems to study the processes that govern the evolution of reproductive barriers, especially at advanced stages of the speciation process. An increase in reproductive isolation resulting from selection against maladaptive hybrids is thought to contribute to reproductive barrier buildup in secondary contact zones. Although such processes have been invoked for many systems, it remains unclear to which extent they influence contact zone dynamics in nature. Here, we study a very narrow contact zone between the butterfly species Erebia cassioides and Erebia tyndarus in the Swiss Alps. We quantified phenotypic traits related to wing shape and reproduction as well as ecology to compare the degree of intra- and interspecific differentiation. Even though only very few first-generation hybrids occur, we find no strong indications for current reinforcing selection, suggesting that if reinforcement occurred in our system, it likely operated in the past. Additionally, we show that both species differ less in their ecological niche at the contact zone than elsewhere, which could explain why coexistence between these butterflies may currently not be possible.

Neurogenomic divergence during speciation by reinforcement of mating behaviors in chorus frogs (Pseudacris)

BACKGROUND

Species interactions can promote mating behavior divergence, particularly when these interactions are costly due to maladaptive hybridization. Selection against hybridization can indirectly cause evolution of reproductive isolation within species, a process termed cascade reinforcement. This process can drive incipient speciation by generating divergent selection pressures among populations that interact with different species assemblages. Theoretical and empirical studies indicate that divergent selection on gene expression networks has the potential to increase reproductive isolation among populations. After identifying candidate synaptic transmission genes derived from neurophysiological studies in anurans, we test for divergence of gene expression in a system undergoing cascade reinforcement, the Upland Chorus Frog (Pseudacris feriarum).

RESULTS

Our analyses identified seven candidate synaptic transmission genes that have diverged between ancestral and reinforced populations of P. feriarum, including five that encode synaptic vesicle proteins. Our gene correlation network analyses revealed four genetic modules that have diverged between these populations, two possessing a significant concentration of neurotransmission enrichment terms: one for synaptic membrane components and the other for metabolism of the neurotransmitter nitric oxide. We also ascertained that a greater number of genes have diverged in expression by geography than by sex. Moreover, we found that more genes have diverged within females as compared to males between populations. Conversely, we observed no difference in the number of differentially-expressed genes within the ancestral compared to the reinforced population between the sexes.

CONCLUSIONS

This work is consistent with the idea that divergent selection on mating behaviors via cascade reinforcement contributed to evolution of gene expression in P. feriarum. Although our study design does not allow us to fully rule out the influence of environment and demography, the fact that more genes diverged in females than males points to a role for cascade reinforcement. Our discoveries of divergent candidate genes and gene networks related to neurotransmission support the idea that neural mechanisms of acoustic mating behaviors have diverged between populations, and agree with previous neurophysiological studies in frogs. Increasing support for this hypothesis, however, will require additional experiments under common garden conditions. Our work points to the importance of future replicated and tissue-specific studies to elucidate the relative contribution of gene expression divergence to the evolution of reproductive isolation during incipient speciation.

Reinforcement and the Proliferation of Species

Adaptive radiations are characterized by the rapid proliferation of species. Explaining how adaptive radiations occur therefore depends, in part, on identifying how populations become reproductively isolated-and ultimately become different species. Such reproductive isolation could arise when populations adapting to novel niches experience selection to avoid interbreeding and, consequently, evolve mating traits that minimize such hybridization via the process of reinforcement. Here, we highlight that a downstream consequence of reinforcement is divergence of conspecific populations, and this further divergence can instigate species proliferation. Moreover, we evaluate when reinforcement will-and will not-promote species proliferation. Finally, we discuss empirical approaches to test what role, if any, reinforcement plays in species proliferation and, consequently, in adaptive radiation. To date, reinforcement's downstream effects on species proliferation remain largely unknown and speculative. Because the ecological and evolutionary contexts in which adaptive radiations occur are conducive to reinforcement and its downstream consequences, adaptive radiations provide an ideal framework in which to evaluate reinforcement's role in diversification.

Female northern grass lizards judge mates by body shape to reinforce local adaptation

Background

Identifying the factors that contribute to divergence among populations in mate preferences is important for understanding of the manner in which premating reproductive isolation might arise and how this isolation may in turn contribute to the evolutionary process of population divergence. Here, we offered female northern grass lizards (Takydromus septentrionalis) a choice of males between their own population and another four populations to test whether the preferences that females display in the mating trials correlate with phenotypic adaptation to local environments, or to the neutral genetic distance measured by divergence of mitochondrial DNA sequence loci.

Results

Females showed a strong preference for native over foreign males. Females that mated with native versus foreign males did not differ from each other in mating latency, or copulation duration. From results of the structural equation modelling we knew that: 1) geographical distance directly contributed to genetic differentiation and environmental dissimilarity; 2) genetic differentiation and environmental dissimilarity indirectly contributed to female mate preference, largely through their effects on morphological divergence; and 3) females judged mates by body shape (appearance) and discriminated more strongly against morphologically less familiar allopatric males.

Conclusions

Local adaptation rather than neutral genetic distance influences female mate preference in T. septentrionalis. The tendency to avoid mating with foreign males may indicate that, in T. septentrionalis, local adaptations are more valuable than genetic novelties. Our results highlight the importance of comprehensive studies integrating ecological, molecular and behavioral approaches to understand population divergence in female mate preferences as the consequence of local adaptations.

To accept or reject heterospecific mates: behavioural decisions underlying premating isolation

Premating isolation in animals involves decision-making processes that affect whether individuals accept or reject heterospecific mates. An integrative understanding of the behavioural processes underlying heterospecific acceptance can clarify the conditions under which premating isolation evolves. As an illustration, we review how Reeve's (Reeve 1989 Am. Nat. 133, 407-435. (doi:10.1086/284926)) acceptance threshold model can help make sense of patterns of premating isolation in nature. This model derives a threshold trait value for acceptance for rejection of recipients of an action (e.g. mating) based on the fitness consequences of these decisions. We show that the maintenance of partial reproductive isolation can be an outcome of optimal acceptance thresholds, even in the face of reinforcement. We also use this model to clarify how the composition of multispecies communities can shape premating isolation. The acceptance threshold model can also be viewed as the behavioural underpinning of reproductive character displacement and cascading reinforcement. Finally, we highlight potential limitations of the acceptance threshold model with respect to investigating the role of sexual selection in speciation, and we propose that integration of behavioural models in speciation research will help us gain a full picture of the mechanisms underlying premating isolation. This article is part of the theme issue 'Signal detection theory in recognition systems: from evolving models to experimental tests'.

Ongoing hybridization obscures phylogenetic relationships in the Drosophila subquinaria species complex

Inferring evolutionary relationships among recently diverged lineages is necessary to understand how isolating barriers produce independent lineages. Here, we investigate the phylogenetic relationships between three incompletely isolated and closely related mushroom-feeding Drosophila species. These species form the Drosophila subquinaria species complex and consist of one Eurasian species (D. transversa) and two widespread North American species (D. subquinaria and D. recens) that are sympatric in central Canada. Although patterns of pre- and post-mating isolation among these species are well characterized, previous work on their phylogenetic relationships is limited and conflicting. In this study, we generated a multi-locus data set of 29 loci from across the genome sequenced in a population sample from each species, and then, we inferred species relationships and patterns of introgression. We find strong statistical support that D. subquinaria is paraphyletic, showing that samples from the geographic region sympatric with D. recens are most closely related to D. recens, whereas samples from the geographic region allopatric with D. recens are most closely related to D. transversa. We present several lines of evidence that both incomplete lineage sorting and gene flow are causing phylogenetic discordance. We suggest that ongoing gene flow primarily from D. recens into D. subquinaria in the sympatric part of their ranges causes phylogenetic uncertainty in the evolutionary history of these species. Our results highlight how population genetic data can be used to disentangle the sources of phylogenetic discordance among closely related species.

Strength of sexual and postmating prezygotic barriers varies between sympatric populations with different histories and species abundances

The impact of different reproductive barriers on species or population isolation may vary in different stages of speciation depending on evolutionary forces acting within species and through species' interactions. Genetic incompatibilities between interacting species are expected to reinforce prezygotic barriers in sympatric populations and lead to cascade reinforcement between conspecific populations living within and outside the areas of sympatry. We tested these predictions and studied whether and how the strength and target of reinforcement between Drosophila montana and Drosophila flavomontana vary between sympatric populations with different histories and species abundances. All barriers between D. montana females and D. flavomontana males were nearly complete, while in the reciprocal cross strong postzygotic isolation was accompanied by prezygotic barriers whose strength varied according to population composition. Sexual isolation between D. flavomontana females and D. montana males was increased in long-established sympatric populations, where D. flavomontana is abundant, while postmating prezygotic (PMPZ) barriers were stronger in populations where this species is a new invader and still rare and where female discrimination against heterospecific males was lower. Strengthening of sexual and PMPZ barriers in this cross also induced cascade reinforcement of respective barriers between D. flavomontana populations, which is a classic signature of reinforcement process.

A phylogenetic examination of host use evolution in the quinaria and testacea groups of Drosophila

Adaptive radiations provide an opportunity to examine complex evolutionary processes such as ecological specialization and speciation. While a well-resolved phylogenetic hypothesis is critical to completing such studies, the rapid rates of evolution in these groups can impede phylogenetic studies. Here we study the quinaria and testacea species groups of the immigrans-tripunctata radiation of Drosophila, which represent a recent adaptive radiation and are a developing model system for ecological genetics. We were especially interested in understanding host use evolution in these species. In order to infer a phylogenetic hypothesis for this group we sampled loci from both the nuclear genome and the mitochondrial DNA to develop a dataset of 43 protein-coding loci for these two groups along with their close relatives in the immigrans-tripunctata radiation. We used this dataset to examine their evolutionary relationships along with the evolution of feeding behavior. Our analysis recovers strong support for the monophyly of the testacea but not the quinaria group. Results from our ancestral state reconstruction analysis suggests that the ancestor of the testacea and quinaria groups exhibited mushroom-feeding. Within the quinaria group, we infer that transition to vegetative feeding occurred twice, and that this transition did not coincide with a genome-wide change in the rate of protein evolution.

Fine-scale geographic patterns of gene flow and reproductive character displacement in Drosophila subquinaria and Drosophila recens

When two species are incompletely isolated, strengthening premating isolation barriers in response to the production of low fitness hybrids may complete the speciation process. Here, we use the sister species Drosophila subquinaria and Drosophila recens to study the conditions under which this reinforcement of species boundaries occurs in natural populations. We first extend the region of known sympatry between these species, and then we conduct a fine-scale geographic survey of mate discrimination coupled with estimates of gene flow within and admixture between species. Within D. subquinaria, reinforcement is extremely effective: we find variation in mate discrimination both against D. recens males and against conspecific allopatric males on the scale of a few kilometres and in the face of gene flow both from conspecific populations and introgression from D. recens. In D. recens, we do not find evidence for increased mate discrimination in sympatry, even where D. recens is rare, consistent with substantial gene flow throughout the species' range. Finally, we find that introgression between species is asymmetric, with more from D. recens into D. subquinaria than vice versa. Within each species, admixture is highest in the geographic region where it is rare relative to the other species, suggesting that when hybrids are produced they are of low fitness. In sum, reinforcement within D. subquinaria is effective at maintaining species boundaries, but even when reinforcing selection is strong it may not always result in a pattern of strong reproductive character displacement due to variation in the frequency of hybridization and gene flow from neighbouring populations.

Male-driven reproductive and agonistic character displacement in darters and its implications for speciation in allopatry

Selection against hybridization can cause mating traits to diverge between species in sympatry via reproductive character displacement (RCD). Additionally, selection against interspecific fighting can cause aggressive traits to diverge between sympatric species via agonistic character displacement (ACD). By directly affecting conspecific recognition traits, RCD and ACD between species can also incidentally cause divergence in mating and fighting traits among populations within a species [termed cascade RCD (CRCD) and cascade ACD]. Here, we demonstrate patterns consistent with male-driven RCD and ACD in 2 groups of darters (orangethroat darter clade Ceasia and rainbow darter Etheostoma caeruleum). In both groups, males that occur in sympatry (between Ceasia and E. caeruleum) have higher levels of preference for mating and fighting with conspecifics over heterospecifics than do males from allopatry. This is consistent with RCD and ACD. We also found patterns consistent with CRCD and cascade ACD among species of Ceasia. Ceasia males that are sympatric to E. caeruleum (but allopatric to one another) also have heightened preferences for mating and fighting with conspecific versus heterospecific Ceasia. In contrast, Ceasia males that are allopatric to E. caeruleum readily mate and fight with heterospecific Ceasia. We suggest that RCD and ACD between Ceasia and E. caeruleum has incidentally led to divergence in mating and fighting traits among Ceasia species. This study is unique in that male preferences evolve via both RCD (male preference for conspecific females) and ACD (male preference to fight conspecific males) which leads to subsequent divergence among allopatric lineages.

Correlated evolution of male and female reproductive traits drive a cascading effect of reinforcement in Drosophila yakuba

Selection against maladaptive hybridization can drive the evolution of reproductive isolation in a process called reinforcement. While the importance of reinforcement in evolution has been historically debated, many examples now exist. Despite these examples, we typically lack a detailed understanding of the mechanisms limiting the spread of reinforced phenotypes throughout a species' range. Here we address this issue in the fruit fly Drosophila yakuba, a species that hybridizes with its sister species D. santomea and is undergoing reinforcement in a well-defined hybrid zone on the island of São Tomé. Within this region, female D. yakuba show increased postmating-prezygotic (gametic) isolation towards D. santomea when compared with females from allopatric populations. We use a combination of natural collections, fertility assays, and experimental evolution to understand why reinforced gametic isolation in D. yakuba is confined to this hybrid zone. We show that, among other traits, D. yakuba males from sympatric populations sire fewer progeny than allopatric males when mated to allopatric D. yakuba females. Our results provide a novel example of reinforcement acting on a postmating-prezygotic trait in males, resulting in a cascade of reproductive isolation among conspecific populations.

Geographic variation in hybridization across a reinforcement contact zone of chorus frogs (Pseudacris)

Reinforcement contact zones, which are secondary contact zones where species are diverging in reproductive behaviors due to selection against hybridization, represent natural laboratories for studying speciation‐in‐action. Here, we examined replicate localities across the entire reinforcement contact zone between North American chorus frogs Pseudacris feriarum and P. nigrita to investigate geographic variation in hybridization frequencies and to assess whether reinforcement may have contributed to increased genetic divergence within species. Previous work indicated these species have undergone reproductive character displacement (RCD) in male acoustic signals and female preferences due to reinforcement. We also examined acoustic signal variation across the contact zone to assess whether signal characteristics reliably predict hybrid index and to elucidate whether the degree of RCD predicts hybridization rate. Using microsatellites, mitochondrial sequences, and acoustic signal information from >1,000 individuals across >50 localities and ten sympatric focal regions, we demonstrate: (1) hybridization occurs and (2) varies substantially across the geographic range of the contact zone, (3) hybridization is asymmetric and in the direction predicted from observed patterns of asymmetric RCD, (4) in one species, genetic distance is higher between conspecific localities where one or both have been reinforced than between nonreinforced localities, after controlling for geographic distance, (5) acoustic signal characters strongly predict hybrid index, and (6) the degree of RCD does not strongly predict admixture levels. By showing that hybridization occurs in all sympatric localities, this study provides the fifth and final line of evidence that reproductive character displacement is due to reinforcement in the chorus frog contact zone. Furthermore, this work suggests that the dual action of cascade reinforcement and partial geographic isolation is promoting genetic diversification within one of the reinforced species.

Poison frog color morphs express assortative mate preferences in allopatry but not sympatry

The concurrent divergence of mating traits and preferences is necessary for the evolution of reproductive isolation via sexual selection, and such coevolution has been demonstrated in diverse lineages. However, the extent to which assortative mate preferences are sufficient to drive reproductive isolation in nature is less clear. Natural contact zones between lineages divergent in traits and preferences provide exceptional opportunities for testing the predicted evolutionary consequences of such divergence. The strawberry poison frog (Oophaga pumilio) displays extreme color polymorphism in and around the young Bocas del Toro archipelago. In a transition zone between red and blue allopatric lineages, we asked whether female preferences diverged along with coloration, and whether any divergent preferences persist in a zone of sympatry. When choosing among red, blue and phenotypically intermediate males, females from monomorphic red and monomorphic blue populations both expressed assortative preferences. However, red, blue, and intermediate females from the contact zone all preferred red males, suggesting that divergent preferences may be insufficient to effect behavioral isolation. Our results highlight the complexity of behavioral isolation, and the need for studies that can reveal the circumstances under which divergent preferences do and do not contribute to speciation.

Population differentiation at a regional scale in spadefoot toads: contributions of distance and divergent selective environments

The causes of population differentiation can provide insight into the origins of early barriers to gene flow. Two key drivers of population differentiation are geographic distance and local adaptation to divergent selective environments. When reproductive isolation arises because some populations of a species are under selection to avoid hybridization while others are not, population differentiation and even speciation can result. Spadefoot toad populations Spea multiplicata that are sympatric with a congener have undergone reinforcement. This reinforcement has resulted not only in increased reproductive isolation from the congener, but also in the evolution of reproductive isolation from nearby and distant conspecific allopatric populations. We used multiple approaches to evaluate the contributions of geographic distance and divergent selective environments to population structure across this regional scale in S. multiplicata, based on genotypes from six nuclear microsatellite markers. We compared groups of populations varying in both geographic location and in the presence of a congener. Hierarchical F-statistics and results from cluster analyses and discriminant analyses of principal components all indicate that geographic distance is the stronger contributor to genetic differentiation among S. multiplicata populations at a regional scale. However, we found evidence that adaptation to divergent selective environments also contributes to population structure. Our findings highlight how variation in the balance of evolutionary forces acting across a species' range can lead to variation in the relative contributions of geographic distance and local adaptation to population differentiation across different spatial scales.

Reinforcement's incidental effects on reproductive isolation between conspecifics

Reinforcement—the process whereby maladaptive hybridization leads to the strengthening of prezygotic isolation between species—has a long history in the study of speciation. Because reinforcement affects traits involved in mate choice and fertility, it can have indirect effects on reproductive isolation between populations within species. Here we review examples of these “cascading effects of reinforcement” (CER) and discuss different mechanisms through which they can arise. We discuss three factors that are predicted to influence the potential occurrence of CER: rates of gene flow among populations, the strength of selection acting on the traits involved in reinforcement, and the genetic basis of those traits. We suggest that CER is likely if (1) the rate of gene flow between conspecific populations is low; (2) divergent selection acts on phenotypes involved in reinforcement between sympatric and allopatric populations; and (3) the genetic response to reinforcement differs among conspecific populations subject to parallel reinforcing selection. Future work continuing to address gene flow, selection, and the genetic basis of the traits involved in the reinforcement will help develop a better understanding of reinforcement as a process driving the production of species diversity, both directly and incidentally.

Is cascade reinforcement likely when sympatric and allopatric populations exchange migrants?

When partially reproductively isolated species come back into secondary contact, these taxa may diverge in mating preferences and sexual cues to avoid maladaptive hybridization, a process known as reinforcement. This phenomenon often leads to reproductive character displacement (RCD) between sympatric and allopatric populations of reinforcing species that differ in their exposure to hybridization. Recent discussions have reinvigorated the idea that RCD may give rise to additional speciation between conspecific sympatric and allopatric populations, dubbing the concept "cascade reinforcement." Despite some empirical studies supporting cascade reinforcement, we still know very little about the conditions for its evolution. In the present article, we address this question by developing an individual-based population genetic model that explicitly simulates cascade reinforcement when one of the hybridizing species is split into sympatric and allopatric populations. Our results show that when sympatric and allopatric populations reside in the same environment and only differ in their exposure to maladaptive hybridization, migration between them generally inhibits the evolution of cascade by spreading the reinforcement alleles from sympatry into allopatry and erasing RCD. Under these conditions, cascade reinforcement only evolved when migration rate between sympatric and allopatric populations was very low. This indicates that stabilizing sexual selection in allopatry is generally ineffective in preventing the spread of reinforcement alleles. Only when sympatric and allopatric populations experienced divergent ecological selection did cascade reinforcement evolve in the presence of substantial migration. These predictions clarify the conditions for cascade reinforcement and facilitate our understanding of existing cases in nature.

Patterns of reproductive isolation in the Drosophila subquinaria complex: can reinforced premating isolation cascade to other species?

The reinforcement of premating barriers due to reduced hybrid fitness in sympatry may cause secondary sexual isolation within a species as a by-product. Consistent with this, in the fly Drosophila subquinaria, females that are sympatric with D. recens mate at very low rates not only with D. recens, but also with conspecific D. subquinaria males from allopatry. Here, we ask if these effects of reinforcement cascade more broadly to affect sexual isolation with other closely related species. We assay reproductive isolation of these species with D. transversa and find that choosy D. subquinaria females from the region sympatric with D. recens discriminate strongly against male D. transversa, whereas D. subquinaria from the allopatric region do not. This increased sexual isolation cannot be explained by natural selection to avoid mating with this species, as they are allopatric in geographic range and we do not identify any intrinsic postzygotic isolation between D. subquinaria and D. transversa. Variation in epicuticular hydrocarbons, which are used as mating signals in D. subquinaria, follow patterns of premating isolation: D. transversa and allopatric D. subquinaria are most similar to each other and differ from sympatric D. subquinaria, and those of D. recens are distinct from the other two species. We suggest that the secondary effects of reinforcement may cascade to strengthen reproductive isolation with other species that were not a target of selection. These effects may enhance the divergence that occurs in allopatry to help explain why some species are already sexually isolated upon secondary contact.