Genomic mechanisms of climate adaptation in polyploid bioenergy switchgrass.

Long-term climate change and periodic environmental extremes threaten food and fuel security1 and global crop productivity2-4. Although molecular and adaptive breeding strategies can buffer the effects of climatic stress and improve crop resilience5, these approaches require sufficient knowledge of the genes that underlie productivity and adaptation6-knowledge that has been limited to a small number of well-studied model systems. Here we present the assembly and annotation of the large and complex genome of the polyploid bioenergy crop switchgrass (Panicum virgatum). Analysis of biomass and survival among 732 resequenced genotypes, which were grown across 10 common gardens that span 1,800 km of latitude, jointly revealed extensive genomic evidence of climate adaptation. Climate-gene-biomass associations were abundant but varied considerably among deeply diverged gene pools. Furthermore, we found that gene flow accelerated climate adaptation during the postglacial colonization of northern habitats through introgression of alleles from a pre-adapted northern gene pool. The polyploid nature of switchgrass also enhanced adaptive potential through the fractionation of gene function, as there was an increased level of heritable genetic diversity on the nondominant subgenome. In addition to investigating patterns of climate adaptation, the genome resources and gene-trait associations developed here provide breeders with the necessary tools to increase switchgrass yield for the sustainable production of bioenergy.

A generalist-specialist trade-off between switchgrass cytotypes impacts climate adaptation and geographic range.

Polyploidy results from whole-genome duplication and is a unique form of heritable variation with pronounced evolutionary implications. Different ploidy levels, or cytotypes, can exist within a single species, and such systems provide an opportunity to assess how ploidy variation alters phenotypic novelty, adaptability, and fitness, which can, in turn, drive the development of unique ecological niches that promote the coexistence of multiple cytotypes. Switchgrass, Panicum virgatum, is a widespread, perennial C4 grass in North America with multiple naturally occurring cytotypes, primarily tetraploids (4×) and octoploids (8×). Using a combination of genomic, quantitative genetic, landscape, and niche modeling approaches, we detect divergent levels of genetic admixture, evidence of niche differentiation, and differential environmental sensitivity between switchgrass cytotypes. Taken together, these findings support a generalist (8×)­specialist (4×) trade-off. Our results indicate that the 8× represent a unique combination of genetic variation that has allowed the expansion of switchgrass' ecological niche and thus putatively represents a valuable breeding resource.

Multiple constraints cause positive and negative feedbacks limiting grassland soil CO2 efflux under CO2 enrichment.

Terrestrial ecosystems are increasingly enriched with resources such as atmospheric CO2 that limit ecosystem processes. The consequences for ecosystem carbon cycling depend on the feedbacks from other limiting resources and plant community change, which remain poorly understood for soil CO2 efflux, JCO2, a primary carbon flux from the biosphere to the atmosphere. We applied a unique CO2 enrichment gradient (250 to 500 µL L-1) for eight years to grassland plant communities on soils from different landscape positions. We identified the trajectory of JCO2 responses and feedbacks from other resources, plant diversity [effective species richness, exp(H)], and community change (plant species turnover). We found linear increases in JCO2 on an alluvial sandy loam and a lowland clay soil, and an asymptotic increase on an upland silty clay soil. Structural equation modeling identified CO2 as the dominant limitation on JCO2 on the clay soil. In contrast with theory predicting limitation from a single limiting factor, the linear JCO2 response on the sandy loam was reinforced by positive feedbacks from aboveground net primary productivity and exp(H), while the asymptotic JCO2 response on the silty clay arose from a net negative feedback among exp(H), species turnover, and soil water potential. These findings support a multiple resource limitation view of the effects of global change drivers on grassland ecosystem carbon cycling and highlight a crucial role for positive or negative feedbacks between limiting resources and plant community structure. Incorporating these feedbacks will improve models of terrestrial carbon sequestration and ecosystem services.

Negative effects of nitrogen override positive effects of phosphorus on grassland legumes worldwide.

