Nutrient and moisture limitations reveal keystone metabolites linking rhizosphere metabolomes and microbiomes.

Plants release a wealth of metabolites into the rhizosphere that can shape the composition and activity of microbial communities in response to environmental stress. The connection between rhizodeposition and rhizosphere microbiome succession has been suggested, particularly under environmental stress conditions, yet definitive evidence is scarce. In this study, we investigated the relationship between rhizosphere chemistry, microbiome dynamics, and abiotic stress in the bioenergy crop switchgrass grown in a marginal soil under nutrient-limited, moisture-limited, and nitrogen (N)-replete, phosphorus (P)-replete, and NP-replete conditions. We combined 16S rRNA amplicon sequencing and LC-MS/MS-based metabolomics to link rhizosphere microbial communities and metabolites. We identified significant changes in rhizosphere metabolite profiles in response to abiotic stress and linked them to changes in microbial communities using network analysis. N-limitation amplified the abundance of aromatic acids, pentoses, and their derivatives in the rhizosphere, and their enhanced availability was linked to the abundance of bacterial lineages from Acidobacteria, Verrucomicrobia, Planctomycetes, and Alphaproteobacteria. Conversely, N-amended conditions increased the availability of N-rich rhizosphere compounds, which coincided with proliferation of Actinobacteria. Treatments with contrasting N availability differed greatly in the abundance of potential keystone metabolites; serotonin and ectoine were particularly abundant in N-replete soils, while chlorogenic, cinnamic, and glucuronic acids were enriched in N-limited soils. Serotonin, the keystone metabolite we identified with the largest number of links to microbial taxa, significantly affected root architecture and growth of rhizosphere microorganisms, highlighting its potential to shape microbial community and mediate rhizosphere plant-microbe interactions.

miR156-PvSPL2 controls culm development by transcriptional repression of switchgrass CYTOKININ OXIDASE/DEHYDROGENASE4.

Culm development in grasses can be controlled by both miR156 and cytokinin. However, the crosstalk between the miR156-SPL module and the cytokinin metabolic pathway remains largely unknown. Here, we found CYTOKININ OXIDASE/DEHYDROGENASE4 (PvCKX4) plays a negative regulatory role in culm development of the bioenergy grass Panicum virgatum (switchgrass). Overexpression of PvCKX4 in switchgrass reduced the internode diameter and length without affecting tiller number. Interestingly, we also found that PvCKX4 was always upregulated in miR156 overexpressing (miR156OE) transgenic switchgrass lines. Additionally, upregulation of either miR156 or PvCKX4 in switchgrass reduced the content of isopentenyl adenine (iP) without affecting trans-zeatin (tZ) accumulation. It is consistent with the evidence that the recombinant PvCKX4 protein exhibited much higher catalytic activity against iP than tZ in vitro. Furthermore, our results showed that miR156-targeted SPL2 bound directly to the promoter of PvCKX4 to repress its expression. Thus, alleviating the SPL2-mediated transcriptional repression of PvCKX4 through miR156 overexpression resulted in a significant increase in cytokinin degradation and impaired culm development in switchgrass. On the contrary, suppressing PvCKX4 in miR156OE transgenic plants restored iP content, internode diameter, and length to wild-type levels. Most strikingly, the double transgenic lines retained the same increased tiller numbers as the miR156OE transgenic line, which yielded more biomass than the wild type. These findings indicate that the miR156-SPL module can control culm development through transcriptional repression of PvCKX4 in switchgrass, which provides a promising target for precise design of shoot architecture to yield more biomass from grasses.

Genome-wide profiling of histone (H3) lysine 4 (K4) tri-methylation (me3) under drought, heat, and combined stresses in switchgrass.

