Patterns of polymorphism resulting from long-range colonization in the Mediterranean conifer Aleppo pine

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

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.

Intraspecific responses to climate reveal nonintuitive warming impacts on a widespread thermophilic conifer

Many ecologically important forest trees from dry areas have been insufficiently investigated for their ability to adapt to the challenges posed by climate change, which hampers the implementation of mitigation policies. We analyzed 14 common-garden experiments across the Mediterranean which studied the widespread thermophilic conifer Pinus halepensis and involved 157 populations categorized into five ecotypes. Ecotype-specific tree height responses to climate were applied to projected climate change (2071-2100 ad), to project potential growth patterns both locally and across the species' range. We found contrasting ecotypic sensitivities to annual precipitation but comparatively uniform responses to mean temperature, while evidence of local adaptation for tree height was limited to mesic ecotypes. We projected intriguing patterns of response range-wide, implying either height inhibition or stimulation of up to 75%, and deduced that the ecotype currently experiencing more favorable (wetter) conditions will show the largest inhibition. Extensive height reductions can be expected for coastal areas of France, Greece, Spain and northern Africa. Our findings underline the fact that intraspecific variations in sensitivity to precipitation must be considered when projecting tree height responses of dry forests to future climate. The ecotype-specific projected performances call for management activities to ensure forest resilience in the Mediterranean through, for example, tailored deployment strategies.

Looking for Local Adaptation: Convergent Microevolution in Aleppo Pine (Pinus halepensis)

Finding outlier loci underlying local adaptation is challenging and is best approached by suitable sampling design and rigorous method selection. In this study, we aimed to detect outlier loci (single nucleotide polymorphisms, SNPs) at the local scale by using Aleppo pine (Pinus halepensis), a drought resistant conifer that has colonized many habitats in the Mediterranean Basin, as the model species. We used a nested sampling approach that considered replicated altitudinal gradients for three contrasting sites. We genotyped samples at 294 SNPs located in genomic regions selected to maximize outlier detection. We then applied three different statistical methodologies-Two Bayesian outlier methods and one latent factor principal component method-To identify outlier loci. No SNP was an outlier for all three methods, while eight SNPs were detected by at least two methods and 17 were detected only by one method. From the intersection of outlier SNPs, only one presented an allelic frequency pattern associated with the elevational gradient across the three sites. In a context of multiple populations under similar selective pressures, our results underline the need for careful examination of outliers detected in genomic scans before considering them as candidates for convergent adaptation.

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.

Resource stability and geographic isolation are associated with genome divergence in western Palearctic crossbills

While many conifers produce annually variable seed crops, serotinous species (which hold seeds in cones for multiple years) represent unusually stable food resources for seed predators. Such stability is conducive to residency and potentially population divergence of consumers as exemplified by the Cassia crossbill (Loxia sinesciuris) in North America. We used genotyping by sequencing (GBS) to test whether three Mediterranean subspecies of common crossbills (L. curvirostra) associated with the serotinous Aleppo pine (Pinus halepensis) were more genetically distinct than European crossbills associated with nonserotinous conifers. We assembled a Cassia crossbill draft genome as a reference for mapping GBS reads and as a first step towards a more contiguous genome assembly. We found clear patterns of genetic divergence for each of the P. halepensis-associated subspecies. Geographic isolation, as promoted by resource stability and residency, is associated with genetic divergence of two of these subspecies. However, geographic isolation cannot account for divergence of L. c. hispana. Instead, resource stability likely contributed to divergence by reducing dispersal and increasing resource competition that may limit breeding by immigrants. In contrast, we found no differentiation among common crossbills associated with less stable resources, and only slight differentiation between common crossbills and parrot crossbills (L. pytyopsittacus). The substantial morphological divergence between common and parrot crossbills has likely originated or been maintained by selection despite gene flow generated by spatiotemporal resource fluctuation. Our results indicate that phenological as well as morphological characteristics of conifers have influenced crossbill diversification, and suggest a possible link between resource stability and population divergence.

