Genomic consequences of multiple speciation processes in a stick insect

Diverse geographical modes and mechanisms of speciation are known, and individual speciation genes have now been identified. Despite this progress, genome-wide outcomes of different evolutionary processes during speciation are less understood. Here, we integrate ecological and spatial information, mating trials, transplantation data and analysis of 86 130 single nucleotide polymorphisms (SNPs) in eight populations (28 pairwise comparisons) of Timema cristinae stick insects to test the effects of different factors on genomic divergence in a system undergoing ecological speciation. We find patterns consistent with effects of numerous factors, including geographical distance, gene flow, divergence in host plant use and climate, and selection against maladaptive hybridization (i.e. reinforcement). For example, the number of highly differentiated ‘outlier loci’, allele-frequency clines and the overall distribution of genomic differentiation were recognizably affected by these factors. Although host use has strong effects on phenotypic divergence and reproductive isolation, its effects on genomic divergence were subtler and other factors had pronounced effects. The results demonstrate how genomic data can provide new insights into speciation and how genomic divergence can be complex, yet predictable. Future work could adopt experimental, mapping and functional approaches to directly test which genetic regions are affected by selection and determine their physical location in the genome.     

Ernst Mayr and the integration of geographic and ecological factors in speciation

Mayr's best recognized scientific contributions include the biological species concept and the theory of geographic speciation. In the latter, reproductive isolation evolves as an incidental by‐product of genetic divergence between allopatric populations. Mayr noted that divergent natural selection could accelerate speciation, but also argued that gene flow so strongly retards divergence that, even with selection, non‐allopatric speciation is unlikely. However, current theory and data demonstrate that substantial divergence, and even speciation, in the face of gene flow is possible. Here, I attempt to connect some opposing views about speciation by integrating Mayr's ideas about the roles of ecology and geography in speciation with current data and theory. My central premise is that the speciation process (i.e. divergence) is often continuous, and that the opposing processes of selection and gene flow interact to determine the degree of divergence (i.e. the degree of progress towards the completion of speciation). I first establish that, in the absence of gene flow, divergent selection often promotes speciation. I then discuss how population differentiation in the face of gene flow is common when divergent selection occurs. However, such population differentiation does not always lead to the evolution of discontinuities, strong reproductive isolation, and thus speciation per se. I therefore explore the genetic and ecological circumstances that facilitate speciation in the face of gene flow. For example, particular genetic architectures or ecological niches may tip the balance between selection and gene flow strongly in favour of selection. The circumstances allowing selection to overcome gene flow to the extent that a discontinuity develops, and how often these circumstances occur, are major remaining questions in speciation research. © 2008 The Linnean Society of London, Biological Journal of the Linnean Society, 2008, 95 , 26–46.     

Mutation‐order divergence by sexual selection: diversification of sexual signals in similar environments as a first step in speciation

The origin of species remains a central question, and recent research focuses on the role of ecological differences in promoting speciation. Ecological differences create opportunities for divergent selection (i.e. ‘ecological’ speciation), a Darwinian hypothesis that hardly requires justification. In contrast, ‘mutation‐order’ speciation proposes that, instead of adapting to different environments, populations find different ways to adapt to similar environments, implying that speciation does not require ecological differences. This distinction is critical as it provides an alternative hypothesis to the prevailing view that ecological differences drive speciation. Speciation by sexual selection lies at the centre of debates about the importance of ecological differences in promoting speciation; here, we present verbal and mathematical models of mutation‐order divergence by sexual selection. We develop three general cases and provide a two‐locus population genetic model for each. Results indicate that alternative secondary sexual traits can fix in populations that initially experience similar natural and sexual selection and that divergent traits and preferences can remain stable in the face of low gene flow. This stable divergence can facilitate subsequent divergence that completes or reinforces speciation. We argue that a mutation‐order process could explain widespread diversity in secondary sexual traits among closely related, allopatric species.     

Evidence for spatial clines and mixed geographic modes of speciation for North American cherry‐infesting Rhagoletis (Diptera: Tephritidae) flies

