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.     

Measuring microbial fitness in a field reciprocal transplant experiment

Microbial fitness is easy to measure in the laboratory, but difficult to measure in the field. Laboratory fitness assays make use of controlled conditions and genetically modified organisms, neither of which are available in the field. Among other applications, fitness assays can help researchers detect adaptation to different habitats or locations. We designed a competitive fitness assay to detect adaptation of Saccharomyces paradoxus isolates to the habitat they were isolated from (oak or larch leaf litter). The assay accurately measures relative fitness by tracking genotype frequency changes in the field using digital droplet PCR (DDPCR). We expected locally adapted S. paradoxus strains to increase in frequency over time when growing on the leaf litter type from which they were isolated. The DDPCR assay successfully detected fitness differences among S. paradoxus strains, but did not find a tendency for strains to be adapted to the habitat they were isolated from. Instead, we found that the natural alleles of the hexose transport gene we used to distinguish S. paradoxus strains had significant effects on fitness. The origin of a strain also affected its fitness: strains isolated from oak litter were generally fitter than strains from larch litter. Our results suggest that dispersal limitation and genetic drift shape S. paradoxus populations in the forest more than local selection does, although further research is needed to confirm this. Tracking genotype frequency changes using DDPCR is a practical and accurate microbial fitness assay for natural environments.     

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.     

The genomics of wild yeast populations sheds light on the domestication of man's best (micro) friend

The domestication of plants, animals and microbes by humans are the longest artificial evolution experiments ever performed. The study of these long‐term experiments can teach us about the genomics of adaptation through the identification of the genetic bases underlying the traits favoured by humans. In laboratory evolution, the characterization of the molecular changes that evolved specifically in some lineages is straightforward because the ancestors are readily available, for instance in the freezer. However, in the case of domesticated species, the ancestor is often missing, which leads to the necessity of going back to nature in order to infer the most likely ancestral state. Significant and relatively recent examples of this approach include wolves as the closest wild relative to domestic dogs (Axelsson et al. 2013) and teosinte as the closest relative to maize (reviewed in Hake & Ross‐Ibarra 2015). In both cases, the joint analysis of domesticated lineages and their wild cousins has been key in reconstructing the molecular history of their domestication. While the identification of closest wild relatives has been done for many plants and animals, these comparisons represent challenges for micro‐organisms. This has been the case for the budding yeast Saccharomyces cerevisiae, whose natural ecological niche is particularly challenging to define. For centuries, this unicellular fungus has been the cellular factory for wine, beer and bread crafting, and currently for bioethanol and drug production. While the recent development of genomics has lead to the identification of many genetic elements associated with important wine characteristics, the historical origin of some of the domesticated wine strains has remained elusive due to the lack of knowledge of their close wild relatives. In this issue of Molecular Ecology, Almeida et al. (2015) identified what is to date the closest known wild population of the wine yeast. This population is found associated with oak trees in Europe, presumably its natural host. Using population genomics analyses, Almeida and colleagues discovered that the initial divergence between natural and domesticated wine yeasts in the Mediterranean region took place around the early days of wine production. Surprisingly, genomic regions that are key to wine production today appeared not to be derived from these natural populations but from genes gained from other yeast species.     

Getting specific: making taxonomic and ecological sense of large sequencing data sets

Eukaryotic microbes comprise a diverse collection of phototrophic and heterotrophic creatures known to play fundamental roles in ecological processes. Some can be identified by light microscopy, generally the largest and with conspicuous shapes, while the smallest can be counted by epifluorescence microscopy or flow cytometry but remain largely unidentified. Microbial diversity studies greatly advanced with the analysis of phylogenetic markers sequenced from natural assemblages. Molecular surveys began in 1990 targeting marine bacterioplankton (Giovannoni et al. 1990 ) and first approached microbial eukaryotes in three studies published in 2001 (Díez et al. 2001 ; López‐García et al. 2001 ; Moon‐van der Staay et al. 2001 ). These seminal studies, based on cloning and Sanger sequencing the complete 18S rDNA, were critical for obtaining broad pictures of microbial diversity in contrasted habitats and for describing novel lineages by robust phylogenies, but were limited by the number of sequences obtained. So, inventories of species richness in a given sample and community comparisons through environmental gradients were very incomplete. These limitations have been overcome with the advent of high‐throughput sequencing (HTS) methods, initially 454‐pyrosequencing, today Illumina and soon others to come. In this issue of Molecular Ecology, Egge et al. ( 2015 ) show a nice example of the use of HTS to study the biodiversity and seasonal succession of a particularly important group of marine microbial eukaryotes, the haptophytes. Temporal changes were analysed first at the community level, then at the clade level, and finally at the lowest rank comparable to species. Interesting and useful ecological insights were obtained at each taxonomic scale. Haptophyte diversity differed along seasons in a systematic manner, with some species showing seasonal preferences and others being always present. Many of these species had no correspondence with known species, pointing out the high level of novelty in microbial assemblages, only accessible by molecular tools. Moreover, the number of species detected was limited, agreeing with a putative scenario of constrained evolutionary diversification in free‐living small eukaryotes. This study illustrates the potential of HTS to address ecological relevant questions in an accessible way by processing large data sets that, nonetheless, need to be treated with a fair understanding of their limitations.     

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.     

