Search

The effect of past environmental changes on the demography and genetic diversity of natural populations remains a contentious issue and has rarely been investigated across multiple, phylogenetically distant species. Here, we perform comparative population genomic analyses and demographic inferences for seven widely distributed and ecologically contrasting European forest tree species based on concerted sampling of 164 populations across their natural ranges. For all seven species, the effective population size, N e , increased or remained stable over many glacial cycles and up to 15 million years in the most extreme cases. Surprisingly, the drastic environmental changes associated with the Pleistocene glacial cycles have had little impact on the level of genetic diversity of dominant forest tree species, despite major shifts in their geographic ranges. Based on their trajectories of N e over time, the seven tree species can be divided into three major groups, highlighting the importance of life history and range size in determining synchronous variation in genetic diversity over time. Altogether, our results indicate that forest trees have been able to retain their evolutionary potential over very long periods of time despite strong environmental changes., (© 2024. The Author(s).)

Closely related species often use the same genes to adapt to similar environments. However, we know little about why such genes possess increased adaptive potential and whether this is conserved across deeper evolutionary lineages. Adaptation to climate presents a natural laboratory to test these ideas, as even distantly related species must contend with similar stresses. Here, we re-analyse genomic data from thousands of individuals from 25 plant species as diverged as lodgepole pine and Arabidopsis (~300 Myr). We test for genetic repeatability based on within-species associations between allele frequencies in genes and variation in 21 climate variables. Our results demonstrate significant statistical evidence for genetic repeatability across deep time that is not expected under randomness, identifying a suite of 108 gene families (orthogroups) and gene functions that repeatedly drive local adaptation to climate. This set includes many orthogroups with well-known functions in abiotic stress response. Using gene co-expression networks to quantify pleiotropy, we find that orthogroups with stronger evidence for repeatability exhibit greater network centrality and broader expression across tissues (higher pleiotropy), contrary to the 'cost of complexity' theory. These gene families may be important in helping wild and crop species cope with future climate change, representing important candidates for future study., (© 2024. The Author(s).)

Genetic and genomic data are collected for a vast array of scientific and applied purposes. Despite mandates for public archiving, data are typically used only by the generating authors. The reuse of genetic and genomic datasets remains uncommon because it is difficult, if not impossible, due to non-standard archiving practices and lack of contextual metadata. But as the new field of macrogenetics is demonstrating, if genetic data and their metadata were more accessible and FAIR (findable, accessible, interoperable and reusable) compliant, they could be reused for many additional purposes. We discuss the main challenges with existing genetic and genomic data archives, and suggest best practices for archiving genetic and genomic data. Recognizing that this is a longstanding issue due to little formal data management training within the fields of ecology and evolution, we highlight steps that research institutions and publishers could take to improve data archiving., (© 2024. Springer Nature Limited.)

Detecting the extrinsic selective pressures shaping genomic variation is critical for a better understanding of adaptation and for forecasting evolutionary responses of natural populations to changing environmental conditions. With increasing availability of geo-referenced environmental data, landscape genomics provides unprecedented insights into how genomic variation and underlying gene functions affect traits potentially under selection. Yet, the robustness of genotype-environment associations used in landscape genomics remains tempered due to various limitations, including the characteristics of environmental data used, sampling designs employed, and statistical frameworks applied. Here, we argue that using complementary or new environmental data sources and well-informed sampling designs may help improve the detection of selective pressures underlying patterns of local adaptation in various organisms and environments., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Background: Plant mitogenomes vary widely in size and genomic architecture. Although hundreds of plant mitogenomes of angiosperm species have already been sequence-characterized, only a few mitogenomes are available from gymnosperms. Silver fir (Abies alba) is an economically important gymnosperm species that is widely distributed in Europe and occupies a large range of environmental conditions. Reference sequences of the nuclear and chloroplast genome of A. alba are available, however, the mitogenome has not yet been assembled and studied., Results: Here, we used paired-end Illumina short reads generated from a single haploid megagametophyte in combination with PacBio long reads from high molecular weight DNA of needles to assemble the first mitogenome sequence of A. alba. Assembly and scaffolding resulted in 11 mitogenome scaffolds, with the largest scaffold being 0.25 Mbp long. Two of the scaffolds displayed a potential circular structure supported by PCR. The total size of the A. alba mitogenome was estimated at 1.43 Mbp, similar to the size (1.33 Mbp) of a draft assembly of the Abies firma mitogenome. In total, 53 distinct genes of known function were annotated in the A. alba mitogenome, comprising 41 protein-coding genes, nine tRNA, and three rRNA genes. The proportion of highly repetitive elements (REs) was 0.168. The mitogenome seems to have a complex and dynamic structure featured by high combinatorial variation, which was specifically confirmed by PCR for the contig with the highest mapping coverage. Comparative analysis of all sequenced mitogenomes of gymnosperms revealed a moderate, but significant positive correlation between mitogenome size and proportion of REs., Conclusions: The A. alba mitogenome provides a basis for new comparative studies and will allow to answer important structural, phylogenetic and other evolutionary questions. Future long-read sequencing with higher coverage of the A. alba mitogenome will be the key to further resolve its physical structure. The observed positive correlation between mitogenome size and proportion of REs will be further validated once available mitogenomes of gymnosperms would become more numerous. To test whether a higher proportion of REs in a mitogenome leads to an increased recombination and higher structural complexity and variability is a prospective avenue for future research., (© 2022. The Author(s).)

