Search

Abstract 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.

Abstract Local adaptation is key for ecotypic differentiation and species evolution. Understanding underlying genomic patterns can allow the prediction of future maladaptation and ecosystem stability. Here, we report the whole-genome resequencing of 874 individuals from 100 range-wide populations of European beech (Fagus sylvatica L.), an important forest tree species in Europe. We show that genetic variation closely mirrors geography with a clear pattern of isolation-by-distance. Genome-wide analyses for genotype-environment associations (GEAs) identify relatively few potentially adaptive variants after correcting for an overwhelming signal of statistically significant but non-causal GEAs. We characterize the single high confidence genomic region and pinpoint a candidate gene possibly involved in winter temperature adaptation via modulation of spring phenology. Surprisingly, allelic variation at this locus does not result in any apparent fitness differences in a common garden. More generally, reciprocal transplant experiments across large climate distances suggest extensive phenotypic plasticity. Nevertheless, we find indications of polygenic adaptation which may be essential in natural ecosystems. This polygenic signal exhibits broad- and fine-scale variation across the landscape, highlighting the relevance of spatial resolution. In summary, our results emphasize the importance, but also exemplify the complexity, of employing natural genetic variation for forest conservation under climate change.

Abstract Effective population size (Ne) is a pivotal evolutionary parameter with crucial implications in conservation practice and policy. Genetic methods to estimate Ne have been preferred over demographic methods because they rely on genetic data rather than time‐consuming ecological monitoring. Methods based on linkage disequilibrium (LD), in particular, have become popular in conservation as they require a single sampling and provide estimates that refer to recent generations. A software program based on the LD method, GONE, looks particularly promising to estimate contemporary and recent‐historical Ne (up to 200 generations in the past). Genomic datasets from non‐model species, especially plants, may present some constraints to the use of GONE, as linkage maps and reference genomes are seldom available, and SNP genotyping is usually based on reduced‐representation methods. In this study, we use empirical datasets from four plant species to explore the limitations of plant genomic datasets when estimating Ne using the algorithm implemented in GONE, in addition to exploring some typical biological limitations that may affect Ne estimation using the LD method, such as the occurrence of population structure. We show how accuracy and precision of Ne estimates potentially change with the following factors: occurrence of missing data, limited number of SNPs/individuals sampled, and lack of information about the location of SNPs on chromosomes, with the latter producing a significant bias, previously unexplored with empirical data. We finally compare the Ne estimates obtained with GONE for the last generations with the contemporary Ne estimates obtained with the programs currentNe and NeEstimator.

Abstract Forest ecosystems are increasingly challenged by extreme events, for example, drought, storms, pest attacks, and fungal pathogen outbreaks, causing severe ecological and economic losses. Understanding the genetic basis of adaptive traits in tree species is of key importance to preserve forest ecosystems, as genetic variation in a trait (i.e., heritability) determines its potential for human‐mediated or evolutionary change. Maritime pine (Pinus pinaster Aiton), a conifer widely distributed in southwestern Europe and northwestern Africa, grows under contrasted environmental conditions promoting local adaptation. Genetic variation at adaptive phenotypes, including height, spring phenology, and susceptibility to two fungal pathogens (Diplodia sapinea and Armillaria ostoyae) and an insect pest (Thaumetopoea pityocampa), was assessed in a range‐wide clonal common garden of maritime pine. Broad‐sense heritability was significant for height (0.219), spring phenology (0.165–0.310), and pathogen susceptibility (necrosis length caused by D. sapinea, 0.152; and by A. ostoyae, 0.021, measured on inoculated, excised branches under controlled conditions), but not for pine processionary moth incidence in the common garden. The correlations of trait variation among populations revealed contrasting trends for pathogen susceptibility to D. sapinea and A. ostoyae with respect to height. Taller trees showed longer necrosis length caused by D. sapinea while shorter trees were more affected by A. ostoyae. Moreover, maritime pine populations from areas with high summer temperatures and frequent droughts were less susceptible to D. sapinea but more susceptible to A. ostoyae. Finally, an association study using 4227 genome‐wide SNPs revealed several loci significantly associated with each trait (range of 3–26), including a possibly disease‐induced translation initiation factor, eIF‐5, associated with needle discoloration caused by D. sapinea. This study provides important insights to develop genetic conservation and breeding strategies integrating species responses to biotic stressors.

