Reduced representation sequencing to understand the evolutionary history of Torrey pine (Pinus torreyana parry) with implications for rare species conservation.

Understanding the contribution of neutral and adaptive evolutionary processes to population differentiation is often necessary for better informed management and conservation of rare species. In this study, we focused on Pinus torreyana Parry (Torrey pine), one of the world's rarest pines, endemic to one island and one mainland population in California. Small population size, low genetic diversity, and susceptibility to abiotic and biotic stresses suggest Torrey pine may benefit from interpopulation genetic rescue to preserve the species' evolutionary potential. We leveraged reduced representation sequencing to tease apart the respective contributions of stochastic and deterministic evolutionary processes to population differentiation. We applied these data to model spatial and temporal demographic changes in effective population sizes and genetic connectivity, to identify loci possibly under selection, and evaluate genetic rescue as a potential conservation strategy. Overall, we observed exceedingly low standing variation within both Torrey pine populations, reflecting consistently low effective population sizes across time, and limited genetic differentiation, suggesting maintenance of gene flow between populations following divergence. However, genome scans identified more than 2000 candidate SNPs potentially under divergent selection. Combined with previous observations indicating population phenotypic differentiation, this indicates natural selection has probably contributed to the evolution of population genetic differences. Thus, while reduced genetic diversity, small effective population size, and genetic connectivity between populations suggest genetic rescue could mitigate the adverse effects of rarity, evidence for adaptive differentiation suggests genetic mixing could disrupt adaptation. Further work evaluating the fitness consequences of inter-population admixture is necessary to empirically evaluate the trade-offs associated with genetic rescue in Torrey pine.

SNPs in a Large Genomic Scaffold Are Strongly Associated with Cr1R, Major Gene for Resistance to White Pine Blister Rust in Range-Wide Samples of Sugar Pine (Pinus lambertiana).

Sugar pine, Pinus lambertiana Douglas, is a keystone species of montane forests from Baja California to southern Oregon. Like other North American white pines, populations of sugar pine have been greatly reduced by the disease white pine blister rust (WPBR) caused by a fungal pathogen, Cronartium ribicola, that was introduced into North America early in the twentieth century. Major gene resistance to WPBR segregating in natural populations has been documented in sugar pine. Indeed, the dominant resistance gene in this species, Cr1, was genetically mapped, although not precisely. Genomic single nucleotide polymorphisms (SNPs) placed in a large scaffold were reported to be associated with the allele for this major gene resistance (Cr1R). Forest restoration efforts often include sugar pine seed derived from the rare resistant individuals (typically Cr1R/Cr1r) identified through an expensive 2-year phenotypic testing program. To validate and geographically characterize the variation in this association and investigate its potential to expedite genetic improvement in forest restoration, we developed a simple PCR-based, diploid genotyping of DNA from needle tissue. By applying this to range-wide samples of susceptible and resistant (Cr1R) trees, we show that the SNPs exhibit a strong, though not complete, association with Cr1R. Paralleling earlier studies of the geographic distribution of Cr1R and the inferred demographic history of sugar pine, the resistance-associated SNPs are marginally more common in southern populations, as is the frequency of Cr1R. Although the strength of the association of the SNPs with Cr1R and thus, their predictive value, also varies with geography, the potential value of this new tool in quickly and efficiently identifying candidate WPBR-resistant seed trees is clear.

Seed morphological traits as a tool to quantify variation maintained in ex situ collections: a case study in Pinus torreyana.

Understanding the within- and among-population distribution of trait variation within seed collections may provide a means to approximate standing genetic variation and inform plant conservation. This study aimed to estimate population- and family-level seed trait variability for existing seed collections of Torrey pine (Pinus torreyana), and to use these data to guide sampling of future collections. We quantified variation in 14 seed morphological traits and seedling emergence within and among Torrey pine populations. Using a simulation-based approach, we used estimates of within-population variance to assess the number of maternal families required to capture 95 % of trait variation within each existing seed collection. Substantial structure was observed both within and among Torrey pine populations, with island and mainland seeds varying in seed size and seed coat thickness. Despite morphological differences, seedling emergence was similar across populations. Simulations revealed that 83 % and 71 % of all maternal families within island and mainland seed collections respectively needed to be resampled to capture 95 % of seed trait variation within existing collections. From a conservation perspective, our results indicate that to optimize genetic diversity captured in Torrey pine seed collections, maximizing the number of maternal families sampled within each population will be necessary.

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.

