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.     

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

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.     

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.     

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.     

Everything old is new again: (linc)RNAs make proteins!

Long non‐coding RNAs have become the focus of considerable interest over the past few years. Intriguing novel functions have been reported for lincRNAs. Three recent papers identify lincRNAs that work in a more conventional way—encoding protein—in each case a small polypeptide with an interesting biological activity (Magny et al, 2013 ; Pauli et al, 2014 ), (Bazzini et al, 2014 ).     

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.     

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.     

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.     

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

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

Predicting the distribution of non‐vagile taxa: a response to McNyset and Blackburn (2006) and re‐evaluation of Stockman et al. (2006)

This paper addresses the issues raised by McNyset and Blackburn (2006 ) in their response to Stockman et al. (2006 ). Re‐evaluation of our published GARP analyses by McNyset and Blackburn showed that a much improved ecological niche model is obtained for predicting the distribution of the trapdoor spider genus Promyrmekiaphila in central/northern California. The improved niche model results in a substantially reduced omission error rate and a predictive model comparable to models obtained using other methods (GLM and BIOCLIM). However, the improved GARP models have a high commission error rate (> 0.75); consequently, the inferences regarding difficulties in modelling non‐vagile taxa drawn by Stockman et al. remain valid. Finally, we discuss other relatively minor criticisms of our study raised by McNyset and Blackburn and issues related to the peer review of our original paper.     

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

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.     

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.