An important criterion for understanding speciation is the geographic context of population divergence. Three major modes of allopatric, parapatric, and sympatric speciation define the extent of spatial overlap and gene flow between diverging populations. However, mixed modes of speciation are also possible, whereby populations experience periods of allopatry, parapatry, and/or sympatry at different times as they diverge. Here, we report clinal patterns of variation for 21 nuclear‐encoded microsatellites and a wing spot phenotype for cherry‐infesting Rhagoletis (Diptera: Tephritidae) across North America consistent with these flies having initially diverged in parapatry followed by a period of allopatric differentiation in the early Holocene. However, mitochondrial DNA (mtDNA) displays a different pattern; cherry flies at the ends of the clines in the eastern USA and Pacific Northwest share identical haplotypes, while centrally located populations in the southwestern USA and Mexico possess a different haplotype. We hypothesize that the mitochondrial difference could be due to lineage sorting but more likely reflects a selective sweep of a favorable mtDNA variant or the spread of an endosymbiont. The estimated divergence time for mtDNA suggests possible past allopatry, secondary contact, and subsequent isolation between USA and Mexican fly populations initiated before the Wisconsin glaciation. Thus, the current genetics of cherry flies may involve different mixed modes of divergence occurring in different portions of the fly''s range. We discuss the need for additional DNA sequencing and quantification of prezygotic and postzygotic reproductive isolation to verify the multiple mixed‐mode hypothesis for cherry flies and draw parallels from other systems to assess the generality that speciation may commonly involve complex biogeographies of varying combinations of allopatric, parapatric, and sympatric divergence.     

Rugged adaptive landscapes shape a complex,sympatric radiation

Strong disruptive ecological selection can initiate speciation, even in the absence of physical isolation of diverging populations. Species evolving under disruptive ecological selection are expected to be ecologically distinct but, at least initially, genetically weakly differentiated. Strong selection and the associated fitness advantages of narrowly adapted individuals, coupled with assortative mating, are predicted to overcome the homogenizing effects of gene flow. Theoretical plausibility is, however, contrasted by limited evidence for the existence of rugged adaptive landscapes in nature. We found evidence for multiple, disruptive ecological selection regimes that have promoted divergence in the sympatric, incipient radiation of ‘sharpfin’ sailfin silverside fishes in ancient Lake Matano (Sulawesi, Indonesia). Various modes of ecological specialization have led to adaptive morphological differences between the species, and differently adapted morphs display significant but incomplete reproductive isolation. Individual fitness and variation in morphological key characters show that disruptive selection shapes a rugged adaptive landscape in this small but complex incipient lake fish radiation.     

The genetic mosaic suggests a new role for hitchhiking in ecological speciation

Early in ecological speciation, the genomically localized effects of divergent selection cause heterogeneity among loci in divergence between incipient species. We call this pattern of genomic variability in divergence the 'genetic mosaic of speciation'. Previous studies have used F(ST) outliers as a way to identify divergently selected genomic regions, but the nature of the relationship between outlier loci and quantitative trait loci (QTL) involved in reproductive isolation has not yet been quantified. Here, we show that F(ST) outliers between a pair of incipient species are significantly clustered around QTL for traits that cause ecologically based reproductive isolation. Around these key QTL, extensive 'divergence hitchhiking' occurs because reduced inter-race mating and negative selection decrease the opportunity for recombination between chromosomes bearing different locally adapted QTL alleles. Divergence hitchhiking is likely to greatly increase the opportunity for speciation in populations that are sympatric, regardless of whether initial divergence was sympatric or allopatric. Early in ecological speciation, analyses of population structure, gene flow or phylogeography based on different random or arbitrarily chosen neutral markers should be expected to conflict--only markers in divergently selected genomic regions will reveal the evolutionary history of adaptive divergence and ecologically based reproductive isolation. Species retain mosaic genomes for a very long time, and gene exchange in hybrid zones can vary dramatically among loci. However, in hybridizing species, the genomic regions that affect ecologically based reproductive isolation are difficult to distinguish from regions that have diverged for other reasons.     

Divergence,gene flow,and speciation in eight lineages of trans‐Beringian birds

Determining how genetic diversity is structured between populations that span the divergence continuum from populations to biological species is key to understanding the generation and maintenance of biodiversity. We investigated genetic divergence and gene flow in eight lineages of birds with a trans‐Beringian distribution, where Asian and North American populations have likely been split and reunited through multiple Pleistocene glacial cycles. Our study transects the speciation process, including eight pairwise comparisons in three orders (ducks, shorebirds and passerines) at population, subspecies and species levels. Using ultraconserved elements (UCEs), we found that these lineages represent conditions from slightly differentiated populations to full biological species. Although allopatric speciation is considered the predominant mode of divergence in birds, all of our best divergence models included gene flow, supporting speciation with gene flow as the predominant mode in Beringia. In our eight lineages, three were best described by a split‐migration model (divergence with gene flow), three best fit a secondary contact scenario (isolation followed by gene flow), and two showed support for both models. The lineages were not evenly distributed across a divergence space defined by gene flow (M) and differentiation (F_ST), instead forming two discontinuous groups: one with relatively shallow divergence, no fixed single nucleotide polymorphisms (SNPs), and high rates of gene flow between populations; and the second with relatively deeply divergent lineages, multiple fixed SNPs, and low gene flow. Our results highlight the important role that gene flow plays in avian divergence in Beringia.     

