Climatic niche and neutral genetic diversity of the six Iberian pine species: a retrospective and prospective view

Stone Pine (Pinus pinea L.) High-Added-Value Genetics: An Overview

Stone pine (Pinus pinea L.) has received limited attention in terms of genetic research. However, genomic techniques hold promise for decoding the stone pine genome and contributing to developing a more resilient bioeconomy. Retrotransposon and specific genetic markers are effective tools for determining population-specific genomic diversity. Studies on the transcriptome and proteome have identified differentially expressed genes PAS1, CLV1, ATAF1, and ACBF involved in shoot bud formation. The stone pine proteome shows variation among populations and shows the industrial potential of the enzyme pinosylvin. Microsatellite studies have revealed low levels of polymorphism and a unique genetic diversity in stone pine, which may contribute to its environmental adaptation. Transcriptomic and proteomic analyses uncover the genetic and molecular responses of stone pine to fungal infections and nematode infestations, elucidating the defense activation, gene regulation, and the potential role of terpenes in pathogen resistance. Transcriptomics associated with carbohydrate metabolism, dehydrins, and transcription factors show promise as targets for improving stone pine's drought stress response and water retention capabilities. Stone pine presents itself as an important model tree for studying climate change adaptation due to its characteristics. While knowledge gaps exist, stone pine's genetic resources hold significant potential, and ongoing advancements in techniques offer prospects for future exploration.

Measuring recent effective gene flow among large populations in Pinus sylvestris: Local pollen shedding does not preclude substantial long-distance pollen immigration

The estimation of recent gene flow rates among vast and often weakly genetically differentiated tree populations remains a great challenge. Yet, empirical information would help understanding the interaction between gene flow and local adaptation in present-day non-equilibrium forests. We investigate here recent gene flow rates between two large native Scots pine (Pinus sylvestris L.) populations in central Iberian Peninsula (Spain), which grow on contrasting edaphic conditions six kilometers apart from each other and show substantial quantitative trait divergence in common garden experiments. Using a sample of 1,200 adult and offspring chloroplast-microsatellite haplotypes and a Bayesian inference model, we estimated substantial male gametic gene flow rates (8 and 21%) between the two natural populations, and even greater estimated immigration rates (42 and 64%) from nearby plantations into the two natural populations. Our results suggest that local pollen shedding within large tree populations does not preclude long-distance pollen immigration from large external sources, supporting the role of gene flow as a homogenizing evolutionary force contributing to low molecular genetic differentiation among populations of widely distributed wind-pollinated species. Our results also indicate the high potential for reproductive connectivity in large fragmented populations of wind-pollinated trees, and draw attention to a potential scenario of adaptive genetic divergence in quantitative traits under high gene flow.

Candidate Genes for the High-Altitude Adaptations of Two Mountain Pine Taxa

Mountain plants, challenged by vegetation time contractions and dynamic changes in environmental conditions, developed adaptations that help them to balance their growth, reproduction, survival, and regeneration. However, knowledge regarding the genetic basis of species adaptation to higher altitudes remain scarce for most plant species. Here, we attempted to identify such corresponding genomic regions of high evolutionary importance in two closely related European pines, Pinus mugo and P. uncinata, contrasting them with a reference lowland relative—P. sylvestris. We genotyped 438 samples at thousands of single nucleotide polymorphism (SNP) markers, tested their genetic differentiation and population structure followed by outlier detection and gene ontology annotations. Markers clearly differentiated the species and uncovered patterns of population structure in two of them. In P. uncinata three Pyrenean sites were grouped together, while two outlying populations constituted a separate cluster. In P. sylvestris, Spanish population appeared distinct from the remaining four European sites. Between mountain pines and the reference species, 35 candidate genes for altitude-dependent selection were identified, including such encoding proteins responsible for photosynthesis, photorespiration and cell redox homeostasis, regulation of transcription, and mRNA processing. In comparison between two mountain pines, 75 outlier SNPs were found in proteins involved mainly in the gene expression and metabolism.

Champions of winter survival: cold acclimation and molecular regulation of cold hardiness in evergreen conifers

