The value of long-term ecological research for evolutionary insights.

Scientists must have an integrative understanding of ecology and evolution across spatial and temporal scales to predict how species will respond to global change. Although comprehensively investigating these processes in nature is challenging, the infrastructure and data from long-term ecological research networks can support cross-disciplinary investigations. We propose using these networks to advance our understanding of fundamental evolutionary processes and responses to global change. For ecologists, we outline how long-term ecological experiments can be expanded for evolutionary inquiry, and for evolutionary biologists, we illustrate how observed long-term ecological patterns may motivate new evolutionary questions. We advocate for collaborative, multi-site investigations and discuss barriers to conducting evolutionary work at network sites. Ultimately, these networks offer valuable information and opportunities to improve predictions of species' responses to global change.

Genomic prediction of regional-scale performance in switchgrass (Panicum virgatum) by accounting for genotype-by-environment variation and yield surrogate traits.

Switchgrass is a potential crop for bioenergy or carbon capture schemes, but further yield improvements through selective breeding are needed to encourage commercialization. To identify promising switchgrass germplasm for future breeding efforts, we conducted multi-site and multi-trait genomic prediction with a diversity panel of 630 genotypes from 4 switchgrass subpopulations (Gulf, Midwest, Coastal, and Texas), which were measured for spaced plant biomass yield across 10 sites. Our study focused on the use of genomic prediction to share information among traits and environments. Specifically, we evaluated the predictive ability of cross-validation (CV) schemes using only genetic data and the training set, (cross validation 1: CV1), a subset of the sites (cross validation 2: CV2), and/or with two yield surrogates (flowering time and fall plant height). We found that genotype-by-environment interactions were largely due to the north-south distribution of sites. The genetic correlations between yield surrogates and biomass yield were generally positive (mean height r=0.85; mean flowering time r=0.45) and did not vary due to subpopulation or growing region (North, Middle, South). Genomic prediction models had cross-validation predictive abilities of -0.02 for individuals using only genetic data (CV1) but 0.55, 0.69, 0.76, 0.81, and 0.84 for individuals with biomass performance data from one, two, three, four and five sites included in the training data (CV2), respectively. To simulate a resource-limited breeding program, we determined the predictive ability of models provided with: one site observation of flowering time (0.39), one site observation of flowering time and fall height (0.51), one site observation of fall height (0.52), one site observation of biomass (0.55), and five site observations of biomass yield (0.84). The ability to share information at a regional scale is very encouraging but further research is required to accurately translate spaced plant biomass to commercial-scale sward biomass performance.

Programmed cell death regulator BAP2 is required for IRE1-mediated unfolded protein response in Arabidopsis.

Environmental and physiological situations can challenge the balance between protein synthesis and folding capacity of the endoplasmic reticulum (ER) and cause ER stress, a potentially lethal condition. The unfolded protein response (UPR) restores ER homeostasis or actuates programmed cell death (PCD) when ER stress is unresolved. The cell fate determination mechanisms of the UPR are not well understood, especially in plants. Here, we integrate genetics and ER stress profiling with natural variation and quantitative trait locus analysis of 350 natural accessions of the model species Arabidopsis thaliana. Our analyses implicate a single nucleotide polymorphism to the loss of function of the general PCD regulator BON-ASSOCIATED PROTEIN2 (BAP2) in UPR outcomes. We establish that ER stress-induced BAP2 expression is antagonistically regulated by the UPR master regulator, inositol-requiring enzyme 1 (IRE1), and that BAP2 controls adaptive UPR amplitude in ER stress and ignites pro-death mechanisms in conditions of UPR insufficiency.

Genetically correlated leaf tensile and morphological traits are driven by growing season length in a widespread perennial grass.

