Genetic diversity, structure, and demographic change in tanoak, Lithocarpus densiflorus (Fagaceae), the most susceptible species to sudden oak death in California

Selection of Restoration Material for Abies koreana Based on Its Genetic Diversity on Mt. Hallasan

The restoration of damaged or disrupted forests with genetically appropriate restoration planting material that can adapt to future environmental conditions will ensure the conservation of forest genetic resources. Abies koreana is endemic to the Republic of Korea, with declining populations under current environmental changes. In this study, we examined the genetic diversity of its largest population growing on Mt. Hallasan to determine the sampling size of planting material from the population that will ensure 95% coverage of alleles in the population. We evaluated the genetic diversity and spatial genetic structure of three subpopulations of A. koreana on Mt. Hallasan. A total of 456 samples were evaluated using 10 microsatellites. The observed heterozygosity and expected heterozygosity were 0.538 and 0.614 at the population level, respectively. The differences among the subpopulations accounted for 4% of the total variance. Intervals between individuals of the sample to be extracted were based on the two-target distance (5 and 10 m) inferred from the spatial genetic structure. Through random sampling methods considering the target distance, we showed that genetic diversity can be captured by obtaining at least 35 individuals in the population of A. koreana on Mt. Hallasan.

Innate Resistance and Phosphite Treatment Affect Both the Pathogen's and Host's Transcriptomes in the Tanoak-Phytophthora ramorum Pathosystem

Phosphites have been used to control Sudden Oak Death; however, their precise mode of action is not fully understood. To study the mechanism of action of phosphites, we conducted an inoculation experiment on two open-pollinated tanoak families, previously found to be partially resistant. Stems of treatment group individuals were sprayed with phosphite, and seven days later, distal leaves were inoculated with the Sudden Oak Death pathogen Phytophthora ramorum. Leaves from treated and untreated control plants were harvested before and seven days after inoculation, and transcriptomes of both host and pathogen were analyzed. We found that tanoak families differed in the presence of innate resistance (resistance displayed by untreated tanoak) and in the response to phosphite treatment. A set of expressed genes associated with innate resistance was found to overlap with an expressed gene set for phosphite-induced resistance. This observation may indicate that phosphite treatment increases the resistance of susceptible host plants. In addition, genes of the pathogen involved in detoxification were upregulated in phosphite-treated plants compared to phosphite-untreated plants. In summary, our RNA-Seq analysis supports a two-fold mode of action of phosphites, including a direct toxic effect on P. ramorum and an indirect enhancement of resistance in the tanoak host.

Variation in Susceptibility of Tanoak to the NA1 and EU1 Lineages of Phytophthora ramorum, the Cause of Sudden Oak Death

Phytophthora ramorum, the cause of sudden oak death (SOD), kills tanoak (Notholithocarpus densiflorus) trees in southwestern Oregon and California. Two lineages of P. ramorum are now found in wildland forests of Oregon (NA1 and EU1). In addition to the management of SOD in forest ecosystems, disease resistance could be used as a way to mitigate the impact of P. ramorum. The objectives of this study were to (i) characterize the variability in resistance of N. densiflorus among families using lesion length; (ii) determine whether lineage, isolate, family, or their interactions significantly affect variation in lesion length; and (iii) determine whether there are differences among isolates and among families in terms of lesion length. The parameters isolate nested within lineage (isolate[lineage]) and family × isolate(lineage) interaction explained the majority of the variation in lesion length. There was no significant difference between the NA1 and EU1 lineages in terms of mean lesion length; however, there were differences among the six isolates. Lesions on seedlings collected from surviving trees at infested sites were smaller, on average, than lesions of seedlings collected from trees at noninfested sites (P = 0.0064). The results indicate that there is potential to establish a breeding program for tanoak resistance to SOD and that several isolates of P. ramorum should be used in an artificial inoculation assay.

Promise and Pitfalls of Endemic Resistance for Cultural Resources Threatened by Phytophthora ramorum

