Seasonality and phenology alter functional leaf traits

JA, Larson A, Townsend PA, Montgomery R, Cavender‐Bares J, Reich PB. All the light we cannot see: Climate manipulations leave short and long-term imprints in spectral reflectance of trees

Anthropogenic climate change, particularly changes in temperature and precipitation, affects plants in multiple ways. Because plants respond dynamically to stress and acclimate to changes in growing conditions, diagnosing quantitative plant-environment relationships is a major challenge. One approach to this problem is to quantify leaf responses using spectral reflectance, which provides rapid, inexpensive, and nondestructive measurements that capture a wealth of information about genotype as well as phenotypic responses to the environment. However, it is unclear how warming and drought affect spectra. To address this gap, we used an open-air field experiment that manipulates temperature and rainfall in 36 plots at two sites in the boreal-temperate ecotone of northern Minnesota, USA. We collected leaf spectral reflectance (400-2400 nm) at the peak of the growing season for three consecutive years on juveniles (two to six years old) of five tree species planted within the experiment. We hypothesized that these mid-season measurements of spectral reflectance capture a snapshot of the leaf phenotype encompassing a suite of physiological, structural, and biochemical responses to both long- and short-time scale environmental conditions. We show that the imprint of environmental conditions experienced by plants hours to weeks before spectral measurements is linked to regions in the spectrum associated with stress, namely the water absorption regions of the near-infrared and short-wave infrared. In contrast, the environmental conditions plants experience during leaf development leave lasting imprints on the spectral profiles of leaves, attributable to leaf structure and chemistry (e.g., pigment content and associated ratios). Our analyses show that after accounting for baseline species spectral differences, spectral responses to the environment do not differ among the species. This suggests that building a general framework for understanding forest responses to climate change through spectral metrics may be possible, likely having broader implications if the common responses among species detected here represent a widespread phenomenon. Consequently, these results demonstrate that examining the entire spectrum of leaf reflectance for environmental imprints in contrast to single features (e.g., indices and traits) improves inferences about plant-environment relationships, which is particularly important in times of unprecedented climate change.

Photosynthetic traits of Phragmites australis along an ecological gradient and developmental stages

Common reed (Phragmites australis) is a cosmopolitan species, though its dieback is a worldwide phenomenon. In order to assess the evolutionary role of phenotypic plasticity in a successful plant, the values and plasticity of photophysiological traits of Phragmites australis were investigated in the Lake Fertő wetlands at 5 sites with different degrees of reed degradation and along a seasonal sequence. On the one hand, along the established ecological degradation gradient, photophysiological traits of Phragmites changed significantly, affecting plant productivity, although no consistent gradient-type trends were observed. Gradual changes within a season in the values of photosynthetic traits were observed that were recorded in both degraded and stable stands, suggesting a universal response to seasonally changing environmental conditions that could not be overridden by the ecological gradient. On the other hand, reed plants exposed to different levels of degradation showed comparable physiological plasticity; there was no difference in trait variability between stable and degraded stands. This relatively uniform plasticity is likely to contribute to the resilience of reed plants by providing a wider range of adaptive traits under different conditions. In contrast, the 150-200% gradual change in photophysiological trait plasticity with senescence in Phragmites was also demonstrated, reflecting a more dynamic response of the photosynthetic apparatus to seasonal changes. Senescence affected the plasticity of plant traits independently of their degradation status, suggesting a more universal nature of seasonal changes. This research shows that under conditions of conservative resource use determined by stressful habitats, trait values respond to conditions, while trait plasticity shows minimal changes. Furthermore, phenological sequence significantly influenced both the values and the plasticity of the photosynthetic traits studied. Our results underline the impact of ecological degradation on reed physiology and highlight the importance of understanding both trait values and plasticity in plant responses to environmental and seasonal change.

Unraveling the drivers and impacts of leaf phenological diversity in a subtropical forest: A fine-scale analysis using PlanetScope CubeSats

Leaf phenology variations within plant communities shape community assemblages and influence ecosystem properties and services. However, questions remain regarding quantification, drivers, and productivity impacts of intra-site leaf phenological diversity. With a 50-ha subtropical forest plot in China's Heishiding Provincial Nature Reserve (part of the global ForestGEO network) as a testbed, we gathered a unique dataset combining ground-derived abiotic (topography, soil) and biotic (taxonomic diversity, functional diversity, functional traits) factors. We investigated drivers underlying leaf phenological diversity extracted from high-resolution PlanetScope data, and its influence on aboveground biomass (AGB) using structural equation modeling (SEM). Our results reveal considerable fine-scale leaf phenological diversity across the subtropical forest landscape. This diversity is directly and indirectly influenced by abiotic and biotic factors (e.g. slope, soil, traits, taxonomic diversity; r2 = 0.43). While a notable bivariate relationship between AGB and leaf phenological diversity was identified (r = -0.24, P < 0.05), this relationship did not hold in SEM analysis after considering interactions with other biotic and abiotic factors (P > 0.05). These findings unveil the underlying mechanism regulating intra-site leaf phenological diversity. While leaf phenology is known to be associated with ecosystem properties, our findings confirm that AGB is primarily influenced by functional trait composition and taxonomic diversity rather than leaf phenological diversity.

Unveiling the transferability of PLSR models for leaf trait estimation: lessons from a comprehensive analysis with a novel global dataset

Leaf traits are essential for understanding many physiological and ecological processes. Partial least squares regression (PLSR) models with leaf spectroscopy are widely applied for trait estimation, but their transferability across space, time, and plant functional types (PFTs) remains unclear. We compiled a novel dataset of paired leaf traits and spectra, with 47 393 records for > 700 species and eight PFTs at 101 globally distributed locations across multiple seasons. Using this dataset, we conducted an unprecedented comprehensive analysis to assess the transferability of PLSR models in estimating leaf traits. While PLSR models demonstrate commendable performance in predicting chlorophyll content, carotenoid, leaf water, and leaf mass per area prediction within their training data space, their efficacy diminishes when extrapolating to new contexts. Specifically, extrapolating to locations, seasons, and PFTs beyond the training data leads to reduced R2 (0.12-0.49, 0.15-0.42, and 0.25-0.56) and increased NRMSE (3.58-18.24%, 6.27-11.55%, and 7.0-33.12%) compared with nonspatial random cross-validation. The results underscore the importance of incorporating greater spectral diversity in model training to boost its transferability. These findings highlight potential errors in estimating leaf traits across large spatial domains, diverse PFTs, and time due to biased validation schemes, and provide guidance for future field sampling strategies and remote sensing applications.

Unraveling the ecological threads: How invasive alien plants influence soil carbon dynamics

