Two are better than one: combining landscape genomics and common gardens for detecting local adaptation in forest trees

Predicting likely species responses to an alteration of their local environment is key to decision‐making in resource management, ecosystem restoration and biodiversity conservation practice in the face of global human‐induced habitat disturbance. This is especially true for forest trees which are a dominant life form on Earth and play a central role in supporting diverse communities and structuring a wide range of ecosystems. In Europe, it is expected that most forest tree species will not be able to migrate North fast enough to follow the estimated temperature isocline shift given current predictions for rapid climate warming. In this context, a topical question for forest genetics research is to quantify the ability for tree species to adapt locally to strongly altered environmental conditions (Kremer et al. 2012 ). Identifying environmental factors driving local adaptation is, however, a major challenge for evolutionary biology and ecology in general but is particularly difficult in trees given their large individual and population size and long generation time. Empirical evaluation of local adaptation in trees has traditionally relied on fastidious long‐term common garden experiments (provenance trials) now supplemented by reference genome sequence analysis for a handful of economically valuable species. However, such resources have been lacking for most tree species despite their ecological importance in supporting whole ecosystems. In this issue of Molecular Ecology, De Kort et al. ( 2014 ) provide original and convincing empirical evidence of local adaptation to temperature in black alder, Alnus glutinosa L. Gaertn, a surprisingly understudied keystone species supporting riparian ecosystems. Here, De Kort et al. ( 2014 ) use an innovative empirical approach complementing state‐of‐the‐art landscape genomics analysis of A. glutinosa populations sampled in natura across a regional climate gradient with phenotypic trait assessment in a common garden experiment (Fig. 1 ). By combining the two methods, De Kort et al. ( 2014 ) were able to detect unequivocal association between temperature and phenotypic traits such as leaf size as well as with genetic loci putatively under divergent selection for temperature. The research by De Kort et al. ( 2014 ) provides valuable insight into adaptive response to temperature variation for an ecologically important species and demonstrates the usefulness of an integrated approach for empirical evaluation of local adaptation in nonmodel species (Sork et al. 2013 ).     

The genetic dimension of pest pressure in the tropical rainforest

Wet tropical forests are among the most diverse ecosystems on Earth and can host several hundreds of tree species per hectare. To maintain such diversity, the community must contain large numbers of relatively rare species rather than be dominated by a few very common trees, as is often the case in temperate forests. Explaining the mechanisms preventing dominance by common species has been a major task of tropical forest ecology. One of the most promising mechanisms is negative density dependence (NDD) of tree abundance driven by pests, including fungal diseases (‘pest pressure’). NDD entails that the chance of survival of a sapling increases with the distance from a mature tree of the same species, thus preventing species from becoming locally dominant. Curiously, the strength of NDD is negatively correlated with abundance, meaning that tree species that are more common generally show weaker NDD (Comita et al. 2010 ). Interactions between plants and soil pathogens have been shown to play an important role in NDD (Klironomos 2002 ), and rare species are apparently more strongly affected (Mangan et al. 2010 ). However, the genetic mechanisms underlying this phenomenon have remained obscure. In this issue of Molecular Ecology, Marden et al. ( 2017 ) suggest that reduced diversity of the genes involved in pathogen recognition (Resistance genes or R genes) could explain why NDD is stronger in locally rare species.     

Homomorphic plant sex chromosomes are coming of age

Sex chromosomes are a very peculiar part of the genome that have evolved independently in many groups of animals and plants (Bull 1983 ). Major research efforts have so far been focused on large heteromorphic sex chromosomes in a few animal and plant species (Chibalina & Filatov 2011 ; Zhou & Bachtrog 2012 ; Bellott et al. 2014 ; Hough et al. 2014 ; Zhou et al. 2014 ), while homomorphic (cytologically indistinguishable) sex chromosomes have largely been neglected. However, this situation is starting to change. In this issue, Geraldes et al. ( 2015 ) describe a small (~100 kb long) sex‐determining region on the homomorphic sex chromosomes of poplars (Populus trichocarpa and related species, Fig.  1 ). All species in Populus and its sister genus Salix are dioecious, suggesting that dioecy and the sex chromosomes, if any, should be relatively old. Contrary to this expectation, Geraldes et al. ( 2015 ) demonstrate that the sex‐determining region in poplars is of very recent origin and probably evolved within the genus Populus only a few million years ago.     

