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.     

Establishment of new mutations under divergence and genome hitchhiking

Theoretical models addressing genome-wide patterns of divergence during speciation are needed to help us understand the evolutionary processes generating empirical patterns. Here, we examine a critical issue concerning speciation-with-gene flow: to what degree does physical linkage (r < 0.5) of new mutations to already diverged genes aid the build-up of genomic islands of differentiation? We used simulation and analytical approaches to partition the probability of establishment for a new divergently selected mutation when the mutation (i) is the first to arise in an undifferentiated genome (the direct effect of selection), (ii) arises unlinked to any selected loci (r = 0.5), but within a genome that has some already diverged genes (the effect of genome-wide reductions in gene flow for facilitating divergence, which we term 'genome hitchhiking'), and (iii) arises in physical linkage to a diverged locus (divergence hitchhiking). We find that the strength of selection acting directly on a new mutation is generally the most important predictor for establishment, with divergence and genomic hitchhiking having smaller effects. We outline the specific conditions under which divergence and genome hitchhiking can aid mutation establishment. The results generate predictions about genome divergence at different points in the speciation process and avenues for further work.     

Widespread yet heterogeneous genomic divergence

Genetic differentiation during adaptive divergence and speciation is heterogeneous among genomic regions. Some regions can be highly differentiated between populations, for example, because they harbour genes under divergent selection or those causing reproductive isolation and thus are resistant to gene flow. Other regions might be homogenized by gene flow and thus weakly differentiated. Debates persist about the number of differentiated regions expected under divergence with gene flow, and their causes, size, and genomic distribution. In this issue of Molecular Ecology, a study of freshwater stickleback used next-generation sequencing to shed novel insight into these issues (Roesti et al. 2012). Many genomic regions distributed across the genome were strongly differentiated, indicating divergence with gene flow can involve a greater number of loci than often thought. Nonetheless, differentiation of some regions, such as those near the centre of chromosomes where recombination is reduced, was strongly accentuated over others. Thus, divergence was widespread yet highly heterogeneous across the genome. Moreover, different population pairs varied in patterns of differentiation, illustrating how genomic divergence builds up across stages of the speciation process. The study demonstrates how variation in different evolutionary processes, such as selection and recombination rate, can combine to result in similar genomic patterns. Future work could focus on teasing apart the contributions of different processes for causing differentiation, a task facilitated by experimental manipulations.     

The genomics of speciation-with-gene-flow

The emerging field of speciation genomics is advancing our understanding of the evolution of reproductive isolation from the individual gene to a whole-genome perspective. In this new view it is important to understand the conditions under which 'divergence hitchhiking' associated with the physical linkage of gene regions, versus 'genome hitchhiking' associated with reductions in genome-wide rates of gene flow caused by selection, can enhance speciation-with-gene-flow. We describe here a theory predicting four phases of speciation, defined by changes in the relative effectiveness of divergence and genome hitchhiking, and review empirical data in light of the theory. We outline future directions, emphasizing the need to couple next-generation sequencing with selection, transplant, functional genomics, and mapping studies. This will permit a natural history of speciation genomics that will help to elucidate the factors responsible for population divergence and the roles that genome structure and different forms of hitchhiking play in facilitating the genesis of new biodiversity.     

Genomics of Divergence along a Continuum of Parapatric Population Differentiation

The patterns of genomic divergence during ecological speciation are shaped by a combination of evolutionary forces. Processes such as genetic drift, local reduction of gene flow around genes causing reproductive isolation, hitchhiking around selected variants, variation in recombination and mutation rates are all factors that can contribute to the heterogeneity of genomic divergence. On the basis of 60 fully sequenced three-spined stickleback genomes, we explore these different mechanisms explaining the heterogeneity of genomic divergence across five parapatric lake and river population pairs varying in their degree of genetic differentiation. We find that divergent regions of the genome are mostly specific for each population pair, while their size and abundance are not correlated with the extent of genome-wide population differentiation. In each pair-wise comparison, an analysis of allele frequency spectra reveals that 25–55% of the divergent regions are consistent with a local restriction of gene flow. Another large proportion of divergent regions (38–75%) appears to be mainly shaped by hitchhiking effects around positively selected variants. We provide empirical evidence that alternative mechanisms determining the evolution of genomic patterns of divergence are not mutually exclusive, but rather act in concert to shape the genome during population differentiation, a first necessary step towards ecological speciation.     

Genome scans and elusive candidate genes: detecting the variation that matters for speciation

