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.

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.

Regional biogeography versus intra-annual dynamics of the root and soil microbiome

BACKGROUND

Root and soil microbial communities constitute the below-ground plant microbiome, are drivers of nutrient cycling, and affect plant productivity. However, our understanding of their spatiotemporal patterns is confounded by exogenous factors that covary spatially, such as changes in host plant species, climate, and edaphic factors. These spatiotemporal patterns likely differ across microbiome domains (bacteria and fungi) and niches (root vs. soil).

RESULTS

To capture spatial patterns at a regional scale, we sampled the below-ground microbiome of switchgrass monocultures of five sites spanning > 3 degrees of latitude within the Great Lakes region. To capture temporal patterns, we sampled the below-ground microbiome across the growing season within a single site. We compared the strength of spatiotemporal factors to nitrogen addition determining the major drivers in our perennial cropping system. All microbial communities were most strongly structured by sampling site, though collection date also had strong effects; in contrast, nitrogen addition had little to no effect on communities. Though all microbial communities were found to have significant spatiotemporal patterns, sampling site and collection date better explained bacterial than fungal community structure, which appeared more defined by stochastic processes. Root communities, especially bacterial, were more temporally structured than soil communities which were more spatially structured, both across and within sampling sites. Finally, we characterized a core set of taxa in the switchgrass microbiome that persists across space and time. These core taxa represented < 6% of total species richness but > 27% of relative abundance, with potential nitrogen fixing bacteria and fungal mutualists dominating the root community and saprotrophs dominating the soil community.

CONCLUSIONS

Our results highlight the dynamic variability of plant microbiome composition and assembly across space and time, even within a single variety of a plant species. Root and soil fungal community compositions appeared spatiotemporally paired, while root and soil bacterial communities showed a temporal lag in compositional similarity suggesting active recruitment of soil bacteria into the root niche throughout the growing season. A better understanding of the drivers of these differential responses to space and time may improve our ability to predict microbial community structure and function under novel conditions.

Prevalence of Deltaproteobacterial sequences in nifH gene pools associated with the rhizosphere of native switchgrass from Tall Grass Prairie (Oklahoma, USA)

The potential nitrogen-fixing bacterial diversity in the rhizospheric soil of the native switchgrass (Panicum virgatum L.) from Tall Grass Prairies of Northern Oklahoma was studied using a partial region of nitrogenase structural gene-nifH. Eleven clone libraries constructed from nifH amplicons gave 407 good-quality sequences. More than 70% of sequences showed similarity of nifH with uncultured bacteria (< 98%). The dominance of sequences affiliated with Deltaproteobacterial nifH was observed, followed by Betaproteobacterial nifH sequences. The nifH gene library was dominated by the genera Geobacter, Rhizobacter, Paenibacillus, and Azoarcus. Sequences affiliated with rhizobia, such as Bradyrhizobium, Methylocystis, Ensifer, etc., were also in the rhizosphere in small numbers. From Deltaproteobacteria, five genera, namely Geobacter, Pelobacter, Geomonas, Desulfovibrio, and Anaeromyxobacter, contributed to 48% of the total sequences suggesting the dominance of group Deltaproteobacteria in the rhizosphere of native switchgrass. Considering the percent similarity of the nifH sequences with cultivated bacteria, this study demonstrated the presence of novel bacterial species in switchgrass rhizospheric soil from Tall Grass Prairie.

Molecular evolution, diversification, and expression assessment of MADS gene family in Setaria italica, Setaria viridis, and Panicum virgatum

KEY MESSAGE

This paper sheds light on the evolution and expression patterns of MADS genes in Setaria and Panicum virgatum. SiMADS51 and SiMADS64 maybe involved in the ABA-dependent pathway of drought response. The MADS gene family is a key regulatory factor family that controls growth, reproduction, and response to abiotic stress in plants. However, the molecular evolution of this family is rarely reported. Here, a total of 265 MADS genes were identified in Setaria italica (foxtail millet), Setaria viridis (green millet), and Panicum virgatum (switchgrass) and analyzed by bioinformatics, including physicochemical characteristics, subcellular localization, chromosomal position and duplicate, motif distribution, genetic structure, genetic evolvement, and expression patterns. Phylogenetic analysis was used to categorize these genes into M and MIKC types. The distribution of motifs and gene structure were similar for the corresponding types. According to a collinearity study, the MADS genes have been mostly conserved during evolution. The principal cause of their expansion is segmental duplication. However, the MADS gene family tends to shrink in foxtail millet, green millet, and switchgrass. The MADS genes were subjected to purifying selection, but several positive selection sites were also identified in three species. And most of the promoters of MADS genes contain cis-elements related to stress and hormonal response. RNA-seq and quantitative Real-time PCR (qRT-PCR) analysis also were examined. SiMADS genes expression levels are considerably changed in reaction to various treatments, following qRT-PCR analysis. This sheds fresh light on the evolution and expansion of the MADS family in foxtail millet, green millet, and switchgrass, and lays the foundation for further research on their functions.

Field performance of switchgrass plants engineered for reduced recalcitrance

Switchgrass (Panicum virgatum L.) is a promising perennial bioenergy crop that achieves high yields with relatively low nutrient and energy inputs. Modification of cell wall composition for reduced recalcitrance can lower the costs of deconstructing biomass to fermentable sugars and other intermediates. We have engineered overexpression of OsAT10, encoding a rice BAHD acyltransferase and QsuB, encoding dehydroshikimate dehydratase from Corynebacterium glutamicum, to enhance saccharification efficiency in switchgrass. These engineering strategies demonstrated low lignin content, low ferulic acid esters, and increased saccharification yield during greenhouse studies in switchgrass and other plant species. In this work, transgenic switchgrass plants overexpressing either OsAT10 or QsuB were tested in the field in Davis, California, USA for three growing seasons. No significant differences in the content of lignin and cell wall-bound p-coumaric acid or ferulic acid were detected in transgenic OsAT10 lines compared with the untransformed Alamo control variety. However, the transgenic overexpressing QsuB lines had increased biomass yield and slightly increased biomass saccharification properties compared to the control plants. This work demonstrates good performance of engineered plants in the field, and also shows that the cell wall changes in the greenhouse were not replicated in the field, emphasizing the need to validate engineered plants under relevant field conditions.

Could a Legume-Switchgrass Sod-Seeding System Increase Forage Productivity?

Nowadays, the lack of cattle feed, particularly green fodder, has become a key limiting factor in the agricultural economy. Switchgrass appears to offer a viable solution to the feed shortage. An improved cultivation practice might be needed to boost switchgrass forage production all season long. This study was conducted to quantify the positive effects of introducing different legume crops (vetch and pea), optimally fertilized, on the production and quality of mixed harvested switchgrass-legumes hay in late spring (May) and switchgrass hay harvested once more in early fall (September). The studied intercropping systems, independently of the legume species used, increased forage productivity (almost threefold), reaching 7.5 t ha^-1 and quality characteristics, with protein content almost rising threefold, reaching 12.5%. The aforementioned practice can assist the perennial crop (switchgrass) in providing a high hay production during the early fall harvest, even without fertilization. The overall annual economic benefit for the farmers may be increased by 90-720 € per ha, depending on the prevailing weather conditions. Overall, it may be concluded that the suggested cropping system produces a significantly higher yield of cattle feed compared to traditional monocultures, improving the agricultural economy while reducing the negative effects of modern agriculture on the environment.

Identification of microRNAs responsive to arbuscular mycorrhizal fungi in Panicum virgatum (switchgrass)

BACKGROUND

MicroRNAs (miRNAs) are important post-transcriptional regulators involved in the control of a range of processes, including symbiotic interactions in plants. MiRNA involvement in arbuscular mycorrhizae (AM) symbiosis has been mainly studied in model species, and our study is the first to analyze global miRNA expression in the roots of AM colonized switchgrass (Panicum virgatum), an emerging biofuel feedstock. AM symbiosis helps plants gain mineral nutrition from the soil and may enhance switchgrass biomass production on marginal lands. Our goals were to identify miRNAs and their corresponding target genes that are controlling AM symbiosis in switchgrass.

RESULTS

