Evolution of rarity and phylogeny determine above- and belowground biomass in plant-plant interactions.

Rare species are often considered inferior competitors due to occupancy of small ranges, specific habitats, and small local populations. However, the phylogenetic relatedness and rarity level (level 1-7 and common) of interacting species in plant-plant interactions are not often considered when predicting the response of rare plants in a biotic context. We used a common garden of 25 species of Tasmanian Eucalyptus, to differentiate non-additive patterns in the biomass of rare versus common species when grown in mixtures varying in phylogenetic relatedness and rarity. We demonstrate that rare species maintain progressively positive non-additive responses in biomass when interacting with phylogenetically intermediate, less rare and common species. This trend is not reflected in common species that out-performed in monocultures compared to mixtures. These results offer predictability as to how rare species' productivity will respond within various plant-plant interactions. However, species-specific interactions, such as those involving E. globulus, yielded a 97% increase in biomass compared to other species-specific interaction outcomes. These results are important because they suggest that plant rarity may also be shaped by biotic interactions, in addition to the known environmental and population factors normally used to describe rarity. Rare species may utilize potentially facilitative interactions with phylogenetically intermediate and common species to escape the effects of limiting similarity. Biotically mediated increases in rare plant biomass may have subsequent effects on the competitive ability and geographic occurrence of rare species, allowing rare species to persist at low abundance across plant communities. Through the consideration of species rarity and evolutionary history, we can more accurately predict plant-plant interaction dynamics to preserve unique ecosystem functions and fundamentally challenge what it means to be "rare".

Above- and belowground fungal biodiversity of Populus trees on a continental scale.

Understanding drivers of terrestrial fungal communities over large scales is an important challenge for predicting the fate of ecosystems under climate change and providing critical ecological context for bioengineering plant-microbe interactions in model systems. We conducted an extensive molecular and microscopy field study across the contiguous United States measuring natural variation in the Populus fungal microbiome among tree species, plant niche compartments and key symbionts. Our results show clear biodiversity hotspots and regional endemism of Populus-associated fungal communities explained by a combination of climate, soil and geographic factors. Modelling climate change impacts showed a deterioration of Populus mycorrhizal associations and an increase in potentially pathogenic foliar endophyte diversity and prevalence. Geographic differences among these symbiont groups in their sensitivity to environmental change are likely to influence broader forest health and ecosystem function. This dataset provides an above- and belowground atlas of Populus fungal biodiversity at a continental scale.

Stacked distribution models predict climate-driven loss of variation in leaf phenology at continental scales.

Climate change is having profound effects on species distributions and is likely altering the distribution of genetic variation across landscapes. Maintaining population genetic diversity is essential for the survival of species facing rapid environmental change, and variation loss will further ecological and evolutionary change. We used trait values of spring foliar leaf-out phenology of 400 genotypes from three geographically isolated populations of Populus angustifolia grown under common conditions, in concert with stacked species distribution modeling, to ask: (a) How will climate change alter phenological variation across the P. angustifolia species-range, and within populations; and (b) will the distribution of phenological variation among and within populations converge (become more similar) in future climatic conditions? Models predicted a net loss of phenological variation in future climate scenarios on 20-25% of the landscape across the species' range, with the trailing edge population losing variation on as much as 47% of the landscape. Our models also predicted that population's phenological trait distributions will become more similar over time. This stacked distribution model approach allows for the identification of areas expected to experience the greatest loss of genetically based functional trait variation and areas that may be priorities to conserve as future genetic climate refugia.

Aggregate population-level models informed by genetics predict more suitable habitat than traditional species-level model across the range of a widespread riparian tree.