Landscape genomics of Quercus lobata reveals genes involved in local climate adaptation at multiple spatial scales.

Understanding how the environment shapes genetic variation provides critical insight about the evolution of local adaptation in natural populations. At multiple spatial scales and multiple geographic contexts within a single species, such information could address a number of fundamental questions about the scale of local adaptation and whether or not the same loci are involved at different spatial scales or geographic contexts. We used landscape genomic approaches from three local elevational transects and rangewide sampling to (a) identify genetic variation underlying local adaptation to environmental gradients in the California endemic oak, Quercus lobata; (b) examine whether putatively adaptive SNPs show signatures of selection at multiple spatial scales; and (c) map putatively adaptive variation to assess the scale and pattern of local adaptation. Of over 10 k single-nucleotide polymorphisms (SNPs) generated with genotyping-by-sequencing, we found signatures of natural selection by climate or local environment at over 600 SNPs (536 loci), some at multiple spatial scales across multiple analyses. Candidate SNPs identified with gene-environment tests (LFMM) at the rangewide scale also showed elevated associations with climate variables compared to the background at both rangewide and elevational transect scales with gradient forest analysis. Some loci overlap with those detected in other oak species, raising the question of whether the same loci might be involved in local climate adaptation in different congeneric species that inhabit different geographic contexts. Mapping landscape patterns of adaptive versus background genetic variation identified regions of marked local adaptation and suggests nonlinear association of candidate SNPs and environmental variables. Taken together, our results offer robust evidence for novel candidate genes for local climate adaptation at multiple spatial scales.

Adaptational lag to temperature in valley oak (Quercus lobata) can be mitigated by genome-informed assisted gene flow.

Climate change over the next century is predicted to cause widespread maladaptation in natural systems. This prediction, as well as many sustainable management and conservation practices, assumes that species are adapted to their current climate. However, this assumption is rarely tested. Using a large-scale common garden experiment combined with genome-wide sequencing, we found that valley oak (Quercus lobata), a foundational tree species in California ecosystems, showed a signature of adaptational lag to temperature, with fastest growth rates occurring at cooler temperatures than populations are currently experiencing. Future warming under realistic emissions scenarios was predicted to lead to further maladaptation to temperature and reduction in growth rates for valley oak. We then identified genotypes predicted to grow relatively fast under warmer temperatures and demonstrated that selecting seed sources based on their genotype has the potential to mitigate predicted negative consequences of future climate warming on growth rates in valley oak. These results illustrate that the belief of local adaptation underlying many management and conservation practices, such as using local seed sources for restoration, may not hold for some species. If contemporary adaptational lag is commonplace, we will need new approaches to help alleviate predicted negative consequences of climate warming on natural systems. We present one such approach, "genome-informed assisted gene flow," which optimally matches individuals to future climates based on genotype-phenotype-environment associations.

Seedling response to water stress in valley oak (Quercus lobata) is shaped by different gene networks across populations.

Drought is a major stress for plants, creating a strong selection pressure for traits that enable plant growth and survival in dry environments. Many drought responses are conserved species-wide responses, while others vary among populations distributed across heterogeneous environments. We tested how six populations of the widely distributed California valley oak (Quercus lobata) sampled from contrasting climates would differ in their response to soil drying relative to well-watered controls in a common environment by measuring ecophysiological traits in 93 individuals and gene expression (RNA-seq) in 42 individuals. Populations did not differ in their adjustment of turgor loss point during soil drying, suggesting a generalized species-wide response. Differential expression analysis identified 689 genes with a common response to treatment across populations and 470 genes with population-specific responses. Weighted gene co-expression network analysis (WGCNA) identified groups of genes with similar expression patterns that may be regulated together (gene modules). Several gene modules responded differently to water stress among populations, suggesting regional differences in gene network regulation. Populations from sites with a high mean annual temperature responded to the imposed water stress with significantly greater changes in gene module expression, indicating that these populations may be locally adapted to respond to drought. We propose that this variation among valley oak populations provides a mechanism for differential tolerance to the increasingly frequent and severe droughts in California.

