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.     

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.     

Aquatic biodiversity assessment for the lazy

The world is covered in DNA. In any ecosystem, extracellular DNA fragments can be found that once formed the genomes of a variety of micro‐ and macroorganisms. A few years ago, it was proposed to use this environmental DNA (eDNA) as a source of information on local vertebrate biodiversity (Ficetola et al. 2008 ; Taberlet et al. 2012 ). This idea offered an elegant solution to take up the gauntlet of rapidly increasing monitoring needs. Coupled with barcoding efforts, it promised to be cost‐efficient in many respects, for example man‐hours and taxonomic expertise. Ecologists and conservation biologists with an interest in aquatic ecosystems have enthusiastically adopted and pioneered this new method, producing dozens of eDNA studies. Most of these studies have, however, focused on a single or a few aquatic species. In this issue of Molecular Ecology, Valentini et al. ( 2016 ) move the field a step further by demonstrating that metabarcoding approaches – which simultaneously target large groups of organisms such as amphibians or fish – can match and sometimes even outperform other inventory methods.     

Will reduced host connectivity curb the spread of a devastating epidemic?

The white‐nose syndrome (WNS), caused by the fungal pathogen Pseudogymnoascus destructans, is threatening the cave‐dwelling bat fauna of North America by killing individuals by the thousands in hibernacula each winter since its appearance in New York State less than ten years ago. Epidemiological models predict that WNS will reach the western coast of the USA by 2035, potentially eliminating most populations of susceptible bat species in its path (Frick et al. 2015; O'Regan et al. 2015). These models were built and validated using distributional data from the early years of the epidemic, which spread throughout eastern North America following a route driven by cave density and winter severity (Maher et al. 2012). In this issue of Molecular Ecology, Wilder et al. (2015) refine these findings by showing that connectivity among host populations, as assessed by population genetic markers, is crucial in determining the spread of the pathogen. Because host connectivity is much reduced in the hitherto disease free western half of North America, Wilder et al. make the reassuring prediction that the disease will spread more slowly west of the Great Plains.     

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.     

Standing and flowing: the complex origins of adaptive variation

A population faced with a new selection pressure can only adapt if appropriate genetic variation is available. This genetic variation might come from new mutations or from gene exchange with other populations or species, or it might already segregate in the population as standing genetic variation (which might itself have arisen from either mutation or gene flow). Understanding the relative importance of these sources of adaptive variation is a fundamental issue in evolutionary genetics (Orr & Betancourt 2001 ; Barrett & Schluter 2008 ; Gladyshev et al. 2008 ) and has practical implications for conservation, plant and animal breeding, biological control and infectious disease prevention (e.g. Robertson 1960 ; Soulé & Wilcox 1980 ; Prentis et al. 2008 ; Pennings 2012 ). In this issue of Molecular Ecology, Roesti et al. ( 2014 ) make an important contribution to this longstanding debate.     

The discovery of Antarctic RNA viruses: a new game changer

Antarctic ecosystems are dominated by micro‐organisms, and viruses play particularly important roles in the food webs. Since the first report in 2009 (López‐Bueno et al. 2009 ), ‘omic’‐based studies have greatly enlightened our understanding of Antarctic aquatic microbial diversity and ecosystem function (Wilkins et al. 2013 ; Cavicchioli 2015 ). This has included the discovery of many new eukaryotic viruses (López‐Bueno et al. 2009 ), virophage predators of algal viruses (Yau et al. 2011 ), bacteria with resistance to phage (Lauro et al. 2011 ) and mechanisms of haloarchaeal evasion, defence and adaptation to viruses (Tschitschko et al. 2015 ). In this issue of Molecular Ecology, López‐Bueno et al. ( 2015 ) report the first discovery of RNA viruses from an Antarctic aquatic environment. High sequence coverage enabled genome variation to be assessed for four positive‐sense single‐stranded RNA viruses from the order Picornavirales. By examining the populations present in the water column and in the lake's catchment area, populations of ‘quasispecies’ were able to be linked to local environmental factors. In view of the importance of viruses in Antarctic ecosystems but lack of data describing them, this study represents a significant advance in the field.     

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.     

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.     

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 ).     

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.     