Invasive alien plants (IAPs) pose significant threats to native ecosystems and biodiversity worldwide. However, the understanding of their precise impact on soil carbon (C) dynamics in invaded ecosystems remains a crucial area of research. This review comprehensively explores the mechanisms through which IAPs influence soil C pools, fluxes, and C budgets, shedding light on their effects and broader consequences. Key mechanisms identified include changes in litter inputs, rates of organic matter decomposition, alterations in soil microbial communities, and shifts in nutrient cycling, all driving the impact of IAPs on soil C dynamics. These mechanisms affect soil C storage, turnover rates, and ecosystem functioning. Moreover, IAPs tend to increase gross primary productivity and net primary productivity leading to the alterations in fluxes and C budgets. The implications of IAP-induced alterations in soil C dynamics are significant and extend to plant-soil interactions, ecosystem structure, and biodiversity. Additionally, they have profound consequences for C sequestration, potentially impacting climate change mitigation. Restoring native plant communities, promoting soil health, and implementing species-specific management are essential measures to significantly mitigate the impacts of IAPs on soil C dynamics. Overall, understanding and mitigating the effects of IAPs on soil C storage, nutrient cycling, and related processes will contribute to the conservation of native biodiversity and complement global C neutrality efforts.

A trait-based investigation into evergreen woody plants for traffic-related air pollution mitigation over time

This study investigated influences of leaf traits on particulate matter (PM) wash-off and (re)capture (i.e., net removal) over time. Leaf samples were taken before and after three rainfall events from a range of 10 evergreen woody plants (including five different leaf types), which were positioned with an optical particle counter alongside a busy road. Scanning electron microscopy was used to quantify the density (no./mm2), mass (μg/cm2), and elemental composition of deposited particles. To enable leaf area comparison between scale-like leaves and other leaf types, a novel metric (FSA: foliage surface area per unit branch length) was developed, which may be utilised by future research. Vehicle-related particles constituted 15 % of total deposition, and there was a notable 50 % decrease in the proportion of tyre wear particles after rainfall. T. baccata presented the lowest proportion (11.1 %) of vehicle-related particle deposition but the most consistent performance in terms of net PM removal. Only four of the 10 plant specimens (C. japonica, C. lawsoniana, J. chinensis, and T. baccata) presented effective PM wash-off across all particle size fractions and rainfall intensities, with a generally positive relationship observed between rainfall intensity and wash-off. Mass deposition was more significantly determined by particle size than number density. Interestingly, larger particles were also less easily washed off than smaller particles. Some traits typically considered to be advantageous (e.g., greater hairiness) may in fact hinder net removal over time due to retention under rainfall. Small leaf area is one trait that may promote both accumulation and wash-off. However, FSA was found to be the most influential trait, with an inverse relationship between FSA and wash-off efficacy. This finding poses trade-offs and opportunities for green infrastructure design, which are discussed. Finally, numerous areas for future research are recommended, underlining the importance of systems approaches in developing vegetation management frameworks.

The leaf economics spectrum of triploid and tetraploid C4 grass Miscanthus x giganteus

The leaf economics spectrum (LES) describes multivariate correlations in leaf structural, physiological and chemical traits, originally based on diverse C₃ species grown under natural ecosystems. However, the specific contribution of C₄ species to the global LES is studied less widely. C₄ species have a CO₂ concentrating mechanism which drives high rates of photosynthesis and improves resource use efficiency, thus potentially pushing them towards the edge of the LES. Here, we measured foliage morphology, structure, photosynthesis, and nutrient content for hundreds of genotypes of the C₄ grass Miscanthus× giganteus grown in two common gardens over two seasons. We show substantial trait variations across M.× giganteus genotypes and robust genotypic trait relationships. Compared to the global LES, M.× giganteus genotypes had higher photosynthetic rates, lower stomatal conductance, and less nitrogen content, indicating greater water and photosynthetic nitrogen use efficiency in the C₄ species. Additionally, tetraploid genotypes produced thicker leaves with greater leaf mass per area and lower leaf density than triploid genotypes. By expanding the LES relationships across C₃ species to include C₄ crops, these findings highlight that M.× giganteus occupies the boundary of the global LES and suggest the potential for ploidy to alter LES traits.

Assisted migration is plausible for a boreal tree species under climate change: A quantitative and population genetics study of trembling aspen (Populus tremuloides Michx.) in western Canada

A novel method was tested for improving tree breeding strategies that integrate quantitative and population genetics based on range-wide reciprocal transplant experiments. Five reciprocal common garden tests of Populus tremuloides were investigated including 6450 trees across western Canada focusing on adaptation traits and growth. Both genetic parameters and home-site transplant models were evaluated. We found a genetic trade-off between growth and early spring leaf flush and late fall senescence. Coefficients of phenotypic variation (CVp) of cell lysis (CL), a measure of freezing injury, shrank from 0.28 to 0.10 during acclimation in the fall, and the CVp slope versus the freezing temperature was significantly different from zero (R 2 = 0.33, p = .02). There was more between-population genetic variation in fall phenology than in spring leaf phenology. We suggest that P. tremuloides demonstrated a discrepancy between the ecological optimum and the physiological optimum minimum winter temperature. The sub-optimal growing condition of P. tremuloides is potentially caused by the warmer ecological optimum than the physiological optimum. Assisted migration and breeding of fast growers to reforest cooler plantation sites can improve productivity. Transferring the study populations to less than 4°C of extreme minimum temperature appears safe for reforestation aligning with the historical recolonization direction of the species. This is equivalent to a 5-10° latitudinal northward movement. Fall frost hardiness is an effective criterion for family selection in the range tested in this study.

Litter quality and decomposition responses to drought in a northeastern US deciduous forest

Even though drought impacts on tree physiology have been identified, whether drought affects leaf litter chemistry that, in turn, influences litter decay rates is still poorly understood. We compared litter quality and decomposition for two cohorts of leaves from five co-occurring seasonally deciduous tree species: Acer saccharum, Tilia americana, Quercus rubra, Quercus alba, and Ostrya virginiana. One cohort experienced a growing-season drought, and the other cohort came from the same trees in the ensuing, post-drought growing season. Leaf litter production was greater for drought litter than post-drought litter for all five species. Specific leaf area and nitrogen concentrations were 20% greater for the drought cohort than the post-drought cohort. Concentrations of non-structural carbohydrates were about 14% greater for the drought cohort, except for greater values for post-drought A. saccharum litter. Pectin in the middle lamella of leaf litter was 31% lower for the drought cohort compared to post-drought cohort. We found few differences in litter decay rates between drought and post-drought cohorts, although Q. rubra litter had more decomposition for the post-drought cohort than the drought cohort, whereas A. saccharum litter had more decomposition for the drought cohort than the post-drought cohort. Leaf litter decay rates for the drought cohort were related to litter nitrogen and lignin concentrations, whereas decay rates for the post-drought cohort were related to litter carbohydrate concentrations. Our findings suggest that the role of drought events on seasonally deciduous forest ecosystems must recognize species-specific, idiosyncratic responses in leaf litter quality and decomposition.

New Insight into Genetic Structure and Diversity of Scots Pine (Pinus sylvestris L.) Populations in Lithuania Based on Nuclear, Chloroplast and Mitochondrial DNA Markers