DYNAMICS OF HYBRID INCOMPATIBILITY IN GENE NETWORKS IN A CONSTANT ENVIRONMENT

After an ancestral population splits into two allopatric populations, different mutations may fix in each. When pairs of mutations are brought together in a hybrid offspring, epistasis may cause reduced fitness. Such pairs are known as Bateson–Dobzhansky–Muller (BDM) incompatibilities. A well-known model of BDM incompatibility due to Orr suggests that the fitness load on hybrids should initially accelerate, and continue to increase as the number of potentially incompatible substitutions increases (the "snowball effect"). In the gene networks model, which violates a key assumption of Orr's model (independence of fixation probabilities), the snowball effect often does not occur. Instead, we describe three possible dynamics in a constant environment: (1) Stabilizing selection can constrain two allopatric populations to remain near-perfectly compatible. (2) Despite constancy of environment, punctuated evolution may obtain; populations may experience rare adaptations asynchronously, permitting incompatibility. (3) Despite stabilizing selection, developmental system drift may permit genetic change, allowing two populations to drift in and out of compatibility. We reinterpret Orr's model in terms of genetic distance. We extend Orr's model to the finite loci case, which can limit incompatibility. Finally, we suggest that neutral evolution of gene regulation in nature, to the point of speciation, is a distinct possibility.     

Molecular correlates of reproductive isolation

Evolution of reproductive isolation as a byproduct of genetic divergence in isolated populations is the dominant (albeit not exclusive) mode of speciation in sexual animals. But little is known about the factors linking speciation to general divergence. Several authors have argued that allopatric speciation should proceed more rapidly if isolated populations also experience divergent selection. Reproductive isolation between allopatric populations is not subject to direct selection; it can accumulate only by random drift or as a fortuitous byproduct of selection on other traits. Here I present a novel analysis of published data, demonstrating that pre- and postmating isolation of Drosophila species are more tightly correlated with allozyme divergence than with silent DNA divergence. Inasmuch as proteins are more subject to the action of natural selection than are silent DNA polymorphisms, this result provides broad support for a model of selection-mediated allopatric speciation.     

Disruptive ecological selection on a mating cue

Adaptation to divergent ecological niches can result in speciation. Traits subject to disruptive selection that also contribute to non-random mating will facilitate speciation with gene flow. Such ‘magic’ or ‘multiple-effect’ traits may be widespread and important for generating biodiversity, but strong empirical evidence is still lacking. Although there is evidence that putative ecological traits are indeed involved in assortative mating, evidence that these same traits are under divergent selection is considerably weaker. Heliconius butterfly wing patterns are subject to positive frequency-dependent selection by predators, owing to aposematism and Müllerian mimicry, and divergent colour patterns are used by closely related species to recognize potential mates. The amenability of colour patterns to experimental manipulation, independent of other traits, presents an excellent opportunity to test their role during speciation. We conducted field experiments with artificial butterflies, designed to match natural butterflies with respect to avian vision. These were complemented with enclosure trials with live birds and real butterflies. Our experiments showed that hybrid colour-pattern phenotypes are attacked more frequently than parental forms. For the first time, we demonstrate disruptive ecological selection on a trait that also acts as a mating cue.     

Adaptive chromosomal divergence driven by mixed geographic mode of evolution

Chromosomal inversions are ubiquitous in nature and of great significance for understanding adaptation and speciation. Inversions were the first markers used to investigate the genetic structure of natural populations, leading to the concept of coadapted gene complexes and theories concerning founder effects and genetic drift in small populations. However, we still lack elements of a general theory accounting for the origins and distribution of inversions in nature. Here, we use computer simulations to show that a "mixed geographic mode" of evolution involving allopatric separation of populations followed by secondary contact and gene flow generates chromosomal divergence by natural selection under wider conditions than previous hypotheses. This occurs because inversions arising in allopatry contain a full complement of locally adapted genes. Once gene flow ensues, reduced recombination within inversions keeps these favorable genotypic combinations intact, resulting in inverted genomic regions being favored over collinear regions. This process allows inversions to establish to high frequencies. Our model can account for several classic patterns in the geographic distribution of inversions and highlights how selection on standing genetic variation allows rapid chromosomal evolution without the waiting time for new mutations. As inversion differences often separate closely related taxa, mixed modes of divergence could be common.     