Anthropogenic nutrient enrichment is driving global biodiversity decline and modifying ecosystem functions. Theory suggests that plant functional types that fix atmospheric nitrogen have a competitive advantage in nitrogen-poor soils, but lose this advantage with increasing nitrogen supply. By contrast, the addition of phosphorus, potassium, and other nutrients may benefit such species in low-nutrient environments by enhancing their nitrogen-fixing capacity. We present a global-scale experiment confirming these predictions for nitrogen-fixing legumes (Fabaceae) across 45 grasslands on six continents. Nitrogen addition reduced legume cover, richness, and biomass, particularly in nitrogen-poor soils, while cover of non-nitrogen-fixing plants increased. The addition of phosphorous, potassium, and other nutrients enhanced legume abundance, but did not mitigate the negative effects of nitrogen addition. Increasing nitrogen supply thus has the potential to decrease the diversity and abundance of grassland legumes worldwide regardless of the availability of other nutrients, with consequences for biodiversity, food webs, ecosystem resilience, and genetic improvement of protein-rich agricultural plant species.

Legacies of precipitation influence primary production in Panicum virgatum.

Precipitation is a key driver of primary production worldwide, but primary production does not always track year-to-year variation in precipitation linearly. Instead, plant responses to changes in precipitation may exhibit time lags, or legacies of past precipitation. Legacies can be driven by multiple mechanisms, including persistent changes in plant physiological and morphological traits and changes to the physical environment, such as plant access to soil water. We used three precipitation manipulation experiments in central Texas, USA to evaluate the magnitude, duration, and potential mechanisms driving precipitation legacies on aboveground primary production of the perennial C4 grass, Panicum virgatum. Specifically, we performed a rainout shelter study, where eight genotypes grew under different precipitation regimes; a transplant study, where plants that had previously grown in a rainout shelter under different precipitation regimes were moved to a common environment; and a mesocosm study, where the effect of swapping precipitation regime was examined with a single genotype. Across these experiments, plants previously grown under wet conditions generally performed better than expected when exposed to drought. Panicum virgatum exhibited stronger productivity legacies of past wet years on current-year responses to drought than of past dry years on current-year responses to wet conditions. Additionally, previous year tiller counts, a proxy for meristem availability, were important in determining legacy effects on aboveground production. As climate changes and precipitation extremes-both dry and wet-become more common, these results suggest that populations of P. virgatum may become less resilient.

Nutrient identity modifies the destabilising effects of eutrophication in grasslands.

Nutrient enrichment can simultaneously increase and destabilise plant biomass production, with co-limitation by multiple nutrients potentially intensifying these effects. Here, we test how factorial additions of nitrogen (N), phosphorus (P) and potassium with essential nutrients (K+) affect the stability (mean/standard deviation) of aboveground biomass in 34 grasslands over 7 years. Destabilisation with fertilisation was prevalent but was driven by single nutrients, not synergistic nutrient interactions. On average, N-based treatments increased mean biomass production by 21-51% but increased its standard deviation by 40-68% and so consistently reduced stability. Adding P increased interannual variability and reduced stability without altering mean biomass, while K+ had no general effects. Declines in stability were largest in the most nutrient-limited grasslands, or where nutrients reduced species richness or intensified species synchrony. We show that nutrients can differentially impact the stability of biomass production, with N and P in particular disproportionately increasing its interannual variability.

Linking changes in species composition and biomass in a globally distributed grassland experiment.

Global change drivers, such as anthropogenic nutrient inputs, are increasing globally. Nutrient deposition simultaneously alters plant biodiversity, species composition and ecosystem processes like aboveground biomass production. These changes are underpinned by species extinction, colonisation and shifting relative abundance. Here, we use the Price equation to quantify and link the contributions of species that are lost, gained or that persist to change in aboveground biomass in 59 experimental grassland sites. Under ambient (control) conditions, compositional and biomass turnover was high, and losses (i.e. local extinctions) were balanced by gains (i.e. colonisation). Under fertilisation, the decline in species richness resulted from increased species loss and decreases in species gained. Biomass increase under fertilisation resulted mostly from species that persist and to a lesser extent from species gained. Drivers of ecological change can interact relatively independently with diversity, composition and ecosystem processes and functions such as aboveground biomass due to the individual contributions of species lost, gained or persisting.

Soil carbon stocks in temperate grasslands differ strongly across sites but are insensitive to decade-long fertilization.

