Landscape Biology of Western White Pine: Implications for Conservation of a Widely-distributed Five-needle Pine at Its Southern Range Limit

Throughout much of the range of western white pine, Pinus monticola Dougl., timber harvesting, fire exclusion and the presence of Cronartium ribicola J. C. Fisch., the white pine blister rust (WPBR) pathogen, have led to negative population and genetic consequences. To address these interactions, we examined population dynamics and genetic diversity in 10 populations of western white pine in upper montane forests of the Lake Tahoe Basin. We documented negative population trends for three of the 10 populations. These populations exhibited low estimated growth rates (λ), moderate to high incidences of WPBR and mountain pine beetle (MPB), and high levels of mortality. In contrast, seven populations appear to be stable (λ ě 1.0), with low to moderate disease and insect incidence, and evidence for genetic resistance to WPBR. Genetic diversity (H E) for a set of 160 single nucleotide polymorphisms was in the range of 0.245–0.272 across populations, and population-specific estimates of F ST ranged from 0.0062 to 0.0244. Allele frequency of the Cr2 gene, which confers complete resistance to C. ribicola in western white pine, was low, averaging 0.009 for all populations sampled. However, a low frequency of pollen receptors (i.e., susceptible maternal parents pollinated by a local resistant parent) was found in nine of 10 populations. A moderate and negative relationship was found between the frequency of pollen receptors in a population and the incidence of WPBR (r 2 = 0.32). In the context of an introduced pathogen, climate driven outbreaks of MPB, fire exclusion, and prolonged drought, conservation and management strategies are warranted for this species in the Lake Tahoe Basin and likely other locations in California. These strategies include gene conservation of western white pine, WPBR resistance screening, and forest restoration treatments.


