Evolutionary history of the mediterranean Pinus halepensis-brutia species complex using gene-resequencing and transcriptomic approaches.

KEY MESSAGE: Complementary gene-resequencing and transcriptomic approaches reveal contrasted evolutionary histories in a species complex. Pinus halepensis and Pinus brutia are closely related species that can intercross, but occupy different geographical ranges and bioclimates. To study the evolution of this species complex and to provide genomic resources for further research, we produce and analyze two new complementary sets of genetic resources: (i) a set of 172 re-sequenced genomic target loci analyzed in 45 individuals, and (ii) a set of 11 transcriptome assemblies. These two datasets provide insights congruent with previous studies: P. brutia displays high level of genetic diversity and no genetic sub-structure, while P. halepensis shows three main genetic clusters, the western Mediterranean and North African clusters displaying much lower genetic diversity than the eastern Mediterranean cluster, the latter cluster having similar genetic diversity to P. brutia. In addition, these datasets provide new insights on the timing of the species-complex history: the two species would have split at the end of the tertiary, and the changing climatic conditions of the Mediterranean region at the end of the Tertiary-beginning of the Quaternary, together with the distinct species tolerance to harsh climatic conditions would have resulted in different geographic distributions, demographic histories and genetic patterns of the two pines. The multiple glacial-interglacial cycles during the Quaternary would have led to the expansion of P. brutia in the Middle East, while P. halepensis would have been through bottlenecks. The last glaciations, from 0.6 Mya on, would have affected further the Western genetic pool of P. halepensis.

Ectomycorrhizal fungal diversity decreases in Mediterranean pine forests adapted to recurrent fires.

Fire is a major disturbance linked to the evolutionary history and climate of Mediterranean ecosystems, where the vegetation has evolved fire-adaptive traits (e.g., serotiny in pines). In Mediterranean forests, mutualistic feedbacks between trees and ectomycorrhizal (ECM) fungi, essential for ecosystem dynamics, might be shaped by recurrent fires. We tested how the structure and function of ECM fungal communities of Pinus pinaster and Pinus halepensis vary among populations subjected to high and low fire recurrence in Mediterranean ecosystems, and analysed the relative contribution of environmental (climate, soil properties) and tree-mediated (serotiny) factors. For both pines, local and regional ECM fungal diversity were lower in areas of high than low fire recurrence, although certain fungal species were favoured in the former. A general decline of ECM root-tip enzymatic activity for P. pinaster was associated with high fire recurrence, but not for P. halepensis. Fire recurrence and fire-related factors such as climate, soil properties or tree phenotype explained these results. In addition to the main influence of climate, the tree fire-adaptive trait serotiny recovered a great portion of the variation in structure and function of ECM fungal communities associated with fire recurrence. Edaphic conditions (especially pH, tightly linked to bedrock type) were an important driver shaping ECM fungal communities, but mainly at the local scale and probably independently of the fire recurrence. Our results show that ECM fungal community shifts are associated with fire recurrence in fire-prone dry Mediterranean forests, and reveal complex feedbacks among trees, mutualistic fungi and the surrounding environment in these ecosystems.

Increased fire frequency promotes stronger spatial genetic structure and natural selection at regional and local scales in Pinus halepensis Mill.

Background and Aims: The recurrence of wildfires is predicted to increase due to global climate change, resulting in severe impacts on biodiversity and ecosystem functioning. Recurrent fires can drive plant adaptation and reduce genetic diversity; however, the underlying population genetic processes have not been studied in detail. In this study, the neutral and adaptive evolutionary effects of contrasting fire regimes were examined in the keystone tree species Pinus halepensis Mill. (Aleppo pine), a fire-adapted conifer. The genetic diversity, demographic history and spatial genetic structure were assessed at local (within-population) and regional scales for populations exposed to different crown fire frequencies. Methods: Eight natural P. halepensis stands were sampled in the east of the Iberian Peninsula, five of them in a region exposed to frequent crown fires (HiFi) and three of them in an adjacent region with a low frequency of crown fires (LoFi). Samples were genotyped at nine neutral simple sequence repeats (SSRs) and at 251 single nucleotide polymorphisms (SNPs) from coding regions, some of them potentially important for fire adaptation. Key Results: Fire regime had no effects on genetic diversity or demographic history. Three high-differentiation outlier SNPs were identified between HiFi and LoFi stands, suggesting fire-related selection at the regional scale. At the local scale, fine-scale spatial genetic structure (SGS) was overall weak as expected for a wind-pollinated and wind-dispersed tree species. HiFi stands displayed a stronger SGS than LoFi stands at SNPs, which probably reflected the simultaneous post-fire recruitment of co-dispersed related seeds. SNPs with exceptionally strong SGS, a proxy for microenvironmental selection, were only reliably identified under the HiFi regime. Conclusions: An increasing fire frequency as predicted due to global change can promote increased SGS with stronger family structures and alter natural selection in P. halepensis and in plants with similar life history traits.