Enhancing soil carbon (C) storage has the potential to offset human-caused increases in atmospheric CO2 . Rising CO2 has occurred concurrently with increasing supply rates of biologically limiting nutrients such as nitrogen (N) and phosphorus (P). However, it is unclear how increased supplies of N and P will alter soil C sequestration, particularly in grasslands, which make up nearly a third of non-agricultural land worldwide. Here, we leverage a globally distributed nutrient addition experiment (the Nutrient Network) to examine how a decade of N and P fertilization (alone and in combination) influenced soil C and N stocks at nine grassland sites spanning the continental United States. We measured changes in bulk soil C and N stocks and in three soil C fractions (light and heavy particulate organic matter, and mineral-associated organic matter fractions). Nutrient amendment had variable effects on soil C and N pools that ranged from strongly positive to strongly negative, while soil C and N pool sizes varied by more than an order of magnitude across sites. Piecewise SEM clarified that small increases in plant C inputs with fertilization did not translate to greater soil C storage. Nevertheless, peak season aboveground plant biomass (but not root biomass or production) was strongly positively related to soil C storage at seven of the nine sites, and across all nine sites, soil C covaried with moisture index and soil mineralogy, regardless of fertilization. Overall, we show that site factors such as moisture index, plant productivity, soil texture, and mineralogy were key predictors of cross-site soil C, while nutrient amendment had weaker and site-specific effects on C sequestration. This suggests that prioritizing the protection of highly productive temperate grasslands is critical for reducing future greenhouse gas losses arising from land use change.

The genetic basis for panicle trait variation in switchgrass (Panicum virgatum).

KEY MESSAGE: We investigate the genetic basis of panicle architecture in switchgrass in two mapping populations across a latitudinal gradient, and find many stable, repeatable genetic effects and limited genetic interactions with the environment. Grass species exhibit large diversity in panicle architecture influenced by genes, the environment, and their interaction. The genetic study of panicle architecture in perennial grasses is limited. In this study, we evaluate the genetic basis of panicle architecture including panicle length, primary branching number, and secondary branching number in an outcrossed switchgrass QTL population grown across ten field sites in the central USA through multi-environment mixed QTL analysis. We also evaluate genetic effects in a diversity panel of switchgrass grown at three of the ten field sites using genome-wide association (GWAS) and multivariate adaptive shrinkage. Furthermore, we search for candidate genes underlying panicle traits in both of these independent mapping populations. Overall, 18 QTL were detected in the QTL mapping population for the three panicle traits, and 146 unlinked genomic regions in the diversity panel affected one or more panicle trait. Twelve of the QTL exhibited consistent effects (i.e., no QTL by environment interactions or no QTL × E), and most (four of six) of the effects with QTL × E exhibited site-specific effects. Most (59.3%) significant partially linked diversity panel SNPs had significant effects in all panicle traits and all field sites and showed pervasive pleiotropy and limited environment interactions. Panicle QTL co-localized with significant SNPs found using GWAS, providing additional power to distinguish between true and false associations in the diversity panel.

Integrative modelling reveals mechanisms linking productivity and plant species richness.

How ecosystem productivity and species richness are interrelated is one of the most debated subjects in the history of ecology. Decades of intensive study have yet to discern the actual mechanisms behind observed global patterns. Here, by integrating the predictions from multiple theories into a single model and using data from 1,126 grassland plots spanning five continents, we detect the clear signals of numerous underlying mechanisms linking productivity and richness. We find that an integrative model has substantially higher explanatory power than traditional bivariate analyses. In addition, the specific results unveil several surprising findings that conflict with classical models. These include the isolation of a strong and consistent enhancement of productivity by richness, an effect in striking contrast with superficial data patterns. Also revealed is a consistent importance of competition across the full range of productivity values, in direct conflict with some (but not all) proposed models. The promotion of local richness by macroecological gradients in climatic favourability, generally seen as a competing hypothesis, is also found to be important in our analysis. The results demonstrate that an integrative modelling approach leads to a major advance in our ability to discern the underlying processes operating in ecological systems.

Addition of multiple limiting resources reduces grassland diversity.

