Divergent selection and heterogeneous migration rates across the range of Sitka spruce (Picea sitchensis)

Metapopulation heterogeneities in host mobility, productivity, and immunocompetency always increase virulence and infectiousness

The epidemiology and evolution of diseases unfold in populations that are rarely homogeneous. Instead, hosts infected by pathogens often form metapopulations, in which local populations connected by the movement of hosts experience different demographic and epidemiological conditions. Here, we develop a general theory of the evolution of pathogens in heterogeneous metapopulations. We reveal the following key insights into the evolution of pathogen virulence and infectiousness: (1) When the mobility (movement rate), productivity (birth rate and carrying capacity), or immunocompetency (immunity-loss rate) differ among local populations, this variance always increases pathogen virulence and infectiousness (2) The increment of pathogen virulence caused by such heterogeneity is approximately proportional to the variance of the corresponding heterogeneous local conditions (3) This increment can be expressed as the covariance between the local selection pressures and the local reproductive values experienced by the pathogen (4) The reason why heterogeneity always increases pathogen virulence is explained by the positive correlation of local selection pressures with reproductive values (5) Combinations of multiple independent heterogeneities further increase virulence and infectiousness, even more so when their covariances are positive. Our key findings robustly hold for different epidemiological frameworks - including SI, SIS, SIR, and SIRS models, with both density- and frequency-dependent transmission as well as with superinfection. They provide insights into the risks of growing pathogen infectiousness in a world in which heterogeneity - caused, e.g., by the concentration of human populations in urban areas - is rising.

Monitoring of species' genetic diversity in Europe varies greatly and overlooks potential climate change impacts

Genetic monitoring of populations currently attracts interest in the context of the Convention on Biological Diversity but needs long-term planning and investments. However, genetic diversity has been largely neglected in biodiversity monitoring, and when addressed, it is treated separately, detached from other conservation issues, such as habitat alteration due to climate change. We report an accounting of efforts to monitor population genetic diversity in Europe (genetic monitoring effort, GME), the evaluation of which can help guide future capacity building and collaboration towards areas most in need of expanded monitoring. Overlaying GME with areas where the ranges of selected species of conservation interest approach current and future climate niche limits helps identify whether GME coincides with anticipated climate change effects on biodiversity. Our analysis suggests that country area, financial resources and conservation policy influence GME, high values of which only partially match species' joint patterns of limits to suitable climatic conditions. Populations at trailing climatic niche margins probably hold genetic diversity that is important for adaptation to changing climate. Our results illuminate the need in Europe for expanded investment in genetic monitoring across climate gradients occupied by focal species, a need arguably greatest in southeastern European countries. This need could be met in part by expanding the European Union's Birds and Habitats Directives to fully address the conservation and monitoring of genetic diversity.

Pattern of Adaptive Divergence in Zingiber kawagoii Hayata (Zingiberaceae) along a Narrow Latitudinal Range

Ecological and evolutionary processes linking adaptation to environment are related to species’ range shifts. In this study, we employed amplified-fragment-length-polymorphism-based genome scan methods to identify candidate loci among Zingiber kawagoii populations inhabiting varying environments distributed at low to middle elevations (143–1488 m) in a narrow latitudinal range (between 21.90 and 25.30° N). Here, we show evidence of selection driving the divergence of Z. kawagoii. Twenty-six F_ST outliers were detected, which were significantly correlated with various environmental variables. The allele frequencies of nine F_ST outliers were either positively or negatively correlated with the population mean F_ST. Using several independent approaches, we found environmental variables act in a combinatorial fashion, best explaining outlier genetic variation. Nonetheless, we found that adaptive divergence was affected mostly by annual temperature range, and it is significantly positively correlated with latitude and significantly negatively correlated with the population mean F_ST. This study addresses a latitudinal pattern of changes in annual temperature range (which ranged from 13.8 °C in the Lanyu population to 18.5 °C in the Wulai population) and emphasizes the pattern of latitudinal population divergence closely linked to the allele frequencies of adaptive loci, acting in a narrow latitudinal range. Our results also indicate environmentally dependent local adaptation for both leading- and trailing-edge populations.

Genetic Diversity and Differentiation of Pedunculate Oak (Quercus robur L.) Populations at the Southern Margin of Its Distribution Range—Implications for Conservation

