Why are triploid quaking aspen (Populus tremuloides) common?

PREMISE: Quaking aspen is a clonal tree species that has mixed ploidy, often with high relative abundance of both diploids and triploids but no haploids or tetraploids. Triploids typically have low fertility, leaving their occurrence apparently unlikely from an evolutionary perspective, unless they provide a "triploid bridge" to generating higher-fitness tetraploids-which are not observed in this species. This study focused on how triploidy can be maintained in quaking aspen. METHODS: A computational model was used to simulate gamete production, sexual reproduction, asexual reproduction, parent survival, and offspring survival in a population. All parameters were assumed to be cytotype-dependent and environment-independent. Sampling methods were used to identify parameter combinations consistent with observed cytotype frequencies. RESULTS: Many processes and parameter values were sufficient to yield a moderate frequency of triploids, and very few were necessary. The most plausible route involved higher triploid survival at the parent or offspring stage and limited unreduced gamete production by either diploid or triploid parents. Triploid fertility was helpful but not necessary. CONCLUSIONS: The coexistence of diploids and triploids in quaking aspen is statistically likely and promoted by the existence of commonly observed, long-lived triploid clones. However, other mechanisms not captured by the model related to environmental variation could also occur. Further empirical data or more complex but difficult-to-parameterize models are needed to gain further insight.

Negative allometry of leaf xylem conduit diameter and double-wall thickness: implications for implosion safety.

Xylem conduits have lignified walls to resist crushing pressures. The thicker the double-wall (T) relative to its diameter (D), the greater the implosion safety. Having safer conduits may incur higher costs and reduced flow, while having less resistant xylem may lead to catastrophic collapse under drought. Although recent studies have shown that conduit implosion commonly occurs in leaves, little is known about how leaf xylem scales T vs D to trade off safety, flow efficiency, mechanical support, and cost. We measured T and D in > 7000 conduits of 122 species to investigate how T vs D scaling varies across clades, habitats, growth forms, leaf, and vein sizes. As conduits become wider, their double-cell walls become proportionally thinner, resulting in a negative allometry between T and D. That is, narrower conduits, which are usually subjected to more negative pressures, are proportionally safer than wider ones. Higher implosion safety (i.e. higher T/D ratios) was found in asterids, arid habitats, shrubs, small leaves, and minor veins. Despite the strong allometry, implosion safety does not clearly trade off with other measured leaf functions, suggesting that implosion safety at whole-leaf level cannot be easily predicted solely by individual conduits' anatomy.

Climate lags and genetics determine phenology in quaking aspen (Populus tremuloides).

Spatiotemporal patterns of phenology may be affected by mosaics of environmental and genetic variation. Environmental drivers may have temporally lagged impacts, but patterns and mechanisms remain poorly known. We combine multiple genomic, remotely sensed, and physically modeled datasets to determine the spatiotemporal patterns and drivers of canopy phenology in quaking aspen, a widespread clonal dioecious tree species with diploid and triploid cytotypes. We show that over 391 km2 of southwestern Colorado: greenup date, greendown date, and growing season length vary by weeks and differ across sexes, cytotypes, and genotypes; phenology has high phenotypic plasticity and heritabilities of 31-61% (interquartile range); and snowmelt date, soil moisture, and air temperature predict phenology, at temporal lags of up to 3 yr. Our study shows that lagged environmental effects are needed to explain phenological variation and that the effect of cytotype on phenology is obscured by its correlation with topography. Phenological patterns are consistent with responses to multiyear accumulation of carbon deficit or hydraulic damage.

Vein network redundancy and mechanical resistance mitigate gas exchange losses under simulated herbivory in desert plants.