Niche dimensionality provides a general theoretical explanation for biodiversity-more niches, defined by more limiting factors, allow for more ways that species can coexist. Because plant species compete for the same set of limiting resources, theory predicts that addition of a limiting resource eliminates potential trade-offs, reducing the number of species that can coexist. Multiple nutrient limitation of plant production is common and therefore fertilization may reduce diversity by reducing the number or dimensionality of belowground limiting factors. At the same time, nutrient addition, by increasing biomass, should ultimately shift competition from belowground nutrients towards a one-dimensional competitive trade-off for light. Here we show that plant species diversity decreased when a greater number of limiting nutrients were added across 45 grassland sites from a multi-continent experimental network. The number of added nutrients predicted diversity loss, even after controlling for effects of plant biomass, and even where biomass production was not nutrient-limited. We found that elevated resource supply reduced niche dimensionality and diversity and increased both productivity and compositional turnover. Our results point to the importance of understanding dimensionality in ecological systems that are undergoing diversity loss in response to multiple global change factors.

QTL × environment interactions underlie adaptive divergence in switchgrass across a large latitudinal gradient.

Local adaptation is the process by which natural selection drives adaptive phenotypic divergence across environmental gradients. Theory suggests that local adaptation results from genetic trade-offs at individual genetic loci, where adaptation to one set of environmental conditions results in a cost to fitness in alternative environments. However, the degree to which there are costs associated with local adaptation is poorly understood because most of these experiments rely on two-site reciprocal transplant experiments. Here, we quantify the benefits and costs of locally adaptive loci across 17° of latitude in a four-grandparent outbred mapping population in outcrossing switchgrass (Panicum virgatum L.), an emerging biofuel crop and dominant tallgrass species. We conducted quantitative trait locus (QTL) mapping across 10 sites, ranging from Texas to South Dakota. This analysis revealed that beneficial biomass (fitness) QTL generally incur minimal costs when transplanted to other field sites distributed over a large climatic gradient over the 2 y of our study. Therefore, locally advantageous alleles could potentially be combined across multiple loci through breeding to create high-yielding regionally adapted cultivars.

Soil properties as key predictors of global grassland production: Have we overlooked micronutrients?

Fertilisation experiments have demonstrated that nutrient availability is a key determinant of biomass production and carbon sequestration in grasslands. However, the influence of nutrients in explaining spatial variation in grassland biomass production has rarely been assessed. Using a global dataset comprising 72 sites on six continents, we investigated which of 16 soil factors that shape nutrient availability associate most strongly with variation in grassland aboveground biomass. Climate and N deposition were also considered. Based on theory-driven structural equation modelling, we found that soil micronutrients (particularly Zn and Fe) were important predictors of biomass and, together with soil physicochemical properties and C:N, they explained more unique variation (32%) than climate and N deposition (24%). However, the association between micronutrients and biomass was absent in grasslands limited by NP. These results highlight soil properties as key predictors of global grassland biomass production and point to serial co-limitation by NP and micronutrients.

Geographic variation in the genetic basis of resistance to leaf rust between locally adapted ecotypes of the biofuel crop switchgrass (Panicum virgatum).

Local adaptation is an important process in plant evolution, which can be impacted by differential pathogen pressures along environmental gradients. However, the degree to which pathogen resistance loci vary in effect across space and time is incompletely described. To understand how the genetic architecture of resistance varies across time and geographic space, we quantified rust (Puccinia spp.) severity in switchgrass (Panicum virgatum) plantings at eight locations across the central USA for 3 yr and conducted quantitative trait locus (QTL) mapping for rust progression. We mapped several variable QTLs, but two large-effect QTLs which we have named Prr1 and Prr2 were consistently associated with rust severity in multiple sites and years, particularly in northern sites. By contrast, there were numerous small-effect QTLs at southern sites, indicating a genotype-by-environment interaction in rust resistance loci. Interestingly, Prr1 and Prr2 had a strong epistatic interaction, which also varied in the strength and direction of effect across space. Our results suggest that abiotic factors covarying with latitude interact with the genetic loci underlying plant resistance to control rust infection severity. Furthermore, our results indicate that segregating genetic variation in epistatically interacting loci may play a key role in determining response to infection across geographic space.