Inferring selection in instances of long-range colonization: The Aleppo pine (Pinus halepensis) in the Mediterranean Basin

Teasing apart the effects of natural selection and demography on current allele frequencies is challenging, due to both processes leaving a similar molecular footprint. In particular, when attempting to identify selection in species that have undergone a recent range expansion, the increase in genetic drift at the edges of range expansions ("allele surfing") can be a confounding factor. To address this potential issue, we first assess the long-range colonization history of the Aleppo pine across the Mediterranean Basin, using molecular markers. We then look for single nucleotide polymorphisms (SNPs) involved in local adaptation using: (a) environmental correlation methods (bayenv2), focusing on bioclimatic variables important for the species' adaptation (i.e., temperature, precipitation and water availability); and (b) FST -related methods (pcadapt). To assess the rate of false positives caused by the allele surfing effect, these results are compared with results from simulated SNP data that mimics the species' past range expansions and the effect of genetic drift, but with no selection. We find that the Aleppo pine shows a previously unsuspected complex genetic structure across its range, as well as evidence of selection acting on SNPs involved with the response to bioclimatic variables such as drought. This study uses an original approach to disentangle the confounding effects of drift and selection in range margin populations. It also contributes to the increased evidence that plant populations are able to adapt to new environments despite the expected accumulation of deleterious mutations that takes place during long-range colonizations.

Rethinking the history of common walnut (Juglans regia L.) in Europe: Its origins and human interactions

Common walnut (Juglans regia L) is an economically important species cultivated worldwide for its high-quality wood and nuts. It is generally accepted that after the last glaciation J. regia survived and grew in almost completely isolated stands in Asia, and that ancient humans dispersed walnuts across Asia and into new habitats via trade and cultural expansion. The history of walnut in Europe is a matter of debate, however. In this study, we estimated the genetic diversity and structure of 91 Eurasian walnut populations using 14 neutral microsatellites. By integrating fossil pollen, cultural, and historical data with population genetics, and approximate Bayesian analysis, we reconstructed the demographic history of walnut and its routes of dispersal across Europe. The genetic data confirmed the presence of walnut in glacial refugia in the Balkans and western Europe. We conclude that human-mediated admixture between Anatolian and Balkan walnut germplasm started in the Early Bronze Age, and between western Europe and the Balkans in eastern Europe during the Roman Empire. A population size expansion and subsequent decline in northeastern and western Europe was detected in the last five centuries. The actual distribution of walnut in Europe resulted from the combined effects of expansion/contraction from multiple refugia after the Last Glacial Maximum and its human exploitation over the last 5,000 years.

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.

Indirect Evidence for Genetic Differentiation in Vulnerability to Embolism in Pinus halepensis

Climate change is increasing mean temperatures and in the eastern Mediterranean is expected to decrease annual precipitation. The resulting increase in aridity may be too rapid for adaptation of tree species unless their gene pool already possesses variation in drought resistance. Vulnerability to embolism, estimated by the pressure inducing 50% loss of xylem hydraulic conductivity (P 50), is strongly associated with drought stress resistance in trees. Yet, previous studies on various tree species reported low intraspecific genetic variation for this trait, and therefore limited adaptive capacities to increasing aridity. Here we quantified differences in hydraulic efficiency (xylem hydraulic conductance) and safety (resistance to embolism) in four contrasting provenances of Pinus halepensis (Aleppo pine) in a provenance trial, which is indirect evidence for genetic differences. Results obtained with three techniques (bench dehydration, centrifugation and X-ray micro-CT) evidenced significant differentiation with similar ranking between provenances. Inter-provenance variation in P 50 correlated with pit anatomical properties (torus overlap and pit aperture size). These results suggest that adaptation of P. halepensis to xeric habitats has been accompanied by modifications of bordered pit function driven by variation in pit aperture. This study thus provides evidence that appropriate exploitation of provenance differences will allow continued forestry with P. halepensis in future climates of the Eastern Mediterranean.

Local adaptation to temperature and precipitation in naturally fragmented populations of Cephalotaxus oliveri, an endangered conifer endemic to China

