Genetic signatures of divergent selection in European beech (Fagus sylvatica L.) are associated with the variation in temperature and precipitation across its distribution range

Rolling down that mountain: microgeographical adaptive divergence during a fast population expansion along a steep environmental gradient in European beech

Forest tree populations harbour high genetic diversity thanks to large effective population sizes and strong gene flow, allowing them to diversify through adaptation to local environmental pressures within dispersal distance. Many tree populations also experienced historical demographic fluctuations, including spatial population contraction or expansions at various temporal scales, which may constrain their ability to adapt to environmental variations. Our aim is to investigate how recent contraction and expansion events interfere with local adaptation, by studying patterns of adaptive divergence between closely related stands undergoing environmentally contrasted conditions, and having or not recently expanded. To investigate genome-wide signatures of local adaptation while accounting for demography, we analysed divergence in a European beech population by testing pairwise differentiation among four tree stands at ~35k Single Nucleotide Polymorphisms from ~9k genomic regions. We applied three divergence outlier search methods resting on different assumptions and targeting either single SNPs or contiguous genomic regions, while accounting for the effect of population size variations on genetic divergence. We found 27 signals of selective signatures in 19 target regions. Putatively adaptive divergence involved all stand pairs. We retrieved signals both when comparing old-growth stands and recently colonised areas and when comparing stands within the old-growth area. Therefore, adaptive divergence processes have taken place both over short time spans, under strong environmental contrasts, and over short ecological gradients, in populations that have been stable in the long term. This suggests that standing genetic variation supports local, microgeographic divergence processes, which can maintain genetic diversity at the landscape level.

A genome-wide genetic association study reveals SNPs significantly associated with environmental variables and specific leaf area in European beech

European beech is negatively affected by climate change and a further growth decline is predicted for large parts of its distribution range. Despite the importance of this species, little is known about its genetic adaptation and especially the genetic basis of its physiological traits. Here, we used genotyping by sequencing to identify SNPs in 43 German European beech populations growing under different environmental conditions. In total, 28 of these populations were located along a precipitation and temperature gradient in northern Germany, and single tree-based hydraulic and morphological traits were available. We obtained a set of 13,493 high-quality SNPs that were used for environmental and SNP-trait association analysis. In total, 22 SNPs were identified that were significantly associated with environmental variables or specific leaf area (SLA). Several SNPs were located in genes related to stress response. The majority of the significant SNPs were located in non-coding (intergenic and intronic) regions. These may be in linkage disequilibrium with the causative coding or regulatory regions. Our study gives insights into the genetic basis of abiotic adaptation in European beech, and provides genetic resources that can be used in future studies on this species. Besides clear patterns of local adaptation to environmental conditions of the investigated populations, the analyzed morphological and hydraulic traits explained most of the explainable genetic variation. Thus, they could successfully be altered in tree breeding programs, which may help to increase the adaptation of European beech to changing environmental conditions in the future.

Climate-driven mitochondrial selection in lacertid lizards

The mitochondrion, which is an intracellular organelle responsible for most of the energy-producing pathways, can have its genome targeted for climate-driven selection. However, climate-driven mitochondrial selection remains a sparsely studied area in reptiles. Here, we reported the complete mitochondrial genome sequence of a lacertid lizard (Takydromus intermedius) and used mitogenomes from 54 species of lacertid lizards to study their phylogenetic relationships and to identify the mitochondrial genes under positive selection by climate. The length of the complete mitochondrial genome sequence of T. intermedius was 17,713 bp, which was within the range of lengths (17,224-18,943) ever reported for Takydromus species. The arrangement of mitochondrial genes in T. intermedius was the same as in other congeneric species. The 54 lacertid species could be divided into three geographically and climatically different clades. We identified three mitochondrial genes (ATP6, ATP8, and ND3) under positive selection by climate, and found that isothermality, temperature seasonality, precipitation of wettest month, and precipitation seasonality were the most important climatic variables contributing to the gene selection.