Global impacts of fertilization and herbivore removal on soil net nitrogen mineralization are modulated by local climate and soil properties.

Soil nitrogen (N) availability is critical for grassland functioning. However, human activities have increased the supply of biologically limiting nutrients, and changed the density and identity of mammalian herbivores. These anthropogenic changes may alter net soil N mineralization (soil net Nmin ), that is, the net balance between N mineralization and immobilization, which could severely impact grassland structure and functioning. Yet, to date, little is known about how fertilization and herbivore removal individually, or jointly, affect soil net Nmin across a wide range of grasslands that vary in soil and climatic properties. Here we collected data from 22 grasslands on five continents, all part of a globally replicated experiment, to assess how fertilization and herbivore removal affected potential (laboratory-based) and realized (field-based) soil net Nmin . Herbivore removal in the absence of fertilization did not alter potential and realized soil net Nmin . However, fertilization alone and in combination with herbivore removal consistently increased potential soil net Nmin. Realized soil net Nmin , in contrast, significantly decreased in fertilized plots where herbivores were removed. Treatment effects on potential and realized soil net Nmin were contingent on site-specific soil and climatic properties. Fertilization effects on potential soil net Nmin were larger at sites with higher mean annual precipitation (MAP) and temperature of the wettest quarter (T.q.wet). Reciprocally, realized soil net Nmin declined most strongly with fertilization and herbivore removal at sites with lower MAP and higher T.q.wet. In summary, our findings show that anthropogenic nutrient enrichment, herbivore exclusion and alterations in future climatic conditions can negatively impact soil net Nmin across global grasslands under realistic field conditions. This is an important context-dependent knowledge for grassland management worldwide.

Eutrophication weakens stabilizing effects of diversity in natural grasslands.

Studies of experimental grassland communities have demonstrated that plant diversity can stabilize productivity through species asynchrony, in which decreases in the biomass of some species are compensated for by increases in others. However, it remains unknown whether these findings are relevant to natural ecosystems, especially those for which species diversity is threatened by anthropogenic global change. Here we analyse diversity-stability relationships from 41 grasslands on five continents and examine how these relationships are affected by chronic fertilization, one of the strongest drivers of species loss globally. Unmanipulated communities with more species had greater species asynchrony, resulting in more stable biomass production, generalizing a result from biodiversity experiments to real-world grasslands. However, fertilization weakened the positive effect of diversity on stability. Contrary to expectations, this was not due to species loss after eutrophication but rather to an increase in the temporal variation of productivity in combination with a decrease in species asynchrony in diverse communities. Our results demonstrate separate and synergistic effects of diversity and eutrophication on stability, emphasizing the need to understand how drivers of global change interactively affect the reliable provisioning of ecosystem services in real-world systems.

CO2 enrichment and soil type additively regulate grassland productivity.

Atmospheric CO2 enrichment usually increases the aboveground net primary productivity (ANPP) of grassland vegetation, but the magnitude of the ANPP-CO2 response differs among ecosystems. Soil properties affect ANPP via multiple mechanisms and vary over topographic to geographic gradients, but have received little attention as potential modifiers of the ANPP-CO2 response. We assessed the effects of three soil types, sandy loam, silty clay and clay, on the ANPP response of perennial C3 /C4 grassland communities to a subambient to elevated CO2 gradient over 10 yr in Texas, USA. We predicted an interactive, rather than additive, effect of CO2 and soil type on ANPP. Contrary to prediction, CO2 and soil additively influenced grassland ANPP. Increasing CO2 by 250 µl l-1 increased ANPP by 170 g m-2 across soil types. Increased clay content from 10% to 50% among soils reduced ANPP by 50 g m-2 . CO2 enrichment increased ANPP via a predominant direct effect, accompanied by a smaller indirect effect mediated by a successional shift to increased dominance of the C4 tallgrass Sorghastrum nutans. Our results indicate a large, positive influence of CO2 enrichment on grassland productivity that resulted from the direct physiological benefits of CO2 augmented by species succession, and was expressed similarly across soils of differing physical properties.

Spatial heterogeneity in species composition constrains plant community responses to herbivory and fertilisation.