PREMISE: Leaf tensile resistance, a leaf's ability to withstand pulling forces, is an important determinant of plant ecological strategies. One potential driver of leaf tensile resistance is growing season length. When growing seasons are long, strong leaves, which often require more time and resources to construct than weak leaves, may be more advantageous than when growing seasons are short. Growing season length and other ecological conditions may also impact the morphological traits that underlie leaf tensile resistance. METHODS: To understand variation in leaf tensile resistance, we measured size-dependent leaf strength and size-independent leaf toughness in diverse genotypes of the widespread perennial grass Panicum virgatum (switchgrass) in a common garden. We then used quantitative genetic approaches to estimate the heritability of leaf tensile resistance and whether there were genetic correlations between leaf tensile resistance and other morphological traits. RESULTS: Leaf tensile resistance was positively associated with aboveground biomass (a proxy for fitness). Moreover, both measures of leaf tensile resistance exhibited high heritability and were positively genetically correlated with leaf lamina thickness and leaf mass per area (LMA). Leaf tensile resistance also increased with the growing season length in the habitat of origin, and this effect was mediated by both LMA and leaf thickness. CONCLUSIONS: Differences in growing season length may promote selection for different leaf lifespans and may explain existing variation in leaf tensile resistance in P. virgatum. In addition, the high heritability of leaf tensile resistance suggests that P. virgatum will be able to respond to climate change as growing seasons lengthen.

Arabidopsis AN3 and OLIGOCELLULA genes link telomere maintenance mechanisms with cell division and expansion control.

Telomeres are conserved chromosomal structures necessary for continued cell division and proliferation. In addition to the classical telomerase pathway, multiple other genes including those involved in ribosome metabolism and chromatin modification contribute to telomere length maintenance. We previously reported that Arabidopsis thaliana ribosome biogenesis genes OLI2/NOP2A, OLI5/RPL5A and OLI7/RPL5B have critical roles in telomere length regulation. These three OLIGOCELLULA genes were also shown to function in cell proliferation and expansion control and to genetically interact with the transcriptional co-activator ANGUSTIFOLIA3 (AN3). Here we show that AN3-deficient plants progressively lose telomeric DNA in early homozygous mutant generations, but ultimately establish a new shorter telomere length setpoint by the fifth mutant generation with a telomere length similar to oli2/nop2a -deficient plants. Analysis of double an3 oli2 mutants indicates that the two genes are epistatic for telomere length control. Telomere shortening in an3 and oli mutants is not caused by telomerase inhibition; wild type levels of telomerase activity are detected in all analyzed mutants in vitro. Late generations of an3 and oli mutants are prone to stem cell damage in the root apical meristem, implying that genes regulating telomere length may have conserved functional roles in stem cell maintenance mechanisms. Multiple instances of anaphase fusions in late generations of oli5 and oli7 mutants were observed, highlighting an unexpected effect of ribosome biogenesis factors on chromosome integrity. Overall, our data implicate AN3 transcription coactivator and OLIGOCELLULA proteins in the establishment of telomere length set point in plants and further suggest that multiple regulators with pleiotropic functions can connect telomere biology with cell proliferation and cell expansion pathways.

Physiological Responses of C4 Perennial Bioenergy Grasses to Climate Change: Causes, Consequences, and Constraints.

C4 perennial bioenergy grasses are an economically and ecologically important group whose responses to climate change will be important to the future bioeconomy. These grasses are highly productive and frequently possess large geographic ranges and broad environmental tolerances, which may contribute to the evolution of ecotypes that differ in physiological acclimation capacity and the evolution of distinct functional strategies. C4 perennial bioenergy grasses are predicted to thrive under climate change-C4 photosynthesis likely evolved to enhance photosynthetic efficiency under stressful conditions of low [CO2], high temperature, and drought-although few studies have examined how these species will respond to combined stresses or to extremes of temperature and precipitation. Important targets for C4 perennial bioenergy production in a changing world, such as sustainability and resilience, can benefit from combining knowledge of C4 physiology with recent advances in crop improvement, especially genomic selection.

Seed traits and recruitment interact with habitats to generate patterns of local adaptation in a perennial grass.

A fundamental challenge in the field of ecology involves understanding the adaptive traits and life history stages regulating the population dynamics of species across diverse habitats. Seed traits and early seedling vigor are thought to be key functional traits in plants, with important consequences for recruitment, establishment, and population persistence. However, little is known about how diverse seed traits interact with seed and microsite availability to impact plant populations. Here, we performed a factorial experiment involving seed addition and surface soil disturbance to explore the combined effects of seed and site availability using genotypes characterized by varying seed mass and dormancy traits. Additionally, we included hybrids that exhibited recombined seed trait relationships compared with natural genotypes, allowing us to assess the impact of specific seed traits on establishment across different sites. We detected a significant three-way interaction between seed addition, site conditions, and soil surface disturbance, influencing both seedling establishment and adult recruitment in Panicum hallii, a perennial grass found in coastal mesic (lowland) and inland xeric (upland) habitats. This establishment/recruitment pattern suggests that mesic and xeric establishment at foreign sites is constrained by the interplay of seed and site limitations. Notably, soil surface disturbance facilitated establishment and recruitment of the xeric genotype while limiting the mesic genotype across all sites. Our results highlight the importance of seed size and dormancy as key factors impacting seedling establishment and adult recruitment, suggesting a potential interactive relationship between these traits.