We studied the genetic differentiation, structure, and diversity of Scots pine populations in Lithuania based on nuclear, chloroplast microsatellite, and mitochondrial DNA markers. We focused on revealing evolutionary history, country-wide geneflow patterns, and structuring among the Scots pine populations. We genotyped 439 Scots pine individuals of mature age from 23 natural Scots pine stands in Lithuania and used the AMOVA and a set of genetic-clustering methods. The among-population differentiation was weak for nuclear microsatellite loci (nSSRs) (FST = 0.005) but much stronger for cpSSRs (PhiST = 0.240). The populations were structured into highland and lowland populations based on cpSSRs and eastern highland versus the rest for nSSRs. We detected two mtDNA mitotypes—the universal type A and northeastern type B, and the latter occurred at a markedly higher frequency in eastern Lithuania. Within-population genetic diversity was higher in large pine-dominated forest tracts in the eastern highlands than in fragmented forests in the western highlands. We concluded that phenology-based genetic networks following the temperature climate gradients have a strong effect on shaping the genetic structure of otherwise rather homogeneous gene pools of Scots pine populations in Lithuania. The possible effects of human interference with forests on genetic diversity of Scots pine populations in Lithuania are discussed.

Forest phenoclusters for Argentina based on vegetation phenology and climate

Forest biodiversity conservation and species distribution modeling greatly benefit from broad-scale forest maps depicting tree species or forest types rather than just presence and absence of forest, or coarse classifications. Ideally, such maps would stem from satellite image classification based on abundant field data for both model training and accuracy assessments, but such field data do not exist in many parts of the globe. However, different forest types and tree species differ in their vegetation phenology, offering an opportunity to map and characterize forests based on the seasonal dynamic of vegetation indices and auxiliary data. Our goal was to map and characterize forests based on both land surface phenology and climate patterns, defined here as forest phenoclusters. We applied our methodology in Argentina (2.8 million km² ), which has a wide variety of forests, from rainforests to cold-temperate forests. We calculated phenology measures after fitting a harmonic curve of the enhanced vegetation index (EVI) time series derived from 30-m Sentinel 2 and Landsat 8 data from 2018-2019. For climate, we calculated land surface temperature (LST) from Band 10 of the thermal infrared sensor (TIRS) of Landsat 8, and precipitation from Worldclim (BIO12). We performed stratified X-means cluster classifications followed by hierarchical clustering. The resulting clusters separated well into 54 forest phenoclusters with unique combinations of vegetation phenology and climate characteristics. The EVI 90th percentile was more important than our climate and other phenology measures in providing separability among different forest phenoclusters. Our results highlight the potential of combining remotely sensed phenology measures and climate data to improve broad-scale forest mapping for different management and conservation goals, capturing functional rather than structural or compositional characteristics between and within tree species. Our approach results in classifications that go beyond simple forest-nonforest in areas where the lack of detailed ecological field data precludes tree species-level classifications, yet conservation needs are high. Our map of forest phenoclusters is a valuable tool for the assessment of natural resources, and the management of the environment at scales relevant for conservation actions.

Predictability of leaf traits with climate and elevation: a case study in Gongga Mountain, China

Leaf mass per area (Ma), nitrogen content per unit leaf area (Narea), maximum carboxylation capacity (Vcmax) and the ratio of leaf-internal to ambient CO2 partial pressure (χ) are important traits related to photosynthetic function, and they show systematic variation along climatic and elevational gradients. Separating the effects of air pressure and climate along elevational gradients is challenging due to the covariation of elevation, pressure and climate. However, recently developed models based on optimality theory offer an independent way to predict leaf traits and thus to separate the contributions of different controls. We apply optimality theory to predict variation in leaf traits across 18 sites in the Gongga Mountain region. We show that the models explain 59% of trait variability on average, without site- or region-specific calibration. Temperature, photosynthetically active radiation, vapor pressure deficit, soil moisture and growing season length are all necessary to explain the observed patterns. The direct effect of air pressure is shown to have a relatively minor impact. These findings contribute to a growing body of research indicating that leaf-level traits vary with the physical environment in predictable ways, suggesting a promising direction for the improvement of terrestrial ecosystem models.

Transcriptome-wide analysis of introgression-resistant regions reveals genetic divergence genes under positive selection in Populus trichocarpa

Comparing gene expression patterns and genetic polymorphisms between populations is of central importance for understanding the origin and maintenance of biodiversity. Based on population-specific gene expression levels and allele frequency differences, we sought to identify population divergence (PD) genes across the introgression-resistant genomic regions of Populus trichocarpa. Genes containing highly diverged loci [i.e., genetic divergence (GD)] or showing expression divergence (ED) between populations were widely distributed in the genome and substantially enriched in functional categories related to stress responses, disease resistance, timing of flowering, cell cycle regulation, plant growth, and development. Nine genomic regions showing evidence of strong positive selection were overlapped with GD genes, which had significant differences between Oregon (a southernmost peripheral deme) and the other demes. However, we did not find evidence that genes under positive selection show an enrichment for ED. PD genes and genes under selection pertained to the same gene classes, such as SERINE/CYSTEINE PROTEASE, ABC TRANSPORTER, GLYCOSYLTRANSFERASE and other transferases. Our analysis also revealed that GD genes were polymorphic within the species (41.9 ± 3.66 biallelic variants per gene), as previously reported in herbaceous plants. By contrast, ED genes contained less genetic variants (10.73 ± 1.14) and were likely highly expressed. In addition, we found that trans- rather than cis-acting variants considerably contribute to the evolution of >90% PD genes. Overall, this study elucidates that cohorts of PD genes agree with the general attributes of known speciation genes and GD genes will provide substrates for positive selection to operate on.

Comparison of the Scaling Relationships of Leaf Biomass versus Surface Area between Spring and Summer for Two Deciduous Tree Species

The scaling relationship between either leaf dry or fresh mass (M) and surface area (A) can reflect the photosynthetic potential and efficiency of light harvesting in different broad-leaved plants. In growing leaves, lamina area expansion is typically finished before the completion of leaf biomass accumulation, thereby affecting the M vs. A scaling relationship at different developmental stages of leaves (e.g., young vs. adult leaves). In addition, growing plants can have different-sized leaves at different plant ages, potentially also changing M vs. A scaling. Furthermore, leaf shape can also change during the course of ontogeny and modify the M vs. A scaling relationship. Indeed, the effect of seasonal changes in leaf shape on M vs. A scaling has not been examined in any previous studies known to us. The study presented here was conducted using two deciduous tree species: Alangium chinense (saplings forming leaves through the growing season) and Liquidambar formosana (adult trees producing only one leaf flush in spring) that both have complex but nearly bilaterally symmetrical leaf shapes. We determined (i) whether leaf shapes differed in spring versus summer; (ii) whether the M vs. A scaling relationship varied over time; and (iii) whether there is a link between leaf shape and the scaling exponent governing the M vs. A scaling relationship. The data indicated that (i) the leaf dissection index in spring was higher than that in summer for both species (i.e., leaf-shape complexity decreased from young to adult leaves); (ii) there was a significant difference in the numerical value of the scaling exponent of leaf perimeter vs. area between leaves sampled at the two dates; (iii) spring leaves had a higher water content than summer leaves, and the scaling exponents of dry mass vs. area and fresh mass vs. area were all greater than unity; (iv) the scaling relationship between fresh mass and area was statistically more robust than that between leaf dry mass and area; (v) the scaling exponents of leaf dry and fresh mass vs. area of A. chinense leaves in spring were greater than those in summer (i.e., leaves in younger plants tend to be larger than leaves in older plants), whereas, for the adult trees of L. formosana, the scaling exponent in spring was smaller than that in summer, indicating increases in leaf dry mass per unit area with increasing leaf age; and (vi) leaf shape appears not to be related to the scaling relationship between either leaf dry or fresh mass and area, but is correlated with the scaling exponent of leaf perimeter vs. area (which tends to be a ½ power function). These trends indicate that studies of leaf morphometrics and scaling relationships must consider the influence of seasonality and plant age in sampling of leaves and the interpretation of data.