BACKGROUND: Switchgrass (Panicum virgatum L.) is a warm-season perennial (C4) grass identified as an important biofuel crop in the United States. It is well adapted to the marginal environment where heat and moisture stresses predominantly affect crop growth. However, the underlying molecular mechanisms associated with heat and drought stress tolerance still need to be fully understood in switchgrass. The methylation of H3K4 is often associated with transcriptional activation of genes, including stress-responsive. Therefore, this study aimed to analyze genome-wide histone H3K4-tri-methylation in switchgrass under heat, drought, and combined stress. RESULTS: In total, ~ 1.3 million H3K4me3 peaks were identified in this study using SICER. Among them, 7,342; 6,510; and 8,536 peaks responded under drought (DT), drought and heat (DTHT), and heat (HT) stresses, respectively. Most DT and DTHT peaks spanned 0 to + 2000 bases from the transcription start site [TSS]. By comparing differentially marked peaks with RNA-Seq data, we identified peaks associated with genes: 155 DT-responsive peaks with 118 DT-responsive genes, 121 DTHT-responsive peaks with 110 DTHT-responsive genes, and 175 HT-responsive peaks with 136 HT-responsive genes. We have identified various transcription factors involved in DT, DTHT, and HT stresses. Gene Ontology analysis using the AgriGO revealed that most genes belonged to biological processes. Most annotated peaks belonged to metabolite interconversion, RNA metabolism, transporter, protein modifying, defense/immunity, membrane traffic protein, transmembrane signal receptor, and transcriptional regulator protein families. Further, we identified significant peaks associated with TFs, hormones, signaling, fatty acid and carbohydrate metabolism, and secondary metabolites. qRT-PCR analysis revealed the relative expressions of six abiotic stress-responsive genes (transketolase, chromatin remodeling factor-CDH3, fatty-acid desaturase A, transmembrane protein 14C, beta-amylase 1, and integrase-type DNA binding protein genes) that were significantly (P < 0.05) marked during drought, heat, and combined stresses by comparing stress-induced against un-stressed and input controls. CONCLUSION: Our study provides a comprehensive and reproducible epigenomic analysis of drought, heat, and combined stress responses in switchgrass. Significant enrichment of H3K4me3 peaks downstream of the TSS of protein-coding genes was observed. In addition, the cost-effective experimental design, modified ChIP-Seq approach, and analyses presented here can serve as a prototype for other non-model plant species for conducting stress studies.

Genome-wide analysis of the switchgrass YABBY family and functional characterization of PvYABBY14 in response to ABA and GA stress in Arabidopsis.

BACKGROUND: The small YABBY plant-specific transcription factor has a prominent role in regulating plant growth progress and responding to abiotic stress. RESULTS: Here, a total of 16 PvYABBYs from switchgrass (Panicum virgatum L.) were identified and classified into four distinct subgroups. Proteins within the same subgroup exhibited similar conserved motifs and gene structures. Synteny analyses indicated that segmental duplication contributed to the expansion of the YABBY gene family in switchgrass and that complex duplication events occurred in rice, maize, soybean, and sorghum. Promoter regions of PvYABBY genes contained numerous cis-elements related to stress responsiveness and plant hormones. Expression profile analysis indicated higher expression levels of many PvYABBY genes during inflorescence development and seed maturation, with lower expression levels during root growth. Real-time quantitative PCR analysis demonstrated the sensitivity of multiple YABBY genes to PEG, NaCl, ABA, and GA treatments. The overexpression of PvYABBY14 in Arabidopsis resulted in increased root length after treatment with GA and ABA compared to wild-type plants. CONCLUSIONS: Taken together, our study provides the first genome-wide overview of the YABBY transcription factor family, laying the groundwork for understanding the molecular basis and regulatory mechanisms of PvYABBY14 in response to ABA and GA responses in switchgrass.

Genome-wide identification and expression analysis of two-component system genes in switchgrass (Panicum virgatum L.).