Evergreen conifers are champions of winter survival, based on their remarkable ability to acclimate to cold and develop cold hardiness. Counterintuitively, autumn cold acclimation is triggered not only by exposure to low temperature, but also by a combination of decreasing temperature, decreasing photoperiod and changes in light quality. These environmental cues control a network of signaling pathways that coordinate cold acclimation and cold hardiness in overwintering conifers, leading to cessation of growth, bud dormancy, freezing tolerance and changes in energy metabolism. Advances in genomic, transcriptomic and metabolomic tools for conifers have improved our understanding of how trees sense and respond to changes in temperature and light during cold acclimation and the development of cold hardiness, but there remain considerable gaps deserving further research in conifers. In the first section of this review, we focus on the physiological mechanisms used by evergreen conifers to adjust metabolism seasonally and to protect overwintering tissues against winter stresses. In the second section, we review how perception of low temperature and photoperiod regulate the induction of cold acclimation. Finally, we explore the evolutionary context of cold acclimation in conifers and evaluate challenges imposed on them by changing climate and discuss emerging areas of research in the field.

The legacy of climate variability over the last century on populations' phenotypic variation in tree height

Phenotypic plasticity and local adaptation are the two main processes underlying trait variability. Under rapid environmental change, phenotypic plasticity, if adaptive, could increase the odds for organisms to persist. However, little is known on how environmental variation has shaped plasticity across species ranges over time. Here, we assess whether the portion of phenotypic variation of tree populations linked to the environment is related to the inter-annual climate variability of the last century and how it varies among populations across species ranges and age. To this aim, we used 372,647 individual tree height measurements of three pine species found in low elevation forests in Europe: Pinus nigra Arnold, P. pinaster Aiton and P. pinea L. Measurements were taken in a network of 38 common gardens established in Europe and North Africa with 315 populations covering the distribution range of the species. We fitted linear mixed-effect models of tree height as a function of age, population, climate and competition effects. Models allowed us to estimate tree height response curves at the population level and indexes of populations' phenotypic variation, as a proxy of phenotypic plasticity, at 4, 8 and 16 years old, and relate these indexes to the inter-annual climate variability of the last century. We found that phenotypic variation in tree height was higher in young trees than in older ones. We also found that P. pinea showed the highest phenotypic variation in tree height compared with P. pinaster and P. nigra. Finally, phenotypic variation in tree height may be partly adaptive, and differently across species, as climate variability during the last century at the origin of the populations explained between 51 and 69% of the current phenotypic variation of P. nigra and P. pinea, almost twice of the levels of P. pinaster. MAIN CONCLUSIONS: Populations' phenotypic variation in tree height is largely explained by the climate variability that the populations experienced during the last century, which we attribute to the genetic diversity among populations.

Seed Sourcing Strategies Considering Climate Change Forecasts: A Practical Test in Scots Pine

Research Highlights: We experimentally tested different seed sourcing strategies (local, predictive, climate-predictive, climate-adjusted, composite and admixture) under a climate change high emissions scenario using a Scots pine multi-site provenance test. Background and Objectives: There is an urgent need to conserve genetic resources and to support resilience of conifer species facing expected changes and threats. Seed sourcing strategies have been proposed to maximize the future adaptation and resilience of our forests. However, these proposals are yet to be tested, especially in long-lived organisms as forest trees, due to methodological constraints. In addition, some methods rely on the transfer of material from populations matching the future conditions of the sites. However, at the rear edge of the species, some specific problems (high fragmentation, high genetic differentiation, role of genetic drift) challenge the theoretical expectations of some of these methods. Materials and Methods: We used a Scots pine multi-site provenance test, consisting of seventeen provenances covering the distribution range of the species in Spain tested in five representative sites. We measured height, diameter and survival at 5, 10 and 15 years after planting. We simulated populations of 50 trees by bootstrapping material of the provenance test after removing the intra-site environmental effects, simulating different seed sourcing strategies. Results: We found that local and predictive methods behaved better than methods based on the selection of future climate-matching strategies (predictive-climate and climate-adjusted) and those combining several seed sources (composite and admixture seed sourcing strategies). Conclusions: Despite the theoretical expectations, for Scots pine, a forest tree species at its rear edge of its distribution, seed-sourcing methods based on climate matching or a combination of seed sources do not perform better than traditional local or predictive methods or they are not feasible because of the lack of future climate-matching populations.

Exome Resequencing Reveals Evolutionary History, Genomic Diversity, and Targets of Selection in the Conifers Pinus taeda and Pinus elliottii

