Genome‐wide evidence for divergent selection between populations of a major agricultural pathogen

The genetic and environmental homogeneity in agricultural ecosystems is thought to impose strong and uniform selection pressures. However, the impact of this selection on plant pathogen genomes remains largely unknown. We aimed to identify the proportion of the genome and the specific gene functions under positive selection in populations of the fungal wheat pathogen Zymoseptoria tritici. First, we performed genome scans in four field populations that were sampled from different continents and on distinct wheat cultivars to test which genomic regions are under recent selection. Based on extended haplotype homozygosity and composite likelihood ratio tests, we identified 384 and 81 selective sweeps affecting 4% and 0.5% of the 35 Mb core genome, respectively. We found differences both in the number and the position of selective sweeps across the genome between populations. Using a XtX‐based outlier detection approach, we identified 51 extremely divergent genomic regions between the allopatric populations, suggesting that divergent selection led to locally adapted pathogen populations. We performed an outlier detection analysis between two sympatric populations infecting two different wheat cultivars to identify evidence for host‐driven selection. Selective sweep regions harboured genes that are likely to play a role in successfully establishing host infections. We also identified secondary metabolite gene clusters and an enrichment in genes encoding transporter and protein localization functions. The latter gene functions mediate responses to environmental stress, including interactions with the host. The distinct gene functions under selection indicate that both local host genotypes and abiotic factors contributed to local adaptation.     

The Park Grass Experiment and next‐generation approaches: local adaptation of sweet vernal grass revisited

Long‐term ecological experiments provide unique opportunities to observe the effects of natural selection. The Park Grass Experiment at Rothamsted Experiment Station in Hertfordshire, UK, is the longest running ecological experiment that incorporates fertilization treatments and has been ongoing since 1856. In the 1970s, local adaptation was observed in the grass Anthoxanthum odoratum to the elevated soil aluminium levels of the fertilized plots. Gould et al. ( 2014 ) have utilized this system to reevaluate the extent of local adaptation, first documented nearly 45 years ago (Snaydon 1970 ), and to use emerging molecular approaches to identify candidate genes for the adaptation. From their work, they identify several plausible candidate loci for aluminium tolerance. This work shows the power of long‐term field‐based trials in a scientific age concentrated on rapidly emerging molecular techniques often utilized in short, narrowly focused laboratory or controlled environment experiments. The current study clearly illustrates the benefits gained by combining these molecular approaches within long‐term monitoring experiments that can be regularly revisited in a changing world and used to address questions on evolutionary scales.     

Genomic evidence for parallel adaptation to cities

Urban evolutionary biology is the study of rapid evolutionary change in response to humans and our uses of land to support city dwellers. Because cities are relatively modern additions to the natural world, research on urban evolution tends to focus on microevolutionary change that has happened across a few to many hundreds of generations. These questions still fall under the broad purview of evolutionary ecology. However, the severity, rapidity and replication of environmental changes that drive evolution in this context make it worthy of specific attention. Urban evolution provides the opportunity to study the earliest stages of evolution in a context that is scientifically interesting and societally important. The newness of urban populations and their proximity to natural populations also creates challenges when trying to detect population genetic change. In a From the Cover article in this issue of Molecular Ecology, Mueller et al. use whole genome resequencing data to address some of these challenges while exploring genetic changes associated with urbanization in three replicate urban‐rural burrowing owl (Athene cunicularia) populations. Combining multiple approaches across these sample sites Mueller et al. find evidence for selection on genes whose function is related to synapses, neuron projections, brain connectivity and cognitive function in general. That selection was parallel suggests that phenotypes related to brain processes were probably particularly important for urban adaptation.     

Rapid adaptation to climate change

In recent years, amid growing concerns that changing climate is affecting species distributions and ecosystems, predicting responses to rapid environmental change has become a major goal. In this issue, Franks and colleagues take a first step towards this objective (Franks et al. 2016). They examine genomewide signatures of selection in populations of Brassica rapa after a severe multiyear drought. Together with other authors, Franks had previously shown that flowering time was reduced after this particular drought and that the reduction was genetically encoded. Now, the authors have sequenced previously stored samples to compare allele frequencies before and after the drought and identify the loci with the most extreme shifts in frequencies. The loci they identify largely differ between populations, suggesting that different genetic variants may be responsible for reduction in flowering time in the two populations.     

Population genomics and the causes of local differentiation

