ABSTRACT
Brachypodium distachyon has recently emerged as a premier model plant for monocot biology, akin to Arabidopsis thaliana. We previously reported genome-wide transcriptomic and alternative splicing changes occurring in Brachypodium during compatible infections with Panicum mosaic virus (PMV) and its satellite virus (SPMV). Here, we dissected the role of Brachypodium phenylalanine ammonia lyase 1 (PAL1), a key enzyme for phenylpropanoid and salicylic acid (SA) biosynthesis and the induction of plant defenses. Targeted metabolomics profiling of PMV-infected and PMV- plus SPMV-infected (PMV/SPMV) Brachypodium plants revealed enhanced levels of multiple defense-related hormones and metabolites such as cinnamic acid, SA, and fatty acids and lignin precursors during disease progression. The virus-induced accumulation of SA and lignin was significantly suppressed upon knockdown of B. distachyon PAL1 (BdPAL1) using RNA interference (RNAi). The compromised SA accumulation in PMV/SPMV-infected BdPAL1 RNAi plants correlated with weaker induction of multiple SA-related defense gene markers (pathogenesis related 1 [PR-1], PR-3, PR-5, and WRKY75) and enhanced susceptibility to PMV/SPMV compared to that of wild-type (WT) plants. Furthermore, exogenous application of SA alleviated the PMV/SPMV necrotic disease phenotypes and delayed plant death caused by single and mixed infections. Together, our results support an antiviral role for BdPAL1 during compatible host-virus interaction, perhaps as a last resort attempt to rescue the infected plant.
IMPORTANCE Although the role of plant defense mechanisms against viruses are relatively well studied in dicots and in incompatible plant-microbe interactions, studies of their roles in compatible interactions and in grasses are lagging behind. In this study, we leveraged the emerging grass model Brachypodium and genetic resources to dissect Panicum mosaic virus (PMV)- and its satellite virus (SPMV)-compatible grass-virus interactions. We found a significant role for PAL1 in the production of salicylic acid (SA) in response to PMV/SPMV infections and that SA is an essential component of the defense response preventing the plant from succumbing to viral infection. Our results suggest a convergent role for the SA defense pathway in both compatible and incompatible plant-virus interactions and underscore the utility of Brachypodium for grass-virus biology.
INTRODUCTION
Plants, under relentless challenge by pests and pathogens, have evolved a variety of defense mechanisms. The basal host response to pathogen- or microbe-associated molecular patterns (PAMPs or MAMPs, respectively) (1–3) is the activation of specialized transmembrane proteins, pattern recognition receptors (PRRs). This first layer of defense known as P/MAMP-triggered immunity (PTI) is, in turn, suppressed by pathogen-harbored effectors. While a large body of knowledge exists for bacterial pathogen-triggered PTI mechanisms, less is known for virus-triggered PTI responses. A recent study showed that brassinosteroid-insensitive 1-associated kinase 1 (BAK1) or BAK1-like 1 (BKK1) induced immunity in Arabidopsis against Turnip crinkle virus (TCV) infection, and loss-of-function mutations in BAK1 or BKK1 enhanced susceptibility to virus infection (4). Another study demonstrated that Arabidopsis bak1 mutants exhibited increased susceptibility to three RNA viruses (TCV, Tobacco mosaic virus [TMV], and Oilseed rape mosaic virus [ORMV]) during compatible interactions, and crude viral extracts of inoculated leaf tissues induced several PTI markers in a BAK1-dependent manner (5). These studies showed that BAK1-dependent PTI contributed to host antiviral responses.
An effector interferes with various cellular defense pathways to help promote virulence. These effectors are either secreted or, in the case of viruses, expressed within the plant cells. To counter the action of effectors, the host cell responds with a second line of defense, effector-triggered immunity (ETI). Generally, ETI is triggered when a host receptor recognizes an effector, typically a nucleotide binding site-leucine rich repeat (NBS-LRR) protein. NBS-LRR proteins are a type of resistance (R) proteins that prompt a cascade of biochemical reactions that result in a localized hypersensitive response (HR) and/or systemic acquired resistance (SAR) to limit pathogen spread (1, 2). With regard to antiviral defenses, both PTI and ETI can activate downstream defense hormone pathways (e.g., salicylic acid [SA], jasmonic acid [JA], ethylene, or nitric oxide) or mediate changes in the cell wall due to altered composition of lignin, cellulose, and fatty acids; such compounds also act as signaling molecules and precursors to defense hormones and affect host gene expression (6). For instance, earlier studies of an antiviral role of SA in tobacco revealed that SA inhibits replication of Tobacco mosaic virus (TMV) and systemic movement of Cucumber mosaic virus (CMV) in tobacco (7–10). SA application also resulted in reduced accumulation of Potato virus X (PVX) in inoculated tobacco (8). Conversely, unrestricted spread of Potato virus Y (PVY) and enhanced disease symptoms were observed in NahG potato plants that were depleted of SA (11, 12). In bean and cowpea plants, systemic movement of Alfalfa mosaic virus (AMV) was reduced by 90% upon exogenous SA application (13). These and other studies suggest that SA limits virus infection by interfering in multiple steps of viral replication and intercellular spread as well as systemic movement, depending on the virus-host system (14–16).
Recent studies of compatible plant-microbe interactions have certainly advanced our understanding of plant immune responses and strategies used by viruses to target host defenses. Typically, these studies have relied on the dicot model plant Arabidopsis thaliana and Nicotiana spp. However, our understanding of the genetic basis of antiviral defenses in monocots has been at a disadvantage due to the lack of a suitable monocot model plant. Within the past decade, Brachypodium distachyon has been developed as a superb model system for plant research and has opened new avenues for the study of the Poaceae (Gramineae) to include host defense responses to microbes, plant biochemistry, and evolutionary biology (17, 18). The Brachypodium distachyon research community has developed multiple open-access genomic and genetic resources. With the recent expansion of pan-genomics resources for related species in the genus Brachypodium (19), it is now possible to evaluate grass host defense responses to pathogens of significant environmental and economic importance.
For the study of genomic and transcriptomic outcomes of virus infections, Brachypodium has shown its merits. We and others have demonstrated this plant species’ utility in studying Panicum mosaic virus (PMV) and its satellite virus (SPMV), Foxtail mosaic virus (FMV), Barley stripe mosaic virus (BSMV), Brome mosaic virus (BMV), Wheat streak mosaic virus (WSMV), Maize mild mottle virus (MMMV), Sugarcane mosaic virus (SCMV), and Sorghum yellow banding virus (SYBV) (20–28). Brachypodium also supports infections by Sugarcane mosaic virus, Bamboo mosaic virus, and Barley yellow dwarf virus (BYDV-GRV) (29–31). Taken together, Brachypodium is well suited for studies of grass-infecting viruses in several genera.
PMV is the type member of the genus Panicovirus within the family Tombusviridae. The 4,326-nucleotide (nt), single-stranded positive-sense RNA is enveloped in 30-nm icosahedral virions. PMV encodes six open reading frames: two from the genomic RNA to express the replicase-associated proteins (P48 and P112) (32) and four from the subgenomic polycistronic RNA to express a 26-kDa capsid protein (CP) and three nonstructural proteins (P8, P6.6, and P15) that facilitate virus movement (20, 32, 33). PMV is unusual in that it supports replication and movement of two single-stranded RNA satellites: an 826-nt satellite virus (SPMV) and two ∼400-nt satellite RNAs (satRNAs) (27, 28, 34, 35). PMV and its satellites are found in naturally occurring infections of St. Augustinegrass and switchgrass (27, 28, 34–36). SPMV RNA encodes a 17-kDa CP to assemble 16-nm spherical T=1 virions (37, 38). In addition, the SPMV CP elicits a nonhost-hypersensitive response in Nicotiana benthamiana (39). The presence of SPMV in coinfection with PMV exacerbates disease symptoms on systemic leaves in millets and Brachypodium (24, 25). Using PMV, SPMV, and Brachypodium, we have focused our studies on the biology of host-virus interactions.
One early marker of the host response to infection, including to some viruses, is the induction of SA (14). Phenylalanine ammonia lyase (PAL) is a key enzyme for the production of SA and various secondary metabolites—such as monolignols, which are polymerized to form lignin, flavonoids, hydroxycinnamic acid, phenolics, and anthocyanin—which play important roles in plant defense (40–42). Lignin is particularly important, serving as a constituent of secondary cell walls and being involved in defense against several pathogens (43–45). The deposition of lignin strengthens the cell wall and provides a physical barrier to pathogen ingress (44–49).
Our previous work with PMV and PMV plus SPMV (PMV+SPMV) infections of Brachypodium identified genome-wide changes in host gene expression profiles and alternative splicing landscapes (24, 26). Here, we have determined the extent of perturbation in defense-related metabolites during PMV or PMV+SPMV infections (PMV/SPMV) and the function of PAL in antiviral defenses. PAL activity was induced during PMV/SPMV infections, paralleling responses observed for other biotic and abiotic stresses (45, 50–53). PMV/SPMV infections of B. distachyon PAL1 (BdPAL1) RNA interference (RNAi) lines showed lower pathogen-induced SA and lignin production and resulted in exacerbation of disease symptoms compared to that in wild-type (WT) plants. These negative effects were abrogated by exogenous application of SA, thus providing evidence that PAL-mediated SA production in response to PMV/SPMV infection is an important host mechanism for antiviral defenses in grasses.
RESULTS AND DISCUSSION
Identification of primary and secondary metabolic pathways altered by virus infection in Brachypodium.Our previous studies identified widespread transcriptome and splicing changes in response to PMV/SPMV infections in Brachypodium (24, 26). Genes encoding enzymes putatively involved in plant metabolic and biosynthetic pathways affected by PMV infection were identified using the differential gene expression data (stage II, ∼10 to 14 days after infection) and MapMan metabolic pathway analyses (54, 55) (Fig. 1; see also Data Set S1 in the supplemental material). This included gene expression changes predicted to affect production of primary metabolites such as carbohydrates (e.g., glucose and sucrose) that provide cellular energy and function as plant defense signaling molecules (56–58). MapMan analysis of transcripts predicted downregulation of primary metabolites related to light reactions in photosynthesis, which are associated with leaf chlorosis and decreased chlorophyll levels in early infection stages (24, 26). Such breakdown of chloroplast function and the associated leaf chlorosis are known outcomes of plant virus and other plant pathogen infections (59–61).
MapMan overview of primary and secondary metabolic pathways perturbed during PMV infection in Brachypodium. The square boxes in green and red indicate reduced and induced level of expression (>2 log2 fold-change, FDR < 0.05), respectively. The fold change expression data are provided in Data Set 1 in the supplemental material. Mito, mitochondrial; TCA, tricarboxylic acid; OPP, oxidative pentose phosphate; CHO, carbohydrates.
DATA SET S1
This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.
Metabolic pathways related to cell wall and lipid biosynthesis were also downregulated in PMV-infected plants (Fig. 1). This was associated with an overall decrease in the biomass and the stature of virus-infected plants at later infection stages (24, 26). In contrast, central metabolic pathways associated with respiration, such as the mitochondrial electron transport chain, Krebs cycle, and the oxidative pentose phosphate (OPP) pathway, were induced. These changes may reflect the high demand for energy production in virus-infected plants to sustain cellular activities related to adaptation and defense. Secondary metabolites are also important for host responses to biotic and abiotic stress (62).
In PMV-infected plants, multiple secondary metabolite pathways were substantially induced. These included pathways involved in the production of lignin, flavonoids, phenylpropanoids, phenolics, waxes, and terpenes, all of which play crucial roles in host defense and adaptation to plant stress (Fig. 1). Several nucleotide and amino acid biosynthetic pathways (related to phenylalanine, tyrosine, tryptophan, glycine, serine, and cysteine) were also upregulated, possibly as a consequence of a reprogramming of the cellular transcriptome and proteome. Alternatively, increased amounts of amino acids and nucleotides could be redirected to biosynthesis of secondary metabolites or serve as cofactors for proteins involved in plant defense signal transduction. This bioinformatic analysis of primary and secondary metabolic pathways provided a framework to explore the host mechanisms resulting in phenotypic changes and symptomology induced by virus infection.
PMV and SPMV infection perturb defense hormone and fatty acid profiles.Salicylic acid (SA) and its derivative methyl-SA (MeSA) are key hormones mediating plant defense responses to pathogens, including viruses (1, 63–65). SA-mediated defenses are induced in both ETI and PTI and during incompatible and compatible interactions. PMV infection modulates the expression profiles of genes in defense hormone pathways such as SA and jasmonic acid (JA) (24, 26). Several pathogenesis-related (PR) genes such as PR-1, PR-3, and PR-5 were found to be upregulated in PMV/SPMV-infected Brachypodium (25).
Here, we aimed to investigate if gene expression changes reflected the expected changes in defense-related hormones and their precursors. Brachypodium plants were challenged with PMV or with PMV plus SPMV (PMV+SPMV), followed by targeted metabolic profiling using gas chromatography mass spectrometry (GC-MS). To determine the temporal dynamics of infection, multiple time points were inspected on inoculated leaves: early stage, no visual symptoms (7 days postinfection [dpi], stage I); mid stage, with mild chlorosis symptoms (14 dpi, stage II); and late stage, with severe chlorosis and necrosis symptoms (21 dpi, stage III), as described previously (24). Significant changes in SA accumulation were found in all three stages upon virus infection. The accumulation of cinnamic acid (CA) in PMV- and PMV+SPMV-infected plants was elevated at 14 dpi but not at 7 dpi. The accumulation of SA was significantly higher (P ≤ 0.05) in PMV-infected (1,747 ng/g [fresh weight]) and PMV+SPMV-infected (2,570 ng/g [fresh weight]) plants by stage II than in mock-inoculated plants (SA, 815 ng/g [fresh weight]), and these levels remained elevated through stage III (Fig. 2A). CA concentration was elevated but not significantly affected at stage II upon PMV/SPMV infections; however, it was significantly induced (P ≤ 0.05) at stage III (Fig. 2B). SA and CA concentrations were consistently higher in PMV+SPMV-inoculated plants at stage III than in PMV-inoculated plants. The higher levels of SA and CA in virus-infected plants correspond with the upregulation of several SA signaling genes, including pathogenesis-related (PR) proteins (24, 25). We previously reported that the PMV+SPMV synergism exacerbated disease symptoms and had an additive effect on several defense-related gene expression profiles compared to that from PMV-only infections (24). Such additive affects were also observed in the metabolite profiles, with a higher accumulation of SA and CA in PMV+SPMV-coinfected plants than in plants infected with PMV alone (Fig. 2A and B).
Quantification of salicylic acid (SA), cinnamic acid (CA), jasmonic acid (JA), and precursors in PMV- and PMV+SPMV-infected Brachypodium. Levels of SA (A), CA (B), trans- and cis-JA (C), and trans- and cis-12-OPDA (D) were determined using GC-MS. The fresh weight (FW) leaf samples were collected at three different stages of disease progression (early stage I, 7 dpi; mid-stage II, 14 dpi; late stage III, 21 dpi), as described previously (24). The y axes represent average contents of respective metabolites from five biological replicates. The error bars represent standard errors of means (n = 3 and n = 4 for early and for mid and late stages, respectively). *, P ≤ 0.05; **, P ≤ 0.01 between mock- and virus-infected samples as determined by two-sample t test (one-tailed). M, mock; P, PMV; PS, PMV+SPMV; OPDA, 12-oxo-phytodienoate.
In a manner similar to that for SA, jasmonic acid (JA) and its derivative methyl-JA (meJA) regulate various physiological and developmental processes as well as responses to biotic stresses, particularly to necrotrophic pathogens, wounding, and herbivores (66–69). Biosynthesis of JA is initiated in the chloroplast with the conversion of linolenic acid to 12-oxo-phytodienoate (12-OPDA), followed by its reduction in the peroxisome by a 12-OPDA reductase through β-oxidation (70, 71). With some exceptions, the JA and SA pathways are antagonistic to each other (25, 72–74). The expression of several JA signaling genes was downregulated in PMV/SPMV-infected Brachypodium, while SA signaling genes were upregulated (24, 25). Metabolite profiling analysis revealed that at stage III, when SA was highest in virus-infected plants (Fig. 2A), cis-JA was significantly lower (P ≤ 0.05) than in mock-infected plants (Fig. 2C), suggesting SA-JA antagonism. Because the cis isoform of JA was highly affected compared to trans, this indicates a downregulation of biosynthesis rather than increased JA conversion or turnover. Interestingly, 12-OPDA, the precursor to JA synthesis, was significantly higher (P ≤ 0.01) in PMV+SPMV-inoculated plants at stage I (Fig. 2D), when both SA and JA were largely unaltered from mock-inoculated plants. This suggests that 12-OPDA or related oxylipins act as signaling molecules, independent of JA, during the early stages of virus infection (75, 76).
Fatty acids (FAs) are important structural and metabolic constituents of membranes and the cell wall and are involved in intracellular signaling and protein modification processes (77–79). For example, 18-carbon mono- and polyunsaturated FAs such as oleic acid (18:1), linoleic acid (18:2), and linolenic acid (18:3) are precursors to oxylipin-derived defense hormones, with key roles in affecting outcomes in plant-microbe interactions (80–86). We next investigated total and free 18-carbon FA levels in Brachypodium leaves during virus infection at stages II and III (Fig. 3A to F). In stage II, there were no significant changes in levels of total FAs except for linoleic acid (18:2) (Fig. 3A to C). However, at stage III, levels of total oleic acid (18:1) and linoleic acid (18:2) were significantly greater (P ≤ 0.01) in virus-infected plants than in mock-infected plants (Fig. 3A to C). In contrast, total linolenic acid (18:3) levels were significantly lower (P ≤ 0.01) in virus-infected leaves (Fig. 3C).
Quantification of 18-carbon unsaturated fatty acid levels in PMV- and PMV+SPMV-infected Brachypodium. Levels of free and total 18:1, 18:2, and 18:3 fatty acids in mock-, PMV-, and PMV+SPMV-infected plants at stages II (14 dpi) and III (21 dpi) of infection. The error bars represent standard errors of means (n = 3 and n = 4 for stages II and III, respectively). *, P ≤ 0.05; **, P ≤ 0.01; ***, P ≤ 0.001 between mock- and virus-infected samples as determined by two-sample t test (one-tailed). M, mock; P, PMV; PS, PMV+SPMV.
The bottleneck in total 18:3 accumulation corresponds to the decreased expression of FAD7 in PMV/SPMV-infected Brachypodium that we previously observed (24). FAD7 encodes a chloroplast-localized enzyme that desaturates 18:2 to 18:3 (87), and the suppressed expression of FAD7 could result in the accumulation of its substrate, 18:2. Plant membranes, especially chloroplast membranes, are also a repository for 18:3, esterified in glycolipids and phospholipids. Upon cell perturbation or damage, phospholipases are activated, which release free 18:3 from membranes (88). Interestingly, levels of free 18:3 were significantly higher (P ≤ 0.01) in virus-infected plants, as early as stage II (Fig. 3F), suggesting that despite the bottleneck in biosynthesis, during PMV/SPMV infection, there is a greater amount of free 18:3 being released from repositories and/or as breakdown products. Given that plastidic FA-derived lipid signals are critical for modulating plant defense signaling (84, 89, 90), it is possible that the free 18:3 or derived lipid signals are important in host responses to PMV/SPMV infections.
BdPAL1 is required for SA and lignin accumulation during PMV and PMV+SPMV infections.SA is derived from chorismate or phenylalanine by two key enzymatic pathways involving isochorismate synthase (ICS) or phenylalanine ammonia lyase (PAL), respectively (91–93). Although the ICS pathway is thought to be dominant in Arabidopsis (94), other plant species have various contributions of either ICS or PAL to SA production during infection. For instance, SA is predominantly derived from the phenylpropanoid pathway in tobacco (95). In contrast, in soybean, both ICS and PAL pathways contribute equally to SA biosynthesis upon pathogen infection (96). It is unclear which pathway of SA biosynthesis is dominant in Brachypodium. PAL catalyzes the conversion of phenylalanine to trans-cinnamate by nonoxidative deamination. This is a key rate-limiting step between primary and secondary metabolism (40, 41, 97, 98). The PAL pathway also leads to production of several other critical defense-related metabolites such as lignin and flavonoids. Brachypodium encodes eight putative BdPAL homologs. BdPAL1 (Bradi3g49250) is predominantly expressed in vegetative tissues, and its downregulation by RNA interference (RNAi) resulted in significant changes in defense-related secondary metabolites and enhanced susceptibility to fungal pathogens (50).
We previously showed that PMV/SPMV infection activated multiple SA-related genes (e.g., PR genes) (24, 99). Expression of a few PAL isoforms was also upregulated (data not shown). Because PMV/SPMV infections affected SA biosynthesis (Fig. 2) along with multiple primary and secondary metabolic pathways, including cell wall-related components (Fig. 1), we further examined the role of PAL in PMV/SPMV-induced SA and phenylpropanoid biosynthesis by leveraging the BdPAL1 RNAi plants (50).
Prior to investigating the host effects of BdPAL1 knockdown on PMV/SPMV infections, we reconfirmed the silencing of BdPAL1 expression and activity in the BdPAL1 RNAi plants compared to those in WT plants (Fig. 4A and B). As anticipated, PMV and PMV+SPMV infections significantly induced PAL activity in WT (P ≤ 0.05) compared to that in mock-inoculated leaves (Fig. 4B). Although PAL activity was induced in PMV- and PMV+SPMV-infected BdPAL1 RNAi plants (P ≤ 0.05), it was consistently lower than in the corresponding WT plants, because of the PAL1 downregulation (Fig. 4B). Next, we quantified defense-related and primary and secondary cell wall-associated compounds in WT and BdPAL1 RNAi plants, in mock-inoculated and PMV/SPMV-infected plants. Among the studied compounds, sucrose levels did not follow any specific pattern of change during virus infection; however, they were lower in healthy BdPAL1 RNAi plants than in WT plants (Fig. 4C). We speculate the decrease in sucrose could be a reflection of PAL function in normal plant growth and development. BdPAL1 RNAi plants do have slightly delayed growth compared to that of the WT (50). Glucose, the end product of photosynthesis, was significantly lower in both WT and BdPAL1 RNAi plants infected with PMV/SPMV than in their healthy counterparts (P ≤ 0.05) (Fig. 4D). The levels of glucose correlated with the decreased level of various enzymes involved in photosynthetic and primary carbohydrate metabolic pathways (Fig. 1) as a result of PMV+SPMV infection, and the observed chlorotic symptomology and stunted growth of the PMV+SPMV-infected plants (Fig. 5C). Although it is difficult to determine if the lowered primary carbohydrates are the cause or consequence of the growth delays and disease phenotypes, they underscore the broader function of PAL1 and their metabolites, including SA, on plant growth, development, and defense processes (50, 100).
Functional analysis of wild-type (WT) and BdPAL1 RNAi plants during PMV and PMV+SPMV infection. Expression of BdPAL1 (A) and PAL activity (B) in mock-, PMV-, and PMV+SPMV-inoculated WT and BdPAL1 RNAi plants at stage II (14 dpi) of infection. Levels of soluble carbohydrates (sucrose and glucose) (C), structural carbohydrates (cellulose and hemicellulose) and total lignin (D), and salicylic acid (E) at stage II (14 dpi) of infection. M, mock; P, PMV; PS, PMV+SPMV. Relative expression of defense-related genes (SA and JA signaling components) in PMV-infected (F) and PMV+SPMV-infected (G) WT and BdPAL1 RNAi plants. The error bars represent standard errors of means (n = 3 for panels A to D, F, and G and n = 5 for panel E). Statistically significant differences in panel A were determined by two-sample t test (one-tailed) and, in panels B to G, were assessed using one or two-way ANOVA followed by the Tukey’s test. Unlike lowercase letters represent significant differences among the group means (P ≤ 0.05). F and t test statistics of ANOVA and two-sample t tests are indicated.
BdPAL1 RNAi plants showed enhanced susceptibility to PMV+SPMV infection. (A) Mock- and PMV+SPMV-infected WT (left) and BdPAL1 RNAi (right) plants at stage III (21 dpi). (B) Closeup of mock- and PMV+SPMV-infected leaves of WT (left) and BdPAL1 RNAi (right) plants at stage III (21 dpi). (C) Percentages of leaves with necrosis of wild-type (WT) and BdPAL1 RNAi plants, mock- and virus-infected plants at stage III (21 dpi). Statistically significant differences were assessed using one-way ANOVA followed by the Tukey’s multiple-comparison test. Unlike lowercase letters represent significant differences among the group means (P ≤ 0.05). F statistics of ANOVA are indicated. (D) RT-qPCR detection of mRNA encoding PMV CP (PCP) and SPMV CP (SPCP) in noninoculated leaves at 14 dpi. M, mock; PS, PMV+SPMV. Statistically significant differences were assessed between mock- and virus-inoculated samples using two-sample t test (one-tailed). *, P ≤ 0.05; ***, P ≤ 0.001. The error bars represent standard errors of the means (n = 3).
Among structural cell wall components, total lignin was significantly induced in WT and BdPAL1 RNAi upon PMV/SPMV infections (P ≤ 0.001), although compared with that in WT, BdPAL1 RNAi plants exhibited low total lignin levels (Fig. 4E). It is possible that the elevated levels of total lignin in WT plants, derived from PAL, may contribute to limiting cell-to-cell movement of PMV and SPMV and/or be a general consequence of tissue stress. Levels of cellulose remained unaltered in healthy and virus-infected WT and BdPAL1 RNAi plants (Fig. 4F). BdPAL1 RNAi lines showed a statistically significant decrease in hemicellulose (P ≤ 0.05), compared to that in WT plants (Fig. 4F). This could be due to the compromised ferulate (FA) production in the BdPAL1 RNAi plants (50). FA is a critical moiety present in grass cell walls, which helps cross-link the cell wall components such as hemicellulose (50). Alternatively, lowered sucrose (Fig. 4C), a carbon source for hemicellulose, may negatively affect its production/levels. There were no significant differences in hemicellulose upon virus infection (Fig. 4F).
Metabolic profiling of SA revealed that virus infection triggered significant increases in SA levels in WT plants (P ≤ 0.05) (Fig. 4G). SA accumulation was slightly higher in virus-infected BdPAL1 RNAi plants than in mock-treated plants but was not statistically significant upon virus infection (Fig. 4G). However, in a manner similar to that for total lignin, virus-triggered SA accumulation was significantly suppressed (P ≤ 0.05) in the BdPAL1 RNAi plants compared to that in WT plants (Fig. 4G). Basal levels of SA in WT and BdPAL1 RNAi plants were largely unaffected (Fig. 4G), suggesting that BdPAL1 primarily contributes to induce SA production upon pathogen infection. The compromised SA levels in PMV+SPMV-infected BdPAL1 RNAi plants correlated with weaker induction of multiple SA-related defense gene markers (PR-1-like, PR-3-like, PR-5-like, and WRKY75-like) (Fig. 4H and I). For instance, PMV/SPMV caused an 18- to 22-fold increase in the PR-1-like gene in WT plants, compared to an ∼2- to 12-fold increase in BdPAL1 RNAi plants. Consistent with previous studies (24, 25, 101–104), we also observed the antagonistic relationship between SA and JA, wherein several JA components (AOS, LOX2, and FAD7) were downregulated in PMV/SPMV-infected WT and BdPAL1 RNAi plants (Fig. 4H to I), which is in agreement with our transcriptomic analyses (24, 25). The suppressed expression of LOX2, AOS, and FAD7 was associated with decreased accumulation of JA in WT plants (Fig. 2C). These findings support the notion that an antagonistic relationship exists between SA and JA during PMV/SPMV infections in Brachypodium.
Inhibition of BdPAL1 results in enhanced susceptibility to PMV and SPMV in Brachypodium.Next, we determined whether the compromised lignin and SA levels and SA-related defense gene activation in the BdPAL1 RNAi plants influence PMV/SPMV disease outcomes. For this, we assessed the disease symptoms (e.g., chlorotic and necrotic lesions) of WT and BdPAL1 RNAi plants. Typical chlorotic symptoms appeared by 12 to 14 dpi of virus infection in both WT and BdPAL1 RNAi plants, and necrotic symptoms appeared by 21 dpi (Fig. 5A and B). At this time point, the percentage necrosis on the upper noninoculated (systemically infected) leaves was quantified. Severe necrotic lesions were observed on greater than 72% of the leaves of PMV+SPMV-infected BdPAL1 RNAi plants, whereas only ∼5% of the WT leaves displayed necrotic lesions (Fig. 5C). The presence of PMV and SPMV RNA in the respective samples was confirmed by reverse transcription-quantitative PCR (RT-qPCR) (Fig. 5D). Levels of SPMV were significantly higher (P ≤ 0.001) in BdPAL1 RNAi plants (Fig. 5D) than in WT plants and paralleled the severe necrotic symptoms observed on the PMV/SPMV-infected BdPAL1 RNAi plants (Fig. 5B and C). We also observed that 14% of the leaves of nonchallenged BdPAL1 RNAi plants also showed necrotic lesions. These symptoms were reminiscent of spontaneous lesions on lesion mimic mutants (LMMs) which also occur in nonchallenged plants. Several genes are linked to LMMs, many of them involved in plant development, metabolism, and defense signaling (105).
Exogenous SA treatment alleviates the hypersensitive phenotype of BdPAL1 RNAi plants.Lastly, we investigated the effects of exogenous supplementation of SA—could it rescue the enhanced PMV/SPMV susceptibility phenotype of BdPAL1 RNAi plants? For this, soil containing WT and BdPAL1 RNAi plants was drenched with 50 ml of water or 100 ppm SA, 24 h prior to inoculation with PMV alone or PMV+SPMV (106). SA was reapplied 24 h postinoculation and then once a week for up to 4 weeks. As expected, on PMV/SPMV-inoculated plants, symptom onset of WT plants, typified by leaf chlorosis, occurred ∼12 to 14 dpi. However, these symptoms were delayed by ∼4 days for SA-treated WT and BdPAL1 RNAi plants. The application of exogenous SA alleviated the necrotic disease phenotype and subsequent plant death at 21 dpi (Fig. 6A) and remained protective through 42 dpi (Fig. 6B). The lower percentage of chlorosis and necrosis on BdPAL1 RNAi plants drenched with SA solution was significant, to levels almost comparable to that for WT plants (P ≤ 0.05) (Fig. 6C and D). It appears the exogenous SA treatment rescued the enhanced susceptibility of BdPAL1 RNAi plants to PMV/SPMV infection. Regarding the possible mechanism of the SA effect on PMV/SPMV, we speculate that the SA applied to the roots, or a mobile SA form (107), could be triggering a systemic acquired resistance (SAR) response. The SAR signal could be perpetuated from roots to distant shoot tissues. Early studies showed that SA not only inhibits virus replication in inoculated tissues but also limits virus cell-to-cell and systemic movement via systemic resistance pathways (7–10, 108, 109). Recent work further bolstered the role of SA-mediated signaling in SAR responses against several biotrophic pathogens, including viruses (11, 12, 107, 110–112). Given that PMV and SPMV accumulate both in shoots and roots (24), it is conceivable that an SA-mediated SAR response originating in the roots could limit PMV/SPMV accumulation on a whole-plant level and attenuate the disease symptoms. Regardless of the exact mechanisms, the results suggest that the induction of the SA pathway in WT plants or supplementation of SA to BdPAL1 knockdown plants can confer tolerance to PMV/SPMV infections.
Exogenous application of salicylic acid (SA) attenuates PMV and SPMV disease symptoms. Representative chlorosis and necrosis at 21 dpi (A) and overall stunting symptoms at 42 dpi (B) of PMV- and PMV+SPMV-infected plants treated with water or salicylic acid (SA; 100 ppm). Quantification of percent chlorotic (C) and necrotic (D) leaves in wild-type (WT) and BdPAL1 RNAi plants treated either with water or SA (100 ppm) at 21 dpi. Statistically significant differences were assessed using one-way ANOVA followed by the Tukey’s test. Unlike lowercase letters represent significant differences among the group means (P ≤ 0.05). F test statistics of ANOVA are indicated. The error bars represent standard errors of means (n = 3). M, mock; P, PMV; PS, PMV+SPMV.
Conclusion.We previously established the emerging model grass Brachypodium as an alternate host to study PMV- and SPMV-host interactions and viral synergism (24) and demonstrated by genomic and transcriptomic approaches that several defense-related hormones and signaling pathways are perturbed in PMV/SPMV interactions with Brachypodium (25). Here, we used Brachypodium and PMV/SPMV as a model pathosystem to gain a broader understanding of the defense-related metabolic pathways perturbed during compatible grass virus infections and the role of BdPAL1 in antiviral defenses. We have shown that PMV/SPMV infections triggered changes in various defense-related metabolites, including the defense hormone salicylic acid (SA) and immune-related genes. Specifically, BdPAL1-induced SA is important for tolerating PMV/SPMV. This metabolic, biochemical, and genetic characterization extended our knowledge of PAL-mediated antiviral defenses in a C3 grass and underscores the utility of Brachypodium as a superb model system for the study of grass-microbe interactions. Further translational studies aimed to enhance the PAL pathway or similar conserved immune features of grasses will likely lead to new measures to develop robust disease resistance in the field.
MATERIALS AND METHODS
Plant material, growth conditions, and virus inoculation.Seeds of Brachypodium wild type (WT, accession Bd21-3) and the BdPAL1 RNAi line were planted in Metro Mix-366 (Sun Gro Horticulture) in 2- by 3-in. pots, stratified at 4°C for 4 days. Information on generation and characterization of BdPAL1 RNAi lines was described previously (50). Plants were maintained in a growth chamber with a photoperiod of 14 h of light (21°C)/10 h of dark (18°C) with light intensity of approximately 250 to 280 mmol/m2s (24). Brachypodium leaves were rub inoculated at the three-leaf stage using in vitro transcribed PMV and SPMV genomic RNA (gRNA), as described previously (24, 113). Briefly, full-length viral cDNA clones of PMV and SPMV in pUC19 were linearized using EcoICR1 and BglII endonucleases, respectively. Linearized infectious viral clones (∼500 ng) were used for in vitro synthesis of gRNA using HiScribe T7 in vitro transcription kit (New England Biolabs). After confirming the size and integrity of RNA by gel electrophoresis, PMV or PMV+SPMV gRNAs were mixed (50% [vol/vol]) in RNA inoculation buffer (50 mM KH2PO4, 50 mM glycine, pH 9.0, 1% bentonite, and 1% Celite) and rub inoculated on the lower-most leaves and cotyledons (∼300 ng of transcript per plant). Throughout this paper, PMV/SPMV indicates results associated with both PMV alone and PMV+SPMV infections. Mock inoculations consisted of RNA inoculation buffer alone. All mock-infected and virus-inoculated plants were covered by a Humi-Dome (Hummert International), incubated in the dark for 24 h, and then moved to the growth chamber.
RNA isolation, RT-qPCR, and molecular diagnostics.The presence of PMV and SPMV RNA in infected plants was confirmed by quantitative RT-PCR detection of noninoculated (systemically infected) leaves as described previously (24). Briefly, total RNA was isolated from 100 mg of leaf tissue using an RNeasy plant kit (Qiagen) with on-column DNase treatment (Invitrogen). One microgram of RNA was used for cDNA synthesis with a SuperScript IV first-strand cDNA synthesis kit (Invitrogen) and random hexamer primers, followed by PCR amplification to detect PMV and SPMV.
For SA and JA marker gene expression analysis, RNA was isolated from mock- and virus-inoculated leaves at 14 dpi (RNeasy plant kit; Qiagen) according to the manufacturer’s instructions. The RNA concentration and quality were validated using a NanoDrop 1000 and agarose gel electrophoresis, respectively. Two micrograms of RNA was used for cDNA synthesis with a SuperScript IV first-strand cDNA synthesis kit (Invitrogen). Gene expression analysis was performed by RT-qPCR using a CFX384 real-time PCR detection system (Bio-Rad) with iTaq universal SYBR green supermix (Bio-Rad), 1 μl of 2-fold diluted cDNA, and 0.2 μM target specific primers, according to the manufacturer’s instructions. Each assay was performed using three biological replicates (composed of RNA isolated from three different plants) and two technical replicates. The primers were designed using QuantPrime (114) and are listed in Table S1 in the supplemental material. The expression data were normalized to the constitutively expressed ubiquitin 18 (Bradi4g00660) gene (115). The relative quantification of gene expression was performed using the comparative threshold cycle (2−ΔΔCT) method (116) and is presented as relative to mock-inoculated plants, which are set to 1.
TABLE S1
This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.
MapMan pathway analysis and GC-MS metabolite profiling.Brachypodium metabolic pathways affected during PMV infection were determined by MapMan metabolic pathway analysis (54, 55, 117) using the differential gene expression data from a previous transcriptome study (stage II, ∼10 to 14 days after infection [24]). MapMan allows correlating gene expression data placement onto schematic maps or pathways of biological processes. The MapMan mapping units for Brachypodium represent orthologous, nonredundant functional gene ontologies retrieved from similarity searches using BLAST (123) (against Arabidopsis TAIR10, Plant SWISS-Prot, and UniRef90), RPS-BLAST (against NCBI Conserved Domain Database [CDD], EuKaryotic Orthologous Groups [KOG]), and Inter-Pro scan. Only mapping classifications with the highest reliability are retained in these mapping units. Briefly, the normalized expression counts (fragments per kilobase per million [FPKM]) of significantly differentially expressed genes (>2 log2 fold up- or downregulated at a false discovery rate [FDR] of <0.05) were overlaid onto the Brachypodium primary and secondary metabolic maps. A transcript-metabolite correlation map was generated, which provided a bird’s eye view of up- and downregulated pathways, indicated by red and green colored squares, respectively (Fig. 1).
For quantification of metabolites, PMV-, PMV+SPMV-, and mock-inoculated Brachypodium leaves were collected at 7, 14, and 21 dpi. Frozen ground tissue samples (100 mg) were solvent extracted and partitioned into duplicate samples. For phytohormones and free fatty acids, one set of the partitioned samples was methylated, collected on a polymeric adsorbent using vapor-phase extraction (VPE), and analyzed using GC/isobutane chemical ionization MS (CI-MS) as previously described (118). For total fatty acid analysis, the second set of partitioned samples was dried down under liquid nitrogen, and then 1 ml of 1:10 3 M KOH to methanol (MeOH) was added and samples were incubated at 80°C for 1 h. Following incubation, samples were acidified with 50 μl of 6 M HCl, 1 ml of hexane was added, and samples were vortexed. Then, 500 μl was transferred to a 4-ml glass vial and the samples were methylated, collected on a polymeric adsorbent using vapor-phase extraction (VPE), and analyzed using GC/isobutane CI-MS as previously described (118). Metabolite quantification was based d5-JA (C/D/N Isotopes, Canada), or U-13C18:3 (Cambridge Isotope Laboratories) as internal standards.
Primary and secondary cell wall composition analysis.Cell wall composition analysis, including that of structural carbohydrates (cellulose and hemicellulose), lignin, and nonstructural carbohydrates (glucose and sucrose), of mock- and virus-inoculated Brachypodium leaves was performed at 14 dpi (stage II), as described previously (119, 120). Briefly, leaf samples were collected, freeze-dried, and pulverized in a mortar with a pestle. The dried sample (7 mg) was extracted twice in 2 ml of water with 2 ml of 95% ethanol, sonicated, and centrifuged each time for 20 min. The supernatants were discarded, and the pellets dried overnight at 65°C and subsequently weighed to calculate the percentage of soluble extractives. Next, the dried pellet was hydrolyzed in two steps: (i) by the addition of 70 μl of 72% H2SO4 and incubating for 60 min at 30°C and then (ii) by diluting the acid solution to 3% with ultrapure water and autoclaving for 60 min at 121°C. After autoclaving, the hydrolysis liquor was centrifuged, and 150 μl of the supernatant, diluted to 50%, was used to measure absorbance on a UV-visible spectrophotometer at 240 nm to calculate the percentage of acid soluble lignin. Finally, 1.5 ml of the remaining hydrolyzed liquor was used to analyze structural carbohydrates (percent cellulose and percent hemicellulose) in the sample, as described previously (119, 120). The structural carbohydrate analysis was performed by high-performance liquid chromatography (HPLC) on a system equipped with a refractive index detector using a Shodex sugar SP0810 column after neutralization of the hydrolysis liquor with calcium carbonate. To analyze the acid insoluble lignin, the insoluble residue remaining after hydrolysis was washed with ultrapure water three times to remove the acid, dried overnight at 65°C, and weighed. The percent acid insoluble lignin was estimated by subtracting the percent extractive free ash from the percent acid-insoluble residue. The percent total lignin was the sum of the percent acid-insoluble lignin and percent acid-soluble lignin (119, 120).
Nonstructural carbohydrates were determined as described previously (121, 122), with minor modification. Briefly, 1 ml of 80% ethanol was added to 20 mg of dried sample and incubated in a water bath at 80°C for 20 min. The ethanol in the tubes was decanted, and the process was repeated two times. The decanted ethanol was evaporated by incubating the tube in a water bath at 80°C for 4 h; the remaining volume in the tube was adjusted to 1 ml with ultrapure water, and 40 μl was used to read absorbance at 490 nm. For glucose determination, 180 μl peroxidase/glucose oxidase (PGO) was added, and the sample was incubated at room temperature for 20 min followed by measurement of absorbance at 490 nm. Glucose content was determined by the difference between the background absorbance (before addition of PGO) and the absorbance with PGO (after incubation) and comparing it with a glucose standard curve. For sucrose determination, 40 μl of invertase solution was added to the previous sample and incubated at 35°C for 60 min followed by an absorbance reading at 490 nm. To calculate the sucrose content, the background and glucose absorbances were subtracted from the final absorbance reading using sucrose standard curve.
Measurement of PAL activity.PAL activity was estimated using a phenylalanine ammonia lyase (PAL) microplate assay kit (Cohesion Biosciences) according to the manufacturer’s instructions. Briefly, 100 mg of leaves from mock- and virus-inoculated WT or BdPAL1 RNAi plants was homogenized with 1 ml of assay buffer in a prechilled mortar and pestle. Homogenized samples were centrifuged at 8,000 × g for 15 min at 4°C. The supernatant was transferred to a new 1.5-ml microcentrifuge tube and kept on ice. Assay samples were prepared in a 96-well plate by adding 10 μl of sample, 130 μl of reaction buffer, and 50 μl of l-phenylalanine (substrate) followed by a 30-min incubation at 30°C. After 30 min, 10 μl of stop solution was added, and PAL activity was quantified by measuring absorption spectra (optical density [OD]) at 290 nm. PAL activity was calculated according to manufacturer’s formula shown below:
Here, one unit (U) is defined as the OD value change of 0.01 per minute, and W is the weight of sample in grams.
Exogenous SA treatment and disease monitoring.Brachypodium plants were grown to the three-leaf stage under the growth conditions mentioned previously. Then the WT and BdPAL1 RNAi plants were treated with 50 ml of 100 ppm SA (Sigma-Aldrich) or sterile water by soil drenching, followed by PMV or PMV+SPMV inoculation after 24 h (106). The SA treatment was repeated 24 h postinoculation. Then, the SA application was made once every 7 days for up to 4 weeks. Throughout the experiment, we monitored the progression of symptoms on noninoculated leaves (chlorosis and necrosis) and the overall plant habit. At 21 dpi, the noninoculated leaves were scored for percentage of leaves with chlorosis and necrosis symptoms.
Data analysis and statistics.Data analysis, statistics, and graphs were prepared using Microsoft Excel (version 2009), Python package bioinfokit v1.0.2 (https://github.com/mandadi-lab/bioinfokit), and R package agricolae v1.3-3. Statistically significant differences were analyzed using two-sample t test (one-tailed) or analysis of variance (ANOVA) followed by Tukey’s test for multiple pairwise comparisons.
ACKNOWLEDGMENTS
This project was supported in part by the USDA Agriculture and Food Research Initiative National Institute of Food and Agriculture (NIFA; 2016-67013-24738, HATCH 1023984, and CSP503406 to K.K.M. and K.-B.G.S.) and the US Department of Energy Great Lakes Bioenergy Research Center (Department of Energy, Biological and Environmental Research, Office of Science, grant no. DE-FC02-07ER64494 to J.C.S.).
We thank Brianna Jacques, Gerleene Acuna, and Rolando Mireles (Texas A&M AgriLife Research, Weslaco, TX) for their technical assistance with plant propagation and cell wall composition analysis.
FOOTNOTES
- Received 14 December 2020
- Accepted 9 January 2021
- Published 16 February 2021
This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.
REFERENCES
- 1.↵
- 2.↵
- 3.↵
- 4.↵
- 5.↵
- 6.↵
- 7.↵
- 8.↵
- 9.↵
- 10.↵
- 11.↵
- 12.↵
- 13.↵
- 14.↵
- 15.↵
- 16.↵
- 17.↵
- 18.↵
- 19.↵
- 20.↵
- 21.↵
- 22.↵
- 23.↵
- 24.↵
- 25.↵
- 26.↵
- 27.↵
- 28.↵
- 29.↵
- 30.↵
- 31.↵
- 32.↵
- 33.↵
- 34.↵
- 35.↵
- 36.↵
- 37.↵
- 38.↵
- 39.↵
- 40.↵
- 41.↵
- 42.↵
- 43.↵
- 44.↵
- 45.↵
- 46.↵
- 47.↵
- 48.↵
- 49.↵
- 50.↵
- 51.↵
- 52.↵
- 53.↵
- 54.↵
- 55.↵
- 56.↵
- 57.↵
- 58.↵
- 59.↵
- 60.↵
- 61.↵
- 62.↵
- 63.↵
- 64.↵
- 65.↵
- 66.↵
- 67.↵
- 68.↵
- 69.↵
- 70.↵
- 71.↵
- 72.↵
- 73.↵
- 74.↵
- 75.↵
- 76.↵
- 77.↵
- 78.↵
- 79.↵
- 80.↵
- 81.↵
- 82.↵
- 83.↵
- 84.↵
- 85.↵
- 86.↵
- 87.↵
- 88.↵
- 89.↵
- 90.↵
- 91.↵
- 92.↵
- 93.↵
- 94.↵
- 95.↵
- 96.↵
- 97.↵
- 98.↵
- 99.↵
- 100.↵
- 101.↵
- 102.↵
- 103.↵
- 104.↵
- 105.↵
- 106.↵
- 107.↵
- 108.↵
- 109.↵
- 110.↵
- 111.↵
- 112.↵
- 113.↵
- 114.↵
- 115.↵
- 116.↵
- 117.↵
- 118.↵
- 119.↵
- 120.↵
- 121.↵
- 122.↵
- 123.↵