Loblolly pine (Pinus taeda) and slash pine (Pinus elliottii) are ecologically and economically important pine species that dominate many forest ecosystems in the southern United States, but like all conifers, the study of their genetic diversity and demographic history has been hampered by their large genome size. A small number of studies mainly based on candidate-gene sequencing have been reported for P. taeda to date, whereas none are available for P. elliottii. Targeted exome resequencing has recently enabled population genomics studies for conifers, approach used here to assess genomic diversity, signatures of selection, population structure, and demographic history of P. elliottii and P. taeda. Extensive similarities were revealed between these species: both species feature rapid linkage disequilibrium decay and high levels of genetic diversity. Moreover, genome-wide positive correlations for measures of genetic diversity between the species were also observed, likely due to shared structural genomic constraints. Also, positive selection appears to be targeting a common set of genes in both pines. Demographic history differs between both species, with only P. taeda being affected by a dramatic bottleneck during the last glacial period. The ability of P. taeda to recover from a dramatic reduction in population size while still retaining high levels of genetic diversity shows promise for other pines facing environmental stressors associated with climate change, indicating that these too may be able to adapt successfully to new future conditions even after a drastic population size contraction.

Adaptation and plasticity in aboveground allometry variation of four pine species along environmental gradients

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.

Contrasting patterns of genetic variation in core and peripheral populations of highly outcrossing and wind pollinated forest tree species

Gene flow tends to have a homogenising effect on a species' background genetic variation over large geographical areas. However, it is usually unknown to what extent the genetic structure of populations is influenced by gene exchange between core and peripheral populations that may represent stands of different evolutionary and demographic history. In this study, we looked at the patterns of population differentiation in Scots pine-a highly outcrossing and wind pollinated conifer species that forms large ecosystems of great ecological and economic importance in Europe and Asia. A set of 13 polymorphic nuclear microsatellite loci was analysed to infer the genetic relationships among 24 populations (676 individuals) from Europe and Asia Minor. The study included specimens from the primary continuous range and from isolated, marginal stands that are considered to be autochthonous populations representative of the species' putative refugial areas. Despite their presumably different histories, a similar level of genetic variation and no evidence of a population bottleneck was found across the populations. Differentiation among populations was relatively low (average FST = 0.035); however, the population structure was not homogenous, which was clearly evident from the allelic frequency spectra and Bayesian assignment analysis. Significant differentiation over short geographical distances was observed between isolated populations within the Iberian and Anatolian Peninsulas (Asia Minor), which contrasted with the absence of genetic differentiation observed between distant populations e.g., between central and northern Europe. The analysed populations were assigned to several groups that corresponded to the geographical regions of their occurrence. These results will be useful in genetics studies in Scots pine that aim to link nucleotide and phenotypic variation across the species distribution range and for development of sustainable breeding and management programs.

Phylogeography of the Spanish Moon Moth Graellsia isabellae (Lepidoptera, Saturniidae)

BACKGROUND

Geographic and demographic factors as well as specialisation to a new host-plant may lead to host-associated differentiation in plant-feeding insects. We explored the phylogeography of a protected moth, Graellsia isabellae, and its two recognised host-plant species (Pinus sylvestris and P. nigra) in order to seek for any concordance useful to disentangle the evolutionary history of this iconic lepidopteran.

RESULTS

DNA variation in one mitochondrial marker and nine nuclear microsatellite loci revealed a strong phylogeographic pattern across 28 populations of G. isabellae studied in Spain and France comprising six groups mostly distributed along different mountain ranges. Reanalysis of a previously published chloroplast microsatellite dataset revealed a three and two-group structure for Spanish P. sylvestris and P. nigra, respectively. Overall, the population groupings of this protected moth did not match the ones of P. sylvestris and P. nigra.

CONCLUSIONS

There was no evidence of host-associated differentiation between populations using P. sylvestris and the ones inhabiting P. nigra. The two major mitochondrial clades of G. isabellae likely diverged before the Last Glacial Maximum and geographically separated the species into a "southern" (Central and Southern Iberian clusters) and a "northern" lineage (Eastern Iberian, Pyrenean and French Alpine clusters). The Eastern Iberian System, where this insect uses both host-plants, harboured the highest level of genetic diversity. Such a group independently colonised the West and East parts of the Pyrenees. Our results point to a native origin for the French populations occurring in the Alps, genetically related to the Eastern Iberian and Pyrenean sites. The Central Iberian group derived from Southern Iberian ancestors. Secondary contacts were inferred between the Southern/Central Iberian populations and Eastern Iberian cluster as well as between the two Pyrenean ones. The mito-nuclear discordance observed with regard to the Eastern Iberian cluster is congruent with a secondary contact after the evolution of mito-nuclear incompatibilities in geographically isolated areas.