Mock communities highlight the diversity of host‐associated eukaryotes

Host‐associated microbes are ubiquitous. Every multicellular eukaryote, and even many unicellular eukaryotes (protists), hosts a diverse community of microbes. High‐throughput sequencing (HTS) tools have illuminated the vast diversity of host‐associated microbes and shown that they have widespread influence on host biology, ecology and evolution (McFall‐Ngai et al. 2013 ). Bacteria receive most of the attention, but protists are also important components of microbial communities associated with humans (Parfrey et al. 2011 ) and other hosts. As HTS tools are increasingly used to study eukaryotes, the presence of numerous and diverse host‐associated eukaryotes is emerging as a common theme across ecosystems. Indeed, HTS studies demonstrate that host‐associated lineages account for between 2 and 12% of overall eukaryotic sequences detected in soil, marine and freshwater data sets, with much higher relative abundances observed in some samples (Ramirez et al. 2014 ; Simon et al. 2015 ; de Vargas et al. 2015 ). Previous studies in soil detected large numbers of predominantly parasitic lineages such as Apicomplexa, but did not delve into their origin [e.g. (Ramirez et al. 2014 )]. In this issue of Molecular Ecology, Geisen et al. ( 2015 ) use mock communities to show that many of the eukaryotic organisms detected by environmental sequencing in soils are potentially associated with animal hosts rather than free‐living. By isolating the host‐associated fraction of soil microbial communities, Geisen and colleagues help explain the surprisingly high diversity of parasitic eukaryotic lineages often detected in soil/terrestrial studies using high‐throughput sequencing (HTS) and reinforce the ubiquity of these host‐associated microbes. It is clear that we can no longer assume that organisms detected in bulk environmental sequencing are free‐living, but instead need to design studies that specifically enumerate the diversity and function of host‐associated eukaryotes. Doing so will allow the field to determine the role host‐associated eukaryotes play in soils and other environments and to evaluate hypotheses on assembly of host‐associated communities, disease ecology and more.     

Host plant richness explains diversity of ectomycorrhizal fungi: Response to the comment of Tedersoo et al. (2014)

Exploring the relationships between the biodiversity of groups of interacting organisms yields insight into ecosystem stability and function (Hooper et al. 2000 ; Wardle 2006 ). We demonstrated positive relationships between host plant richness and ectomycorrhizal (EM) fungal diversity both in a field study in subtropical China (Gutianshan) and in a meta‐analysis of temperate and tropical studies (Gao et al. 2013 ). However, based on re‐evaluation of our data sets, Tedersoo et al. ( 2014 ) argue that the observed positive correlation between EM fungal richness and EM plant richness at Gutianshan and also in our metastudies was based mainly from (i) a sampling design with inconsistent species pool and (ii) poor data compilation for the meta‐analysis. Accordingly, we checked our data sets and repeated the analysis performed by Tedersoo et al. ( 2014 ). In contrast to Tedersoo et al. ( 2014 ), our re‐analysis still confirms a positive effect of plant richness on EM fungal diversity in Gutianshan, temperate and tropical ecosystems, respectively.     

Using RAD‐seq to recognize sex‐specific markers and sex chromosome systems