Plant physical defenses contribute to a latitudinal gradient in resistance to insect herbivory within a widespread perennial grass.

PREMISE: Herbivore pressure can vary across the range of a species, resulting in different defensive strategies. If herbivory is greater at lower latitudes, plants may be better defended there, potentially driving a latitudinal gradient in defense. However, relationships that manifest across the entire range of a species may be confounded by differences within genetic subpopulations, which may obscure the drivers of these latitudinal gradients. METHODS: We grew plants of the widespread perennial grass Panicum virgatum in a common garden that included genotypes from three genetic subpopulations spanning an 18.5° latitudinal gradient. We then assessed defensive strategies of these plants by measuring two physical resistance traits-leaf mass per area (LMA) and leaf ash, a proxy for silica-and multiple measures of herbivory by caterpillars of the generalist herbivore fall armyworm (Spodoptera frugiperda). RESULTS: Across all genetic subpopulations, low-latitude plants experienced less herbivory than high-latitude plants. Within genetic subpopulations, however, this relationship was inconsistent-the most widely distributed and phenotypically variable subpopulation (Atlantic) exhibited more consistent latitudinal trends than either of the other two subpopulations. The two physical resistance traits, LMA and leaf ash, were both highly heritable and positively associated with resistance to different measures of herbivory across all subpopulations, indicating their importance in defense against herbivores. Again, however, these relationships were inconsistent within subpopulations. CONCLUSIONS: Defensive gradients that occur across the entire species range may not arise within localized subpopulations. Thus, identifying the drivers of latitudinal gradients in herbivory defense may depend on adequately sampling the diversity within a species.

Arabidopsis AN3 and OLIGOCELLULA genes link telomere maintenance mechanisms with cell division and expansion control.

Telomeres are conserved chromosomal structures necessary for continued cell division and proliferation. In addition to the classical telomerase pathway, multiple other genes including those involved in ribosome metabolism and chromatin modification contribute to telomere length maintenance. We previously reported that Arabidopsis thaliana ribosome biogenesis genes OLI2/NOP2A, OLI5/RPL5A and OLI7/RPL5B have critical roles in telomere length regulation. These three OLIGOCELLULA genes were also shown to function in cell proliferation and expansion control and to genetically interact with the transcriptional co-activator ANGUSTIFOLIA3 (AN3). Here we show that AN3-deficient plants progressively lose telomeric DNA in early homozygous mutant generations, but ultimately establish a new shorter telomere length setpoint by the fifth mutant generation with a telomere length similar to oli2/nop2a - deficient plants. Analysis of double an3 oli2 mutants indicates that the two genes are epistatic for telomere length control. Telomere shortening in an3 and oli mutants is not caused by telomerase inhibition; wild type levels of telomerase activity are detected in all analyzed mutants in vitro. Late generations of an3 and oli mutants are prone to stem cell damage in the root apical meristem, implying that genes regulating telomere length may have conserved functional roles in stem cell maintenance mechanisms. Multiple instances of anaphase fusions in late generations of oli5 and oli7 mutants were observed, highlighting an unexpected effect of ribosome biogenesis factors on chromosome integrity. Overall, our data implicate AN3 transcription coactivator and OLIGOCELLULA proteins in the establishment of telomere length set point in plants and further suggest that multiple regulators with pleiotropic functions can connect telomere biology with cell proliferation and cell expansion pathways.

Evolutionary Analyses of Gene Expression Divergence in Panicum hallii: Exploring Constitutive and Plastic Responses Using Reciprocal Transplants.