Identifying and predicting how species ranges will shift in response to climate change is paramount for conservation and restoration. Ecological niche models are the most common method used to estimate potential distributions of species; however, they traditionally omit knowledge of intraspecific variation that can allow populations to respond uniquely to change. Here, we aim to test how population X environment relationships influence predicted suitable geographic distributions by comparing aggregated population-level models with species-level model predictions of suitable habitat within population ranges and across the species' range. We also test the effect of two variable selection methods on these predictions-both addressing the possibility of local adaptation: Models were built with (a) a common set, and number, of predictors and, (b) a unique combination and number of predictors specific to each group's training extent. Our study addresses the overarching hypothesis that populations have unique environmental niches, and specifically that (1) species-level models predict more suitable habitat within the ranges of genetic populations than individual models built from those groups, particularly when compared models are built with the same set of environmental predictors; and (2) aggregated genetic population models predict more suitable habitat across the species' range than the species-level model, an = d this difference will increase when models are trained with individualized predictors. We found the species models predicted more habitat within population ranges for two of three genetic groups regardless of variable selection, and that aggregated population models predicted more habitat than species' models, but that individualized predictors increased this difference. Our study emphasizes the extent to which changes to model predictions depend on the inclusion of genetic information and on the type and selection of predictors. Results from these modeling decisions can have broad implications for predicting population-level ecological and evolutionary responses to climate change.

Natural soil microbiome variation affects spring foliar phenology with consequences for plant productivity and climate-driven range shifts.

Identifying the potential for natural soil microbial communities to predictably affect complex plant traits is an important frontier in climate change research. Plant phenology varies with environmental and genetic factors, but few studies have examined whether the soil microbiome interacts with plant population differentiation to affect phenology and ecosystem function. We compared soil microbial variation in a widespread tree species (Populus angustifolia) with different soil inoculum treatments in a common garden environment to test how the soil microbiome affects spring foliar phenology and subsequent biomass growth. We hypothesized and show that soil bacterial and fungal communities vary with tree conditioning from different populations and elevations, that this soil community variation influences patterns of foliar phenology and plant growth across populations and elevation gradients, and that transferring lower elevation plant genotypes to higher elevation soil communities delayed foliar phenology, thereby shortening the growing season and reducing annual biomass production. Our findings show the importance of plant-soil interactions that help shape the timing of tree foliar phenology and productivity. These geographic patterns in plant population × microbiome interactions also broaden our understanding of how soil communities impact plant phenotypic variation across key climate change gradients, with consequences for ecosystem functioning.

Climate-driven divergence in plant-microbiome interactions generates range-wide variation in bud break phenology.

Soil microbiomes are rapidly becoming known as an important driver of plant phenotypic variation and may mediate plant responses to environmental factors. However, integrating spatial scales relevant to climate change with plant intraspecific genetic variation and soil microbial ecology is difficult, making studies of broad inference rare. Here we hypothesize and show: 1) the degree to which tree genotypes condition their soil microbiomes varies by population across the geographic distribution of a widespread riparian tree, Populus angustifolia; 2) geographic dissimilarity in soil microbiomes among populations is influenced by both abiotic and biotic environmental variation; and 3) soil microbiomes that vary in response to abiotic and biotic factors can change plant foliar phenology. We show soil microbiomes respond to intraspecific variation at the tree genotype and population level, and geographic variation in soil characteristics and climate. Using a fully reciprocal plant population by soil location feedback experiment, we identified a climate-based soil microbiome effect that advanced and delayed bud break phenology by approximately 10 days. These results demonstrate a landscape-level feedback between tree populations and associated soil microbial communities and suggest soil microbes may play important roles in mediating and buffering bud break phenology with climate warming, with whole ecosystem implications.

Intraspecific trait variation across elevation predicts a widespread tree species' climate niche and range limits.

Global change is widely altering environmental conditions which makes accurately predicting species range limits across natural landscapes critical for conservation and management decisions. If climate pressures along elevation gradients influence the distribution of phenotypic and genetic variation of plant functional traits, then such trait variation may be informative of the selective mechanisms and adaptations that help define climatic niche limits. Using extensive field surveys along 16 elevation transects and a large common garden experiment, we tested whether functional trait variation could predict the climatic niche of a widespread tree species (Populus angustifolia) with a double quantile regression approach. We show that intraspecific variation in plant size, growth, and leaf morphology corresponds with the species' total climate range and certain climatic limits related to temperature and moisture extremes. Moreover, we find evidence of genetic clines and phenotypic plasticity at environmental boundaries, which we use to create geographic predictions of trait variation and maximum values due to climatic constraints across the western US. Overall, our findings show the utility of double quantile regressions for connecting species distributions and climate gradients through trait-based mechanisms. We highlight how new approaches like ours that incorporate genetic variation in functional traits and their response to climate gradients will lead to a better understanding of plant distributions as well as identifying populations anticipated to be maladapted to future environments.