Through genome-wide analysis of next-generation miRNA sequencing reads generated from switchgrass roots, we identified 122 mature miRNAs, including 28 novel miRNAs. By comparing miRNA expression profiles of AM-inoculated and control switchgrass roots, we identified 15 AM-responsive miRNAs across lowland accession "Alamo", upland accession "Dacotah", and two upland/lowland F₁ hybrids. We used degradome sequencing to identify target genes of the AM-responsive miRNAs revealing targets of miRNAs residing on both K and N subgenomes. Notably, genes involved in copper ion binding were targeted by downregulated miRNAs, while upregulated miRNAs mainly targeted GRAS family transcription factors.

CONCLUSION

Through miRNA analysis and degradome sequencing, we revealed that both upland and lowland switchgrass genotypes as well as upland-lowland hybrids respond to AM by altering miRNA expression. We demonstrated complex GRAS transcription factor regulation by the miR171 family, with some miR171 family members being AM responsive while others remained static. Copper miRNA downregulation was common amongst the genotypes tested and we identified superoxide dismutases and laccases as targets, suggesting that these Cu-miRNAs are likely involved in ROS detoxification and lignin deposition, respectively. Other prominent targets of the Cu miRNAs were blue copper proteins. Overall, the potential effect of AM colonization on lignin deposition pathways in this biofuel crop highlights the importance of considering AM and miRNA in future biofuel crop development strategies.

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.

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

Polyploidy, which occurs in roughly half of all flowering plants and an even higher percentage of grasses, is thought to be a major driver of adaptation. Higher numbers of copies of each gene in polyploid genomes can increase genetic diversity, which could drive shifts in habitat preference, adaptability, and fitness. To test the effects of increased ploidy, we compared genomic diversity, environmental niche, and fitness responses across climatic gradients between tetraploid and octoploid switchgrass. We found that the octoploids contained novel combinations of the ancestral tetraploid genetic diversity, which was linked to the expansion of switchgrass into unsuitable habitats for tetraploid populations. Our experiments revealed evidence of niche divergence, differential fitness, and a generalist–specialist trade-off between cytotypes.

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.

Global analysis of switchgrass (Panicum virgatum L.) transcriptomes in response to interactive effects of drought and heat stresses

BACKGROUND

Sustainable production of high-quality feedstock has been of great interest in bioenergy research. Despite the economic importance, high temperatures and water deficit are limiting factors for the successful cultivation of switchgrass in semi-arid areas. There are limited reports on the molecular basis of combined abiotic stress tolerance in switchgrass, particularly the combination of drought and heat stress. We used transcriptomic approaches to elucidate the changes in the response of switchgrass to drought and high temperature simultaneously.

RESULTS

We conducted solely drought treatment in switchgrass plant Alamo AP13 by withholding water after 45 days of growing. For the combination of drought and heat effect, heat treatment (35 °C/25 °C day/night) was imposed after 72 h of the initiation of drought. Samples were collected at 0 h, 72 h, 96 h, 120 h, 144 h, and 168 h after treatment imposition, total RNA was extracted, and RNA-Seq conducted. Out of a total of 32,190 genes, we identified 3912, as drought (DT) responsive genes, 2339 and 4635 as, heat (HT) and drought and heat (DTHT) responsive genes, respectively. There were 209, 106, and 220 transcription factors (TFs) differentially expressed under DT, HT and DTHT respectively. Gene ontology annotation identified the metabolic process as the significant term enriched in DTHT genes. Other biological processes identified in DTHT responsive genes included: response to water, photosynthesis, oxidation-reduction processes, and response to stress. KEGG pathway enrichment analysis on DT and DTHT responsive genes revealed that TFs and genes controlling phenylpropanoid pathways were important for individual as well as combined stress response. For example, hydroxycinnamoyl-CoA shikimate/quinate hydroxycinnamoyl transferase (HCT) from the phenylpropanoid pathway was induced by single DT and combinations of DTHT stress.

CONCLUSION

Through RNA-Seq analysis, we have identified unique and overlapping genes in response to DT and combined DTHT stress in switchgrass. The combination of DT and HT stress may affect the photosynthetic machinery and phenylpropanoid pathway of switchgrass which negatively impacts lignin synthesis and biomass production of switchgrass. The biological function of genes identified particularly in response to DTHT stress could further be confirmed by techniques such as single point mutation or RNAi.

Genome-wide analysis and characterization of GRAS family in switchgrass

Panicum virgatum, a model plant of cellulosic ethanol conversion, not only has high large biomass and strong adaptability to soil, but also grows well in marginal soil and has the advantage of improving saline-alkali soil. GRAS transcription factor gene family play important roles in individual environment adaption, and these vital functions has been proved in several plants, however, the research of GRAS in the development of switchgrass (Panicum virgatum) were limited. A comprehensive study was investigated to explore the relationship between GRAS gene family and resistance. According to the phylogenetic analysis, a total of 144 GRAS genes were identified and renamed which were classified into eight subfamilies. Chromosome distribution, tandem and segmental repeats analysis indicated that gene duplication events contributed a lot to the expansion of GRAS genes in the switchgrass genome. Sixty-six GRAS genes in switchgrass were identified as having orthologous genes with rice through gene duplication analysis. Most of these GRAS genes contained zero or one intron, and closely related genes in evolution shared similar motif composition. Interaction networks were analyzed including DELLA and ten interaction proteins that were primarily involved in gibberellin acid mediated signaling. Notably, online analysis indicated that the promoter regions of the identified PvGRAS genes contained many cis-elements including light responsive elements, suggesting that PvGRAS might involve in light signal cross-talking. This work provides key insights into resistance and bioavailability in switchgrass and would be helpful to further study the function of GRAS and GRAS-mediated signal transduction pathways.

Contributions of environmental and maternal transmission to the assembly of leaf fungal endophyte communities

Leaf fungal endophytes (LFEs) contribute to plant growth and responses to stress. Fungi colonize leaves through maternal transmission, e.g. via the seed, and through environmental transmission, e.g. via aerial dispersal. The relative importance of these two pathways in assembly and function of the LFE community is poorly understood. We used amplicon sequencing to track switchgrass (Panicum virgatum) LFEs in a greenhouse and field experiment as communities assembled from seed endophytes and rain fungi (integration of wet and dry aerial dispersal) in germinating seeds, seedlings, and adult plants. Rain fungi varied temporally and hosted a greater portion of switchgrass LFE richness (greater than 65%) than were found in seed endophytes (greater than 25%). Exposure of germinating seeds to rain inoculum increased dissimilarity between LFE communities and seed endophytes, increasing the abundance of rain-derived taxa, but did not change diversity. In the field, seedling LFE composition changed more over time, with a decline in seed-derived taxa and an increase in richness, in response to environmental transmission than LFEs of adult plants. We show that environmental transmission is an important driver of LFE assembly, and likely plant growth, but its influence depends on both the conditions at the time of colonization and plant life stage.

Phytoremediation potential of three terrestrial plant species for removal of atrazine, azoxystrobin, and imidacloprid

Pesticides can move off-target resulting in contamination of sensitive water bodies and causing adverse effects on inhabiting species. Through best management practices, such as the implementation of vegetative buffer strips, off-target movement of pesticides can be decreased, and compound degradation can be increased via phytoremediation. In this study, blueflag iris (Iris versicolor), broomsedge (Andropogon virginicus) and switchgrass (Panicum virgatum) were planted in soil treated with one of three commonly used pesticides. At 28, 56 and 112 days after treatment (DAT), plants were destructively harvested and analyzed for pesticide residue in soil and above-ground and below-ground vegetation using high-performance liquid chromatography (HPLC). Relative to the amount of pesticide found in planted pots compared to non-planted pots, I. versicolor was found to reduce greater atrazine in soil compared to non-planted pots at 112 DAT by 58.7%. I. versicolor was also the most capable of reducing azoxystrobin, by 86.9% compared to non-planted pots, from the soil at 112 DAT. At the same sampling time, I. versicolor and P. virgatum reduced greater imidacloprid from soil by 62.5% and 64.3% compared to non-planted pots, respectively. This information supports the recommendation for establishment of diverse plant species for optimization of phytoremediation capacities. Novelty statement While research has found that plants can absorb and remediate synthetic chemicals, this practice is only sustainable if used with native plants that require low maintenance and are tolerant to the applied substances. Various previous studies observe plants that are fast-growing, tolerant to environmental conditions, require low-maintenance, and are hardy. However, these plant species are not always suitable for any location and are often considered invasive and/or weed-like. The present research initiates a list of plant species which can be used within the southeastern United States and similar areas to phytoremediate commonly used pesticides atrazine, azoxystrobin, and imidacloprid and prevent off-target movement.