Consequences of divergence and introgression for speciation in Andean cloud forest birds

Divergence with gene flow is well documented and reveals the influence of ecological adaptation on speciation. Yet, it remains intuitive that gene exchange inhibits speciation in many scenarios, particularly among ecologically similar populations. The influence of gene flow on the divergence of populations facing similar selection pressures has received less empirical attention than scenarios where differentiation is coupled with local environmental adaptation. I used a paired study design to test the influence of genomic divergence and introgression on plumage differentiation between ecologically similar allopatric replacements of Andean cloud forest birds. Through analyses of short‐read genome‐wide sequences from over 160 individuals in 16 codistributed lineages, I found that plumage divergence is associated with deep genetic divergence, implicating a prominent role of geographic isolation in speciation. By contrast, lineages that lack plumage divergence across the same geographic barrier are more recently isolated or exhibit a signature of secondary genetic introgression, indicating a negative relationship between gene flow and divergence in phenotypic traits important to speciation. My results suggest that the evolutionary outcomes of cycles of isolation and divergence in this important theatre of biotic diversification are sensitive to time spent in the absence of gene flow.     

Hybrid fitness in a locally adapted parasite

The parasite (Red Queen) hypothesis for the maintenance of sexual reproduction and genetic diversity assumes that host-parasite interactions result from tight genetic specificity. Hence, hybridization between divergent parasite populations would be expected to disrupt adaptive gene combinations, leading to reduced infectivity on exposure to parental sympatric hosts, as long as gene effects are nonadditive. In contrast, hybridization would not cause reduced infectivity on allopatric hosts unless the divergent parasite populations possess alleles that are intrinsically incompatible when they are combined. In three different experiments, we compared the infectivity of locally adapted parasite (trematode) populations with that of F(1) hybrid parasites when exposed to host (snail) populations that were sympatric to one of the two parasite populations. We tested for intrinsic genetic incompatibilities in two experiments by including one host population that was allopatric to both parasite populations. As predicted, when the target host populations were sympatric to the parasite populations, the hybrids were significantly less infective than the parental average, while hybrid parasites on allopatric hosts were not, thereby ruling out intrinsic genetic incompatibilities. The results are consistent with nonadditive gene effects and tightly specific host-driven selection underlying parasite divergence, as envisioned by coevolutionary theory and the Red Queen hypothesis.     

基因流存在条件下的物种形成研究述评:生殖隔离机制进化

物种形成过程是生物多样性形成的基础, 长期以来一直是进化生物学的中心议题之一。传统的异域物种形成理论认为, 地理隔离是物种分化的主要决定因子, 物种形成只有在种群之间存在地理隔离的情况下才能发生。近年来, 随着种群基因组学的发展和溯祖理论分析方法的完善, 种群间存在基因流情况下的物种形成成为进化生物学领域新的研究焦点。物种形成过程中是否有基因流的发生?基因流如何影响物种的形成与分化?基因流存在条件下物种形成的生殖隔离机制是什么?根据已发表的相关文献资料, 作者综述了当前物种形成研究中基因流的时间和空间分布模式、基因流对物种分化的影响以及生殖隔离机制形成等问题, 指出基因流存在条件下的物种形成可能是自然界普遍发生的一种模式。     

Divergent host-plant use promotes reproductive isolation among cynipid gall wasp populations

Ecological speciation occurs when reproductive isolation evolves as a consequence of divergent natural selection among environments. A direct prediction of this process is that ecologically divergent pairs of populations will exhibit greater reproductive isolation than ecologically similar pairs of populations. By comparing allopatric populations of the cynipid gall wasp Belonocnema treatae infesting Quercus virginiana and Quercus geminata, we tested the role that divergent host use plays in generating ecological divergence and sexual isolation. We found differences in body size and gall structure associated with divergent host use, but no difference in neutral genetic divergence between populations on the same or different host plant. We observed significant assortative mating between populations from alternative host plants but not between allopatric populations on the same host plant. Thus, we provide evidence that divergent host use promotes speciation among gall wasp populations.     

Vicariance divergence and gene flow among islet populations of an endemic lizard

