Identification, characterization, and gene expression analysis of nucleotide binding site (NB)-type resistance gene homologues in switchgrass

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.

Genome-wide identification and analysis of the NLR gene family in Medicago ruthenica

Medicago ruthenica, important forage in the legume family, possesses high nutritional value and carries abundant tolerance genes. This study used whole-genome data of M. ruthenica to perform a genome-wide analysis of the nucleotide-binding site-leucine-rich repeat receptor (NLR) gene family, which is the largest family of plant disease resistance genes (R genes). A total of 338 NLR genes were identified in the M. ruthenica genome, including 160 typical genes that contained 80 coiled-coil (CC)-NBS-LRR (CNL) genes, 76 toll/interleukin-1 receptor (TIR)-NBS-LRR (TNL) genes, four resistance to powdery mildew 8 (RPW8)-NBS-LRR (RNL) subclass genes, and 178 atypical NLR genes encoding proteins without at least one important domain. Among its eight chromosomes, M. ruthenica chromosomes 3 and 8 contained most of the NLR genes. More than 40% of all NLR genes were located on these two chromosomes, mainly in multigene clusters. The NLR proteins of M. ruthenica had six highly conserved motifs: P-loop, GLPL, RNBS-D, kinase-2, RNBS-C, and MHDV. Phylogenetic analysis revealed that the NLR genes of M. ruthenica formed three deeply separated clades according to the N-terminal domain of the proteins encoded by these genes. Gene duplication and syntenic analysis suggested four gene duplication types in the NLR genes of M. ruthenica, namely, tandem, proximal, dispersed, and segmental duplicates, which involved 189, 49, 59, and 41 genes, respectively. A total of 41 segmental duplication genes formed 23 NLR gene pairs located on syntenic chromosomal blocks mainly between chromosomes 6 and 7. In addition, syntenic analysis between M. truncatula and M. ruthenica revealed 193 gene pairs located on syntenic chromosomal blocks of the two species. The expression analysis of M. ruthenica NLR genes showed that 303 (89.6%) of the NLR genes were expressed in different varieties. Overall, this study described the full NLR profile of the M. ruthenica genome to provide an important resource for mining disease-resistant genes and disease-resistant breeding.

Effect of cultivar and temperature on the synergistic interaction between panicum mosaic virus and satellite panicum mosaic virus in switchgrass

Panicum mosaic virus (PMV), the type member of the genus Panicovirus in the family Tombusviridae, naturally infects switchgrass (Panicum virgatum L.). PMV and its molecular partner, satellite panicum mosaic virus (SPMV), interact synergistically in coinfected millets to exacerbate the disease phenotype and increase the accumulation of PMV compared to plants infected with PMV alone. In this study, we examined the reaction of switchgrass cvs. Summer and Kanlow to PMV and PMV+SPMV infections at 24°C and 32°C. Switchgrass cv. Summer was susceptible to PMV at both temperatures. In contrast, cv. Kanlow was tolerant to PMV at 24°C, but not at 32°C, suggesting that Kanlow harbors temperature-sensitive resistance to PMV. At 24°C, PMV was readily detected in inoculated leaves, but not in upper uninoculated leaves of Kanlow, suggesting that resistance to PMV was likely mediated by abrogation of long-distance virus transport. Coinfection by PMV and SPMV at 24°C and 32°C in cv. Summer, but not in Kanlow, caused increased symptomatic systemic infection and mild disease synergism with slightly increased PMV accumulation compared to plants infected with PMV alone. These data suggest that the interaction between PMV and SPMV in switchgrass is cultivar-dependent, manifested in Summer but not in Kanlow. However, co-inoculation of cv. Kanlow with PMV+SPMV caused an enhanced asymptomatic infection, suggesting a role of SPMV in enhancement of symptomless infection in a tolerant cultivar. These data suggest that enhanced asymptomatic infections in a virus-tolerant switchgrass cultivar could serve as a source of virus spread and play an important role in panicum mosaic disease epidemiology under field conditions. Our data reveal that the cultivar, coinfection with SPMV, and temperature influence the severity of symptoms elicited by PMV in switchgrass.

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.