De novo transcriptome in roots of switchgrass (Panicum virgatum L.) reveals gene expression dynamic and act network under alkaline salt stress

Background

Soil salinization is a major limiting factor for crop cultivation. Switchgrass is a perennial rhizomatous bunchgrass that is considered an ideal plant for marginal lands, including sites with saline soil. Here we investigated the physiological responses and transcriptome changes in the roots of Alamo (alkaline-tolerant genotype) and AM-314/MS-155 (alkaline-sensitive genotype) under alkaline salt stress.

Results

Alkaline salt stress significantly affected the membrane, osmotic adjustment and antioxidant systems in switchgrass roots, and the ASTTI values between Alamo and AM-314/MS-155 were divergent at different time points. A total of 108,319 unigenes were obtained after reassembly, including 73,636 unigenes in AM-314/MS-155 and 65,492 unigenes in Alamo. A total of 10,219 DEGs were identified, and the number of upregulated genes in Alamo was much greater than that in AM-314/MS-155 in both the early and late stages of alkaline salt stress. The DEGs in AM-314/MS-155 were mainly concentrated in the early stage, while Alamo showed greater advantages in the late stage. These DEGs were mainly enriched in plant-pathogen interactions, ubiquitin-mediated proteolysis and glycolysis/gluconeogenesis pathways. We characterized 1480 TF genes into 64 TF families, and the most abundant TF family was the C2H2 family, followed by the bZIP and bHLH families. A total of 1718 PKs were predicted, including CaMK, CDPK, MAPK and RLK. WGCNA revealed that the DEGs in the blue, brown, dark magenta and light steel blue 1 modules were associated with the physiological changes in roots of switchgrass under alkaline salt stress. The consistency between the qRT-PCR and RNA-Seq results confirmed the reliability of the RNA-seq sequencing data. A molecular regulatory network of the switchgrass response to alkaline salt stress was preliminarily constructed on the basis of transcriptional regulation and functional genes.

Conclusions

Alkaline salt tolerance of switchgrass may be achieved by the regulation of ion homeostasis, transport proteins, detoxification, heat shock proteins, dehydration and sugar metabolism. These findings provide a comprehensive analysis of gene expression dynamic and act network induced by alkaline salt stress in two switchgrass genotypes and contribute to the understanding of the alkaline salt tolerance mechanism of switchgrass and the improvement of switchgrass germplasm.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07368-w.

Genetic Transformation of Recalcitrant Upland Switchgrass Using Morphogenic Genes

Switchgrass (Panicum virgatum) is an excellent feedstock for biofuel production. While genetic transformation is routinely done in lowland switchgrass, upland cultivars remain recalcitrant to genetic transformation. Here we report the establishment of an efficient and reproducible transformation protocol for two upland cultivars, 'Summer' and 'Blackwell', by ectopic overexpression of morphogenic genes, Baby boom (Bbm) and Wuschel2 (Wus2). Two auxotrophic Agrobacterium strains, LBA4404Thy- and EHA105Thy-, each harboring the same construct containing ZmBbm, ZmWus2, and a green fluorescence protein (GFP) gene, ZsGreen1, were used to infect immature leaf segments derived from in vitro grown seedlings. The Agrobacterium strains also contain a transformation helper plasmid that carry additional copies of Agrobacterium virulence genes. GFP-expressing calli were identified and selected for regeneration. The highest transformation efficiency of 6% was obtained for the tetraploid cultivar Summer when LBA4404Thy- was used for infection, which is twice of that for the octoploid cultivar Blackwell. LBA4404Thy- consistently outperformed EHA105Thy- on transformation frequency across the two cultivars. Fifteen randomly selected putative transgenic plants of Summer and Blackwell, representing independent callus events, were confirmed as transgenic by the presence of the transgene, ZmAls, and the absence of AtuFtsZ, a chromosomal gene specific to the Agrobacterium strain LBA4404 using polymerase chain reaction. Transgene integration and expression was further confirmed by the detection of GFP in roots, and the resistance to herbicide injury to leaves of selected putative transgenic plants. The ZmBbm and ZmWus2 genes were successfully removed from 40 to 33.3% of the transgenic plants of Summer and Blackwell, respectively, via the Cre-Lox recombination system upon heat treatment of GFP-expressing embryogenic calli. Our successful transformation of recalcitrant upland switchgrass provides a method for gene function analysis and germplasm enhancement via biotechnology.

Linkage mapping evidence for a syntenic QTL associated with flowering time in perennial C4 rhizomatous grasses Miscanthus and switchgrass

Flowering in perennial species is directed via complex signalling pathways that adjust to developmental regulations and environmental cues. Synchronized flowering in certain environments is a prerequisite to commercial seed production, and so the elucidation of the genetic architecture of flowering time in Miscanthus and switchgrass could aid breeding in these underdeveloped species. In this context, we assessed a mapping population in Miscanthus and two ecologically diverse switchgrass mapping populations over 3 years from planting. Multiple flowering time quantitative trait loci (QTL) were identified in both species. Remarkably, the most significant Miscanthus and switchgrass QTL proved to be syntenic, located on linkage groups 4 and 2, with logarithm of odds scores of 17.05 and 21.8 respectively. These QTL regions contained three flowering time transcription factors: Squamosa Promoter-binding protein-Like, MADS-box SEPELLATA2 and gibberellin-responsive bHLH137. The former is emerging as a key component of the age-related flowering time pathway.

Differential Defense Responses of Upland and Lowland Switchgrass Cultivars to a Cereal Aphid Pest

Yellow sugarcane aphid (YSA) (Sipha flava, Forbes) is a damaging pest on many grasses. Switchgrass (Panicum virgatum L.), a perennial C4 grass, has been selected as a bioenergy feedstock because of its perceived resilience to abiotic and biotic stresses. Aphid infestation on switchgrass has the potential to reduce the yields and biomass quantity. Here, the global defense response of switchgrass cultivars Summer and Kanlow to YSA feeding was analyzed by RNA-seq and metabolite analysis at 5, 10, and 15 days after infestation. Genes upregulated by infestation were more common in both cultivars compared to downregulated genes. In total, a higher number of differentially expressed genes (DEGs) were found in the YSA susceptible cultivar (Summer), and fewer DEGs were observed in the YSA resistant cultivar (Kanlow). Interestingly, no downregulated genes were found in common between each time point or between the two switchgrass cultivars. Gene co-expression analysis revealed upregulated genes in Kanlow were associated with functions such as flavonoid, oxidation-response to chemical, or wax composition. Downregulated genes for the cultivar Summer were found in co-expression modules with gene functions related to plant defense mechanisms or cell wall composition. Global analysis of defense networks of the two cultivars uncovered differential mechanisms associated with resistance or susceptibility of switchgrass in response to YSA infestation. Several gene co-expression modules and transcription factors correlated with these differential defense responses. Overall, the YSA-resistant Kanlow plants have an enhanced defense even under aphid uninfested conditions.

Soil precipitation legacies influence intraspecific plant-soil feedback

Feedbacks between plants and soil microbial communities can play an important role in structuring plant communities. However, little is known about how soil legacies caused by environmental disturbances such as drought and extreme precipitation events may affect plant-soil feedback or whether plant-soil feedback operates within species as it does between species. If soil legacies alter plant-soil feedback among genotypes within a plant species, then soil legacies may alter the diversity within plant populations. We conducted a fully factorial pairwise plant-soil feedback experiment to test how precipitation legacies influenced intraspecific plant-soil feedbacks among three genotypes of a dominant grass species, Panicum virgatum. Panicum virgatum experienced negative intraspecific plant-soil feedback, i.e., genotypes generally performed worse on soil from the same genotype than different genotypes. Soil precipitation legacies reversed the rank order of the strength of negative feedback among the genotypes. Feedback is often positively correlated with plant relative abundance. Therefore, our results suggest that soil precipitation legacies may alter the genotypic composition of P. virgatum populations, favoring genotypes that develop less negative feedback. Changes in intraspecific diversity will likely further affect community structure and ecosystem functioning, and may constrain the ability of populations to respond to future changes in climate.

Aphid-Responsive Defense Networks in Hybrid Switchgrass