The two-component system (TCS) consists of histidine kinase (HK), histidine phosphate transfer protein (HP), and response regulatory factor (RR). It is one of the most crucial components of signal transduction in plants, playing a significant role in regulating plant growth, development, and responses to various abiotic stresses. Although TCS genes have been extensively identified in a variety of plants, the genome-wide recognition and examination of TCS in switchgrass remain unreported. Accordingly, this study identified a total of 87 TCS members in the genome of switchgrass, comprising 20 HK(L)s, 10 HPs, and 57 RRs. Detailed analyses were also conducted on their gene structures, conserved domains, and phylogenetic relationships. Moreover, this study analysed the gene expression profiles across diverse organs and investigated their response patterns to adverse environmental stresses. Results revealed that 87 TCS genes were distributed across 18 chromosomes, with uneven distribution. Expansion of these genes in switchgrass was achieved through both fragment and tandem duplication. PvTCS members are relatively conservative in the evolutionary process, but the gene structure varies significantly. Various cis-acting elements, varying in types and amounts, are present in the promoter region of PvTCSs, all related to plant growth, development, and abiotic stress, due to the TCS gene structure. Protein-protein interaction and microRNA prediction suggest complex interactions and transcriptional regulation among TCS members. Additionally, most TCS members are expressed in roots and stems, with some genes showing organ-specific expression at different stages of leaf and inflorescence development. Under conditions of abiotic stress such as drought, low temperature, high temperature, and salt stress, as well as exogenous abscisic acid (ABA), the expression of most TCS genes is either stimulated or inhibited. Our systematic analysis could offer insight into the characterization of the TCS genes, and further the growth of functional studies in switchgrass.

Blocking miR528 function promotes tillering and regrowth in switchgrass.

MiRNAs have been reported to be the key regulators involving a wide range of biological processes in diverse plant species, but their functions in switchgrass, an important biofuel and forage crop, are largely unknown. Here, we reported the novel function of miR528, which has expanded to four copies in switchgrass, in controlling biomass trait of tillering number and regrowth rate after mowing. Blocking miR528 activity by expressing short tandem target mimic (STTM) increased tiller number and regrowth rate after mowing. The quadruple pvmir528 mutant lines derived from genome editing also showed such improved traits. Degradome and RNA-seq analysis, combined with in situ hybridization assay revealed that up-regulation of two miR528 targets coding for Cu/Zn-SOD enzymes, might be responsible for the improved traits of tillering and regrowth in pvmir528 mutant. Additionally, natural variations in the miR528-SOD interaction exist in C3 and C4 monocot species, implying the distinct regulatory strength of the miR528-SOD module during monocot evolution. Overall, our data illuminated a novel role of miR528 in controlling biomass traits and provided a new target for genetic manipulation-mediated crop improvement.

UDP-glycosyltransferase UGT96C10 functions as a novel detoxification factor for conjugating the activated dinitrotoluene sulfonate in switchgrass.

Dinitrotoluene sulfonates (DNTSes) are highly toxic hazards regulated by the Resource Conservation and Recovery Act (RCRA) in the United States. The trinitrotoluene (TNT) red water formed during the TNT purification process consists mainly of DNTSes. Certain plants, including switchgrass, reed and alfalfa, can detoxify low concentrations of DNTS in TNT red water-contaminated soils. However, the precise mechanism by which these plants detoxify DNTS remains unknown. In order to aid in the development of phytoremediation resources with high DNTS removal rates, we identified and characterized 1-hydroxymethyl-2,4-dinitrobenzene sulfonic acid (HMDNBS) and its glycosylated product HMDNBS O-glucoside as the degradation products of 2,4-DNT-3-SO3Na, the major isoform of DNTS in TNT red water-contaminated soils, in switchgrass via LC-MS/MS- and NMR-based metabolite analyses. Transcriptomic analysis revealed that 15 UDP-glycosyltransferase genes were dramatically upregulated in switchgrass plants following 2,4-DNT-3-SO3Na treatment. We expressed, purified and assayed the activity of recombinant UGT proteins in vitro and identified PvUGT96C10 as the enzyme responsible for the glycosylation of HMDNBS in switchgrass. Overexpression of PvUGT96C10 in switchgrass significantly alleviated 2,4-DNT-3-SO3Na-induced plant growth inhibition. Notably, PvUGT96C10-overexpressing transgenic switchgrass plants removed 83.1% of 2,4-DNT-3-SO3Na in liquid medium after 28 days, representing a 3.2-fold higher removal rate than that of control plants. This work clarifies the DNTS detoxification mechanism in plants for the first time, suggesting that PvUGT96C10 is crucial for DNTS degradation. Our results indicate that PvUGT96C10-overexpressing plants may hold great potential for the phytoremediation of TNT red water-contaminated soils.

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.

