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.     

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

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

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.     

A hot topic: the genetics of adaptation to geothermal vents in Mimulus guttatus

Identifying the individual loci and mutations that underlie adaptation to extreme environments has long been a goal of evolutionary biology. However, finding the genes that underlie adaptive traits is difficult for several reasons. First, because many traits and genes evolve simultaneously as populations diverge, it is difficult to disentangle adaptation from neutral demographic processes. Second, finding the individual loci involved in any trait is challenging given the respective limitations of quantitative and population genetic methods. In this issue of Molecular Ecology, Hendrick et al. (2016) overcome these difficulties and determine the genetic basis of microgeographic adaptation between geothermal vent and nonthermal populations of Mimulus guttatus in Yellowstone National Park. The authors accomplish this by combining population and quantitative genetic techniques, a powerful, but labour‐intensive, strategy for identifying individual causative adaptive loci that few studies have used (Stinchcombe & Hoekstra 2008 ). In a previous common garden experiment (Lekberg et al. 2012), thermal M. guttatus populations were found to differ from their closely related nonthermal neighbours in various adaptive phenotypes including trichome density. Hendrick et al. (2016) combine quantitative trait loci (QTL) mapping, population genomic scans for selection and admixture mapping to identify a single genetic locus underlying differences in trichome density between thermal and nonthermal M. guttatus. The candidate gene, R2R3 MYB, is homologous to genes involved in trichome development across flowering plants. The major trichome QTL, Tr14, is also involved in trichome density differences in an independent M. guttatus population comparison (Holeski et al. 2010) making this an example of parallel genetic evolution.     

Mapping the genomic architecture of ecological speciation in the wild: does linkage disequilibrium hold the key?

The hunt for the genes underlying ecological speciation has now closed in on a number of candidates, but making the link from genotype to phenotype continues to pose a significant challenge. This is partly because genetic studies in many systems remain impeded by long generation times or an inability to perform controlled crosses. Now, in this issue of Molecular Ecology, Malek et al. (2012) demonstrate the utility of a novel admixture mapping approach that can be used to identify genomic regions contributing to adaptive trait divergence between natural populations. Remarkably, they validate their approach by mapping traits associated with mate choice in a wild limnetic and benthic threespine stickleback (Gasterosteus aculeatus) species pair, finding several loci associated with male nuptial coloration and shape. While this study benefited from tried‐and‐true microsatellites in a well‐characterized species with a detailed genetic map (and genome sequence), the field is quickly moving towards the use of next‐generation sequencing, especially for nonmodel systems. The ability to characterize molecular polymorphisms for any system suggests that molecular ecologists working on virtually any species may benefit from applying Malek et al.'s approach, if naturally admixed populations are available.     

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

Selection on QTL and complex traits in complex environments

Understanding genetic variation for complex traits in heterogeneous environments is a fundamental problem in biology. In this issue of Molecular Ecology, Fournier‐Level et al. ( 2013 ) analyse quantitative trait loci (QTL) influencing ecologically important phenotypes in mapping populations of Arabidopsis thaliana grown in four habitats across its native European range. They used causal modelling to quantify the selective consequences of life history and morphological traits and QTL on components of fitness. They found phenology QTL colocalizing with known flowering time genes as well as novel loci. Most QTL influenced fitness via life history and size traits, rather than QTL having direct effects on fitness. Comparison of phenotypes among environments found no evidence for genetic trade‐offs for phenology or growth traits, but genetic trade‐offs for fitness resulted because flowering time had opposite fitness effects in different environments. These changes in QTL effects and selective consequences may maintain genetic variation among populations.     

The interaction of sexually and naturally selected traits in the adaptive radiations of cichlid fishes

The question of how genetic variation translates into organismal diversity has puzzled biologists for decades. Despite recent advances in evolutionary and developmental genetics, the mechanisms that underlie adaptation, diversification and evolutionary innovation remain largely unknown. The exceptionally diverse species flocks of cichlid fishes are textbook examples of adaptive radiation and explosive speciation and emerge as powerful model systems to study the genetic basis of animal diversification. East Africa's hundreds of endemic cichlid species are akin to a natural mutagenesis screen and differ greatly not only in ecologically relevant (hence naturally selected) characters such as mouth morphology and body shape, but also in sexually selected traits such as coloration. One of the most fascinating aspects of cichlid evolution is the frequent occurrence of evolutionary parallelisms, which has led to the question whether selection alone is sufficient to produce these parallel morphologies, or whether a developmental or genetic bias has influenced the direction of diversification. Here, I review fitness-relevant traits that could be responsible for the cichlids' evolutionary success and assess whether these were shaped by sexual or natural selection. I then focus on the interaction and the relative importance of sexual vs. natural selection in cichlid evolution. Finally, I discuss what is currently known about the genes underlying the morphogenesis of adaptively relevant traits and highlight the importance of the forthcoming cichlid genomes in the quest of the genetic basis of diversification in this group.     

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.     

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

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.     

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.     

Genomics and museum specimens