Understanding intraspecific genetic variation is one of the principal requirements for the evaluation of tree species capacity to cope with intensive climatic changes, as well as designing long-term conservation programs. Herein, we evaluated the genetic diversity and genetic structure of seven pedunculate oak (Quercus robur L.) populations, located at the southern margin of its distribution range on the Balkan Peninsula (Serbia). The objective of the study was to propose future in situ conservation measures aimed at protection of pedunculate oak adaptive and neutral genetic diversity at the species rear-edge. Genetic diversity and structure were estimated using twelve highly polymorphic simple sequence repeat (SSR) markers. The mean expected heterozygosity (He) was 0.769, allelic richness (AR) 9.63, and private allelic richness (pAR) 0.79, indicating high genetic diversity in the studied populations. Genetic differentiation among the populations was low (Fst = 0.032). Structure analysis, the unweighted pair group method with arithmetic mean (UPGMA) showed the existence of two gene pools unrelated to the populations’ area of occurrence. Taking into consideration the results of the current study and previous conservation activities on the pedunculate oak in Serbia, as well as the importance of rear-edge populations in the long-term conservation of the species genetic diversity, we suggested establishing three additional gene conservation units for securing long-term sustainability of the species.

Endemic Juniperus Montane Species Facing Extinction Risk under Climate Change in Southwest China: Integrative Approach for Conservation Assessment and Prioritization

Simple Summary

Climate change is one of the most significant drivers of habitat loss and species extinction, particularly montane endemic species such as Juniper trees, which are restricted to unique habitats. Therefore, assessing the impact of climate change on the extinction risk of species is a promising tool or guide for species conservation planning. The loss in species habitat due to global warming indicates the level of extinction or endangerment. Predictions of suitable habitats are outputs from assessment analysis. This will help conservationists discover new populations of endemic species and help raise the awareness of local people to save and rescue these endangered species.

Abstract

Climate change is an important driver of biodiversity loss and extinction of endemic montane species. In China, three endemic Juniperus spp. (Juniperuspingii var. pingii, J.tibetica, and J.komarovii) are threatened and subjected to the risk of extinction. This study aimed to predict the potential distribution of these three Juniperus species under climate change and dispersal scenarios, to identify critical drivers explaining their potential distributions, to assess the extinction risk by estimating the loss percentage in their area of occupancy (AOO), and to identify priority areas for their conservation in China. We used ensemble modeling to evaluate the impact of climate change and project AOO. Our results revealed that the projected AOOs followed a similar trend in the three Juniperus species, which predicted an entire loss of their suitable habitats under both climate and dispersal scenarios. Temperature annual range and isothermality were the most critical key variables explaining the potential distribution of these three Juniperus species; they contribute by 16–56.1% and 20.4–38.3%, respectively. Accounting for the use of different thresholds provides a balanced approach for species distribution models’ applications in conservation assessment when the goal is to assess potential climatic suitability in new geographical areas. Therefore, south Sichuan and north Yunnan could be considered important priority conservation areas for in situ conservation and search for unknown populations of these three Juniperus species.

Modern Strategies to Assess and Breed Forest Tree Adaptation to Changing Climate

Studying the genetics of adaptation to new environments in ecologically and industrially important tree species is currently a major research line in the fields of plant science and genetic improvement for tolerance to abiotic stress. Specifically, exploring the genomic basis of local adaptation is imperative for assessing the conditions under which trees will successfully adapt in situ to global climate change. However, this knowledge has scarcely been used in conservation and forest tree improvement because woody perennials face major research limitations such as their outcrossing reproductive systems, long juvenile phase, and huge genome sizes. Therefore, in this review we discuss predictive genomic approaches that promise increasing adaptive selection accuracy and shortening generation intervals. They may also assist the detection of novel allelic variants from tree germplasm, and disclose the genomic potential of adaptation to different environments. For instance, natural populations of tree species invite using tools from the population genomics field to study the signatures of local adaptation. Conventional genetic markers and whole genome sequencing both help identifying genes and markers that diverge between local populations more than expected under neutrality, and that exhibit unique signatures of diversity indicative of "selective sweeps." Ultimately, these efforts inform the conservation and breeding status capable of pivoting forest health, ecosystem services, and sustainable production. Key long-term perspectives include understanding how trees' phylogeographic history may affect the adaptive relevant genetic variation available for adaptation to environmental change. Encouraging "big data" approaches (machine learning-ML) capable of comprehensively merging heterogeneous genomic and ecological datasets is becoming imperative, too.