A population study of killer viruses reveals different evolutionary histories of two closely related Saccharomyces sensu stricto yeasts

Microbes have evolved ways of interference competition to gain advantage over their ecological competitors. The use of secreted killer toxins by yeast cells through acquiring double‐stranded RNA viruses is one such prominent example. Although the killer behaviour has been well studied in laboratory yeast strains, our knowledge regarding how killer viruses are spread and maintained in nature and how yeast cells co‐evolve with viruses remains limited. We investigated these issues using a panel of 81 yeast populations belonging to three Saccharomyces sensu stricto species isolated from diverse ecological niches and geographic locations. We found that killer strains are rare among all three species. In contrast, killer toxin resistance is widespread in Saccharomyces paradoxus populations, but not in Saccharomyces cerevisiae or Saccharomyces eubayanus populations. Genetic analyses revealed that toxin resistance in S. paradoxus is often caused by dominant alleles that have independently evolved in different populations. Molecular typing identified one M28 and two types of M1 killer viruses in those killer strains. We further showed that killer viruses of the same type could lead to distinct killer phenotypes under different host backgrounds, suggesting co‐evolution between the viruses and hosts in different populations. Taken together, our data suggest that killer viruses vary in their evolutionary histories even within closely related yeast species.     

Is colour polymorphism advantageous to populations and species?

I am writing in response to an article by Bolton, Rollins and Griffith (2015) entitled ‘The danger within: the role of genetic, behavioural and ecological factors in population persistence of colour polymorphic species’ that was recently published as an Opinion under the NEWS AND VIEWS section in Molecular Ecology. Bolton et al. (Molecular Ecology, 2015, 24 , 2907) argue that colour polymorphism may reduce population fitness and increase extinction risk and emphasize that this is contrary to predictions put forward by Forsman et al. (Ecology, 89 , 2008, 34) and Wennersten & Forsman (Biological Reviews 87 , 2012, 756) that the existence of multiple colour morphs with co‐adapted gene complexes and associated trait values may increase the ecological and evolutionary success of polymorphic populations and species. Bolton et al. (Molecular Ecology, 2015, 24 , 2907) further state that there is no clear evidence from studies of ‘true polymorphic species’ that polymorphism promotes population persistence. In response, I (i) challenge their classifications of polymorphisms and revisit the traditional definitions recognizing the dynamic nature of polymorphisms, (ii) review empirical studies that have examined whether and how polymorphism is associated with extinction risk, (iii) discuss the roles of trait correlations between colour pattern and other phenotypic dimensions for population fitness and (iv) highlight that the causes and mechanisms that influence the composition and maintenance of polymorphisms are different from the consequences of the polymorphic condition and how it may impact on aspects of ecological success and long‐term persistence of populations and species.     

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.     

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.     

It's a family act: the geminin triplets take center stage in motile ciliogenesis

The balance between proliferation and differentiation is a fundamental aspect of multicellular life. Perhaps nowhere is this delicate balance more palpable than in the multiciliated cells (MCCs) that line the respiratory tract, the ependyma, and the oviduct. These cells contain dozens to hundreds of motile cilia that beat in a concerted fashion to generate directed fluid flow over the tissue surface. Although MCCs have exited the cell cycle, remarkably, they retain the ability to duplicate their centrioles and to mature those centrioles into ciliary basal bodies—two features, which are known to be normally under strict cell cycle control (Firat‐Karalar & Stearns, 2014 ). How post‐mitotic MCCs retain this ability, remains unclear. In the past several months, four research articles, including one from Terré et al in this issue of The EMBO Journal, have described a vital role for the geminin coiled‐coil domain‐containing protein (Gemc1) in the MCC gene expression program in multiple tissues and organisms, that bring us closer to understanding this question (Kyrousi et al, 2015 ; Zhou et al, 2015 ; Arbi et al, 2016 ; Terré et al, 2016 ).     

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.     