Invasive forest pathogens can harm cultural, economic, and ecological resources. Here, we demonstrate the potential of endemic tree pathogen resistance in forest disease management using Phytophthora ramorum, cause of sudden oak death, in the context of management of tanoak (Notholithocarpus densiflorus), an ecologically unique and highly valued tree within Native American communities of northern California and southern Oregon in the United States. We surveyed resistance to P. ramorum on the Hoopa Valley Indian Reservation and Yurok Indian Reservation in a set of study sites with variable management intensities. Variation in resistance was found at all sites with similar mean and variation across stands, and resistance tended to have a random spatial distribution within stands but was not associated with previous stand management (thinning or prescribed fire) or structural characteristics such as tree density, basal area, or pairwise relatedness among study trees. These results did not suggest host, genetic, management, or environment interactions that could be easily leveraged into treatments to increase the prevalence of resistant trees. We applied epidemiological models to assess the potential application of endemic resistance in this system and to examine our assumption that in planta differences in lesion size-our measure of resistance-reflect linkages between mortality and transmission (resistance) versus reduced mortality with no change in transmission (tolerance). This assumption strongly influenced infection dynamics but changes in host populations-our conservation focus-was dependent on community-level variation in transmission. For P. ramorum, slowing mortality rates (whether by resistance or tolerance) conserves host resources when a second source of inoculum is present; these results are likely generalizable to pathogens with a broader host range. However, when the focal host is the sole source of inoculum, increasing tolerant individuals led to the greatest stand-level pathogen accumulation in our model. When seeking to use variation in mortality rates to affect conservation strategies, it is important to understand how these traits are linked with transmission because tolerance will be more useful for management in mixed-host stands that are already invaded, compared with single-host stands with low or no pathogen presence, where resistance will have the greatest conservation benefits.

Population-level genetic variation and climate change in a biodiversity hotspot

INTRODUCTION

Estimated future climate scenarios can be used to predict where hotspots of endemism may occur over the next century, but life history, ecological and genetic traits will be important in informing the varying responses within myriad taxa. Essential to predicting the consequences of climate change to individual species will be an understanding of the factors that drive genetic structure within and among populations. Here, I review the factors that influence the genetic structure of plant species in California, but are applicable elsewhere; existing levels of genetic variation, life history and ecological characteristics will affect the ability of an individual taxon to persist in the presence of anthropogenic change.

FACTORS INFLUENCING THE DISTRIBUTION OF GENETIC VARIATION

Persistence in the face of climate change is likely determined by life history characteristics: dispersal ability, generation time, reproductive ability, degree of habitat specialization, plant-insect interactions, existing genetic diversity and availability of habitat or migration corridors. Existing levels of genetic diversity in plant populations vary based on a number of evolutionary scenarios that include endemism, expansion since the last glacial maximum, breeding system and current range sizes.

REGIONAL PRIORITIES AND EXAMPLES

A number of well-documented examples are provided from the California Floristic Province. Some predictions can be made for the responses of plant taxa to rapid environmental changes based on geographic position, evolutionary history, existing genetic variation, and ecological amplitude.

CONCLUSIONS, SOLUTIONS AND RECOMMENDATIONS

The prediction of how species will respond to climate change will require a synthesis drawing from population genetics, geography, palaeontology and ecology. The important integration of the historical factors that have shaped the distribution and existing genetic structure of California's plant taxa will enable us to predict and prioritize the conservation of species and areas most likely to be impacted by rapid climate change, human disturbance and invasive species.

Divergence maintained by climatic selection despite recurrent gene flow: a case study of Castanopsis carlesii (Fagaceae)

Local adaptation to different environments has the potential to maintain divergence between populations despite recurrent gene flow and is an important driver for generating biological diversity. In this study, we investigate the role of adaptation in the maintenance of two parapatric varieties of a forest tree. We used sequence variation of chloroplastic DNA and restriction site-associated DNA to investigate the genetic structure of two varieties of Castanopsis carlesii in subtropical China and relate it to climatic variation. We used niche reconstruction methods to investigate niche differentiation between the two varieties and to estimate the past distribution of this species. A deep divergence was observed between the two varieties, but evidence of introgression and genetic admixture was detected in two phenotypically and geographically intermediate populations. Niche reconstruction suggests that the distribution of the two varieties was disjunct during periods of global cooling and that the two varieties occupy significantly different niches. The genetic structure was mainly driven by environmental factors, and 13 outlier loci under divergent selection were correlated with climatic variation. These results suggest that the two varieties evolved in allopatry and came back into secondary contact after the last glacial maximum and that they are an evolutionary example of divergence maintained by climatic selection despite recurrent gene flow.

Post-glacial phylogeography and evolution of a wide-ranging highly-exploited keystone forest tree, eastern white pine (Pinus strobus) in North America: single refugium, multiple routes

Background