Next‐generation phylogenetics takes root

It has been a tumultuous 5 years in phylogeography and phylogenetics during which both fields have struggled to harness the power of next‐generation sequencing (NGS) (Ekblom & Galindo 2010 ; McCormack et al. 2012a ). Fortunately, several methodological approaches appear to be taking root. In this issue of Molecular Ecology, O'Neill et al. 2013 ) employ one such method – parallel tagged sequencing (PTS) – to elucidate the phylogeography of a tiger salamander (Ambystoma tigrinum) species complex. This study demonstrates a practical application of NGS on a scale appropriate (and not overkill) for most biologists interested in phylogeography (~100 loci for ~100 individuals), and their results highlight several analytical challenges that lie ahead for researchers employing NGS techniques.     

A tale of two lineages: unexpected,long‐term persistence of the amphibian‐killing fungus in Brazil

For the past 17 years, scientists have been compiling a list of amphibian species susceptible to infection by the amphibian‐killing chytrid fungus, Batrachochytrium dendrobatidis (Bd), all over the world, with >500 species infected on every continent except Antarctica (Olson et al. 2013 ). Where Bd has been found, the impacts on amphibians has been one of two types: either Bd arrives into a naïve amphibian population followed by a mass die‐off and population declines (e.g. Lips et al. 2006 ), or Bd is present at some moderate prevalence, usually infecting many species but at apparently nonlethal intensities for a long time. In this issue of Molecular Ecology, Rodriguez et al. ( 2014 ) discover that the Atlantic Coastal Forest of Brazil is home to two Bd lineages: the Global Pandemic Lineage (Bd‐GPL) – the strain responsible for mass die‐offs and population declines – and a lineage endemic to Brazil (Bd‐Bz). Even more surprising was that both lineages have been present in this area for the past 100 years, making these the oldest records of Bd infecting amphibians. The team also described a moderate but steady prevalence of ~20% across all sampled anuran families for over 100 years, indicating that Brazil has been in an enzootic disease state for over a century. Most amphibians were infected with Bd‐GPL, suggesting this lineage may be a better competitor than Bd‐Bz or may be replacing the Bd‐Bz lineage. Rodriguez et al. ( 2014 ) also detected likely hybridization of the two Bd lineages, as originally described by Schloegel et al. ( 2012 ).     

Enalikter is not an annelid: homology,autapomorphies and the interpretation of problematic fossils

A megacheiran arthropod, Enalikter aphson, was recently described by Siveter et al. (2014) from the mid‐Silurian (late Wenlock) of Herefordshire. Previously, megacheirans had only been recognized from the Cambrian. Struck et al. (2015) considered the body plan of Enalikter to be incompatible with this affinity, arguing that many of the arthropod features were either not present or misinterpreted. Instead, they compared Enalikter to polychaete annelids, identifying characters from numerous polychaete lineages which they considered to be present in Enalikter. A reply to this critique by Siveter et al. (2015) reaffirmed arthropod affinities for Enalikter by presenting additional evidence for key arthropod features, such as arthropodized appendages. Here, we augment Siveter et al. by critically addressing the putative annelid characters of Enalikter presented by Struck et al. and additionally explore the morphological and phylogenetic implications of their hypothesis. We conclude that similarities between Enalikter and polychaetes are superficial and that character combinations proposed by Struck et al. are not present in any annelid, living or extinct. This taxon highlights the importance of using a phylogenetic framework for interpreting fossils that present unusual morphologies, such that proposed shared characters are hypotheses of homology rather than merely phenotypic similarities. Crucially, we argue that autapomorphic characters of subgroups of large taxa (like families or classes within phyla) should not be used to diagnose problematic fossils.     

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.     