Ecological genomics in full colour

Colour patterns in animals have long offered an opportunity to observe adaptive traits in natural populations. Colour plays myriad roles in interactions within and among species, from reproductive signalling to predator avoidance, leading to multiple targets of natural and sexual selection and opportunities for diversification. Understanding the genetic and developmental underpinnings of variation in colour promises a fuller understanding of these evolutionary processes, but the path to unravelling these connections can be arduous. The advent of genomic techniques suitable for nonmodel organisms is now beginning to light the way. Two new studies in this issue of Molecular Ecology use genomic sequencing of laboratory crosses to map colour traits in cichlid fishes, a remarkably diverse group in which coloration has played a major role in diversification. They illustrate how genomic approaches, specifically RAD sequencing, can rapidly identify both simple and more complex genetic variation underlying ecologically important traits. In the first, Henning et al. ( 2014 ) detect a single locus that appears to control in a Mendelian fashion the presence of horizontal stripes, a trait that has evolved in numerous cichlid lineages. In the second, Albertson et al. ( 2014 ) identify several genes and epistatic interactions affecting multiple colour traits, as well as a novel metric describing integration across colour traits. Albertson et al. ( 2014 ) go further, by quantifying differential expression of parental alleles at a candidate locus and by relating differentiation among natural populations at mapped loci to trait divergence. Herein lies the promise of ecological genomics – efficiently integrating genetic mapping of phenotypes with population genomic data to both identify functional genes and unravel their evolutionary history. These studies offer guidance on how genomic techniques can be tailored to a research question or study system, and they also add to the growing body of empirical examples addressing basic questions about how ecologically important traits evolve in natural populations.     

Trade‐offs and utility of alternative RADseq methods: Reply to Puritz et al.

Puritz et al. provide a review of several RADseq methodological approaches in response to our ‘Population Genomic Data Analysis’ workshop (Sept 2013) review (Andrews & Luikart 2014). We agree with Puritz et al. on the importance for researchers to thoroughly understand RADseq library preparation and data analysis when choosing an approach for answering their research questions. Some of us are currently using multiple RADseq protocols, and we agree that the different methods may offer advantages in different cases. Our workshop review did not intend to provide a thorough review of RADseq because the workshop covered a broad range of topics within the field of population genomics. Similarly, neither the response of Puritz et al. nor our comments here provide sufficient space to thoroughly review RADseq. Nonetheless, here we address some key points that we find unclear or potentially misleading in their evaluation of techniques.     

The interaction of Saccharomyces paradoxus with its natural competitors on oak bark

The natural history of the model yeast Saccharomyces cerevisiae is poorly understood and confounded by domestication. In nature, S. cerevisiae and its undomesticated relative S. paradoxus are usually found on the bark of oak trees, a habitat very different from wine or other human fermentations. It is unclear whether the oak trees are really the primary habitat for wild yeast, or whether this apparent association is due to biased sampling. We use culturing and high‐throughput environmental sequencing to show that S. paradoxus is a very rare member of the oak bark microbial community. We find that S. paradoxus can grow well on sterile medium made from oak bark, but that its growth is strongly suppressed when the other members of the community are present. We purified a set of twelve common fungal and bacterial species from the oak bark community and tested how each affected the growth of S. paradoxus in direct competition on oak bark medium at summer and winter temperatures, identifying both positive and negative interactions. One Pseudomonas species produces a diffusible toxin that suppresses S. paradoxus as effectively as either the whole set of twelve species together or the complete community present in nonsterilized oak medium. Conversely, one of the twelve species, Mucilaginibacter sp., had the opposite effect and promoted S. paradoxus growth at low temperatures. We conclude that, in its natural oak tree habitat, S. paradoxus is a rare species whose success depends on the much more abundant microbial species surrounding it.     

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

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.     

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.     

Prmt5: a guardian of the germline protects future generations

Primordial germ cells (PGCs) are the embryonic precursors of the germ cell lineage that form sperm and egg cells. It is of great importance to preserve the germline from DNA damage and potentially from epimutations in order to ensure the survival of future generations. Recent research highlights the role of the protein arginine methyltransferase 5 (PRMT5) as an important player in DNA protection during germline development in the mouse (Kim et al, 2014 & Li et al, 2015 ).