More precisely biased: increasing the number of markers is not a silver bullet in genetic bottleneck testing

In response to our review of the use of genetic bottleneck tests in the conservation literature (Peery et al. 2012, Molecular Ecology, 21 , 3403–3418), Hoban et al. (2013, Molecular Ecology, in press) conducted population genetic simulations to show that the statistical power of genetic bottleneck tests can be increased substantially by sampling large numbers of microsatellite loci, as they suggest is now possible in the age of genomics. While we agree with Hoban and co‐workers in principle, sampling large numbers of microsatellite loci can dramatically increase the probability of committing type 1 errors (i.e. detecting a bottleneck in a stable population) when the mutation model is incorrectly assumed. Using conservative values for mutation model parameters can reduce the probability of committing type 1 errors, but doing so can result in significant losses in statistical power. Moreover, we believe that practical limitations associated with developing large numbers of high‐quality microsatellite loci continue to constrain sample sizes, a belief supported by a literature review of recent studies using next generation sequencing methods to develop microsatellite libraries. conclusion, we maintain that researchers employing genetic bottleneck tests should proceed with caution and carefully assess both statistical power and type 1 error rates associated with their study design.     

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.     

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.     

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.     

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.     

Novel aspects in the life cycle and biotrophic interactions in Pezizomycetes (Ascomycota,Fungi)

The ascomycete class Pezizomycetes (single order Pezizales) is known for its cup‐shaped fruit bodies and the evolution of edible truffles and morels, but little is known about the ontogeny and ecology of this large and ecologically diverse fungal group. In this issue of Molecular Ecology, Healy et al. ( 2013 ) make a great leap forward by describing and identifying asexual, anamorphic structures that produce mitotic spores in many ectomycorrhiza‐forming truffle and nontruffle species on soil surfaces worldwide (Fig.  1 ). Although such anamorphic forms have been reported sporadically from certain ectomycorrhizal and saprotrophic Pezizomycetes (e.g. Warcup 1990 ), Healy et al. ( 2013 ) demonstrate that these terricolous asexual forms are both taxonomically and geographically more widespread and, in fact, much more common than previously understood. We anticipate that deeper insight into other substrates, provided by molecular analyses of materials such as dead wood and seeds, is likely to reveal numerous anamorphs of saprotrophic and pathogenic Pezizomycetes as well (see Marek et al. 2009 ).     

Ready,Set…Poised!: Polycomb target genes are bound by poised RNA polymerase II throughout differentiation

In embryonic stem cells (ESCs), silent genes with major developmental functions display a unique epigenetic state in which strong and broad binding by Polycomb repressive complexes (PRCs) is accompanied by the presence of poised RNA polymerase II (RNAPII) and activating histone marks (e.g. H3K4me3) (Azuara et al, 2006 ; Bernstein et al, 2006 ; Stock et al, 2007 ; Brookes et al, 2012 ). It has been suggested that the plasticity and broad differentiation potential of pluripotent cells might rely, at least partly, on this unique epigenetic state (Bernstein et al, 2006 ; Stock et al, 2007 ). In their recent study, Pombo and colleagues (Ferrai et al, 2017 ) show that a similar epigenetic state can be found at a subset of major developmental genes throughout the differentiation of ESCs into neurons, providing novel and exciting insights into the molecular basis of cellular plasticity in differentiated cells.     

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

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

Bacteriophage exclusion,a new defense system

The ability to withstand viral predation is critical for survival of most microbes. Accordingly, a plethora of phage resistance systems has been identified in bacterial genomes (Labrie et al, 2010 ), including restriction‐modification systems (R‐M) (Tock & Dryden, 2005 ), abortive infection (Abi) (Chopin et al, 2005 ), Argonaute‐based interference (Swarts et al, 2014 ), as well as clustered regularly interspaced short palindromic repeats (CRISPR) and associated protein (Cas) adaptive immune system (CRISPR‐Cas) (Barrangou & Marraffini, 2014 ; Van der Oost et al, 2014 ). Predictably, the dark matter of bacterial genomes contains a wealth of genetic gold. A study published in this issue of The EMBO Journal by Goldfarb et al ( 2015 ) unveils bacteriophage exclusion (BREX) as a novel, widespread bacteriophage resistance system that provides innate immunity against virulent and temperate phage in bacteria.     

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.     

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.     

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.     

Artificial intelligence for microbial biotechnology: beyond the hype

It has been a landmark year for artificial intelligence (AI) and biotechnology. Perhaps the most noteworthy of these advances was Google DeepMind’s AlphaFold2 algorithm which smashed records in protein structure prediction (Jumper et al., 2021, Nature, 596, 583) complemented by progress made by other research groups around the globe (Baek et al., 2021, Science, 373, 871; Zheng et al., 2021, Proteins). For the first time in history, AI achieved protein structure models rivalling the accuracy of experimentally determined structures. The power of accurate protein structure prediction at our fingertips has countless implications for drug discovery, de novo protein design and fundamental research in chemical biology. While acknowledging the significance of these breakthroughs, this perspective aims to cut through the hype and examine some key limitations using AlphaFold2 as a lens to consider the broader implications of AI for microbial biotechnology for the next 15 years and beyond.     

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.     

No sex please,we're (in)breeding

A genetic trick allows induction of haploid maize plants by a process known as gynogenesis, which is a useful tool for breeders. In this issue of The EMBO Journal, Gilles et al ( 2017 ) show that loss of function of a patatin-like phospholipase A underlies the induction of gynogenesis, findings that were also made in two other recent studies (Kelliher et al, 2017 ; Liu et al, 2017 ).     

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.     

Population genomic analyses reveal possible drivers of population divergence

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