Rapid human-induced environmental changes like climate warming represent a challenge for forest ecosystems. Due to their biological complexity and the long generation time of their keystone tree species, genetic adaptation in these ecosystems might not be fast enough to keep track with conditions changing at such a fast pace. The study of adaptation to environmental change and its genetic mechanisms is therefore key for ensuring a sustainable support and management of forests. The 4-day conference of the European Research Group EvolTree (https://www.evoltree.eu) on the topic of "Genomics and Adaptation in Forest Ecosystems" brought together over 130 scientists to present and discuss the latest developments and findings in forest evolutionary research. Genomic studies in forest trees have long been hampered by the lack of high-quality genomics resources and affordable genotyping methods. This has dramatically changed in the last few years; the conference impressively showed how such tools are now being applied to study past demography, adaptation and interactions with associated organisms. Moreover, genomic studies are now finally also entering the world of conservation and forest management, for example by measuring the value or cost of interspecific hybridization and introgression, assessing the vulnerability of species and populations to future change, or accurately delineating evolutionary significant units. The newly launched conference series of EvolTree will hopefully play a key role in the exchange and synthesis of such important investigations., Supplementary Information: The online version contains supplementary material available at 10.1007/s11295-022-01542-1., Competing Interests: Conflict of interestThe authors declare no competing interests., (© The Author(s) 2022.)

Many forest tree species have characteristics that make the study of their evolutionary ecology complex. For example, they are long-lived and thus have long generation times, and their often large, complex genomes have hampered establishing genomic resources. One way to tackle this challenge is to access multiple complementary data sources and analytical approaches when attempting to infer patterns of adaptive evolution. In the cover article of this issue of Molecular Ecology, Depardieu et al. (2021) combine large amounts of environmental, genomic, dendrochronological, and gene expression data in a common garden to explore the polygenic basis of drought resistance in white spruce (Picea glauca), a long-lived conifer. They identify candidate genes involved in growth and resistance to extreme drought events and show how multiple data sets may deliver complementary evidence to circumvent the manifold challenges of the research field., (© 2021 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Warmer climate and more frequent extreme droughts will pose major threats to forest ecosystems. Past demography processes due to post-glacial recolonization and adaptation to local environmental conditions are among the main contributors to genetic differentiation processes among provenances. Assessing the intra-specific variability of tree growth responses to such changes is crucial to explore a species' potential to cope with climate warming. We combined growth-related traits derived from tree-ring width series with neutral genetic information of 18 European provenances of silver fir (Abies alba Mill.) growing in two common garden experiments in Switzerland. Analyses based on neutral single nucleotide polymorphisms revealed that the studied provenances grouped into three longitudinal clusters. These three genetic clusters showed differences in growth traits (height and DBH), with the provenances from the eastern cluster exhibiting the highest growth. The Pyrenees cluster showed significantly lower recovery and resilience to the extreme drought of 2003 as well as lower values of growth autocorrelation. Q ST -F ST and correlation analyses with climate of provenance origin suggest that the differences among provenances found in some traits result from natural selection. Our study suggests that the last post-glacial re-colonization and natural selection are the major drivers explaining the intra-specific variability in growth of silver fir across Europe. These findings highlight the importance of combining dendroecology and genetic analyses on fitness-related traits to assess the potential of a species to cope with global environmental change and provide insights to support assisted gene flow to ensure the persistence of the species in European forests., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2021 The Authors. Published by Elsevier B.V. All rights reserved.)

In nature conservation, there is keen interest in predicting how populations will respond to environmental changes such as climate change. These predictions can help determine whether a population can be self-sustaining under future alterations of its habitat or whether it may require human intervention such as protection, restoration, or assisted migration. An increasingly popular approach in this respect is the concept of genomic offset, which combines genomic and environmental data from different time points and/or locations to assess the degree of possible maladaptation to new environmental conditions. Here, we argue that the concept of genomic offset holds great potential, but an exploration of its risks and limitations is needed to use it for recommendations in conservation or assisted migration. After briefly describing the concept, we list important issues to consider (e.g., statistical frameworks, population genetic structure, migration, independent evidence) when using genomic offset or developing these methods further. We conclude that genomic offset is an area of development that still lacks some important features and should be used in combination with other approaches to inform conservation measures., Competing Interests: The authors declare no conflict of interest., (© 2021 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

Background: Progress in the field of evolutionary forest ecology has been hampered by the huge challenge of phenotyping trees across their ranges in their natural environments, and the limitation in high-resolution environmental information., Findings: The GenTree Platform contains phenotypic and environmental data from 4,959 trees from 12 ecologically and economically important European forest tree species: Abies alba Mill. (silver fir), Betula pendula Roth. (silver birch), Fagus sylvatica L. (European beech), Picea abies (L.) H. Karst (Norway spruce), Pinus cembra L. (Swiss stone pine), Pinus halepensis Mill. (Aleppo pine), Pinus nigra Arnold (European black pine), Pinus pinaster Aiton (maritime pine), Pinus sylvestris L. (Scots pine), Populus nigra L. (European black poplar), Taxus baccata L. (English yew), and Quercus petraea (Matt.) Liebl. (sessile oak). Phenotypic (height, diameter at breast height, crown size, bark thickness, biomass, straightness, forking, branch angle, fructification), regeneration, environmental in situ measurements (soil depth, vegetation cover, competition indices), and environmental modeling data extracted by using bilinear interpolation accounting for surrounding conditions of each tree (precipitation, temperature, insolation, drought indices) were obtained from trees in 194 sites covering the species' geographic ranges and reflecting local environmental gradients., Conclusion: The GenTree Platform is a new resource for investigating ecological and evolutionary processes in forest trees. The coherent phenotyping and environmental characterization across 12 species in their European ranges allow for a wide range of analyses from forest ecologists, conservationists, and macro-ecologists. Also, the data here presented can be linked to the GenTree Dendroecological collection, the GenTree Leaf Trait collection, and the GenTree Genomic collection presented elsewhere, which together build the largest evolutionary forest ecology data collection available., (© The Author(s) 2021. Published by Oxford University Press GigaScience.)

The ongoing increase in global temperature affects biodiversity, especially in mountain regions where climate change is exacerbated. As sessile, long-lived organisms, trees are especially challenged in terms of adapting to rapid climate change. Here, we show that low rates of allele frequency shifts in Swiss stone pine (Pinus cembra) occurring near the treeline result in high genomic vulnerability to future climate warming, presumably due to the species' long generation time. Using exome sequencing data from adult and juvenile cohorts in the Swiss Alps, we found an average rate of allele frequency shift of 1.23 × 10 -2 /generation (i.e. 40 years) at presumably neutral loci, with similar rates for putatively adaptive loci associated with temperature (0.96 × 10 -2 /generation) and precipitation (0.91 × 10 -2 /generation). These recent shifts were corroborated by forward-in-time simulations at neutral and adaptive loci. Additionally, in juvenile trees at the colonisation front we detected alleles putatively beneficial under a future warmer and drier climate. Notably, the observed past rate of allele frequency shift in temperature-associated loci was decidedly lower than the estimated average rate of 6.29 × 10 -2 /generation needed to match a moderate future climate scenario (RCP4.5). Our findings suggest that species with long generation times may have difficulty keeping up with the rapid climate change occurring in high mountain areas and thus are prone to local extinction in their current main elevation range., (© 2020 John Wiley & Sons Ltd.)

It has long been discussed to what extent related species develop similar genetic mechanisms to adapt to similar environments. Most studies documenting such convergence have either used different lineages within species or surveyed only a limited portion of the genome. Here, we investigated whether similar or different sets of orthologous genes were involved in genetic adaptation of natural populations of three related plant species to similar environmental gradients in the Alps. We used whole-genome pooled population sequencing to study genome-wide SNP variation in 18 natural populations of three Brassicaceae (Arabis alpina, Arabidopsis halleri, and Cardamine resedifolia) from the Swiss Alps. We first de novo assembled draft reference genomes for all three species. We then ran population and landscape genomic analyses with ~3 million SNPs per species to look for shared genomic signatures of selection and adaptation in response to similar environmental gradients acting on these species. Genes with a signature of convergent adaptation were found at significantly higher numbers than expected by chance. The most closely related species pair showed the highest relative over-representation of shared adaptation signatures. Moreover, the identified genes of convergent adaptation were enriched for nonsynonymous mutations, suggesting functional relevance of these genes, even though many of the identified candidate genes have hitherto unknown or poorly described functions based on comparison with Arabidopsis thaliana. We conclude that adaptation to heterogeneous Alpine environments in related species is partly driven by convergent evolution, but that most of the genomic signatures of adaptation remain species-specific., (© 2020 John Wiley & Sons Ltd.)

It is generally accepted that the spatial distribution of neutral genetic diversity within a species' native range mostly depends on effective population size, demographic history, and geographic position. However, it is unclear how genetic diversity at adaptive loci correlates with geographic peripherality or with habitat suitability within the ecological niche. Using exome-wide genomic data and distribution maps of the Alpine range, we first tested whether geographic peripherality correlates with four measures of population genetic diversity at > 17,000 SNP loci in 24 Alpine populations (480 individuals) of Swiss stone pine (Pinus cembra) from Switzerland. To distinguish between neutral and adaptive SNP sets, we used four approaches (two gene diversity estimates, F ST outlier test, and environmental association analysis) that search for signatures of selection. Second, we established ecological niche models for P. cembra in the study range and investigated how habitat suitability correlates with genetic diversity at neutral and adaptive loci. All estimates of neutral genetic diversity decreased with geographic peripherality, but were uncorrelated with habitat suitability. However, heterozygosity (H e ) at adaptive loci based on Tajima's D declined significantly with increasingly suitable conditions. No other diversity estimates at adaptive loci were correlated with habitat suitability. Our findings suggest that populations at the edge of a species' geographic distribution harbour limited neutral genetic diversity due to demographic properties. Moreover, we argue that populations from suitable habitats went through strong selection processes, are thus well adapted to local conditions, and therefore exhibit reduced genetic diversity at adaptive loci compared to populations at niche margins., (© 2020 John Wiley & Sons Ltd.)

Long generation times have been suggested to hamper rapid genetic adaptation of organisms to changing environmental conditions. We examined if environmental memory of the parental Scots pines (Pinus sylvestris L.) drive offspring survival and growth. We used seeds from trees growing under naturally dry conditions (control), irrigated trees (irrigated from 2003 to 2016), and formerly irrigated trees ("irrigation stop"; irrigated from 2003-2013; control condition since 2014). We performed two experiments, one under controlled greenhouse conditions and one at the experimental field site. In the greenhouse, the offspring from control trees exposed regularly to drought were more tolerant to hot-drought conditions than the offspring from irrigated trees and showed lower mortality even though there was no genetic difference. However, under optimal conditions (high water supply and full sunlight), these offspring showed lower growth and were outperformed by the offspring of the irrigated trees. This different offspring growth, with the offspring of the "irrigation-stop" trees showing intermediate responses, points to the important role of transgenerational memory for the long-term acclimation of trees. Such memory effects, however, may be overridden by climatic extremes during germination and early growth stages such as the European 2018 mega-drought that impacted our field experiment., (© 2020 The Authors. Plant, Cell & Environment published by John Wiley & Sons Ltd.)

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Background and Aims: Hybridization and introgression play an important role in the evolution and diversification of plants. To assess the degree of past and current hybridization, the level of genetic admixture in populations needs to be investigated. Ongoing hybridization and blurred species separation have made it challenging to assign European white oak taxa based on leaf morphology and/or genetic markers and to assess the level of admixture. Therefore, there is a need for powerful markers that differentiate between taxa. Here, we established a condensed set of single-nucleotide polymorphism (SNP) markers to reliably differentiate between the three most common oak species in temperate European forests (Quercus robur, Q. petraea, Q. pubescens) and to assess the degree of admixture in a large set of selected Swiss populations., Methods: A training set of 194 presumably pure reference samples from Switzerland and Europe was used to assign 633 test individuals with two different approaches (population genetic-based/Bayesian vs. assumption-free/discriminative classifier) using 58 selected SNPs from coding regions. Admixture was calculated at the individual and population level with the Shannon diversity index based on individual assignment probabilities., Key Results: Depending on the approach, 97.5-100 % of training individuals were assigned correctly, and additional analyses showed that the established SNP set could be further reduced while maintaining its discriminatory power. The two assignment approaches showed high overlap (99 %) in assigning training individuals and slightly less overlap in test individuals (84 %). Levels of admixture varied widely among populations. Mixed stands of Q. petraea and Q. pubescens revealed much higher degrees of admixture than mixed stands of the other two taxon pairs, accentuating high levels of gene flow between these two taxa in Switzerland., Conclusions: Our set of SNPs warrants reliable taxon discrimination with great potential for further applications. We show that the three European white oak taxa have largely retained their species integrity in Switzerland despite high levels of admixture., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Annals of Botany Company. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Global climate change (GCC) increasingly threatens biodiversity through the loss of species, and the transformation of entire ecosystems. Many species are challenged by the pace of GCC because they might not be able to respond fast enough to changing biotic and abiotic conditions. Species can respond either by shifting their range, or by persisting in their local habitat. If populations persist, they can tolerate climatic changes through phenotypic plasticity, or genetically adapt to changing conditions depending on their genetic variability and census population size to allow for de novo mutations. Otherwise, populations will experience demographic collapses and species may go extinct. Current approaches to predicting species responses to GCC begin to combine ecological and evolutionary information for species distribution modelling. Including an evolutionary dimension will substantially improve species distribution projections which have not accounted for key processes such as dispersal, adaptive genetic change, demography, or species interactions. However, eco-evolutionary models require new data and methods for the estimation of a species' adaptive potential, which have so far only been available for a small number of model species. To represent global biodiversity, we need to devise large-scale data collection strategies to define the ecology and evolutionary potential of a broad range of species, especially of keystone species of ecosystems. We also need standardized and replicable modelling approaches that integrate these new data to account for eco-evolutionary processes when predicting the impact of GCC on species' survival. Here, we discuss different genomic approaches that can be used to investigate and predict species responses to GCC. This can serve as guidance for researchers looking for the appropriate experimental setup for their particular system. We furthermore highlight future directions for moving forward in the field and allocating available resources more effectively, to implement mitigation measures before species go extinct and ecosystems lose important functions., (© 2019 The Authors. Evolution Letters published by Wiley Periodicals, Inc. on behalf of Society for the Study of Evolution (SSE) and European Society for Evolutionary Biology (ESEB).)

The dataset presented here was collected by the GenTree project (EU-Horizon 2020), which aims to improve the use of forest genetic resources across Europe by better understanding how trees adapt to their local environment. This dataset of individual tree-core characteristics including ring-width series and whole-core wood density was collected for seven ecologically and economically important European tree species: silver birch (Betula pendula), European beech (Fagus sylvatica), Norway spruce (Picea abies), European black poplar (Populus nigra), maritime pine (Pinus pinaster), Scots pine (Pinus sylvestris), and sessile oak (Quercus petraea). Tree-ring width measurements were obtained from 3600 trees in 142 populations and whole-core wood density was measured for 3098 trees in 125 populations. This dataset covers most of the geographical and climatic range occupied by the selected species. The potential use of it will be highly valuable for assessing ecological and evolutionary responses to environmental conditions as well as for model development and parameterization, to predict adaptability under climate change scenarios.

Silver fir ( Abies alba Mill.) is a keystone conifer of European montane forest ecosystems that has experienced large fluctuations in population size during during the Quaternary and, more recently, due to land-use change. To forecast the species' future distribution and survival, it is important to investigate the genetic basis of adaptation to environmental change, notably to extreme events. For this purpose, we here provide a first draft genome assembly and annotation of the silver fir genome, established through a community-based initiative. DNA obtained from haploid megagametophyte and diploid needle tissue was used to construct and sequence Illumina paired-end and mate-pair libraries, respectively, to high depth. The assembled A. alba genome sequence accounted for over 37 million scaffolds corresponding to 18.16 Gb, with a scaffold N50 of 14,051 bp. Despite the fragmented nature of the assembly, a total of 50,757 full-length genes were functionally annotated in the nuclear genome. The chloroplast genome was also assembled into a single scaffold (120,908 bp) that shows a high collinearity with both the A. koreana and A. sibirica complete chloroplast genomes. This first genome assembly of silver fir is an important genomic resource that is now publicly available in support of a new generation of research. By genome-enabling this important conifer, this resource will open the gate for new research and more precise genetic monitoring of European silver fir forests., (Copyright © 2019 Mosca et al.)

Despite decreasing sequencing costs, whole-genome sequencing for population-based genome scans for selection is still prohibitively expensive for organisms with large genomes. Moreover, the repetitive nature of large genomes often represents a challenge in bioinformatic and downstream analyses. Here, we use in-depth transcriptome sequencing to design probes for exome capture in Swiss stone pine (Pinus cembra), a conifer with an estimated genome size of 29.3 Gbp and no reference genome available. We successfully applied around 55,000 self-designed probes, targeting 25,000 contigs, to DNA pools of seven populations from the Swiss Alps and identified >160,000 SNPs in around 15,000 contigs. The probes performed equally well in pools of the closely related species Pinus sibirica; in both species, more than 70% of the targeted contigs were sequenced at a depth ≥40× (number of haplotypes in the pool). However, a thorough analysis of individually sequenced P. cembra samples indicated that a majority of the contigs (63%) represented multi-copy genes. We therefore removed paralogous contigs based on heterozygote excess and deviation from allele balance. Without putatively paralogous contigs, allele frequencies of population pools represented accurate estimates of individually determined allele frequencies. We show that inferences of neutral and adaptive genetic variation may be biased when not accounting for such multi-copy genes. Without individual genotype data, it would have been nearly impossible to recognize and deal with the problem of multi-copy contigs. We advocate to put more emphasis on identifying paralogous loci, which will be facilitated by the establishment of additional high-quality reference genomes., (© 2018 John Wiley & Sons Ltd.)

When plants adapt to local environments, strong signatures of selection are expected in the genome, particularly in high-stress environments such as trace metal element enriched (metalliferous) soils. Using Arabidopsis halleri, a model species for metal homeostasis and adaptation to extreme environments, we identifid genes, gene variants, and pathways that are associated with soil properties and may thus contribute to adaptation to high concentrations of trace metal elements. We analysed whole-genome Pool-seq data from two metallicolous (from metalliferous soils) and two non-metallicolous populations (in total 119 individuals) and associated allele frequencies of the identified single-nucleotide polymorphisms (SNPs) with soil variables measured on site. Additionally, we accounted for polygenic adaptation by searching for gene pathways showing enrichment of signatures of selection. Out of >2.5 million SNPs, we identified 57 SNPs in 19 genes that were significantly associated with soil variables and are members of three enriched pathways. At least three of these candidate genes and pathways are involved in transmembrane transport and/or associated with responses to various stresses such as oxidative stress. We conclude that both allocation and detoxification processes play a crucial role in A. halleri for coping with these unfavourable conditions.

Over the last decade, the genomic revolution has offered the possibility to generate tremendous amounts of data that contain valuable information on the genetic basis of phenotypic traits, such as those linked to human diseases or those that allow for species to adapt to a changing environment. Most ecologically relevant traits are controlled by a large number of genes with small individual effects on trait variation, but that are connected with one another through complex developmental, metabolic and biochemical networks. As a result, it has recently been suggested that most adaptation events in natural populations are reached via correlated changes at multiple genes at a time, for which the name polygenic adaptation has been coined. The current challenge is to develop methods to extract the relevant information from genomic data to detect the signature of polygenic evolutionary change. The symposium entitled "Detecting the Genomic Signal of Polygenic Adaptation and the Role of Epistasis in Evolution" held in 2017 at the University of Zürich aimed at reviewing our current state of knowledge. In this review, we use the talks of the invited speakers to summarize some of the most recent developments in this field., (© 2018 John Wiley & Sons Ltd.)

Background: Microsatellite markers are widely used for estimating genetic diversity within and differentiation among populations. However, it has rarely been tested whether such estimates are useful proxies for genome-wide patterns of variation and differentiation. Here, we compared microsatellite variation with genome-wide single nucleotide polymorphisms (SNPs) to assess and quantify potential marker-specific biases and derive recommendations for future studies. Overall, we genotyped 180 Arabidopsis halleri individuals from nine populations using 20 microsatellite markers. Twelve of these markers were originally developed for Arabidopsis thaliana (cross-species markers) and eight for A. halleri (species-specific markers). We further characterized 2 million SNPs across the genome with a pooled whole-genome re-sequencing approach (Pool-Seq)., Results: Our analyses revealed that estimates of genetic diversity and differentiation derived from cross-species and species-specific microsatellites differed substantially and that expected microsatellite heterozygosity (SSR-H e ) was not significantly correlated with genome-wide SNP diversity estimates (SNP-H e and θ Watterson ) in A. halleri. Instead, microsatellite allelic richness (A r ) was a better proxy for genome-wide SNP diversity. Estimates of genetic differentiation among populations (F ST ) based on both marker types were correlated, but microsatellite-based estimates were significantly larger than those from SNPs. Possible causes include the limited number of microsatellite markers used, marker ascertainment bias, as well as the high variance in microsatellite-derived estimates. In contrast, genome-wide SNP data provided unbiased estimates of genetic diversity independent of whether genome- or only exome-wide SNPs were used. Further, we inferred that a few thousand random SNPs are sufficient to reliably estimate genome-wide diversity and to distinguish among populations differing in genetic variation., Conclusions: We recommend that future analyses of genetic diversity within and differentiation among populations use randomly selected high-throughput sequencing-based SNP data to draw conclusions on genome-wide diversity patterns. In species comparable to A. halleri, a few thousand SNPs are sufficient to achieve this goal.

Testing how populations are locally adapted and predicting their response to their future environment is of key importance in view of climate change. Landscape genomics is a powerful approach to investigate genes and environmental factors involved in local adaptation. In a pooled amplicon sequencing approach of 94 genes in 71 populations, we tested whether >3500 single nucleotide polymorphisms (SNPs) in the three most common oak species in Switzerland (Quercus petraea, Q. pubescens, Q. robur) show an association with abiotic factors related to local topography, historical climate and soil characteristics. In the analysis including all species, the most frequently associated environmental factors were those best describing the habitats of the species. In the species-specific analyses, the most important environmental factors and associated SNPs greatly differed among species. However, we identified one SNP and seven genes that were associated with the same environmental factor across all species. We finally used regressions of allele frequencies of the most strongly associated SNPs along environmental gradients to predict the risk of nonadaptedness (RONA), which represents the average change in allele frequency at climate-associated loci theoretically required to match future climatic conditions. RONA is considerable for some populations and species (up to 48% in single populations) and strongly differs among species. Given the long generation time of oaks, some of the required allele frequency changes might not be realistic to achieve based on standing genetic variation. Hence, future adaptedness requires gene flow or planting of individuals carrying beneficial alleles from habitats currently matching future climatic conditions., (© 2016 John Wiley & Sons Ltd.)

In nature, host individuals are commonly simultaneously infected with more than one genotype of the same parasite species. These co-infecting parasites often interact, which can affect their fitness and shape host-parasite ecology and evolution. Many of such interactions take place through competition for limited host resources. Therefore, variation in ecological factors modifying the host resource level could be important in determining the intensity of competition and the outcome of co-infections. We tested this hypothesis by measuring the relative reproductive success of co-infecting genotypes of the trematode parasite Diplostomum pseudospathaceum in its snail host Lymnaea stagnalis while experimentally manipulating snail resource level using contrasting feeding treatments (ad libitum food supply, no food). We found that food deprivation constrained the overall parasite within-host reproduction as the release of parasite transmission stages (cercariae) was reduced. This indicates intensified competition among the parasite genotypes. The genotypic composition of the released cercariae, however, was not affected by the feeding treatments. This suggests that in this system, the relative reproductive success of co-infecting parasite genotypes, which is an important component determining their fitness, is robust to variation in ecological factors modifying the strength of resource competition., (Copyright © 2015 Elsevier B.V. All rights reserved.)

Landscape genomics is an emerging research field that aims to identify the environmental factors that shape adaptive genetic variation and the gene variants that drive local adaptation. Its development has been facilitated by next-generation sequencing, which allows for screening thousands to millions of single nucleotide polymorphisms in many individuals and populations at reasonable costs. In parallel, data sets describing environmental factors have greatly improved and increasingly become publicly accessible. Accordingly, numerous analytical methods for environmental association studies have been developed. Environmental association analysis identifies genetic variants associated with particular environmental factors and has the potential to uncover adaptive patterns that are not discovered by traditional tests for the detection of outlier loci based on population genetic differentiation. We review methods for conducting environmental association analysis including categorical tests, logistic regressions, matrix correlations, general linear models and mixed effects models. We discuss the advantages and disadvantages of different approaches, provide a list of dedicated software packages and their specific properties, and stress the importance of incorporating neutral genetic structure in the analysis. We also touch on additional important aspects such as sampling design, environmental data preparation, pooled and reduced-representation sequencing, candidate-gene approaches, linearity of allele-environment associations and the combination of environmental association analyses with traditional outlier detection tests. We conclude by summarizing expected future directions in the field, such as the extension of statistical approaches, environmental association analysis for ecological gene annotation, and the need for replication and post hoc validation studies., (© 2015 John Wiley & Sons Ltd.)

Sequencing of pooled samples (Pool-Seq) using next-generation sequencing technologies has become increasingly popular, because it represents a rapid and cost-effective method to determine allele frequencies for single nucleotide polymorphisms (SNPs) in population pools. Validation of allele frequencies determined by Pool-Seq has been attempted using an individual genotyping approach, but these studies tend to use samples from existing model organism databases or DNA stores, and do not validate a realistic setup for sampling natural populations. Here we used pyrosequencing to validate allele frequencies determined by Pool-Seq in three natural populations of Arabidopsis halleri (Brassicaceae). The allele frequency estimates of the pooled population samples (consisting of 20 individual plant DNA samples) were determined after mapping Illumina reads to (i) the publicly available, high-quality reference genome of a closely related species (Arabidopsis thaliana) and (ii) our own de novo draft genome assembly of A. halleri. We then pyrosequenced nine selected SNPs using the same individuals from each population, resulting in a total of 540 samples. Our results show a highly significant and accurate relationship between pooled and individually determined allele frequencies, irrespective of the reference genome used. Allele frequencies differed on average by less than 4%. There was no tendency that either the Pool-Seq or the individual-based approach resulted in higher or lower estimates of allele frequencies. Moreover, the rather high coverage in the mapping to the two reference genomes, ranging from 55 to 284x, had no significant effect on the accuracy of the Pool-Seq. A resampling analysis showed that only very low coverage values (below 10-20x) would substantially reduce the precision of the method. We therefore conclude that a pooled re-sequencing approach is well suited for analyses of genetic variation in natural populations.

Natural genetic variation is essential for the adaptation of organisms to their local environment and to changing environmental conditions. Here, we examine genomewide patterns of nucleotide variation in natural populations of the outcrossing herb Arabidopsis halleri and associations with climatic variation among populations in the Alps. Using a pooled population sequencing (Pool-Seq) approach, we discovered more than two million SNPs in five natural populations and identified highly differentiated genomic regions and SNPs using FST -based analyses. We tested only the most strongly differentiated SNPs for associations with a nonredundant set of environmental factors using partial Mantel tests to identify topo-climatic factors that may underlie the observed footprints of selection. Possible functions of genes showing signatures of selection were identified by Gene Ontology analysis. We found 175 genes to be highly associated with one or more of the five tested topo-climatic factors. Of these, 23.4% had unknown functions. Genetic variation in four candidate genes was strongly associated with site water balance and solar radiation, and functional annotations were congruent with these environmental factors. Our results provide a genomewide perspective on the distribution of adaptive genetic variation in natural plant populations from a highly diverse and heterogeneous alpine environment., (© 2013 John Wiley & Sons Ltd.)

Vertebrate hosts often defend themselves against several co-infecting parasite genotypes simultaneously. This has important implications for the ecological dynamics and the evolution of host defence systems and parasite strategies. For example, it can drive the specificity of the adaptive immune system towards high genotype-specificity or cross-reactivity against several parasite genotypes depending on the sequence and probability of re-infections. However, to date, there is very little evidence on these interactions outside mammalian disease literature. In this study we asked whether genotype-specific or cross-reactive responses dominate in the adaptive immune system of a fish host towards a common macroparasite. In other words, we investigated if the infection success of a parasite genotype is influenced by the immunization genotype. We reciprocally immunized and re-exposed rainbow trout (Oncorhynchus mykiss) to a range of genotypes of the trematode eye fluke Diplostomum pseudospathaceum, and measured infection success of the parasite. We found that the infection success of the parasite genotypes in the re-exposure did not depend on the immunization genotype. While immunization reduced average infection success by 31%, the reduction was not larger against the initial immunization genotype. Our results suggest significant cross-reactivity, which may be advantageous for the host in genetically diverse re-exposures and have significant evolutionary implications for parasite strategies. Overall, our study is among the first to demonstrate cross-reactivity of adaptive immunity against genetically diverse macroparasites with complex life cycles.

Characterizing genetic variation in parasite transmission traits and its contribution to parasite vigor is essential for understanding the evolution of parasite life-history traits. We measured genetic variation in output, activity, survival, and infection success of clonal transmission stages (cercaria larvae) of a complex life cycle parasite (Diplostomum pseudospathaceum). We further tested if variation in host nutritional stage had an effect on these traits by keeping hosts on limited or ad libitum diet. The traits we measured were highly variable among parasite genotypes indicating significant genetic variation in these life-history traits. Traits were also phenotypically variable, for example, there was significant variation in the measured traits over time within each genotype. However, host nutritional stage had no effect on the parasite traits suggesting that a short-term reduction in host resources was not limiting the cercarial output or performance. Overall, these results suggest significant interclonal and phenotypic variation in parasite transmission traits that are not affected by host nutritional status.

Parasite communities of fishes are known to respond directly to the abiotic environment of the host, for example, to water quality and water temperature. Biotic factors are also important as they affect the exposure profile through heterogeneities in parasite distribution in the environment. Parasites in a particular environment may pose a strong selection on fish. For example, ecological differences in selection by parasites have been hypothesized to facilitate evolutionary differentiation of freshwater fish morphs specializing on different food types. However, as parasites may also respond directly to abiotic environment the parasite risk does not depend only on biotic features of the host environment. It is possible that different morphs experience specific selection gradients by parasites but it is not clear how consistent the selection is when abiotic factors change. We examined parasite pressure in sympatric morphs of threespine stickleback (Gasterosteus aculeatus) across a temperature gradient in two large Icelandic lakes, Myvatn and Thingvallavatn. Habitat-specific temperature gradients in these lakes are opposite. Myvatn lava rock morph lives in a warm environment, while the mud morph lives in the cold. In Thingvallavatn, the lava rock morph lives in a cold environment and the mud morph in a warm habitat. We found more parasites in fish living in higher temperature in both lakes, independent of the fish morph, and this pattern was similar for the two dominating parasite taxa, trematodes and cestodes. However, at the same time, we also found higher parasite abundance in a third morph living in deep cold-water habitat in Thingvallavatn compared to the cold-water lava morph, indicating strong effect of habitat-specific biotic factors. Our results suggest complex interactions between water temperature and biotic factors in determining the parasite community structure, a pattern that may have implications for differentiation of stickleback morphs.

In nature, parasites commonly share hosts with other conspecific parasite genotypes. While adult parasites typically show aggregated distribution in their final hosts, aggregation of clonal parasite genotypes in intermediate hosts, such as those of trematodes in molluscs, is not generally known. However, infection of a host by multiple parasite genotypes has significant implications for evolution of virulence and host-parasite coevolution. Aggregated distribution of the clonal stages can increase host mortality and reduce larval output of each infecting genotype through interclonal competition, and therefore have significant implications for parasite epidemiology. The aim of this study was (i) to find out how common multiple genotype infections (MGIs) are in aquatic snails serving as intermediate hosts for different trematode species; (ii) to find out whether the prevalence of infection could be used to predict MGI frequencies and (iii) to use the relationship to infer whether MGIs aggregate in molluscan hosts. We determined the prevalence of trematode (Diplostomum pseudospathaceum) infections and the frequency of MGIs in snail (Lymnaea stagnalis) host populations as well as compiled corresponding literature data from a range of snail-trematode systems. We used Bayesian simulations to explore the relationship between prevalence of infection and MGI frequency in these data, and tested whether genotypes aggregate in snails by comparing the simulated relations with null model (Poisson and demographic Poisson) expectations. Our results show that MGIs are common in aquatic snails with up to 90% of the infected snails carrying MGIs. Parasite prevalence is a good predictor of MGI frequencies at a range of prevailing natural prevalences of infection (0-50%). The frequency of MGIs was higher than expected by both null models, indicating parasite aggregation at genotype level. These findings are in sharp contrast with the absence of multiple infections in snails at level of trematode species, suggesting that co-infections by multiple species and multiple genotypes of one species are controlled by different biological processes. Aggregation of MGIs in snail hosts appears to be widespread across different snail-trematode systems., (© 2012 The Authors. Journal of Animal Ecology © 2012 British Ecological Society.)

Understanding genetic specificity in factors determining the outcome of host-parasite interactions is especially important as it contributes to parasite epidemiology, virulence, and maintenance of genetic variation. Such specificity, however, is still generally poorly understood. We examined genetic specificity in interactions among coinfecting parasites. In natural populations, individual hosts are often simultaneously infected by multiple parasite species and genotypes that interact. Such interactions could maintain genetic variation in parasite populations if they are genetically specific so that the relative fitness of parasite genotypes varies across host individuals depending on (1) the presence/absence of coinfections and/or (2) the genetic composition of the coinfecting parasite community. We tested these predictions using clones of fish eye flukes Diplostomum pseudospathaceum and Diplostomum gasterostei. We found that interactions among parasites had a strong genetic basis and that this modified genetic variation in infection success of D. pseudospathaceum between single and multiple infections as well as across multiply infected host individuals depending on the genetic identity of the coinfecting D. gasterostei. The relative magnitude of these effects, however, depended on the exposure dose, suggesting that ecological factors can modify genetic interactions between parasites.

Co-infecting parasite genotypes typically compete for host resources limiting their fitness. The intensity of such competition depends on whether parasites are reproducing in a host, or using it primarily as a transmission vehicle while not multiplying in host tissues (referred to as 'competition hypothesis'). Alternatively, simultaneous attack and co-infection by several parasite genotypes might facilitate parasite infection because such a diverse attack could present an additional challenge to host immune defence (referred to as 'facilitation hypothesis'). We tested the competition hypothesis by comparing the production of transmission stages (cercariae) from snails infected with one or two genotypes of the trematode Diplostomum pseudospathaceum. We found that cercarial production did not differ between the two groups of snails, suggesting lower per genotype production in double infections, and competition for host resources. Second, we tested the facilitation hypothesis by comparing parasite infection success on fishes (proportion of parasites establishing in the host) using cercariae originating from single-infected snails, double-infected snails and artificial mixtures of the single genotypes. In both cases, we found higher infection success when fishes were challenged with two parasite genotypes instead of one, supporting the facilitation hypothesis. Our results suggest that constraints defining the success of multiple genotype infections in parasites with multiple host life cycles include both between-genotype resource competition in the host and performance of host immune defences against a diverse parasite challenge.

Infections by multiple parasite species are common in nature and have important consequences for between-species interactions and coevolutionary dynamics with the host populations. For example, ecological and evolutionary factors underlying the structure of parasite communities determine the range of hosts a parasite can infect and set the basis for both evolution of host defences and parasite virulence, as well as management of diseases. Studies investigating these factors have been facilitated in the recent past by genetic methods, which surmount difficulties of traditional morphological taxonomy in identifying individual parasite species. Here we take a step further and present a novel methodological approach to analyze parasite communities as a whole. We determined the relative frequencies of interspecific SNP alleles by pyrosequencing naturally pooled samples of closely related Diplostomum spp. trematodes infecting eye lenses of freshwater fish. Pyrosequencing allowed us to use naturally pooled community samples (lenses) to increase the sample size and statistical power, without sequencing single parasite specimens. In a case study, we applied this method to analyze the community structure of the eye flukes across fish host species of a freshwater system in Finland. We found that the fish species harboured significantly different parasite communities and that multiple species infections were common. Our study provides an example of how quantitative pyrosequencing can be used to answer evolutionary and ecological questions in natural communities of parasites., (Copyright © 2011 Elsevier B.V. All rights reserved.)

Due to their particular way of life, dispersal of parasites is often mediated by their host's biology. Dispersal distance is relevant for parasites because high degree of dispersal leads to high gene flow, which counters the rate of parasite local adaptation in the host populations. Parasites with complex life cycles need to exploit sequentially more than one host species to complete their life cycle. Most trematode parasites have such complex life cycles involving invertebrate and vertebrate hosts. The spatial scales of invertebrate and vertebrate host populations are often different, which may decrease the probability that the parasite cycles locally in the intermediate host population. We used neutral microsatellite markers to determine genetic structure in Diplostomum pseudospathaceum parasites collected from local populations of freshwater snails (Lymnaea stagnalis). D. pseudospathaceum is a trematode that has two intermediate hosts (snail and fish) and a highly motile definitive host (bird). We found that the parasite population infecting the local snail populations showed no genetic structure over a large geographic range (over 300km). We also did not detect evidence for isolation by distance in the parasite. We conclude that dispersal in the motile definitive host is likely to prevent emergence of local population genetic structure in the parasite. Our results suggest that parasite dispersal in the definitive host may limit local cycling of the parasites in the intermediate host populations., (Copyright © 2010 Elsevier B.V. All rights reserved.)

To quantify the effect of the filter feeder Daphnia on the aggregation of mineral particles, temporal changes in the particle size distribution of inorganic colloids were experimentally determined both in the presence and in the absence of Daphnia in water samples of Lake Brienz, Switzerland, an oligotrophic lake rich in suspended inorganic colloids. The results obtained show that daphnids favour the aggregation of mineral colloids, but only for particle sizes above the Daphnia filter mesh size. However, the number concentration of particles smaller than the Daphnia filter mesh size simultaneously increases in the presence of the filter feeder, suggesting either the break-down of existing aggregates or the aggregation of particles with initial sizes below the measured size range. The density of daphnids in this lake is currently too low to have any significant effect on the fate of inorganic colloidal particles as compared with aggregation due to physical processes of particle collision. However, in more productive water bodies where Daphnia is more abundant, they may play a significant role.