Functional outcomes of fungal community shifts driven by tree genotype and spatial-temporal factors in Mediterranean pine forests.

Fungi provide relevant ecosystem services contributing to primary productivity and the cycling of nutrients in forests. These fungal inputs can be decisive for the resilience of Mediterranean forests under global change scenarios, making necessary an in-deep knowledge about how fungal communities operate in these ecosystems. By using high-throughput sequencing and enzymatic approaches, we studied the fungal communities associated with three genotypic variants of Pinus pinaster trees, in 45-year-old common garden plantations. We aimed to determine the impact of biotic (i.e., tree genotype) and abiotic (i.e., season, site) factors on the fungal community structure, and to explore whether structural shifts triggered functional responses affecting relevant ecosystem processes. Tree genotype and spatial-temporal factors were pivotal structuring fungal communities, mainly by influencing their assemblage and selecting certain fungi. Diversity variations of total fungal community and of that of specific fungal guilds, together with edaphic properties and tree's productivity, explained relevant ecosystem services such as processes involved in carbon turnover and phosphorous mobilization. A mechanistic model integrating relations of these variables and ecosystem functional outcomes is provided. Our results highlight the importance of structural shifts in fungal communities because they may have functional consequences for key ecosystem processes in Mediterranean forests.

Adaptive evolution of Mediterranean pines.

Mediterranean pines represent an extremely heterogeneous assembly. Although they have evolved under similar environmental conditions, they diversified long ago, ca. 10 Mya, and present distinct biogeographic and demographic histories. Therefore, it is of special interest to understand whether and to what extent they have developed specific strategies of adaptive evolution through time and space. To explore evolutionary patterns, the Mediterranean pines' phylogeny was first reconstructed analyzing a new set of 21 low-copy nuclear genes with multilocus Bayesian tree reconstruction methods. Secondly, a phylogenetic approach was used to search for footprints of natural selection and to examine the evolution of multiple phenotypic traits. We identified two genes (involved in pines' defense and stress responses) that have likely played a role in the adaptation of Mediterranean pines to their environment. Moreover, few life-history traits showed historical or evolutionary adaptive convergence in Mediterranean lineages, while patterns of character evolution revealed various evolutionary trade-offs linking growth-development, reproduction and fire-related traits. Assessing the evolutionary path of important life-history traits, as well as the genomic basis of adaptive variation is central to understanding the past evolutionary success of Mediterranean pines and their future response to environmental changes.

Molecular footprints of local adaptation in two Mediterranean conifers.

This study combines neutrality tests and environmental correlations to identify nonneutral patterns of evolution in candidate genes related to drought stress in two closely related Mediterranean conifers, Pinus pinaster Ait. and P. halepensis Mill. Based on previous studies, we selected twelve amplicons covering six candidate genes that were sequenced in a large sample spanning the full range of these two species. Neutrality tests relatively robust to demography (DHEW compound test and maximum likelihood multilocus Hudson-Kreitman-Aguadé test) were used to detect selection events at different temporal scales. Environmental associations between variation at candidate genes and climatic variables were also examined. These combined approaches detected distinct genes that may be targeted by selection, most of them specific to only one of the two conifers, despite their recent divergence (<10 Ma). An exception was 4-coumarate: CoA ligase, a gene involved in the production of various important secondary products that appeared to play a role in local adaptation processes of both pines. Another remarkable result was that all significant environmental correlations involved temperature indices, highlighting the importance of this climatic factor as a selective driver on Mediterranean pines. The ability to detect natural selection at the DNA sequence level depends on the nature and the strength of the selection events, on the timescale at which they occurred, and on the sensitivity of the methods to other evolutionary forces that can mimic selection (e.g., demography and population structure). Using complementary approaches can help to capture different aspects of the evolutionary processes that govern molecular variation at both intra- and interspecific levels.