Marginal and peripheral populations are important for biodiversity conservation. Their original situation in a species’ geographic and ecological space often confers them genetic diversity and traits of high adaptive value. Yet theoretical hypotheses related to marginality are difficult to test because of confounding factors that influence marginality, namely environment, geography, and history. There is an urgent need to develop metrics to disentangle these confounding factors. We designed nine quantitative indices of marginality and peripherality that define where margins lie within species distributions, from a geographical, an environmental and a historical perspective. Using the distribution maps of eight European forest tree species, we assessed whether these indices were idiosyncratic or whether they conveyed redundant information. Using a database on marginal and peripheral populations based on expert knowledge, we assessed the capacity of the indices to predict the marginality status of a population. There was no consistent pattern of correlation between indices across species, confirming that the indices conveyed different information related to the specific geometry of the species distributions. Contrasting with this heterogeneity of correlation patterns across species, the relative importance of the indices to predict the marginality status of populations was consistent across species. However, there was still a significant country effect in the marginality status, showing a variation in expert opinion of marginality vis-á-vis the species distribution. The marginality indices that we developed are entirely based on distribution maps and can be used for any species. They pave the way for testing hypotheses related to marginality and peripherality, with important implications in quantitative ecology, genetics, and biodiversity conservation.

Abstract European forests are expanding and becoming denser following the widespread abandonment of farmland and rural areas. Spontaneous forest regrowth provides a cost‐effective opportunity to restore ecosystems, enhance multifunctionality and sustainability and mitigate climate change. Yet, little is known about the goods and services that such forests provide to people. We assessed the changes in nature's contributions to people (NCP) from spontaneous forest regrowth, i.e. forest expansion and densification, in South‐West Europe. We investigated 65 forest plots in four different landscapes with contrasting ecological and societal contexts. Two landscapes are located in rural areas undergoing human exodus and forest expansion and densification; the other two, in peri‐urban areas with intense land use and forest densification but negligible expansion. For each forest plot, we estimated variables related to ten out of the 18 main NCP defined by the Intergovernmental Science‐Policy Platform on Biodiversity and Ecosystem Services (IPBES). Regulating and material NCP were addressed using variables measured in the field as proxies. Non‐material NCP were studied through stakeholder interviews. Our results show across the cases that forest expansion and densification are generally associated with greater climate regulation and energy provision. Changes in other NCP, especially in non‐material ones, were strongly context‐dependent. The social perception of spontaneous forest regrowth was primarily negative in rural areas and more positive in peri‐urban landscapes. Passive restoration through spontaneous forest expansion and densification can enhance regulating and material NCP, especially when adaptive management is applied. To optimise NCP and to increase the societal awareness of and interest in spontaneous forest regrowth, the effects of this process should be analysed in close coordination with local stakeholders to unveil and quantify the many and complex trade‐offs involved in rural or peri‐urban social perceptions. A free Plain Language Summary can be found within the Supporting Information of this article.