Plant genetic variation drives geographic differences in atmosphere-plant-ecosystem feedbacks.

The objective of this study was to understand how genetic variation in a riparian species, Populus angustifolia, affects mass and energy exchange between the land and atmosphere across ~1,700 km of latitude of the western United States. To examine the potential for large-scale land-atmosphere feedbacks in hydrologic processes driven by geographic differences in plant population traits, we use a physical hydrology model, paired field, and greenhouse observations of plant traits, and stable isotope compositions of soil, stem, and leaf water of P. angustifolia populations. Populations show patterns of local adaptation in traits related to landscape hydrologic functioning-a 47% difference in stomatal density in greenhouse conditions and a 74% difference in stomatal ratio in the field. Trait and stable isotope differences reveal that populations use water differently which is related to historical landscape hydrologic functioning (evapotranspiration and streamflow). Overall, results suggest that populations from landscapes with different hydrologic histories will differ in their ability to maintain favorable water balance with changing atmospheric demands for water, with ecosystem consequences.

Climate-driven reduction of genetic variation in plant phenology alters soil communities and nutrient pools.

We examined the hypothesis that climate-driven evolution of plant traits will influence associated soil microbiomes and ecosystem function across the landscape. Using a foundation tree species, Populus angustifolia, observational and common garden approaches, and a base population genetic collection that spans 17 river systems in the western United States, from AZ to MT, we show that (a) as mean annual temperature (MAT) increases, genetic and phenotypic variation for bud break phenology decline; (b) soil microbiomes, soil nitrogen (N), and soil carbon (C) vary in response to MAT and conditioning by trees; and (c) with losses of genetic variation due to warming, population-level regulation of community and ecosystem functions strengthen. These results demonstrate a relationship between the potential evolutionary response of populations and subsequent shifts in ecosystem function along a large temperature gradient.

The ecological importance of intraspecific variation.

Human activity is causing wild populations to experience rapid trait change and local extirpation. The resulting effects on intraspecific variation could have substantial consequences for ecological processes and ecosystem services. Although researchers have long acknowledged that variation among species influences the surrounding environment, only recently has evidence accumulated for the ecological importance of variation within species. We conducted a meta-analysis comparing the ecological effects of variation within a species (intraspecific effects) with the effects of replacement or removal of that species (species effects). We evaluated direct and indirect ecological responses, including changes in abundance (or biomass), rates of ecological processes and changes in community composition. Our results show that intraspecific effects are often comparable to, and sometimes stronger than, species effects. Species effects tend to be larger for direct ecological responses (for example, through consumption), whereas intraspecific effects and species effects tend to be similar for indirect responses (for example, through trophic cascades). Intraspecific effects are especially strong when indirect interactions alter community composition. Our results summarize data from the first generation of studies examining the relative ecological effects of intraspecific variation. Our conclusions can help inform the design of future experiments and the formulation of strategies to quantify and conserve biodiversity.

Divergent plant-soil feedbacks could alter future elevation ranges and ecosystem dynamics.

Plant-soil feedbacks (PSF) are important interactions that may influence range dynamics in a changing world. What remains largely unknown is the generality of plant-soil biotic interactions across populations and the potential role of specific soil biota, both of which are key for understanding how PSF might change future communities and ecosystems. We combined landscape-level field observations and experimental soil treatments to test whether a dominant tree alters soil environments to impact its own performance and range shifts towards higher elevations. We show: (1) soil conditioning by trees varies with elevation, (2) soil biota relate to PSF, (3) under simulated conditions, biotic PSF constrain range shifts at lower elevations but allow for expansions at higher elevations, and (4) differences in soil conditioning predict feedback outcomes in specific range-shift scenarios. These results suggest that variable plant-soil biotic interactions may influence the migration and fragmentation of tree species, and that models incorporating soil parameters will more accurately predict future species distributions.