Applying landscape genomic tools to forest management and restoration of Hawaiian koa (Acacia koa) in a changing environment.

Identifying and quantifying the importance of environmental variables in structuring population genetic variation can help inform management decisions for conservation, restoration, or reforestation purposes, in both current and future environmental conditions. Landscape genomics offers a powerful approach for understanding the environmental factors that currently associate with genetic variation, and given those associations, where populations may be most vulnerable under future environmental change. Here, we applied genotyping by sequencing to generate over 11,000 single nucleotide polymorphisms from 311 trees and then used nonlinear, multivariate environmental association methods to examine spatial genetic structure and its association with environmental variation in an ecologically and economically important tree species endemic to Hawaii, Acacia koa. Admixture and principal components analyses showed that trees from different islands are genetically distinct in general, with the exception of some genotypes that match other islands, likely as the result of recent translocations. Gradient forest and generalized dissimilarity models both revealed a strong association between genetic structure and mean annual rainfall. Utilizing a model for projected future climate on the island of Hawaii, we show that predicted changes in rainfall patterns may result in genetic offset, such that trees no longer may be genetically matched to their environment. These findings indicate that knowledge of current and future rainfall gradients can provide valuable information for the conservation of existing populations and also help refine seed transfer guidelines for reforestation or replanting of koa throughout the state.

Genetic conservation and management of the California endemic, Torrey pine (Pinus torreyana Parry): Implications of genetic rescue in a genetically depauperate species.

Rare species present a challenge under changing environmental conditions as the genetic consequences of rarity may limit species ability to adapt to environmental change. To evaluate the evolutionary potential of a rare species, we assessed variation in traits important to plant fitness using multigenerational common garden experiments. Torrey pine, Pinus torreyana Parry, is one of the rarest pines in the world, restricted to one mainland and one island population. Morphological differentiation between island and mainland populations suggests adaptation to local environments may have contributed to trait variation. The distribution of phenotypic variances within the common garden suggests distinct population-specific growth trajectories underlay genetic differences, with the island population exhibiting substantially reduced genetic variance for growth relative to the mainland population. Furthermore, F1 hybrids, representing a cross between mainland and island trees, exhibit increased height accumulation and fecundity relative to mainland and island parents. This may indicate genetic rescue via intraspecific hybridization could provide the necessary genetic variation to persist in environments modified as a result of climate change. Long-term common garden experiments, such as these, provide invaluable resources to assess the distribution of genetic variance that may inform conservation strategies to preserve evolutionary potential of rare species, including genetic rescue.

Whole-transcriptome response to water stress in a California endemic oak, Quercus lobata.

Reduced water availability during drought can create major stress for many plant species. Within a species, populations with a history of seasonal drought may have evolved the ability to tolerate drought more than those in areas of high precipitation and low seasonality. In this study, we assessed response to water stress in a California oak species, Quercus lobata Née, by measuring changes in gene expression profiles before and after a simulated drought stress treatment through water deprivation of seedlings in a greenhouse setting. Using whole-transcriptome sequencing from nine samples from three collection localities, we identified which genes are involved in response to drought stress and tested the hypothesis that seedlings sampled from climatically different regions of the species range respond to water stress differently. We observed a surprisingly massive transcriptional response to drought: 35,347 of 68,434 contigs (52%) were differentially expressed before versus after drought treatment, of which 18,111 were down-regulated and 17,236 were up-regulated. Genes functionally associated with abiotic stresses and death were enriched among the up-regulated genes, whereas metabolic and cell part-related genes were enriched among the down-regulated. We found 56 contigs that exhibited significantly different expression responses to the drought treatment among the three populations (treatment × population interaction), suggesting that those genes may be involved in local adaptation to drought stress. These genes have stress response (e.g., WRKY DNA-binding protein 51 and HSP20-like chaperones superfamily protein), metabolic (e.g., phosphoglycerate kinase and protein kinase superfamily protein), transport/transfer (e.g., cationic amino acid transporter 7 and K+ transporter) and regulatory functions (e.g., WRKY51 and Homeodomain-like transcriptional regulator). Baseline expression levels of 1310 unique contigs also differed among pairs of populations, and they were enriched for metabolic and cell part-related genes. Out of the large fraction of the transcriptome that was differentially expressed in response to our drought treatment, we identified several novel genes that are candidates for involvement in local adaptation to drought.

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.

