Advancing research on ectomycorrhizal fungal adaptation with landscape genomics.

Ectomycorrhizal (ECM) fungi serve key functions in forest ecosystems by supplying water and nutrients to tree hosts, yet mutualistic plant-fungi interactions are jeopardised by environmental alterations. Here, we discuss the great potential and current limitations of landscape genomics in investigating signatures of local adaptation in natural populations of ECM fungi.

Re-thinking the environment in landscape genomics.

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.

The fall of the summer truffle: Recurring hot, dry summers result in declining fruitbody production of Tuber aestivum in Central Europe.

Global warming is pushing populations outside their range of physiological tolerance. According to the environmental envelope framework, the most vulnerable populations occur near the climatic edge of their species' distributions. In contrast, populations from the climatic center of the species range should be relatively buffered against climate warming. We tested this latter prediction using a combination of linear mixed effects and machine learning algorithms on an extensive, citizen-scientist generated dataset on the fruitbody productivity of the Burgundy (aka summer) truffle (Tuber aestivum Vittad.), a keystone, ectomycorrhizal tree-symbiont occurring on a wide range of temperate climates. T. aestivum's fruitbody productivity was monitored at 3-week resolution over up to 8 continuous years at 20 sites distributed in the climatic center of its European distribution in southwest Germany and Switzerland. We found that T. aestivum fruitbody production is more sensitive to summer drought than would be expected from the breadth of its species' climatic niche. The monitored populations occurring nearly 5°C colder than the edge of their species' climatic distribution. However, interannual fruitbody productivity (truffle mass year-1 ) fell by a median loss of 22% for every 1°C increase in summer temperature over a site's 30-year mean. Among the most productive monitored populations, the temperature sensitivity was even higher, with single summer temperature anomalies of 3°C sufficient to stop fruitbody production altogether. Interannual truffle productivity was also related to the phenology of host trees, with ~22 g less truffle mass for each 1-day reduction in the length of the tree growing season. Increasing summer drought extremes are therefore likely to reduce fruiting among summer truffle populations throughout Central Europe. Our results suggest that European T. aestivum may be a mosaic of vulnerable populations, sensitive to climate-driven declines at lower thresholds than implied by its species distribution model.

Genetic basis of growth reaction to drought stress differs in contrasting high-latitude treeline ecotones of a widespread conifer.

Climate change is increasing the frequency and intensity of drought events in many boreal forests. Trees are sessile organisms with a long generation time, which makes them vulnerable to fast climate change and hinders fast adaptations. Therefore, it is important to know how forests cope with drought stress and to explore the genetic basis of these reactions. We investigated three natural populations of white spruce (Picea glauca) in Alaska, located at one drought-limited and two cold-limited treelines with a paired plot design of one forest and one treeline plot. We obtained individual increment cores from 458 trees and climate data to assess dendrophenotypes, in particular the growth reaction to drought stress. To explore the genetic basis of these dendrophenotypes, we genotyped the individual trees at 3000 single nucleotide polymorphisms in candidate genes and performed genotype-phenotype association analysis using linear mixed models and Bayesian sparse linear mixed models. Growth reaction to drought stress differed in contrasting treeline populations. Therefore, the populations are likely to be unevenly affected by climate change. We identified 40 genes associated with dendrophenotypic traits that differed among the treeline populations. Most genes were identified in the drought-limited site, indicating comparatively strong selection pressure of drought-tolerant phenotypes. Contrasting patterns of drought-associated genes among sampled sites and in comparison to Canadian populations in a previous study suggest that drought adaptation acts on a local scale. Our results highlight genes that are associated with wood traits which in turn are critical for the establishment and persistence of future forests under climate change.

Hidden fairy rings and males-Genetic patterns of natural Burgundy truffle (Tuber aestivum Vittad.) populations reveal new insights into its life cycle.

Burgundy truffles are heterothallic ascomycetes that grow in symbiosis with trees. Despite their esteemed belowground fruitbodies, the species' complex lifecycle is still not fully understood. Here, we present the genetic patterns in three natural Burgundy truffle populations based on genotyped fruitbodies, ascospore extracts and ectomycorrhizal root tips using microsatellites and the mating-type locus. Distinct genetic structures with high relatedness in close vicinity were found for females (forming the fruitbodies) and males (fertilizing partner as inferred from ascospore extracts), with high genotypic diversity and annual turnover of males, suggesting that ephemeral male mating partners are germinating ascospores from decaying fruitbodies. The presence of hermaphrodites and the interannual persistence of a few males suggest that persistent mycelia may sporadically also act as males. Only female or hermaphroditic individuals were detected on root tips. At one site, fruitbodies grew in a fairy ring formed by a large female individual that showed an outward growth rate of 30 cm per year, with the mycelium decaying within the ring and being fertilized by over 50 male individuals. While fairy ring structures have never been shown for truffles, the genetics of Burgundy truffle populations support a similar reproductive biology as those of other highly prized truffles.