Abstract Detecting the molecular basis of local adaptation and identifying selective drivers is still challenging in nonmodel species. The use of purely population genetic approaches is limited by some characteristics of genetic systems, such as pleiotropy and polygenic control, and parallel evidence from phenotypic‐based experimental comparisons is required. In long‐lived organisms, the detection of selective pressures might also be precluded by evolutionary lag times in response to the environment. Here, we used the English yew to showcase an example of a multiscale integrative approach in a nonmodel species with limited plant and genomic resources. We combined information from two independent sources, phenotypes in a common environment and genomic data in natural populations, to investigate the signature of selection. Growth differences among populations in a common environment, and phenological patterns of both shoot elongation and male strobili maturation, were associated with climate clines, providing evidence for local adaptation and guiding us in the selection of populations for genomic analyses. We used information on over 25,000 SNPs from c. 1,200 genes to infer the demographic history and to test for molecular signatures of selection at different levels: SNP, gene, and biological pathway. Our results confirmed an overall demographic history of population decline, but we also found evidence for putative local adaptation at the molecular level. We identified or confirmed several candidate genes for positive and negative selection in forest trees, including the pseudo‐response regulator 7 (PRR7), an essential component of the circadian clock in plants. In addition, we successfully tested an approach to detect polygenic adaptation in biological pathways, allowing us to identify the flavonoid biosynthesis pathway as a candidate stress‐response pathway that deserves further attention in other plants. Finally, our study contributes to the emerging view that explaining contemporary standing genetic variation requires considering adaptation to past climates, especially for long‐lived trees.

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.

Abstract Plant species aboveground allometry can be viewed as a functional trait that reflects the evolutionary trade‐off between above‐ and belowground resources. In forest trees, allometry is related to productivity and resilience in different environments, and it is tightly connected with a compromise between efficiency‐safety and competitive ability. A better understanding on how this trait varies within and across species is critical to determine the potential of a species/population to perform along environmental gradients. We followed a hierarchical framework to assess tree height‐diameter allometry variation within and across four common European Pinus species. Tree height‐diameter allometry variation was a function of solely genetic components –approximated by either population effects or clinal geographic responses of the population's site of origin– and differential genetic plastic responses –approximated by the interaction between populations and two climatic variables of the growing sites (temperature and precipitation)–. Our results suggest that, at the species level, climate of the growing sites set the tree height‐diameter allometry of xeric and mesic species (Pinus halepensis, P. pinaster and P. nigra) apart from the boreal species (P. sylvestris), suggesting a weak signal of their phylogenies in the tree height‐diameter allometry variation. Moreover, accounting for interpopulation variability within species for the four pine species aided to: (1) detect genetic differences among populations in allometry variation, which in P. nigra and P. pinaster were linked to gene pools –genetic diversity measurements–; (2) reveal the presence of differential genetic variation in plastic responses along two climatic gradients in tree allometry variation. In P. sylvestris and P. nigra, genetic variation was the result of adaptive patterns to climate, while in P. pinaster and P. halepensis, this signal was either weaker or absent, respectively; and (3) detect local adaptation in the exponent of the tree height‐diameter allometry relationship in two of the four species (P. sylvestris and P. nigra), as it was a function of populations' latitude and altitude variables. Our findings suggest that the four species have been subjected to different historical and climatic constraints that might have driven their aboveground allometry and promoted different life strategies.

Aim of study: Our main goal is to determine the relationship between cone production and radial growth in Pinus pinaster Ait. under different climatic conditions across the Iberian Peninsula. Area of study: Coca Intensive Sampling Plateau, Northern Plateau (Spain). Material and methods: Cone counts were conducted on an intensive monitoring plot in Coca (North-Central Spain) during the years 2000, 2006 and 2007. A ZIP (zero-inflated Poisson) model was adjusted for simultaneously estimating the probability of obtaining crop cones and its amount. The Northern Atlantic Oscillation (NAO) index was used as explanatory variable, together with a wide variety of tree and local stand variables. Climate (as evaluated by NAO), local stand density (here estimated from the six nearest trees), tree size and vigor, competition and growth efficiency significantly influenced both occurrence and intensity of cone production. Main results: ZIP models for predicting reproductive effort seems an adequate tool to predict reproductive responses to climatic fluctuations and the resulting future species distribution in the face of climate change, as well as to identify silviculture actions that would promote reproductive success in naturally-regenerated stands, list and discuss relevant results (including numeric values of experimental results). Research highlights: Climate, stand density and tree conditions (size and vigor, competition and growth efficiency) influence significantly both cone occurrence and intensity of fruiting as shown by a ZIP model. As the climate variables included in the model (based on Northern Atlantic Oscillation, NAO) are general and easily obtained, the proposed model has practical applicability to predicting Pinus pinaster cone production in the Iberian Peninsula.