Cephalotaxus oliveri is an endangered tertiary relict conifer endemic to China. The species survives in a wide range from west to east with heterogeneous climatic conditions. Precipitation and temperature are main restrictive factors for distribution of C. oliveri. In order to comprehend the mechanism of adaptive evolution to climate variation, we employed ISSR markers to detect adaptive evolution loci, to identify the association between variation in temperature and precipitation and adaptive loci, and to investigate the genetic structure for 22 C. oliveri natural populations. In total, 14 outlier loci were identified, of which five were associated with temperature and precipitation. Among outlier loci, linkage disequilibrium (LD) was high (42.86%), which also provided strong evidence for selection. In addition, C. oliveri possessed high genetic variation (93.31%) and population differentiation, which may provide raw material to evolution and accelerate local adaptation, respectively. Ecological niche modeling showed that global warming will cause a shift for populations of C. oliveri from south to north with a shrinkage of southern areas. Our results contribute to understand the potential response of conifers to climatic changes, and provide new insights for conifer resource management and conservation strategies.

Intraspecific variation in the use of water sources by the circum-Mediterranean conifer Pinus halepensis

The relevance of interspecific variation in the use of plant water sources has been recognized in drought-prone environments. By contrast, the characterization of intraspecific differences in water uptake patterns remains elusive, although preferential access to particular soil layers may be an important adaptive response for species along aridity gradients. Stable water isotopes were analysed in soil and xylem samples of 56 populations of the drought-avoidant conifer Pinus halepensis grown in a common garden test. We found that most populations reverted to deep soil layers as the main plant water source during seasonal summer droughts. More specifically, we detected a clear geographical differentiation among populations in water uptake patterns even under relatively mild drought conditions (early autumn), with populations originating from more arid regions taking up more water from deep soil layers. However, the preferential access to deep soil water was largely independent of aboveground growth. Our findings highlight the high plasticity and adaptive relevance of the differential access to soil water pools among Aleppo pine populations. The observed ecotypic patterns point to the adaptive relevance of resource investment in deep roots as a strategy towards securing a source of water in dry environments for P. halepensis.

Climatic adaptation and ecological divergence between two closely related pine species in Southeast China

Climate is one of the most important drivers for adaptive evolution in forest trees. Climatic selection contributes greatly to local adaptation and intraspecific differentiation, but this kind of selection could also have promoted interspecific divergence through ecological speciation. To test this hypothesis, we examined intra- and interspecific genetic variation at 25 climate-related candidate genes and 12 reference loci in two closely related pine species, Pinus massoniana Lamb. and Pinus hwangshanensis Hisa, using population genetic and landscape genetic approaches. These two species occur in Southeast China but have contrasting ecological preferences in terms of several environmental variables, notably altitude, although hybrids form where their distributions overlap. One or more robust tests detected signals of recent and/or ancient selection at two-thirds (17) of the 25 candidate genes, at varying evolutionary timescales, but only three of the 12 reference loci. The signals of recent selection were species specific, but signals of ancient selection were mostly shared by the two species likely because of the shared evolutionary history. FST outlier analysis identified six SNPs in five climate-related candidate genes under divergent selection between the two species. In addition, a total of 24 candidate SNPs representing nine candidate genes showed significant correlation with altitudinal divergence in the two species based on the covariance matrix of population history derived from reference SNPs. Genetic differentiation between these two species was higher at the candidate genes than at the reference loci. Moreover, analysis using the isolation-with-migration model indicated that gene flow between the species has been more restricted for climate-related candidate genes than the reference loci, in both directions. Taken together, our results suggest that species-specific and divergent climatic selection at the candidate genes might have counteracted interspecific gene flow and played a key role in the ecological divergence of these two closely related pine species.

First insights into the transcriptome and development of new genomic tools of a widespread circum-Mediterranean tree species, Pinus halepensis Mill