Long-distance pollen dispersal during recent colonization favors a rapid but partial recovery of genetic diversity in Picea sitchensis

Tree species in the northern hemisphere are currently subject to rapid anthropogenic climate change and are shifting their ranges in response. This prompts questions about the mechanisms allowing tree populations to respond quickly to selection pressures when establishing into new areas. Focusing on the northern expanding range edge of Picea sitchensis, a widespread conifer of western North America, we ask how genetic structure and diversity develop during colonization, and assess the role of demographic history in shaping the evolutionary trajectory of an establishing population. We combined 500 yr of tree-ring and genotyping-by-sequencing data in 639 trees at the expansion front on the Kodiak Archipelago. We show that alleles accumulated rapidly during an increase in recruitment rate in the early 1700s. A shift from foreign to local pollen flow subsequently homogenized genetic structure at the expansion front. Taking advantage of the exceptional longevity of conifers, we highlight the major role of long-distance pollen dispersal in the rapid but incomplete recovery of genetic diversity during the initial stages of colonization. We also warn that slow initial population growth as well as long-lasting dominance of local gene flow by early founders could increase evolutionary load under a rapidly changing climate.

Gene flow improves fitness at a range edge under climate change

Populations at the margins of a species' geographic range are often thought to be poorly adapted to their environment. According to theoretical predictions, gene flow can inhibit these range edge populations if it disrupts adaptation to local conditions. Alternatively, if range edge populations are small or isolated, gene flow can provide beneficial genetic variation and may facilitate adaptation to environmental change. We tested these competing predictions in the annual wildflower Clarkia pulchella using greenhouse crosses to simulate gene flow from sources across the geographic range into two populations at the northern range margin. We planted these between-population hybrids in common gardens at the range edge and evaluated how genetic differentiation and climatic differences between edge populations and gene flow sources affected lifetime fitness. During an anomalously warm study year, gene flow from populations occupying historically warm sites improved fitness at the range edge and plants with one or both parents from warm populations performed best. The effects of the temperature provenance of gene flow sources were most apparent at early life history stages, but precipitation provenance also affected reproduction. We also found benefits of gene flow that were independent of climate: after climate was controlled for, plants with parents from different populations performed better at later lifestages than those with parents from the same population, indicating that gene flow may improve fitness via relieving homozygosity. Further supporting this result, we found that increasing genetic differentiation of parental populations had positive effects on fitness of hybrid seeds. Gene flow from warmer populations, when it occurs, is likely to contribute adaptive genetic variation to populations at the northern range edge as the climate warms. On heterogeneous landscapes, climate of origin may be a better predictor of gene flow effects than geographic proximity.

Invoking adaptation to decipher the genetic legacy of past climate change

Persistence of natural populations during periods of climate change is likely to depend on migration (range shifts) or adaptation. These responses were traditionally considered discrete processes and conceptually divided into the realms of ecology and evolution. In a milestone paper, Davis and Shaw (2001) Science 292:673 argued that the interplay of adaptation and migration was central to biotic responses to Quaternary climate, but since then there has been no synthesis of efforts made to set up this research program. Here we review some of the salient findings from molecular genetic studies assessing ecological and evolutionary responses to Quaternary climate change. These studies have revolutionized our understanding of population processes associated with past species migration. However, knowledge remains limited about the role of natural selection for local adaptation of populations to Quaternary environmental fluctuations and associated range shifts, and for the footprints this might have left on extant populations. Next-generation sequencing technologies, high-resolution paleoclimate analyses, and advances in population genetic theory offer an unprecedented opportunity to test hypotheses about adaptation through time. Recent population genomics studies have greatly improved our understanding of the role of contemporary adaptation to local environments in shaping spatial patterns of genetic diversity across modern-day landscapes. Advances in this burgeoning field provide important conceptual and methodological bases to decipher the historical role of natural selection and assess adaptation to past environmental variation. We suggest that a process called "temporal conditional neutrality" has taken place: some alleles favored in glacial environments become selectively neutral in modern-day conditions, whereas some alleles that had been neutral during glacial periods become under selection in modern environments. Building on this view, we present a new integrative framework for addressing the interplay of demographic and adaptive evolutionary responses to Quaternary climate dynamics, the research agenda initially envisioned by Davis and Shaw (2001) Science 292:673.