Exactly 50 years ago, a revolution in empirical population genetics began with the introduction of methods for detecting allelic variation using protein electrophoresis (Throckmorton 1962; Hubby 1963; Lewontin & Hubby 1966). These pioneering scientists showed that populations are chock‐full of genetic variation. This variation was a surprise that required a re‐thinking of evolutionary genetic heuristics. Understanding the causes for the maintenance of this variation became and remains a major area of research. In the process of addressing the causes, this same group of scientists documented geographical genetic structure (Prakash et al. 1969), spawning the continued accumulation of what is now a huge case study catalogue of geographical differentiation (e.g. Loveless & Hamrick 1984; Linhart & Grant 1996). Geographical differentiation is clearly quite common. Yet, a truly general understanding of the patterns in and causes of spatial genetic structure across the genome remains elusive. To what extent is spatial structure driven by drift and phylogeography vs. geographical differences in environmental sources of selection? What proportion of the genome participates? A general understanding requires range‐wide data on spatial patterning of variation across the entire genome. In this issue of Molecular Ecology, Lasky et al. (2012) make important strides towards addressing these issues, taking advantage of three contemporary revolutions in evolutionary biology. Two are technological: high‐throughput sequencing and burgeoning computational power. One is cultural: open access to data from the community of scientists and especially data sets that result from large collaborative efforts. Together, these developments may at last put answers within reach.     

Paths to selection on life history loci in different natural environments across the native range of Arabidopsis thaliana

Selection on quantitative trait loci (QTL) may vary among natural environments due to differences in the genetic architecture of traits, environment‐specific allelic effects or changes in the direction and magnitude of selection on specific traits. To dissect the environmental differences in selection on life history QTL across climatic regions, we grew a panel of interconnected recombinant inbred lines (RILs) of Arabidopsis thaliana in four field sites across its native European range. For each environment, we mapped QTL for growth, reproductive timing and development. Several QTL were pleiotropic across environments, three colocalizing with known functional polymorphisms in flowering time genes (CRY2, FRI and MAF2‐5), but major QTL differed across field sites, showing conditional neutrality. We used structural equation models to trace selection paths from QTL to lifetime fitness in each environment. Only three QTL directly affected fruit number, measuring fitness. Most QTL had an indirect effect on fitness through their effect on bolting time or leaf length. Influence of life history traits on fitness differed dramatically across sites, resulting in different patterns of selection on reproductive timing and underlying QTL. In two oceanic field sites with high prereproductive mortality, QTL alleles contributing to early reproduction resulted in greater fruit production, conferring selective advantage, whereas alleles contributing to later reproduction resulted in larger size and higher fitness in a continental site. This demonstrates how environmental variation leads to change in both QTL effect sizes and direction of selection on traits, justifying the persistence of allelic polymorphism at life history QTL across the species range.     

Genome‐wide single‐generation signatures of local selection in the panmictic European eel

Next‐generation sequencing and the collection of genome‐wide data allow identifying adaptive variation and footprints of directional selection. Using a large SNP data set from 259 RAD‐sequenced European eel individuals (glass eels) from eight locations between 34 and 64^oN, we examined the patterns of genome‐wide genetic diversity across locations. We tested for local selection by searching for increased population differentiation using F_ST‐based outlier tests and by testing for significant associations between allele frequencies and environmental variables. The overall low genetic differentiation found (F_ST = 0.0007) indicates that most of the genome is homogenized by gene flow, providing further evidence for genomic panmixia in the European eel. The lack of genetic substructuring was consistent at both nuclear and mitochondrial SNPs. Using an extensive number of diagnostic SNPs, results showed a low occurrence of hybrids between European and American eel, mainly limited to Iceland (5.9%), although individuals with signatures of introgression several generations back in time were found in mainland Europe. Despite panmixia, a small set of SNPs showed high genetic differentiation consistent with single‐generation signatures of spatially varying selection acting on glass eels. After screening 50 354 SNPs, a total of 754 potentially locally selected SNPs were identified. Candidate genes for local selection constituted a wide array of functions, including calcium signalling, neuroactive ligand–receptor interaction and circadian rhythm. Remarkably, one of the candidate genes identified is PERIOD, possibly related to differences in local photoperiod associated with the >30° difference in latitude between locations. Genes under selection were spread across the genome, and there were no large regions of increased differentiation as expected when selection occurs within just a single generation due to panmixia. This supports the conclusion that most of the genome is homogenized by gene flow that removes any effects of diversifying selection from each new generation.     