Trait variation along elevation gradients in a dominant woody shrub is population-specific and driven by plasticity.

Elevation gradients are frequently used as space-for-time substitutions to infer species' trait responses to climate change. However, studies rarely investigate whether trait responses to elevation are widespread or population-specific within a species, and the relative genetic and plastic contributions to such trait responses may not be well understood. Here, we examine plant trait variation in the dominant woody shrub, Rhododendron maximum, along elevation gradients in three populations in the South Central Appalachian Mountains, USA, in both field and common garden environments. We ask the following: (i) do plant traits vary along elevation? (ii) do trait responses to elevation differ across populations, and if so, why? and (iii) does genetic differentiation or phenotypic plasticity drive trait variation within and among populations? We found that internode length, shoot length, leaf dry mass, and leaf area varied along elevation, but that these responses were generally unique to one population, suggesting that trait responses to environmental gradients are population-specific. A common garden experiment identified no genetic basis to variation along elevation or among populations in any trait, suggesting that plasticity drives local and regional trait variation and may play a key role in the persistence of plant species such as R. maximum with contemporary climate change. Overall, our findings highlight the importance of examining multiple locations in future elevation studies and indicate that, for a given plant species, the magnitude of trait responses to global climate change may vary by location.

Risk assessment of imidacloprid use in forest settings on the aquatic macroinvertebrate community.

The isolated effects of a single insecticide can be difficult to assess in natural settings because of the presence of numerous pollutants in many watersheds. Imidacloprid use for suppressing hemlock woolly adelgid, Adelges tsugae (Annand) (Hemiptera: Adelgidae), in forests offers a rare opportunity to assess potential impacts on aquatic macroinvertebrates in relatively pristine landscapes. Aquatic macroinvertebrate communities were assessed in 9 streams in Great Smoky Mountains National Park (southern Appalachian Mountains, USA). The streams flow through hemlock conservation areas where imidacloprid soil drench treatments were applied for hemlock woolly adelgid suppression. Sites were located upstream and downstream of the imidacloprid treatments. Baseline species presence data (pre-imidacloprid treatment) were available from previous sample collections at downstream sites. Downstream and upstream sites did not vary in numerous community measures. Although comparisons of paired upstream and downstream sites showed differences in diversity in 7 streams, higher diversity was found more often in downstream sites. Macroinvertebrate functional feeding groups and life habits were similar between downstream and upstream sites. Downstream and baseline stream samples were similar. While some functional feeding group and life habit species richness categories varied, variations did not indicate poorer quality downstream communities. Imidacloprid treatments applied according to US Environmental Protection Agency federal restrictions did not result in negative effects to aquatic macroinvertebrate communities, which indicates that risks of imidacloprid use in forest settings are low. Environ Toxicol Chem 2017;36:3108-3119. © 2017 SETAC.

Phylogeny is a powerful tool for predicting plant biomass responses to nitrogen enrichment.

