@article{12514, abstract = {The concept of a “speciation continuum” has gained popularity in recent decades. It emphasizes speciation as a continuous process that may be studied by comparing contemporary population pairs that show differing levels of divergence. In their recent perspective article in Evolution, Stankowski and Ravinet provided a valuable service by formally defining the speciation continuum as a continuum of reproductive isolation, based on opinions gathered from a survey of speciation researchers. While we agree that the speciation continuum has been a useful concept to advance the understanding of the speciation process, some intrinsic limitations exist. Here, we advocate for a multivariate extension, the speciation hypercube, first proposed by Dieckmann et al. in 2004, but rarely used since. We extend the idea of the speciation cube and suggest it has strong conceptual and practical advantages over a one-dimensional model. We illustrate how the speciation hypercube can be used to visualize and compare different speciation trajectories, providing new insights into the processes and mechanisms of speciation. A key strength of the speciation hypercube is that it provides a unifying framework for speciation research, as it allows questions from apparently disparate subfields to be addressed in a single conceptual model.}, author = {Bolnick, Daniel I. and Hund, Amanda K. and Nosil, Patrik and Peng, Foen and Ravinet, Mark and Stankowski, Sean and Subramanian, Swapna and Wolf, Jochen B.W. and Yukilevich, Roman}, issn = {1558-5646}, journal = {Evolution: International journal of organic evolution}, number = {1}, pages = {318--328}, publisher = {Oxford University Press}, title = {{A multivariate view of the speciation continuum}}, doi = {10.1093/evolut/qpac004}, volume = {77}, year = {2023}, } @article{14604, abstract = {Sex chromosomes have evolved independently multiple times, but why some are conserved for more than 100 million years whereas others turnover rapidly remains an open question. Here, we examine the homology of sex chromosomes across nine orders of insects, plus the outgroup springtails. We find that the X chromosome is likely homologous across insects and springtails; the only exception is in the Lepidoptera, which has lost the X and now has a ZZ/ZW sex-chromosome system. These results suggest the ancestral insect X chromosome has persisted for more than 450 million years—the oldest known sex chromosome to date. Further, we propose that the shrinking of gene content the dipteran X chromosome has allowed for a burst of sex-chromosome turnover that is absent from other speciose insect orders.}, author = {Toups, Melissa A and Vicoso, Beatriz}, issn = {1558-5646}, journal = {Evolution}, number = {11}, pages = {2504--2511}, publisher = {Oxford University Press}, title = {{The X chromosome of insects likely predates the origin of class Insecta}}, doi = {10.1093/evolut/qpad169}, volume = {77}, year = {2023}, } @article{11334, abstract = {Hybridization is a common evolutionary process with multiple possible outcomes. In vertebrates, interspecific hybridization has repeatedly generated parthenogenetic hybrid species. However, it is unknown whether the generation of parthenogenetic hybrids is a rare outcome of frequent hybridization between sexual species within a genus or the typical outcome of rare hybridization events. Darevskia is a genus of rock lizards with both hybrid parthenogenetic and sexual species. Using capture sequencing, we estimate phylogenetic relationships and gene flow among the sexual species, to determine how introgressive hybridization relates to the origins of parthenogenetic hybrids. We find evidence for widespread hybridization with gene flow, both between recently diverged species and deep branches. Surprisingly, we find no signal of gene flow between parental species of the parthenogenetic hybrids, suggesting that the parental pairs were either reproductively or geographically isolated early in their divergence. The generation of parthenogenetic hybrids in Darevskia is, then, a rare outcome of the total occurrence of hybridization within the genus, but the typical outcome when specific species pairs hybridize. Our results question the conventional view that parthenogenetic lineages are generated by hybridization in a window of divergence. Instead, they suggest that some lineages possess specific properties that underpin successful parthenogenetic reproduction.}, author = {Freitas, Susana and Westram, Anja M and Schwander, Tanja and Arakelyan, Marine and Ilgaz, Çetin and Kumlutas, Yusuf and Harris, David James and Carretero, Miguel A. and Butlin, Roger K.}, issn = {1558-5646}, journal = {Evolution}, number = {5}, pages = {899--914}, publisher = {Wiley}, title = {{Parthenogenesis in Darevskia lizards: A rare outcome of common hybridization, not a common outcome of rare hybridization}}, doi = {10.1111/evo.14462}, volume = {76}, year = {2022}, } @article{12234, abstract = {Hybrid speciation—the origin of new species resulting from the hybridization of genetically divergent lineages—was once considered rare, but genomic data suggest that it may occur more often than once thought. In this study, Noguerales and Ortego found genomic evidence supporting the hybrid origin of a grasshopper that is able to exploit a broader range of host plants than either of its putative parents.}, author = {Stankowski, Sean}, issn = {1558-5646}, journal = {Evolution}, keywords = {General Agricultural and Biological Sciences, Genetics, Ecology, Evolution, Behavior and Systematics}, number = {11}, pages = {2784--2785}, publisher = {Wiley}, title = {{Digest: On the origin of a possible hybrid species}}, doi = {10.1111/evo.14632}, volume = {76}, year = {2022}, } @article{12247, abstract = {Chromosomal inversions have been shown to play a major role in a local adaptation by suppressing recombination between alternative arrangements and maintaining beneficial allele combinations. However, so far, their importance relative to the remaining genome remains largely unknown. Understanding the genetic architecture of adaptation requires better estimates of how loci of different effect sizes contribute to phenotypic variation. Here, we used three Swedish islands where the marine snail Littorina saxatilis has repeatedly evolved into two distinct ecotypes along a habitat transition. We estimated the contribution of inversion polymorphisms to phenotypic divergence while controlling for polygenic effects in the remaining genome using a quantitative genetics framework. We confirmed the importance of inversions but showed that contributions of loci outside inversions are of similar magnitude, with variable proportions dependent on the trait and the population. Some inversions showed consistent effects across all sites, whereas others exhibited site-specific effects, indicating that the genomic basis for replicated phenotypic divergence is only partly shared. The contributions of sexual dimorphism as well as environmental factors to phenotypic variation were significant but minor compared to inversions and polygenic background. Overall, this integrated approach provides insight into the multiple mechanisms contributing to parallel phenotypic divergence.}, author = {Koch, Eva L. and Ravinet, Mark and Westram, Anja M and Johannesson, Kerstin and Butlin, Roger K.}, issn = {1558-5646}, journal = {Evolution}, keywords = {General Agricultural and Biological Sciences, Genetics, Ecology, Evolution, Behavior and Systematics}, number = {10}, pages = {2332--2346}, publisher = {Wiley}, title = {{Genetic architecture of repeated phenotypic divergence in Littorina saxatilis evolution}}, doi = {10.1111/evo.14602}, volume = {76}, year = {2022}, } @article{8743, abstract = {Montane cloud forests are areas of high endemism, and are one of the more vulnerable terrestrial ecosystems to climate change. Thus, understanding how they both contribute to the generation of biodiversity, and will respond to ongoing climate change, are important and related challenges. The widely accepted model for montane cloud forest dynamics involves upslope forcing of their range limits with global climate warming. However, limited climate data provides some support for an alternative model, where range limits are forced downslope with climate warming. Testing between these two models is challenging, due to the inherent limitations of climate and pollen records. We overcome this with an alternative source of historical information, testing between competing model predictions using genomic data and demographic analyses for a species of beetle tightly associated to an oceanic island cloud forest. Results unequivocally support the alternative model: populations that were isolated at higher elevation peaks during the Last Glacial Maximum are now in contact and hybridizing at lower elevations. Our results suggest that genomic data are a rich source of information to further understand how montane cloud forest biodiversity originates, and how it is likely to be impacted by ongoing climate change.}, author = {Salces-Castellano, Antonia and Stankowski, Sean and Arribas, Paula and Patino, Jairo and Karger, Dirk N. and Butlin, Roger and Emerson, Brent C.}, issn = {1558-5646}, journal = {Evolution}, number = {2}, pages = {231--244}, publisher = {Wiley}, title = {{Long-term cloud forest response to climate warming revealed by insect speciation history}}, doi = {10.1111/evo.14111}, volume = {75}, year = {2021}, } @article{9252, abstract = {This paper analyses the conditions for local adaptation in a metapopulation with infinitely many islands under a model of hard selection, where population size depends on local fitness. Each island belongs to one of two distinct ecological niches or habitats. Fitness is influenced by an additive trait which is under habitat‐dependent directional selection. Our analysis is based on the diffusion approximation and accounts for both genetic drift and demographic stochasticity. By neglecting linkage disequilibria, it yields the joint distribution of allele frequencies and population size on each island. We find that under hard selection, the conditions for local adaptation in a rare habitat are more restrictive for more polygenic traits: even moderate migration load per locus at very many loci is sufficient for population sizes to decline. This further reduces the efficacy of selection at individual loci due to increased drift and because smaller populations are more prone to swamping due to migration, causing a positive feedback between increasing maladaptation and declining population sizes. Our analysis also highlights the importance of demographic stochasticity, which exacerbates the decline in numbers of maladapted populations, leading to population collapse in the rare habitat at significantly lower migration than predicted by deterministic arguments.}, author = {Szep, Eniko and Sachdeva, Himani and Barton, Nicholas H}, issn = {1558-5646}, journal = {Evolution}, keywords = {Genetics, Ecology, Evolution, Behavior and Systematics, General Agricultural and Biological Sciences}, number = {5}, pages = {1030--1045}, publisher = {Wiley}, title = {{Polygenic local adaptation in metapopulations: A stochastic eco‐evolutionary model}}, doi = {10.1111/evo.14210}, volume = {75}, year = {2021}, } @article{9374, abstract = {If there are no constraints on the process of speciation, then the number of species might be expected to match the number of available niches and this number might be indefinitely large. One possible constraint is the opportunity for allopatric divergence. In 1981, Felsenstein used a simple and elegant model to ask if there might also be genetic constraints. He showed that progress towards speciation could be described by the build‐up of linkage disequilibrium among divergently selected loci and between these loci and those contributing to other forms of reproductive isolation. Therefore, speciation is opposed by recombination, because it tends to break down linkage disequilibria. Felsenstein then introduced a crucial distinction between “two‐allele” models, which are subject to this effect, and “one‐allele” models, which are free from the recombination constraint. These fundamentally important insights have been the foundation for both empirical and theoretical studies of speciation ever since.}, author = {Butlin, Roger K. and Servedio, Maria R. and Smadja, Carole M. and Bank, Claudia and Barton, Nicholas H and Flaxman, Samuel M. and Giraud, Tatiana and Hopkins, Robin and Larson, Erica L. and Maan, Martine E. and Meier, Joana and Merrill, Richard and Noor, Mohamed A. F. and Ortiz‐Barrientos, Daniel and Qvarnström, Anna}, issn = {1558-5646}, journal = {Evolution}, keywords = {Genetics, Ecology, Evolution, Behavior and Systematics, General Agricultural and Biological Sciences}, number = {5}, pages = {978--988}, publisher = {Wiley}, title = {{Homage to Felsenstein 1981, or why are there so few/many species?}}, doi = {10.1111/evo.14235}, volume = {75}, year = {2021}, } @article{9383, abstract = {A primary roadblock to our understanding of speciation is that it usually occurs over a timeframe that is too long to study from start to finish. The idea of a speciation continuum provides something of a solution to this problem; rather than observing the entire process, we can simply reconstruct it from the multitude of speciation events that surround us. But what do we really mean when we talk about the speciation continuum, and can it really help us understand speciation? We explored these questions using a literature review and online survey of speciation researchers. Although most researchers were familiar with the concept and thought it was useful, our survey revealed extensive disagreement about what the speciation continuum actually tells us. This is due partly to the lack of a clear definition. Here, we provide an explicit definition that is compatible with the Biological Species Concept. That is, the speciation continuum is a continuum of reproductive isolation. After outlining the logic of the definition in light of alternatives, we explain why attempts to reconstruct the speciation process from present‐day populations will ultimately fail. We then outline how we think the speciation continuum concept can continue to act as a foundation for understanding the continuum of reproductive isolation that surrounds us.}, author = {Stankowski, Sean and Ravinet, Mark}, issn = {1558-5646}, journal = {Evolution}, number = {6}, pages = {1256--1273}, publisher = {Oxford University Press}, title = {{Defining the speciation continuum}}, doi = {10.1111/evo.14215}, volume = {75}, year = {2021}, } @article{6680, abstract = {This paper analyzes how partial selfing in a large source population influences its ability to colonize a new habitat via the introduction of a few founder individuals. Founders experience inbreeding depression due to partially recessive deleterious alleles as well as maladaptation to the new environment due to selection on a large number of additive loci. I first introduce a simplified version of the Inbreeding History Model (Kelly, 2007) in order to characterize mutation‐selection balance in a large, partially selfing source population under selection involving multiple non‐identical loci. I then use individual‐based simulations to study the eco‐evolutionary dynamics of founders establishing in the new habitat under a model of hard selection. The study explores how selfing rate shapes establishment probabilities of founders via effects on both inbreeding depression and adaptability to the new environment, and also distinguishes the effects of selfing on the initial fitness of founders from its effects on the long‐term adaptive response of the populations they found. A high rate of (but not complete) selfing is found to aid establishment over a wide range of parameters, even in the absence of mate limitation. The sensitivity of the results to assumptions about the nature of polygenic selection are discussed.}, author = {Sachdeva, Himani}, issn = {1558-5646}, journal = {Evolution}, number = {9}, pages = {1729--1745}, publisher = {Wiley}, title = {{Effect of partial selfing and polygenic selection on establishment in a new habitat}}, doi = {10.1111/evo.13812}, volume = {73}, year = {2019}, } @article{1851, abstract = {We consider mating strategies for females who search for males sequentially during a season of limited length. We show that the best strategy rejects a given male type if encountered before a time-threshold but accepts him after. For frequency-independent benefits, we obtain the optimal time-thresholds explicitly for both discrete and continuous distributions of males, and allow for mistakes being made in assessing the correct male type. When the benefits are indirect (genes for the offspring) and the population is under frequency-dependent ecological selection, the benefits depend on the mating strategy of other females as well. This case is particularly relevant to speciation models that seek to explore the stability of reproductive isolation by assortative mating under frequency-dependent ecological selection. We show that the indirect benefits are to be quantified by the reproductive values of couples, and describe how the evolutionarily stable time-thresholds can be found. We conclude with an example based on the Levene model, in which we analyze the evolutionarily stable assortative mating strategies and the strength of reproductive isolation provided by them.}, author = {Priklopil, Tadeas and Kisdi, Eva and Gyllenberg, Mats}, issn = {1558-5646}, journal = {Evolution}, number = {4}, pages = {1015 -- 1026}, publisher = {Wiley}, title = {{Evolutionarily stable mating decisions for sequentially searching females and the stability of reproductive isolation by assortative mating}}, doi = {10.1111/evo.12618}, volume = {69}, year = {2015}, } @article{9931, abstract = {Gene duplication is important in evolution, because it provides new raw material for evolutionary adaptations. Several existing hypotheses about the causes of duplicate retention and diversification differ in their emphasis on gene dosage, subfunctionalization, and neofunctionalization. Little experimental data exist on the relative importance of gene expression changes and changes in coding regions for the evolution of duplicate genes. Furthermore, we do not know how strongly the environment could affect this importance. To address these questions, we performed evolution experiments with the TEM-1 beta lactamase gene in Escherichia coli to study the initial stages of duplicate gene evolution in the laboratory. We mimicked tandem duplication by inserting two copies of the TEM-1 gene on the same plasmid. We then subjected these copies to repeated cycles of mutagenesis and selection in various environments that contained antibiotics in different combinations and concentrations. Our experiments showed that gene dosage is the most important factor in the initial stages of duplicate gene evolution, and overshadows the importance of point mutations in the coding region.}, author = {Dhar, Riddhiman and Bergmiller, Tobias and Wagner, Andreas}, issn = {1558-5646}, journal = {Evolution}, number = {6}, pages = {1775--1791}, publisher = {Wiley}, title = {{Increased gene dosage plays a predominant role in the initial stages of evolution of duplicate TEM-1 beta lactamase genes}}, doi = {10.1111/evo.12373}, volume = {68}, year = {2014}, } @article{4305, abstract = {The common shrew (Sorex araneus) is subdivided into several karyotypic races in Britain. Two of these races meet near Oxford to form the "Oxford-Hermitage" hybrid zone. We present a model which describes this system as a "tension zone," i.e., a set of clines maintained by a balance between dispersal and selection against chromosomal heterozygotes. The Oxford and Hermitage races differ by Robertsonian fusions with monobrachial homology (kq, no versus ko), and so F1 hybrids between them would have low fertility. However, the acrocentric karyotype is found at high frequency within the hybrid zone, so that complex Robertsonian heterozygotes (kq no/q ko n) are replaced by more fertile combinations, such as (kq no/k q n o). This suggests that the hybrid zone has been modified so as to increase hybrid fitness. Mathematical analysis and simulation show that, if selection against complex heterozygotes is sufficiently strong relative to selection against simple heterozygotes, acrocentrics increase, and displace the clines for kq and no from the cline for ko. Superimposed on this separation is a tendency for the hybrid zone to move m favor of the Oxford (kq no) race. We compare the model with estimates of linkage disequilibrium and cline shape made from field data.}, author = {Hatfield, Todd and Barton, Nicholas H and Searle, Jeremy}, issn = {1558-5646}, journal = {Evolution; International Journal of Organic Evolution}, number = {4}, pages = {1129 -- 1145}, publisher = {Wiley-Blackwell}, title = {{A model of a hybrid zone between two chromosomal races of the common shrew (Sorex araneus)}}, doi = {10.1111/j.1558-5646.1992.tb00624.x}, volume = {46}, year = {1992}, } @article{4308, author = {Barton, Nicholas H}, issn = {1558-5646}, journal = {Evolution; International Journal of Organic Evolution}, number = {2}, pages = {551 -- 557}, publisher = {Wiley-Blackwell}, title = {{On the spread of new gene combinations in the third phase of Wright's shifting balance}}, volume = {46}, year = {1992}, } @article{3646, abstract = {We compare the pattern of morphological and electrophoretic variation in the hybrid zone between Bombina bombina and B. variegata across two transects: one near Cracow and one 200 km away, near Przemysl in southeastern Poland. Morphological variation across the Przemysl transect had been surveyed more than 50 years ago; though we found a significant shift at one site, there is no evidence for gross movement over this period. Morphological and electrophoretic changes coincide, and the average shape of the clines is the same across both transects. At the center, most of the change in frequency of six diagnostic allozymes occurs within w = 6.05 km (2-unit support limits 5.56-6.54 km). These steep gradients are generated not by selection on the allozymes themselves, but by associations with other loci: though these markers are unlinked, they are in strong linkage disequilibrium with each other [R = D/ = 0.22 (0.15-0.29) at the center]. Disequilibria are broken up as alleles diffuse away from the zone and flow into the new genetic background. The net barrier to the flow of genes from bombina into variegata, which is generated by these disequilibria, is B = 51 (22-81) km. The fitness of hybrids must be substantially reduced to produce such a barrier [W̄H/W̄P = 0.58 (0.54-0.68)], and this selection must be spread over many loci [N = 55 (26-88)]. Alleles introgress significantly less far than would be expected from the age of the zone and the estimated dispersal rate [σ = 0.99 (0.82-1.14) km gen.-1/2]: this implies selection of se = 0.37 (0.15-0.58)% on the enzymes themselves. There is weak but significant linkage disequilibrium well away from the center of the zone; this, together with the presence of parental and F1 genotypes, suggests some long-range migration. However, such migration is not likely to cause significant introgression. }, author = {Szymura, Jacek and Barton, Nicholas H}, issn = {1558-5646}, journal = {Evolution}, number = {2}, pages = {237 -- 261}, publisher = {Wiley-Blackwell}, title = {{The genetic structure of the hybrid zone between the fire-bellied toads Bombina bombina and B. variegata: comparisons between transects and between loci}}, doi = {10.1111/j.1558-5646.1991.tb04400.x}, volume = {45}, year = {1991}, } @article{3648, abstract = {We investigate the probability of fixation of a chromosome rearrangement in a subdivided population, concentrating on the limit where migration is so large relative to selection (m ≫ s) that the population can be thought of as being continuously distributed. We study two demes, and one- and two-dimensional populations. For two demes, the probability of fixation in the limit of high migration approximates that of a population with twice the size of a single deme: migration therefore greatly reduces the fixation probability. However, this behavior does not extend to a large array of demes. Then, the fixation probability depends primarily on neighborhood size (Nb), and may be appreciable even with strong selection and free gene flow (≈exp(-B·Nb) in one dimension, ≈exp(-B\cdotNb) in two dimensions). Our results are close to those for the more tractable case of a polygenic character under disruptive selection.}, author = {Barton, Nicholas H and Rouhani, Shahin}, issn = {1558-5646}, journal = {Evolution}, number = {3}, pages = {499 -- 517}, publisher = {Wiley-Blackwell}, title = {{The probability of fixation of a new karyotype in a continuous population}}, doi = {10.1111/j.1558-5646.1991.tb04326.x}, volume = {45}, year = {1991}, } @article{3653, abstract = {Frequency-dependent selection on warning color can maintain narrow hybrid zones between unpalatable prey taxa. To measure such selection, we transferred marked Heliconius erato (Lepidoptera: Nymphalidae) in both directions across a 10-km-wide hybrid zone between Peruvian races differing in color pattern. These experimental H. erato were released at four sites, along with control H. erato of the phenotype native to each site. Survival of experimental butterflies was significantly lower than that of controls at two sites and overall. Most selection, measured as differences in survival, occurred soon after release. Selection against foreign morphs was 52% (confidence limits: 25-71%) and was probably due to bird attacks on unusual warning-color morphs (more than 10% of the recaptures had beak marks). Since only three major loci determine the color-pattern differences, this suggests an average selection coefficient of 0.17 per locus, sufficient to maintain the narrow clines in H. erato.}, author = {Mallet, James and Barton, Nicholas H}, issn = {1558-5646}, journal = {Evolution}, number = {2}, pages = {421 -- 431}, publisher = {Wiley-Blackwell}, title = {{Strong natural selection in a warning color hybrid zone}}, doi = {10.2307/2409217 }, volume = {43}, year = {1989}, } @article{4309, abstract = {Three methods for estimating the average level of gene flow in natural population are discussed and compared. The three methods are FST, rare alleles, and maximum likelihood. All three methods yield estimates of the combination of parameters (the number of migrants [Nm] in a demic model or the neighborhood size [4πDσ2] in a continuum model) that determines the relative importance of gene flow and genetic drift. We review the theory underlying these methods and derive new analytic results for the expectation of FST in stepping-stone and continuum models when small sets of samples are taken. We also compare the effectiveness of the different methods using a variety of simulated data. We found that the FST and rare-alleles methods yield comparable estimates under a wide variety of conditions when the population being sampled is demographically stable. They are roughly equally sensitive to selection and to variation in population structure, and they approach their equilibrium values at approximately the same rate. We found that two different maximum-likelihood methods tend to yield biased estimates when relatively small numbers of locations are sampled but more accurate estimates when larger numbers are sampled. Our conclusion is that, although FST and rare-alleles methods are expected to be equally effective in analyzing ideal data, practical problems in estimating the frequencies of rare alleles in electrophoretic studies suggest that FST is likely to be more useful under realistic conditions.}, author = {Slatkin, Montgomery and Barton, Nicholas H}, issn = {1558-5646}, journal = {Evolution; International Journal of Organic Evolution}, number = {7}, pages = {1349 -- 1368}, publisher = {Wiley-Blackwell}, title = {{A comparison of three methods for estimating average levels of gene flow}}, doi = {10.1111/j.1558-5646.1989.tb02587.x }, volume = {43}, year = {1989}, } @article{4321, abstract = {The fire-bellied toads Bombina bombina and B. variegata differ extensively in biochemistry, morphology, and behavior. We use a survey of five diagnostic enzyme loci across the hybrid zone near Cracow in Southern Poland to estimate the dispersal rate, selection pressures, and numbers of loci which maintain this zone. The enzyme clines coincide closely with each other and with morphological and mitochondrial DNA clines. Although the zone lies on a broad transition between environments suitable for bombina and variegata, the close concordance of diverse characters, together with increased aberrations and mortality in hybrids, suggest that the zone is maintained largely by selection against hybrids. There are strong “linkage disequilibria” between each pair of (unlinked) enzyme loci (R̄ = 0.129 [2-unit support limits: 0.119–0.139]). These are probably caused by gene flow into the zone, and they give an estimate of dispersal (σ = 890 [790–940] m gen−½). The clines are sharply stepped, with most of the change occurring within 6.15 (5.45–6.45) km, but with long tails of introgression on either side. This implies that the effective selection pressure on each enzyme marker (due largely to disequilibrium with other loci) is s* = 0.17 (0.159–0.181) at the center but that the selection acting directly on the enzyme loci is weak or zero (se < 0.0038). The stepped pattern implies a barrier to gene flow of 220 (48–415) km. This would substantially delay neutral introgression but would have little effect on advantageous alleles; the two taxa need not evolve independently. Strong selection is needed to maintain such a barrier: hybrid populations must have their mean fitness reduced by a factor of 0.65 (0.60–0.77). This selection must be spread over a large number of loci to account for the concordant patterns and the observed cline widths (N = 300 [80–2,000]).}, author = {Szymura, Jacek and Barton, Nicholas H}, issn = {1558-5646}, journal = {Evolution; International Journal of Organic Evolution}, pages = {1141 -- 1159}, publisher = {Society for the Study of Evolution}, title = {{Genetic analysis of a hybrid zone between the fire-bellied toads Bombina bombina and B. variegata, near Cracow in Southern Poland}}, doi = {10.1111/j.1558-5646.1986.tb05740.x}, volume = {40}, year = {1986}, } @article{3668, abstract = {When two populations which differ at many loci meet, the degree of introgression of alleles across the boundary will depend on the selection acting on each locus (s), the rate of recombination between adjacent loci (r), and the number of loci involved (n). Simple scaling arguments suggest that the behavior of the system should depend on the ratio of selection to recombination (θ = s/r), and on n. This is borne out by mathematical analysis of two demes which exchange individuals at a low rate; when selection is stronger than recombination (θ > 1), the effective selection on each locus is comparable to the total selection over the whole genome (s* ∼ ns). When selection is weaker than recombination (θ < 1), the effective selection is much weaker, but is still stronger than the selection on each locus alone (s* \sim sn20 for small θ). When n is very large, these two regimes are separated by a sharp threshold at θ = 1. The results are extended to two taxa which meet in a continuous habitat; the effective selection pressure, which determines the width of the hybrid zone, behaves in the same way as for the simpler case above. Even when selection is weak compared to recombination, multilocus clines have a sharp step at their center, flanked by tails of introgression in which the alleles behave independently of each other. The set of clines acts as a barrier to gene flow, and it is shown that the barrier is strongest when selection is spread over many loci. The implications of the results for divergence and speciation are discussed.}, author = {Barton, Nicholas H}, issn = {1558-5646}, journal = {Evolution; International Journal of Organic Evolution}, number = {3}, pages = {454 -- 471}, publisher = {Society for the Study of Evolution}, title = {{Multilocus clines}}, doi = {10.2307/2408260}, volume = {37}, year = {1983}, }