Next‐generation sequencing methods have initiated a revolution in molecular ecology and evolution (Tautz et al. 2010 ). Among the most impressive of these sequencing innovations is restriction site‐associated DNA sequencing or RAD‐seq (Baird et al. 2008 ; Andrews et al. 2016 ). RAD‐seq uses the Illumina sequencing platform to sequence fragments of DNA cut by a specific restriction enzyme and can generate tens of thousands of molecular genetic markers for analysis. One of the many uses of RAD‐seq data has been to identify sex‐specific genetic markers, markers found in one sex but not the other (Baxter et al. 2011 ; Gamble & Zarkower 2014 ). Sex‐specific markers are a powerful tool for biologists. At their most basic, they can be used to identify the sex of an individual via PCR. This is useful in cases where a species lacks obvious sexual dimorphism at some or all life history stages. For example, such tests have been important for studying sex differences in life history (Sheldon 1998 ; Mossman & Waser 1999 ), the management and breeding of endangered species (Taberlet et al. 1993 ; Griffiths & Tiwari 1995 ; Robertson et al. 2006 ) and sexing embryonic material (Hacker et al. 1995 ; Smith et al. 1999 ). Furthermore, sex‐specific markers allow recognition of the sex chromosome system in cases where standard cytogenetic methods fail (Charlesworth & Mank 2010 ; Gamble & Zarkower 2014 ). Thus, species with male‐specific markers have male heterogamety (XY) while species with female‐specific markers have female heterogamety (ZW). In this issue, Fowler & Buonaccorsi ( 2016 ) illustrate the ease by which RAD‐seq data can generate sex‐specific genetic markers in rockfish (Sebastes). Moreover, by examining RAD‐seq data from two closely related rockfish species, Sebastes chrysomelas and Sebastes carnatus (Fig.  1 ), Fowler & Buonaccorsi ( 2016 ) uncover shared sex‐specific markers and a conserved sex chromosome system.     

The importance of genomic novelty in social evolution

Insect societies dominate the natural world: They mould landscapes, sculpt habitats, pollinate plants, sow seeds and control pests. The secret to their success lies in the evolution of queen (reproductive) and worker (provisioner and carer) castes (Oster & Wilson 1978 ). A major problem in evolutionary biology is explaining the evolution of insect castes, particularly the workers (Darwin 1859 ). Next‐generation sequencing technologies now make it possible to understand how genomic material is born, lost and reorganized in the evolution of alternative phenotypes. Such analyses are revealing a general role for novel (e.g. taxonomically restricted) genes in phenotypic innovations across the animal kingdom (Chen et al. 2013). In this issue of molecular ecology, Feldmeyer et al. (2014) provide overwhelming evidence for the importance of novel genes in caste evolution in an ant. Feldmeyer et al.'s study is important and exciting because it cements the role of genomic novelty, as well as conservation, firmly into the molecular jigsaw of social evolution. Evolution is eclectic in its exploitation of both old and new genomic material to generate replicated phenotypic innovations across the tree of life.     

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.     

Tree diversity in relation to maximum tree height: evidence for the harshness hypothesis of species diversity gradients

Marks et al. (Ecology Letters, 19, 2016, 743) showed tree species richness correlates with maximum tree height, and interpret this as evidence that the environmental stressors that limit tree height also act as ecological filters on species richness. Here, we strengthen these arguments by further addressing the roles of environmental covariates and beta diversity.     

Opportunities and challenges of next‐generation sequencing applications in ecological epigenetics

Evolutionary theory posits that adaptation can result when populations harbour heritable phenotypic variation for traits that increase tolerance to local conditions. However, the actual mechanisms that underlie heritable phenotypic variation are not completely understood (Keller 2014 ). Recently, the potential role of epigenetic mechanisms in the process of adaptive evolution has been the subject of much debate (Pigliucci & Finkelman 2014 ). Studies of variation in DNA methylation in particular have shown that natural populations harbour high amounts of epigenetic variation, which can be inherited across generations and can cause heritable trait variation independently of genetic variation (Kilvitis et al. 2014 ). While we have made some progress addressing the importance of epigenetics in ecology and evolution using methylation‐sensitive AFLP (MS‐AFLP), this approach provides relatively few anonymous and dominant markers per individual. MS‐AFLP are difficult to link to functional genomic elements or phenotype and are difficult to compare directly to genetic variation, which has limited the insights drawn from studies of epigenetic variation in natural nonmodel populations (Schrey et al. 2013 ). In this issue, Platt et al. provide an example of a promising approach to address this problem by applying a reduced representation bisulphite sequencing (RRBS) approach based on next‐generation sequencing methods in an ecological context.     

Why some parasites are widespread and abundant while others are local and rare?