The evolution of gene expression is thought to be an important mechanism of local adaptation and ecological speciation. Gene expression divergence occurs through the evolution of cis- polymorphisms and through more widespread effects driven by trans-regulatory factors. Here, we explore expression and sequence divergence in a large sample of Panicum hallii accessions encompassing the species range using a reciprocal transplantation experiment. We observed widespread genotype and transplant site drivers of expression divergence, with a limited number of genes exhibiting genotype-by-site interactions. We used a modified FST-QST outlier approach (QPC analysis) to detect local adaptation. We identified 514 genes with constitutive expression divergence above and beyond the levels expected under neutral processes. However, no plastic expression responses met our multiple testing correction as QPC outliers. Constitutive QPC outlier genes were involved in a number of developmental processes and responses to abiotic environments. Leveraging earlier expression quantitative trait loci results, we found a strong enrichment of expression divergence, including for QPC outliers, in genes previously identified with cis and cis-environment interactions but found no patterns related to trans-factors. Population genetic analyses detected elevated sequence divergence of promoters and coding sequence of constitutive expression outliers but little evidence for positive selection on these proteins. Our results are consistent with a hypothesis of cis-regulatory divergence as a primary driver of expression divergence in P. hallii.

Seed size variation impacts local adaptation and life-history strategies in a perennial grass.

Seed mass is an ecologically important trait that often differs considerably among ecotypes. Yet, because few studies examine the impacts of seed mass on adult life-history traits, its role in local adaptation is unclear. In this study, using accessions of Panicum hallii that spanned the two major ecotypes, we examined whether covariation between seed mass, seedling and reproductive traits impacts ecotypic divergence and local adaptation. The perennial grass P. hallii has two distinct ecotypes-a large-seeded upland ecotype adapted to xeric environments and a small-seeded lowland ecotype adapted to mesic environments. In the greenhouse, seed mass varied greatly across P. hallii genotypes in a manner consistent with ecotypic divergence. Seed mass covaried significantly with several seedling and reproductive traits. At field sites representing the habitats of the two ecotypes, seed mass had different impacts on seedling and adult recruitment: selection favoured large seeds in upland habitat and small seeds in lowland habitat, which was consistent with local adaptation. By demonstrating the central role of seed mass in ecotypic differences in P. hallii and its importance to seedling and adult recruitment under field conditions, these studies show that early life-history traits can promote local adaptation and potentially explain ecotype formation.

Genetic determinants of switchgrass-root-associated microbiota in field sites spanning its natural range.

A fundamental goal in plant microbiome research is to determine the relative impacts of host and environmental effects on root microbiota composition, particularly how host genotype impacts bacterial community composition. Most studies characterizing the effect of plant genotype on root microbiota undersample host genetic diversity and grow plants outside of their native ranges, making the associations between host and microbes difficult to interpret. Here, we characterized the root microbiota of a large diversity panel of switchgrass, a North American native C4 bioenergy crop, in three field locations spanning its native range. Our data, composed of 1,961 samples, suggest that field location is the primary determinant of microbiome composition; however, substantial heritable variation is widespread across bacterial taxa, especially those in the Sphingomonadaceae family. Despite diverse compositions, relatively few highly prevalent taxa make up the majority of the switchgrass root microbiota, a large fraction of which is shared across sites. Local genotypes preferentially recruit/filter for local microbes, supporting the idea of affinity between local plants and their microbiota. Using genome-wide association, we identified loci impacting the abundance of >400 microbial strains and found an enrichment of genes involved in immune responses, signaling pathways, and secondary metabolism. We found loci associated with over half of the core microbiota (i.e., microbes in >80% of samples), regardless of field location. Finally, we show a genetic relationship between a basal plant immunity pathway and relative abundances of root microbiota. This study brings us closer to harnessing and manipulating beneficial microbial associations via host genetics.

Transcriptome and DNA methylome divergence of inflorescence development between 2 ecotypes in Panicum hallii.

