Landscape drivers of genomic diversity and divergence in woodland Eucalyptus

Spatial genetic patterns are influenced by numerous factors, and they can vary even among coexisting, closely related species due to differences in dispersal and selection. Eucalyptus (L'Héritier 1789; the “eucalypts”) are foundation tree species that provide essential habitat and modulate ecosystem services throughout Australia. Here we present a study of landscape genomic variation in two woodland eucalypt species, using whole‐genome sequencing of 388 individuals of Eucalyptus albens and Eucalyptus sideroxylon. We found exceptionally high genetic diversity (π ≈ 0.05) and low genome‐wide, interspecific differentiation (F_ST = 0.15) and intraspecific differentiation between localities (F_ST ≈ 0.01–0.02). We found no support for strong, discrete population structure, but found substantial support for isolation by geographic distance (IBD) in both species. Using generalized dissimilarity modelling, we identified additional isolation by environment (IBE). Eucalyptus albens showed moderate IBD, and environmental variables have a small but significant amount of additional predictive power (i.e. IBE). Eucalyptus sideroxylon showed much stronger IBD and moderate IBE. These results highlight the vast adaptive potential of these species and set the stage for testing evolutionary hypotheses of interspecific adaptive differentiation across environments.     

The genetic basis of evolution, adaptation and speciation

A primary question in biology concerns the genetic basis of the evolution of novel traits, often in response to environmental changes, and how this can subsequently cause species isolation. This topic was the focus of the symposium on the Genetics of Speciation and Evolution at the annual meeting of the Canadian Society for Ecology and Evolution, held in Banff in May 2011. The presentations revealed some of the rapid advances being made in understanding the genetic basis of adaptation and speciation, as well as the elegant interplay between an organism's genetic complement and the environment that organism experiences.     

Making inferences about speciation using sophisticated statistical genetics methods: look before you leap

Understanding speciation is a fundamental aim of evolutionary biology and a very challenging one. Speciation can be viewed as the dynamics of genetic differentiation between populations resulting in substantial reproductive isolation (Gavrilets 2003). It was generally accepted that very small levels of migration prevent genetic differentiation among populations and, therefore, speciation. However, recent theoretical work showed that sympatric speciation is possible (Gavrilets 2003). Nevertheless, providing empirical evidence that gene flow occurred during speciation is challenging because several gene flow scenarios can explain observed patterns of genetic differentiation. Positive migration rate estimates alone do not prove ongoing gene flow during divergence. We also need to know whether migration took place before, during or after speciation. There is no statistical method specifically aimed at estimating gene flow timing, but several studies used the isolation with migration model (Hey & Nielsen 2004, 2007; Hey 2010b) to estimate this parameter and make inferences about speciation scenarios. It is tempting to use statistical methods to estimate important evolutionary parameters even if they do not appear explicitly in the inference model. Nevertheless, before doing so, we need to determine whether they can provide reliable results. In this issue of Molecular Ecology, Strasburg and Rieseberg (2011) present a simulation study that examines the degree to which gene flow timing estimates obtained from IMa2 (Hey 2010b) can be used to make inferences about speciation mode. Their results are sobering; gene flow timing estimates obtained in this way are not reliable and cannot be used to unequivocally establish if gene flow was ongoing during speciation.     

Divergence hitchhiking and the spread of genomic isolation during ecological speciation-with-gene-flow

In allopatric populations, geographical separation simultaneously isolates the entire genome, allowing genetic divergence to accumulate virtually anywhere in the genome. In sympatric populations, however, the strong divergent selection required to overcome migration produces a genetic mosaic of divergent and non-divergent genomic regions. In some recent genome scans, each divergent genomic region has been interpreted as an independent incidence of migration/selection balance, such that the reduction of gene exchange is restricted to a few kilobases around each divergently selected gene. I propose an alternative mechanism, 'divergence hitchhiking' (DH), in which divergent selection can reduce gene exchange for several megabases around a gene under strong divergent selection. Not all genes/markers within a DH region are divergently selected, yet the entire region is protected to some degree from gene exchange, permitting genetic divergence from mechanisms other than divergent selection to accumulate secondarily. After contrasting DH and multilocus migration/selection balance (MM/SB), I outline a model in which genomic isolation at a given genomic location is jointly determined by DH and genome-wide effects of the progressive reduction in realized migration, then illustrate DH using data from several pairs of incipient species in the wild.     

History repeats itself: genomic divergence in copepods

Press stop, erase everything from now till some arbitrary time in the past and start recording life as it evolves once again. Would you see the same tape of life playing itself over and over, or would a different story unfold every time? The late Steven Jay Gould called this experiment replaying the tape of life and argued that any replay of the tape would lead evolution down a pathway radically different from the road actually taken (Gould 1989). This thought experiment has puzzled evolutionary biologists for a long time: how repeatable are evolutionary events? And if history does indeed repeat itself, what are the factors that may help us predict the path taken? A powerful means to address these questions at a small evolutionary scale is to study closely related populations that have evolved independently, under similar environmental conditions. This is precisely what Pereira et al. ( 2016 ) set out to do using marine copepods Tigriopus californicus, and present their results in this issue of Molecular Ecology. They show that evolution can be repeatable and even partly predictable, at least at the molecular level. As expected from theory, patterns of divergence were shaped by natural selection. At the same time, strong genetic drift due to small population sizes also constrained evolution down a similar evolutionary road, and probably contributed to repeatable patterns of genomic divergence.     

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.     

Evolutionary potential in the wild: more than meets the eye

The genus Aquilegia consists of 60–70 perennial plant species widely distributed throughout the northern hemisphere. Its flowers have a delicate and ornamental appearance that makes them a favourite of gardeners. In this genus, adaptive radiations for both floral and vegetative traits have occurred. These adaptive radiations, and the key phylogenetic placement of Aquilegia between Arabidopsis and rice, make this genus a ‘model system’ for plant evolution ( Kramer 2009 ). In this issue, Castellanos et al. (2011) use a marker‐based method to infer heritability for floral and vegetative traits in two Aquilegia species. Layered on top of this are estimates of the strength of natural selection. This novel joint estimation of heritability and selection in the wild showed that vegetative traits, compared to floral traits, have the highest evolutionarily potential. Evolutionary potential is the most important quantity to measure in wild populations. It combines inheritance and strength of selection and predicts the potential for populations to adapt to changing environments. The combination of molecular techniques with species in natural environments makes this work a model for molecular ecological investigations.     

Butterfly genomics sheds light on the process of hybrid speciation

How common is hybridization between species and what effect does it have on the evolutionary process? Can hybridization generate new species and what indeed is a species? In this issue, Gompert et al. (2014) show how massive, genome‐scale data sets can be used to shed light on these questions. They focus on the Lycaeides butterflies, and in particular, several populations from the western USA, which have characteristics suggesting that they may contain hybrids of two or more different species (Gompert et al. 2006). They demonstrate that these populations do contain mosaic genomes made up of components from different parental species. However, this appears to have been largely driven by historical admixture, with more recent processes appearing to be isolating the populations from each other. Therefore, these populations are on their way to becoming distinct species (if they are not already) but have arisen following extensive hybridization between other distinct populations or species (Fig.  1 ).

Figure 1 Open in figure viewer PowerPoint There has been extensive historical admixture between Lycaeides species with some new species arising from hybrid populations. (Photo credits: Lauren Lucas, Chris Nice, and James Fordyce).

Their data set must be one of the largest outside of humans, with over one and half thousand butterflies genotyped at over 45 thousand variable nucleotide positions. It is this sheer amount of data that has allowed the researchers to distinguish between historical and more recent evolutionary and demographic processes. This is because it has allowed them to partition the data into common and rare genetic variants and perform separate analyses on these. Common genetic variants are likely to be older while rare variants are more likely to be due to recent mutations. Therefore, by splitting the genetic variation into these components, the researchers were able to show more admixture among common variants, while rare variants showed less admixture and clear separation of the populations. The extensive geographic sampling of individuals, including overlapping distributions of several of the putative species, also allowed the authors to rule out the possibility that the separation of the populations was simply due to geographical distance. The authors have developed a new programme for detecting population structure and admixture, which does the same job as STRUCTURE (Pritchard et al. 2000 ), identifying genetically distinct populations and admixture between these populations, but is designed to be used with next generation sequence data. They use the output of this model for another promising new method to distinguish between contemporary and historical admixtures. They fixed the number of source populations in the model at two and estimated the proportion of each individual's genome coming from these two populations. Therefore, an individual can either be purely population 1, or population 2 or some mixture of the two (they call this value q, the same parameter exists in STRUCTURE). They then compared this to the level of heterozygosity coming from the two source populations in the individual's genome. If an individual is an F1 hybrid of two source populations, then it would have a q of 0.5 and also be heterozygous at all loci that distinguish the parental populations. On the other hand, if it is a member of a stable hybrid lineage, it might also have a q of 0.5 but would not be expected to be heterozygous at these loci, because over time the population would become fixed for one or other of the source population states either by drift or selection (Fig.  2 ). This is indeed what they find in the hybrid populations. They tend to have intermediate q values, but the level of heterozygosity coming from the source populations (which they call Q₁₂) was consistently lower than expected.

Figure 2 Open in figure viewer PowerPoint The Q‐matrix analysis used by Gompert et al. ( 2014 ) to distinguish between contemporary (hybrid swarm) and historical (stable hybrid lineage) admixture.

Overall, the results support several of the populations as being stable hybrid lineages. Nevertheless, the strictest definitions of hybrid species specify that the process of hybridization between the parental species must be instrumental in driving the reproductive isolation of the new species from both parental populations (Abbott et al. 2013 ). This is extremely hard to demonstrate conclusively because it requires us to first of all identify the isolating mechanisms that operated in the early evolution of the species and then to show that these were caused by the hybridization event itself. One advantage of the Lycaeides system is that the species appear to be in the early stages of divergence, so barriers to gene flow that are operating currently are likely to be those that are driving the species divergence. While there is some evidence that hybridization gave rise to traits that allowed the new populations to colonize new environments (Gompert et al. 2006 ; Lucas et al. 2008 ), there is clearly further work to be carried out in this direction. One of the rare examples of homoploid hybrid speciation (hybrid speciation without a change in chromosome number) where the reproductive isolation criterion has been demonstrated, comes from the Heliconius butterflies. In this case, hybridization of two species has been shown to give rise to a new colour pattern that instantly becomes reproductively isolated from the parental species due to mate preference for that pattern (Mavárez et al. 2006 ). However, while this has become a widely accepted example (Abbott et al. 2013 ), the naturally occurring ‘hybrid species’ in fact has derived most of its genome from one of the parental species, with largely just the colour pattern controlling locus coming from the other parent, a process that has been termed ‘hybrid trait speciation’ (Salazar et al. 2010 ). This distinction is an important one in terms of our understanding of the organization of biological diversity. While hybrid trait speciation will still largely fit the model of a neatly branching evolutionary tree, with perhaps only the region surrounding the single introgressed gene deviating from this model, hybrid species that end up with mosaic genomes, like Lycaeides, will not fit this model when considering the genome as a whole. This distinction also more broadly applies when comparing the patterns of divergence between Heliconius and Lycaeides. These two butterfly genera have been driving forward our understanding of the prevalence and importance of hybridization at the genomic level, but they reveal different ways in which hybridization can influence the organization of biological diversity. Recent work in Heliconius has shown that admixture is extensive and has been ongoing over a large portion of the evolutionary history of species (Martin et al. 2013 ; Nadeau et al. 2013 ). Nevertheless, this has not obscured the clear and robust pattern of a bifurcating evolutionary tree when considering the genome as a whole (Nadeau et al. 2013 ). In contrast in Lycaeides, the genome‐wide phylogeny clearly does not fit a bifurcating tree, resembling more of a messy shrub, with hybrid taxa falling at intermediate positions on the phylogeny (Gompert et al. 2014 ). The extent to which we need to rethink the way we describe and organize biological diversity will depend on the relative prevalence of these different outcomes of hybridization. We are likely to see many more of these types of large sequence data sets for ecologically interesting organisms. Gompert et al. ( 2014 ) show that these data need not only be a quantitative advance, but can also qualitatively change our understanding of the evolutionary history of these organisms. In particular, analysing common and rare genetic variants separately may provide information that would otherwise be missed. The emerging field of ‘speciation genomics’ (Seehausen et al. 2014 ) should follow this lead in developing new ways of making the most of the flood of genomic data that is being generated, but also improve methods for integrating this with field observations and experiments to identify the sources and targets of selection and divergence.

References

• Abbott R , Albach D , Ansell S et al. (2013 ) Hybridization and speciation . Journal of Evolutionary Biology, 26 , 229 – 246 . Wiley Online Library CAS PubMed Web of Science® Google Scholar
• Gompert Z , Fordyce JA , Forister ML , Shapiro AM , Nice CC (2006 ) Homoploid hybrid speciation in an extreme habitat . Science, 314 , 1923 – 1925 . Crossref CAS PubMed Web of Science® Google Scholar
• Gompert Z , Lucas LK , Buerkle CA et al. (2014 ) Admixture and the organization of genetic diversity in a butterfly species complex revealed through common and rare genetic variants . Molecular Ecology, 23 , 4555 – 4573 . Wiley Online Library CAS PubMed Web of Science® Google Scholar
• Lucas LK , Fordyce JA , Nice CC (2008 ) Patterns of genitalic morphology around suture zones in North American Lycaeides (Lepidoptera: Lycaenidae): implications for taxonomy and historical biogeography . Annals of the Entomological Society of America, 101 , 172 – 180 . Crossref Web of Science® Google Scholar
• Martin SH , Dasmahapatra KK , Nadeau NJ et al. (2013 ) Genome‐wide evidence for speciation with gene flow in Heliconius butterflies . Genome Research, 23 , 1817 – 1828 . Crossref CAS PubMed Web of Science® Google Scholar
• Mavárez J , Salazar CA , Bermingham E et al. (2006 ) Speciation by hybridization in Heliconius butterflies . Nature, 441 , 868 – 871 . Crossref CAS PubMed Web of Science® Google Scholar
• Nadeau NJ , Martin SH , Kozak KM et al. (2013 ) Genome‐wide patterns of divergence and gene flow across a butterfly radiation . Molecular Ecology, 22 , 814 – 826 . Wiley Online Library CAS PubMed Web of Science® Google Scholar
• Pritchard JK , Stephens M , Donnelly P (2000 ) Inference of population structure using multilocus genotype data . Genetics, 155 , 945 – 959 . Wiley Online Library CAS PubMed Web of Science® Google Scholar
• Salazar C , Baxter SW , Pardo‐Diaz C et al. (2010 ) Genetic evidence for hybrid trait speciation in Heliconius butterflies . PLoS Genetics, 6 , e1000930 . Crossref CAS PubMed Web of Science® Google Scholar
• Seehausen O , Butlin RK , Keller I et al. (2014 ) Genomics and the origin of species . Nature Reviews Genetics, 15 , 176 – 192 . Crossref CAS PubMed Web of Science® Google Scholar

This article was written and figures prepared by N.N. except as specified in the text (photo credits).

Citing Literature

Number of times cited according to CrossRef: 4

• V. Alex Sotola, David S. Ruppel, Timothy H. Bonner, Chris C. Nice, Noland H. Martin, Asymmetric introgression between fishes in the Red River basin of Texas is associated with variation in water quality, Ecology and Evolution, 10.1002/ece3.4901, 9 , 4, (2083-2095), (2019). Wiley Online Library
• Matej Bocek, Dominik Kusy, Michal Motyka, Ladislav Bocak, Persistence of multiple patterns and intraspecific polymorphism in multi-species Müllerian communities of net-winged beetles, Frontiers in Zoology, 10.1186/s12983-019-0335-8, 16 , 1, (2019). Crossref
• Nicola J. Nadeau, Takeshi Kawakami, Population Genomics of Speciation and Admixture, , 10.1007/13836_2018_24, (2018). Crossref
• Amanda Roe, Julian Dupuis, Felix Sperling, Molecular Dimensions of Insect Taxonomy in the Genomics Era, Insect Biodiversity, 10.1002/9781118945568, (547-573), (2017). Wiley Online Library

Volume 23 , Issue 18 September 2014

Pages 4441-4443     

Advances in ecological genomics in forest trees and applications to genetic resources conservation and breeding

Forest trees are an unparalleled group of organisms in their combined ecological, economic and societal importance. With widespread distributions, predominantly random mating systems and large population sizes, most tree species harbour extensive genetic variation both within and among populations. At the same time, demographic processes associated with Pleistocene climate oscillations and land‐use change have affected contemporary range‐wide diversity and may impinge on the potential for future adaptation. Understanding how these adaptive and neutral processes have shaped the genomes of trees species is therefore central to their management and conservation. As for many other taxa, the advent of high‐throughput sequencing methods is expected to yield an understanding of the interplay between the genome and environment at a level of detail and depth not possible only a few years ago. An international conference entitled ‘Genomics and Forest Tree Genetics’ was held in May 2016, in Arcachon (France), and brought together forest geneticists with a wide range of research interests to disseminate recent efforts that leverage contemporary genomic tools to probe the population, quantitative and evolutionary genomics of trees. An important goal of the conference was to discuss how such data can be applied to both genome‐enabled breeding and the conservation of forest genetic resources under land use and climate change. Here, we report discoveries presented at the meeting and discuss how the ecological genomic toolkit can be used to address both basic and applied questions in tree biology.     

A population genetic assessment of taxonomic species: The case of Lake Malawi cichlid fishes

Organisms sampled for population‐level research are typically assigned to species by morphological criteria. However, if those criteria are limited to one sex or life stage, or the organisms come from a complex of closely related forms, the species assignments may misdirect analyses. The impact of such sampling can be assessed from the correspondence of genetic clusters, identified only from patterns of genetic variation, to the species identified using only phenotypic criteria. We undertook this protocol with the rock‐dwelling mbuna cichlids of Lake Malawi, for which species within genera are usually identified using adult male coloration patterns. Given high local endemism of male colour patterns, and considerable allele sharing among species, there persists considerable taxonomic uncertainty in these fishes. Over 700 individuals from a single transect were photographed, genotyped and separately assigned: (a) to morphospecies using photographs; and (b) to genetic clusters using five widely used methods. Overall, the correspondence between clustering methods was strong for larger clusters, but methods varied widely in estimated number of clusters. The correspondence between morphospecies and genetic clusters was also strong for larger clusters, as well as some smaller clusters for some methods. These analyses generally affirm (a) adult male‐limited sampling and (b) the taxonomic status of Lake Malawi mbuna, as the species in our study largely appear to be well‐demarcated genetic entities. More generally, our analyses highlight the challenges for clustering methods when the number of populations is unknown, especially in cases of highly uneven sample sizes.     

A new multiplex SNP genotyping assay for detecting hybridization and introgression between the M and S molecular forms of Anopheles gambiae

The M and S forms of Anopheles gambiae have been the subject of intense study, but are morphologically indistinguishable and can only be identified using molecular techniques. PCR‐based assays to distinguish the two forms have been designed and applied widely. However, the application of these assays towards identifying hybrids between the two forms, and backcrossed hybrids in particular, has been problematic as the currently available diagnostic assays are based on single locus and/or are located within a multicopy gene. Here, we present an alternative genotyping method for detecting hybridization and introgression between M and S molecular forms based on a multilocus panel of single‐nucleotide polymorphisms (SNPs) fixed between the M and S forms. The panel of SNPs employed is located in so‐called islands of divergence leading us to describe this method as the ‘Divergence Island SNP’ (DIS) assay. We show this multilocus SNP genotyping approach can robustly and accurately detect F1 hybrids as well as backcrossed individuals.     

Female mate choice of male signals is unlikely to promote ecological adaptation in Enchenopa treehoppers (Hemiptera: Membracidae)

A key question in speciation research is how ecological and sexual divergence arise and interact. We tested the hypothesis that mate choice causes local adaptation and ecological divergence using the rationale that the performance~signal trait relationship should parallel the attractiveness~signal trait relationship. We used female fecundity as a measure of ecological performance. We used a species in the Enchenopa binotata treehopper complex, wherein speciation involves adaptation to novel environments and divergence in sexual communication. We used a full‐sibling, split‐family rearing design to estimate genetic correlations (r_G) between fecundity and signal traits, and compared those relationships against population‐level mate preferences for the signal traits. Animal model estimates for r_G between female fecundity and male signal traits overlapped zero—rejecting the hypothesis—but could reflect sample size limitations. The magnitude of r_G correlated with the strength of the mate preferences for the corresponding signal traits, especially for signal frequency, which has the strongest mate preference and the most divergence in the complex. However, signal frequencies favored by the population‐level mate preference are not associated with high fecundity. Therefore, mate preferences do not appear to have been selected to favor high‐performance genotypes. Our findings suggest that ecological and sexual divergence may arise separately, but reinforce each other, during speciation.     

The genomic and ecological context of hybridization affects the probability that symmetrical incompatibilities drive hybrid speciation

Despite examples of homoploid hybrid species, theoretical work describing when, where, and how we expect homoploid hybrid speciation to occur remains relatively rare. Here, I explore the probability of homoploid hybrid speciation due to “symmetrical incompatibilities” under different selective and genetic scenarios. Through simulation, I test how genetic architecture and selection acting on traits that do not themselves generate incompatibilities interact to affect the probability that hybrids evolve symmetrical incompatibilities with their parent species. Unsurprisingly, selection against admixture at “adaptive” loci that are linked to loci that generate incompatibilities tends to reduce the probability of evolving symmetrical incompatibilities. By contrast, selection that favors admixed genotypes at adaptive loci can promote the evolution of symmetrical incompatibilities. The magnitude of these outcomes is affected by the strength of selection, aspects of genetic architecture such as linkage relationships and the linear arrangement of loci along a chromosome, and the amount of hybridization following the formation of a hybrid zone. These results highlight how understanding the nature of selection, aspects of the genetics of traits affecting fitness, and the strength of reproductive isolation between hybridizing taxa can all be used to inform when we expect to observe homoploid hybrid speciation due to symmetrical incompatibilities.     

The role of structural genomic variants in population differentiation and ecotype formation in Timema cristinae walking sticks

Theory predicts that structural genomic variants such as inversions can promote adaptive diversification and speciation. Despite increasing empirical evidence that adaptive divergence can be triggered by one or a few large inversions, the degree to which widespread genomic regions under divergent selection are associated with structural variants remains unclear. Here we test for an association between structural variants and genomic regions that underlie parallel host‐plant‐associated ecotype formation in Timema cristinae stick insects. Using mate‐pair resequencing of 20 new whole genomes we find that moderately sized structural variants such as inversions, deletions and duplications are widespread across the genome, being retained as standing variation within and among populations. Using 160 previously published, standard‐orientation whole genome sequences we find little to no evidence that the DNA sequences within inversions exhibit accentuated differentiation between ecotypes. In contrast, a formerly described large region of reduced recombination that harbours genes controlling colour‐pattern exhibits evidence for accentuated differentiation between ecotypes, which is consistent with differences in the frequency of colour‐pattern morphs between host‐associated ecotypes. Our results suggest that some types of structural variants (e.g., large inversions) are more likely to underlie adaptive divergence than others, and that structural variants are not required for subtle yet genome‐wide genetic differentiation with gene flow.     

Comparing geographical genetic differentiation between candidate and noncandidate loci for adaptation strengthens support for parallel ecological divergence in the marine snail Littorina saxatilis

The Galician sympatric ecotypes of Littorina saxatilis have been proposed as a model system for studying parallel ecological speciation. Such a model system makes a clear prediction: candidate loci (for divergent adaptation) should present a higher level of geographical differentiation than noncandidate (neutral) loci. We used 2356 amplified fragment length polymorphisms (AFLPs) and four microsatellite loci to identify candidate loci for ecological adaptation using the F _ST outlier method. Three per cent of the studied AFLP loci were identified as candidate loci associated with adaptation, after multitest adjustments, thus contributing to ecotype differentiation (candidate loci were not detected within ecotypes). Candidate and noncandidate loci were analysed separately at four different F _ST partitions: differences between ecotypes (overall and local), differences between localities and micro-geographical differences within ecotypes. The magnitude of F _ST differed between candidate and noncandidate loci for all partitions except in the case of microgeographical differentiation within ecotypes, and the microsatellites (putatively neutral) showed an identical pattern to noncandidate loci. Thus, variation in candidate loci is determined partially independent by divergent natural selection (in addition to stochastic forces) at each locality, while noncandidate loci are exclusively driven by stochastic forces. These results support the evolutionary history described for these particular populations, considered to be a clear example of incomplete sympatric ecological speciation.     

Evaluation of genetic isolation within an island flora reveals unusually widespread local adaptation and supports sympatric speciation

It is now recognized that speciation can proceed even when divergent natural selection is opposed by gene flow. Understanding the extent to which environmental gradients and geographical distance can limit gene flow within species can shed light on the relative roles of selection and dispersal limitation during the early stages of population divergence and speciation. On the remote Lord Howe Island (Australia), ecological speciation with gene flow is thought to have taken place in several plant genera. The aim of this study was to establish the contributions of isolation by environment (IBE) and isolation by community (IBC) to the genetic structure of 19 plant species, from a number of distantly related families, which have been subjected to similar environmental pressures over comparable time scales. We applied an individual-based, multivariate, model averaging approach to quantify IBE and IBC, while controlling for isolation by distance (IBD). Our analyses demonstrated that all species experienced some degree of ecologically driven isolation, whereas only 12 of 19 species were subjected to IBD. The prevalence of IBE within these plant species indicates that divergent selection in plants frequently produces local adaptation and supports hypotheses that ecological divergence can drive speciation in sympatry.     

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.     

Hybridization at an ecotone: ecological and genetic barriers between three Iberian vipers

The formation of stable genetic boundaries between emerging species is often diagnosed by reduced hybrid fitness relative to parental taxa. This reduced fitness can arise from endogenous and/or exogenous barriers to gene flow. Although detecting exogenous barriers in nature is difficult, we can estimate the role of ecological divergence in driving species boundaries by integrating molecular and ecological niche modelling tools. Here, we focus on a three‐way secondary contact zone between three viper species (Vipera aspis, V. latastei and V. seoanei) to test for the contribution of ecological divergence to the development of reproductive barriers at several species traits (morphology, nuclear DNA and mitochondrial DNA). Both the nuclear and mitochondrial data show that all taxa are genetically distinct and that the sister species V. aspis and V. latastei hybridize frequently and backcross over several generations. We find that the three taxa have diverged ecologically and meet at a hybrid zone coincident with a steep ecotone between the Atlantic and Mediterranean biogeographical provinces. Integrating landscape and genetic approaches, we show that hybridization is spatially restricted to habitats that are suboptimal for parental taxa. Together, these results suggest that niche separation and adaptation to an ecological gradient confer an important barrier to gene flow among taxa that have not achieved complete reproductive isolation.