Adaptive differentiation in seedling traits in a hybrid pine species complex, Pinus densata and its parental species, on the Tibetan Plateau

Evidence from molecular genetics demonstrates that Pinus densata is a natural homoploid hybrid originating from the parent species Pinus tabuliformis and Pinus yunnanensis, and ecological selection may have played a role in the speciation of P. densata. However, data on differentiation in adaptive traits in the species complex are scarce. In this study, we performed a common garden test on 16 seedling traits to examine the differences between P. densata and its parental species in a high altitude environment. We found that among the 16 analyzed traits, 15 were significantly different among the species. Pinus tabuliformis had much earlier bud set and a relatively higher bud set ratio but poorer seedling growth, and P. yunnanensis had opposite responses for the same traits. P. densata had the greatest fitness with higher viability and growth rates than the parents. The relatively high genetic contribution of seedling traits among populations suggested that within each species the evolutionary background is complex. The correlations between the seedling traits of a population within a species and the environmental factors indicated different impacts of the environment on species evolution. The winter temperature is among the most important climate factors that affected the fitness of the three pine species. Our investigation provides empirical evidence on adaptive differentiation among this pine species complex at seedling stages.

Plasticity in dendroclimatic response across the distribution range of Aleppo pine (Pinus halepensis)

We investigated the variability of the climate-growth relationship of Aleppo pine across its distribution range in the Mediterranean Basin. We constructed a network of tree-ring index chronologies from 63 sites across the region. Correlation function analysis identified the relationships of tree-ring index to climate factors for each site. We also estimated the dominant climatic gradients of the region using principal component analysis of monthly, seasonal, and annual mean temperature and total precipitation from 1,068 climatic gridpoints. Variation in ring width index was primarily related to precipitation and secondarily to temperature. However, we found that the dendroclimatic relationship depended on the position of the site along the climatic gradient. In the southern part of the distribution range, where temperature was generally higher and precipitation lower than the regional average, reduced growth was also associated with warm and dry conditions. In the northern part, where the average temperature was lower and the precipitation more abundant than the regional average, reduced growth was associated with cool conditions. Thus, our study highlights the substantial plasticity of Aleppo pine in response to different climatic conditions. These results do not resolve the source of response variability as being due to either genetic variation in provenance, to phenotypic plasticity, or a combination of factors. However, as current growth responses to inter-annual climate variability vary spatially across existing climate gradients, future climate-growth relationships will also likely be determined by differential adaptation and/or acclimation responses to spatial climatic variation. The contribution of local adaptation and/or phenotypic plasticity across populations to the persistence of species under global warming could be decisive for prediction of climate change impacts across populations. In this sense, a more complex forest dynamics modeling approach that includes the contribution of genetic variation and phenotypic plasticity can improve the reliability of the ecological inferences derived from the climate-growth relationships.

Genetic variation of drought tolerance in Pinus pinaster at three hierarchical levels: a comparison of induced osmotic stress and field testing

Understanding the survival capacity of forest trees to periods of severe water stress could improve knowledge of the adaptive potential of different species under future climatic scenarios. In long lived organisms, like forest trees, the combination of induced osmotic stress treatments and field testing can elucidate the role of drought tolerance during the early stages of establishment, the most critical in the life of the species. We performed a Polyethylene glycol-osmotic induced stress experiment and evaluated two common garden experiments (xeric and mesic sites) to test for survival and growth of a wide range clonal collection of Maritime pine. This study demonstrates the importance of additive vs non additive effects for drought tolerance traits in Pinus pinaster, and shows differences in parameters determining the adaptive trajectories of populations and family and clones within populations. The results show that osmotic adjustment plays an important role in population variation, while biomass allocation and hydric content greatly influence survival at population level. Survival in the induced osmotic stress experiment presented significant correlations with survival in the xeric site, and height growth at the mesic site, at population level, indicating constraints of adaptation for those traits, while at the within population level no significant correlation existed. These results demonstrate that population differentiation and within population genetic variation for drought tolerance follow different patterns.

The intraspecific genetic variability of siliceous and calcareous Gentiana species is shaped by contrasting demographic and re-colonization processes