Monitoring Plant Functional Diversity Using the Reflectance and Echo from Space

Plant functional diversity (FD) is an important component of biodiversity. Evidence shows that FD strongly determines ecosystem functioning and stability and also regulates various ecosystem services that underpin human well-being. Given the importance of FD, it is critical to monitor its variations in an explicit manner across space and time, a highly demanding task that cannot be resolved solely by field data. Today, high hopes are placed on satellite-based observations to complement field plot data. The promise is that multiscale monitoring of plant FD, ecosystem functioning, and their services is now possible at global scales in near real-time. However, non-trivial scale challenges remain to be overcome before plant ecology can capitalize on the latest advances in Earth Observation (EO). Here, we articulate the existing scale challenges in linking field and satellite data and further elaborated in detail how to address these challenges via the latest innovations in optical and radar sensor technologies and image analysis algorithms. Addressing these challenges not only requires novel remote sensing theories and algorithms but also urges more effective communication between remote sensing scientists and field ecologists to foster mutual understanding of the existing challenges. Only through a collaborative approach can we achieve the global plant functional diversity monitoring goal.

Differences in growth and gas exchange between southern and northern provenances of silver birch (Betula pendula Roth) in northern Europe

Due to its ubiquity across northern latitudes, silver birch (Betula pendula Roth) is an attractive model species for studying geographical trait variation and acclimation capacity. Six birch provenances from 60 to 67°N across Finland were grown in a common garden and studied for provenance and genotype variation. We looked for differences in height growth, photosynthetic gas exchange and chlorophyll content index (CCI) and compared the gas exchange of early and late leaves on short and long shoots, respectively. The provenances stratified into southern and northern groups. Northern provenances attained less height growth increment and had higher stomatal conductance (gs) and lower intrinsic water-use efficiency (WUE, Anet/gs) than southern provenances, whereas net photosynthesis (Anet) or CCI did not show clear grouping. Short shoot leaves had lower gs and higher WUE than long shoot leaves in all provenances, but there was no difference in Anet between shoot types. The separation of the provenances into two groups according to their physiological responses might reflect the evolutionary history of B. pendula. Latitudinal differences in gas exchange and water use traits can have plausible consequences for global carbon and water fluxes in a warming climate.

Phenotypic plasticity of natural Populus trichocarpa populations in response to temporally environmental change in a common garden

BACKGROUND

Natural selection on fitness-related traits can be temporally heterogeneous among populations. As climate changes, understanding population-level responses is of scientific and practical importance. We examined 18 phenotypic traits associated with phenology, biomass, and ecophysiology in 403 individuals of natural Populus trichocarpa populations, growing in a common garden.

RESULTS

Compared with tree origin settings, propagules likely underwent drought exposures in the common garden due to significantly low rainfall during the years of measurement. All study traits showed population differentiation reflecting adaptive responses due to local genetic adaptation. Phenology and biomass traits were strongly under selection and showed plastic responses between years, co-varying with latitude. While phenological events (e.g., bud set and growth period) and biomass were under positive directional selection, post-bud set period, particularly from final bud set to the onset of leaf drop, was selected against. With one exception to water-use efficiency, ecophysiology traits were under negative directional selection. Moreover, extended phenological events jointly evolved with source niches under increased temperature and decreased rainfall exposures. High biomass coevolved with climatic niches of high temperature; low rainfall promoted high photosynthetic rates evolution.

CONCLUSIONS

This work underpins that P. trichocarpa is likely to experience increased fitness (height gain) by evolving toward extended bud set and growth period, abbreviated post-bud set period, and increased drought resistance, potentially constituting a powerful mechanism for long-lived tree species in surviving unpredictably environmental extremes (e.g., drought).

Visible and near-infrared hyperspectral indices explain more variation in lower-crown leaf nitrogen concentrations in autumn than in summer

Autumn canopy phenological transitions are increasing in length as a consequence of climate change. Here, we assess how well hyperspectral indices in the visible and near-infrared (NIR) wavelengths predict nitrogen (N) concentrations in lower-canopy leaves in the autumn phenological transition as they are generally understudied in leaf trait research. Using a Bayesian framework, we tested how well published indices are able to predict N concentrations in Fagus grandifolia Ehrh., Liriodendron tulipifera L., and Betula lenta L. from mid-summer through senescence, and how related the indices are to autumn phenological change. No indices were able to determine a trend in differences in N in mid-summer leaves. Indices that included wavelengths in the green and NIR ranges were the first indices able to detect a trend and had among the highest correlations with N concentration in both the last green collection and the senescing collection. Models were unique when indices were fit to data from different phenophases. Indices that focused on only the red edge (i.e., the sharp increase in reflectance between the red and NIR wavelengths) had the strongest explanatory power across the autumn phenological transition, but had less explanatory power for individual collections. These indices, as well as those that have been correlated with chlorophyll (CCI) and carotenoids (PRI), were the strongest descriptors of autumn progression. This study provides insights on challenges and capabilities to monitor a leaf's N concentration throughout and across canopy senescence.

A role for SPEECHLESS in the integration of leaf stomatal patterning with the growth vs disease trade-off in poplar

Occurrence of stomata on both leaf surfaces (amphistomaty) promotes higher stomatal conductance and photosynthesis while simultaneously increasing exposure to potential disease agents in black cottonwood (Populus trichocarpa). A genome-wide association study (GWAS) with 2.2M single nucleotide polymorphisms generated through whole-genome sequencing found 280 loci associated with variation in adaxial stomatal traits, implicating genes regulating stomatal development and behavior. Strikingly, numerous loci regulating plant growth and response to biotic and abiotic stresses were also identified. The most significant locus was a poplar homologue of SPEECHLESS (PtSPCH1). Individuals possessing PtSPCH1 alleles associated with greater adaxial stomatal density originated primarily from environments with shorter growing seasons (e.g. northern latitudes, high elevations) or with less precipitation. PtSPCH1 was expressed in developing leaves but not developing stem xylem. In developing leaves, RNA sequencing showed patterns of coordinated expression between PtSPCH1 and other GWAS-identified genes. The breadth of our GWAS results suggests that the evolution of amphistomaty is part of a larger, complex response in plants. Suites of genes underpin this response, retrieved through genetic association to adaxial stomata, and show coordinated expression during development. We propose that the occurrence of amphistomaty in P. trichocarpa involves PtSPCH1 and reflects selection for supporting rapid growth over investment in immunity.