Transcriptome analysis of Pinus halepensis under drought stress and during recovery

Forest trees use various strategies to cope with drought stress and these strategies involve complex molecular mechanisms. Pinus halepensis Miller (Aleppo pine) is found throughout the Mediterranean basin and is one of the most drought-tolerant pine species. In order to decipher the molecular mechanisms that P. halepensis uses to withstand drought, we performed large-scale physiological and transcriptome analyses. We selected a mature tree from a semi-arid area with suboptimal growth conditions for clonal propagation through cuttings. We then used a high-throughput experimental system to continuously monitor whole-plant transpiration rates, stomatal conductance and the vapor pressure deficit. The transcriptomes of plants were examined at six physiological stages: pre-stomatal response, partial stomatal closure, minimum transpiration, post-irrigation, partial recovery and full recovery. At each stage, data from plants exposed to the drought treatment were compared with data collected from well-irrigated control plants. A drought-stressed P. halepensis transcriptome was created using paired-end RNA-seq. In total, ~6000 differentially expressed, non-redundant transcripts were identified between drought-treated and control trees. Cluster analysis has revealed stress-induced down-regulation of transcripts related to photosynthesis, reactive oxygen species (ROS)-scavenging through the ascorbic acid (AsA)-glutathione cycle, fatty acid and cell wall biosynthesis, stomatal activity, and the biosynthesis of flavonoids and terpenoids. Up-regulated processes included chlorophyll degradation, ROS-scavenging through AsA-independent thiol-mediated pathways, abscisic acid response and accumulation of heat shock proteins, thaumatin and exordium. Recovery from drought induced strong transcription of retrotransposons, especially the retrovirus-related transposon Tnt1-94. The drought-related transcriptome illustrates this species' dynamic response to drought and recovery and unravels novel mechanisms.

The role of niche divergence and geographic arrangement in the speciation of Eared Pheasants (Crossoptilon, Hodgson 1938)

One of the most contentious theories in current ecology is the ecological niche conservatism, which is defined as conservatism among closely related species; however, the ecological niche can also be shifted, as documented in several cases. Genetic drift and ecological divergent selection may cause ecological niche divergence. The current study aims to test whether the ecological niche is conserved or divergent and to determine the main factor that drives ecological niche divergence or conservation. We analyzed the phylogenetic relationship, ecological niche model (ENM) and demographic history of Eared Pheasants in the genus Crossoptilon (Galliformes: Phasianidae) to test niche conservatism with respect to different geographically distributed patterns. The phylogenetic relationship was reconstructed using ∗BEAST with mitochondrial cytochrome b (cyt b) and 44 unlinked autosomal exonic loci, and ENMs were reconstructed in MAXENT using an average of 41 occurrence sites in each species and 22 bioclimatic variables. A background similarity test was used to detect whether the ecological niche is conserved. Demographic history was estimated using the isolation with migration (IM) model. We found that there was asymmetric gene flow between the allopatric sister species Crossoptilon mantchuricum and C. auritum and the parapatric sister species C. harmani and C. crossoptilon. We found that ecological niches were divergent, not conserved, between C. mantchuricum and C. auritum, which began to diverge at approximately 0.3 million years ago. However, the ecological niches were conserved between C. crossoptilon and C. harmani, which gradually diverged approximately half a million years ago. Ecological niches can be either conserved or divergent, and ecological divergent selection for local adaptation is probably an important factor that promotes and maintains niche divergence in the face of gene flow. This study provides a better understanding of the role that divergent selection has in the initial speciation process. The platform combined demographic processes and ecological niches to offer new insights into the mechanism of biogeography patterns.

Genetic Structure in the Northern Range Margins of Common Ash, Fraxinus excelsior L