Herbivory can impact gas exchange, but the causes of interspecific variation in response remain poorly understood. We aimed to determine (1) what effects does experimental herbivory damage to leaf midveins have on leaf gas exchange and, (2) whether changes in leaf gas exchange after damage was predicted by leaf mechanical or venation traits. We hypothesized that herbivory-driven impacts on leaf gas exchange would be mediated by (1a/1b) venation networks, either by more vein resistance, or possibly trading off with other structural defenses; (2a/2b) or more reticulation (resilience, providing more alternate flow pathways after damage) or less reticulation (sectoriality, preventing spread of reduced functionality after damage). We simulated herbivory by damaging the midveins of four leaves from each of nine Sonoran Desert species. We then measured the percent change in photosynthesis (ΔAn%), transpiration (ΔEt%) and stomatal conductance (Δgsw%) between treated and control leaves. We assessed the relationship of each with leaf venation traits and other mechanical traits. ΔAn% varied between +10 % and -55%, similar to ΔEt% (+27%, -54%) and Δgsw% (+36%, -53%). There was no tradeoff between venation and other structural defenses. Increased damage resilience (reduced ΔAn%, ΔEt%, Δgsw%) was marginally associated with lower force-to-tear (P < 0.05), and higher minor vein density (P < 0.10) but not major vein density or reticulation. Leaf venation networks may thus partially mitigate the response of gas exchange to herbivory and other types of vein damage through either resistance or resilience.

Navigation between initial and desired community states using shortcuts.

Ecological management problems often involve navigating from an initial to a desired community state. We ask whether navigation without brute-force additions and deletions of species is possible via: adding/deleting a small number of individuals of a species, changing the environment, and waiting. Navigation can yield direct paths (single sequence of actions) or shortcut paths (multiple sequences of actions with lower cost than a direct path). We ask (1) when is non-brute-force navigation possible?; (2) do shortcuts exist and what are their properties?; and (3) what heuristics predict shortcut existence? Using a state diagram framework applied to several empirical datasets, we show that (1) non-brute-force navigation is only possible between some state pairs, (2) shortcuts exist between many state pairs; and (3) changes in abundance and richness are the strongest predictors of shortcut existence, independent of dataset and algorithm choices. State diagrams thus unveil hidden strategies for manipulating species coexistence and efficiently navigating between states.

Plant water use theory should incorporate hypotheses about extreme environments, population ecology, and community ecology.

Plant water use theory has largely been developed within a plant-performance paradigm that conceptualizes water use in terms of value for carbon gain and that sits within a neoclassical economic framework. This theory works very well in many contexts but does not consider other values of water to plants that could impact their fitness. Here, we survey a range of alternative hypotheses for drivers of water use and stomatal regulation. These hypotheses are organized around relevance to extreme environments, population ecology, and community ecology. Most of these hypotheses are not yet empirically tested and some are controversial (e.g. requiring more agency and behavior than is commonly believed possible for plants). Some hypotheses, especially those focused around using water to avoid thermal stress, using water to promote reproduction instead of growth, and using water to hoard it, may be useful to incorporate into theory or to implement in Earth System Models.

Does plant ecosystem thermoregulation occur? An extratropical assessment at different spatial and temporal scales.

To what degree plant ecosystems thermoregulate their canopy temperature (Tc ) is critical to assess ecosystems' metabolisms and resilience with climate change, but remains controversial, with opinions from no to moderate thermoregulation capability. With global datasets of Tc , air temperature (Ta ), and other environmental and biotic variables from FLUXNET and satellites, we tested the 'limited homeothermy' hypothesis (indicated by Tc & Ta regression slope < 1 or Tc < Ta around midday) across global extratropics, including temporal and spatial dimensions. Across daily to weekly and monthly timescales, over 80% of sites/ecosystems have slopes ≥1 or Tc > Ta around midday, rejecting the above hypothesis. For those sites unsupporting the hypothesis, their Tc -Ta difference (ΔT) exhibits considerable seasonality that shows negative, partial correlations with leaf area index, implying a certain degree of thermoregulation capability. Spatially, site-mean ΔT exhibits larger variations than the slope indicator, suggesting ΔT is a more sensitive indicator for detecting thermoregulatory differences across biomes. Furthermore, this large spatial-wide ΔT variation (0-6°C) is primarily explained by environmental variables (38%) and secondarily by biotic factors (15%). These results demonstrate diverse thermoregulation patterns across global extratropics, with most ecosystems negating the 'limited homeothermy' hypothesis, but their thermoregulation still occurs, implying that slope < 1 or Tc < Ta are not necessary conditions for plant thermoregulation.

The global spectrum of plant form and function: enhanced species-level trait dataset.