Aphid herbivory elicits plant defense-related networks that are influenced by host genetics. Plants of the upland switchgrass (Panicum virgatum) cultivar Summer can be a suitable host for greenbug aphids (Schizaphis graminum; GB), and yellow sugarcane aphids (Sipha flava, YSA), whereas the lowland cultivar Kanlow exhibited multi-species resistance that curtails aphid reproduction. However, stabilized hybrids of Summer (♀) x Kanlow (♂) (SxK) with improved agronomics can be damaged by both aphids. Here, hormone and metabolite analyses, coupled with RNA-Seq analysis of plant transcriptomes, were utilized to delineate defense networks induced by aphid feeding in SxK switchgrass and pinpoint plant transcription factors (TFs), such as WRKYs that potentially regulate these responses. Abscisic acid (ABA) levels were significantly higher in GB infested plants at 5 and 10 days after infestation (DAI). ABA levels were highest at 15DAI in YSA infested plants. Jasmonic acid levels were significantly elevated under GB infestation, while salicylic acid levels were signifi40cantly elevated only at 15 DAI in YSA infested plants. Similarly, levels of several metabolites were altered in common or specifically to each aphid. YSA infestation induced a significant enrichment of flavonoids consistent with an upregulation of many genes associated with flavonoid biosynthesis at 15DAI. Gene co-expression modules that responded singly to either aphid or in common to both aphids were differentiated and linked to specific TFs. Together, these data provide important clues into the interplay of metabolism and transcriptional remodeling accompanying defense responses to aphid herbivory in hybrid switchgrass.

Extensions of BLUP Models for Genomic Prediction in Heterogeneous Populations: Application in a Diverse Switchgrass Sample

Genomic prediction is a useful tool to accelerate genetic gain in selection using DNA marker information. However, this technology typically relies on standard prediction procedures, such as genomic BLUP, that are not designed to accommodate population heterogeneity resulting from differences in marker effects across populations. In this study, we assayed different prediction procedures to capture marker-by-population interactions in genomic prediction models. Prediction procedures included genomic BLUP and two kernel-based extensions of genomic BLUP which explicitly accounted for population heterogeneity. To model population heterogeneity, dissemblance between populations was either depicted by a unique coefficient (as previously reported), or a more flexible function of genetic distance between populations (proposed herein). Models under investigation were applied in a diverse switchgrass sample under two validation schemes: whole-sample calibration, where all individuals except selection candidates are included in the calibration set, and cross-population calibration, where the target population is entirely excluded from the calibration set. First, we showed that using fixed effects, from principal components or putative population groups, appeared detrimental to prediction accuracy, especially in cross-population calibration. Then we showed that modeling population heterogeneity by our proposed procedure resulted in highly significant improvements in model fit. In such cases, gains in accuracy were often positive. These results suggest that population heterogeneity may be parsimoniously captured by kernel methods. However, in cases where improvement in model fit by our proposed procedure is null-to-moderate, ignoring heterogeneity should probably be preferred due to the robustness and simplicity of the standard genomic BLUP model.

Breeding progress and preparedness for mass-scale deployment of perennial lignocellulosic biomass crops switchgrass, miscanthus, willow and poplar

Genetic improvement through breeding is one of the key approaches to increasing biomass supply. This paper documents the breeding progress to date for four perennial biomass crops (PBCs) that have high output-input energy ratios: namely Panicum virgatum (switchgrass), species of the genera Miscanthus (miscanthus), Salix (willow) and Populus (poplar). For each crop, we report on the size of germplasm collections, the efforts to date to phenotype and genotype, the diversity available for breeding and on the scale of breeding work as indicated by number of attempted crosses. We also report on the development of faster and more precise breeding using molecular breeding techniques. Poplar is the model tree for genetic studies and is furthest ahead in terms of biological knowledge and genetic resources. Linkage maps, transgenesis and genome editing methods are now being used in commercially focused poplar breeding. These are in development in switchgrass, miscanthus and willow generating large genetic and phenotypic data sets requiring concomitant efforts in informatics to create summaries that can be accessed and used by practical breeders. Cultivars of switchgrass and miscanthus can be seed-based synthetic populations, semihybrids or clones. Willow and poplar cultivars are commercially deployed as clones. At local and regional level, the most advanced cultivars in each crop are at technology readiness levels which could be scaled to planting rates of thousands of hectares per year in about 5 years with existing commercial developers. Investment in further development of better cultivars is subject to current market failure and the long breeding cycles. We conclude that sustained public investment in breeding plays a key role in delivering future mass-scale deployment of PBCs.

A dynamic model of lignin biosynthesis in Brachypodium distachyon

BACKGROUND

Lignin is a crucial molecule for terrestrial plants, as it offers structural support and permits the transport of water over long distances. The hardness of lignin reduces plant digestibility by cattle and sheep; it also makes inedible plant materials recalcitrant toward the enzymatic fermentation of cellulose, which is a potentially valuable substrate for sustainable biofuels. Targeted attempts to change the amount or composition of lignin in relevant plant species have been hampered by the fact that the lignin biosynthetic pathway is difficult to understand, because it uses several enzymes for the same substrates, is regulated in an ill-characterized manner, may operate in different locations within cells, and contains metabolic channels, which the plant may use to funnel initial substrates into specific monolignols.

RESULTS

We propose a dynamic mathematical model that integrates various datasets and other information regarding the lignin pathway in Brachypodium distachyon and permits explanations for some counterintuitive observations. The model predicts the lignin composition and label distribution in a BdPTAL knockdown strain, with results that are quite similar to experimental data.

CONCLUSION

Given the present scarcity of available data, the model resulting from our analysis is presumably not final. However, it offers proof of concept for how one may design integrative pathway models of this type, which are necessary tools for predicting the consequences of genomic or other alterations toward plants with lignin features that are more desirable than in their wild-type counterparts.

Genome-Wide Association Study in Pseudo-F2 Populations of Switchgrass Identifies Genetic Loci Affecting Heading and Anthesis Dates

Switchgrass (Panicum virgatum) is a native prairie grass and valuable bio-energy crop. The physiological change from juvenile to reproductive adult can draw important resources away from growth into producing reproductive structures, thereby limiting the growth potential of early flowering plants. Delaying the flowering of switchgrass is one approach by which to increase total biomass. The objective of this research was to identify genetic variants and candidate genes for controlling heading and anthesis in segregating switchgrass populations. Four pseudo-F₂ populations (two pairs of reciprocal crosses) were developed from lowland (late flowering) and upland (early flowering) ecotypes, and heading and anthesis dates of these populations were collected in Lafayette, IN and DeKalb, IL in 2015 and 2016. Across 2 years, there was a 34- and 73-day difference in heading and a 52- and 75-day difference in anthesis at the Lafayette and DeKalb locations, respectively. A total of 37,901 single nucleotide polymorphisms obtained by exome capture sequencing of the populations were used in a genome-wide association study (GWAS) that identified five significant signals at three loci for heading and two loci for anthesis. Among them, a homolog of FLOWERING LOCUS T on chromosome 5b associated with heading date was identified at the Lafayette location across 2 years. A homolog of ARABIDOPSIS PSEUDO-RESPONSE REGULATOR 5, a light modulator in the circadian clock associated with heading date was detected on chromosome 8a across locations and years. These results demonstrate that genetic variants related to floral development could lend themselves to a long-term goal of developing late flowering varieties of switchgrass with high biomass yield.

Dual Utilization of Medicinal and Aromatic Crops as Bioenergy Feedstocks

Dual production of biofuels and chemicals can increase the economic value of lignocellulosic bioenergy feedstocks. We compared the bioenergy potential of several essential oil (EO) crops with switchgrass ( Panicum virgatum L.), a crop chosen to benchmark biomass and lignocellulosic biofuel production. The EO crops of interest were peppermint ( Mentha × piperita L.), "Scotch" spearmint ( Mentha × gracilis Sole), Japanese cornmint ( Mentha canadensis L.), and sweet sagewort ( Artemisia annua L.). We also assessed each crop for EO production in a marginal production environment in Wyoming, USA, with irrigation and nitrogen (N) rates using a split-plot experimental design. Oil content ranged from 0.31 to 0.4% for Japanese cornmint, 0.23 to 0.26% for peppermint, 0.38 to 0.5% for spearmint, and the overall mean of sweet sagewort was 0.34%. Oil yields ranged from (in kg ha^-1) 34 to 165 in Japanese cornmint, 25 to 108 in peppermint, 29.3 to 126 in spearmint, and 39.7 in sweet sagewort. EO production, but not composition, was sensitive to N fertilization. The alternative bioenergy crops and switchgrass produced similar amounts of ethanol from bench-scale simultaneous saccharification and fermentation assays. Value-added incomes from the EO proceeds were estimated to be between $1055 and $5132 ha^-1 from peppermint, $1309 and $5580 ha^-1 from spearmint, $510 and $2460 ha^-1 from Japanese cornmint, and $3613 ha^-1 from sweet sagewort under Wyoming growth conditions. The advantage of the proposed crops over traditional lignocellulosic species is the production of high-value natural products in addition to lignocellulosic biofuel production.