First genome-wide data from Italian European beech (Fagus sylvatica L.): Strong and ancient differentiation between Alps and Apennines

The European beech (Fagus sylvatica L.) is one of the most widespread forest trees in Europe whose distribution and intraspecific diversity has been largely shaped by repeated glacial cycles. Previous studies, mainly based on palaeobotanical evidence and a limited set of chloroplast and nuclear genetic markers, highlighted a complex phylogeographic scenario, with southern and western Europe characterized by a rather heterogeneous genetic structure, as a result of recolonization from different glacial refugia. Despite its ecological and economic importance, the genome of this broad-leaved tree has only recently been assembled, and its intra-species genomic diversity is still largely unexplored. Here, we performed whole-genome resequencing of nine Italian beech individuals sampled from two stands located in the Alpine and Apennine mountain ranges. We investigated patterns of genetic diversity at chloroplast, mitochondrial and nuclear genomes and we used chloroplast genomes to reconstruct a temporally-resolved phylogeny. Results allowed us to test European beech differentiation on a whole-genome level and to accurately date their divergence time. Our results showed comparable, relatively high levels of genomic diversity in the two populations and highlighted a clear differentiation at chloroplast, mitochondrial and nuclear genomes. The molecular clock analysis indicated an ancient split between the Alpine and Apennine populations, occurred between the Günz and the Riss glaciations (approximately 660 kyrs ago), suggesting a long history of separation for the two gene pools. This information has important conservation implications in the context of adaptation to ongoing climate changes.

The establishment of plants following long-distance dispersal

Long-distance dispersal (LDD) beyond the range of a species is an important driver of ecological and evolutionary patterns, but insufficient attention has been given to postdispersal establishment. In this review, we summarize current knowledge of the post-LDD establishment phase in plant colonization, identify six key determinants of establishment success, develop a general quantitative framework for post-LDD establishment, and address the major challenges and opportunities in future research. These include improving detection and understanding of LDD using novel approaches, investigating mechanisms determining post-LDD establishment success using mechanistic modeling and inference, and comparison of establishment between past and present. By addressing current knowledge gaps, we aim to further our understanding of how LDD affects plant distributions, and the long-term consequences of LDD events.

Transcriptome-based analyses of adaptive divergence between two closely related spruce species on the Qinghai-Tibet plateau and adjacent regions

Speciation among populations connected by gene flow is driven by adaptation to different environments, but underlying gene-environment associations remain largely unknown. Here, 162 individuals from 32 populations were sampled to obtain 191,648 independent single nucleotide polymorphisms (SNPs) across the genomes of two closely related spruce species, Picea asperata and Picea crassifolia, which occur on the Qinghai-Tibet Plateau and in surrounding regions. Using the SNP data set, genotype-environment associations and demographic modelling were used to examine local adaptation and genetic divergence between these two species. While morphologically similar, the two Picea species were genetically differentiated in multiple analyses. These species diverged despite continuous gene flow, and their initial divergence was dated back to the late Quaternary. The effective population sizes of both species have expanded since their divergence, as confirmed by niche distribution simulations. A total of 6365 genes were associated with the tested environmental variables; of these, 41 were positively selected in P. asperata and were mainly associated with temperature, while 83 were positively selected in P. crassifolia and were primarily associated with precipitation. These results deepen our understanding of the adaptive divergence and demographic histories of these two spruce species and highlight the importance of genomic data in deciphering the environmental selection underlying Quaternary interspecific divergence.

A novel synthesis of two decades of microsatellite studies on European beech reveals decreasing genetic diversity from glacial refugia