Leaf age effects on the spectral predictability of leaf traits in Amazonian canopy trees

Recent work has shown that leaf traits and spectral properties change through time and/or seasonally as leaves age. Current field and hyperspectral methods used to estimate canopy leaf traits could, therefore, be significantly biased by variation in leaf age. To explore the magnitude of this effect, we used a phenological dataset comprised of leaves of different leaf age groups -developmental, mature, senescent and mixed-age- from canopy and emergent tropical trees in southern Peru. We tested the performance of partial least squares regression models developed from these different age groups when predicting traits for leaves of different ages on both a mass and area basis. Overall, area-based models outperformed mass-based models with a striking improvement in prediction observed for area-based leaf carbon (C_area) estimates. We observed trait-specific age effects in all mass-based models while area-based models displayed age effects in mixed-age leaf groups for P_area and N_area. Spectral coefficients and variable importance in projection (VIPs) also reflected age effects. Both mass- and area-based models for all five leaf traits displayed age/temporal sensitivity when we tested their ability to predict the traits of leaves of other age groups. Importantly, mass-based mature models displayed the worst overall performance when predicting the traits of leaves from other age groups. These results indicate that the widely adopted approach of using fully expanded mature leaves to calibrate models that estimate remotely-sensed tree canopy traits introduces error that can bias results depending on the phenological stage of canopy leaves. To achieve temporally stable models, spectroscopic studies should consider producing area-based estimates as well as calibrating models with leaves of different age groups as they present themselves through the growing season. We discuss the implications of this for surveys of canopies with synchronised and unsynchronised leaf phenology.

Intraspecific variation in soy across the leaf economics spectrum

Background and Aims

Intraspecific trait variation (ITV) is an important dimension of plant ecological diversity, particularly in agroecosystems, where phenotypic ITV (within crop genotypes) is an important correlate of key agroecosystem processes including yield. There are few studies that have evaluated whether plants of the same genotype vary along well-defined axes of biological variation, such as the leaf economics spectrum (LES). There is even less information disentangling environmental and ontogenetic determinants of crop ITV along an intraspecific LES, and whether or not a plant's position along an intraspecific LES is correlated with reproductive output.

Methods

We sought to capture the extent of phenotypic ITV within a single cultivar of soy (Glycine max) - the world's most commonly cultivated legume - using a data set of nine leaf traits measured on 402 leaves, sampled from 134 plants in both agroforestry and monoculture management systems, across three distinct whole-plant ontogenetic stages (while holding leaf age and canopy position stable).

Key Results

Leaf traits covaried strongly along an intraspecific LES, in patterns that were largely statistically indistinguishable from the 'universal LES' observed across non-domesticated plants. Whole-plant ontogenetic stage explained the highest proportion of phenotypic ITV in LES traits, with plants progressively expressing more 'resource-conservative' LES syndromes throughout development. Within ontogenetic stages, leaf traits differed systematically across management systems, with plants growing in monoculture expressing more 'resource-conservative' trait syndromes: trends largely owing to an approximately ≥50% increases in leaf mass per area (LMA) in high-light monoculture vs. shaded agroforestry systems. Certain traits, particularly LMA, leaf area and maximum photosynthetic rates, correlated closely with plant-level reproductive output.

Conclusions

Phenotypic ITV in soy is governed by constraints in trait trade-offs along an intraspecific LES, which in turn (1) underpins plant responses to managed environmental gradients, and (2) reflects shifts in plant functional biology and resource allocation that occur throughout whole-plant ontogeny.

Ecological genomics of variation in bud-break phenology and mechanisms of response to climate warming in Populus trichocarpa

Spring bud-break phenology is a critical adaptive feature common to temperate perennial woody plants. Understanding the molecular underpinnings of variation in bud-break is important for elucidating adaptive evolution and predicting outcomes relating to climate change. Field and controlled growth chamber tests were used to assess population-wide patterns in bud-break from wild-sourced black cottonwood (Populus trichocarpa) genotypes. We conducted a genome-wide association study (GWAS) with single nucleotide polymorphisms (SNPs) derived from whole genome sequencing to test for loci underlying variation in bud-break. Bud-break had a quadratic relationship with latitude, where southern- and northern-most provenances generally broke bud earlier than those from central parts of the species' range. Reduced winter chilling increased population-wide variation in bud-break, whereas greater chilling decreased variation. GWAS uncovered 16 loci associated with bud-break. Phenotypic changes connected with allelic variation were replicated in an independent set of P. trichocarpa trees. Despite phenotypic similarities, genetic profiles between southern- and northern-most genotypes were dissimilar based on our GWAS-identified SNPs. We propose that the GWAS-identified loci underpin the geographical pattern in P. trichocarpa and that variation in bud-break reflects different selection for winter chilling and heat sum accumulation, both of which can be affected by climate warming.

Sexual homomorphism in dioecious trees: extensive tests fail to detect sexual dimorphism in Populus †

The evolution of sexual dimorphism and expansion of sex chromosomes are both driven through sexual conflict, arising from differing fitness optima between males and females. Here, we pair work in poplar (Populus) describing one of the smallest sex-determining regions known thus far in complex eukaryotes (~100 kbp) with comprehensive tests for sexual dimorphism using >1300 individuals from two Populus species and assessing 96 non-reproductive functional traits. Against expectation, we found sexual homomorphism (no non-reproductive trait differences between the sexes), suggesting that gender is functionally neutral with respect to non-reproductive features that affect plant survival and fitness. Combined with a small sex-determining region, we infer that sexual conflict may be effectively stymied or non-existent within these taxa. Both sexual homomorphism and the small sex-determining region occur against a background of strong environmental selection and local adaptation in Populus. This presents a powerful hypothesis for the evolution of dioecious species. Here, we suggest that environmental selection may be sufficient to suppress and stymy sexual conflict if it acts orthogonal to sexual selection, thereby placing limitations on the evolution of sexual dimorphism and genomic expansion of sex chromosomes.

Convergence in relationships between leaf traits, spectra and age across diverse canopy environments and two contrasting tropical forests

Leaf age structures the phenology and development of plants, as well as the evolution of leaf traits over life histories. However, a general method for efficiently estimating leaf age across forests and canopy environments is lacking. Here, we explored the potential for a statistical model, previously developed for Peruvian sunlit leaves, to consistently predict leaf ages from leaf reflectance spectra across two contrasting forests in Peru and Brazil and across diverse canopy environments. The model performed well for independent Brazilian sunlit and shade canopy leaves (R²  = 0.75-0.78), suggesting that canopy leaves (and their associated spectra) follow constrained developmental trajectories even in contrasting forests. The model did not perform as well for mid-canopy and understory leaves (R²  = 0.27-0.29), because leaves in different environments have distinct traits and trait developmental trajectories. When we accounted for distinct environment-trait linkages - either by explicitly including traits and environments in the model, or, even better, by re-parameterizing the spectra-only model to implicitly capture distinct trait-trajectories in different environments - we achieved a more general model that well-predicted leaf age across forests and environments (R²  = 0.79). Fundamental rules, linked to leaf environments, constrain the development of leaf traits and allow for general prediction of leaf age from spectra across species, sites and canopy environments.