The genetic basis of the root economics spectrum in a perennial grass.

Construction economics of plant roots exhibit predictable relationships with root growth, death, and nutrient uptake strategies. Plant taxa with inexpensively constructed roots tend to more precisely explore nutrient hotspots than do those with costly constructed roots but at the price of more frequent tissue turnover. This trade-off underlies an acquisitive to conservative continuum in resource investment, described as the "root economics spectrum (RES)." Yet the adaptive role and genetic basis of RES remain largely unclear. Different ecotypes of switchgrass (Panicum virgatum) display root features exemplifying the RES, with costly constructed roots in southern lowland and inexpensively constructed roots in northern upland ecotypes. We used an outbred genetic mapping population derived from lowland and upland switchgrass ecotypes to examine the genetic architecture of the RES. We found that absorptive roots (distal first and second orders) were often "deciduous" in winter. The percentage of overwintering absorptive roots was decreased by northern upland alleles compared with southern lowland alleles, suggesting a locally-adapted conservative strategy in warmer and acquisitive strategy in colder regions. Relative turnover of absorptive roots was genetically negatively correlated with their biomass investment per unit root length, suggesting that the key trade-off in framing RES is genetically facilitated. We also detected strong genetic correlations among root morphology, root productivity, and shoot size. Overall, our results reveal the genetic architecture of multiple traits that likely impacts the evolution of RES and plant aboveground-belowground organization. In practice, we provide genetic evidence that increasing switchgrass yield for bioenergy does not directly conflict with enhancing its root-derived carbon sequestration.

The effect of AMF combined with biochar on plant growth and soil quality under saline-alkali stress: Insights from microbial community analysis.

Arbuscular mycorrhizal fungi (AMF) and biochar application individually can enhance plant tolerance to saline-alkali stress and promote plant growth efficiency. However, little is known about the potential synergistic effects of their combination on improving plant growth and soil quality under saline-alkali stress. This experiment adopted the potted method to explore the effects of four treatments on switchgrass growth and soil quality: biochar (BC), Rhizophagus irregularis (Ri), biochar + Ri (BR) and a control without biochar or Ri (CK). Compared to the CK treatment, the switchgrass biomass increased by 92.4 %, 148.6 %, and 177.3 % in the BC, Ri, and BR treatment groups, respectively. Similarly, the rhizosphere soil quality index increased by 29.33 %, 22.7 %, and 49.1 % in the respective treatment groups. The BR treatment significantly altered the rhizosphere soil microbial composition and diversity. Notably, compared to the other treatments, the archaeal α-diversity in the BR group showed a significant decrease. BR treatment significantly increased the relative abundance of bacteria, fungi and archaea at the genus level (e.g., Bacillus, Trichome and candidatus_methanopenens). Network analysis showed that the complexity and closeness of interactions between different microbial taxa were stronger in the BC, Ri and BR treatments than in the CK treatment, with BR being the more prominent. In summary, biochar combined with Ri has a better effect on promoting the growth of switchgrass under saline-alkali stress, improving the quality of saline-alkali soil, and increasing soil microbial diversity. This study provides a new approach for the efficient development and utilization of saline-alkali land.

First report of ergot caused by Claviceps clavispora in switchgrass (Panicum virgatum L.) in Georgia, USA.