Characterization of Nucleotide Binding -Site-Encoding Genes in Sweetpotato, Ipomoea batatas(L.) Lam., and Their Response to Biotic and Abiotic Stresses

Sweetpotato, Ipomoea batatas (L.) Lam., is an important and widely grown crop, yet its production is affected severely by biotic and abiotic stresses. The nucleotide binding site (NBS)-encoding genes have been shown to improve stress tolerance in several plant species. However, the characterization of NBS-encoding genes in sweetpotato is not well-documented to date. In this study, a comprehensive analysis of NBS-encoding genes has been conducted on this species by using bioinformatics and molecular biology methods. A total of 315 NBS-encoding genes were identified, and 260 of them contained all essential conserved domains while 55 genes were truncated. Based on domain architectures, the 260 NBS-encoding genes were grouped into 6 distinct categories. Phylogenetic analysis grouped these genes into 3 classes: TIR, CC (I), and CC (II). Chromosome location analysis revealed that the distribution of NBS-encoding genes in chromosomes was uneven, with a number ranging from 1 to 34. Multiple stress-related regulatory elements were detected in the promoters, and the NBS-encoding genes' expression profiles under biotic and abiotic stresses were obtained. According to the bioinformatics analysis, 9 genes were selected for RT-qPCR analysis. The results revealed that IbNBS75, IbNBS219, and IbNBS256 respond to stem nematode infection; Ib-NBS240, IbNBS90, and IbNBS80 respond to cold stress, while IbNBS208, IbNBS71, and IbNBS159 respond to 30% PEG treatment. We hope these results will provide new insights into the evolution of NBS-encoding genes in the sweetpotato genome and contribute to the molecular breeding of sweetpotato in the future.

Genome-wide analysis of the abiotic stress-related bZIP family in switchgrass

The large basic leucine zipper (bZIP) transcription factor family is conserved in plants. These proteins regulate growth, development, and stress response. Here, we conducted a genome-wide analysis to identify the bZIP genes associated with stress resistance in switchgrass (Panicum virgatum L.). We identified 178 PvbZIPs unevenly distributed on 18 switchgrass chromosomes. An evolutionary analysis segregated them into 10 subfamilies. Gene structure and conserved motif analyses indicated that the same subfamily members shared similar intron-exon modes and motif compositions. This finding corroborated the proposed PvbZIP family grouping. A promoter analysis showed that PvbZIP genes participate in various stress responses. Phylogenetic and synteny analyses characterized 111 switchgrass bZIPs as orthologs of 70 rice bZIPs. A protein interaction network analysis revealed that 22 proteins are involved in salt and drought tolerance. An expression atlas disclosed that the expression patterns of several PvbZIPs differ among various tissues and developmental stages. Online data demonstrated that 16 PvbZIPs were significantly downregulated and five were significantly upregulated in response to heat stress. Other PvbZIPs participated in responses to abiotic stress such as salt, drought, cold, and heat. Our genome-wide analysis and identification of the switchgrass bZIP family characterized multiple candidate PvbZIPs that regulate growth and stress response. This study lays theoretical and empirical foundations for future functional investigations into other transcription factors.

Transcriptome divergence during leaf development in two contrasting switchgrass (Panicum virgatum L.) cultivars