Natural acetylation impacts carbohydrate recovery during deconstruction of Populus trichocarpa wood

BACKGROUND

Significant variation in the inherent degree of acetylation naturally exists in the xylem cell walls of Populus trichocarpa. During pretreatment, endogenous acetate hydrolyzes to acetic acid that can subsequently catalyze the breakdown of poplar wood, increasing the efficiency of biomass pretreatment.

RESULTS

Poplar genotypes varying in cell wall composition were pretreated in 0.3% H₂SO₄ in non-isothermal batch reactors. Acetic acid released from the wood was positively related to sugar release during pretreatment (R ≥ 0.9), and inversely proportional to the lignin content of the poplar wood (R = 0.6).

CONCLUSION

There is significant variation in wood chemistry among P. trichocarpa genotypes. This study elucidated patterns of cell wall deconstruction and clearly links carbohydrate solubilization to acetate release. Tailoring biomass feedstocks for acetate release could enhance pretreatment efficiencies.

Variation in fluxes estimated from nitrogen isotope discrimination corresponds with independent measures of nitrogen flux in Populus balsamifera L

Acquisition of mineral nitrogen by roots from the surrounding environment is often not completely efficient, in which a variable amount of leakage (efflux) relative to gross uptake (influx) occurs. The efflux/influx ratio (E/I) is, therefore, inversely related to the efficiency of nutrient uptake at the root level. Time-integrated estimates of E/I and other nitrogen-use traits may be obtainable from variation in stable isotope ratios or through compartmental analysis of tracer efflux (CATE) using radioactive or stable isotopes. To compare these two methods, Populus balsamifera L. genotypes were selected, a priori, for high or low nitrogen isotope discrimination. Vegetative cuttings were grown hydroponically, and E/I was calculated using an isotope mass balance model (IMB) and compared to E/I calculated using (15) N CATE. Both methods indicated that plants grown with ammonium had greater E/I than nitrate-grown plants. Genotypes with high or low E/I using CATE also had similarly high or low estimates of E/I using IMB, respectively. Genotype-specific means were linearly correlated (r = 0.77; P = 0.0065). Discrepancies in E/I between methods may reflect uncertainties in discrimination factors for the assimilatory enzymes, or temporal differences in uptake patterns. By utilizing genotypes with known variation in nitrogen isotope discrimination, a relationship between nitrogen isotope discrimination and bidirectional nitrogen fluxes at the root level was observed.

Phenological variation of leaf functional traits within species

A basic assumption of the trait-based approach in plant ecology is that differences in functional trait values are greater between species than within species. We questioned this assumption by assessing (1) the relative extent of inter- and intraspecific leaf trait variation throughout a complete growing season (phenological variation) in a group of deciduous and evergreen woody species, and (2) whether species rankings based on leaf traits were maintained across the growing season. We analysed leaf mass per area (LMA) and leaf nutrient concentrations (C, N, P), including the C:N and N:P ratios. Intraspecific trait variation (ITV) due to phenology was significantly greater than interspecific variation for leaf N concentration on a mass basis (Nm; 68.90 %) and for the leaf C:N ratio (60.60 %), whereas interspecific variation was significantly higher than ITV for LMA (62.30 %) and for leaf C concentration on a mass (Cm) and area (Ca) basis (Cm 70.40 %; Ca 65.30 %). ITV was particularly low for LMA (<20 %). Species rankings were highly modified by phenology for a number of leaf traits (Pm, N:P ratio) but were relatively well conserved throughout the growing season for others (LMA, Nm). Patterns of ITV across the growing season differed significantly between deciduous and evergreen species for all traits except leaf P but did not vary between native and exotic species. Overall, our results show that intraspecific phenological variation in leaf traits may be similar to or greater than interspecific variation and that temporal patterns of ITV vary considerably among traits and species, especially for leaf nutrient concentrations, factors which can potentially affect quantitative interspecific relationships.

Evolutionary Quantitative Genomics of Populus trichocarpa

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

Higher photosynthetic capacity from higher latitude: foliar characteristics and gas exchange of southern, central and northern populations of Populus angustifolia

Narrowleaf cottonwood (Populus angustifolia James) is an obligate riparian poplar that is a foundation species in river valleys along the Rocky Mountains, spanning 16° of latitude from southern Arizona, USA to southern Alberta, Canada. Its current distribution is fragmented, and genetic variation shows regional population structure consistent with the effects of geographic barriers and past climate. It is thus very well-suited for investigating ecophysiological adaptation associated with latitude. In other section Tacamahaca poplar species, genotypes from higher latitudes show evidence of short-season adaptation with foliar traits that contribute to higher photosynthetic capacity. We tested for similar adaptation in three populations of narrowleaf cottonwoods: from Arizona (south), Alberta (north) and Utah, near the centre of the latitudinal distribution. We propagated 20 genotypes from each population in a common garden in Alberta, and measured foliar and physiological traits after 3 years. Leaves of genotypes from the northern population had higher leaf mass per area (LMA), increased nitrogen (N) content and higher carotenoid and chlorophyll content, and these were associated with higher light-saturated net photosynthesis (Asat). In leaves of all populations the majority of stomata were abaxial, with the proportion of abaxial stomata highest in the southern population. Stomatal conductance (gs) and transpiration rates were higher in the northern population but water-use efficiency (Asat/gs) and leaf carbon isotope composition (δ(13)C) did not differ across the populations. These results (i) establish links between Asat and gs, N, chlorophyll and LMA among populations within this species, (ii) are consistent with the discrimination of populations from prior investigation of genetic variation and (iii) support the concept of latitudinal adaptation, whereby deciduous trees from higher latitudes display higher photosynthetic capacity, possibly compensating for a shorter and cooler growth season and reduced insolation.

Population genetics of freeze tolerance among natural populations of Populus balsamifera across the growing season

Protection against freeze damage during the growing season influences the northern range limits of plants. Freeze tolerance and freeze avoidance are the two major freeze resistance strategies. Winter survival strategies have been extensively studied in perennials, but few have addressed them and their genetic basis during the growing season. We examined intraspecific phenotypic variation in freeze resistance of Populus balsamifera across latitude and the growing season. To investigate the molecular basis of this variation, we surveyed nucleotide diversity and examined patterns of gene expression in the poplar C-repeat binding factor (CBF) gene family. Foliar freeze tolerance exhibited latitudinal and seasonal variation indicative of natural genotypic variation. CBF6 showed signatures of recent selective sweep. Of the 46 SNPs surveyed across the six CBF homologs, only CBF2_619 exhibited latitudinal differences consistent with increased freeze tolerance in the north. All six CBF genes were cold inducible, but showed varying patterns of expression across the growing season. Some Poplar CBF homologs exhibited patterns consistent with historical selection and clinal variation in freeze tolerance documented here. However, the CBF genes accounted for only a small amount of the variation, indicating that other genes in this and other molecular pathways likely play significant roles in nature.

Comparative physiology of allopatric Populus species: geographic clines in photosynthesis, height growth, and carbon isotope discrimination in common gardens