Abundances and distributions of species are usually associated. This implies that as a species declines in abundance so does the number of sites it occupies. Conversely, when there is an increase in a species' range size, it is usually followed by an increase in population size (Gaston et al. 2000 ). This ecological phenomenon, also known as the abundance–occupancy relationship (AOR), is well documented in several species of animals and plants (Gaston et al. 2000 ) but has been little investigated in parasites. In this issue of Molecular Ecology, Drovetski et al. ( 2014 ) investigated the AOR in avian haemosporidians (vector‐borne blood parasites) using data from four well‐sampled bird communities. In support of the AOR, the research group found that the abundance of parasite cytochrome b lineages (a commonly used proxy for species identification within this group of parasites) was positively linked with the abundance of susceptible avian host species and that the most abundant haemospordian lineages were those with larger ranges. Drovetski et al. ( 2014 ) also found evidence for both hypotheses proposed to explain the AOR in parasites: the trade‐off hypothesis (TOH) and the niche‐breadth hypothesis (NBH). Interestingly, the main predictor of the AOR was the number of susceptible hosts (i.e. number of infected birds) and not the number of host species the parasites were able to exploit.     

Demystifying the RAD fad

We are writing in response to the population and phylogenomics meeting review by Andrews & Luikart ( 2014 ) entitled ‘Recent novel approaches for population genomics data analysis’. Restriction‐site‐associated DNA (RAD) sequencing has become a powerful and useful approach in molecular ecology, with several different published methods now available to molecular ecologists, none of which can be considered the best option in all situations. A&L report that the original RAD protocol of Miller et al. ( 2007 ) and Baird et al. ( 2008 ) is superior to all other RAD variants because putative PCR duplicates can be identified (see Baxter et al. 2011 ), thereby reducing the impact of PCR artefacts on allele frequency estimates (Andrews & Luikart 2014 ). In response, we (i) challenge the assertion that the original RAD protocol minimizes the impact of PCR artefacts relative to that of other RAD protocols, (ii) present additional biases in RADseq that are at least as important as PCR artefacts in selecting a RAD protocol and (iii) highlight the strengths and weaknesses of four different approaches to RADseq which are a representative sample of all RAD variants: the original RAD protocol (mbRAD, Miller et al. 2007 ; Baird et al. 2008 ), double digest RAD (ddRAD, Peterson et al. 2012 ), ezRAD (Toonen et al. 2013 ) and 2bRAD (Wang et al. 2012 ). With an understanding of the strengths and weaknesses of different RAD protocols, researchers can make a more informed decision when selecting a RAD protocol.     

The genomics of adaptation,divergence and speciation: a congealing theory

In this issue, Flaxman et al. ( 2014 ) report the results of sophisticated whole‐genome simulations of speciation with gene flow, enhancing our understanding of the process by building on previous single‐locus, multilocus and analytical works. Their findings provide us with new insights about how genomes can diverge and the importance of statistical and chromosomal linkage in facilitating reproductive isolation. The authors characterize the conditions under which, even with high gene flow and weak divergent selection, reproductive isolation between populations can occur due to the emergent stochastic process of genomewide congealing, where numerous statistically or physically linked loci of small effect allow selection to limit effective migration rates. The initial congealing event can occur within a broad range conditions, and once initiated, the self‐reinforcing process leads to rapid divergence and ultimately two reproductively isolated populations. Flaxman et al.'s ( 2014 ) work is a valuable contribution to our understanding of speciation with gene flow and in making a more predictive field of evolutionary genomics and speciation.     

Marine island biogeography. Response to comment on ‘Island biogeography: patterns of marine shallow‐water organisms’

In this response we have incorporated data on gastropod and seaweed biodiversity referred to by Ávila et al. (2016, Journal of Biogeography, doi: 10.1111/jbi.12816 ) to allow an updated analysis on marine shallow‐water biogeography patterns. When compared to the biogeography patterns reported in Hachich et al. (2015, Journal of Biogeography, 42 , 1871–1882), we find (1) no differences in the patterns originally reported for reef fish or seaweeds, (2) minor differences in gastropod species–area and species–age patterns and (3) a significant difference for the gastropod species‐isolation pattern. In our original work, we reported that there was limited evidence that gastropod species richness was influenced by island isolation; however, our new analysis reveals a power‐model relationship between these variables. Thus, we are now able to conclude that gastropod species diversity, whose dispersal capacity is intermediate between seaweeds (lowest) and reef fish (highest), is also influenced by island isolation.     

Tending a complex microbiota requires major immune complexity