The genetics and responses to biotic stressors of tetraploid switchgrass (Panicum virgatum L.) lowland cultivar 'Kanlow' and upland cultivar Summer are distinct and can be exploited for trait improvement. In general, there is a paucity of data on the basal differences in transcription across tissue developmental times for switchgrass cultivars. Here, the changes in basal and temporal expression of genes related to leaf functions were evaluated for greenhouse grown 'Kanlow', and 'Summer' plants. Three biological replicates of the 4th leaf pooled from 15 plants per replicate were harvested at regular intervals beginning from leaf emergence through senescence. Increases and decreases in leaf chlorophyll and N content were similar for both cultivars. Likewise, multidimensional scaling (MDS) analysis indicated both cultivar-independent and cultivar-specific gene expression. Cultivar-independent genes and gene-networks included those associated with leaf function, such as growth/senescence, carbon/nitrogen assimilation, photosynthesis, chlorophyll biosynthesis, and chlorophyll degradation. However, many genes encoding nucleotide-binding leucine rich repeat (NB-LRRs) proteins and wall-bound kinases associated with detecting and responding to environmental signals were differentially expressed. Several of these belonged to unique cultivar-specific gene co-expression networks. Analysis of genomic resequencing data provided several examples of NB-LRRs genes that were not expressed and/or apparently absent in the genomes of Summer plants. It is plausible that cultivar (ecotype)-specific genes and gene-networks could be one of the drivers for the documented differences in responses to leaf-borne pathogens between these two cultivars. Incorporating broad resistance to plant pathogens in elite switchgrass germplasm could improve sustainability of biomass production under low-input conditions.

Fall armyworm (Spodoptera frugiperda Smith) feeding elicits differential defense responses in upland and lowland switchgrass

Switchgrass (Panicum virgatum L.) is a low input, high biomass perennial grass being developed for the bioenergy sector. Upland and lowland cultivars can differ in their responses to insect herbivory. Fall armyworm [FAW; Spodoptera frugiperda JE Smith (Lepidoptera: Noctuidae)] is a generalist pest of many plant species and can feed on switchgrass as well. Here, in two different trials, FAW larval mass were significantly reduced when fed on lowland cultivar Kanlow relative to larvae fed on upland cultivar Summer plants after 10 days. Hormone content of plants indicated elevated levels of the plant defense hormone jasmonic acid (JA) and its bioactive conjugate JA-Ile although significant differences were not observed. Conversely, the precursor to JA, 12-oxo-phytodienoic acid (OPDA) levels were significantly different between FAW fed Summer and Kanlow plants raising the possibility of differential signaling by OPDA in the two cultivars. Global transcriptome analysis revealed a stronger response in Kanlow plant relative to Summer plants. Among these changes were a preferential upregulation of several branches of terpenoid and phenylpropanoid biosynthesis in Kanlow plants suggesting that enhanced biosynthesis or accumulation of antifeedants could have negatively impacted FAW larval mass gain on Kanlow plants relative to Summer plants. A comparison of the switchgrass-FAW RNA-Seq dataset to those from maize-FAW and switchgrass-aphid interactions revealed that key components of plant responses to herbivory, including induction of JA biosynthesis, key transcription factors and JA-inducible genes were apparently conserved in switchgrass and maize. In addition, these data affirm earlier studies with FAW and aphids that the cultivar Kanlow can provide useful genetics for the breeding of switchgrass germplasm with improved insect resistance.

Genome-wide characterization revealed role of NBS-LRR genes during powdery mildew infection in Vitis vinifera

NBS-LRR comprises a large class of disease resistance (R) proteins that play a widespread role in plant protection against pathogens. In grapevine, powdery mildew cause significant losses in its productivity and efforts are being directed towards finding of resistance loci or genes imparting resistance/tolerance against such fungal diseases. In the present study, we performed genome-wide analysis of NBS-LRR genes during PM infection in grapevine. We identified 18, 23, 12, 16, 10, 10, 9, 20 and 14 differentially expressed NBS-LRR genes in response to PM infection in seven partially PM-resistant (DVIT3351.27, Husseine, Karadzhandal, Khalchili, Late vavilov, O34-16, Sochal) and 2 PM-susceptible (Carignan and Thompson seedless) V. vinifera accessions. Further, the identified sequences were characterized based on chromosomal locations, physicochemical properties, gene structure and motif analysis, and functional annotation by Gene Ontology (GO) mapping. The NBS-LRR genes responsive to powdery mildew could potentially be exploited to improve resistance in grapes.

Genome-Wide Identification of NBS-Encoding Resistance Genes in Sunflower (Helianthus annuus L.)