Herbivores and edaphic factors constrain the realized niche of a native plant.

Biotic interactions, such as competition and herbivory, can limit plant species ranges to a subset of edaphically suitable habitats, termed the realized niche. Here we explored the role that herbivores play in restricting the niche of serpentine ecotypes of the native California annual Collinsia sparsiflora. We planted seeds from four populations into a range of natural field environments that varied in the presence/absence of naturally occurring C. sparsiflora and in predicted suitability for growth and survival of the serpentine ecotype of C. sparsiflora. Path analysis was then used to model the direct and herbivore-mediated indirect effects of environmental variables on the survival of C. sparsiflora serpentine ecotypes. We found that C. sparsiflora received more herbivory when planted into areas where serpentine ecotypes of C. sparsiflora were not predicted to persist, and that increased herbivory was associated with decreased survival, suggesting that herbivores may limit the distribution of C. sparsiflora serpentine ecotypes. Additionally, we demonstrated that edaphic environmental variables impacted the survival of C. sparsiflora serpentine ecotypes both directly and indirectly, by altering interactions with herbivores. These indirect effects were probably trait-mediated and probably occurred because edaphic factors may influence plant traits that, in turn, alter attractiveness to herbivores. Although the magnitude of direct effects exceeded the magnitude of indirect effects, many strong herbivore-mediated indirect effects were detected. Thus, interactions between the abiotic environment and insect herbivory contributed to restricting the niche of C. sparsiflora serpentine ecotypes to a subset of available habitat.

Collinsia sparsiflora in serpentine and nonserpentine habitats: using F2 hybrids to detect the potential role of selection in ecotypic differentiation.

Here we document phenotypic differences between serpentine and nonserpentine ecotypes of Collinsia sparsiflora, as well as patterns of selection in these contrasting soil habitats. We transplanted the two parental ecotypes and experimental F2 hybrids into six field sites, and collected morphological, phenological and fitness data on emergent plants. To focus on edaphically mediated selection, rather than on pollinator-mediated selection, we used pollinator-exclusion cages. Transplanted parentals of the two ecotypes showed genetic differentiation for floral traits, but not for phenological traits or cotyledon size. For the F2 hybrids growing on serpentine soils, there was significant directional selection on cotyledon size, flower size and flower shape. However, the pattern of selection did not differ significantly between serpentine and nonserpentine sites. Overall, we did not see evidence for divergent selection on the two soil types. We conclude that differences in floral traits between the ecotypes do not reflect adaptation to physical conditions associated with soil type, and that there are unmeasured traits that must be contributing to ecotypic differentiation.

Experimental verification of ecological niche modeling in a heterogeneous environment.

The current range of ecological habitats occupied by a species reflects a combination of the ecological tolerance of the species, dispersal limitation, and competition. Whether the current distribution of a species accurately reflects its niche has important consequences for the role of ecological niche modeling in predicting changes in species ranges as the result of biological invasions and climate change. We employed a detailed data set of species occurrence and spatial variation in biotic and abiotic attributes to model the niche of a native California annual plant, Collinsia sparsiflora. We tested the robustness of our model for both the realized and fundamental niche by planting seeds collected from four populations, representing two ecotypes, into plots that fully represented the five-dimensional niche space described by our model. The model successfully predicted which habitats allowed for C. sparsiflora persistence, but only for one of the two source ecotypes. Our results show that substantial niche divergence has occurred in our sample of four study populations, illustrating the importance of adequately sampling and describing within-species variation in niche modeling.