Aleppo pine (Pinus halepensis Mill.) is a relevant conifer species for studying adaptive responses to drought and fire regimes in the Mediterranean region. In this study, we performed Illumina next-generation sequencing of two phenotypically divergent Aleppo pine accessions with the aims of (i) characterizing the transcriptome through Illumina RNA-Seq on trees phenotypically divergent for adaptive traits linked to fire adaptation and drought, (ii) performing a functional annotation of the assembled transcriptome, (iii) identifying genes with accelerated evolutionary rates, (iv) studying the expression levels of the annotated genes and (v) developing gene-based markers for population genomic and association genetic studies. The assembled transcriptome consisted of 48,629 contigs and covered about 54.6 Mbp. The comparison of Aleppo pine transcripts to Picea sitchensis protein-coding sequences resulted in the detection of 34,014 SNPs across species, with a Ka /Ks average value of 0.216, suggesting that the majority of the assembled genes are under negative selection. Several genes were differentially expressed across the two pine accessions with contrasted phenotypes, including a glutathione-s-transferase, a cellulose synthase and a cobra-like protein. A large number of new markers (3334 amplifiable SSRs and 28,236 SNPs) have been identified which should facilitate future population genomics and association genetics in this species. A 384-SNP Oligo Pool Assay for genotyping with the Illumina VeraCode technology has been designed which showed an high overall SNP conversion rate (76.6%). Our results showed that Illumina next-generation sequencing is a valuable technology to obtain an extensive overview on whole transcriptomes of nonmodel species with large genomes.

Decoding the massive genome of loblolly pine using haploid DNA and novel assembly strategies

Background

The size and complexity of conifer genomes has, until now, prevented full genome sequencing and assembly. The large research community and economic importance of loblolly pine, Pinus taeda L., made it an early candidate for reference sequence determination.

Results

We develop a novel strategy to sequence the genome of loblolly pine that combines unique aspects of pine reproductive biology and genome assembly methodology. We use a whole genome shotgun approach relying primarily on next generation sequence generated from a single haploid seed megagametophyte from a loblolly pine tree, 20-1010, that has been used in industrial forest tree breeding. The resulting sequence and assembly was used to generate a draft genome spanning 23.2 Gbp and containing 20.1 Gbp with an N50 scaffold size of 66.9 kbp, making it a significant improvement over available conifer genomes. The long scaffold lengths allow the annotation of 50,172 gene models with intron lengths averaging over 2.7 kbp and sometimes exceeding 100 kbp in length. Analysis of orthologous gene sets identifies gene families that may be unique to conifers. We further characterize and expand the existing repeat library based on the de novo analysis of the repetitive content, estimated to encompass 82% of the genome.

Conclusions

In addition to its value as a resource for researchers and breeders, the loblolly pine genome sequence and assembly reported here demonstrates a novel approach to sequencing the large and complex genomes of this important group of plants that can now be widely applied.

The number of markers and samples needed for detecting bottlenecks under realistic scenarios, with and without recovery: a simulation-based study

Detecting bottlenecks is a common task in molecular ecology. While several bottleneck detection method sexist, evaluations of their power have focused only on severe bottlenecks (e.g. to Ne ~10). As a component of a recent review, Peery et al. (2012) analysed the power of two approaches, the M-ratio and heterozygote excess tests, to detect moderate bottlenecks (e.g. to Ne ~100),which is realistic for many conservation situations. In this Comment, we address three important points relevant to but not considered in Peery et al. Under moderate bottleneck scenarios, we test the (i) relative advantage of sampling more markers vs. more individuals, (ii) potential power to detect the bottleneck when utilizing dozens of microsatellites (a realistic possibility for contemporary studies) and (iii) reduction in power when post bottle neck recovery has occurred. For the realistic situations examined,we show that (i) doubling the number of loci shows equal or better power than tripling the number of individuals,(ii) increasing the number of markers (up to 100) results in continued additive gains in power, and (iii)recovery after a moderate amount of time or gradual change in size reduces power, by up to one-half. Our results provide a practical supplement to Peery et al. and encourage the continued use of bottleneck detection methods in the genomic age, but also emphasize that the power under different sampling schemes should be estimated,using simulation modelling, as a routine component of molecular ecology studies.