Next‐generation sequencing is providing us with vast amounts of genetic data, yet we are currently struggling in our attempts to make sense of them. In particular, it has proven difficult to link phenotypic divergence and speciation to genome level divergence. In the current issue of Molecular Ecology, Ruegg et al. ( 2014 ) present new empirical results from two closely related bird taxa. They use a promising approach combining genome scan and candidate gene analysis. Their results suggest that we may have been looking in vain for candidate speciation genes by focusing only on genes found within genomic islands of divergence. This is because genes important in divergence and speciation may not be detected by genome scans and because features of the genomic architecture per se may have a large effect on the pattern of genome divergence.     

Mutation and selection processes regulating short tandem repeats give rise to genetic and phenotypic diversity across species

Short tandem repeats (STRs) are units of 1–6 bp that repeat in a tandem fashion in DNA. Along with single nucleotide polymorphisms and large structural variations, they are among the major genomic variants underlying genetic, and likely phenotypic, divergence. STRs experience mutation rates that are orders of magnitude higher than other well-studied genotypic variants. Frequent copy number changes result in a wide range of alleles, and provide unique opportunities for modulating complex phenotypes through variation in repeat length. While classical studies have identified key roles of individual STR loci, the advent of improved sequencing technology, high-quality genome assemblies for diverse species, and bioinformatics methods for genome-wide STR analysis now enable more systematic study of STR variation across wide evolutionary ranges. In this review, we explore mutation and selection processes that affect STR copy number evolution, and how these processes give rise to varying STR patterns both within and across species. Finally, we review recent examples of functional and adaptive changes linked to STRs.     

Population genomic analyses reveal possible drivers of population divergence

Recent advances in sequencing technology and efficiency enable new and improved methods to investigate how populations diverge and species evolve. Fungi have relatively small and simple genomes and can often be cultured in the laboratory. Fungal populations can thus be sequenced for a relatively low cost, which makes them ideal for population genomic analyses. In several recent population genomic studies, wild populations of fungal model organisms and human pathogens have been analysed, for example Neurospora crassa (Ellison et al. 2011 ), Saccharomyces uvarum (Almeida et al. 2014 ), Coccidioides spp. (Neafsey et al. 2010 ) and Cryptococcus gatti (Engelthaler et al. 2014 ). In this issue of Molecular Ecology, Branco et al. ( 2015 ) apply population genomic tools to understand population divergence and adaptation in a symbiotic (mycorrhizal) fungus. This study exemplifies the possibilities of diving deeper into the genomic features involved in population divergence and speciation, also for nonmodel organisms, and how molecular and analytical tools will improve our understanding of the patterns and mechanisms that underlie adaptation to habitats, population divergence and dispersal limitation of fungi.     

Genome-wide patterns of divergence during speciation: the lake whitefish case study

The nature, size and distribution of the genomic regions underlying divergence and promoting reproductive isolation remain largely unknown. Here, we summarize ongoing efforts using young (12 000 yr BP) species pairs of lake whitefish (Coregonus clupeaformis) to expand our understanding of the initial genomic patterns of divergence observed during speciation. Our results confirmed the predictions that: (i) on average, phenotypic quantitative trait loci (pQTL) show higher F(ST) values and are more likely to be outliers (and therefore candidates for being targets of divergent selection) than non-pQTL markers; (ii) large islands of divergence rather than small independent regions under selection characterize the early stages of adaptive divergence of lake whitefish; and (iii) there is a general trend towards an increase in terms of numbers and size of genomic regions of divergence from the least (East L.) to the most differentiated species pair (Cliff L.). This is consistent with previous estimates of reproductive isolation between these species pairs being driven by the same selective forces responsible for environment specialization. Altogether, dwarf and normal whitefish species pairs represent a continuum of both morphological and genomic differentiation contributing to ecological speciation. Admittedly, much progress is still required to more finely map and circumscribe genomic islands of speciation. This will be achieved through the use of next generation sequencing data but also through a better quantification of phenotypic traits moulded by selection as organisms adapt to new environmental conditions.     

The genes underlying the process of speciation

The long-standing goal of finding genes causing reproductive isolation is being achieved. To better link the genetics with the process of speciation, we propose that 'speciation gene' be defined as any gene contributing to the evolution of reproductive isolation. Characterizing a speciation gene involves establishing that the gene affects a component of reproductive isolation; demonstrating that divergence at the locus occurred before completion of speciation; and quantifying the effect size of the gene (i.e. the increase in total reproductive isolation caused by its divergence). Review of a sample of candidate speciation genes found that few meet these criteria. Improved characterization of speciation genes will clarify how numerous they are, their properties and how they affect genome-wide patterns of divergence.     

