Early aerial senescence in switchgrass (Panicum virgatum) can significantly limit biomass yields. WRKY transcription factors that can regulate senescence could be used to reprogram senescence and enhance biomass yields. All potential WRKY genes present in the version 1.0 of the switchgrass genome were identified and curated using manual and bioinformatic methods. Expression profiles of WRKY genes in switchgrass flag leaf RNA-Seq datasets were analyzed using clustering and network analyses tools to identify both WRKY and WRKY-associated gene co-expression networks during leaf development and senescence onset. We identified 240 switchgrass WRKY genes including members of the RW5 and RW6 families of resistance proteins. Weighted gene co-expression network analysis of the flag leaf transcriptomes across development readily separated clusters of co-expressed genes into thirteen modules. A visualization highlighted separation of modules associated with the early and senescence-onset phases of flag leaf growth. The senescence-associated module contained 3000 genes including 23 WRKYs. Putative promoter regions of senescence-associated WRKY genes contained several cis-element-like sequences suggestive of responsiveness to both senescence and stress signaling pathways. A phylogenetic comparison of senescence-associated WRKY genes from switchgrass flag leaf with senescence-associated WRKY genes from other plants revealed notable hotspots in Group I, IIb, and IIe of the phylogenetic tree. We have identified and named 240 WRKY genes in the switchgrass genome. Twenty three of these genes show elevated mRNA levels during the onset of flag leaf senescence. Eleven of the WRKY genes were found in hotspots of related senescence-associated genes from multiple species and thus represent promising targets for future switchgrass genetic improvement. Overall, individual WRKY gene expression profiles could be readily linked to developmental stages of flag leaves.

Aphid infestation of switchgrass (Panicum virgatum) has the potential to reduce yields and biomass quality. Although switchgrass-greenbug (Schizaphis graminum; GB) interactions have been studied at the whole plant level, little information is available on plant defense responses at the molecular level. The global transcriptomic response of switchgrass cv Summer to GB was monitored by RNA-Seq in infested and control (uninfested) plants harvested at 5, 10, and 15 days after infestation (DAI). Differentially expressed genes (DEGs) in infested plants were analyzed relative to control uninfested plants at each time point. DEGs in GB-infested plants induced by 5-DAI included an upregulation of reactive burst oxidases and several cell wall receptors. Expression changes in genes linked to redox metabolism, cell wall structure, and hormone biosynthesis were also observed by 5-DAI. At 10-DAI, network analysis indicated a massive upregulation of defense-associated genes, including NAC, WRKY, and MYB classes of transcription factors and potential ancillary signaling molecules such as leucine aminopeptidases. Molecular evidence for loss of chloroplastic functions was also detected at this time point. Supporting these molecular changes, chlorophyll content was significantly decreased, and ROS levels were elevated in infested plants 10-DAI. Total peroxidase and laccase activities were elevated in infested plants at 10-DAI relative to control uninfested plants. The net result appeared to be a broad scale defensive response that led to an apparent reduction in C and N assimilation and a potential redirection of nutrients away from GB and towards the production of defensive compounds, such as pipecolic acid, chlorogenic acid, and trehalose by 10-DAI. By 15-DAI, evidence of recovery in primary metabolism was noted based on transcript abundances for genes associated with carbon, nitrogen, and nutrient assimilation. Extensive remodeling of the plant transcriptome and the production of ROS and several defensive metabolites in an upland switchgrass cultivar were observed in response to GB feeding. The early loss and apparent recovery in primary metabolism by 15-DAI would suggest that these transcriptional changes in later stages of GB infestation could underlie the recovery response categorized for this switchgrass cultivar. These results can be exploited to develop switchgrass lines with more durable resistance to GB and potentially other aphids.

Switchgrass breeders need to improve the rates of genetic gain in many bioenergy-related traits in order to create improved cultivars that are higher yielding and have optimal biomass composition. One way to achieve this is through genomic selection. However, the heritability of traits needs to be determined as well as the accuracy of prediction in order to determine if efficient selection is possible. Using five distinct switchgrass populations comprised of three lowland, one upland and one hybrid accession, the accuracy of genomic predictions under different cross-validation strategies and prediction methods was investigated. Individual genotypes were collected using GBS while kin-BLUP, partial least squares, sparse partial least squares, and BayesB methods were employed to predict yield, morphological, and NIRS-based compositional data collected in 2012–2013 from a replicated Nebraska field trial. Population structure was assessed by F statistics which ranged from 0.3952 between lowland and upland accessions to 0.0131 among the lowland accessions. Prediction accuracy ranged from 0.57–0.52 for cell wall soluble glucose and fructose respectively, to insignificant for traits with low repeatability. Ratios of heritability across to within-population ranged from 15 to 0.6. Accuracy was significantly affected by both cross-validation strategy and trait. Accounting for population structure with a cross-validation strategy constrained by accession resulted in accuracies that were 69% lower than apparent accuracies using unconstrained cross-validation. Less accurate genomic selection is anticipated when most of the phenotypic variation exists between populations such as with spring regreening and yield phenotypes.