Increasing rates of anthropogenic nitrogen (N) enrichment to soils often lead to the dominance of nitrophilic plant species and reduce plant diversity in natural ecosystems. Yet, we lack a framework to predict which species will be winners or losers in soil N enrichment scenarios, a framework that current literature suggests should integrate plant phylogeny, functional tradeoffs, and nutrient co-limitation. Using a controlled fertilization experiment, we quantified biomass responses to N enrichment for 23 forest tree species within the genus Eucalyptus that are native to Tasmania, Australia. Based on previous work with these species' responses to global change factors and theory on the evolution of plant resource-use strategies, we hypothesized that (1) growth responses to N enrichment are phylogenetically structured, (2) species with more resource-acquisitive functional traits have greater growth responses to N enrichment, and (3) phosphorus (P) limits growth responses to N enrichment differentially across species, wherein P enrichment increases growth responses to N enrichment more in some species than others. We built a hierarchical Bayesian model estimating effects of functional traits (specific leaf area, specific stem density, and specific root length) and P fertilization on species' biomass responses to N, which we then compared between lineages to determine whether phylogeny explains variation in responses to N. In concordance with literature on N limitation, a majority of species responded strongly and positively to N enrichment. Mean responses ranged three-fold, from 6.21 (E. pulchella) to 16.87 (E. delegatensis) percent increases in biomass per g N·m-2 ·yr-1 added. We identified a strong difference in responses to N between two phylogenetic lineages in the Eucalyptus subgenus Symphyomyrtus, suggesting that shared ancestry explains variation in N limitation. However, our model indicated that after controlling for phylogenetic non-independence, eucalypt responses to N were not associated with functional traits (although post-hoc analyses show a phylogenetic pattern in specific root length similar to that of responses to N), nor were responses differentially limited by P. Overall, our model results suggest that phylogeny is a powerful predictor of winners and losers in anthropogenic N enrichment scenarios in Tasmanian eucalypts, which may have implications for other species.

Genetics-based interactions of foundation species affect community diversity, stability and network structure.

We examined the hypothesis that genetics-based interactions between strongly interacting foundation species, the tree Populus angustifolia and the aphid Pemphigus betae, affect arthropod community diversity, stability and species interaction networks of which little is known. In a 2-year experimental manipulation of the tree and its aphid herbivore four major findings emerged: (i) the interactions of these two species determined the composition of an arthropod community of 139 species; (ii) both tree genotype and aphid presence significantly predicted community diversity; (iii) the presence of aphids on genetically susceptible trees increased the stability of arthropod communities across years; and (iv) the experimental removal of aphids affected community network structure (network degree, modularity and tree genotype contribution to modularity). These findings demonstrate that the interactions of foundation species are genetically based, which in turn significantly contributes to community diversity, stability and species interaction networks. These experiments provide an important step in understanding the evolution of Darwin's 'entangled bank', a metaphor that characterizes the complexity and interconnectedness of communities in the wild.

Phylogeny Explains Variation in The Root Chemistry of Eucalyptus Species.

Plants are dependent on their root systems for survival, and thus are defended from belowground enemies by a range of strategies, including plant secondary metabolites (PSMs). These compounds vary among species, and an understanding of this variation may provide generality in predicting the susceptibility of forest trees to belowground enemies and the quality of their organic matter input to soil. Here, we investigated phylogenetic patterns in the root chemistry of species within the genus Eucalyptus. Given the known diversity of PSMs in eucalypt foliage, we hypothesized that (i) the range and concentrations of PSMs and carbohydrates in roots vary among Eucalyptus species, and (ii) that phylogenetic relationships explain a significant component of this variation. To test for interspecific variation in root chemistry and the influence of tree phylogeny, we grew 24 Eucalyptus species representing two subgenera (Eucalyptus and Symphyomyrtus) in a common garden for two years. Fine root samples were collected from each species and analyzed for total phenolics, condensed tannins, carbohydrates, terpenes, and formylated phloroglucinol compounds. Compounds displaying significant interspecific variation were mapped onto a molecular phylogeny and tested for phylogenetic signal. Although all targeted groups of compounds were present, we found that phenolics dominated root defenses and that all phenolic traits displayed significant interspecific variation. Further, these compounds displayed a significant phylogenetic signal. Overall, our results suggest that within these representatives of genus Eucalyptus, more closely related species have more similar root chemistry, which may influence their susceptibility to belowground enemies and soil organic matter accrual.

Evolutionary history and novel biotic interactions determine plant responses to elevated CO2 and nitrogen fertilization.