Here we provide the 'Global Spectrum of Plant Form and Function Dataset', containing species mean values for six vascular plant traits. Together, these traits -plant height, stem specific density, leaf area, leaf mass per area, leaf nitrogen content per dry mass, and diaspore (seed or spore) mass - define the primary axes of variation in plant form and function. The dataset is based on ca. 1 million trait records received via the TRY database (representing ca. 2,500 original publications) and additional unpublished data. It provides 92,159 species mean values for the six traits, covering 46,047 species. The data are complemented by higher-level taxonomic classification and six categorical traits (woodiness, growth form, succulence, adaptation to terrestrial or aquatic habitats, nutrition type and leaf type). Data quality management is based on a probabilistic approach combined with comprehensive validation against expert knowledge and external information. Intense data acquisition and thorough quality control produced the largest and, to our knowledge, most accurate compilation of empirically observed vascular plant species mean traits to date.

Gas exchange analysers exhibit large measurement error driven by internal thermal gradients.

Portable gas exchange analysers provide critical data for understanding plant-atmosphere carbon and water fluxes, and for parameterising Earth system models that forecast climate change effects and feedbacks. We characterised temperature measurement errors in the Li-Cor LI-6400XT and LI-6800, and estimated downstream errors in derived quantities, including stomatal conductance (gsw ) and leaf intercellular CO2 concentration (Ci ). The LI-6400XT exhibited air temperature errors (differences between reported air temperature and air temperature measured near the leaf) up to 7.2°C, leaf temperature errors up to 5.3°C, and relative errors in gsw and Ci that increased as temperatures departed from ambient. This caused errors in leaf-to-air temperature relationships, assimilation-temperature curves and CO2 response curves. Temperature dependencies of maximum Rubisco carboxylation rate (Vcmax ) and maximum RuBP regeneration rate (Jmax ) showed errors of 12% and 35%, respectively. These errors are likely to be idiosyncratic and may differ among machines and environmental conditions. The LI-6800 exhibited much smaller errors. Earth system model predictions may be erroneous, as much of their parametrisation data were measured on the LI-6400XT system, depending on the methods used. We make recommendations for minimising errors and correcting data in the LI-6400XT. We also recommend transitioning to the LI-6800 for future data collection.

Remote sensing of cytotype and its consequences for canopy damage in quaking aspen.

Mapping geographic mosaics of genetic variation and their consequences via genotype x environment interactions at large extents and high resolution has been limited by the scalability of DNA sequencing. Here, we address this challenge for cytotype (chromosome copy number) variation in quaking aspen, a drought-impacted foundation tree species. We integrate airborne imaging spectroscopy data with ground-based DNA sequencing data and canopy damage data in 391 km2 of southwestern Colorado. We show that (1) aspen cover and cytotype can be remotely sensed at 1 m spatial resolution, (2) the geographic mosaic of cytotypes is heterogeneous and interdigitated, (3) triploids have higher leaf nitrogen, canopy water content, and carbon isotope shifts (δ13 C) than diploids, and (4) canopy damage varies among cytotypes and depends on interactions with topography, canopy height, and trait variables. Triploids are at higher risk in hotter and drier conditions.

Process-explicit models reveal pathway to extinction for woolly mammoth using pattern-oriented validation.

Pathways to extinction start long before the death of the last individual. However, causes of early stage population declines and the susceptibility of small residual populations to extirpation are typically studied in isolation. Using validated process-explicit models, we disentangle the ecological mechanisms and threats that were integral in the initial decline and later extinction of the woolly mammoth. We show that reconciling ancient DNA data on woolly mammoth population decline with fossil evidence of location and timing of extinction requires process-explicit models with specific demographic and niche constraints, and a constrained synergy of climatic change and human impacts. Validated models needed humans to hasten climate-driven population declines by many millennia, and to allow woolly mammoths to persist in mainland Arctic refugia until the mid-Holocene. Our results show that the role of humans in the extinction dynamics of woolly mammoth began well before the Holocene, exerting lasting effects on the spatial pattern and timing of its range-wide extinction.

Cytotype and genotype predict mortality and recruitment in Colorado quaking aspen (Populus tremuloides).