Localizing FST outliers on a QTL map reveals evidence for large genomic regions of reduced gene exchange during speciation‐with‐gene‐flow

Populations that maintain phenotypic divergence in sympatry typically show a mosaic pattern of genomic divergence, requiring a corresponding mosaic of genomic isolation (reduced gene flow). However, mechanisms that could produce the genomic isolation required for divergence‐with‐gene‐flow have barely been explored, apart from the traditional localized effects of selection and reduced recombination near centromeres or inversions. By localizing F_ST outliers from a genome scan of wild pea aphid host races on a Quantitative Trait Locus (QTL) map of key traits, we test the hypothesis that between‐population recombination and gene exchange are reduced over large ‘divergence hitchhiking’ (DH) regions. As expected under divergence hitchhiking, our map confirms that QTL and divergent markers cluster together in multiple large genomic regions. Under divergence hitchhiking, the nonoutlier markers within these regions should show signs of reduced gene exchange relative to nonoutlier markers in genomic regions where ongoing gene flow is expected. We use this predicted difference among nonoutliers to perform a critical test of divergence hitchhiking. Results show that nonoutlier markers within clusters of F_ST outliers and QTL resolve the genetic population structure of the two host races nearly as well as the outliers themselves, while nonoutliers outside DH regions reveal no population structure, as expected if they experience more gene flow. These results provide clear evidence for divergence hitchhiking, a mechanism that may dramatically facilitate the process of speciation‐with‐gene‐flow. They also show the power of integrating genome scans with genetic analyses of the phenotypic traits involved in local adaptation and population divergence.     

Divergent selection on flowering time contributes to local adaptation in Mimulus guttatus populations

The timing of when to initiate reproduction is an important transition in any organism's life cycle. There is much variation in flowering time among populations, but we do not know to what degree this variation contributes to local adaptation. Here we use a reciprocal transplant experiment to examine the presence of divergent natural selection for flowering time and local adaptation between two distinct populations of Mimulus guttatus. We plant both parents and hybrids (to tease apart differences in suites of associated parental traits) between these two populations into each of the two native environments and measure floral, vegetative, life-history, and fitness characters to assess which traits are under selection at each site. Analysis of fitness components indicates that each of these plant populations is locally adapted. We obtain striking evidence for divergent natural selection on date of first flower production at these two sites. Early flowering is favored at the montane site, which is inhabited by annual plants and characterized by dry soils in midsummer, whereas intermediate (though later) flowering dates are selectively favored at the temperate coastal site, which is inhabited by perennial plants and is almost continually moist. Divergent selection on flowering time contributes to local adaptation between these two populations of M. guttatus, suggesting that genetic differentiation in the timing of reproduction may also serve as a partial reproductive isolating barrier to gene flow among populations.     

Rapid genome‐wide evolution in Brassica rapa populations following drought revealed by sequencing of ancestral and descendant gene pools