Impacts of biomass production at civil airports on grassland bird conservation and aviation strike risk

Growing concerns about climate change, foreign oil dependency, and environmental quality have fostered interest in perennial native grasses (e.g., switchgrass [Panicum virgatum]) for bioenergy production while also maintaining biodiversity and ecosystem function. However, biomass cultivation in marginal landscapes such as airport grasslands may have detrimental effects on aviation safety as well as conservation efforts for grassland birds. In 2011-2013, we investigated effects of vegetation composition and harvest frequency on seasonal species richness and habitat use of grassland birds and modeled relative abundance, aviation risk, and conservation value of birds associated with biomass crops. Avian relative abundance was greater in switchgrass monoculture plots during the winter months, whereas Native Warm-Season Grass (NWSG) mixed species plantings were favored by species during the breeding season. Conversely, treatment differences in aviation risk and conservation value were not biologically significant. Only 2.6% of observations included avian species of high hazard to aircraft, providing support for semi-natural grasslands as a feasible landcover option at civil airports. Additionally, varied harvest frequencies across a mosaic of switchgrass monocultures and NWSG plots allows for biomass production with multiple vegetation structure options for grassland birds to increase seasonal avian biodiversity and habitat use.

Mathematical models of lignin biosynthesis

BACKGROUND

Lignin is a natural polymer that is interwoven with cellulose and hemicellulose within plant cell walls. Due to this molecular arrangement, lignin is a major contributor to the recalcitrance of plant materials with respect to the extraction of sugars and their fermentation into ethanol, butanol, and other potential bioenergy crops. The lignin biosynthetic pathway is similar, but not identical in different plant species. It is in each case comprised of a moderate number of enzymatic steps, but its responses to manipulations, such as gene knock-downs, are complicated by the fact that several of the key enzymes are involved in several reaction steps. This feature poses a challenge to bioenergy production, as it renders it difficult to select the most promising combinations of genetic manipulations for the optimization of lignin composition and amount.

RESULTS

Here, we present several computational models than can aid in the analysis of data characterizing lignin biosynthesis. While minimizing technical details, we focus on the questions of what types of data are particularly useful for modeling and what genuine benefits the biofuel researcher may gain from the resulting models. We demonstrate our analysis with mathematical models for black cottonwood (Populus trichocarpa), alfalfa (Medicago truncatula), switchgrass (Panicum virgatum) and the grass Brachypodium distachyon.

CONCLUSIONS

Despite commonality in pathway structure, different plant species show different regulatory features and distinct spatial and topological characteristics. The putative lignin biosynthes pathway is not able to explain the plant specific laboratory data, and the necessity of plant specific modeling should be heeded.

Non-targeted Colonization by the Endomycorrhizal Fungus, Serendipita vermifera, in Three Weeds Typically Co-occurring with Switchgrass

Serendipita vermifera (=Sebacina vermifera; isolate MAFF305830) is a mycorrhizal fungus originally isolated from the roots of an Australian orchid that we have previously shown to be beneficial in enhancing biomass yield and drought tolerance in switchgrass, an important bioenergy crop for cellulosic ethanol production in the United States. However, almost nothing is known about how this root-associated fungus proliferates and grows through the soil matrix. Such information is critical to evaluate the possibility of non-target effects, such as unintended spread to weedy plants growing near a colonized switchgrass plant in a field environment. A microcosm experiment was conducted to study movement of vegetative mycelia of S. vermifera between intentionally inoculated switchgrass (Panicum virgatum L.) and nearby weeds. We constructed size-exclusion microcosms to test three different common weeds, large crabgrass (Digitaria sanguinalis L.), Texas panicum (Panicum texanum L.), and Broadleaf signalgrass (Brachiaria platyphylla L.), all species that typically co-occur in Southern Oklahoma and potentially compete with switchgrass. We report that such colonization of non-target plants by S. vermifera can indeed occur, seemingly via co-mingled root systems. As a consequence of colonization, significant enhancement of growth was noted in signalgrass, while a mild increase (albeit not significant) was evident in crabgrass. Migration of the fungus seems unlikely in root-free bulk soil, as we failed to see transmission when the roots were kept separate. This research is the first documentation of non-targeted colonization of this unique root symbiotic fungus and highlights the need for such assessments prior to deployment of biological organisms in the field.

Analyses of methylomes of upland and lowland switchgrass (Panicum virgatum) ecotypes using MeDIP-seq and BS-seq

Background

Switchgrass is a crop with many desirable traits for bioenergy production. Plant genomes have high DNA methylation levels throughout genes and transposable elements and DNA methylation is known to play a role in silencing transposable elements. Here we analyzed methylomes in two switchgrass genotypes AP13 and VS16. AP13 is derived from a lowland ecotype and VS16, typically considered drought-tolerant, is derived from an upland ecotype, both genotypes are tetraploid (2n = 4× = 36).

Results

Methylated DNA immunoprecipitation-sequencing (MeDIP-seq) and bisulfite-sequencing (BS-seq) were used to profile DNA methylation in genomic features of AP13 and VS16. The methylation patterns in genes and transposable elements were similar to other plants, however, overall CHH methylation levels were comparatively low. Differentially methylated regions (DMRs) were assessed and a total of 1777 CG-DMRs, 573 CHG-DMRs, and 3 CHH-DMRs were detected between the two genotypes. TEs and their flanking regions were higher than that of genic regions. Different types of TEs had different methylation patterns, but the two LTRs (Copia and Gypsy) were similarly methylated, while LINEs and DNA transposons typically had different methylation patterns. MeDIP-seq data was compared to BS-seq data and most of the peaks generated by MeDIP-seq were confirmed to be highly methylated by BS-seq.

Conclusions

DNA methylation in switchgrass genotypes obtained from the two ecotypes were found similar. Collinear gene pairs in two subgenomes (A and B) were not significantly differentially methylated. Both BS-seq and MeDIP-seq methodologies were found effective. Methylation levels were highest at CG and least in CHH. Increased DNA methylation was seen in TEs compared to genic regions. Exploitation of TE methylations can be a viable option in future crop improvement.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4218-0) contains supplementary material, which is available to authorized users.

Crop-associated virus reduces the rooting depth of non-crop perennial native grass more than non-crop-associated virus with known viral suppressor of RNA silencing (VSR)

As agricultural acreage expanded and came to dominate landscapes across the world, viruses gained opportunities to move between crop and wild native plants. In the Midwestern USA, virus exchange currently occurs between widespread annual Poaceae crops and remnant native perennial prairie grasses now under consideration as bioenergy feedstocks. In this region, the common aphid species Rhopalosiphum padi L. (the bird cherry-oat aphid) transmits several virus species in the family Luteoviridae, including Barley yellow dwarf virus (BYDV-PAV, genus Luteovirus) and Cereal yellow dwarf virus (CYDV-RPV and -RPS, genus Polerovirus). The yellow dwarf virus (YDV) species in these two genera share genetic similarities in their 3'-ends, but diverge in the 5'-regions. Most notably, CYDVs encode a P0 viral suppressor of RNA silencing (VSR) absent in BYDV-PAV. Because BYDV-PAV has been reported more frequently in annual cereals and CYDVs in perennial non-crop grasses, we examine the hypothesis that the viruses' genetic differences reflect different affinities for crop and non-crop hosts. Specifically, we ask (i) whether CYDVs might persist within and affect a native non-crop grass more strongly than BYDV-PAV, on the grounds that the polerovirus VSR could better moderate the defenses of a well-defended perennial, and (ii) whether the opposite pattern of effects might occur in a less defended annual crop. Because previous work found that the VSR of CYDV-RPS possessed greater silencing suppressor efficiency than that of CYDV-RPV, we further explored (iii) whether a novel grass-associated CYDV-RPS isolate would influence a native non-crop grass more strongly than a comparable CYDV-RPV isolate. In growth chamber studies, we found support for this hypothesis: only grass-associated CYDV-RPS stunted the shoots and crowns of Panicum virgatum L. (switchgrass), a perennial native North American prairie grass, whereas crop-associated BYDV-PAV (and coinfection with BYDV-PAV and CYDV-RPS) most stunted annual Avena sativa L. (oats). These findings suggest that some of the diversity in grass-infecting Luteoviridae reflects viral capacity to modulate defenses in different host types. Intriguingly, while all virus treatments also reduced root production in both host species, only crop-associated BYDV-PAV (or co-infection) reduced rooting depths. Such root effects may increase host susceptibility to drought, and indicate that BYDV-PAV pathogenicity is determined by something other than a P0 VSR. These findings contribute to growing evidence that pathogenic crop-associated viruses may harm native species as well as crops. Critical next questions include the extent to which crop-associated selection pressures drive viral pathogenesis.