Fire structures pine serotiny at different scales

PREMISE OF THE STUDY

Serotiny (delayed seed release with the consequent accumulation of a canopy seedbank) confers fitness benefits in environments with crown-fire regimes. Thus, we predicted that serotiny level should be higher in populations recurrently subjected to crown-fires than in populations where crown-fires are rare. In addition, under a high frequency of fires, space and resources are recurrently available, permitting recruitment around each mother to follow the seed rain shadow. Thus, we also predicted spatial aggregation of serotiny within populations.

METHODS

We compared serotiny, considering both the proportion and the age of serotinous cones, in populations living in contrasting fire regimes for two iconic Mediterranean pine species (Pinus halepensis, P. pinaster). We framed our results by quantitatively comparing the strength of the fire-serotiny relationship with previous studies worldwide.

KEY RESULTS

For the two species, populations living under high crown-fire recurrence regimes had a higher serotiny level than those populations where the recurrence of crown-fires was low. For P. halepensis (the species with higher serotiny), populations in high fire recurrence regimes had higher fine-scale spatial aggregation of serotiny than those inhabiting low fire recurrence systems. The strength of the observed fire-serotiny relationship in P. halepensis is among the highest in published literature.

CONCLUSIONS

Fire regime shapes serotiny level among populations, and in populations with high serotiny, recurrent fires maintain a significant spatial structure for this trait. Consequently, fire has long-term evolutionary implications at different scales, emphasizing its prominent role in shaping the ecology of pines.

Novel conserved segments are associated with differential expression patterns for Pinaceae dehydrins

Dehydrins are thought to play an essential role in the response, acclimation and tolerance to different abiotic stresses, such as cold and drought. These proteins have been classified into five groups according to the presence of conserved and repeated motifs in their amino acid sequence. Due to their putative functions in the response to stress, dehydrins have been often used as candidate genes in studies on population variability and local adaptation to environmental conditions. However, little is still known regarding the differential role played by such groups or the mechanism underlying their function. Based on the sequences corresponding to dehydrins available in public databases we have isolated eight different dehydrins from cDNA of Pinus pinaster. We have obtained also their genomic sequences and identified their intron/exon structure. Quantitative RT-PCR analysis of their expression pattern in needles, stems and roots during a severe and prolonged drought stress, similar to the ones trees must face in nature, is also reported. Additionally, we have identified two amino acid motifs highly conserved and repeated in Pinaceae dehydrins and absent in angiosperms, presumably related to the divergent expression profiles observed.

Contrasting patterns of nucleotide diversity for four conifers of Alpine European forests

A candidate gene approach was used to identify levels of nucleotide diversity and to identify genes departing from neutral expectations in coniferous species of the Alpine European forest. Twelve samples were collected from four species that dominate montane and subalpine forests throughout Europe: Abies alba Mill, Larix decidua Mill, Pinus cembra L., and Pinus mugo Turra. A total of 800 genes, originally resequenced in Pinus taeda L., were resequenced across 12 independent trees for each of the four species. Genes were assigned to two categories, candidate and control, defined through homology-based searches to Arabidopsis. Estimates of nucleotide diversity per site varied greatly between polymorphic candidate genes (range: 0.0004–0.1295) and among species (range: 0.0024–0.0082), but were within the previously established ranges for conifers. Tests of neutrality using stringent significance thresholds, performed under the standard neutral model, revealed one to seven outlier loci for each species. Some of these outliers encode proteins that are involved with plant stress responses and form the basis for further evolutionary enquiries.

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.

Effects of forest plantations on the genetic composition of conspecific native Aleppo pine populations