Switchgrass (Panicum virgatum L.) is a North American grass species with biofuel potential. Claviceps spp. is known to infect the florets of various grass species, initially characterized by a sticky honeydew exudate and later as sclerotium replacing the seed in the infected ovary (Tanaka et al 2023). Since 2019, from July to October, switchgrass panicles in Georgia have been observed with honeydew and black sclerotia. The disease was first noted on some of the 285 accessions of the genome-wide association study (GWAS) panels at the University of Georgia's Iron Horse Farm in Watkinsville and Gibbs Farm in Tifton, GA. In Watkinsville, GA, ergot incidence was 5%, 6%, 65% and 54% in 2019, 2020, 2021 and 2022, respectively. Symptomatic panicles with honeydew and sclerotia were collected in 2021 (sample Scl) and 2023 (samples Cla_M and ATH20cl) from Watkinsville, GA. Under microcopy, panicles with honeydew symptoms had mycelium and conidia (9 µm long [range of 5-13 µm] and 4 µm wide [3-5 µm]) consistent with Claviceps spp. (Tooley et al. 2001). Sclerotia were 1.5 mm long (range 1-3 mm). Sclerotia were surface sterilized for 3 minutes in 5% NaOCl, followed by 70% ethanol then rinsed three times in distilled water. Sterilized sclerotia were plated on potato dextrose agar and placed on bench top with a 12-hour day/night cycle at room temperature (22°C) (Singh 1976). After 2 months, the sclerotia produced sterile apothecia. For molecular identification, genomic DNA was extracted from three honeydew samples following the protocol of Doyle and Doyle (1987). The internal transcribed spacer (ITS) region and RNA polymerase second largest subunit (RPB2) gene were amplified using ITS4/ITS5 (White et al. 1990) and 5F2/7CR (Liu et al. 1999) primer sets. The ITS region of ATH20cl, Cla_M, and Scl (GenBank nos. PP546317- PP546319) showed 92.83-97.48% identity to C. clavispora (NR_163506.1). The RPB2 region of ATH20cl, Cla_M, and Scl (GenBank nos. PP573916- PP573918) showed 97.46-97.72% identity to C. clavispora (LT216566.1). The maximum likelihood tree constructed in MEGA-X (Kumar et al. 2018) using concatenated ITS (539 bp) and RPB2 (792 bp) gene sequences from this study and eleven reference sequences from Tanaka et al. (2023), revealed close relatedness of ATH20cl, Cla_M, and Scl to C. clavispora under section Pusillae. The pathogenicity test for samples ATH20cl and Cla_M was conducted in the greenhouse on switchgrass cultivar 'Alamo' grown in injection molded pots containing Sungro professional growing mix. Three replicates plants at reproductive (R3) growth stage were inoculated by immersing panicles in 105 spores/ml suspension for 5 minutes and bagged for 3 days (Tooley et al. 2001). Control plants were immersed in distilled water. Honeydew symptoms and sclerotia appeared within 7- and 90-days post-inoculation, respectively, whereas control plants remained symptom-free. The honeydew collected from the infected Alamo panicles were reconfirmed to produce similar spores to those collected from the field. This is the first report of ergot caused by C. clavispora in switchgrass. Ergot was also observed since 2019 on some accessions from the same GWAS panel planted at the University of Tennessee Plant Sciences Farm in Knoxville, TN. This information will be beneficial to determine the effect of ergot on biomass production and seed quality of switchgrass in the southeastern U.S.

Consolidated bioprocessing of lignocellulosic wastes in Northwest China for D-glucaric acid production by an artificial microbial consortium.

D-glucaric acid is a platform chemical of great importance and the consolidated bioprocessing (CBP) of lignocellulose by the microbial consortium of Trichoderma reesei C10 and Saccharomyces cerevisiae LGA-1C3S2 features prospects in biomanufacturing it. Here we compared some representative lignocelluloses in Northwest China including corn stover, wheat straw and switchgrass, and the leading pretreatments including steam explosion, subcritical water pretreatment, sodium hydroxide pretreatment, aqueous ammonia pretreatment, lime pretreatment, and diluted sulfuric acid pretreatment. It was found that sodium hydroxide pretreated switchgrass (SHPSG) was the best substrate for D-glucaric acid production, resulting in the highest D-glucaric acid titers, 11.69 ± 0.73 g/L in shake flask and 15.71 ± 0.80 g/L in 10L airlift fermenter, respectively. To the best of our knowledge, this is the highest D-glucaric acid production titer from lignocellulosic biomass. This work offers a paradigm of producing low-cost D-glucaric acid for low-carbon polyethylene 2,5-furandicarboxylate (PEF) and a reference on developing biorefinery in Northwest China.