Transcriptome-wide SNPs for Botrychium lunaria ferns enable fine-grained analysis of ploidy and population structure.

Ferns are the second most diverse group of land plants after angiosperms. Extant species occupy a wide range of habitats and contribute significantly to ecosystem functioning. Despite the importance of ferns, most taxa are poorly covered by genomic resources and within-species studies based on high-resolution markers are entirely lacking. The genus Botrychium belongs to the family Ophioglossaceae, which includes species with very large genomes and chromosome numbers (e.g., Ophioglossum reticulatum 2n = 1520). The genus has a cosmopolitan distribution with 35 species, half of which are polyploids. Here, we establish a transcriptome for Botrychium lunaria (L.) Sw., a diploid species with an extremely large genome of about ~19.0-23.7 Gb. We assembled 25,677 high-quality transcripts with an average length of 1,333 bp based on deep RNA-sequencing of a single individual. We sequenced 11 additional transcriptomes of individuals from two populations in Switzerland, including the population of the reference individual. Based on read mapping to reference transcript sequences, we identified 374,463 single nucleotide polymorphisms (SNPs) segregating among individuals for an average density of 14 SNPs per kilobase. We found that all 12 transcriptomes were most likely from diploid individuals. The transcriptome-wide markers provided unprecedented resolution of the population genetic structure, revealing substantial variation in heterozygosity among individuals. We also constructed a phylogenomic tree of 92 taxa representing all fern orders to ascertain the placement of the genus Botrychium. High-quality transcriptomic resources and SNP sets constitute powerful population genomic resources to investigate the ecology, and evolution of fern populations.

A dormant resource for genome size estimation in ferns: C-value inference of the Ophioglossaceae using herbarium specimen spores.

PREMISE: The great variation of genome size (C-value) across land plants is linked to various adaptative features. Flow cytometry (FCM), the standard approach to estimating C-values, relies mostly on fresh materials, performing poorly when used with herbarium materials. No fern C-value reports have been derived from herbarium specimens; however, the herbarium spores of some ferns remain highly viable for decades and are thus promising for further investigation. To explore this possibility, we evaluated herbarium spore collections of Ophioglossaceae ferns using FCM. METHODS: Flow cytometry was conducted on 24 spore samples, representing eight of the 12 genera of the Ophioglossaceae, using specimens ranging in age from 2.6 to 111 years obtained from five herbaria. RESULTS: Regardless of the genus or the source herbarium, high-quality C-value data were generated from 17 samples, with the oldest being 26 years old. Estimates of the C-values from sporophytic tissues of known ploidy did not reveal any evidence of apomixis for the species surveyed here. We also detected a pronounced genome downsizing in Sceptridium polyploids. DISCUSSION: The recent success of FCM for C-value estimation using spores provides a much more convenient method of utilizing "dry" refrigerated materials. We demonstrate here that herbarium spores of some ferns are also promising for this use, even for older specimens.

Cryptic genetic structure and copy-number variation in the ubiquitous forest symbiotic fungus Cenococcum geophilum.

Ectomycorrhizal (ECM) fungi associated with plants constitute one of the most successful symbiotic interactions in forest ecosystems. ECM support trophic exchanges with host plants and are important factors for the survival and stress resilience of trees. However, ECM clades often harbour morpho-species and cryptic lineages, with weak morphological differentiation. How this relates to intraspecific genome variability and ecological functioning is poorly known. Here, we analysed 16 European isolates of the ascomycete Cenococcum geophilum, an extremely ubiquitous forest symbiotic fungus with no known sexual or asexual spore-forming structures but with a massively enlarged genome. We carried out whole-genome sequencing to identify single-nucleotide polymorphisms. We found no geographic structure at the European scale but divergent lineages within sampling sites. Evidence for recombination was restricted to specific cryptic lineages. Lineage differentiation was supported by extensive copy-number variation. Finally, we confirmed heterothallism with a single MAT1 idiomorph per genome. Synteny analyses of the MAT1 locus revealed substantial rearrangements and a pseudogene of the opposite MAT1 idiomorph. Our study provides the first evidence for substantial genome-wide structural variation, lineage-specific recombination and low continent-wide genetic differentiation in C. geophilum. Our study provides a foundation for targeted analyses of intra-specific functional variation in this major symbiosis.

Prospects and limitations of genomic offset in conservation management.

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.

The GenTree Platform: growth traits and tree-level environmental data in 12 European forest tree species.

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.

Genomic vulnerability to rapid climate warming in a tree species with a long generation time.

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.

Disentangling the effects of geographic peripherality and habitat suitability on neutral and adaptive genetic variation in Swiss stone pine.

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, FST 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 (He ) 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.

Author Correction: The GenTree Dendroecological Collection, tree-ring and wood density data from seven tree species across Europe.

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

Biotic and abiotic factors shape arbuscular mycorrhizal fungal communities associated with the roots of the widespread fern Botrychium lunaria (Ophioglossaceae).