The Ciminalis section of Gentiana comprises seven species, two of them growing on siliceous substrates (G. alpina and G. acaulis), the other ones being calcareous taxa (G. clusii, G. angustifolia, G. ligustica, G. occidentalis and G. dinarica). A total of 515 individuals from 183 populations over the entire Ciminalis distribution range was analyzed using four chloroplast loci (trnH-psbA, matK, rpoB and rpoC1) and the nuclear ribosomal marker ITS2. The siliceous species display only two chloroplast haplotypes each and are both characterized by patterns of range expansions all over the Alps. Conversely, the calcareous species are on average more diverse (two to 13 haplotypes per species) with strong patterns of local structuring. We suggest that the occurrence of many calcareous refugia at the periphery of the Alps must have led to local adaptation and morphological diversification, and helped preserving intraspecific diversities during the last glaciations for the associated taxa. ITS2 was more efficient in delineating species boundaries than the chloroplast markers for which several haplotypes are shared among species. This might be either due to chloroplast capture among species and/or to recent divergence. Species adapted to the same substrate are generally only distantly related when they co-occur in the same place. For both types of markers, G. clusii is found genetically distant from all other species.

Potential for evolutionary responses to climate change - evidence from tree populations

Evolutionary responses are required for tree populations to be able to track climate change. Results of 250 years of common garden experiments show that most forest trees have evolved local adaptation, as evidenced by the adaptive differentiation of populations in quantitative traits, reflecting environmental conditions of population origins. On the basis of the patterns of quantitative variation for 19 adaptation-related traits studied in 59 tree species (mostly temperate and boreal species from the Northern hemisphere), we found that genetic differentiation between populations and clinal variation along environmental gradients were very common (respectively, 90% and 78% of cases). Thus, responding to climate change will likely require that the quantitative traits of populations again match their environments. We examine what kind of information is needed for evaluating the potential to respond, and what information is already available. We review the genetic models related to selection responses, and what is known currently about the genetic basis of the traits. We address special problems to be found at the range margins, and highlight the need for more modeling to understand specific issues at southern and northern margins. We need new common garden experiments for less known species. For extensively studied species, new experiments are needed outside the current ranges. Improving genomic information will allow better prediction of responses. Competitive and other interactions within species and interactions between species deserve more consideration. Despite the long generation times, the strong background in quantitative genetics and growing genomic resources make forest trees useful species for climate change research. The greatest adaptive response is expected when populations are large, have high genetic variability, selection is strong, and there is ecological opportunity for establishment of better adapted genotypes.

Large-scale longitudinal gradients of genetic diversity: a meta-analysis across six phyla in the Mediterranean basin

Biodiversity is the diversity of life at all scales, from genes to ecosystems. Predicting its patterns of variation across the globe is a fundamental issue in ecology and evolution. Diversity within species, that is, genetic diversity, is of prime importance for understanding past and present evolutionary patterns, and highlighting areas where conservation might be a priority. Using published data on the genetic diversity of species whose populations occur in the Mediterranean basin, we calculated a coefficient of correlation between within-population genetic diversity indices and longitude. Using a meta-analysis framework, we estimated the role of biological, ecological, biogeographic, and marker type factors on the strength and magnitude of this correlation in six phylla. Overall, genetic diversity increases from west to east in the Mediterranean basin. This correlation is significant for both animals and plants, but is not uniformly expressed for all groups. It is stronger in the southern than in the northern Mediterranean, in true Mediterranean plants than in plants found at higher elevations, in trees than in other plants, and in bi-parentally and paternally than in maternally inherited DNA makers. Overall, this correlation between genetic diversity and longitude, and its patterns across biological and ecological traits, suggests the role of two non-mutually exclusive major processes that shaped the genetic diversity in the Mediterranean during and after the cold periods of the Pleistocene: east-west recolonization during the Holocene and population size contraction under local Last Glacial Maximum climate in resident western and low elevation Mediterranean populations.

Nuclear microsatellites for Pinus pinea (Pinaceae), a genetically depauperate tree, and their transferability to P. halepensis

PREMISE OF THE STUDY

Pinus pinea is one of the few widespread plant species that are also genetically depauperate. It is also an important commercial species with high market value seeds. A deeper knowledge of the existing population genetic variation was needed.

METHODS AND RESULTS

Twelve nuclear microsatellites were isolated from genomic and cDNA sequences and screened for variability in 729 individuals from 33 natural populations. Low level of genetic variability was confirmed with average expected heterozygosity of 0.11. Hardy-Weinberg equilibrium expectations were not met in only ∼10% of the possible locus/population combinations. All loci were in linkage equilibrium, and the frequency of null alleles was very low (≤1% in 332 out of 396 locus/population combinations). Nine out of the 12 microsatellites were successfully transferred to P. halepensis.

CONCLUSIONS