Animals maintain complex microbial communities within their guts that fill important roles in the health and development of the host. To what degree a host's genetic background influences the establishment and maintenance of its gut microbial communities is still an open question. We know from studies in mice and humans that external factors, such as diet and environmental sources of microbes, and host immune factors play an important role in shaping the microbial communities (Costello et al. 2012 ). In this issue of Molecular Ecology, Bolnick et al. ( 2014a ) sample the gut microbial community from 150 genetically diverse stickleback isolated from a single lake to provide evidence that another part of the adaptive immune response, the major histocompatibility complex class II (MHCII) receptors of antigen‐presenting cells, may play a role in shaping the gut microbiota of the threespine stickleback, Gasterosteus aculeatus (Bolnick et al. 2014a ). Bolnick et al. ( 2014a ) provide insight into natural, interindividual variation in the diversity of both stickleback MHCII alleles and their gut microbial communities and correlate changes in the diversity of MHCII receptor alleles with changes in the microbiota.     

A robust tool highlights the influence of bird migration on influenza A virus evolution

One of the fundamental unknowns in the field of influenza biology is a panoramic understanding of the role wild birds play in the global maintenance and spread of influenza A viruses. Wild aquatic birds are considered a reservoir host for all lowly pathogenic avian influenza A viruses (AIV) and thus serve as a potential source of zoonotic AIV, such as Australasian‐origin H5N1 responsible for morbidity and mortality in both poultry and humans, as well as genes that may contribute to the emergence of pandemic viruses. Years of broad, in‐depth wild bird AIV surveillance have helped to decipher key observations and ideas regarding AIV evolution and viral ecology including the trending of viral lineages, patterns of gene flow within and between migratory flyways and the role of geographic boundaries in shaping viral evolution (Bahl et al. 2009 ; Lam et al. 2012 ). While these generally ‘virus‐centric’ studies have ultimately advanced our broader understanding of AIV dynamics, recent studies have been more host‐focused, directed at determining the potential impact of host behaviour on AIV, specifically, the influence of bird migration upon AIV maintenance and transmission. A large number of surveillance studies have taken place in Alaska, United States—a region where several global flyways overlap—with the aim of detecting the introduction of novel, Australasian‐origin highly pathogenic H5N1 AIV into North America. By targeting bird species with known migration habits, long‐distance migrators were determined to be involved in the intercontinental movement of individual AIV gene segments, but not entire viruses, between the Australasian and North American flyways (Koehler et al. 2008 ; Pearce et al. 2010 ). Yet, bird movement is not solely limited to long‐distance migration, and the relationship of resident or nonmigratory and intermediate‐distance migrant populations with AIV ecology has only recently been explored by Hill et al. ( 2012 ) in this issue of Molecular Ecology. Applying a uniquely refined, multidimensional approach, Hill et al. validate the innovative use of stable isotope assays for qualifying migration status of wild mallards within the Pacific flyway. The authors reveal that AIV prevalence and diversity did not differ in wintering mallard ducks with different migration strategies, and while migrant mallards do indeed introduce AIV, these viruses do not circulate as the predominant viruses in resident birds. On the other hand, resident mallards from more temperate regions act as reservoirs, possibly contributing to the unseasonal circulation and extended transmission period of AIV. This study highlights the impact of animal behaviour on shaping viral evolution, and the unique observations made will help inform prospective AIV surveillance efforts in wild birds.     

Geographical genetics and the conservation of forest trees

Trees are key ecosystem engineers. Many analyses of the genetic diversity of forest trees over substantial parts of their distributional ranges are appearing. These studies are of relevance for forest and landscape management, the inventory of botanical genetic resources and the conservation biology of rare, endemic, relictual, and endangered tree species. This review focuses on (i) recent investigations of the influence of human disturbances, (ii) comparative analyses of closely related and hybridizing species, (iii) reconstructions of refugia and of the spread of tree populations during the postglacial, (iv) studies of both range-wide and range-edge genetic patterns, and (v) assessments of the role of tree genetic diversity in the face of future climate warming. There is a need to include more tropical and austral trees in genetic analyses, as most studies have dealt with the relatively species-poor Palaearctic and Nearctic regions. Further studies are also needed on the role of tree genetic diversity in variations in phenology, resistance to insect defoliators and fungal pathogens, reactions to increased CO₂ and ozone concentrations, growth, mortality rates and other traits. Most macroecological and scaling patterns of species richness still need to be studied for genetic diversity. Open research questions in this rapidly evolving field involve invasion biology, island biogeography, and urban ecology. There is a need for more knowledge transfer from the many studies of tree genetic diversity to the day-to-day management of trees and forests.     