It's time to listen: there is much to be learned from the sounds of tropical ecosystems

Knowledge that can be gained from acoustic data collection in tropical ecosystems is low‐hanging fruit. There is every reason to record and with every day, there are fewer excuses not to do it. In recent years, the cost of acoustic recorders has decreased substantially (some can be purchased for under US$50, e.g., Hill et al. 2018) and the technology needed to store and analyze acoustic data is continuously improving (e.g., Corrada Bravo et al. 2017, Xie et al. 2017). Soundscape recordings provide a permanent record of a site at a given time and contain a wealth of invaluable and irreplaceable information. Although challenges remain, failure to collect acoustic data now in tropical ecosystems would represent a failure to future generations of tropical researchers and the citizens that benefit from ecological research. In this commentary, we (1) argue for the need to increase acoustic monitoring in tropical systems; (2) describe the types of research questions and conservation issues that can be addressed with passive acoustic monitoring (PAM) using both short‐ and long‐term data in terrestrial and freshwater habitats; and (3) present an initial plan for establishing a global repository of tropical recordings.     

High‐throughput SNPs for all: genotyping‐in‐thousands

Understanding the genetic structure of species is essential for conservation. It is only with this information that managers, academics, user groups and land‐use planners can understand the spatial scale of migration and local adaptation, source‐sink dynamics and effective population size. Such information is essential for a multitude of applications including delineating management units, balancing management priorities, discovering cryptic species and implementing captive breeding programmes. Species can range from locally adapted by hundreds of metres (Pavey et al. 2010 ) to complete species panmixia (Côté et al. 2013 ). Even more remarkable is that this essential information can be obtained without fully sequenced or annotated genomes, but from mere (putatively) nonfunctional variants. First with allozymes, then microsatellites and now SNPs, this neutral genetic variation carries a wealth of information about migration and drift. For many of us, it may be somewhat difficult to remember our understanding of species conservation before the widespread usage of these useful tools. However most species on earth have yet to give us that ‘peek under the curtain’. With the current diversity on earth estimated to be nearly 9 million species (Mora et al. 2011 ), we have a long way to go for a comprehensive meta‐phylogeographic understanding. A method presented in this issue by Campbell and colleagues (Campbell et al. 2015 ) is a tool that will accelerate the pace in this area. Genotyping‐in‐thousands (GT‐seq) leverages recent advancements in sequencing technology to save many hours and dollars over previous methods to generate this important neutral genetic information.     

Don't throw the baby out with the bathwater: identifying and mapping paralogs in salmonids

Many eukaryotic genomes contain a large fraction of gene duplicates (or paralogs) as a result of ancient or recent whole‐genome duplications (Ohno 1970 ; Jaillon et al. 2004 ; Kellis et al. 2004 ). Identifying paralogs with NGS data is a pervasive problem in both ancient polyploids and neopolyploids. Likewise, paralogs are often treated as a nuisance that has to be detected and removed (Everett et al. 2012 ). In this issue of Molecular Ecology Resources, Waples et al. ( 2015 ) show that exclusion might not be necessary and how we may miss out on important genomic information in doing so. They present a novel statistical approach to detect paralogs based on the segregation of RAD loci in haploid offspring and test their method by constructing linkage maps with and without these duplicated loci in chum salmon, Oncorhynchus keta (Fig.  1 ). Their linkage map including the resolved paralogs shows that these are mostly located in the distal regions of several linkage groups. Particularly intriguing is their finding that these homoeologous regions appear impoverished in transposable elements (TE). Given the role that TE play in genome remodelling, it is noteworthy that these elements are of low abundance in regions showing residual tetrasomic inheritance. This raises the question whether re‐diploidization is constrained in these regions and whether they might have a role to play in salmonid speciation. This study provides an original approach to identifying duplicated loci in species with a pedigree, as well as providing a dense linkage map for chum salmon, and interesting insights into the retention of gene duplicates in an ancient polyploid.