Can genomic data alone tell us whether speciation happened with gene flow?

The allopatric model, which requires a period of geographical isolation for speciation to complete, has been the standard model in the modern era. Recently, “speciation with gene flow” has been widely discussed in relation to the model of “strict allopatry” and the level of DNA divergence across genomic regions. We wish to caution that genomic data by themselves may only permit the rejection of the simplest form of allopatry. Even a slightly more complex and realistic model that starts with subdivided populations would be impossible to reject by the genomic data alone. To resolve this central issue of speciation, other forms of observations such as the sequencing of reproductive isolation genes or the identification of geographical barrier(s) will be necessary.     

The shape of things to come in the study of the origin of species?

Perhaps Darwin would agree that speciation is no longer the mystery of mysteries that it used to be. It is now generally accepted that evolution by natural selection can contribute to ecological adaptation, resulting in the evolution of reproductive barriers and, hence, to the evolution of new species (Schluter & Conte 2009 ; Meyer 2011 ; Nosil 2012 ). From genes that encode silencing proteins that cause infertility in hybrid mice (Mihola et al. 2009 ), to segregation distorters linked to speciation in fruit flies (Phadnis & Orr 2009 ), or pollinator‐mediated selection on flower colour alleles driving reinforcement in Texan wildflowers (Hopkins & Rausher 2012 ), characterization of the genes that drive speciation is providing clues to the origin of species (Nosil & Schluter 2011 ). It is becoming apparent that, while recent work continues to overturn historical ideas about sympatric speciation (e.g. Barluenga et al. 2006 ), ecological circumstances strongly influence patterns of genomic divergence, and ultimately the establishment of reproductive isolation when gene flow is present (Elmer & Meyer 2011 ). Less clear, however, are the genetic mechanisms that cause speciation, particularly when ongoing gene flow is occurring. Now, in this issue, Franchini et al. ( 2014 ) employ a classic genetic mapping approach augmented with new genomic tools to elucidate the genomic architecture of ecologically divergent body shapes in a pair of sympatric crater lake cichlid fishes. From over 450 segregating SNPs in an F2 cross, 72 SNPs were linked to 11 QTL associated with external morphology measured by means of traditional and geometric morphometrics. Annotation of two highly supported QTL further pointed to genes that might contribute to ecological divergence in body shape in Midas cichlids, overall supporting the hypothesis that genomic regions of large phenotypic effect may be contributing to early‐stage divergence in Midas cichlids.     

Selection on Plant Male Function Genes Identifies Candidates for Reproductive Isolation of Yellow Monkeyflowers

Understanding the genetic basis of reproductive isolation promises insight into speciation and the origins of biological diversity. While progress has been made in identifying genes underlying barriers to reproduction that function after fertilization (post-zygotic isolation), we know much less about earlier acting pre-zygotic barriers. Of particular interest are barriers involved in mating and fertilization that can evolve extremely rapidly under sexual selection, suggesting they may play a prominent role in the initial stages of reproductive isolation. A significant challenge to the field of speciation genetics is developing new approaches for identification of candidate genes underlying these barriers, particularly among non-traditional model systems. We employ powerful proteomic and genomic strategies to study the genetic basis of conspecific pollen precedence, an important component of pre-zygotic reproductive isolation among yellow monkeyflowers (Mimulus spp.) resulting from male pollen competition. We use isotopic labeling in combination with shotgun proteomics to identify more than 2,000 male function (pollen tube) proteins within maternal reproductive structures (styles) of M. guttatus flowers where pollen competition occurs. We then sequence array-captured pollen tube exomes from a large outcrossing population of M. guttatus, and identify those genes with evidence of selective sweeps or balancing selection consistent with their role in pollen competition. We also test for evidence of positive selection on these genes more broadly across yellow monkeyflowers, because a signal of adaptive divergence is a common feature of genes causing reproductive isolation. Together the molecular evolution studies identify 159 pollen tube proteins that are candidate genes for conspecific pollen precedence. Our work demonstrates how powerful proteomic and genomic tools can be readily adapted to non-traditional model systems, allowing for genome-wide screens towards the goal of identifying the molecular basis of genetically complex traits.     

Utility of closely related taxa for genetic studies of adaptive variation and speciation: Current state and perspectives in plants with focus on forest tree species