A single amino acid change led to structural and functional differentiation of PvHd1 to control flowering in switchgrass.

Switchgrass, a forage and bioenergy crop, occurs as two main ecotypes with different but overlapping ranges of adaptation. The two ecotypes differ in a range of characteristics, including flowering time. Flowering time determines the duration of vegetative development and therefore biomass accumulation, a key trait in bioenergy crops. No causal variants for flowering time differences between switchgrass ecotypes have, as yet, been identified. In this study, we mapped a robust flowering time quantitative trait locus (QTL) on chromosome 4K in a biparental F2 population and characterized the flowering-associated transcription factor gene PvHd1, an ortholog of CONSTANS in Arabidopsis and Heading date 1 in rice, as the underlying causal gene. Protein modeling predicted that a serine to glycine substitution at position 35 (p.S35G) in B-Box domain 1 greatly altered the global structure of the PvHd1 protein. The predicted variation in protein compactness was supported in vitro by a 4 °C shift in denaturation temperature. Overexpressing the PvHd1-p.35S allele in a late-flowering CONSTANS-null Arabidopsis mutant rescued earlier flowering, whereas PvHd1-p.35G had a reduced ability to promote flowering, demonstrating that the structural variation led to functional divergence. Our findings provide us with a tool to manipulate the timing of floral transition in switchgrass cultivars and, potentially, expand their cultivation range.

Plant Host Traits Mediated by Foliar Fungal Symbionts and Secondary Metabolites.

Fungal symbionts living inside plant leaves ("endophytes") can vary from beneficial to parasitic, but the mechanisms by which the fungi affect the plant host phenotype remain poorly understood. Chemical interactions are likely the proximal mechanism of interaction between foliar endophytes and the plant, as individual fungal strains are often exploited for their diverse secondary metabolite production. Here, we go beyond single strains to examine commonalities in how 16 fungal endophytes shift plant phenotypic traits such as growth and physiology, and how those relate to plant metabolomics profiles. We inoculated individual fungi on switchgrass, Panicum virgatum L. This created a limited range of plant growth and physiology (2-370% of fungus-free controls on average), but effects of most fungi overlapped, indicating functional similarities in unstressed conditions. Overall plant metabolomics profiles included almost 2000 metabolites, which were broadly correlated with plant traits across all the fungal treatments. Terpenoid-rich samples were associated with larger, more physiologically active plants and phenolic-rich samples were associated with smaller, less active plants. Only 47 metabolites were enriched in plants inoculated with fungi relative to fungus-free controls, and of these, Lasso regression identified 12 metabolites that explained from 14 to 43% of plant trait variation. Fungal long-chain fatty acids and sterol precursors were positively associated with plant photosynthesis, conductance, and shoot biomass, but negatively associated with survival. The phytohormone gibberellin, in contrast, was negatively associated with plant physiology and biomass. These results can inform ongoing efforts to develop metabolites as crop management tools, either by direct application or via breeding, by identifying how associations with more beneficial components of the microbiome may be affected.

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.

Genome-wide characterization of LBD transcription factors in switchgrass (Panicum virgatum L.) and the involvement of PvLBD12 in salt tolerance.