Nucleotide Binding Site—Leucine-Rich Repeat (NBS-LRR) genes encode disease resistance proteins involved in plants’ defense against their pathogens. Although sunflower is affected by many diseases, only a few molecular details have been uncovered regarding pathogenesis and resistance mechanisms. Recent availability of sunflower whole genome sequences in publicly accessible databases allowed us to accomplish a genome-wide identification of Toll-interleukin-1 receptor-like Nucleotide-binding site Leucine-rich repeat (TNL), Coiled Coil (CC)-NBS-LRR (CNL), Resistance to powdery mildew8 (RPW8)-NBS-LRR (RNL) and NBS-LRR (NL) protein encoding genes. Hidden Markov Model (HMM) profiling of 52,243 putative protein sequences from sunflower resulted in 352 NBS-encoding genes, among which 100 genes belong to CNL group including 64 genes with RX_CC like domain, 77 to TNL, 13 to RNL, and 162 belong to NL group. We also identified signal peptides and nuclear localization signals present in the identified genes and their homologs. We found that NBS genes were located on all chromosomes and formed 75 gene clusters, one-third of which were located on chromosome 13. Phylogenetic analyses between sunflower and Arabidopsis NBS genes revealed a clade-specific nesting pattern in CNLs, with RNLs nested in the CNL-A clade, and species-specific nesting pattern for TNLs. Surprisingly, we found a moderate bootstrap support (BS = 50%) for CNL-A clade being nested within TNL clade making both the CNL and TNL clades paraphyletic. Arabidopsis and sunflower showed 87 syntenic blocks with 1049 high synteny hits between chromosome 5 of Arabidopsis and chromosome 6 of sunflower. Expression data revealed functional divergence of the NBS genes with basal level tissue-specific expression. This study represents the first genome-wide identification of NBS genes in sunflower paving avenues for functional characterization and potential crop improvement.

Data on the genome-wide identification of CNL R-genes in Setaria italica (L.) P. Beauv

We report data associated with the identification of 242 disease resistance genes (R-genes) in the genome of Setaria italica as presented in “Genetic diversity of disease resistance genes in foxtail millet (Setaria italica L.)” (Andersen and Nepal, 2017) [1]. Our data describe the structure and evolution of the Coiled-coil, Nucleotide-binding site, Leucine-rich repeat (CNL) R-genes in foxtail millet. The CNL genes were identified through rigorous extraction and analysis of recently available plant genome sequences using cutting-edge analytical software. Data visualization includes gene structure diagrams, chromosomal syntenic maps, a chromosomal density plot, and a maximum-likelihood phylogenetic tree comparing Sorghum bicolor, Panicum virgatum, Setaria italica, and Arabidopsis thaliana. Compilation of InterProScan annotations, Gene Ontology (GO) annotations, and Basic Local Alignment Search Tool (BLAST) results for the 242 R-genes identified in the foxtail millet genome are also included in tabular format.

Genome-wide survey of switchgrass NACs family provides new insights into motif and structure arrangements and reveals stress-related and tissue-specific NACs

NAC proteins comprise of a plant-specific transcription factor (TF) family and play important roles in plant development and stress responses. Switchgrass (Panicum virgatum) is the prime candidate and model bioenergy grass across the world. Excavating agronomically valuable genes is important for switchgrass molecular breeding. In this study, a total of 251 switchgrass NAC (PvNACs) family genes clustered into 19 subgroups were analyzed, and those potentially involved in stress response or tissue-specific expression patterns were pinpointed. Specifically, 27 PvNACs were considered as abiotic stress-related including four membrane-associated ones. Among 40 tissue-specific PvNACs expression patterns eight factors were identified that might be relevant for lignin biosynthesis and/or secondary cell wall formation. Conserved functional domains and motifs were also identified among the PvNACs and potential association between these motifs and their predicted functions were proposed, that might encourage experimental studies to use PvNACs as possible targets to improve biomass production and abiotic stress tolerance.