Species responses to climate change depend on environment, genetics, and interactions among these factors. Intraspecific cytotype (ploidy level) variation is a common type of genetic variation in many species. However, the importance of intraspecific cytotype variation in determining demography across environments is poorly known. We studied quaking aspen (Populus tremuloides), which occurs in diploid and triploid cytotypes. This widespread tree species is experiencing contractions in its western range, which could potentially be linked to cytotype-dependent drought tolerance. We found that interactions between cytotype and environment drive mortality and recruitment across 503 plots in Colorado. Triploids were more vulnerable to mortality relative to diploids and had reduced recruitment on more drought-prone and disturbed plots relative to diploids. Furthermore, there was substantial genotype-dependent variation in demography. Thus, cytotype and genotype variation are associated with decline in this foundation species. Future assessment of demographic responses to climate change will benefit from knowledge of how genetic and environmental mosaics interact to determine species' ecophysiology and demography.

Automated and accurate segmentation of leaf venation networks via deep learning.

Leaf vein network geometry can predict levels of resource transport, defence and mechanical support that operate at different spatial scales. However, it is challenging to quantify network architecture across scales due to the difficulties both in segmenting networks from images and in extracting multiscale statistics from subsequent network graph representations. Here we developed deep learning algorithms using convolutional neural networks (CNNs) to automatically segment leaf vein networks. Thirty-eight CNNs were trained on subsets of manually defined ground-truth regions from >700 leaves representing 50 southeast Asian plant families. Ensembles of six independently trained CNNs were used to segment networks from larger leaf regions (c. 100 mm2 ). Segmented networks were analysed using hierarchical loop decomposition to extract a range of statistics describing scale transitions in vein and areole geometry. The CNN approach gave a precision-recall harmonic mean of 94.5% ± 6%, outperforming other current network extraction methods, and accurately described the widths, angles and connectivity of veins. Multiscale statistics then enabled the identification of previously undescribed variation in network architecture across species. We provide a LeafVeinCNN software package to enable multiscale quantification of leaf vein networks, facilitating the comparison across species and the exploration of the functional significance of different leaf vein architectures.

Linking functional traits to multiscale statistics of leaf venation networks.

Leaf venation networks evolved along several functional axes, including resource transport, damage resistance, mechanical strength, and construction cost. Because functions may depend on architectural features at different scales, network architecture may vary across spatial scales to satisfy functional tradeoffs. We develop a framework for quantifying network architecture with multiscale statistics describing elongation ratios, circularity ratios, vein density, and minimum spanning tree ratios. We quantify vein networks for leaves of 260 southeast Asian tree species in samples of up to 2 cm2 , pairing multiscale statistics with traits representing axes of resource transport, damage resistance, mechanical strength, and cost. We show that these multiscale statistics clearly differentiate species' architecture and delineate a phenotype space that shifts at larger scales; functional linkages vary with scale and are weak, with vein density, minimum spanning tree ratio, and circularity ratio linked to mechanical strength (measured by force to punch) and elongation ratio and circularity ratio linked to damage resistance (measured by tannins); and phylogenetic conservatism of network architecture is low but scale-dependent. This work provides tools to quantify the function and evolution of venation networks. Future studies including primary and secondary veins may uncover additional insights.

High water use in desert plants exposed to extreme heat.

Many plant water use models predict leaves maximize carbon assimilation while minimizing water loss via transpiration. Alternate scenarios may occur at high temperature, including heat avoidance, where leaves increase water loss to evaporatively cool regardless of carbon uptake; or heat failure, where leaves non-adaptively lose water also regardless of carbon uptake. We hypothesized that these alternative scenarios are common in species exposed to hot environments, with heat avoidance more common in species with high construction cost leaves. Diurnal measurements of leaf temperature and gas exchange for 11 Sonoran Desert species revealed that 37% of these species increased transpiration in the absence of increased carbon uptake. High leaf mass per area partially predicted this behaviour (r2  = 0.39). These data are consistent with heat avoidance and heat failure, but failure is less likely given the ecological dominance of the focal species. These behaviours are not yet captured in any extant plant water use model.

Community Assembly and Climate Mismatch in Late Quaternary Eastern North American Pollen Assemblages.