KEY MESSAGE: PvLBD12 enhanced the salt tolerance by increasing proline accumulation, improving K+ accumulation, and decreasing reactive oxygen species level in switchgrass. Abiotic stresses are the serious factors which limit plant development and productivity and restrict the agricultural economy. It is important, therefore, to understand the mechanism of abiotic tolerance in plants. Lateral organ boundaries domain (LBD) proteins as plant-specific transcription factors play important function in plant lateral organ development, plant regeneration, and abiotic stress. In our study, we identify 69 LBD members from switchgrass genome-wide sequences and classify them based on their homology with LBD proteins in Arabidopsis. RT-qPCR showed that PvLBD genes had different expression patterns under abiotic stress conditions, indicating that they play important roles in various stress. PvLBD12 was selected as a candidate gene for further functional analysis because it had the highest expression level under salt stress. Overexpression of PvLBD12 enhanced salt tolerance by altering a wide range of physiological responses (like increased proline accumulation, reduced malondialdehyde production, improved K+ accumulation, and reduced Na+ absorption) in switchgrass. Some stress response genes such as proline biosynthesis gene PvP5CS1, vacuolar Na+(K+)/H+ antiporter gene PvNHX1, two key ROS-scavenging enzyme genes PvCAT and PvSOD were all upregulated in PvLBD12 overexpression lines. Taken together, PvLBD12 plays a pivotal role in response to salt stress by increasing proline accumulation, improving K+ accumulation, reducing Na+ absorption, and decreasing reactive oxygen species level. It will be better to understand the potential biological functions of LBD genes in other plants.

Ectopic Expression of PvHMA2.1 Enhances Cadmium Tolerance in Arabidopsis thaliana.

Cadmium (Cd) in soil inhibits plant growth and development and even harms human health through food chain transmission. Switchgrass (Panicum virgatum L.), a perennial C4 biofuel crop, is considered an ideal plant for phytoremediation due to its high efficiency in removing Cd and other heavy metals from contaminated soil. The key to understanding the mechanisms of switchgrass Cd tolerance is to identify the genes involved in Cd transport. Heavy-metal ATPases (HMAs) play pivotal roles in heavy metal transport, including Cd, in Arabidopsis thaliana and Oryza sativa, but little is known about the functions of their orthologs in switchgrass. Therefore, we identified 22 HMAs in switchgrass, which were distributed on 12 chromosomes and divided into 4 groups using a phylogenetic analysis. Then, we focused on PvHMA2.1, which is one of the orthologs of the rice Cd transporter OsHMA2. We found that PvHMA2.1 was widely expressed in roots, internodes, leaves, spikelets, and inflorescences, and was significantly induced in the shoots of switchgrass under Cd treatment. Moreover, PvHMA2.1 was found to have seven transmembrane domains and localized at the cell plasma membrane, indicating that it is a potential transporter. The ectopic expression of PvHMA2.1 alleviated the reduction in primary root length and the loss of fresh weight of Arabidopsis seedlings under Cd treatment, suggesting that PvHMA2.1 enhanced Cd tolerance in Arabidopsis. The higher levels of relative water content and chlorophyll content of the transgenic lines under Cd treatment reflected that PvHMA2.1 maintained water retention capacity and alleviated photosynthesis inhibition under Cd stress in Arabidopsis. The roots of the PvHMA2.1 ectopically expressed lines accumulated less Cd compared to the WT, while no significant differences were found in the Cd contents of the shoots between the transgenic lines and the WT under Cd treatment, suggesting that PvHMA2.1 reduced Cd absorption from the environment through the roots in Arabidopsis. Taken together, our results showed that PvHMA2.1 enhanced Cd tolerance in Arabidopsis, providing a promising target that could be engineered in switchgrass to repair Cd-contaminated soil.

Dynamic Reconfiguration of Switchgrass Proteomes in Response to Rust (Puccinia novopanici) Infection.