Environmental change can result in substantial shifts in community composition. The associated immigration and extinction events are likely constrained by the spatial distribution of species. Still, studies on environmental change typically quantify biotic responses at single spatial (time series within a single plot) or temporal (spatial beta diversity at single time points) scales, ignoring their potential interdependence. Here, we use data from a global network of grassland experiments to determine how turnover responses to two major forms of environmental change - fertilisation and herbivore loss - are affected by species pool size and spatial compositional heterogeneity. Fertilisation led to higher rates of local extinction, whereas turnover in herbivore exclusion plots was driven by species replacement. Overall, sites with more spatially heterogeneous composition showed significantly higher rates of annual turnover, independent of species pool size and treatment. Taking into account spatial biodiversity aspects will therefore improve our understanding of consequences of global and anthropogenic change on community dynamics.

Bacterial community response to a preindustrial-to-future CO2 gradient is limited and soil specific in Texas Prairie grassland.

Rising atmospheric CO2 concentration directly stimulates plant productivity and affects nutrient dynamics in the soil. However, the influence of CO2 enrichment on soil bacterial communities remains elusive, likely due to their complex interactions with a wide range of plant and soil properties. Here, we investigated the bacterial community response to a decade long preindustrial-to-future CO2 gradient (250-500 ppm) among three contrasting soil types using 16S rRNA gene amplicon sequencing. In addition, we examined the effect of seasonal variation and plant species composition on bacterial communities. We found that Shannon index (H') and Faith's phylogenetic diversity (PD) did not change in response to the CO2 gradient (R2  = 0.01, p > 0.05). CO2 gradient and season also had a negligible effect on overall community structure, although silty clay soil communities were better structured on a CO2 gradient (p < 0.001) among three soils. Similarly, CO2 gradient had no significant effect on the relative abundance of different phyla. However, we observed soil-specific variation of CO2 effects in a few individual families. For example, the abundance of Pirellulaceae family decreased linearly with CO2 gradient, but only in sandy loam soils. Conversely, the abundance of Micromonosporaceae and Gaillaceae families increased with CO2 gradient in clay soils. Soil water content (SWC) and nutrient properties were the key environmental constraints shaping bacterial community structure, one manifestation of which was a decline in bacterial diversity with increasing SWC. Furthermore, the impact of plant species composition on community structure was secondary to the strong influence of soil properties. Taken together, our findings indicate that bacterial communities may be largely unresponsive to indirect effects of CO2 enrichment through plants. Instead, bacterial communities are strongly regulated by edaphic conditions, presumably because soil differences create distinct environmental niches for bacteria.

Flowering in grassland predicted by CO2 and resource effects on species aboveground biomass.

Continuing enrichment of atmospheric CO2 may change plant community composition, in part by altering the availability of other limiting resources including soil water, nutrients, or light. The combined effects of CO2 enrichment and altered resource availability on species flowering remain poorly understood. We quantified flowering culm and ramet production and biomass allocation to flowering culms/ramets for 10 years in C4 -dominated grassland communities on contrasting soils along a CO2 concentration gradient spanning pre-industrial to expected mid-21st century levels (250-500 µl/L). CO2 enrichment explained up to 77% of the variation in flowering culm count across soils for three of the five species, and was correlated with flowering culm count on at least one soil for four of five species. In contrast, allocation to flowering culms was only weakly correlated with CO2 enrichment for two species. Flowering culm counts were strongly correlated with species aboveground biomass (AGB; R2  = .34-.74), a measure of species abundance. CO2 enrichment also increased soil moisture and decreased light levels within the canopy but did not affect soil inorganic nitrogen availability. Structural equation models fit across the soils suggested species-specific controls on flowering in two general forms: (1) CO2 effects on flowering culm count mediated by canopy light level and relative species AGB (species AGB/total AGB) or by soil moisture effects on flowering culm count; (2) effects of canopy light level or soil inorganic nitrogen on flowering and/or relative species AGB, but with no significant CO2 effect. Understanding the heterogeneity in species responses to CO2 enrichment in plant communities across soils in edaphically variable landscapes is critical to predict CO2 effects on flowering and other plant fitness components, and species potential to adapt to future environmental changes.