Despite low polymorphism, these new markers will be useful to resolve population structure and hold potential for seed origin identification and traceability.

Wind pollination over mesoscale distances: an investigation with Scots pine

There is a gap between the order of magnitude of maximum documented distances of airborne tree pollen transport (up to 10(2)-10(3) km) and effective wind pollination (up to 10(1) km), which may partly derive from greater difficulties in detecting the latter. This study aims to assess wind pollination over scales closer to the maximum observed physical pollen transport distances. The origin of effective pollen immigrants into a strongly isolated Iberian Pinus sylvestris remnant was investigated using paternally inherited microsatellite markers and maximum-likelihood estimation combined with Monte Carlo assessment of parameter uncertainty. The results revealed significant effective pollen flow (up to 4.4%) from a large population located 100 km away, suggesting that the well-known mesoscale airborne transport of viable pine pollen can result in successful pollination over larger scales than previously reported for wind-pollinated tree species. This study supports the view that the gap between documented potential and effective wind pollen dispersal scales might not accurately reflect biological reality. Expanding the expected range of effective wind pollination has an impact on the assessment of a wide range of ecological and evolutionary processes, including reproductive assurance on fragmentation or colonization, metapopulation connectivity and interactions with local adaptation in heterogeneous habitats.

The Strait of Gibraltar as a major biogeographic barrier in Mediterranean conifers: a comparative phylogeographic survey

The Strait of Gibraltar (SG) is reputed for being both a bridge and a geographic barrier to biological exchanges between Europe and Africa. Major genetic breaks associated with this strait have been identified in various taxa, but it is unknown whether these disjunctions have been produced simultaneously or by independent biogeographic processes. Here, the genetic structure of five conifers distributed on both sides of the SG was investigated using mitochondrial (nad1 b/c, nad5-1, nad5-4 and nad7-1) and chloroplast (Pt1254, Pt15169, Pt30204, Pt36480, Pt71936 and Pt87268) DNA markers. The distribution of genetic variation was partially congruent between types of markers within the same species. Across taxa, there was a significant overlapping between the SG and the genetic breaks detected, especially for the four Tertiary species surveyed (Abies pinsapo complex, Pinus nigra, Pinus pinaster and Taxus baccata). For most of these taxa, the divergence of populations across the SG could date back to long before the Pleistocene glaciations. However, their strongly different cpDNA G(ST) and R(ST) values point out that they have had dissimilar population histories, which might include contrasting amounts of pollen-driven gene flow since their initial establishment in the region. The fifth species, Pinus halepensis, was genetically depauperated and homogenous on both sides of the SG. A further analysis of nuclear DNA sequences with coalescent-based isolation with migration models suggests a Pleistocene divergence of P. halepensis populations across the SG, which is in sharp contrast with the pre-Pleistocene divergence dates obtained for P. pinaster. Altogether, these results indicate that the genetic breaks observed across this putative biogeographical barrier have been produced by independent evolutionary processes related to the biological history of each individual species instead of a common vicariant phenomenon.

Back to nature: ecological genomics of loblolly pine (Pinus taeda, Pinaceae)

Genetic variation is often arrayed in latitudinal or altitudinal clines, reflecting either adaptation along environmental gradients, migratory routes, or both. For forest trees, climate is one of the most important drivers of adaptive phenotypic traits. Correlations of single and multilocus genotypes with environmental gradients have been identified for a variety of forest trees. These correlations are interpreted normally as evidence of natural selection. Here, we use a genome-wide dataset of single nucleotide polymorphisms (SNPs) typed from 1730 loci in 682 loblolly pine (Pinus taeda L.) trees sampled from 54 local populations covering the full-range of the species to examine allelic correlations to five multivariate measures of climate. Applications of a Bayesian generalized linear mixed model, where the climate variable was a fixed effect and an estimated variance-covariance matrix controlled random effects due to shared population history, identified several well-supported SNPs associating to principal components corresponding to geography, temperature, growing degree-days, precipitation and aridity. Functional annotation of those genes with putative orthologs in Arabidopsis revealed a diverse set of abiotic stress response genes ranging from transmembrane proteins to proteins involved in sugar metabolism. Many of these SNPs also had large allele frequency differences among populations (F(ST) = 0.10-0.35). These results illustrate a first step towards a ecosystem perspective of population genomics for non-model organisms, but also highlight the need for further integration of the methodologies employed in spatial statistics, population genetics and climate modeling during scans for signatures of natural selection from genomic data.