Simultaneous Downregulation of MTHFR and COMT in Switchgrass Affects Plant Performance and Induces Lesion-Mimic Cell Death

Switchgrass (Panicum virgatum) has been developed into a model lignocellulosic bioenergy crop. Downregulation of caffeic acid O-methyltransferase (COMT), a key enzyme in lignin biosynthesis, has been shown to alter lignification and increase biofuel yield in switchgrass. Methylenetetrahydrofolate reductase (MTHFR) mediates C1 metabolism and provides methyl units consumed by COMT. It was predicted that co-silencing of MTHFR and COMT would impact lignification even more than either of the single genes. However, our results showed that strong downregulation of MTHFR in a COMT-deficient background led to altered plant growth and development, but no significant change in lignin content or composition was found when compared with COMT plants. Another unexpected finding was that the double MTHFR/COMT downregulated plants showed a novel lesion-mimic leaf phenotype. Molecular analyses revealed that the lesion-mimic phenotype was caused by the synergistic effect of MTHFR and COMT genes, with MTHFR playing a predominant role. Microarray analysis showed significant induction of genes related to oxidative and defense responses. The results demonstrated the lack of additive effects of MTHFR and COMT on lignification. Furthermore, this research revealed an unexpected role of the two genes in the modulation of lesion-mimic cell death as well as their synergistic effects on agronomic performance.

Growth and fecundity of fertile Miscanthus × giganteus ("PowerCane") compared to feral and ornamental Miscanthus sinensis in a common garden experiment: Implications for invasion

Perennial grasses are promising candidates for bioenergy crops, but species that can escape cultivation and establish self-sustaining naturalized populations (feral) may have the potential to become invasive. Fertile Miscanthus × giganteus, known as "PowerCane," is a new potential biofuel crop. Its parent species are ornamental, non-native Miscanthus species that establish feral populations and are sometimes invasive in the USA. As a first step toward assessing the potential for "PowerCane" to become invasive, we documented its growth and fecundity relative to one of its parent species (Miscanthus sinensis) in competition with native and invasive grasses in common garden experiments located in Columbus, Ohio and Ames, Iowa, within the targeted range of biofuel cultivation. We conducted a 2-year experiment to compare growth and reproduction among three Miscanthus biotypes-"PowerCane," ornamental M. sinensis, and feral M. sinensis-at two locations. Single Miscanthus plants were subjected to competition with a native grass (Panicum virgatum), a weedy grass (Bromus inermis), or no competition. Response variables were aboveground biomass, number of shoots, basal area, and seed set. In Iowa, all Miscanthus plants died after the first winter, which was unusually cold, so no further results are reported from the Iowa site. In Ohio, we found significant differences among biotypes in growth and fecundity, as well as significant effects of competition. Interactions between these treatments were not significant. "PowerCane" performed as well or better than ornamental or feral M. sinensis in vegetative traits, but had much lower seed production, perhaps due to pollen limitation. In general, ornamental M. sinensis performed somewhat better than feral M. sinensis. Our findings suggest that feral populations of "PowerCane" could become established adjacent to biofuel production areas. Fertile Miscanthus × giganteus should be studied further to assess its potential to spread via seed production in large, sexually compatible populations.

Improved tissue culture conditions for the emerging C4 model Panicum hallii

BACKGROUND