Afforestation is a common and widespread management practice throughout the world, yet its implications for the genetic diversity of native populations are still poorly understood. We examined the effect of Aleppo pine (Pinus halepensis) plantations on the genetic composition of nearby conspecific native populations. We focused on two native populations in Israel with different levels of isolation from the surrounding plantations and compared the genetic diversity of naturally established young trees within the native populations with that of local native adults, using nine nuclear microsatellite markers. We found that the genetic composition of the recruits was significantly different from that of local adults in both populations, with allelic frequency changes between generations that could not be ascribed to random drift, but rather to substantial gene flow from the surrounding planted Aleppo pine populations. The more isolated population experienced a lower gene-flow level (22%) than the less isolated population (49%). The genetic divergence between native populations at the adult-tree stage (F(st) = 0.32) was more than twice as high as that of the young trees naturally established around native adults (F(st) = 0.15). Our findings provide evidence for a rapid genetic homogenization process of native populations following the massive planting efforts in the last decades. These findings have important implications for forest management and nature conservation and constitute a warning sign for the risk of translocation of biota for local biodiversity.

Temperature-related natural selection in a wind-pollinated tree across regional and continental scales

Adaptive genetic variation is a key factor in evolutionary biology, but the detection of signatures of natural selection remains challenging in nonmodel organisms. We used a genome scan approach to detect signals of natural selection in the Black alder (Alnus glutinosa), a widespread wind-pollinated tree. Gene flow through pollen dispersal is believed to be high in this species, and we therefore expected to find a clear response to natural selection. In combination with two different landscape genetic approaches, we determined which environmental variables were most associated with the inferred selection. This analysis was performed on a regional scale (northern Belgium) and on a continental scale (Europe). Because climate-related differences are much more pronounced at the continental scale, we expected to find more selection-sensitive genetic markers across Europe than across northern Belgium. At both spatial scales, a substantial number of genetic loci were considered outliers, with respect to neutral expectations, and were therefore identified as selective. Based on results from our combined approach, four putative selective loci (or 2.5%) were recovered with high statistical support. Although these loci seemed to be associated with different environmental variables, they were mainly temperature-related. Our study demonstrates that the use of complementary methods in landscape genetics allows the discovery of selective loci which otherwise might stay hidden. In combination with a genome scan, the selective loci can be verified and the nature of the selection pressure can be identified.

High genetic diversity at the extreme range edge: nucleotide variation at nuclear loci in Scots pine (Pinus sylvestris L.) in Scotland