Genetic diversity influences the evolutionary potential of forest trees under changing environmental conditions, thus indirectly the ecosystem services that forests provide. European beech (Fagus sylvatica L.) is a dominant European forest tree species that increasingly suffers from climate change-related die-back. Here, we conducted a systematic literature review of neutral genetic diversity in European beech and created a meta-data set of expected heterozygosity (He) from all past studies providing nuclear microsatellite data. We propose a novel approach, based on population genetic theory and a min-max scaling to make past studies comparable. Using a new microsatellite data set with unprecedented geographic coverage and various re-sampling schemes to mimic common sampling biases, we show the potential and limitations of the scaling approach. The scaled meta-dataset reveals the expected trend of decreasing genetic diversity from glacial refugia across the species range and also supports the hypothesis that different lineages met and admixed north of the European mountain ranges. As a result, we present a map of genetic diversity across the range of European beech which could help to identify seed source populations harboring greater diversity and guide sampling strategies for future genome-wide and functional investigations of genetic variation. Our approach illustrates how to combine information from several nuclear microsatellite data sets to describe patterns of genetic diversity extending beyond the geographic scale or mean number of loci used in each individual study, and thus is a proof-of-concept for synthesizing knowledge from existing studies also in other species.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11295-022-01577-4.

Integrating transcriptomic and metabolomic analysis in roots of wild soybean seedlings in response to low-phosphorus stress

INTRODUCTION

Plants undergo divergent adaptations to form different ecotypes when exposed to different habitats. Ecotypes with ecological adaptation advantages are excellent germplasm resources for crop improvement.

METHODS

his study comprehensively compared the differences in morphology and physiological mechanisms in the roots of two different ecotypes of wild soybean (Glycine soja) seedlings under artificially simulated low-phosphorus (LP) stress.

RESULT

The seedlings of barren-tolerant wild soybean (GS2) suffered less damage than common wild soybean (GS1). GS2 absorbed more phosphorus (P) by increasing root length. In-depth integrated analyses of transcriptomics and metabolomics revealed the formation process of the ecological adaptability of the two different ecotypes wild soybean from the perspective of gene expression and metabolic changes. This study revealed the adaptation process of GS2 from the perspective of the adaptation of structural and molecular metabolism, mainly including: (1) Enhancing the metabolism of phenolic compounds, lignin, and organic acid metabolism could activate unavailable soil P; (2) Up-regulating genes encoding pectinesterase and phospholipase C (PLC) specifically could promote the reuse of structural P; (3) Some factors could reduce the oxidative damage to the membranes caused by LP stress, such as accumulating the metabolites putrescine and ascorbate significantly, up-regulating the genes encoding SQD2 (the key enzyme of sulfolipid substitution of phospholipids) substantially and enhancing the synthesis of secondary antioxidant metabolite anthocyanins and the AsA-GSH cycle; (4) enhancing the uptake of soil P by upregulating inorganic phosphate transporter, acid phosphatase ACP1, and purple acid phosphatase genes; (5) HSFA6b and MYB61 are the key TFs to resist LP stress.

DISCUSSION

In general, GS2 could resist LP stress by activating unavailable soil P, reusing plant structural P, rebuilding membrane lipids, and enhancing the antioxidant membrane protection system. Our study provides a new perspective for the study of divergent adaptation of plants.

A set of nuclear SNP loci derived from single sample double digest RAD and from pool sequencing for large-scale genetic studies in the European beech Fagus sylvatica

The large-scale spatial genetic structure of European beech, Fagus sylvatica, has been until now poorly studied. We conducted double digest RAD sequencing (ddRADseq) on 54 beech individuals stemming from 36 provenances to discover spatially informative nuclear SNP loci. In addition, two pools derived from 14 early and 14 late flushing individuals each were sequenced with Illumina HiSeq. From an initial amount of 5,464 loci detected by ddRADseq, we selected 559 informative loci. Further 27 additional loci showing significant allelic differences among early and late flushing individuals could be identified after a genotyping on 95 test individuals. The final selection of 578 loci was submitted to probe design for targeted genotyping by sequencing, which yielded 543 loci. The new set of SNP loci should be, after validation on a larger sample size, useful for large-scale genetic studies in this economically-important species.