Allopatry and allopatric speciation can arise through two different mechanisms: vicariance or colonization through dispersal. Distinguishing between these different allopatric mechanisms is difficult and one of the major challenges in biogeographical research. Here, we address whether allopatric isolation in an endemic island lizard is the result of vicariance or dispersal. We estimated the amount and direction of gene flow during the divergence of isolated islet populations and subspecies of the endemic Skyros wall lizard Podarcis gaigeae, a phenotypically variable species that inhabits a major island and small islets in the Greek archipelago. We applied isolation-with-migration models to estimate population divergence times, population sizes and gene flow between islet-mainland population pairs. Divergence times were significantly correlated with independently estimated geological divergence times. This correlation strongly supports a vicariance scenario where islet populations have sequentially become isolated from the major island. We did not find evidence for significant gene flow within P. g. gaigeae. However, gene-flow estimates from the islet to the mainland populations were positively affected by islet area and negatively by distance between the islet and mainland. We also found evidence for gene flow from one subspecies (P. g. weigandi) into another (P. g. gaigeae), but not in the other direction. Ongoing gene flow between the subspecies suggests that even in this geographically allopatric scenario with the sea posing a strong barrier to dispersal, divergence with some gene flow is still feasible.     

Parasite-mediated selection drives an immunogenetic trade-off in plains zebras (Equus quagga)

Pathogen evasion of the host immune system is a key force driving extreme polymorphism in genes of the major histocompatibility complex (MHC). Although this gene family is well characterized in structure and function, there is still much debate surrounding the mechanisms by which MHC diversity is selectively maintained. Many studies have investigated relationships between MHC variation and specific pathogens, and have found mixed support for and against the hypotheses of heterozygote advantage, frequency-dependent or fluctuating selection. Few, however, have focused on the selective effects of multiple parasite types on host immunogenetic patterns. Here, we examined relationships between variation in the equine MHC gene, ELA-DRA, and both gastrointestinal (GI) and ectoparasitism in plains zebras (Equus quagga). Specific alleles present at opposing population frequencies had antagonistic effects, with rare alleles associated with increased GI parasitism and common alleles with increased tick burdens. These results support a frequency-dependent mechanism, but are also consistent with fluctuating selection. Maladaptive GI parasite ‘susceptibility alleles’ were reduced in frequency, suggesting that these parasites may play a greater selective role at this locus. Heterozygote advantage, in terms of allele mutational divergence, also predicted decreased GI parasite burden in genotypes with a common allele. We conclude that an immunogenetic trade-off affects resistance/susceptibility to parasites in this system. Because GI and ectoparasites do not directly interact within hosts, our results uniquely show that antagonistic parasite interactions can be indirectly modulated through the host immune system. This study highlights the importance of investigating the role of multiple parasites in shaping patterns of host immunogenetic variation.     

Pervasive genetic associations between traits causing reproductive isolation in Heliconius butterflies

Ecological speciation proceeds through the accumulation of divergent traits that contribute to reproductive isolation, but in the face of gene flow traits that characterize incipient species may become disassociated through recombination. Heliconius butterflies are well known for bright mimetic warning patterns that are also used in mate recognition and cause both pre- and post-mating isolation between divergent taxa. Sympatric sister taxa representing the final stages of speciation, such as Heliconius cydno and Heliconius melpomene, also differ in ecology and hybrid fertility. We examine mate preference and sterility among offspring of crosses between these species and demonstrate the clustering of Mendelian colour pattern loci and behavioural loci that contribute to reproductive isolation. In particular, male preference for red patterns is associated with the locus responsible for the red forewing band. Two further colour pattern loci are associated, respectively, with female mating outcome and hybrid sterility. This genetic architecture in which ‘speciation genes’ are clustered in the genome can facilitate two controversial models of speciation, namely divergence in the face of gene flow and hybrid speciation.     

A neutral view of the evolving genomic architecture of speciation

Analyses of genomewide polymorphism data have begun to shed light on speciation and adaptation. Genome scans to identify regions of the genome that are unusually different between populations or species, possibly due to divergent natural or sexual selection, are widespread in speciation genomics. Theoretical and empirical work suggests that such outlier regions may grow faster than linearly during speciation with gene flow due to a rapid transition between low and high reproductive isolation. We investigate whether this pattern could be attributed to neutral processes by simulating genomes under neutral evolution with varying amounts and timing of gene flow. Under both neutral evolution and divergent selection, simulations with little or no gene flow, or with a long allopatric period after its cessation, resulted in faster than linear growth of the proportion of the genome lying in outlier regions. Without selection, higher recent gene flow erased differentiation; with divergent selection, these same scenarios produced nonlinear growth to a plateau. Our results suggest that, given a history of gene flow, the growth of the divergent genome is informative about selection during divergence, but that in many scenarios, this pattern does not easily distinguish neutral and non‐neutral processes during speciation with gene flow.