Knowledge of the historical distribution and postglacial phylogeography and evolution of a species is important to better understand its current distribution and population structure and potential fate in the future, especially under climate change conditions, and conservation of its genetic resources. We have addressed this issue in a wide-ranging and heavily exploited keystone forest tree species of eastern North America, eastern white pine (Pinus strobus). We examined the range-wide population genetic structure, tested various hypothetical population history and evolutionary scenarios and inferred the location of glacial refugium and post-glacial recolonization routes. Our hypothesis was that eastern white pine survived in a single glacial refugium and expanded through multiple post-glacial recolonization routes.

Results

We studied the range-wide genetic diversity and population structure of 33 eastern white pine populations using 12 nuclear and 3 chloroplast microsatellite DNA markers. We used Approximate Bayesian Computation approach to test various evolutionary scenarios. We observed high levels of genetic diversity, and significant genetic differentiation (F_ST = 0.104) and population structure among eastern white pine populations across its range. A south to north trend of declining genetic diversity existed, consistent with repeated founder effects during post-glaciation migration northwards. We observed broad consensus from nuclear and chloroplast genetic markers supporting the presence of two main post-glacial recolonization routes that originated from a single southern refugium in the mid-Atlantic plain. One route gave rise to populations at the western margin of the species’ range in Minnesota and western Ontario. The second route gave rise to central-eastern populations, which branched into two subgroups: central and eastern. We observed minimal sharing of chloroplast haplotypes between recolonization routes but there was evidence of admixture between the western and west-central populations.

Conclusions

Our study reveals a single southern refugium, two recolonization routes and three genetically distinguishable lineages in eastern white pine that we suggest to be treated as separate Evolutionarily Significant Units. Like many wide-ranging North American species, eastern white pine retains the genetic signatures of post-glacial recolonization and evolution, and its contemporary population genetic structure reflects not just the modern distribution and effects of heavy exploitation but also routes northward from its glacial refugium.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0624-1) contains supplementary material, which is available to authorized users.

Intervarietal and intravarietal genetic structure in Douglas-fir: nuclear SSRs bring novel insights into past population demographic processes, phylogeography, and intervarietal hybridization

Douglas-fir (Pseudotsuga menziesii) is one of numerous wide-range forest tree species represented by subspecies/varieties, which hybridize in contact zones. This study examined the genetic structure of this North American conifer and its two hybridizing varieties, coastal and Rocky Mountain, at intervarietal and intravarietal level. The genetic structure was subsequently associated with the Pleistocene refugial history, postglacial migration and intervarietal hybridization/introgression. Thirty-eight populations from the USA and Canada were genotyped for 13 nuclear SSRs and analyzed with simulations and traditional population genetic structuring methods. Eight genetic clusters were identified. The coastal clusters embodied five refugial populations originating from five distinct refugia. Four coastal refugial populations, three from California and one from western Canada, diverged during the Pleistocene (56.9–40.1 ka). The three Rocky Mountain clusters reflected distinct refugial populations of three glacial refugia. For Canada, ice covered during the Last Glacial Maximum, we present the following three findings. (1) One refugial population of each variety was revealed in the north of the distribution range. Additional research including paleodata is required to support and determine whether both northern populations originated from cryptic refugia situated south or north of the ice-covered area. (2) An interplay between intravarietal gene flow of different refugial populations and intervarietal gene flow by hybridization and introgression was identified. (3) The Canadian hybrid zone displayed predominantly introgressants of the Rocky Mountain into the coastal variety. This study provides new insights into the complex Quaternary dynamics of this conifer essential for understanding its evolution (outside and inside the native range), adaptation to future climates and for forest management.

Scaling up from greenhouse resistance to fitness in the field for a host of an emerging forest disease

Forest systems are increasingly threatened by emergent, exotic diseases, yet management strategies for forest trees may be hindered by long generation times and scant background knowledge. We tested whether nursery disease resistance and growth traits have predictive value for the conservation of Notholithocarpus densiflorus, the host most susceptible to sudden oak death. We established three experimental populations to assess nursery growth and resistance to Phytophthora ramorum, and correlations between nursery-derived breeding values with seedling survival in a field disease trial. Estimates of nursery traits' heritability were low to moderate, with lowest estimates for resistance traits. Within the field trial, survival likelihood was increased in larger seedlings and decreased with the development of disease symptoms. The seed-parent family wide likelihood of survival was likewise correlated with family predictors for size and resistance to disease in 2nd year laboratory assays, though not resistance in 1st year leaf assays. We identified traits and seedling families with increased survivorship in planted tanoaks, and a framework to further identify seed parents favored for restoration. The additive genetic variation and seedling disease dynamics we describe hold promise to refine current disease models and expand the understanding of evolutionary dynamics of emergent infectious diseases in highly susceptible hosts.

