Search

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.

[Motivation] Trait variation within species can reveal plastic and/or genetic responses to environmental gradients, and may indicate where local adaptation has occurred. Here, we present a dataset of rangewide variation in leaf traits from seven of the most ecologically and economically important tree species in Europe. Sample collection and trait assessment are embedded in the GenTree project (EU-Horizon 2020), which aims at characterizing the genetic and phenotypic variability of forest tree species to optimize the management and sustainable use of forest genetic resources. Our dataset captures substantial intra- and interspecific leaf phenotypic variability, and provides valuable information for studying the relationship between ecosystem functioning and trait variability of individuals, and the response and resilience of species to environmental changes., [Main types of variable contained] We chose morphological and chemical characters linked to trade-offs between acquisition and conservation of resources and water use, namely specific leaf area, leaf size, carbon and nitrogen content and their ratio, and the isotopic signature of stable isotope 13C and 15N in leaves., [Spatial location and grain] We surveyed between 18 and 22 populations per species, 141 in total, across Europe., [Time period] Leaf sampling took place between 2016 and 2017., [Major taxa and level of measurement] We sampled at least 25 individuals in each population, 3,569 trees in total, and measured traits in 35,755 leaves from seven European tree species, i.e. the conifers Picea abies, Pinus pinaster and Pinus sylvestris, and the broadleaves Betula pendula, Fagus sylvatica, Populus nigra and Quercus petraea., [Software format] The data files are in ASCII text, tab delimited, not compressed., All authors acknowledge support from the European Union’s Horizon 2020 Research and Innovation Programme under grant agreement no. 676876 (Project GenTree). Laboratory work was also supported by other projects: COMEDIAS (CGL2017- 83170-R, Spanish Ministry of Science, Innovation and Universities), REMEDINAL IV (TE-CM S2018/EMT-4338, Project of the Autonomous Community of Madrid, Spain), SPONFOREST (BIODIVERSA-3, PCIN-2016-055), and International Laboratory of Global Change (LINCGlobal 4540-143AP). The Swiss contribution was supported by the Swiss Secretariat for Education, Research and Innovation (SERI) under contract no. 6.0032. BC was funded by a Scholarship from Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - CAPES (Brazil).

Premise of research. Vascular plants depend on sexual recombination for generating new genetic variability to meet environmental needs. Nevertheless, members of the early-divergent fern genus Botry...

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 similar to 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., Molecular Ecology Resources, 22 (1), ISSN:1755-098X, ISSN:1471-8278, ISSN:1755-0998, ISSN:1471-8286