Arbuscular mycorrhizal fungi (AMF) play central roles in terrestrial ecosystems by interacting with both above and belowground communities as well as by influencing edaphic properties. The AMF communities associated with the roots of the fern Botrychium lunaria (Ophioglossaceae) were sampled in four transects at 2400 m a.s.l. in the Swiss Alps and analyzed using metabarcoding. Members of five Glomeromycota genera were identified across the 71 samples. Our analyses revealed the existence of a core microbiome composed of four abundant Glomus operational taxonomic units (OTUs), as well as a low OTU turnover between samples. The AMF communities were not spatially structured, which contrasts with most studies on AMF associated with angiosperms. pH, microbial connectivity and humus cover significantly shaped AMF beta diversity but only explained a minor fraction of variation in beta diversity. AMF OTUs associations were found to be significant by both cohesion and co-occurrence analyses, suggesting a role for fungus-fungus interactions in AMF community assembly. In particular, OTU co-occurrences were more frequent between different genera than among the same genus, rising the hypothesis of functional complementarity among the AMF associated to B. lunaria. Altogether, our results provide new insights into the ecology of fern symbionts in alpine grasslands.

The GenTree Dendroecological Collection, tree-ring and wood density data from seven tree species across Europe.

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.

Using transcriptome sequencing and pooled exome capture to study local adaptation in the giga-genome of Pinus cembra.

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.

Rapid allopolyploid radiation of moonwort ferns (Botrychium; Ophioglossaceae) revealed by PacBio sequencing of homologous and homeologous nuclear regions.

Polyploidy is a major speciation process in vascular plants, and is postulated to be particularly important in shaping the diversity of extant ferns. However, limitations in the availability of bi-parental markers for ferns have greatly limited phylogenetic investigation of polyploidy in this group. With a large number of allopolyploid species, the genus Botrychium is a classic example in ferns where recurrent polyploidy is postulated to have driven frequent speciation events. Here, we use PacBio sequencing and the PURC bioinformatics pipeline to capture all homeologous or allelic copies of four long (∼1 kb) low-copy nuclear regions from a sample of 45 specimens (25 diploids and 20 polyploids) representing 37 Botrychium taxa, and three outgroups. This sample includes most currently recognized Botrychium species in Europe and North America, and the majority of our specimens were genotyped with co-dominant nuclear allozymes to ensure species identification. We analyzed the sequence data using maximum likelihood (ML) and Bayesian inference (BI) concatenated-data ("gene tree") approaches to explore the relationships among Botrychium species. Finally, we estimated divergence times among Botrychium lineages and inferred the multi-labeled polyploid species tree showing the origins of the polyploid taxa, and their relationships to each other and to their diploid progenitors. We found strong support for the monophyly of the major lineages within Botrychium and identified most of the parental donors of the polyploids; these results largely corroborate earlier morphological and allozyme-based investigations. Each polyploid had at least two distinct homeologs, indicating that all sampled polyploids are likely allopolyploids (rather than autopolyploids). Our divergence-time analyses revealed that these allopolyploid lineages originated recently-within the last two million years-and thus that the genus has undergone a recent radiation, correlated with multiple independent allopolyploidizations across the phylogeny. Also, we found strong parental biases in the formation of allopolyploids, with individual diploid species participating multiple times as either the maternal or paternal donor (but not both). Finally, we discuss the role of polyploidy in the evolutionary history of Botrychium and the interspecific reproductive barriers possibly involved in these parental biases.

Molecular phylogenetics supports widespread cryptic species in moonworts (Botrychium s.s., Ophioglossaceae).

PREMISE OF THE STUDY: Previous phylogenetic studies of moonworts (Botrychium sensu stricto (s.s.)) included few taxa from outside of North America. This low geographical representation limited interpretations of relationships of this group rich in cryptic species. With 18 out of 30 species in the genus being polyploid, understanding their evolutionary history remains a major challenge. METHODS: A new molecular phylogeny was reconstructed using Maximum Likelihood (ML) and Bayesian Inference (BI) analyses based on multiple accessions of the most wide-ranging Arctic taxa of Botrychium in North America and Europe using three noncoding plastid DNA regions (psbA-trnH(GUG), trnL(UAA)-trnF(GAA) intergenic spacer, and rpL16 intron). KEY RESULTS: The new phylogeny confirms the identity of several recently described species and proposed new taxa. Nine subclades are newly identified within the two major clades in Botrychium s.s.: Lanceolatum and Lunaria. Chloroplast DNA was variable enough to separate morphologically cryptic species in the Lunaria clade. On the contrary, much less variation is seen within the morphologically variable Lanceolatum clade despite sampling over the same broad geographic range. The chloroplast region psbA-trnH(GUG) is identified as an efficient DNA barcode for the identification of cryptic taxa in Botrychium s.s. CONCLUSIONS: The combined increase in species representation, samples from throughout the geographic range of each species, and sequencing of multiple plastid DNA regions supports morphologically cryptic species in Botrychium s.s.