There is increasing evidence that evolution can occur rapidly in response to selection. Recent advances in sequencing suggest the possibility of documenting genetic changes as they occur in populations, thus uncovering the genetic basis of evolution, particularly if samples are available from both before and after selection. Here, we had a unique opportunity to directly assess genetic changes in natural populations following an evolutionary response to a fluctuation in climate. We analysed genome‐wide differences between ancestors and descendants of natural populations of Brassica rapa plants from two locations that rapidly evolved changes in multiple phenotypic traits, including flowering time, following a multiyear late‐season drought in California. These ancestor‐descendant comparisons revealed evolutionary shifts in allele frequencies in many genes. Some genes showing evolutionary shifts have functions related to drought stress and flowering time, consistent with an adaptive response to selection. Loci differentiated between ancestors and descendants (F_ST outliers) were generally different from those showing signatures of selection based on site frequency spectrum analysis (Tajima's D), indicating that the loci that evolved in response to the recent drought and those under historical selection were generally distinct. Very few genes showed similar evolutionary responses between two geographically distinct populations, suggesting independent genetic trajectories of evolution yielding parallel phenotypic changes. The results show that selection can result in rapid genome‐wide evolutionary shifts in allele frequencies in natural populations, and highlight the usefulness of combining resurrection experiments in natural populations with genomics for studying the genetic basis of adaptive evolution.     

Integration of Random Forest with population‐based outlier analyses provides insight on the genomic basis and evolution of run timing in Chinook salmon (Oncorhynchus tshawytscha)

Anadromous Chinook salmon populations vary in the period of river entry at the initiation of adult freshwater migration, facilitating optimal arrival at natal spawning. Run timing is a polygenic trait that shows evidence of rapid parallel evolution in some lineages, signifying a key role for this phenotype in the ecological divergence between populations. Studying the genetic basis of local adaptation in quantitative traits is often impractical in wild populations. Therefore, we used a novel approach, Random Forest, to detect markers linked to run timing across 14 populations from contrasting environments in the Columbia River and Puget Sound, USA. The approach permits detection of loci of small effect on the phenotype. Divergence between populations at these loci was then examined using both principle component analysis and F_ST outlier analyses, to determine whether shared genetic changes resulted in similar phenotypes across different lineages. Sequencing of 9107 RAD markers in 414 individuals identified 33 predictor loci explaining 79.2% of trait variance. Discriminant analysis of principal components of the predictors revealed both shared and unique evolutionary pathways in the trait across different lineages, characterized by minor allele frequency changes. However, genome mapping of predictor loci also identified positional overlap with two genomic outlier regions, consistent with selection on loci of large effect. Therefore, the results suggest selective sweeps on few loci and minor changes in loci that were detected by this study. Use of a polygenic framework has provided initial insight into how divergence in a trait has occurred in the wild.     

Fine‐scale temporal adaptation within a salmonid population: mechanism and consequences

We demonstrate a clear example of local adaptation of seasonal timing of spawning and embryo development. The consequence is a population of pink salmon that is segmented into spawning groups that use the same limited habitat. We synthesize published observations with results of new analyses to demonstrate that genetic variation of these traits results in survival differentials related to that variation, and that density‐dependent embryo mortality and seasonally variable juvenile mortality are a mechanism of selection. Most examples of local adaptation in natural systems depend on observed correlations between environments and fitness traits, but do not fully demonstrate local adaptation: that the trait is genetically determined, exhibits different fitness in common environments or across different environments, and its variation is mechanistically connected to fitness differences. The geographic or temporal scales of local adaptation often remain obscure. Here, we show that heritable, fine‐scale differences of timing of reproductive migration in a pink salmon (Oncorhynchus gorbuscha) resulted in temporal structure that persisted several generations; the differences enable a density‐dependent population to pack more spawners into limited spawning habitat, that is, enhance its fitness. A balanced trade‐off of survivals results because embryos from early‐migrating fish have a lower freshwater survival (harsh early physical conditions and disturbance by late spawners), but emigrant fry from late‐migrating fish have lower marine survivals (timing of their vernal emergence into the estuarine environment). Such fine‐scale local adaptations increase the genetic portfolio of the populations and may provide a buffer against the impacts of climate change.     

The genomic and physiological basis of life history variation in a butterfly metapopulation

Unravelling the mechanisms underlying variation in life history traits is of fundamental importance for our understanding of adaptation by natural selection. While progress has been made in mapping fitness-related phenotypes to genotypes, mainly in a handful of model organisms, functional genomic studies of life history adaptations are still in their infancy. In particular, despite a few notable exceptions, the genomic basis of life history variation in natural populations remains poorly understood. This is especially true for the genetic underpinnings of life history phenotypes subject to diversifying selection driven by ecological dynamics in patchy environments--as opposed to adaptations involving strong directional selection owing to major environmental changes, such as latitudinal gradients, extreme climatic events or transitions from salt to freshwater. In this issue of Molecular Ecology,Wheat et al. (2011) now make a significant leap forward by applying the tools of functional genomics to dispersal-related life history variation in a butterfly metapopulation. Using a combination of microarrays, quantitative PCR and physiological measurements, the authors uncover several metabolic and endocrine factors that likely contribute to the observed life history phenotypes. By identifying molecular candidate mechanisms of fitness variation maintained by dispersal dynamics in a heterogeneous environment,they also begin to address fascinating interactions between the levels of physiology, ecology and evolution.     

Landscape genomics and a common garden trial reveal adaptive differentiation to temperature across Europe in the tree species Alnus glutinosa

The adaptive potential of tree species to cope with climate change has important ecological and economic implications. Many temperate tree species experience a wide range of environmental conditions, suggesting high adaptability to new environmental conditions. We investigated adaptation to regional climate in the drought‐sensitive tree species Alnus glutinosa (Black alder), using a complementary approach that integrates genomic, phenotypic and landscape data. A total of 24 European populations were studied in a common garden and through landscape genomic approaches. Genotyping‐by‐sequencing was used to identify SNPs across the genome, resulting in 1990 SNPs. Although a relatively low percentage of putative adaptive SNPs was detected (2.86% outlier SNPs), we observed clear associations among outlier allele frequencies, temperature and plant traits. In line with the typical drought avoiding nature of A. glutinosa, leaf size varied according to a temperature gradient and significant associations with multiple outlier loci were observed, corroborating the ecological relevance of the observed outlier SNPs. Moreover, the lack of isolation by distance, the very low genetic differentiation among populations and the high intrapopulation genetic variation all support the notion that high gene exchange combined with strong environmental selection promotes adaptation to environmental cues.     

Into the wild—WAMBAM goes to Canada

The sixth Wild Animal Models Bi‐Annual Meeting was held in July 2017 in Québec, with 42 participants. This report documents the evolution of questions asked and approaches used in evolutionary quantitative genetic studies of wild populations in recent decades, and how these questions and approaches were represented at the recent meeting. We explore how ideas from previous meetings in this series have developed to their present states, and consider how the format of the meetings may be particularly useful at fostering the rapid development and proliferation of ideas and approaches.     

Unique evolutionary trajectories in repeated adaptation to hydrogen sulphide‐toxic habitats of a neotropical fish (Poecilia mexicana)

Replicated ecological gradients are prime systems to study processes of molecular evolution underlying ecological divergence. Here, we investigated the repeated adaptation of the neotropical fish Poecilia mexicana to habitats containing toxic hydrogen sulphide (H₂S) and compared two population pairs of sulphide‐adapted and ancestral fish by sequencing population pools of >200 individuals (Pool‐Seq). We inferred the evolutionary processes shaping divergence and tested the hypothesis of increase of parallelism from SNPs to molecular pathways. Coalescence analyses showed that the divergence occurred in the face of substantial bidirectional gene flow. Population divergence involved many short, widely dispersed regions across the genome. Analyses of allele frequency spectra suggest that differentiation at most loci was driven by divergent selection, followed by a selection‐mediated reduction of gene flow. Reconstructing allelic state changes suggested that selection acted mainly upon de novo mutations in the sulphide‐adapted populations. Using a corrected Jaccard index to quantify parallel evolution, we found a negligible proportion of statistically significant parallel evolution of J_corr = 0.0032 at the level of SNPs, divergent genome regions (J_corr = 0.0061) and genes therein (J_corr = 0.0091). At the level of metabolic pathways, the overlap was J_corr = 0.2545, indicating increasing parallelism with increasing level of biological integration. The majority of pathways contained positively selected genes in both sulphide populations. Hence, adaptation to sulphidic habitats necessitated adjustments throughout the genome. The largely unique evolutionary trajectories may be explained by a high proportion of de novo mutations driving the divergence. Our findings favour Gould's view that evolution is often the unrepeatable result of stochastic events with highly contingent effects.     

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.     

Widespread hybridization and bidirectional introgression in sympatric species of coral reef fish

Coral reefs are highly diverse ecosystems, where numerous closely related species often coexist. How new species arise and are maintained in these high geneflow environments have been long‐standing conundrums. Hybridization and patterns of introgression between sympatric species provide a unique insight into the mechanisms of speciation and the maintenance of species boundaries. In this study, we investigate the extent of hybridization between two closely related species of coral reef fish: the common coral trout (Plectropomus leopardus) and the bar‐cheek coral trout (Plectropomus maculatus). Using a complementary set of 25 microsatellite loci, we distinguish pure genotype classes from first‐ and later‐generation hybrids, identifying 124 interspecific hybrids from a collection of 2,991 coral trout sampled in inshore and mid‐shelf reefs of the southern Great Barrier Reef. Hybrids were ubiquitous among reefs, fertile and spanned multiple generations suggesting both ecological and evolutionary processes are acting to maintain species barriers. We elaborate on these finding to investigate the extent of genomic introgression and admixture from 2,271 SNP loci recovered from a ddRAD library of pure and hybrid individuals. An analysis of genomic clines on recovered loci indicates that 261 SNP loci deviate from a model of neutral introgression, of which 132 indicate a pattern of introgression consistent with selection favouring both hybrid and parental genotypes. Our findings indicate genome‐wide, bidirectional introgression between two sympatric species of coral reef fishes and provide further support to a growing body of evidence for the role of hybridization in the evolution of coral reef fishes.