During post glacial colonization, loss of genetic diversity due to leading edge effects may be attenuated in forest trees because of their prolonged juvenile phase, allowing many migrants to reach the colonizing front before populations become reproductive. The northern range margins of temperate tree taxa in Europe are particularly suitable to study the genetic processes that follow colonization because they have been little affected by northern refugia. Here we examined how post glacial range dynamics have shaped the genetic structure of common ash (Fraxinus excelsior L.) in its northern range compared to its central range in Europe. We used four chloroplast and six nuclear microsatellites to screen 42 populations (1099 trees), half of which corresponded to newly sampled populations in the northern range and half of which represented reference populations from the central range obtained from previously studies. We found that northern range populations of common ash have the same chloroplast haplotypes as south-eastern European populations, suggesting that colonization of the northern range took place along a single migration route, a result confirmed by the structure at the nuclear microsatellites. Along this route, diversity strongly decreased only in the northern range, concomitantly with increasing population differentiation and complex population substructures, a pattern consistent with a leading edge colonization model. Our study highlights that while diversity is maintained in the central range of common ash due to broad colonizing fronts and high levels of gene flow, it profoundly decreases in the northern range, where colonization was unidirectional and probably involved repeated founder events and population fluctuations. Currently, common ash is threatened by ash dieback, and our results on northern populations will be valuable for developing gene conservation strategies.

Genome-environment association study suggests local adaptation to climate at the regional scale in Fagus sylvatica

The evolutionary potential of long-lived species, such as forest trees, is fundamental for their local persistence under climate change (CC). Genome-environment association (GEA) analyses reveal if species in heterogeneous environments at the regional scale are under differential selection resulting in populations with potential preadaptation to CC within this area. In 79 natural Fagus sylvatica populations, neutral genetic patterns were characterized using 12 simple sequence repeat (SSR) markers, and genomic variation (144 single nucleotide polymorphisms (SNPs) out of 52 candidate genes) was related to 87 environmental predictors in the latent factor mixed model, logistic regressions and isolation by distance/environmental (IBD/IBE) tests. SSR diversity revealed relatedness at up to 150 m intertree distance but an absence of large-scale spatial genetic structure and IBE. In the GEA analyses, 16 SNPs in 10 genes responded to one or several environmental predictors and IBE, corrected for IBD, was confirmed. The GEA often reflected the proposed gene functions, including indications for adaptation to water availability and temperature. Genomic divergence and the lack of large-scale neutral genetic patterns suggest that gene flow allows the spread of advantageous alleles in adaptive genes. Thereby, adaptation processes are likely to take place in species occurring in heterogeneous environments, which might reduce their regional extinction risk under CC.

Time to get moving: assisted gene flow of forest trees

Geographic variation in trees has been investigated since the mid-18th century. Similar patterns of clinal variation have been observed along latitudinal and elevational gradients in common garden experiments for many temperate and boreal species. These studies convinced forest managers that a 'local is best' seed source policy was usually safest for reforestation. In recent decades, experimental design, phenotyping methods, climatic data and statistical analyses have improved greatly and refined but not radically changed knowledge of clines. The maintenance of local adaptation despite high gene flow suggests selection for local adaptation to climate is strong. Concerns over maladaptation resulting from climate change have motivated many new genecological and population genomics studies; however, few jurisdictions have implemented assisted gene flow (AGF), the translocation of pre-adapted individuals to facilitate adaptation of planted forests to climate change. Here, we provide evidence that temperate tree species show clines along climatic gradients sufficiently similar for average patterns or climate models to guide AGF in the absence of species-specific knowledge. Composite provenancing of multiple seed sources can be used to increase diversity and buffer against future climate uncertainty. New knowledge will continue to refine and improve AGF as climates warm further.

Evolutionary Quantitative Genomics of Populus trichocarpa