Panicum hallii Vasey (Hall's panicgrass) is a compact, perennial C4 grass in the family Poaceae, which has potential to enable bioenergy research for switchgrass (Panicum virgatum L.). Unlike P. hallii, switchgrass has a large genome, allopolyploidy, self-incompatibility, a long life cycle, and large stature-all suboptimal traits for rapid genetics research. Herein we improved tissue culture methodologies for two inbred P. hallii populations: FIL2 and HAL2, to enable further development of P. hallii as a model C4 plant.

RESULTS

The optimal seed-derived callus induction medium was determined to be Murashige and Skoog (MS) medium supplemented with 40 mg L-1 L-cysteine, 300 mg L-1 L-proline, 3% sucrose, 1 g L-1 casein hydrolysate, 3 mg L-1 2,4-dichlorophenoxyacetic acid (2,4-D), and 45 μg L-1 6-benzylaminopurine (BAP), which resulted in callus induction of 51 ± 29% for FIL2 and 81 ± 19% for HAL2. The optimal inflorescence-derived callus induction was observed on MP medium (MS medium supplemented with 2 g L-1 L-proline, 3% maltose, 5 mg L-1 2,4-D, and 500 μg L-1 BAP), resulting in callus induction of 100 ± 0.0% for FIL2 and 84 ± 2.4% for HAL2. Shoot regeneration rates of 11.5 ± 0.8 shoots/gram for FIL2 and 11.3 ± 0.6 shoots/gram for HAL2 were achieved using seed-induced callus, whereas shoot regeneration rates of 26.2 ± 2.6 shoots/gram for FIL2 and 29.3 ± 3.6 shoots/gram for HAL2 were achieved from inflorescence-induced callus. Further, cell suspension cultures of P. hallii were established from seed-derived callus, providing faster generation of callus tissue compared with culture using solidified media (1.41-fold increase for FIL2 and 3.00-fold increase for HAL2).

CONCLUSIONS

Aside from abbreviated tissue culture times from callus induction to plant regeneration for HAL2, we noted no apparent differences between FIL2 and HAL2 populations in tissue culture performance. For both populations, the cell suspension cultures outperformed tissue cultures on solidified media. Using the methods developed in this work, P. hallii callus was induced from seeds immediately after harvest in a shorter time and with higher frequencies than switchgrass. For clonal propagation, P. hallii callus was established from R1 inflorescences, similar to switchgrass, which further strengthens the potential of this plant as a C4 model for genetic studies. The rapid cycling (seed-to-seed time) and ease of culture, further demonstrate the potential utility of P. hallii as a C4 model plant.

Defined tetra-allelic gene disruption of the 4-coumarate:coenzyme A ligase 1 (Pv4CL1) gene by CRISPR/Cas9 in switchgrass results in lignin reduction and improved sugar release

BACKGROUND

The development of genome editing technologies offers new prospects in improving bioenergy crops like switchgrass (Panicum virgatum). Switchgrass is an outcrossing species with an allotetraploid genome (2n = 4x = 36), a complexity which forms an impediment to generating homozygous knock-out plants. Lignin, a major component of the plant cell wall and a contributor to cellulosic feedstock's recalcitrance to decomposition, stands as a barrier to efficient biofuel production by limiting enzyme access to cell wall polymers during the fermentation process.

RESULTS

We developed a CRISPR/Cas9 genome editing system in switchgrass to target a key enzyme involved in the early steps of monolignol biosynthesis, 4-Coumarate:coenzyme A ligase (4CL). Three 4CL genes, Pv4CL1, Pv4CL2, and Pv4CL3, were identified in switchgrass. Expression analysis revealed that Pv4CL1 transcripts were more abundant in the stem than in the leaf, while Pv4CL2 transcripts were barely detectable and Pv4CL3 was mainly expressed in the leaf. Pv4CL1 was selected as the target for CRISPR/Cas9 editing because of its preferential expression in highly lignified stem tissues. Specific guide RNA was constructed to target Pv4CL1. After introducing the construct into switchgrass calli, 39 transgenic plants were regenerated. Using two rounds of PCR screening and sequencing, four plants were confirmed to have tetra-allelic mutations simultaneously. The Pv4CL1 knock-out plants had reduced cell wall thickness, an 8-30% reduction in total lignin content, a 7-11% increase in glucose release, and a 23-32% increase in xylose release.

CONCLUSION

This study established a successful CRISPR/Cas9 system in switchgrass with mutation efficiency reaching 10%. The system allows the precise targeting of the selected Pv4CL1 gene to create switchgrass knock-out mutant plants with decreased lignin content and reduced recalcitrance.

Switchgrass (Panicum virgatum) Intercropping within Managed Loblolly Pine (Pinus taeda) Does Not Affect Wild Bee Communities

Intensively-managed pine (Pinus spp.) have been shown to support diverse vertebrate communities, but their ability to support invertebrate communities, such as wild bees, has not been well-studied. Recently, researchers have examined intercropping switchgrass (Panicum virgatum), a native perennial, within intensively managed loblolly pine (P. taeda) plantations as a potential source for cellulosic biofuels. To better understand potential effects of intercropping on bee communities, we investigated visitation of bees within three replicates of four treatments of loblolly pine in Mississippi, U.S.A.: 3-4 year old pine plantations and 9-10 year old pine plantations with and without intercropped switchgrass. We used colored pan traps to capture bees during the growing seasons of 2013 and 2014. We captured 2507 bees comprised of 18 different genera during the two-year study, with Lasioglossum and Ceratina being the most common genera captured. Overall, bee abundances were dependent on plantation age and not presence of intercropping. Our data suggests that switchgrass does not negatively impact or promote bee communities within intensively-managed loblolly pine plantations.

Detection and phylogenetic relationships of Puccinia emaculata and Uromyces graminicola (Pucciniales) on switchgrass in New York State using rDNA sequence information

The species of rust fungi (Pucciniales) inciting disease on switchgrass (Panicum virgatum) grown in bioenergy feedstock systems across the north-central and eastern United States remain unclear. In the present study, the species number and phylogenetic relationships of rust species affecting switchgrass were examined in 2011-2013 at two sites in New York State as well as selected sites in Alabama, Iowa, Nebraska, Pennsylvania, South Dakota, and West Virginia using ribosomal RNA gene data (partial internal transcribed spacer [ITS] 1, complete 5.8 subunit [S] and ITS2, and partial 28S). Uredinial group and teliospore morphology were also utilized to delimit taxa in collection years 2012 and 2013. Maximum likelihood, maximum parsimony, and Bayesian analyses demonstrated two monophyletic clades. Clade I consisted of Puccinia emaculata and included the majority of single-sorus samples across sites, whereas, Clade II included multiple samples from Iowa, Nebraska, and South Dakota. Single-telial samples for Clade I possessed only two-celled teliospores while Clade II samples possessed only one-celled teliospores, and hence, were readily diagnosed morphologically to P. emaculata and Uromyces graminicola, respectively. No U. graminicola sequences exist in GenBank to compare with our Clade II samples; however, based on teliospore morphology, the identity of Clade II taxa is U. graminicola.

Diversity and population structure of northern switchgrass as revealed through exome capture sequencing

Panicum virgatum L. (switchgrass) is a polyploid, perennial grass species that is native to North America, and is being developed as a future biofuel feedstock crop. Switchgrass is present primarily in two ecotypes: a northern upland ecotype, composed of tetraploid and octoploid accessions, and a southern lowland ecotype, composed of primarily tetraploid accessions. We employed high-coverage exome capture sequencing (~2.4 Tb) to genotype 537 individuals from 45 upland and 21 lowland populations. From these data, we identified ~27 million single-nucleotide polymorphisms (SNPs), of which 1 590 653 high-confidence SNPs were used in downstream analyses of diversity within and between the populations. From the 66 populations, we identified five primary population groups within the upland and lowland ecotypes, a result that was further supported through genetic distance analysis. We identified conserved, ecotype-restricted, non-synonymous SNPs that are predicted to affect the protein function of CONSTANS (CO) and EARLY HEADING DATE 1 (EHD1), key genes involved in flowering, which may contribute to the phenotypic differences between the two ecotypes. We also identified, relative to the near-reference Kanlow population, 17 228 genes present in more copies than in the reference genome (up-CNVs), 112 630 genes present in fewer copies than in the reference genome (down-CNVs) and 14 430 presence/absence variants (PAVs), affecting a total of 9979 genes, including two upland-specific CNV clusters. In total, 45 719 genes were affected by an SNP, CNV, or PAV across the panel, providing a firm foundation to identify functional variation associated with phenotypic traits of interest for biofuel feedstock production.

Computational inference of the structure and regulation of the lignin pathway in Panicum virgatum

BACKGROUND

Switchgrass is a prime target for biofuel production from inedible plant parts and has been the subject of numerous investigations in recent years. Yet, one of the main obstacles to effective biofuel production remains to be the major problem of recalcitrance. Recalcitrance emerges in part from the 3-D structure of lignin as a polymer in the secondary cell wall. Lignin limits accessibility of the sugars in the cellulose and hemicellulose polymers to enzymes and ultimately decreases ethanol yield. Monolignols, the building blocks of lignin polymers, are synthesized in the cytosol and translocated to the plant cell wall, where they undergo polymerization. The biosynthetic pathway leading to monolignols in switchgrass is not completely known, and difficulties associated with in vivo measurements of these intermediates pose a challenge for a true understanding of the functioning of the pathway.

RESULTS

In this study, a systems biological modeling approach is used to address this challenge and to elucidate the structure and regulation of the lignin pathway through a computational characterization of alternate candidate topologies. The analysis is based on experimental data characterizing stem and tiller tissue of four transgenic lines (knock-downs of genes coding for key enzymes in the pathway) as well as wild-type switchgrass plants. These data consist of the observed content and composition of monolignols. The possibility of a G-lignin specific metabolic channel associated with the production and degradation of coniferaldehyde is examined, and the results support previous findings from another plant species. The computational analysis suggests regulatory mechanisms of product inhibition and enzyme competition, which are well known in biochemistry, but so far had not been reported in switchgrass. By including these mechanisms, the pathway model is able to represent all observations.

CONCLUSIONS

The results show that the presence of the coniferaldehyde channel is necessary and that product inhibition and competition over cinnamoyl-CoA-reductase (CCR1) are essential for matching the model to observed increases in H-lignin levels in 4-coumarate:CoA-ligase (4CL) knockdowns. Moreover, competition for 4-coumarate:CoA-ligase (4CL) is essential for matching the model to observed increases in the pathway metabolites in caffeic acid O-methyltransferase (COMT) knockdowns. As far as possible, the model was validated with independent data.

C4 bioenergy crops for cool climates, with special emphasis on perennial C4 grasses

There is much interest in cultivating C4 perennial plants in northern climates where there is an abundance of land and a potential large market for biofuels. C4 feedstocks can exhibit superior yields to C3 alternatives during the long warm days of summer at high latitude, but their summer success depends on an ability to tolerate deep winter cold, spring frosts, and early growth-season chill. Here, we review cold tolerance limits in C4 perennial grasses. Dozens of C4 species are known from high latitudes to 63 °N and elevations up to 5200 m, demonstrating that C4 plants can adapt to cold climates. Of the three leading C4 grasses being considered for bioenergy production in cold climates--Miscanthus spp., switchgrass (Panicum virgatum), and prairie cordgrass (Spartina pectinata)--all are tolerant of cool temperatures (10-15 °C), but only cordgrass tolerates hard spring frosts. All three species overwinter as dormant rhizomes. In the productive Miscanthus×giganteus hybrids, exposure to temperatures below -3 °C to -7 °C will kill overwintering rhizomes, while for upland switchgrass and cordgrass, rhizomes survive exposure to temperatures above -20 °C to -24 °C. Cordgrass emerges earlier than switchgrass and M. giganteus genotypes, but lacks the Miscanthus growth potential once warmer days of late spring arrive. To enable C4-based bioenergy production in colder climates, breeding priorities should emphasize improved cold tolerance of M.×giganteus, and enhanced productivity of switchgrass and cordgrass. This should be feasible in the near future, because wild populations of each species exhibit a diverse range of cold tolerance and growth capabilities.

Assessment of drought tolerance of 49 switchgrass (Panicum virgatum) genotypes using physiological and morphological parameters

BACKGROUND

Switchgrass (Panicum virgatum L.) is a warm-season C4 grass that is a target lignocellulosic biofuel species. In many regions, drought stress is one of the major limiting factors for switchgrass growth. The objective of this study was to evaluate the drought tolerance of 49 switchgrass genotypes. The relative drought stress tolerance was determined based on a set of parameters including plant height, leaf length, leaf width, leaf sheath length, leaf relative water content (RWC), electrolyte leakage (EL), photosynthetic rate (Pn), stomatal conductance (g s), transpiration rate (Tr), intercellular CO2 concentration (Ci), and water use efficiency (WUE).

RESULTS

SRAP marker analysis determined that the selected 49 switchgrass genotypes represent a diverse genetic pool of switchgrass germplasm. Principal component analysis (PCA) and drought stress indexes (DSI) of each physiological parameter showed significant differences in the drought stress tolerance among the 49 genotypes. Heatmap and PCA data revealed that physiological parameters are more sensitive than morphological parameters in distinguishing the control and drought treatments. Metabolite profiling data found that under drought stress, the five best drought-tolerant genotypes tended to have higher levels of abscisic acid (ABA), spermine, trehalose, and fructose in comparison to the five most drought-sensitive genotypes.

CONCLUSION

Based on PCA ranking value, the genotypes TEM-SEC, TEM-LoDorm, BN-13645-64, Alamo, BN-10860-61, BN-12323-69, TEM-SLC, T-2086, T-2100, T-2101, Caddo, and Blackwell-1 had relatively higher ranking values, indicating that they are more tolerant to drought. In contrast, the genotypes Grif Nebraska 28, Grenville-2, Central Iowa Germplasm, Cave-in-Rock, Dacotah, and Nebraska 28 were found to be relatively sensitive to drought stress. By analyzing physiological response parameters and different metabolic profiles, the methods utilized in this study identified drought-tolerant and drought-sensitive switchgrass genotypes. These results provide a foundation for future research directed at understanding the molecular mechanisms underlying switchgrass tolerance to drought.

Nucleotide polymorphism and copy number variant detection using exome capture and next-generation sequencing in the polyploid grass Panicum virgatum

Switchgrass (Panicum virgatum) is a polyploid, outcrossing grass species native to North America and has recently been recognized as a potential biofuel feedstock crop. Significant phenotypic variation including ploidy is present across the two primary ecotypes of switchgrass, referred to as upland and lowland switchgrass. The tetraploid switchgrass genome is approximately 1400 Mbp, split between two subgenomes, with significant repetitive sequence content limiting the efficiency of re-sequencing approaches for determining genome diversity. To characterize genetic diversity in upland and lowland switchgrass as a first step in linking genotype to phenotype, we designed an exome capture probe set based on transcript assemblies that represent approximately 50 Mb of annotated switchgrass exome sequences. We then evaluated and optimized the probe set using solid phase comparative genome hybridization and liquid phase exome capture followed by next-generation sequencing. Using the optimized probe set, we assessed variation in the exomes of eight switchgrass genotypes representing tetraploid lowland and octoploid upland cultivars to benchmark our exome capture probe set design. We identified ample variation in the switchgrass genome including 1,395,501 single nucleotide polymorphisms (SNPs), 8173 putative copy number variants and 3336 presence/absence variants. While the majority of the SNPs (84%) detected was bi-allelic, a substantial number was tri-allelic with limited occurrence of tetra-allelic polymorphisms consistent with the heterozygous and polyploid nature of the switchgrass genome. Collectively, these data demonstrate the efficacy of exome capture for discovery of genome variation in a polyploid species with a large, repetitive and heterozygous genome.

Global changes in mineral transporters in tetraploid switchgrasses (Panicum virgatum L.)

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.

Abolishing activity against ascorbate in a cytosolic ascorbate peroxidase from switchgrass

Switchgrass (Panicum virgatum L.) is being developed as a bioenergy species. Recently an early version of its genome has been released permitting a route to the cloning and analysis of key proteins. Ascorbate peroxidases (APx) are an important part of the antioxidant defense system of plant cells and present a well studied model to understand structure-function relationships. Analysis of the genome indicates that switchgrass encodes several cytosolic ascorbate peroxidases with apparent varying levels of tissue expression. A major cytosolic ascorbate peroxidase was thus selected for further studies. This gene was cloned and expressed in Escherichia coli cells to obtain purified active protein. Full heme incorporation of the enzyme was achieved utilizing slow growth and supplementing the media with 5-aminolevulinic acid. The enzyme was observed to be monomeric in solution via size exclusion chromatography. Activity toward ascorbate was observed that was non-Michaelis-Menten in nature. A site-directed mutant, R172S, was made in an attempt to differentiate activity against ascorbate versus other substrates. The R172S protein exhibited negligible ascorbate peroxidase activity, but showed near wild type activity toward other aromatic substrates.

Metagenomes of tropical soil-derived anaerobic switchgrass-adapted consortia with and without iron

Tropical forest soils decompose litter rapidly with frequent episodes of anoxia, making it likely that bacteria using alternate terminal electron acceptors (TEAs) such as iron play a large role in supporting decomposition under these conditions. The prevalence of many types of metabolism in litter deconstruction makes these soils useful templates for improving biofuel production. To investigate how iron availability affects decomposition, we cultivated feedstock-adapted consortia (FACs) derived from iron-rich tropical forest soils accustomed to experiencing frequent episodes of anaerobic conditions and frequently fluctuating redox. One consortium was propagated under fermenting conditions, with switchgrass as the sole carbon source in minimal media (SG only FACs), and the other consortium was treated the same way but received poorly crystalline iron as an additional terminal electron acceptor (SG + Fe FACs). We sequenced the metagenomes of both consortia to a depth of about 150 Mb each, resulting in a coverage of 26× for the more diverse SG + Fe FACs, and 81× for the relatively less diverse SG only FACs. Both consortia were able to quickly grow on switchgrass, and the iron-amended consortium exhibited significantly higher microbial diversity than the unamended consortium. We found evidence of higher stress in the unamended FACs and increased sugar transport and utilization in the iron-amended FACs. This work provides metagenomic evidence that supplementation of alternative TEAs may improve feedstock deconstruction in biofuel production.

Effects of switchgrass (Panicum virgatum) rotations with peanut (Arachis hypogaea L.) on nematode populations and soil microflora

A 3-year field rotation study was conducted to assess the potential of switchgrass (Panicum virgatum) to suppress root-knot nematodes (Meloidogyne arenaria), southern blight (Sclerotium rolfsii), and aflatoxigenic fungi (Aspergillus sp.) in peanut (Arachis hypogaea L.) and to assess shifts in microbial populations following crop rotation. Switchgrass did not support populations of root-knot nematodes but supported high populations of nonparasitic nematodes. Peanut with no nematicide applied and following 2 years of switchgrass had the same nematode populations as continuous peanut plus nematicide. Neither previous crop nor nematicide significantly reduced the incidence of pods infected with Aspergillus. However, pod invasion by A. flavus was highest in plots previously planted with peanut and not treated with nematicide. Peanut with nematicide applied at planting following 2 years of switchgrass had significantly less incidence of southern blight than either continuous peanut without nematicide application or peanut without nematicide following 2 years of cotton. Peanut yield did not differ among rotations in either sample year. Effects of crop rotation on the microbial community structure associated with peanut were examined using indices for diversity, richness, and similarity derived from culture-based analyses. Continuous peanut supported a distinctly different rhizosphere bacterial microflora compared to peanut following 1 year of switchgrass, or continuous switchgrass. Richness and diversity indices for continuous peanut rhizosphere and geocarposphere were not consistently different from peanut following switchgrass, but always differed in the specific genera present. These shifts in community structure were associated with changes in parasitic nematode populations.

Effects of bison grazing on Andropogon gerardii and Panicum virgatum in burned and unburned tallgrass paririe

Responses to clipping and bison grazing in different environmental contexts were examined in two perennial grass species, Andropogon gerardii and Panicum virgatum, on the Konza Prairie in northeastern Kansas. Grazed tillers had lower relative growth rates (RGR) than clipped tillers following defoliation but this difference was transient and final biomass was not affected by mode of defoliation. Grazed tillers of both species had higher RGR throughout the season than ungrazed tillers, resulting in exact compensation for tissue lost to defoliation. However, A. gerardii tillers which had been grazed repeatedly the previous year (1988) had reduced relative growth rates, tiller biomass and tiller survival in 1989. This suggests that the short-term increase in aboveground relative growth rates after defoliation had a cost to future plant growth and tiller survival.In general, the two species had similar responses to defoliation but their responses were altered differentially by fire. The increase in RGR following defoliation of A. gerardii was relatively greater on unburned than burned prairie, and was influenced by topographic position. P. virgatum responses to defoliation were similar in burned and unburned prairie. Thus grazing, fire, and topographical position all interact to influence tiller growth dynamics and these two species respond differently to the fire and grazing interaction. In addition, fire may interact with grazing pattern to influence a plants' grazing history and thus its long-term performance.