SummaryEctomycorrhizal (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.Originality-Significance StatementWe provide the first report on the genetic structure and copy-number variation of the globally ubiquitous and key forest symbiotic fungus Cenococcum geophilum using whole-genome sequencing data. We found divergent lineages within sampling sites, while closely related lineages appear over large geographic distances on a continental scale. Even though no sexual spore forming structures have been reported to date, we provide evidence of recombination in a specific lineage suggesting mating activity. Our findings help explain the high genetic diversity occurring within populations and their resilience to changing and adverse environmental conditions. Furthermore, we identify a single MAT1 idiomorph per genome, confirming heterothallism, and discover that major genomic rearrangements are found in their flanking regions based on chromosomal synteny analysis. Intriguingly, a pseudogene of the opposite functional idiomorph has been characterised in each genome, suggesting a common homothallic ancestor to the species. As Cenococcum geophilum is a pivotal mycorrhizal associate of a broad range of trees and shrubs providing nutrition and water supply to their hosts, we highlight and discuss the potential role of the large genome-wide structural variations in environmental selection.

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.

Ferns are the second most dominant group of land plants after angiosperms. Extant species occupy an extensive range of habitats and contribute significantly to ecosystem functioning. Despite the importance of ferns, most taxa are poorly covered by genomic resources. The genusBotrychiumbelongs to the family Ophioglossaceae, one of the earliest divergent lineages of vascular plants, and has a cosmopolitan distribution with 35 species, half of which are polyploids. Here, we establish a transcriptome forBotrychium lunaria, a diploid species with an extremely large genome with a 1C value of 12.10 pg. We assembled 25,701 high-quality transcripts with an average length of 1,332 bp based on deep RNA-sequencing of a single individual. We sequenced an additional 11 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,510 single nucleotide polymorphisms (SNPs) segregating among individuals for an average density of 14 SNPs per kb. 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 90 taxa representing all fern orders to ascertain the placement of the genusBotrychium. The high-quality transcriptomic resources enable powerful population and phylogenomic studies in an important group of ferns.Significance statementFerns pose substantial puzzles in terms of lifestyles, genome organization and population structure. Progress has been significantly hampered by the lack of genomic resources. Here, we present a transcriptome forBotrychium lunaria, a phylogenetically early-branching fern with an extremely large genome. We show that the new transcriptome improves phylogenetic resolution among early-branching ferns. Based on an additional 11 transcriptomes of the same species, we identify unexpected variation in population-level heterozygosity.

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. Measurement(s) growth ring • wood density Technology Type(s) measuring table • calculation Factor Type(s) tree species Sample Characteristic - Organism Betula pendula • Fagus sylvatica • Picea abies • Populus nigra • Pinus pinaster • Pinus sylvestris • Quercus petraea Sample Characteristic - Location Europe Published

The moonwort genus, Botrychium s. s., includes diploid and polyploid taxa that occur primarily in the northern hemisphere. Their evolutionary history, morphologically cryptic taxa and deep divergence of the family in the phylogeny of ferns has long fascinated pteridologists. Previous molecular studies did not include a complete taxonomic sampling of the taxa in the genus, nor multiple specimens from throughout the known geographical range of each taxon. Therefore, to investigate evolutionary relationships of the major clades of Botrychium s. s., we increased both taxonomic representativeness (multiple accessions per taxa), as well as phylogenetic resolution by including additional new chloroplast markers. To confirm identification and provide evidence from both maternal and paternal parentage of allopolyploids, we also included specimens that have been characterized by allozyme profiles determined by electrophoretic analysis of 20 nuclear enzyme loci for each taxon. We analyzed four chloroplast regions (matK intron, trnHGUG —psbA, andtrnLUAA —trnFGAA intergenic spacers, and rpL16 intron region) of 365 specimens fromAsia, Europe, North America, Oceania, and South America, sampling the geographical range of 34 of 35 accepted Botrychium s. s. taxa and thirteen putatively new taxa. We conducted a phylogenetic analysis of maternal lineages based on 2,385 aligned nucleotides using maximum likelihood and Bayesian inference to explore genetic diversity and phylogenetic relationships among taxa. We found strong support for themonophyly of three major clades: Lanceolatum, Lunaria, and Simplex-Campestre, and resolved 15 subclades. Our results suggest multiple origins for at least four polyploid taxa (B. boreale, B. michiganense, B. yaaxudakeit, and B. watertonense). The Simplex-Campestre clade had the largest number of species, despite having a similar total number of haplotypes as the Lunaria clade (62 and 59, respectively), which has the broadest worldwide distribution. In total, our new molecular phylogeny comprises 47 taxa, of which thirteen are discussed for possible taxonomic recognition.

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.

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.

Recent cytological and molecular studies have investigated genome size variation and evolution in the homosporous ferns, but representatives of the Ophioglossaceae were largely overlooked, despite their evolutionary singularity. Flow cytometry analysis was performed on 41 individuals of eight species of the genera Botrychium (B.), Botrypus, and Sceptridium to estimate their ploidy level variation. In a subset of individuals, we also estimated the absolute genome size and corresponding C values. Additionally, a classical chromosome count was made on the recently described species B. alaskense. Ploidy level and new genome size records were determined for Botrychium alaskense, B. boreale, B. lanceolatum, B. “neolunaria” ined., B. pinnatum, Botrypus virginianus and Sceptridium multifidum. In addition, we confirmed the genome size of B. matricariifolium, B. minganense and B. lunaria. Two of the three major sub-clades of Botrychium differ slightly in their averaged homoploid genome size (subclade Lanceolatum, 24.72 ± 0.40 pg; subclade Lunaria, 27.51 ± 0.47 pg). Flow cytometry and chromosome counting confirmed that B. alaskense is a tetraploid. A new hexaploid cytotype, putatively formed through an autopolyploidization from the sympatric tetraploid cytotype, was detected in a single individual of B. boreale. This is only the second report of hexaploidy in the genus Botrychium and our data highlight the potential to find other ploidy levels within other Botrychium species. Interestingly, no difference within the monoploid genome sizes was detected between ploidy levels, thus supporting the hypothesis of genome size stability after polyploidization and rejecting the scenario of genome downsizing.

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 > 140,000 SNPs in around 13,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 ≥ 40x, i.e. the 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. Using population genetic and landscape genomic methods, we show that inferences of neutral and adaptive genetic variation may be biased when not accounting for such multi-copy genes. Future studies should therefore put more emphasis on identifying paralogous loci, which will be facilitated by the establishment of additional high-quality reference genomes.

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.

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.