Aim of study: To characterize and compare the genetic resources and gain some insights into the evolutionary history of Aleppo pine (Pinus halepensis Mill.) and Brutia pine (Pinus brutia Ten.) species which are both distributed across more than 8 million hectares of area in the Mediterranean Basin. Study area: Fifty-six populations from eight Mediterranean basin countries where P. halepensis and P. brutia species are located. Materials and methods: We analyzed 1344 seeds belonging to 56 populations using five cpSSR primers (Pt15169, Pt30204, Pt41093, Pt87268, and Pt110048). Main results: The analysis of molecular variance (AMOVA) revealed that the genetic diversity among the Brutia pine populations was slightly higher than that of Aleppo pine (27.06% and 24.27%, respectively). The Aleppo pine populations separately displayed a clear east-west differentiation across the Mediterranean Basin, confirming previous results using other markers. Although the Brutia pine populations showed no spatial genetic pattern, geographically close populations and/or populations from their continual distribution range were genetically closer than the fragmented and/or ecologically marginal populations. Research highlights: The seven Aleppo pine populations from the eastern range (Türkiye, Greece, and Italy) were more than two-fold diverse than the 13 populations from the western range (Spain and Morocco). The eastern range of Aleppo pine and Brutia pine populations had similar levels of genetic diversity parameters. These results suggested that the Eastern Mediterranean Basin is a possible genetic diversity center for the two pine species.

The analysis of fine-scale spatial genetic structure (FSGS) within populations can provide insights into eco-evolutionary processes. Restricted dispersal and locally occurring genetic drift are the primary causes for FSGS at equilibrium, as described in the isolation by distance (IBD) model. Beyond IBD expectations, spatial, environmental or historical factors can affect FSGS. We examined FSGS in seven African and Neotropical populations of the late-successional rain forest tree Symphonia globulifera L. f. (Clusiaceae) to discriminate the influence of drift-dispersal vs. landscape/ecological features and historical processes on FSGS. We used spatial principal component analysis and Bayesian clustering to assess spatial genetic heterogeneity at SSRs and examined its association with plastid DNA and habitat features. African populations (from Cameroon and São Tomé) displayed a stronger FSGS than Neotropical populations at both marker types (mean Sp = 0.025 vs. Sp = 0.008 at SSRs) and had a stronger spatial genetic heterogeneity. All three African populations occurred in pronounced altitudinal gradients, possibly restricting animal-mediated seed dispersal. Cyto-nuclear disequilibria in Cameroonian populations also suggested a legacy of biogeographic history to explain these genetic patterns. Conversely, Neotropical populations exhibited a weaker FSGS, which may reflect more efficient wide-ranging seed dispersal by Neotropical bats and other dispersers. The population from French Guiana displayed an association of plastid haplotypes with two morphotypes characterized by differential habitat preferences. Our results highlight the importance of the microenvironment for eco-evolutionary processes within persistent tropical tree populations.

Understanding tree recruitment is needed to forecast future forest distribution. Many studies have reported the relevant ecological factors that affect recruitment success in trees, but the potential for genetic-based differences in recruitment has often been neglected. In this study, we established a semi-natural reciprocal sowing experiment to test for local adaptation and microenvironment effects (evaluated here by canopy cover) in the emergence and early survival of maritime pine (Pinus pinaster Aiton), an emblematic Mediterranean forest tree. A novel application of molecular markers was also developed to test for family selection and, thus, for potential genetic change over generations. Overall, we did not find evidence to support local adaptation at the recruitment stage in our semi-natural experiment. Moreover, only weak family selection (if any) was found, suggesting that in stressful environments with low survival, stochastic processes and among-year climate variability may drive recruitment. Nevertheless, our study revealed that, at early stages of recruitment, microenvironments may favor the population with the best adapted life strategy, irrespectively of its (local or non-local) origin. We also found that emergence time is a key factor for seedling survival in stressful Mediterranean environments. Our study highlights the complexity of the factors influencing the early stages of establishment of maritime pine and provides insights into possible management actions aimed at environmental change impact mitigation. In particular, we found that the high stochasticity of the recruitment process in stressful environments and the differences in population-specific adaptive strategies may difficult assisted migration schemes.