Nucleotide polymorphism at 12 nuclear loci was studied in Scots pine populations across an environmental gradient in Scotland, to evaluate the impacts of demographic history and selection on genetic diversity. At eight loci, diversity patterns were compared between Scottish and continental European populations. At these loci, a similar level of diversity (θ(sil)= ~0.01) was found in Scottish vs mainland European populations, contrary to expectations for recent colonization, however, less rapid decay of linkage disequilibrium was observed in the former (ρ=0.0086±0.0009, ρ=0.0245±0.0022, respectively). Scottish populations also showed a deficit of rare nucleotide variants (multi-locus Tajima's D=0.316 vs D=-0.379) and differed significantly from mainland populations in allelic frequency and/or haplotype structure at several loci. Within Scotland, western populations showed slightly reduced nucleotide diversity (π(tot)=0.0068) compared with those from the south and east (0.0079 and 0.0083, respectively) and about three times higher recombination to diversity ratio (ρ/θ=0.71 vs 0.15 and 0.18, respectively). By comparison with results from coalescent simulations, the observed allelic frequency spectrum in the western populations was compatible with a relatively recent bottleneck (0.00175 × 4N(e) generations) that reduced the population to about 2% of the present size. However, heterogeneity in the allelic frequency distribution among geographical regions in Scotland suggests that subsequent admixture of populations with different demographic histories may also have played a role.

Forest tree genomics: growing resources and applications

Over the past two decades, research in forest tree genomics has lagged behind that of model and agricultural systems. However, genomic research in forest trees is poised to enter into an important and productive phase owing to the advent of next-generation sequencing technologies, the enormous genetic diversity in forest trees and the need to mitigate the effects of climate change. Research on long-lived woody perennials is extending our molecular knowledge of complex life histories and adaptations to the environment - enriching a field that has traditionally drawn biological inference from a few short-lived herbaceous species.

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.

Gene movement and genetic association with regional climate gradients in California valley oak (Quercus lobata Née) in the face of climate change

Rapid climate change jeopardizes tree populations by shifting current climate zones. To avoid extinction, tree populations must tolerate, adapt, or migrate. Here we investigate geographic patterns of genetic variation in valley oak, Quercus lobata Née, to assess how underlying genetic structure of populations might influence this species' ability to survive climate change. First, to understand how genetic lineages shape spatial genetic patterns, we examine historical patterns of colonization. Second, we examine the correlation between multivariate nuclear genetic variation and climatic variation. Third, to illustrate how geographic genetic variation could interact with regional patterns of 21st Century climate change, we produce region-specific bioclimatic distributions of valley oak using Maximum Entropy (MAXENT) models based on downscaled historical (1971-2000) and future (2070-2100) climate grids. Future climatologies are based on a moderate-high (A2) carbon emission scenario and two different global climate models. Chloroplast markers indicate historical range-wide connectivity via colonization, especially in the north. Multivariate nuclear genotypes show a strong association with climate variation that provides opportunity for local adaptation to the conditions within their climatic envelope. Comparison of regional current and projected patterns of climate suitability indicates that valley oaks grow in distinctly different climate conditions in different parts of their range. Our models predict widely different regional outcomes from local displacement of a few kilometres to hundreds of kilometres. We conclude that the relative importance of migration, adaptation, and tolerance are likely to vary widely for populations among regions, and that late 21st Century conditions could lead to regional extinctions.

Aquaporins in the wild: natural genetic diversity and selective pressure in the PIP gene family in five Neotropical tree species

Background

Tropical trees undergo severe stress through seasonal drought and flooding, and the ability of these species to respond may be a major factor in their survival in tropical ecosystems, particularly in relation to global climate change. Aquaporins are involved in the regulation of water flow and have been shown to be involved in drought response; they may therefore play a major adaptive role in these species. We describe genetic diversity in the PIP sub-family of the widespread gene family of Aquaporins in five Neotropical tree species covering four botanical families.

Results

PIP Aquaporin subfamily genes were isolated, and their DNA sequence polymorphisms characterised in natural populations. Sequence data were analysed with statistical tests of standard neutral equilibrium and demographic scenarios simulated to compare with the observed results. Chloroplast SSRs were also used to test demographic transitions. Most gene fragments are highly polymorphic and display signatures of balancing selection or bottlenecks; chloroplast SSR markers have significant statistics that do not conform to expectations for population bottlenecks. Although not incompatible with a purely demographic scenario, the combination of all tests tends to favour a selective interpretation of extant gene diversity.

Conclusions

Tropical tree PIP genes may generally undergo balancing selection, which may maintain high levels of genetic diversity at these loci. Genetic variation at PIP genes may represent a response to variable environmental conditions.

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

Quaternary climatic fluctuations have left contrasting historical footprints on the neutral genetic diversity patterns of existing populations of different tree species. We should expect the demography, and consequently the neutral genetic structure, of taxa less tolerant to particular climatic extremes to be more sensitive to long-term climate fluctuations. We explore this hypothesis here by sampling all six pine species found in the Iberian Peninsula (2464 individuals, 105 populations), using a common set of chloroplast microsatellite markers, and by looking at the association between neutral genetic diversity and species-specific climatic requirements. We found large variation in neutral genetic diversity and structure among Iberian pines, with cold-enduring mountain species (Pinus uncinata, P. sylvestris and P. nigra) showing substantially greater diversity than thermophilous taxa (P. pinea and P. halepensis). Within species, we observed a significant positive correlation between population genetic diversity and summer precipitation for some of the mountain pines. The observed pattern is consistent with the hypotheses that: (i) more thermophilous species have been subjected to stronger demographic fluctuations in the past, as a consequence of their maladaptation to recurrent glacial cold stages; and (ii) altitudinal migrations have allowed the maintenance of large effective population sizes and genetic variation in cold-tolerant species, especially in more humid regions. In the light of these results and hypotheses, we discuss some potential genetic consequences of impending climate change.