Switchgrass (Panicum virgatum L.) can be infected by the rust pathogen (Puccinia novopanici) and results in lowering biomass yields and quality. Label-free quantitative proteomics was conducted on leaf extracts harvested from non-infected and infected plants from a susceptible cultivar (Summer) at 7, 11, and 18 days after inoculation (DAI) to follow the progression of disease and evaluate any plant compensatory mechanisms to infection. Some pustules were evident at 7 DAI, and their numbers increased with time. However, fungal DNA loads did not appreciably change over the course of this experiment in the infected plants. In total, 3830 proteins were identified at 1% false discovery rate, with 3632 mapped to the switchgrass proteome and 198 proteins mapped to different Puccinia proteomes. Across all comparisons, 1825 differentially accumulated switchgrass proteins were identified and subjected to a STRING analysis using Arabidopsis (A. thaliana L.) orthologs to deduce switchgrass cellular pathways impacted by rust infection. Proteins associated with plastid functions and primary metabolism were diminished in infected Summer plants at all harvest dates, whereas proteins associated with immunity, chaperone functions, and phenylpropanoid biosynthesis were significantly enriched. At 18 DAI, 1105 and 151 proteins were significantly enriched or diminished, respectively. Many of the enriched proteins were associated with mitigation of cellular stress and defense.

Activity of Cytosolic Ascorbate Peroxidase (APX) from Panicum virgatum against Ascorbate and Phenylpropanoids.

APX is a key antioxidant enzyme in higher plants, scavenging H2O2 with ascorbate in several cellular compartments. Here, we report the crystal structures of cytosolic ascorbate peroxidase from switchgrass (Panicum virgatum L., Pvi), a strategic feedstock plant with several end uses. The overall structure of PviAPX was similar to the structures of other APX family members, with a bound ascorbate molecule at the ɣ-heme edge pocket as in other APXs. Our results indicated that the H2O2-dependent oxidation of ascorbate displayed positive cooperativity. Significantly, our study suggested that PviAPX can oxidize a broad range of phenylpropanoids with δ-meso site in a rather similar efficiency, which reflects its role in the fortification of cell walls in response to insect feeding. Based on detailed structural and kinetic analyses and molecular docking, as well as that of closely related APX enzymes, the critical residues in each substrate-binding site of PviAPX are proposed. Taken together, these observations shed new light on the function and catalysis of PviAPX, and potentially benefit efforts improve plant health and biomass quality in bioenergy and forage crops.

Structural Similarities and Overlapping Activities among Dihydroflavonol 4-Reductase, Flavanone 4-Reductase, and Anthocyanidin Reductase Offer Metabolic Flexibility in the Flavonoid Pathway.

Flavonoids are potent antioxidants that play a role in defense against pathogens, UV-radiation, and the detoxification of reactive oxygen species. Dihydroflavonol 4-reductase (DFR) and flavanone 4-reductase (FNR) reduce dihydroflavonols and flavanones, respectively, using NAD(P)H to produce flavan-(3)-4-(di)ols in flavonoid biosynthesis. Anthocyanidin reductase (ANR) reduces anthocyanidins to flavan-3-ols. In addition to their sequences, the 3D structures of recombinant DFR, FNR and ANR from sorghum and switchgrass showed a high level of similarity. The catalytic mechanism, substrate-specificity and key residues of three reductases were deduced from crystal structures, site-directed mutagenesis, molecular docking, kinetics, and thermodynamic ana-lyses. Although DFR displayed its highest activity against dihydroflavonols, it also showed activity against flavanones and anthocyanidins. It was inhibited by the flavonol quercetin and high concentrations of dihydroflavonols/flavonones. SbFNR1 and SbFNR2 did not show any activity against dihydroflavonols. However, SbFNR1 displayed activity against flavanones and ANR activity against two anthocyanidins, cyanidin and pelargonidin. Therefore, SbFNR1 and SbFNR2 could be specific ANR isozymes without delphinidin activity. Sorghum has high concentrations of 3-deoxyanthocyanidins in vivo, supporting the observed high activity of SbDFR against flavonols. Mining of expression data indicated substantial induction of these three reductase genes in both switchgrass and sorghum in response to biotic stress. Key signature sequences for proper DFR/ANR classification are proposed and could form the basis for future metabolic engineering of flavonoid metabolism.