Herbaceous perennial plants selected as potential biofuel feedstocks had been understudied at the genomic and functional genomic levels. Recent investments, primarily by the U.S. Department of Energy, have led to the development of a number of molecular resources for bioenergy grasses, such as the partially annotated genome for switchgrass (Panicum virgatum L.), and some related diploid species. In its current version, the switchgrass genome contains 65,878 gene models arising from the A and B genomes of this tetraploid grass. The availability of these gene sequences provides a framework to exploit transcriptomic data obtained from next generation sequencing platforms to address questions of biological importance. One such question pertains to discovery of genes and proteins important for biotic and abiotic stress responses, and how these components might affect biomass quality and stress response in plants engineered for a specific end purpose. It can be expected that production of switchgrass on marginal lands will expose plants to diverse stresses, including herbivory by insects. Class III plant peroxidases have been implicated in many developmental responses such as lignification and in the adaptive responses of plants to insect feeding. Here, we have analyzed the class III peroxidases encoded by the switchgrass genome, and have mined available transcriptomic datasets to develop a first understanding of the expression profiles of the class III peroxidases in different plant tissues. Lastly, we have identified switchgrass peroxidases that appear to be orthologs of enzymes shown to play key roles in lignification and plant defense responses to hemipterans.

Switchgrass (Panicum virgatum L.) is a warm-season perennial grass that can be used as a second generation bioenergy crop. However, foliar fungal pathogens, like switchgrass rust, have the potential to significantly reduce switchgrass biomass yield. Despite its importance as a prominent bioenergy crop, a genome-wide comprehensive analysis of NB-LRR disease resistance genes has yet to be performed in switchgrass. In this study, we used a homology-based computational approach to identify 1011 potential NB-LRR resistance gene homologs (RGHs) in the switchgrass genome (v 1.1). In addition, we identified 40 RGHs that potentially contain unique domains including major sperm protein domain, jacalin-like binding domain, calmodulin-like binding, and thioredoxin. RNA-sequencing analysis of leaf tissue from ‘Alamo’, a rust-resistant switchgrass cultivar, and ‘Dacotah’, a rust-susceptible switchgrass cultivar, identified 2634 high quality variants in the RGHs between the two cultivars. RNA-sequencing data from field-grown cultivar ‘Summer’ plants indicated that the expression of some of these RGHs was developmentally regulated. Our results provide useful insight into the molecular structure, distribution, and expression patterns of members of the NB-LRR gene family in switchgrass. These results also provide a foundation for future work aimed at elucidating the molecular mechanisms underlying disease resistance in this important bioenergy crop.

Switchgrass (Panicum virgatum L) is perennial, C4 grass with great potential as a biofuel crop. An in-depth understanding of the mechanisms that control mineral uptake, distribution and remobilization will benefit sustainable production. Nutrients are mobilized from aerial portions to below-ground crowns and rhizomes as a natural accompaniment to above-ground senescence post seed-set. Mineral uptake and remobilization is dependent on transporters, however, little if any information is available about the specific transporters that are needed and how their relative expression changes over a growing season. Using well-defined classes of mineral transporters, we identified 520 genes belonging to 40 different transporter classes in the tetraploid switchgrass genome. Expression patterns were determined for many of these genes using publically available transcriptomic datasets obtained from both greenhouse and field grown plants. Certain transporters showed strong temporal patterns of expression in distinct developmental stages of the plant. Gene-expression was verified for selected transporters using qRT-PCR. By and large these analyses confirmed the developmental stage-specific expression of these genes. Mineral analyses indicated that K, Fe, Mg, Co and As had a similar pattern of accumulation with apparent limited remobilization at the end of the growing season. These initial analyses will serve as a foundation for more detailed examination of the nutrient biology of switchgrass.