Terra incognita of soil biodiversity: unseen invasions under our feet

Whilst cartographers of the 19th century endeavoured to chart the last unknown lands, the great challenge for biologists in the 21st century is to fill the gaps on the biodiversity map of the Earth. And one of the largest gaps concerns the biodiversity of soils, a terra incognita right under our feet. The study of soil biodiversity, and particularly the complex communities of small invertebrates, has suffered from a severe ‘taxonomic impediment’ (Decaëns 2010 ) leading to great uncertainties about total species richness, phylogenetic diversity, geographical structure, temporal dynamics of soil organisms, and consequently about their role on ecosystem function (Bardgett & van der Putten 2014 ). However, the revolution in high‐throughput sequencing is now revealing the hidden biodiversity of the soil with unprecedented detail (e.g. Arribas et al. 2016 ). In a noteworthy from the Cover article in this issue of Molecular Ecology, Cicconardi et al. ( 2017 ) apply these new tools to study soil communities of Collembola in three distant oceanic islands of volcanic origin, obtaining a striking result: only 38 of 70 species (54%) are exclusively found in a single island, with the remaining shared among islands or with other distant regions, suggesting a massive recent introduction of soil species, whose impact is entirely unknown.     

A tree of tree frogs around the Black Sea

Speciation, the process by which one species evolves into two or more, is a major focus of ongoing debate, particularly regarding the geographic context in which it occurs. Geographic models of speciation tend to fall into discrete categories, typically referred to as allopatric, parapatric and sympatric speciation, according to whether two groups evolve reproductive isolation while geographically isolated, differentiated but connected by gene flow, or completely co‐occurring. Yet molecular studies indicate that full development of reproductive isolation can take very long compared with the timescale at which climatic oscillations occur, such that the geographic context of differentiating forms might change often during the long process to full species. Studies of genetic relationships across the ranges of organisms with low‐dispersal distances have the potential to reveal these complex histories. In a particularly elegant example in this issue, Dufresnes et al. ( 2016 ) use genetic variation and ecological niche modelling to show that a ring of populations of the eastern tree frog (Hyla orientalis) surrounding the Black Sea had a complex history of geographic differentiation. Alternating phases of geographic fragmentation and phases of gene flow between neighbouring populations have produced a pattern of gradual genetic change connecting the western, southern and eastern sides of the ring, with the northwestern and northeastern forms being most differentiated. In the north, a population in Crimea appears to have been produced through mixture of the two extreme forms. The overall genetic relationships are reminiscent of those found in ring species, which have been used as prime demonstrations of the process of speciation. The difference, however, is that the terminal forms appear to have mixed rather than be reproductively isolated, although more research is needed to infer whether there might be some reproductive isolation on the northern side of the ring.     

Coalescing molecular evolution and DNA barcoding

The DNA barcoding concept (Woese et al. 1990 ; Hebert et al. 2003 ) has considerably boosted taxonomy research by facilitating the identification of specimens and discovery of new species. Used alone or in combination with DNA metabarcoding on environmental samples (Taberlet et al. 2012 ), the approach is becoming a standard for basic and applied research in ecology, evolution and conservation across taxa, communities and ecosystems (Scheffers et al. 2012 ; Kress et al. 2015 ). However, DNA barcoding suffers from several shortcomings that still remain overlooked, especially when it comes to species delineation (Collins & Cruickshank 2012 ). In this issue of Molecular Ecology, Barley & Thomson ( 2016 ) demonstrate that the choice of models of sequence evolution has substantial impacts on inferred genetic distances, with a propensity of the widely used Kimura 2‐parameter model to lead to underestimated species richness. While DNA barcoding has been and will continue to be a powerful tool for specimen identification and preliminary taxonomic sorting, this work calls for a systematic assessment of substitution models fit on barcoding data used for species delineation and reopens the debate on the limitation of this approach.