Plant community response to climate change ranges from synchronous tracking to strong mismatch. Explaining this variation in climate change response is critical for accurate global change modeling. Here we quantify how closely assemblages track changes in climate (match/mismatch) and how broadly climate niches are spread within assemblages (narrow/broad ecological tolerance, or "filtering") using data for the past 21,000 years for 531 eastern North American fossil pollen assemblages. Although climate matching has been strong over the last 21 millennia, mismatch increased in 30% of assemblages during the rapid climate shifts between 14.5 and 10 ka. Assemblage matching rebounded toward the present day in 10%-20% of assemblages. Climate-assemblage mismatch was greater in tree-dominated and high-latitude assemblages, consistent with persisting populations, slower dispersal rates, and glacial retreat. In contrast, climate matching was greater for assemblages comprising taxa with higher median seed mass. More than half of the assemblages were climatically filtered at any given time, with peak filtering occurring at 8.5 ka for nearly 80% of assemblages. Thus, vegetation assemblages have highly variable rates of climate mismatch and filtering over millennial scales. These climate responses can be partially predicted by species' traits and life histories. These findings help constrain predictions for plant community response to contemporary climate change.

Leaf venation networks of Bornean trees: images and hand-traced segmentations.

The data set contains images of leaf venation networks obtained from tree species in Malaysian Borneo. The data set contains 726 leaves from 295 species comprising 50 families, sampled from eight forest plots in Sabah. Image extents are approximately 1 × 1 cm, or 50 megapixels. All images contain a region of interest in which all veins have been hand traced. The complete data set includes over 30 billion pixels, of which more than 600 million have been validated by hand tracing. These images are suitable for morphological characterization of these species, as well as for training of machine-learning algorithms that segment biological networks from images. Data are made available under the Open Data Commons Attribution License. You are free to copy, distribute, and use the database; to produce works from the database; and to modify, transform, and build upon the database. You must attribute any public use of the database, or works produced from the database, in the manner specified in the license. For any use or redistribution of the database, or works produced from it, you must make clear to others the license of the database and keep intact any notices on the original database.

Tropical forest leaves may darken in response to climate change.

Tropical forest leaf albedo (reflectance) greatly impacts how much energy the planet absorbs; however; little is known about how it might be impacted by climate change. Here, we measure leaf traits and leaf albedo at ten 1-ha plots along a 3,200-m elevation gradient in Peru. Leaf mass per area (LMA) decreased with warmer temperatures along the elevation gradient; the distribution of LMA was positively skewed at all sites indicating a shift in LMA towards a warmer climate and future reduced tropical LMA. Reduced LMA was significantly (P < 0.0001) correlated with reduced leaf near-infrared (NIR) albedo; community-weighted mean NIR albedo significantly (P < 0.01) decreased as temperature increased. A potential future 2 °C increase in tropical temperatures could reduce lowland tropical leaf LMA by 6-7 g m-2 (5-6%) and reduce leaf NIR albedo by 0.0015-0.002 units. Reduced NIR albedo means that leaves are darker and absorb more of the Sun's energy. Climate simulations indicate this increased absorbed energy will warm tropical forests more at high CO2 conditions with proportionately more energy going towards heating and less towards evapotranspiration and cloud formation.

Global trait-environment relationships of plant communities.

Plant functional traits directly affect ecosystem functions. At the species level, trait combinations depend on trade-offs representing different ecological strategies, but at the community level trait combinations are expected to be decoupled from these trade-offs because different strategies can facilitate co-existence within communities. A key question is to what extent community-level trait composition is globally filtered and how well it is related to global versus local environmental drivers. Here, we perform a global, plot-level analysis of trait-environment relationships, using a database with more than 1.1 million vegetation plots and 26,632 plant species with trait information. Although we found a strong filtering of 17 functional traits, similar climate and soil conditions support communities differing greatly in mean trait values. The two main community trait axes that capture half of the global trait variation (plant stature and resource acquisitiveness) reflect the trade-offs at the species level but are weakly associated with climate and soil conditions at the global scale. Similarly, within-plot trait variation does not vary systematically with macro-environment. Our results indicate that, at fine spatial grain, macro-environmental drivers are much less important for functional trait composition than has been assumed from floristic analyses restricted to co-occurrence in large grid cells. Instead, trait combinations seem to be predominantly filtered by local-scale factors such as disturbance, fine-scale soil conditions, niche partitioning and biotic interactions.