Introduction
Five-needle white pine ecosystems in California, and throughout western North America, are challenged by the interacting threats of an introduced fungal pathogen, Cronartium ribicola J. C. Fisch.(the cause of white pine blister rust, WPBR), climate-driven outbreaks of mountain pine beetle (MPB) (Dendroctonus ponderosae Hopkins), climatic warming, historical logging, and a century of fire suppression [1][2][3][4][5][6][7].Unique to California is the diversity of five-needle white pine species (subgenus Strobus section Quinquefoliae) [8] with six species of the nine North American species distributed throughout the State.Many of the five-needle white pine species in California are at or near their southern range limits.In western North America, sugar (Pinus lambertiana) and western white (P.monticola Dougl.ex D. Don) pines have been of great economic, as well as ecological, importance.Much is known about the population biology, fire ecology, prevalence and resistance to white pine blister rust, and historical disturbances for sugar pine [4,6], but much less is known for western white pine in California.
Western white pine (Pinus monticola Dougl.ex D. Don) is an important co-dominant species of upper montane forests in California.Over the years this species has been severely affected by C. ribicola infection, logging, and fire suppression [1,3,[9][10][11][12][13].In California, western white pine generally occurs in high-elevation forests ranging from 1219 meters on the north coast to 3200 m in the Sierra Nevada.Historically, foresters considered the upper montane to be in the "alpine zone", ranging in elevation from 150 to 2400 m in the Cascades and North Coast ranges, to over 2700 m in the Sierra Nevada; this zone was generally deemed "too inaccessible to be logged profitably" [14].Today, the combined effects of white pine blister rust, logging, fire suppression, and climate-driven outbreaks of mountain pine beetle are affecting populations of western white pine significantly throughout its entire range [2,7,13,15].Thus, western white pine forests merit attention today because of their ecosystem values and the consequences of benign neglect.
Cronartium ribicola is a primary management concern in all five-needle white pine ecosystems in western North America.Cronartium ribicola was introduced into British Columbia around 1910 and arrived in northern California around 1929, spreading southward into the southern Sierra Nevada by 1961 [16].While all size classes of western white pine are susceptible to WPBR, individuals most often infected with stem cankers (as opposed to branch cankers) are generally less than 20 cm in diameter at breast height, making stem cankers lethal to small individuals [16]; this can dramatically alter population structure-by reducing or eliminating small-to intermediate-sized individuals within stands.In a statewide survey of high elevation white pines in California, Maloney [13] found C. ribicola infecting 42% (range: 24%-74%) of western white pine trees in the North Coast, 18% (range: 0%-72%) in the Klamath region and 14% (range: 0%-90%) in the Sierra Nevada.Although WPBR can be fatal to all five-needle white pine species, a simply inherited gene in western white pine confers resistance to C. ribicola [17].This gene, designated Cr2, controls a hypersensitive response in needles that prevents fungal growth [17] and is analogous to Cr1, which confers resistance to C. ribicola in sugar pine [18].Cr1 and Cr2 are the foundation for naturally occurring resistance to C. ribicola in sugar pine and western white pine in California.
Fire is an important disturbance agent in mixed conifer forests throughout California, and while fire return intervals in upper montane forests are likely less frequent than in lower-elevation forests, the effects of fire suppression (e.g., stand densification) are evident in these western white pine forest communities [3,[19][20][21][22][23].Studies from upper montane/red fir zone (Abies magnifica Andr.Murray) and high-elevation forests in California report fire return intervals between 15-76 years (range: five-175 years) [3,[19][20][21][22]24].Given fire suppression policies and estimates of fire return intervals for upper montane and comparable forest types, we assume that many upper montane forests throughout California have missed one-three intervals of fire.
These disturbances (WPBR, MPB, fire exclusion, and logging) and their interactions can have negative population-and genetic-level consequences throughout the range of western white pine in California.At the 1300 km 2 scale of the Lake Tahoe Basin, we have taken a landscape approach to understanding the population dynamics, genetic structure, and diversity of western white pine.Our objectives were to determine: (i) the current population structure and trends of western white pine in the Lake Tahoe Basin; (ii) the genetic structure and diversity of extant western white pine Forests 2016, 7, 93 3 of 17 populations; and (iii) the landscape-scale frequency of WPBR infection and resistance to WPBR (Cr2) in western white pine.We present our results as a summary of important ecological and genetic factors that should be considered to provide guidance for conservation and management strategies appropriate to a forest tree species at its southern range limit.

Study Area
The Lake Tahoe Basin is located in the north-central Sierra Nevada at 120 ˝W and 39 ˝N within the States of California and Nevada.The Basin is flanked to the west by the Sierra Nevada crest and to the east by the Carson Range, and elevations range from 1890 to 3292 m.The climate is Mediterranean and characterized by warm, dry summers and cool wet winters.Most precipitation falls as snow between the months of November and April, with a strong gradient from west to east.East-side locations are influenced by the rainshadow effect of the Sierra Nevada crest (Table 1).Geology of the region is dominated by igneous intrusive rocks, typically granodiorite, and igneous extrusive rocks, typically andesitic lahar, with small amounts of metamorphic rock [25].
In California and western Nevada, western white pine grows in mixed-species stands and is an important associate of upper montane forests at elevations ranging from 2100 m to 2900 m, but can sometimes be found up to 3300 m in the subalpine zone [26].In the Lake Tahoe Basin several forest types occur in the upper montane, but the most dominant type is the red fir forest, with a mix of red fir (A.magnifica Andr.Murray), western white pine, lodgepole pine (P.contorta Dougl.ex.Loud.) and occasionally Jeffrey pine (P.jeffreyi Grev.& Balf.), covering approximately 5000 ha [27].
During the summer of 2008, we selected 10 study populations, with 3 permanent demographic plots per population (sampling area within a population = 4 ha), for a total of 30 plots.For our study, a population is defined as an area in which gene flow occurs (pollen and seed dispersal), which for some conifers has been estimated at about 4 hectares [28,29].Our long-term plot network was established on National Forest System and State Park lands surrounding Lake Tahoe (Lake Tahoe Basin Management Unit/LTBMU and California and Nevada State Parks).Each of the 10 populations of western white pine was located in a distinct watershed and distributed around the Basin to capture variation in the physical environment (e.g., climate, geology, topography).For each demographic plot, annual precipitation from the period 1971-2000 was provided by FHTET (USDA FS Forest Health Technology Enterprise Team, Fort Collins, CO, USA) using the PRISM climatic model [30].Soil survey data were provided by the South Lake Tahoe office of the USDA Natural Resources Conservation Service (NRCS).For a complete description of field methodology as well as insect and pathogen identification, see Maloney et al. [6].
Table 1.Biological and environmental summaries for 10 populations of Pinus monticola in the Lake Tahoe Basin.Stands are listed clockwise in orientation in the Lake Tahoe Basin from north, east, south, and west.d.b.h.(cm) (diameter at breast height); % WPBR (white pine blister rust) = % incidence of WPBR; % MPB (mountain pine beetle) = % incidence of MPB; Elev.= elevation above sea level/a.s.l.(m); Ann.ppt.= total annual precipitation in millimeters (mm); T min ( ˝C) = minimum January temperature; T max ( ˝C) = maximum July temperature; WC 15 bar = soil water capacity at 15 bar; % sand = percentage sand content (sand:silt:clay); CEC = cation exchange capacity.See Methods for sources of climate and soil survey data.

Current Population Structure and Trends
Population trends were assessed with transition matrix models for each of the 10 populations of western white pine.We used transition matrices to describe and summarize current trends in survivorship, fecundity and, to a lesser extent, radial growth.Transition matrix models of populations follow the Lefkovitch [31] model.Entries in the transition matrix represent the contribution that each size class makes to every other size class during a specified time interval.We created sized-based models from six size classes: For growth probabilities, individuals transitioned into the next size class when they were in the upper d.b.h.limit of the size class, or of height for seedlings/saplings.For example, we assumed that an individual with a d.b.h. of 4.9-5.0cm would grow into size class 3 in the next time step.Where mortality occurred, survival was calculated using the number of standing dead trees (years dead ranged from 1 to 12 years) divided by current live and dead stems for each size class.In demographic studies of forest trees where mortality was not observed, mortality was assumed to be 0% or 2% [4,32].In our populations where no mortality was observed, we assumed a minimal value of 1% mortality, an average of these previously published estimates [4,32].Fecundity was estimated from existing recruitment and cone production data for size classes 3, 4, 5, and 6 [6].The population growth rate (λ), was estimated using the dominant eigenvalue of the transition matrix [33].We computed 95% confidence intervals for λ by bootstrapping (n = 10,000) across survivorship, growth and fecundity estimates comprising the transition matrix.All calculations and bootstrapping were performed in Matlab [34].Stepwise multiple linear regression was used to identify important biological and environmental variables explaining western white pine mortality in the Lake Tahoe Basin [35].The stepwise criterion was run with forward direction with p ď 0.10 for variable entry.A check of collinearity in regression models was performed using leverage plots and bivariate scatterplots.Data for water use efficiency, a variable used in the multiple linear regression model, is reported in Maloney et al., (In revision).

Genetic Structure and Diversity
Within demographic plots, needle samples were taken from 25 adult western white pine individuals from each of the 10 populations in the Lake Tahoe Basin (n = 250).Sampled trees were approximately 30 m apart.Needle samples were stored in vials containing three 0.5 g desiccant packets.DNA was extracted in the laboratory from 50 mg of tissue using DNeasy ® Plant 96 Kits following the manufacturer's protocol (Qiagen, Valencia, CA, USA).DNA from sampled trees was genotyped for 384 single nucleotide polymorphisms (SNPs), which were used as genetic markers.These SNPs were discovered through resequencing of a diversity panel comprising 12 haploid megagametophytes collected across the natural range of sugar pine [36].All genotyping was carried out at the DNA Technologies Core at the UC Davis Genome Center.Genotypes were called at each SNP locus for each tree using GenomeStudio (Illumina) software.Quality thresholds for inclusion of SNPs were the same as those given in Eckert et al. [37].
For each SNP, we calculated observed (H O ) and expected (H E ) heterozygosity, Wright's inbreeding coefficient (F IS ) and hierarchical F-statistics among plots nested within populations.Significance of multilocus hierarchical fixation indices was assessed using 99% bootstrap confidence intervals (CIs) (n = 10,000 replicates sampled across SNPs).We examined allelic correlations (r 2 ) between all possible pairs of SNPs to validate our assumption that SNPs were unlinked.Analyses were conducted using the GENETICS and HIERFSTAT [38] packages in R [39].
We assessed population structure using the approach outlined by Nicholson et al. [40] as implemented in the POPGEN library of R [39].This approach models a set of populations that diverges instantaneously from an ancestral population.Allele frequencies in those diverging populations drift away from those in the ancestral population, and are modeled using a binomial-beta hierarchical structure [41].Here, we used the 10 sampled populations of western white pine.There are three parameters of interest in this model: α ilj , the frequency in population i at locus l of allele j; π lj , the ancestral frequency of locus l of allele j; and c i , the variance parameter for population i describing the magnitude of drift away from ancestral allele frequencies across loci.Since the c i parameters are defined in terms of variance in allele frequencies across loci, they are analogous to population-specific estimates of multilocus F ST .Thus, a larger value for an estimate of the c i parameter indicates more drift away from ancestral allele frequencies, meaning that current allele frequencies within population i (α ilj ) are less similar, in terms of variance across loci, to those in the ancestral population (π lj ).As with F ST , there are no natural breaks to define easily what is high and what is low without assuming some form of equilibrium model of gene flow (e.g., island model).A comprehensive review of F ST estimates for conifers is given elsewhere [42,43].
Parameters were estimated in a Bayesian framework using Markov chain Monte Carlo (MCMC), where uniform priors were placed on the c i and π lj parameters and a Metropolis-Hastings algorithm was used to sample from the posterior distribution of each parameter.The α ilj parameters were integrated out of the likelihood function and thus were not sampled (cf.[40]).Markov chains were sampled 10,000 times post burn-in (n = 10,000) to estimate posterior distributions, with point estimates reported as posterior means.Convergence was assessed by comparing posterior distributions of the c i parameters across two replicated runs of the MCMC sampler, using Kolmogorov-Smirnov tests.

Evaluating Disease Resistance and Cr2 Allele Frequency
One-hundred and thirty one families (none of which exhibited symptoms of WPBR infection in field evaluations) were screened for presence of the Cr2 allele at the Placerville Nursery's Rust Resistance Screening Program facility, Camino, CA.Western white pine seed was sown in 2009 (see Maloney et al. [6] for a complete description of culturing and screening protocols).In October 2010 (96 families) and October 2011 (35 families), test seedlings were inoculated by basidiospore casts from infected leaves of Ribes nigrum L., an alternate host to the WPBR pathogen.Ribes nigrum plants had been inoculated with uredospores derived from stored aeciospores of Cronartium ribicola collected and bulked from forested sites on the Eldorado National Forest by the Placerville Nursery, Camino, CA, USA.Seedlings with infected Ribes leaves suspended over them were incubated in dew chambers for 62 h at 17 ˝C and 100% relative humidity (J.Gleason, pers.comm., USDA Forest Service).In March of 2011 and 2012, the frequency of Cr2 in the western white pine family was evaluated.When challenged by C. ribicola, resistant host needles respond to infection by developing hypersensitive spots (small necrotic flecks) that prevent further fungal growth; these flecks contrast markedly with the bright yellow, expanding spots that develop on the foliage of susceptible Cr2/Cr2 genotypes [17].Since the Cr2 allele is dominant and segregates in a 1:1 ratio of resistant to susceptible reactions, its frequency can be estimated for wind-pollinated families by counting the number of seedlings expressing the hypersensitive reaction [17].Parents of families that segregate 1:1 resistant:susceptible are heterozygous for the Cr2 allele (genotype, Cr2Cr2).Cr2 allele frequency in a population is calculated by dividing the number of Cr2 alleles identified per family by twice the total number of families tested.Correlation analyses were conducted to determine relationships between WPBR incidence and the frequency of Cr2 and of pollen receptors (susceptible families that exhibited a very low frequency of resistant offspring, presumably derived from rare Cr2 alleles in the pollen cloud) [35].

Western White Pine and Forest Conditions
Western white pine densities within plots ranged from a low of 56 trees¨ha ´1 to a high of 111 trees¨ha ´1 (mean = 82.5,Table 1).Basal area also varied, ranging from 5.27 m 2 ¨ha ´1 to 21.9 m 2 ¨ha ´1 (mean = 13.7,Table 1).White pine blister rust was found in 90% of the populations and the frequency of infected trees by population ranged from 0% to 45% (Table 1).The highest frequencies of WPBR infection were found at Blackwood Canyon (45%) and Mt Watson (23%).No rust was found Forests 2016, 7, 93 7 of 17 on western white pine in the Heavenly population.Mountain pine beetle was found in 100% of the populations and ranged from 6% to 40% of trees infested (Table 1).High incidence of MPB was found in populations at Heavenly (40%), Flume Trail (23%), Incline Lake (17%), Armstrong Pass (17%), and Blackwood Canyon (17%) (Table 1).Western white pine mortality ranged from 2% to 17%, with the highest cumulative mortality at Heavenly and Blackwood Canyon (Table 1).
Strong environmental gradients exist within the Basin in precipitation, geology and soil type (Table 1).Mean annual precipitation is 1265 mm in west-side and 963 mm in east-side locations, respectively (Table 1).WPBR incidence is highest at Blackwood Canyon on the west-side of the Basin, which was also the most mesic location, receiving, on average, 1472 mm of precipitation annually (Table 1, Figure 1).Moist conditions are required by C. ribicola for successful infection on both hosts; all spore stages require relative humidity >90% and temperatures in the range of 15-18 ˝C [44][45][46].Soil type can influence site characteristics and carrying capacity of western white pine and other forest associates (Table 1).Populations growing on decomposed granite have relatively lower densities (mean = 75 inds.¨ha´1) than neighboring populations growing on volcanic or andesitic soils, which generally have higher available soil moisture (WC 15 bar mean = 7.25 volcanic/andesite versus 3.8 granitic ), and thus support greater tree densities (mean = 94 inds.¨ha´1) (Table 1).
Strong environmental gradients exist within the Basin in precipitation, geology and soil type (Table 1).Mean annual precipitation is 1265 mm in west-side and 963 mm in east-side locations, respectively (Table 1).WPBR incidence is highest at Blackwood Canyon on the west-side of the Basin, which was also the most mesic location, receiving, on average, 1472 mm of precipitation annually (Table 1, Figure 1).Moist conditions are required by C. ribicola for successful infection on both hosts; all spore stages require relative humidity >90% and temperatures in the range of 15-18 °C [44][45][46].Soil type can influence site characteristics and carrying capacity of western white pine and other forest associates (Table 1).Populations growing on decomposed granite have relatively lower densities (mean = 75 inds.•ha−1 ) than neighboring populations growing on volcanic or andesitic soils, which generally have higher available soil moisture (WC 15 bar mean = 7.25volcanic/andesite versus 3.8granitic), and thus support greater tree densities (mean = 94 inds.•ha−1 ) (Table 1).

Population Structure and Trends
Multi-size structure was found for western white pine, with good representation in most size classes, particularly juveniles in the 0.1-5.0cm diameter size class (Figure 2).High proportions of large diameter trees (>40.1 cm) were found at Incline Lake, Flume Trail, Heavenly, and Jakes (Figure 2).Relatively high numbers of seedlings and saplings were found in all populations with the greatest recruitment numbers at Montreal Canyon and Mt Watson, which also corresponds with logging activity within the last 20 years (Figure 2).

Population Structure and Trends
Multi-size structure was found for western white pine, with good representation in most size classes, particularly juveniles in the 0.1-5.0cm diameter size class (Figure 2).High proportions of large diameter trees (>40.1 cm) were found at Incline Lake, Flume Trail, Heavenly, and Jakes (Figure 2).Relatively high numbers of seedlings and saplings were found in all populations with the greatest recruitment numbers at Montreal Canyon and Mt Watson, which also corresponds with logging activity within the last 20 years (Figure 2).Variation in fecundity, which is a function of cone production and the number of recruits that establish successfully, was found among populations (Table 2).Fecundity was greatest at Mt Watson (1.437) and Montreal Canyon (1.036), the two locations that have been most recently logged.Low fecundity was found at Jakes Peak (0.076), Echo Lake (0.112), and Flume Trail (0.150) (Table 2).Survivorship ranged from 0.828 to 0.983 with a mean of 0.918 (Table 2).Estimated growth rates (λ) for 7 of the 10 western white pine populations were ≥1.0, indicating that these populations are apparently stable (Table 2).Blackwood Canyon (0.946), Incline Lake (0.991), and Heavenly (0.997) had λ's <1.0, indicative of populations that may be in decline (Table 2).Because population growth, λ, was calculated from a one-time sampling, our estimates may not reflect intrinsic variation in rates of survivorship, mortality, growth, and fecundity.Our intent was not to predict future population growth but to describe current population conditions.In addition, the confidence intervals are very  Variation in fecundity, which is a function of cone production and the number of recruits that establish successfully, was found among populations (Table 2).Fecundity was greatest at Mt Watson (1.437) and Montreal Canyon (1.036), the two locations that have been most recently logged.Low fecundity was found at Jakes Peak (0.076), Echo Lake (0.112), and Flume Trail (0.150) (Table 2).Survivorship ranged from 0.828 to 0.983 with a mean of 0.918 (Table 2).Estimated growth rates (λ) for 7 of the 10 western white pine populations were ě1.0, indicating that these populations are apparently stable (Table 2).Blackwood Canyon (0.946), Incline Lake (0.991), and Heavenly (0.997) had λ's <1.0, indicative of populations that may be in decline (Table 2).Because population growth, λ, was calculated from a one-time sampling, our estimates may not reflect intrinsic variation in rates of survivorship, mortality, growth, and fecundity.Our intent was not to predict future population growth but to describe current population conditions.In addition, the confidence intervals are very large for λ, which is a cause for concern, but it is difficult to assess the magnitude of this concern, Forests 2016, 7, 93 9 of 17 as λ was estimated from a one-time sampling that likely resulted in the large variances around the point estimate.
Table 2. Demographic and genetic summaries for P. monticola populations from the Lake Tahoe Basin.Mean fecundity, survival (cumulative mortality in bold and parentheses below), and population growth rate (λ along with 2.5% and 97.5% confidence intervals) from transition matrices.Genetic and diversity patterns from 160 marker loci with mean and variance (in parentheses) for H O , H E , F IS , c i , with 95% credible intervals for c i (in parentheses), Cr2 allele frequency, and mean water-use efficiency, δ 13 C. Multiple linear regression of western white pine mortality employed 3 biotic and abiotic variables explaining 88% of the variation (model R 2 = 0.88, F 3,6 = 15.19,p = 0.0029).Mountain pine beetle contributed significant predictive power to the regression model, followed by CEC (cation exchange capacity, a soil property indicative of soil fertility), and water-use efficiency, δ 13 C (Figure 3).was estimated from a one-time sampling that likely resulted in the large variances around the point estimate.
Table 2. Demographic and genetic summaries for P. monticola populations from the Lake Tahoe Basin.Mean fecundity, survival (cumulative mortality in bold and parentheses below), and population growth rate (λ along with 2.5% and 97.5% confidence intervals) from transition matrices.Genetic and diversity patterns from 160 marker loci with mean and variance (in parentheses) for HO, HE, FIS, ci, with 95% credible intervals for ci (in parentheses), Cr2 allele frequency, and mean water-use efficiency, δ 13 C. Multiple linear regression of western white pine mortality employed 3 biotic and abiotic variables explaining 88% of the variation (model R 2 = 0.88, F3,6 = 15.19,p = 0.0029).Mountain pine beetle contributed significant predictive power to the regression model, followed by CEC (cation exchange capacity, a soil property indicative of soil fertility), and water-use efficiency, δ 13 C (Figure 3).

Genetic Structure and Diversity
Of the 384 SNP marker loci assessed, 160 were of high quality and segregated in our sample of 250 trees.Genetic diversity measures (HO and HE) ranged from 0.231 to 0.259, and from 0.245 to 0.272, respectively (Table 2 and Figure 4).Expected heterozygosity was similar to values reported by Kim et al. [47], derived from a limited number of samples of western white pine in California (mean HE = 0.269).The inbreeding coefficient (FIS) ranged from 0.025 to 0.146, with all populations having positive values (Table 2 and Figure 4).High positive values of FIS often indicate inbreeding, with fewer heterozygotes than expected, or within population substructure (e.g., Wahlund effect).Three populations had FIS values >0.100: Blackwood Canyon (0.146), Jakes Peak (0.112), and Heavenly (0.109).The remaining populations were within range of zero and are consistent with the Hardy-Weinberg equilibrium (Table 2).
The drift parameter (ci) ranged from 0.0062 to 0.0244 (Table 2 and Figure 4), consistent with the multilocus estimate of FST across all populations (FST = 0.0124, 99% CI: 0.0077-0.0185).Relatively low values of FST could indicate high rates of gene flow via pollen dispersal.The highest value of ci was at Incline Lake (0.0244), followed by Meiss Meadow (0.0216) (Table 2 and Figure 4).A moderate and negative relationship exists between western white pine fecundity and ci (r 2 = 0.21, F1,9 = 2.19, p = 0.17).Maloney et al. [6] found a similar pattern of high genetic drift for two sugar pine populations in the Lake Tahoe Basin-likely a consequence of historical logging in the region, in both lower and upper

Genetic Structure and Diversity
Of the 384 SNP marker loci assessed, 160 were of high quality and segregated in our sample of 250 trees.Genetic diversity measures (H O and H E ) ranged from 0.231 to 0.259, and from 0.245 to 0.272, respectively (Table 2 and Figure 4).Expected heterozygosity was similar to values reported by Kim et al. [47], derived from a limited number of samples of western white pine in California (mean H E = 0.269).The inbreeding coefficient (F IS ) ranged from 0.025 to 0.146, with all populations having positive values (Table 2 and Figure 4).High positive values of F IS often indicate inbreeding, with fewer heterozygotes than expected, or within population substructure (e.g., Wahlund effect).Three populations had F IS values >0.100: Blackwood Canyon (0.146), Jakes Peak (0.112), and Heavenly (0.109).The remaining populations were within range of zero and are consistent with the Hardy-Weinberg equilibrium (Table 2).
The drift parameter (c i ) ranged from 0.0062 to 0.0244 (Table 2 and Figure 4), consistent with the multilocus estimate of F ST across all populations (F ST = 0.0124, 99% CI: 0.0077-0.0185).Relatively low values of F ST could indicate high rates of gene flow via pollen dispersal.The highest value of c i was at Incline Lake (0.0244), followed by Meiss Meadow (0.0216) (Table 2 and Figure 4).A moderate and negative relationship exists between western white pine fecundity and c i (r 2 = 0.21, F 1,9 = 2.19, p = 0.17).Maloney et al. [6] found a similar pattern of high genetic drift for two sugar pine populations in the Lake Tahoe Basin-likely a consequence of historical logging in the region, in both lower and upper.

Frequency of Disease Resistance
Of the 131 families (maternal parents) screened for Cr2, 4 were heterozygous for the resistance gene (Cr2Cr2), yielding an overall Cr2 allele frequency of 0.023 for sampled western white pine in the Lake Tahoe Basin.Kinloch et al. [17] reported Cr2 frequencies in the Sierra Nevada to be between 0.0046 and 0.0081 (this is because he was reporting gametic frequencies, not zygotic frequencies), nearly 3-fold lower than what we found in the Lake Tahoe Basin.Cr2 was found in 2 of 10 western white pine populations (Table 2), with zygotic frequencies of 0.07 at Meiss Meadow and 0.02 at Armstrong Pass (Table 2, Figure 1).No evidence for Cr2 was found in the remaining populations of western white pine; however, pollen receptors were identified in 9 of the 10 western white pine populations.Pollen receptors (Cr2Cr2) are susceptible female parents that conferred resistance to some of their progeny from rare R alleles in the local pollen cloud.We found a moderate and negative relationship between the frequency of pollen receptors in a population and the incidence of WPBR, r 2 = 0.32, F1,9 = 3.85, p = 0.08 (Figure 5).

Discussion
We evaluated population and genetic characteristics of western white pine in the Lake Tahoe Basin as they relate to stand conditions, land-use, disease, and insect prevalence.Current forest

Frequency of Disease Resistance
Of the 131 families (maternal parents) screened for Cr2, 4 were heterozygous for the resistance gene (Cr2Cr2), yielding an overall Cr2 allele frequency of 0.023 for sampled western white pine in the Lake Tahoe Basin.Kinloch et al. [17] reported Cr2 frequencies in the Sierra Nevada to be between 0.0046 and 0.0081 (this is because he was reporting gametic frequencies, not zygotic frequencies), nearly 3-fold lower than what we found in the Lake Tahoe Basin.Cr2 was found in 2 of 10 western white pine populations (Table 2), with zygotic frequencies of 0.07 at Meiss Meadow and 0.02 at Armstrong Pass (Table 2, Figure 1).No evidence for Cr2 was found in the remaining populations of western white pine; however, pollen receptors were identified in 9 of the 10 western white pine populations.Pollen receptors (Cr2Cr2) are susceptible female parents that conferred resistance to some of their progeny from rare R alleles in the local pollen cloud.We found a moderate and negative relationship between the frequency of pollen receptors in a population and the incidence of WPBR, r 2 = 0.32, F 1,9 = 3.85, p = 0.08 (Figure 5).

Frequency of Disease Resistance
Of the 131 families (maternal parents) screened for Cr2, 4 were heterozygous for the resistance gene (Cr2Cr2), yielding an overall Cr2 allele frequency of 0.023 for sampled western white pine in the Lake Tahoe Basin.Kinloch et al. [17] reported Cr2 frequencies in the Sierra Nevada to be between 0.0046 and 0.0081 (this is because he was reporting gametic frequencies, not zygotic frequencies), nearly 3-fold lower than what we found in the Lake Tahoe Basin.Cr2 was found in 2 of 10 western white pine populations (Table 2), with zygotic frequencies of 0.07 at Meiss Meadow and 0.02 at Armstrong Pass (Table 2, Figure 1).No evidence for Cr2 was found in the remaining populations of western white pine; however, pollen receptors were identified in 9 of the 10 western white pine populations.Pollen receptors (Cr2Cr2) are susceptible female parents that conferred resistance to some of their progeny from rare R alleles in the local pollen cloud.We found a moderate and negative relationship between the frequency of pollen receptors in a population and the incidence of WPBR, r 2 = 0.32, F1,9 = 3.85, p = 0.08 (Figure 5).

Discussion
We evaluated population and genetic characteristics of western white pine in the Lake Tahoe Basin as they relate to stand conditions, land-use, disease, and insect prevalence.Current forest

Discussion
We evaluated population and genetic characteristics of western white pine in the Lake Tahoe Basin as they relate to stand conditions, land-use, disease, and insect prevalence.Current forest conditions (45%) (see Figure 4).Cr2 frequencies in the Lake Tahoe Basin were nearly three-fold higher than what Kinloch et al. [17] reported.

Conclusions
Upper montane forests provide significant biological and hydrological functions and services, including biodiversity, high carbon storage, and deep and long-lasting snowpacks [55][56][57].In the context of fire suppression, an introduced pathogen, prolonged droughts, and catastrophic wildfires, conservation and management strategies appear to be warranted for this important species and forest type.As in lower elevation forests, restoration treatments that mitigate fire and drought may be appropriate for some locations in the upper montane [58,59].Similar to Taylor [3], recommendations to reintroduce fire appear necessary if forest conditions allow, to promote natural western white pine regeneration.Without management action, catastrophic wildfires, infection by WPBR, and elevated MPB activity in this forest type may cause significant ecological (e.g., reduced population sizes), environmental (e.g., greater fuels, carbon loss), and genetic losses, including loss of naturally occurring resistance to WPBR, already present in western white pine.Rangewide seed collections and screening for Cr2 are needed for western white pine in California.Such efforts will be needed to enhance collections of rust-resistant western white pine as well as provide a local and diverse seed source for post-wildfire restoration and post-MPB outbreak recovery.

Figure 1 .Figure 1 .
Figure 1.Location of study area and western white pine plots in the Lake Tahoe Basin.Pinus monticola distribution in California is shown in insert (distribution map source USGS 1999 Digital

Figure 2 .
Figure 2. Seedling and sapling abundance and size structure for ten population of western white pine in the Lake Basin.Population structure is represented by stems in d.b.h.size classes 0.1-5.0cm, 5.1-10.0cm, 10.1-20.0cm, 20.1-4.0 cm, and ≥40.1 cm.Arrows denote populations that were recently logged.

Figure 2 .
Figure 2. Seedling and sapling abundance and size structure for ten population of western white pine in the Lake Basin.Population structure is represented by stems in d.b.h.size classes 0.1-5.0cm, 5.1-10.0cm, 10.1-20.0cm, 20.1-4.0 cm, and ě40.1 cm.Arrows denote populations that were recently logged.

Figure 3 .
Figure 3. Multiple linear regression of western white pine mortality employed three biotic and abiotic variables explaining 88% of the variation (model R 2 = 0.88, F 3,6 = 15.19,p = 0.0029).Mountain pine beetle contributed a significant amount of predictive power to the regression model (A), followed by CEC (B), and water-use efficiency, δ 13 C (C).

Figure 4 .
Figure 4. Genetic and trait patterns for western white pine populations from the Lake Tahoe Basin.

Figure 5 .
Figure 5. Moderate and negative relationship between the incidence of WPBR and frequency of pollen receptors in a population.Blackwood Canyon is denoted by the arrow.

Figure 4 .
Figure 4. Genetic and trait patterns for western white pine populations from the Lake Tahoe Basin.

Figure 4 .
Figure 4. Genetic and trait patterns for western white pine populations from the Lake Tahoe Basin.

Figure 5 .
Figure 5. Moderate and negative relationship between the incidence of WPBR and frequency of pollen receptors in a population.Blackwood Canyon is denoted by the arrow.

Figure 5 .
Figure 5. Moderate and negative relationship between the incidence of WPBR and frequency of pollen receptors in a population.Blackwood Canyon is denoted by the arrow.
Atlas of United States Trees'' by Elbert L. Little, Jr.).Populations with Cr2 are noted and asterisks indicate populations where pollen receptors are present.