Forest trees generally show high levels of local adaptation and efforts focusing on understanding adaptation to climate will be crucial for species survival and management. Here, we address fundamental questions regarding the molecular basis of adaptation in undomesticated forest tree populations to past climatic environments by employing an integrative quantitative genetics and landscape genomics approach. Using this comprehensive approach, we studied the molecular basis of climate adaptation in 433 Populus trichocarpa (black cottonwood) genotypes originating across western North America. Variation in 74 field-assessed traits (growth, ecophysiology, phenology, leaf stomata, wood, and disease resistance) was investigated for signatures of selection (comparing Q_ST -F_ST) using clustering of individuals by climate of origin (temperature and precipitation). 29,354 SNPs were investigated employing three different outlier detection methods and marker-inferred relatedness was estimated to obtain the narrow-sense estimate of population differentiation in wild populations. In addition, we compared our results with previously assessed selection of candidate SNPs using the 25 topographical units (drainages) across the P. trichocarpa sampling range as population groupings. Narrow-sense Q_ST for 53% of distinct field traits was significantly divergent from expectations of neutrality (indicating adaptive trait variation); 2,855 SNPs showed signals of diversifying selection and of these, 118 SNPs (within 81 genes) were associated with adaptive traits (based on significant Q_ST). Many SNPs were putatively pleiotropic for functionally uncorrelated adaptive traits, such as autumn phenology, height, and disease resistance. Evolutionary quantitative genomics in P. trichocarpa provides an enhanced understanding regarding the molecular basis of climate-driven selection in forest trees and we highlight that important loci underlying adaptive trait variation also show relationship to climate of origin. We consider our approach the most comprehensive, as it uncovers the molecular mechanisms of adaptation using multiple methods and tests. We also provide a detailed outline of the required analyses for studying adaptation to the environment in a population genomics context to better understand the species’ potential adaptive capacity to future climatic scenarios.

Genomic patterns of species diversity and divergence in Eucalyptus

We examined genome-wide patterns of DNA sequence diversity and divergence among six species of the important tree genus Eucalyptus and investigated their relationship with genomic architecture. Using c. 90 range-wide individuals of each Eucalyptus species (E. grandis, E. urophylla, E. globulus, E. nitens, E. dunnii and E. camaldulensis), genetic diversity and divergence were estimated from 2840 polymorphic diversity arrays technology markers covering the 11 chromosomes. Species differentiating markers (SDMs) identified in each of 15 pairwise species comparisons, along with species diversity (HHW ) and divergence (FST ), were projected onto the E. grandis reference genome. Across all species comparisons, SDMs totalled 1.1-5.3% of markers and were widely distributed throughout the genome. Marker divergence (FST and SDMs) and diversity differed among and within chromosomes. Patterns of diversity and divergence were broadly conserved across species and significantly associated with genomic features, including the proximity of markers to genes, the relative number of clusters of tandem duplications, and gene density within or among chromosomes. These results suggest that genomic architecture influences patterns of species diversity and divergence in the genus. This influence is evident across the six species, encompassing diverse phylogenetic lineages, geography and ecology.

Holocene re-colonisation, central-marginal distribution and habitat specialisation shape population genetic patterns within an Atlantic European grass species

Corynephorus canescens (L.) P.Beauv. is an outbreeding, short-lived and wind-dispersed grass species, highly specialised on scattered and disturbance-dependent habitats of open sandy sites. Its distribution ranges from the Iberian Peninsula over Atlantic regions of Western and Central Europe, but excludes the two other classical European glacial refuge regions on the Apennine and Balkan Peninsulas. To investigate genetic patterns of this uncommon combination of ecological and biogeographic species characteristics, we analysed AFLP variation among 49 populations throughout the European distribution range, expecting (i) patterns of SW European glacial refugia and post-glacial expansion to the NE; (ii) decreasing genetic diversity from central to marginal populations; and (iii) interacting effects of high gene flow and disturbance-driven genetic drift. Decreasing genetic diversity from SW to NE and distinct gene pool clustering imply refugia on the Iberian Peninsula and in western France, from where range expansion originated towards the NE. High genetic diversity within and moderate genetic differentiation among populations, and a significant pattern of isolation-by-distance indicate a gene flow drift equilibrium within C. canescens, probably due to its restriction to scattered and dynamic habitats and limited dispersal distances. These features, as well as the re-colonisation history, were found to affect genetic diversity gradients from central to marginal populations. Our study emphasises the need for including the specific ecology into analyses of species (re-)colonisation histories and range centre-margin analyses. To account for discontinuous distributions, new indices of marginality were tested for their suitability in studies of centre-periphery gradients.

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.

Evolutionary responses of tree phenology to the combined effects of assortative mating, gene flow and divergent selection

