Trevor Carter and Alison Agneray advised by Dr.Elizabeth Leger
University of Nevada, Reno.
Citation: T. Carter and A. Agneray, “Variability of whitebark pine (Pinus albicaulis Engelm.) leaf traits in the Great Basin” Nevada State Undergraduate Research Journal. V6:I1 Spring-2020. (2020). http://dx.doi.org/10.15629/126.96.36.199.5_6-1_S-2020_1
Whitebark pine are valued for the ecosystem services they provide in subalpine forests of the western United States and have been declining across their range. This project quantiﬁes two leaf traits of contemporary and historical populations (sites) of whitebark pine within the Great Basin. Characteristics of historical herbarium specimens were compared against samples collected in the 2018 ﬁeld season of four diﬀerent populations of whitebark pine within the Great Basin and the eastern Sierra Nevada. We asked how these populations diﬀered from each other. Little change was observed through time for any of the sites, but leaf trait values were diﬀerent among populations. The Jarbidge site within the Great Basin showed the most diﬀerent leaf trait values, with smaller leaf mass per area and fewer leaves than other sites, in both historical and contemporary samples. Our research suggests there are diﬀerences among populations that may reﬂect important diﬀerences among growing conditions, genetic variation, or a combination of these factors. Additional research is needed to determine what is driving variation among whitebark pine within the Great Basin.
Whitebark pine is valued for its ecosystem services and the wildlife habitat it provides (1). It has an expansive range in high elevation western North American forests (2), but whitebark pine has been declining across its range due to a variety of interconnected mortality agents including white pine blister rust (3), the mountain pine beetle (4, 5), and drought (6). The inﬂuence of drought is particularly crucial because it can cause direct mortality (6), lead to decreased tree vigor, and exacerbate other sources of mortality (7, 8). The eﬀects of drought are especially felt in environments where resources are limited, such as sub-alpine ecosystems (9). However, the inﬂuence of drought on whitebark pine adaptation has not been fully examined.
The eﬀects of drought on plants are measured by changes in their ecophysiology (10, 11). These processes can be inferred through leaf traits such as leaf quantity (12), and leaf mass per area (13). Previous studies have shown that leaf traits vary through time in both annual (e.g., Ref. (14)) and perennial species (e.g., Ref. (15)), and leaf morphology is known to respond to drought (12). Understanding the variability of a species across environmental gradients is a key ﬁrst step to understanding whether populations are adapted to their local environments (16) and measuring the variation in traits over time can indicate whether morphology is changing in response to factors such as drought or harvest (17). Leaf traits are easy to measure on herbarium specimens, allowing for analysis through time, even of long-lived perennials.
Determining the distribution of trait variation across populations can better guide restoration eﬀorts (18). Whitebark pine populations are locally adapted across geographic areas with traits related to cold adaptation, and current restoration practices for whitebark pine utilize seed transfer guidelines based on climatic factors (19). These practices include moving genetic stock across environments that do not diﬀer by more than 1.0◦C in mean temperature (19). Whitebark pine habitats also diﬀer in a variety of other environmental conditions such as vapor pressure deﬁcit and precipitation (20), which haven’t been studied for this species. Work to-date on trait distribution in whitebark pine has not fully incorporated populations within the Great Basin, and additional drought-speciﬁc traits may be of interest. For example, leaf number is known to decrease as a way for trees to manage water stress (21), while leaf mass per area (LMA) is shown to be higher in more drought-resistant plants (22). The relationship between LMA and drought may be more complicated. For example, LMA was documented to decrease in Scots pine (Pinus sylvestris) during drought-induced mortality events (23).
The objectives of this research were to ask the following questions: (1) Do leaf quantity and LMA diﬀer within contemporary whitebark pine populations in the Great Basin and the Sierra Nevada mountains? And (2) have leaf quantity, and LMA changed through time? We predicted that populations would diﬀer and that regions with greater water limitation would have lower leaf quantity and higher LMA. We expected that leaf quantity would decrease over time in areas experiencing climate warming, and would increase in areas getting wetter over time, with LMA showing the opposite pattern. We were also interested in whether populations diﬀered in their ability to change in these leaf traits over time, which might indicate a higher degree of plasticity within certain populations and better predict success after seed transfer (16).
Material and methods
We studied four populations of whitebark pine: Mount Rose (MR), Eastern Sierra (ES), Ruby Mountains (RM), and Jarbidge (Jar). We chose populations based on the availability of historical specimens at the University of Nevada, Reno herbarium (Appendix Table 2). Samples were identiﬁed using the Intermountain Region Herbarium Network (24) and selected based on specimen location. Each population, deﬁned as a diﬀerent mountain range, contained a minimum of three historical samples collected at least ten years before sampling. The date of collection for specimens ranged from 1912 – 2006. While these historic specimens varied greatly in age, we included all of them due to the limited data available for historic measurements. Geographic locations were determined based on the collection information provided in the descriptions of the specimens. These locations created a starting point for identifying the contemporary population in the same locality.
To determine the target tree branch size for ﬁeld sampling, we measured the stem diameter of herbarium specimens at the nadir of each sample, near the point where the sample was removed from the tree. We then averaged the diameters within populations to create a target diameter for contemporary sampling, which were collected at each of the four mountain ranges based on this measure. Averaging the stem diameter was done to reduce sampling biases based on branch age, and to conduct the most accurate comparison of modern samples to historical specimens. Trees (n = 30 per site) were chosen randomly within a 5-km radius. Two samples were collected per tree based on the stem diameter at varying locations of the tree. Global Positioning System (GPS) points were recorded at every replicate tree (Appendix Table 1). Samples were then stored in a plant press for at least one week to preserve leaf tissue.
Leaf quantity and LMA were estimated based on methods used by Abrams et al. (10). We counted every fascicle for contemporary specimens and multiplied the ﬁnal count by ﬁve, assuming every fascicle had ﬁve needles. Individual needles were counted for herbarium specimens due to fascicles of herbarium specimens not being readily visible without destructive sampling. The needle count of herbarium specimens was rounded up to the next multiple of ﬁve, assuming that needle count would be a multiple of ﬁve. LMA measurements were conducted on both contemporary and historical specimens. Needles were randomly pulled from a bag of contemporary samples to ensure random selection. The herbarium director permitted us to remove a single needle from the specimen mounts; this was done for all but one delicate specimen (Appendix Table 2). If able, needles were selected at random from a small pocket at the bottom of the herbarium sheet where residual material was stored. All leaves were oven-dried for two days at 40◦C before measuring LMA, similar to work by Hultine & Marshall (25).
Precipitation and vapor pressure deﬁcit (VPd) data were gathered from the PRISM Climate Group (20) for each location and used as a proxy for drought stress. The trends of 30-year average precipitation were used to assess change over time for each population.
All of our analyses were performed using R statistical software (26). ANOVA was used to compare contemporary populations based on LMA and leaf quantity metrics. Ttests were used to compare LMA and leaf quantity through time within the same population. For the purpose of a pairwise statistical test, age was categorized as either “contemporary” or “historical.” These data were compared against the leaf trait values for the corresponding population.
Do leaf quantity and LMA diﬀer within four contemporary populations of the Great Basin and eastern Sierra? Precipitation was variable across the four populations (Table 2). The Mount Rose population had the highest 30-year average precipitation, and Jarbidge had the lowest average rainfall. However, the VPd was less variable and showed no apparent diﬀerences among populations (Table 2). Precipitation was increasing at three sites, and decreasing at one.
Contemporary populations diﬀered in leaf traits, with the Jarbidge population varying the most other populations (Figure 1). Average leaf quantity per replicate at Jarbidge (n = 613) was signiﬁcantly lower than other populations (Jar-ES p = 0.0110; Jar-MR p < 0.0001; Jar-RM p = 0.0001). The other populations did not diﬀer signiﬁcantly from each other in average leaf quantity (Figure 2). Similarly, the LMA measurements for Jarbidge (0.0229 g/cm2) were lower than the other populations (Jar-ES p < 0.0001; Jar-MR p < 0.0001; Jar-RM p = 0.0491). Additionally, the Ruby Mountain population had slightly lower LMA than the Eastern Sierra population and was signiﬁcantly smaller than the Mount Rose LMA measurements (RM-MR p = 0.0146).
Do leaf quantity and LMA change through time within four populations of the Great Basin and eastern Sierra? Of the four populations and two traits, only one leaf trait in one population exhibited change over time (Figure 3). The Mount Rose population showed a signiﬁcant increase in leaf quantity (mean historical = 826.7, mean contemporary = 1168.6; p = 0.0312). The other populations showed no signiﬁcant change of leaf quantity (Appendix Table 3). Leaf mass per area did not show a signiﬁcant change over time within any population. Overall, the contemporary samples trended towards having more leaves than the historical samples, though this diﬀerence was not signiﬁcant at the p = 0.05 level.
Understanding how leaf traits vary among populations is an important ﬁrst step towards understanding how species respond to their environments and for understanding the degree of variation in traits across landscapes. In our study, we measured two leaf traits across multiple populations that diﬀered in climatic conditions and how traits change across time. We found that the driest whitebark pine site, Jarbidge, had both lower leaf quantity and LMA compared to other populations. Lower leaf quantities decrease surface area for evapotranspiration, making plants more water-use eﬃcient (27). However, decreased leaf mass per area has an inverse eﬀect. Decreased LMA is associated with lower densities of mesophyll tissue within the leaves, making plants less drought-tolerant (28, 29). This is of concern for Jarbidge given its notably low precipitation of the four sites and because similar LMA decreases in Scots pine were observed during drought-induced mortality (23).
The only population to show a change between historical and contemporary samples was Mount Rose, which showed an increase in leaf quantity over 105 years. This change in leaf quantity was not driven solely by the range of collection years within the Mount Rose population, as the exclusions of the oldest and most recent specimens yield comparable results. The exclusion of the most recent specimen does not inﬂuence the interpretation of the analysis for any of the populations. Mount Rose had the highest precipitation (Table 2) which has been increasing over time (20). Climate change may be driving the increasing precipitation in this area and lengthening the growing season, allowing for trees to accumulate more biomass. Increases in rainfall caused by climate change occur through the increased water holding capacity of a warming atmosphere (30). However, the interactions between temperature, precipitation, and plant traits are not always easily predictable (31). Although some interactions may extend the period for growth due to more favorable climatic conditions (32), climate islands may be forcing subalpine plants further up mountain tops. As trees increase in altitude there is less available land surface area, decreasing suitable habitat (33). The increase in leaf quantity we observed at Mount Rose may have ecological implications as it relates to leaf area index (LAI = total crown leaf area/ground area). Increased leaf quantity increases the total leaf area, which raises LAI. Increases in LAI have been shown to interact with hydrologic functioning by intercepting and sublimating snowpack (34). Increased canopy structure may also improve habitat for a variety of bird species (35), which are often used as surrogates for ecosystem diversity (36).
Our study indicates that there is variation in the leaf traits of contemporary populations of whitebark pine. Speciﬁcally, the number of leaves and the LMA from individuals at Jarbidge are lower than other populations. Further research is needed to understand the source of this variability, which could be due to local adaptation or genetic drift, as this population is relatively isolated and may experience limited gene ﬂow. Additionally, the Jarbidge population has experienced high beetle mortality since 2008 (37), making the future of these trees uncertain. These unique diﬀerences at Jarbidge, along with the history of drought and mountain pine beetles, make this isolated population a research priority. A reciprocal transplant or common garden experiment would be able to measure the p of local adaptation within this species. Understanding how whitebark pine is changing across its range will be critical for future restoration success.
We would like to thank Drs.Sarah Bisbing and Peter Weisberg at the University of Nevada, Reno for their advice and comments; Dr.Michele Slaton at the US Forest Service for her instrumental assistance in the early stages of this project. This project was supported by funding from the University of Nevada, Reno, and the Bristlecone Chapter of the California Native Plant Society (Alder2018).
 R. E. Keane et al., Gen. Tech. Rep. RMRS-GTR279. Fort Collins, CO: US Department of Agriculture, Forest Service, Rocky Mountain Research Station. 108, 279 (2012).
 S. F. Arno, R. J. Hoﬀ, Gen. Tech. Rep. INT-GTR253. Ogden, UT: US Department of Agriculture, Forest Service, Intermountain Research Station 11, 253 (1989).
 D. F. Tomback, P. Achuﬀ, Forest Pathology 40, 186– 225 (2010).
 K. Gibson et al., US Department of Agriculture Forest Service, Northern Region, Missoula, Montana 20, 8–20 (2008).
 R. E. Keane, P. Morgan, J. P. Menakis, Northwest Science 68 (1994).
 C. I. Millar et al., Canadian Journal of Forest Research 42, 749–765 (2012).
 N. McDowell et al., New Phytologist 178, 719–739 (2008).
 C. D. Allen et al., Forest Ecology and Management 259, 660–684 (2010).
 R. A. Andrus, B. J. Harvey, K. C. Rodman, S. J. Hart, T. T. Veblen, Ecology 99, 567–575 (2018).
 M. D. Abrams, M. E. Kubiske, K. C. Steiner, Tree Physiology 6, 305–315 (1990).
 S. A. Anjum et al., African Journal of Agricultural Research 6, 2026–2032 (2011).
 H. M. Poulos, G. P. Berlyn, The Journal of the Torrey Botanical Society 134, 281–289 (2007).
 J. L. Funk et al., Biological Reviews 92, 1156–1173 (2017).
 E. A. Leger, Global Change Biology 19, 2229–2239 (2013).
 G. R. Guerin, H. Wen, A. J. Lowe, Biology Letters 8, 882–886 (2012).
 L. Benomar et al., Frontiers in Plant Science 7, 48 (2016).
 C. A. Jones, C. C. Daehler, PeerJ 6, e4576 (2018).
 J. K. McKay, C. E. Christian, S. Harrison, K. J. Rice, Restoration Ecology 13, 432–440 (2005).
 A. D. Bower, S. N. Aitken, American Journal of Botany 95, 66–76 (2008).
 PRISM Climate Group, Oregon State University http://prism.oregonstate.edu, created 12 Apr 2019.
 P. A. Vesk, M. Westoby, The New Phytologist 160, 7–14 (2003).
 I. J. Wright, P. B. Reich, M. Westoby, et al., Nature 428, 821–827 (2004).
 R. Poyatos, D. Aguad´e, L. Galiano, M. Mencuccini, J. Martı´enez-Vilalta, New Phytologist 200, 388–401 (2013).
 Intermountain Regional Herbarium Network. 2019. http://intermountainbiota.org/portal/. Downloaded on 30 August 2018.
 K. R. Hultine, J. D. Marshall, Oecologia 123, 32–40 (2000).
 R Core Team. 2019. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. http://www.Rproject.org.
 N. Br´eda, R. Huc, A. Granier, E. Dreyer, Annals of Forest Science 63, 625–644 (2006).
 H. Poorter, U. Niinemets, L. Poorter, I. J. Wright, R. Villar, New Phytologist 182, 565–588 (2009).
 E. G. de la Riva, M. Olmo, H. Poorter, J. L. Ubera, R. Villar, PLoS ONE 11, e0148788 (2016).
 K. E. Trenberth, Climate Research 47, 123–138 (2011).
 M. Westoby, D. S. Falster, A. T. Moles, P. A. Vesk, I. J. Wright, Annual Review of Ecology and Systematics 33, 125–159 (2002).
 T. Hwang et al., Water Resources Research 54, 5359– 5375 (2018).
 W. Romme, M. Turner, Conservation Biology 5, 373– 386 (1991).
 A. N. Gelfan, J. W. Pomeroy, L. S. Kuchment, Journal of Hydrometeorology 5, 785–803 (2004).
 C. E. Swift, K. T. Vierling, A. T. Hudak, L. A. Vierling, Canadian Journal of Remote Sensing 43, 231–243 (2017).
 F. Larsen, J. Bladt, A. Balmford, C. Rahbek, Journal of Applied Ecology 49, 349–356 (2012).
 G. Durham et al., Forest Pest Conditions in Nevada 2010 (Nevada Division of Forestry, 2011).