Studies of adaptation and speciation have greatly benefited from rapid progress of DNA sequencing and genotyping technologies. Comparative genomics of closely related taxa has great potential to advance evolutionary research on genetic architecture of adaptive traits and reproductive isolation. Such studies that utilized closely related plant species and ecotypes have already provided some important insights into genomic regions and/or genes that are potentially involved in local adaptation and speciation. The choice of an appropriate species model for such research is crucial. The paper discusses current approaches used to reveal the patterns of intra‐ and interspecific divergence due to natural selection. Its outcomes in herbaceous plants and forest trees are briefly summarized and compared to reveal general regularities concerning evolutionary processes. We then highlight the importance of multispecies studies and discuss the utility of several related pine taxa as fine candidates for evolutionary inferences. Genetically similar but ecologically and phenotypically diverged taxa seem a promising study system to search for genomic patterns of speciation and adaptive variation.     

Extensive linkage disequilibrium and parallel adaptive divergence across threespine stickleback genomes

Population genomic studies are beginning to provide a more comprehensive view of dynamic genome-scale processes in evolution. Patterns of genomic architecture, such as genomic islands of increased divergence, may be important for adaptive population differentiation and speciation. We used next-generation sequencing data to examine the patterns of local and long-distance linkage disequilibrium (LD) across oceanic and freshwater populations of threespine stickleback, a useful model for studies of evolution and speciation. We looked for associations between LD and signatures of divergent selection, and assessed the role of recombination rate variation in generating LD patterns. As predicted under the traditional biogeographic model of unidirectional gene flow from ancestral oceanic to derived freshwater stickleback populations, we found extensive local and long-distance LD in fresh water. Surprisingly, oceanic populations showed similar patterns of elevated LD, notably between large genomic regions previously implicated in adaptation to fresh water. These results support an alternative biogeographic model for the stickleback radiation, one of a metapopulation with appreciable bi-directional gene flow combined with strong divergent selection between oceanic and freshwater populations. As predicted by theory, these processes can maintain LD within and among genomic islands of divergence. These findings suggest that the genomic architecture in oceanic stickleback populations may provide a mechanism for the rapid re-assembly and evolution of multi-locus genotypes in newly colonized freshwater habitats, and may help explain genetic mapping of parallel phenotypic variation to similar loci across independent freshwater populations.     

Mapping the genomic architecture of ecological speciation in the wild: does linkage disequilibrium hold the key?

The hunt for the genes underlying ecological speciation has now closed in on a number of candidates, but making the link from genotype to phenotype continues to pose a significant challenge. This is partly because genetic studies in many systems remain impeded by long generation times or an inability to perform controlled crosses. Now, in this issue of Molecular Ecology, Malek et al. (2012) demonstrate the utility of a novel admixture mapping approach that can be used to identify genomic regions contributing to adaptive trait divergence between natural populations. Remarkably, they validate their approach by mapping traits associated with mate choice in a wild limnetic and benthic threespine stickleback (Gasterosteus aculeatus) species pair, finding several loci associated with male nuptial coloration and shape. While this study benefited from tried‐and‐true microsatellites in a well‐characterized species with a detailed genetic map (and genome sequence), the field is quickly moving towards the use of next‐generation sequencing, especially for nonmodel systems. The ability to characterize molecular polymorphisms for any system suggests that molecular ecologists working on virtually any species may benefit from applying Malek et al.'s approach, if naturally admixed populations are available.     

Constraints on speciation suggested by comparing lake-stream stickleback divergence across two continents

Adaptation to ecologically distinct environments can coincide with the emergence of reproductive barriers. The outcome of this process is highly variable and can range along a continuum from weak population differentiation all the way to complete, genome-wide divergence. The factors determining how far diverging taxa will move along this continuum remain poorly understood but are most profitably investigated in taxa under replicate divergence. Here, we explore determinants of progress towards speciation by comparing phenotypic and molecular divergence within young (<150 years) lake-stream stickleback pairs from Central Europe to divergence in older (thousands of years) archetypal lake-stream pairs from Vancouver Island, Canada. We generally find relatively weak divergence in most aspects of foraging morphology (gill raker number, body shape) in the European pairs, although substantial adaptive divergence is seen in gill raker length. Combined with striking overall phenotypic differences between the continents, this argues for genetic and time constraints on adaptive divergence in the European pairs. The European lake-stream pairs also do not display the strong habitat-related differentiation in neutral (microsatellite) markers seen in the Canadian watersheds. This indicates either the lack of strong reproductive barriers owing to weak adaptive divergence, or alternatively that neutral markers are poorly suited for detecting reproductive barriers if these emerge rapidly. Overall, our comparative approach suggests constraints on speciation due to genetic architecture and limited time for divergence. The relative importance of these factors remains to be quantified by future investigation.