The timing of bud burst (TBB) in temperate trees is a key adaptive trait, the expression of which is triggered by temperature gradients across the landscape. TBB is strongly correlated with flowering time and is therefore probably mediated by assortative mating. We derived theoretical predictions and realized numerical simulations of evolutionary changes in TBB in response to divergent selection and gene flow in a metapopulation. We showed that the combination of the environmental gradient of TBB and assortative mating creates contrasting genetic clines, depending on the direction of divergent selection. If divergent selection acts in the same direction as the environmental gradient (cogradient settings), genetic clines are established and inflated by assortative mating. Conversely, under divergent selection of the same strength but acting in the opposite direction (countergradient selection), genetic clines are slightly constrained. We explored the consequences of these dynamics for population maladaptation, by monitoring pollen swamping. Depending on the direction of divergent selection with respect to the environmental gradient, pollen filtering owing to assortative mating either facilitates or impedes adaptation in peripheral populations.

Genetic divergence and signatures of natural selection in marginal populations of a keystone, long-lived conifer, Eastern White Pine (Pinus strobus) from Northern Ontario

Marginal populations are expected to provide the frontiers for adaptation, evolution and range shifts of plant species under the anticipated climate change conditions. Marginal populations are predicted to show genetic divergence from central populations due to their isolation, and divergent natural selection and genetic drift operating therein. Marginal populations are also expected to have lower genetic diversity and effective population size (Ne) and higher genetic differentiation than central populations. We tested these hypotheses using eastern white pine (Pinus strobus) as a model for keystone, long-lived widely-distributed plants. All 614 eastern white pine trees, in a complete census of two populations each of marginal old-growth, central old-growth, and central second-growth, were genotyped at 11 microsatellite loci. The central populations had significantly higher allelic and genotypic diversity, latent genetic potential (LGP) and Ne than the marginal populations. However, heterozygosity and fixation index were similar between them. The marginal populations were genetically diverged from the central populations. Model testing suggested predominant north to south gene flow in the study area with curtailed gene flow to northern marginal populations. Signatures of natural selection were detected at three loci in the marginal populations; two showing divergent selection with directional change in allele frequencies, and one balancing selection. Contrary to the general belief, no significant differences were observed in genetic diversity, differentiation, LGP, and Ne between old-growth and second-growth populations. Our study provides information on the dynamics of migration, genetic drift and selection in central versus marginal populations of a keystone long-lived plant species and has broad evolutionary, conservation and adaptation significance.

High-throughput transcriptome sequencing and preliminary functional analysis in four Neotropical tree species

Background

The Amazonian rainforest is predicted to suffer from ongoing environmental changes. Despite the need to evaluate the impact of such changes on tree genetic diversity, we almost entirely lack genomic resources.

Results

In this study, we analysed the transcriptome of four tropical tree species (Carapa guianensis, Eperua falcata, Symphonia globulifera and Virola michelii) with contrasting ecological features, belonging to four widespread botanical families (respectively Meliaceae, Fabaceae, Clusiaceae and Myristicaceae). We sequenced cDNA libraries from three organs (leaves, stems, and roots) using 454 pyrosequencing. We have developed an R and bioperl-based bioinformatic procedure for de novo assembly, gene functional annotation and marker discovery. Mismatch identification takes into account single-base quality values as well as the likelihood of false variants as a function of contig depth and number of sequenced chromosomes. Between 17103 (for Symphonia globulifera) and 23390 (for Eperua falcata) contigs were assembled. Organs varied in the numbers of unigenes they apparently express, with higher number in roots. Patterns of gene expression were similar across species, with metabolism of aromatic compounds standing out as an overrepresented gene function. Transcripts corresponding to several gene functions were found to be over- or underrepresented in each organ. We identified between 4434 (for Symphonia globulifera) and 9076 (for Virola surinamensis) well-supported mismatches. The resulting overall mismatch density was comprised between 0.89 (S. globulifera) and 1.05 (V. surinamensis) mismatches/100 bp in variation-containing contigs.

Conclusion

The relative representation of gene functions in the four transcriptomes suggests that secondary metabolism may be particularly important in tropical trees. The differential representation of transcripts among tissues suggests differential gene expression, which opens the way to functional studies in these non-model, ecologically important species. We found substantial amounts of mismatches in the four species. These newly identified putative variants are a first step towards acquiring much needed genomic resources for tropical tree species.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-238) contains supplementary material, which is available to authorized users.