Populus species with wide geographic ranges display strong adaptation to local environments. We studied the clinal patterns in phenology and ecophysiology in allopatric Populus species adapted to similar environments on different continents under common garden settings. As a result of climatic adaptation, both Populus tremula L. and Populus balsamifera L. display latitudinal clines in photosynthetic rates (A), whereby high-latitude trees of P. tremula had higher A compared to low-latitude trees and nearly so in P. balsamifera (p = 0.06). Stomatal conductance (g s) and chlorophyll content index (CCI) follow similar latitudinal trends. However, foliar nitrogen was positively correlated with latitude in P. balsamifera and negatively correlated in P. tremula. No significant trends in carbon isotope composition of the leaf tissue (δ(13)C) were observed for both species; but, intrinsic water-use efficiency (WUEi) was negatively correlated with the latitude of origin in P. balsamifera. In spite of intrinsically higher A, high-latitude trees in both common gardens accomplished less height gain as a result of early bud set. Thus, shoot biomass was determined by height elongation duration (HED), which was well approximated by the number of days available for free growth between bud flush and bud set. We highlight the shortcoming of unreplicated outdoor common gardens for tree improvement and the crucial role of photoperiod in limiting height growth, further complicating interpretation of other secondary effects.

Introducing a sensor to measure budburst and its environmental drivers

Budburst is a key adaptive trait that can help us understand how plants respond to a changing climate from the molecular to landscape scale. Despite this, acquisition of budburst data is constrained by a lack of information at the plant scale on the environmental stimuli associated with the release of bud dormancy. Additionally, to date, little effort has been devoted to phenotyping plants in natural populations due to the challenge of accounting for the effect of environmental variation. Nonetheless, natural selection operates on natural populations, and investigation of adaptive phenotypes in situ is warranted and can validate results from controlled laboratory experiments. To identify genomic effects on individual plant phenotypes in nature, environmental drivers must be concurrently measured, and characterized. Here, we designed and evaluated a sensor to meet these requirements for temperate woody plants. It was designed for use on a tree branch to measure the timing of budburst together with its key environmental drivers; temperature, and photoperiod. Specifically, we evaluated the sensor through independent corroboration with time-lapse photography and a suite of environmental sampling instruments. We also tested whether the presence of the device on a branch influenced the timing of budburst. Our results indicated the following: the temperatures measured by the budburst sensor's digital thermometer closely approximated the temperatures measured using a thermocouple touching plant tissue; the photoperiod detector measured ambient light with the same accuracy as did time lapse photography; the budburst sensor accurately detected the timing of budburst; and the sensor itself did not influence the budburst timing of Populus clones. Among other potential applications, future use of the sensor may provide plant phenotyping at the landscape level for integration with landscape genomics.

Using Populus as a lignocellulosic feedstock for bioethanol

Populus species along with species from the sister genus Salix will provide valuable feedstock resources for advanced second-generation biofuels. Their inherent fast growth characteristics can particularly be exploited for short rotation management, a time and energy saving cultivation alternative for lignocellulosic feedstock supply. Salicaceae possess inherent cell wall characteristics with favorable cellulose to lignin ratios for utilization as bioethanol crop. We review economically important traits relevant for intensively managed biofuel crop plantations, genomic and phenotypic resources available for Populus, breeding strategies for forest trees dedicated to bioenergy provision, and bioprocesses and downstream applications related to opportunities using Salicaceae as a renewable resource. Challenges need to be resolved for every single step of the conversion process chain, i.e., starting from tree domestication for improved performance as a bioenergy crop, bioconversion process, policy development for land use changes associated with advanced biofuels, and harvest and supply logistics associated with industrial-scale biorefinery plants using Populus as feedstock. Significant hurdles towards cost and energy efficiency, environmental friendliness, and yield maximization with regards to biomass pretreatment, saccharification, and fermentation of celluloses and the sustainability of biorefineries as a whole still need to be overcome.

Patterns and variability in seedling carbon assimilation: implications for tree recruitment under climate change

Predicting future forests' structure and functioning is a critical goal for ecologists, thus information on seedling recruitment will be crucial in determining the composition and structure of future forest ecosystems. In particular, seedlings' photosynthetic response to a changing environment will be a key component determining whether particular species establish enough individuals to maintain populations, as growth is a major determinant of survival. We quantified photosynthetic responses of sugar maple (Acer saccharum Marsh.), pignut hickory (Carya glabra Mill.), northern red oak (Quercus rubra L.) and eastern black oak (Quercus velutina Lam.) seedlings to environmental conditions including light habitat, temperature, soil moisture and vapor pressure deficit (VPD) using extensive in situ gas exchange measurements spanning an entire growing season. We estimated the parameters in a hierarchical Bayesian version of the Farquhar model of photosynthesis, additionally informed by soil moisture and VPD, and found that maximum Rubisco carboxylation (V(cmax)) and electron transport (J(max)) rates showed significant seasonal variation, but not the peaked patterns observed in studies of adult trees. Vapor pressure deficit and soil moisture limited J(max) and V(cmax) for all four species. Predictions indicate large declines in summer carbon assimilation rates under a 3 °C increase in mean annual temperature projected by climate models, while spring and fall assimilation rates may increase. Our model predicts decreases in summer assimilation rates in gap habitats with at least 90% probability, and with 20-99.9% probability in understory habitats depending on species. Predictions also show 70% probability of increases in fall and 52% probability in spring in understory habitats. All species were impacted, but our findings suggest that oak species may be favored in northeastern North America under projected increases in temperature due to superior assimilation rates under these conditions, though as growing seasons become longer, the effects of climate change on seedling photosynthesis may be complex.

Association genetics, geography and ecophysiology link stomatal patterning in Populus trichocarpa with carbon gain and disease resistance trade-offs

Stomata are essential for diffusive entry of gases to support photosynthesis, but may also expose internal leaf tissues to pathogens. To uncover trade-offs in range-wide adaptation relating to stomata, we investigated the underlying genetics of stomatal traits and linked variability in these traits with geoclimate, ecophysiology, condensed foliar tannins and pathogen susceptibility in black cottonwood (Populus trichocarpa). Upper (adaxial) and lower (abaxial) leaf stomatal traits were measured from 454 accessions collected throughout much of the species range. We calculated broad-sense heritability (H(2) ) of stomatal traits and, using SNP data from a 34K Populus SNP array, performed a genome-wide association studies (GWAS) to uncover genes underlying stomatal trait variation. H(2) values for stomatal traits were moderate (average H(2) = 0.33). GWAS identified genes associated primarily with adaxial stomata, including polarity genes (PHABULOSA), stomatal development genes (BRASSINOSTEROID-INSENSITIVE 2) and disease/wound-response genes (GLUTAMATE-CYSTEINE LIGASE). Stomatal traits correlated with latitude, gas exchange, condensed tannins and leaf rust (Melampsora) infection. Latitudinal trends of greater adaxial stomata numbers and guard cell pore size corresponded with higher stomatal conductance (gs ) and photosynthesis (Amax ), faster shoot elongation, lower foliar tannins and greater Melampsora susceptibility. This suggests an evolutionary trade-off related to differing selection pressures across the species range. In northern environments, more adaxial stomata and larger pore sizes reflect selection for rapid carbon gain and growth. By contrast, southern genotypes have fewer adaxial stomata, smaller pore sizes and higher levels of condensed tannins, possibly linked to greater pressure from natural leaf pathogens, which are less significant in northern ecosystems.