A major frontier in global change research is predicting how multiple agents of global change will alter plant productivity, a critical component of the carbon cycle. Recent research has shown that plant responses to climate change are phylogenetically conserved such that species within some lineages are more productive than those within other lineages in changing environments. However, it remains unclear how phylogenetic patterns in plant responses to changing abiotic conditions may be altered by another agent of global change, the introduction of non-native species. Using a system of 28 native Tasmanian Eucalyptus species belonging to two subgenera, Symphyomyrtus and Eucalyptus, we hypothesized that productivity responses to abiotic agents of global change (elevated CO2 and increased soil N) are unique to lineages, but that novel interactions with a non-native species mediate these responses. We tested this hypothesis by examining productivity of 1) native species monocultures and 2) mixtures of native species with an introduced hardwood plantation species, Eucalyptus nitens, to experimentally manipulated soil N and atmospheric CO2. Consistent with past research, we found that N limits productivity overall, especially in elevated CO2 conditions. However, monocultures of species within the Symphyomyrtus subgenus showed the strongest response to N (gained 127% more total biomass) in elevated CO2 conditions, whereas those within the Eucalyptus subgenus did not respond to N. Root:shoot ratio (an indicator of resource use) was on average greater in species pairs containing Symphyomyrtus species, suggesting that functional traits important for resource uptake are phylogenetically conserved and explaining the phylogenetic pattern in plant response to changing environmental conditions. Yet, native species mixtures with E. nitens exhibited responses to CO2 and N that differed from those of monocultures, supporting our hypothesis and highlighting that both plant evolutionary history and introduced species will shape community productivity in a changing world.

Shifts in species interactions due to the evolution of functional differences between endemics and non-endemics: an endemic syndrome hypothesis.

Species ranges have been shifting since the Pleistocene, whereby fragmentation, isolation, and the subsequent reduction in gene flow have resulted in local adaptation of novel genotypes and the repeated evolution of endemic species. While there is a wide body of literature focused on understanding endemic species, very few studies empirically test whether or not the evolution of endemics results in unique function or ecological differences relative to their widespread congeners; in particular while controlling for environmental variation. Using a common garden composed of 15 Eucalyptus species within the subgenus Symphyomyrtus (9 endemic to Tasmania, 6 non-endemic), here we hypothesize and show that endemic species are functionally and ecologically different from non-endemics. Compared to non-endemics, endemic Eucalyptus species have a unique suite of functional plant traits that have extended effects on herbivores. We found that while endemics occupy many diverse habitats, they share similar functional traits potentially resulting in an endemic syndrome of traits. This study provides one of the first empirical datasets analyzing the functional differences between endemics and non-endemics in a common garden setting, and establishes a foundation for additional studies of endemic/non-endemic dynamics that will be essential for understanding global biodiversity in the midst of rapid species extinctions and range shifts as a consequence of global change.

Climate relicts and their associated communities as natural ecology and evolution laboratories.

Climate relicts, marginal populations that have become isolated via climate-driven range shifts, preserve ecological and evolutionary histories that can span millennia. Studies point to climate relicts as 'natural laboratories' for investigating how long-term environmental change impacts species and populations. However, we propose that such research should be expanded to reveal how climate change affects 'interacting' species in ways that reshape community composition and evolution. Biotic interactions and their community and ecosystem effects are often genetically based and driven by associations with foundation species. We discuss evolution in climate relicts within the context of the emerging fields of community and ecosystem genetics, exploring the idea that foundation relicts are also natural community and ecosystem laboratories and windows to future landscapes.

Evolutionary history determines how plant productivity responds to phylogenetic diversity and species richness.

The relationship between biodiversity and ecosystem function has received a great deal of attention in ecological research and recent results, from re-analyses, suggest that ecosystem function improves with increases in phylogenetic diversity. However, many of these results have been generalized across a range of different species and clades, and plants with different evolutionary histories could display different relationships between biodiversity and ecosystem function. To experimentally test this hypothesis, we manipulated species richness and phylogenetic diversity using 26 species from two subgenera of the genus Eucalyptus (subgenus Eucalyptus and subgenus Symphyomyrtus). We found that plant biomass (a measurement of ecosystem function) sometimes, but not always, responded to increases in species richness and phylogenetic diversity. Specifically, Symphyomyrtus plants showed a positive response while no comparable effect was observed for Eucalyptus plants, showing that responses to biodiversity can vary across different phylogenetic groups. Our results show that the impacts of evolutionary history may complicate the relationship between the diversity of plant communities and plant biomass.