Nearly 25 years ago, Allan Wilson and colleagues isolated DNA sequences from museum specimens of kangaroo rats (Dipodomys panamintinus) and compared these sequences with those from freshly collected animals (Thomas et al. 1990 ). The museum specimens had been collected up to 78 years earlier, so the two samples provided a direct temporal comparison of patterns of genetic variation. This was not the first time DNA sequences had been isolated from preserved material, but it was the first time it had been carried out with a population sample. Population geneticists often try to make inferences about the influence of historical processes such as selection, drift, mutation and migration on patterns of genetic variation in the present. The work of Wilson and colleagues was important in part because it suggested a way in which population geneticists could actually study genetic change in natural populations through time, much the same way that experimentalists can do with artificial populations in the laboratory. Indeed, the work of Thomas et al. ( 1990 ) spawned dozens of studies in which museum specimens were used to compare historical and present‐day genetic diversity (reviewed in Wandeler et al. 2007 ). All of these studies, however, were limited by the same fundamental problem: old DNA is degraded into short fragments. As a consequence, these studies mostly involved PCR amplification of short templates, usually short stretches of mitochondrial DNA or microsatellites. In this issue, Bi et al. ( 2013 ) report a breakthrough that should open the door to studies of genomic variation in museum specimens. They used target enrichment (exon capture) and next‐generation (Illumina) sequencing to compare patterns of genetic variation in historic and present‐day population samples of alpine chipmunks (Tamias alpinus) (Fig. 1). The historic samples came from specimens collected in 1915, so the temporal span of this comparison is nearly 100 years.     

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.     

Bayesian parentage analysis reliably controls the number of false assignments in natural populations

Parentage analysis in natural populations is a powerful tool for addressing a wide range of ecological and evolutionary questions. However, identifying parent–offspring pairs in samples collected from natural populations is often more challenging than simply resolving the Mendelian pattern of shared alleles. For example, large numbers of pairwise comparisons and limited numbers of genetic markers can contribute to incorrect assignments, whereby unrelated individuals are falsely identified as parent–offspring pairs. Determining which parentage methods are the least susceptible to making false assignments is an important challenge facing molecular ecologists. In a recent paper, Harrison et al. (2013a) address this challenge by comparing three commonly used parentage methods, including a Bayesian approach, in order to explore the effects of varied proportions of sampled parents on the accuracy of parentage assignments. Unfortunately, Harrison et al. made a simple error in using the Bayesian approach, which led them to incorrectly conclude that this method could not control the rate of false assignment. Here, I briefly outline the basic principles behind the Bayesian approach, identify the error made by Harrison et al., and provide detailed guidelines as to how the method should be correctly applied. Furthermore, using the exact data from Harrison et al., I show that the Bayesian approach actually provides greater control over the number of false assignments than either of the other tested methods. Lastly, I conclude with a brief introduction to solomon , a recently updated version of the Bayesian approach that can account for genotyping error, missing data and false matching.     

Magadi tilapia ecological specialization: filling the early gap in the speciation continuum

Cichlid fish are well known for their high speciation rates, which are usually accompanied by spectacular and rapid diversification in eco‐morphological and secondary sexual traits. This is best illustrated by the famous repeated explosive radiations in the African Great Lakes Tanganyika, Malawi and Victoria, each lake harbouring several hundreds of mostly endemic species. Correspondingly, cichlids diversified very rapidly in many other lakes across their range. Although the larger radiations, unparalleled in vertebrates, are certainly the most intriguing, they are also the most intricate and difficult to address because of their complex nature. This is where smaller, simpler systems may prove to be the most useful. In this issue of Molecular Ecology, Kavembe et al. ( 2016 ) report very recent genetic diversification accompanied by ecological specialization in cichlids of the small and ecologically extreme Lake Magadi, in Kenya. Combining geometric morphometrics, stable isotope analysis, population genomics using RADSeq data and coalescent‐based modelling techniques, the authors characterize the eco‐morphological differences between genetically distinct populations of Magadi tilapia (Alcolapia grahami), which are consistent with the different environmental conditions they experience, and infer their history of divergence. The simplicity of the focal system and the use of a multidisciplinary approach make this work particularly important for our understanding of the early stages of speciation, in both cichlids and other organisms.     

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.     

Admixture increases diversity in managed honey bees: Reply to De la Rúa et al. (2013)

De la Rúa et al. (2013) express some concerns about the conclusions of our recent study showing that management increases genetic diversity of honey bees (Apis mellifera) by promoting admixture (Harpur et al. 2012). We provide a brief review of the literature on the population genetics of A. mellifera and show that we utilized appropriate sampling methods to estimate genetic diversity in the focal populations. Our finding of higher genetic diversity in two managed A. mellifera populations on two different continents is expected to be the norm given the large number of studies documenting admixture in honey bees. Our study focused on elucidating how management affects genetic diversity in honey bees, not on how to best manage bee colonies. We do not endorse the intentional admixture of honey bee populations, and we agree with De la Rúa et al. (2013) that native honey bee subspecies should be conserved.     

Pike (Esox lucius) could have been an exclusive human introduction to Ireland after all: a comment on Pedreschi et al. (2014)

A recent publication (Pedreschi et al., 2014, Journal of Biogeography, 41 , 548–560) casts doubt over the status of pike (Esox lucius) as a non‐native species in Ireland by reporting two distinct genetic groups of pike present: one a human introduction in the Middle Ages, the other hypothesized to result from natural colonization after the Last Glacial Maximum (LGM). While the existence of two groups is not disputed, the hypothesized natural colonization scenario conflicts with the sequence in which the islands of Britain and Ireland became isolated from Europe after the LGM. An alternative natural colonization scenario raised herein was rejected, leaving an earlier, two‐phase, human introduction of pike from Britain or Europe to Ireland as a realistic alternative hypothesis explaining the results of Pedreschi et al. (2014). This leaves the debates on human introduction versus natural colonization, introduced versus native species status, and pike management in Ireland wide open.     

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.