Landscape genomics of Populus trichocarpa: the role of hybridization, limited gene flow, and natural selection in shaping patterns of population structure

Populus trichocarpa is an ecologically important tree across western North America. We used a large population sample of 498 accessions over a wide geographical area genotyped with a 34K Populus SNP array to quantify geographical patterns of genetic variation in this species (landscape genomics). We present evidence that three processes contribute to the observed patterns: (1) introgression from the sister species P. balsamifera, (2) isolation by distance (IBD), and (3) natural selection. Introgression was detected only at the margins of the species' distribution. IBD was significant across the sampled area as a whole, but no evidence of restricted gene flow was detected in a core of drainages from southern British Columbia (BC). We identified a large number of FST outliers. Gene Ontology analyses revealed that FST outliers are overrepresented in genes involved in circadian rhythm and response to red/far-red light when the entire dataset is considered, whereas in southern BC heat response genes are overrepresented. We also identified strong correlations between geoclimate variables and allele frequencies at FST outlier loci that provide clues regarding the selective pressures acting at these loci.

Genome-wide association implicates numerous genes underlying ecological trait variation in natural populations of Populus trichocarpa

In order to uncover the genetic basis of phenotypic trait variation, we used 448 unrelated wild accessions of black cottonwood (Populus trichocarpa) from much of its range in western North America. Extensive data from large-scale trait phenotyping (with spatial and temporal replications within a common garden) and genotyping (with a 34 K Populus single nucleotide polymorphism (SNP) array) of all accessions were used for gene discovery in a genome-wide association study (GWAS). We performed GWAS with 40 biomass, ecophysiology and phenology traits and 29,355 filtered SNPs representing 3518 genes. The association analyses were carried out using a Unified Mixed Model accounting for population structure effects among accessions. We uncovered 410 significant SNPs using a Bonferroni-corrected threshold (P<1.7×10(-6)). Markers were found across 19 chromosomes, explained 1-13% of trait variation, and implicated 275 unique genes in trait associations. Phenology had the largest number of associated genes (240 genes), followed by biomass (53 genes) and ecophysiology traits (25 genes). The GWAS results propose numerous loci for further investigation. Many traits had significant associations with multiple genes, underscoring their genetic complexity. Genes were also identified with multiple trait associations within and/or across trait categories. In some cases, traits were genetically correlated while in others they were not.

Geographical and environmental gradients shape phenotypic trait variation and genetic structure in Populus trichocarpa

• Populus trichocarpa is widespread across western North America spanning extensive variation in photoperiod, growing season and climate. We investigated trait variation in P. trichocarpa using over 2000 trees from a common garden at Vancouver, Canada, representing replicate plantings of 461 genotypes originating from 136 provenance localities. • We measured 40 traits encompassing phenological events, biomass accumulation, growth rates, and leaf, isotope and gas exchange-based ecophysiology traits. With replicated plantings and 29,354 single nucleotide polymorphisms (SNPs) from 3518 genes, we estimated both broad-sense trait heritability (H(2)) and overall population genetic structure from principal component analysis. • Populus trichocarpa had high phenotypic variation and moderate/high H(2) for many traits. H(2) ranged from 0.3 to 0.9 in phenology, 0.3 to 0.8 in biomass and 0.1 to 0.8 in ecophysiology traits. Most traits correlated strongly with latitude, maximum daylength and temperature of tree origin, but not necessarily with elevation, precipitation or heat : moisture indices. Trait H(2) values reflected trait correlation strength with geoclimate variables. The population genetic structure had one significant principal component (PC1) which correlated with daylength and showed enrichment for genes relating to circadian rhythm and photoperiod. • Robust relationships between traits, population structure and geoclimate in P. trichocarpa reflect patterns which suggest that range-wide geographical and environment gradients have shaped its genotypic and phenotypic variability.

Extensive functional pleiotropy of REVOLUTA substantiated through forward genetics

In plants, genes may sustain extensive pleiotropic functional properties by individually affecting multiple, distinct traits. We discuss results from three genome-wide association studies of approximately 400 natural poplar (Populus trichocarpa) accessions phenotyped for 60 ecological/biomass, wood quality, and rust fungus resistance traits. Single-nucleotide polymorphisms (SNPs) in the poplar ortholog of the class III homeodomain-leucine zipper transcription factor gene REVOLUTA (PtREV) were significantly associated with three specific traits. Based on SNP associations with fungal resistance, leaf drop, and cellulose content, the PtREV gene contains three potential regulatory sites within noncoding regions at the gene's 3' end, where alternative splicing and messenger RNA processing actively occur. The polymorphisms in this region associated with leaf abscission and cellulose content are suggested to represent more recent variants, whereas the SNP associated with leaf rust resistance may be more ancient, consistent with REV's primary role in auxin signaling and its functional evolution in supporting fundamental processes of vascular plant development.

Tracking forest phenology and seasonal physiology using digital repeat photography: a critical assessment

Digital repeat photography is becoming widely used for near-surface remote sensing of vegetation. Canopy greenness, which has been used extensively for phenological applications, can be readily quantified from camera images. Important questions remain, however, as to whether the observed changes in canopy greenness are directly related to changes in leaf-level traits, changes in canopy structure, or some combination thereof. We investigated relationships between canopy greenness and various metrics of canopy structure and function, using five years (2008–2012) of automated digital imagery, ground observations of phenological transitions, leaf area index (LAI) measurements, and eddy covariance estimates of gross ecosystem photosynthesis from the Harvard Forest, a temperate deciduous forest in the northeastern United States. Additionally, we sampled canopy sunlit leaves on a weekly basis throughout the growing season of 2011. We measured physiological and morphological traits including leaf size, mass (wet/dry), nitrogen content, chlorophyll fluorescence, and spectral reflectance and characterized individual leaf color with flatbed scanner imagery. Our results show that observed spring and autumn phenological transition dates are well captured by information extracted from digital repeat photography. However, spring development of both LAI and the measured physiological and morphological traits are shown to lag behind spring increases in canopy greenness, which rises very quickly to its maximum value before leaves are even half their final size. Based on the hypothesis that changes in canopy greenness represent the aggregate effect of changes in both leaf-level properties (specifically, leaf color) and changes in canopy structure (specifically, LAI), we developed a two end-member mixing model. With just a single free parameter, the model was able to reproduce the observed seasonal trajectory of canopy greenness. This analysis shows that canopy greenness is relatively insensitive to changes in LAI at high LAI levels, which we further demonstrate by assessing the impact of an ice storm on both LAI and canopy greenness. Our study provides new insights into the mechanisms driving seasonal changes in canopy greenness retrieved from digital camera imagery. The nonlinear relationship between canopy greenness and canopy LAI has important implications both for phenological research applications and for assessing responses of vegetation to disturbances.