Clonal growth and fine-scale genetic structure in tanoak (Notholithocarpus densiflorus: Fagaceae)

The combination of sprouting and reproduction by seed can have important consequences on fine-scale spatial distribution of genetic structure (SGS). SGS is an important consideration for species' restoration because it determines the minimum distance among seed trees to maximize genetic diversity while not prejudicing locally adapted genotypes. Local environmental conditions can be expected to influence levels of clonal spread and SGS, particularly in the case of disturbance regimes such as fire. Here, we characterize fine-scale genetic structure and clonal spread in tanoak from drier upland sites and more mesic lowland woodlands. Clonal spread was a significant mode of stand development, but spread was limited on average to about 5-6 m. Gene dispersal was decomposed into clonal and sexual components. The latter varied according to whether it was estimated from all ramets with the clonal component removed or for a single ramet per genet. We used the difference in these 2 estimates of gene dispersal as a measure of the effect of clonality on effective population size in this species. Although upland sites had a greater number of ramets per genet, most of the other indices computed were not significantly different. However, they tended to show greater heterozygote excess and shorter gene dispersal distances than the lowland sites. The average distance among inferred sibships on upland sites was approximately at the scale of maximum clonal range. This was not the case on lowland sites, where sibs were more dispersed. We recommend minimum distances among seed trees to avoid selecting clones and to maximize genetic diversity for restoration.

Isolation by elevation: genetic structure at neutral and putatively non-neutral loci in a dominant tree of subtropical forests, Castanopsis eyrei

Background

The distribution of genetic diversity among plant populations growing along elevational gradients can be affected by neutral as well as selective processes. Molecular markers used to study these patterns usually target neutral processes only, but may also be affected by selection. In this study, the effects of elevation and successional stage on genetic diversity of a dominant tree species were investigated controlling for neutrality of the microsatellite loci used.

Methodology/Principal Findings

Diversity and differentiation among 24 populations of Castanopsis eyrei from different elevations (251–920 m) and successional stages were analysed by eight microsatellite loci. We found that one of the loci (Ccu97H18) strongly deviated from a neutral model of differentiation among populations due to either divergent selection or hitchhiking with an unknown selected locus. The analysis showed that C. eyrei populations had a high level of genetic diversity within populations (A_R = 7.6, H_E = 0.82). Genetic variation increased with elevation for both the putatively selected locus Ccu97H18 and the neutral loci. At locus Ccu97H18 one allele was dominant at low elevations, which was replaced at higher elevations by an increasing number of other alleles. The level of genetic differentiation at neutral loci was similar to that of other Fagaceae species (F_ST = 0.032,  = 0.15). Population differentiation followed a model of isolation by distance but additionally, strongly significant isolation by elevation was found, both for neutral loci and the putatively selected locus.

Conclusions/Significance

The results indicate higher gene flow among similar elevational levels than across different elevational levels and suggest a selective influence of elevation on the distribution of genetic diversity in C. eyrei. The study underlines the importance to check the selective neutrality of marker loci in analyses of population structure.

Inferring population decline and expansion from microsatellite data: a simulation-based evaluation of the Msvar method

Reconstructing the demographic history of populations is a central issue in evolutionary biology. Using likelihood-based methods coupled with Monte Carlo simulations, it is now possible to reconstruct past changes in population size from genetic data. Using simulated data sets under various demographic scenarios, we evaluate the statistical performance of Msvar, a full-likelihood Bayesian method that infers past demographic change from microsatellite data. Our simulation tests show that Msvar is very efficient at detecting population declines and expansions, provided the event is neither too weak nor too recent. We further show that Msvar outperforms two moment-based methods (the M-ratio test and Bottleneck) for detecting population size changes, whatever the time and the severity of the event. The same trend emerges from a compilation of empirical studies. The latest version of Msvar provides estimates of the current and the ancestral population size and the time since the population started changing in size. We show that, in the absence of prior knowledge, Msvar provides little information on the mutation rate, which results in biased estimates and/or wide credibility intervals for each of the demographic parameters. However, scaling the population size parameters with the mutation rate and scaling the time with current population size, as coalescent theory requires, significantly improves the quality of the estimates for contraction but not for expansion scenarios. Finally, our results suggest that Msvar is robust to moderate departures from a strict stepwise mutation model.