The morphological diversity of the inflorescence determines flower and seed production, which is critical for plant adaptation. Hall's panicgrass (Panicum hallii, P. hallii) is a wild perennial grass that has been developed as a model to study perennial grass biology and adaptive evolution. Highly divergent inflorescences have evolved between the 2 major ecotypes in P. hallii, the upland ecotype (P. hallii var hallii, HAL2 genotype) with compact inflorescence and large seed and the lowland ecotype (P. hallii var filipes, FIL2 genotype) with an open inflorescence and small seed. Here we conducted a comparative analysis of the transcriptome and DNA methylome, an epigenetic mark that influences gene expression regulation, across different stages of inflorescence development using genomic references for each ecotype. Global transcriptome analysis of differentially expressed genes (DEGs) and co-expression modules underlying the inflorescence divergence revealed the potential role of cytokinin signaling in heterochronic changes. Comparing DNA methylome profiles revealed a remarkable level of differential DNA methylation associated with the evolution of P. hallii inflorescence. We found that a large proportion of differentially methylated regions (DMRs) were located in the flanking regulatory regions of genes. Intriguingly, we observed a substantial bias of CHH hypermethylation in the promoters of FIL2 genes. The integration of DEGs, DMRs, and Ka/Ks ratio results characterized the evolutionary features of DMR-associated DEGs that contribute to the divergence of the P. hallii inflorescence. This study provides insights into the transcriptome and epigenetic landscape of inflorescence divergence in P. hallii and a genomic resource for perennial grass biology.

Quantitative genetic-by-soil microbiome interactions in a perennial grass affect functional traits.

Plants interact with diverse microbiomes that can impact plant growth and performance. Recent studies highlight the potential beneficial aspects of plant microbiomes, including the possibility that microbes facilitate the process of local adaptation in their host plants. Microbially mediated local adaptation in plants occurs when local host genotypes have higher fitness than foreign genotypes because of their affiliation with locally beneficial microbes. Here, plant adaptation results from genetic interactions of the host with locally beneficial microbes (e.g. host genotype-by-microbiome interactions). We used a recombinant inbred line (RIL) mapping population derived from upland and lowland ecotypes of the diploid C4 perennial bunch grass Panicum hallii to explore quantitative genetic responses to soil microbiomes focusing on functional root and shoot traits involved in ecotypic divergence. We show that the growth and development of ecotypes and their trait divergence depends on soil microbiomes. Moreover, we find that the genetic architecture is modified by soil microbiomes, revealing important plant genotype-by-microbiome interactions for quantitative traits. We detected a number of quantitative trait loci (QTL) that interact with the soil microbiome. Our results highlight the importance of microbial interactions in ecotypic divergence and trait genetic architecture in C4 perennial grasses.

Burning questions for a warming and changing world: 15 unknowns in plant abiotic stress.

We present unresolved questions in plant abiotic stress biology as posed by 15 research groups with expertise spanning eco-physiology to cell and molecular biology. Common themes of these questions include the need to better understand how plants detect water availability, temperature, salinity, and rising carbon dioxide (CO2) levels; how environmental signals interface with endogenous signaling and development (e.g. circadian clock and flowering time); and how this integrated signaling controls downstream responses (e.g. stomatal regulation, proline metabolism, and growth versus defense balance). The plasma membrane comes up frequently as a site of key signaling and transport events (e.g. mechanosensing and lipid-derived signaling, aquaporins). Adaptation to water extremes and rising CO2 affects hydraulic architecture and transpiration, as well as root and shoot growth and morphology, in ways not fully understood. Environmental adaptation involves tradeoffs that limit ecological distribution and crop resilience in the face of changing and increasingly unpredictable environments. Exploration of plant diversity within and among species can help us know which of these tradeoffs represent fundamental limits and which ones can be circumvented by bringing new trait combinations together. Better defining what constitutes beneficial stress resistance in different contexts and making connections between genes and phenotypes, and between laboratory and field observations, are overarching challenges.

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.

Gene-by-environment interactions in plants: Molecular mechanisms, environmental drivers, and adaptive plasticity.

Plants demonstrate a broad range of responses to environmental shifts. One of the most remarkable responses is plasticity, which is the ability of a single plant genotype to produce different phenotypes in response to environmental stimuli. As with all traits, the ability of plasticity to evolve depends on the presence of underlying genetic diversity within a population. A common approach for evaluating the role of genetic variation in driving differences in plasticity has been to study genotype-by-environment interactions (G × E). G × E occurs when genotypes produce different phenotypic trait values in response to different environments. In this review, we highlight progress and promising methods for identifying the key environmental and genetic drivers of G × E. Specifically, methodological advances in using algorithmic and multivariate approaches to understand key environmental drivers combined with new genomic innovations can greatly increase our understanding about molecular responses to environmental stimuli. These developing approaches can be applied to proliferating common garden networks that capture broad natural environmental gradients to unravel the underlying mechanisms of G × E. An increased understanding of G × E can be used to enhance the resilience and productivity of agronomic systems.