Analysing pollen movement is a key to understanding the reproductive system of plant species and how it is influenced by the spatial distribution of potential mating partners in fragmented populations. Here we infer parameters related to levels of pollen movement and diversity of the effective pollen cloud for the wind-pollinated shrub Pistacia lentiscus across a highly disturbed landscape using microsatellite loci. Paternity analysis and the indirect KinDist and Mixed Effect Mating models were used to assess mating patterns, the pollen dispersal kernel, the effective number of males (N(ep)) and their relative individual fertility, as well as the existence of fine-scale spatial genetic structure in adult plants. All methods showed extensive pollen movement, with high rates of pollen flow from outside the study site (up to 73-93%), fat-tailed dispersal kernels and large average pollination distances (δ = 229-412 m). However, they also agreed in detecting very few pollen donors (N(ep) = 4.3-10.2) and a large variance in their reproductive success: 70% of males did not sire any offspring among the studied female plants and 5.5% of males were responsible for 50% of pollinations. Although we did not find reduced levels of genetic diversity, the adult population showed high levels of biparental inbreeding (14%) and strong spatial genetic structure (S(p) = 0.012), probably due to restricted seed dispersal and scarce safe sites for recruitment. Overall, limited seed dispersal and the scarcity of successful pollen donors can be contributing to generate local pedigrees and to increase inbreeding, the prelude of genetic impoverishment.

BACKGROUND: Host-pathogen interactions that may lead to a competitive co-evolution of virulence and resistance mechanisms present an attractive system to study molecular evolution because strong, recent (or even current) selective pressure is expected at many genomic loci. However, it is unclear whether these selective forces would act to preserve existing diversity, promote novel diversity, or reduce linked neutral diversity during rapid fixation of advantageous alleles. In plants, the lack of adaptive immunity places a larger burden on genetic diversity to ensure survival of plant populations. This burden is even greater if the generation time of the plant is much longer than the generation time of the pathogen. METHODOLOGY/PRINCIPAL FINDINGS: Here, we present nucleotide polymorphism and substitution data for 41 candidate genes from the long-lived forest tree loblolly pine, selected primarily for their prospective influences on host-pathogen interactions. This dataset is analyzed together with 15 drought-tolerance and 13 wood-quality genes from previous studies. A wide range of neutrality tests were performed and tested against expectations from realistic demographic models. CONCLUSIONS/SIGNIFICANCE: Collectively, our analyses found that axr (auxin response factor), caf1 (chromatin assembly factor) and gatabp1 (gata binding protein 1) candidate genes carry patterns consistent with directional selection and erd3 (early response to drought 3) displays patterns suggestive of a selective sweep, both of which are consistent with the arm-race model of disease response evolution. Furthermore, we have identified patterns consistent with diversifying selection at erf1-like (ethylene responsive factor 1), ccoaoemt (caffeoyl-CoA-O-methyltransferase), cyp450-like (cytochrome p450-like) and pr4.3 (pathogen response 4.3), expected under the trench-warfare evolution model. Finally, a drought-tolerance candidate related to the plant cell wall, lp5, displayed patterns consistent with balancing selection. In conclusion, both arms-race and trench-warfare models seem compatible with patterns of polymorphism found in different disease-response candidate genes, indicating a mixed strategy of disease tolerance evolution for loblolly pine, a major tree crop in southeastern United States.