Evolutionary and plastic responses to climate change in terrestrial plant populations

As climate change progresses, we are observing widespread changes in phenotypes in many plant populations. Whether these phenotypic changes are directly caused by climate change, and whether they result from phenotypic plasticity or evolution, are active areas of investigation. Here, we review terrestrial plant studies addressing these questions. Plastic and evolutionary responses to climate change are clearly occurring. Of the 38 studies that met our criteria for inclusion, all found plastic or evolutionary responses, with 26 studies showing both. These responses, however, may be insufficient to keep pace with climate change, as indicated by eight of 12 studies that examined this directly. There is also mixed evidence for whether evolutionary responses are adaptive, and whether they are directly caused by contemporary climatic changes. We discuss factors that will likely influence the extent of plastic and evolutionary responses, including patterns of environmental changes, species' life history characteristics including generation time and breeding system, and degree and direction of gene flow. Future studies with standardized methodologies, especially those that use direct approaches assessing responses to climate change over time, and sharing of data through public databases, will facilitate better predictions of the capacity for plant populations to respond to rapid climate change.

Genetic and morphological structure of a spruce hybrid (Picea sitchensis x P. glauca) zone along a climatic gradient

PREMISE OF THE STUDY

Historic colonization and contemporary evolutionary processes contribute to patterns of genetic variation and differentiation among populations. However, separating the respective influences of these processes remains a challenge, particularly for natural hybrid zones, where standing genetic variation may result from evolutionary processes both preceding and following contact, influencing the evolutionary trajectory of hybrid populations. Where adaptation to novel environments may be facilitated by interspecific hybridization, teasing apart these processes will have practical implications for forest management in changing environments.

METHODS

We evaluated the neutral genetic architecture of the Picea sitchensis (Sitka spruce) × P. glauca (white spruce) hybrid zone along the Nass and Skeena river valleys in northwestern British Columbia using chloroplast, mitochondrial, and nuclear microsatellite markers, in combination with cone morphological traits.

KEY RESULTS

Sitka spruce mitotype "capture", evidenced by this species dominating the maternal lineage, is consistent with earlier colonization of the region by Sitka spruce. This "capture" differs from the spatial distribution of chloroplast haplotypes, indicating pollen dispersal and its contribution to geographic structure. Genetic ancestry, based on nuclear markers, was strongly influenced by climate and geography. Highly parallel results for replicate transects along environmental gradients provide support for the bounded hybrid superiority model of hybrid zone maintenance. •

CONCLUSIONS

This broad-scale analysis of neutral genetic structure indicates the importance of historic and contemporary gene flow, environmental selection, and their interaction in shaping neutral genetic variation within this hybrid zone, informative to seed transfer development and reforestation for future climates.

Historical connectivity, contemporary isolation and local adaptation in a widespread but discontinuously distributed species endemic to Taiwan, Rhododendron oldhamii (Ericaceae)

Elucidation of the evolutionary processes that constrain or facilitate adaptive divergence is a central goal in evolutionary biology, especially in non-model organisms. We tested whether changes in dynamics of gene flow (historical vs contemporary) caused population isolation and examined local adaptation in response to environmental selective forces in fragmented Rhododendron oldhamii populations. Variation in 26 expressed sequence tag-simple sequence repeat loci from 18 populations in Taiwan was investigated by examining patterns of genetic diversity, inbreeding, geographic structure, recent bottlenecks, and historical and contemporary gene flow. Selection associated with environmental variables was also examined. Bayesian clustering analysis revealed four regional population groups of north, central, south and southeast with significant genetic differentiation. Historical bottlenecks beginning 9168-13,092 years ago and ending 1584-3504 years ago were revealed by estimates using approximate Bayesian computation for all four regional samples analyzed. Recent migration within and across geographic regions was limited. However, major dispersal sources were found within geographic regions. Altitudinal clines of allelic frequencies of environmentally associated positively selected outliers were found, indicating adaptive divergence. Our results point to a transition from historical population connectivity toward contemporary population isolation and divergence on a regional scale. Spatial and temporal dispersal differences may have resulted in regional population divergence and local adaptation associated with environmental variables, which may have played roles as selective forces at a regional scale.