Megan McCaghey1†
Dandan Shao1†
Jake Kurcezewski1Ally Lindstrom1
Ashish Ranjan2
Steven A. Whitham3Shawn P. Conley4
Brett Williams5
Damon L. Smith1*
Mehdi Kabbage1*- 1Department of Plant Pathology, University of Wisconsin-Madison, Madison, WI, United States
- 2Department of Plant Pathology, University of Minnesota, St. Paul, MN, United States
- 3Department of Plant Pathology and Microbiology, Iowa State University, Ames, IA, United States
- 4Department of Agronomy, University of Wisconsin-Madison, Madison, WI, United States
- 5Centre for Tropical Crops and Biocommodities, Queensland University of Technology, Brisbane, QLD, Australia
A lack of complete resistance in the current germplasm complicates the management of Sclerotinia stem rot (SSR) caused by Sclerotinia sclerotiorum in soybean. In this study, we used bean pod mottle virus (BPMV) as a vehicle to down-regulate expression of a key enzyme in the production of an important virulence factor in S. sclerotiorum, oxalic acid (OA). Specifically, we targeted a gene encoding oxaloacetate acetylhydrolase (Ssoah1), because Ssoah1 deletion mutants are OA deficient and non-pathogenic on soybean. We first established that S. sclerotiorum can uptake environmental RNAs by monitoring the translocation of Cy3-labeled double-stranded and small interfering RNA (ds/siRNAs) into fungal hyphae using fluorescent confocal microscopy. This translocation led to a significant decrease in Ssoah1 transcript levels in vitro. Inoculation of soybean plants with BPMV vectors targeting Ssoah1 (pBPMV-OA) also led to decreased expression of Ssoah1. Importantly, pBPMV-OA inoculated plants showed enhanced resistance to S. sclerotiorum compared to empty-vector control plants. Our combined results provide evidence supporting the use of HIGS and exogenous applications of ds/siRNAs targeting virulence factors such as OA as viable strategies for the control of SSR in soybean and as discovery tools that can be used to identify previously unknown virulence factors.
Introduction
Sclerotinia sclerotiorum is a broad host range fungal pathogen that infects many dicotyledonous plants worldwide. It is the causal agent of Sclerotinia stem rot (SSR; or white mold) on soybean, a challenging and significant yield-limiting disease (Allen et al., 2017). SSR development is heavily influenced by the weather, and disease onset is favored by cool and wet conditions during flowering (Workneh and Yang, 2000; Peltier et al., 2012; Willbur et al., 2019). SSR causes substantial yield losses to soybean globally (Wrather et al., 2010), and between 2010 and 2014 SSR resulted in yield losses between 150 million and 1.1 billion kgs in the United States, annually (Allen et al., 2017). Current SSR management strategies have limited efficacy and comprise of integrated cultural, chemical, and biological control practices (Peltier et al., 2012; Willbur et al., 2019). Disease control is complicated by the lack of complete genetic resistance to SSR, however, breeding efforts have identified partially resistant soybean genotypes in the laboratory and through field trials (Boland and Hall, 1986; Chun et al., 1987; Kim and Diers, 2000; Iquira et al., 2015; McCaghey et al., 2017).
RNA interference (RNAi) is a gene silencing process that utilizes complementary RNA particles to target mRNA and inhibit gene expression and translation. RNAi has emerged as a promising technology to control pests and pathogens within agricultural systems and has been demonstrated to reduce the aggressiveness of plant pathogenic filamentous fungi both through host-induced gene silencing (HIGS) and exogenous applications of double-stranded (dsRNA) onto the surface of plants (Yin et al., 2010; Andrade et al., 2016; Koch et al., 2016; Wang et al., 2016; Tiwari et al., 2017; Gilbert et al., 2018; Wagner et al., 2020). Spray-induced gene silencing (SIGS), involving ectopic application of dsRNAs targeting pathogen genes, has the added benefit of circumventing the genetic manipulation of plants, and thus further avoiding the regulatory expenses and public concern surrounding genetically modified organisms (GMOs). SIGS has been tested on different plant-pathosystems and was effective against various groups of plant pests and pathogens (Dubrovina and Kiselev, 2019) including viruses (Carbonell et al., 2008), insects (Li et al., 2015; Gogoi et al., 2017), and fungi (Koch et al., 2016; Wang et al., 2016; McLoughlin et al., 2018; Hu et al., 2020). Furthermore, applied dsRNAs targeting pathogen genes can be taken up and processed by plants (Koch et al., 2016; Mitter et al., 2017). For instance, spraying dsRNAs targeting a gene involved in Fusarium graminearum ergosterol biosynthesis on detached barley leaves inhibits fungal growth both at the application site and on distal parts of leaves (Koch et al., 2016). Due to the wide array of target options, the specificity required for RNA binding, and advances in protective materials such as clay nanosheets to adhere and protect dsRNAs on the plant surface (Mitter et al., 2017), HIGS and exogenous applications of dsRNA constitute promising crop protection technologies.
The efficacy of HIGS and SIGS approaches against fungal pathogens is dependent, in part, on the nature of the selected target sequence. Genes encoding pathogenicity determinants can constitute ideal targets of such approaches, because they are essential for pathogens to cause disease. In previous studies, silencing of pathogenicity genes by HIGS in multiple pathogens suppressed disease development (Panwar et al., 2013; Pliego et al., 2013; Yin et al., 2014, 2015; Shanmugam et al., 2017). In S. sclerotiorum, oxalic acid (OA) is a key factor governing its pathogenic success. OA can contribute to pathogenesis in a variety of ways, including the manipulation of the redox status of the host (Williams et al., 2011; Ranjan et al., 2018), the induction of programmed cell death (Kim et al., 2008; Kabbage et al., 2013), and the activation of degradative enzymes (Dutton and Evans, 1996). Thus, the use of HIGS to target OA production could provide a valuable tool to manage SSR and reduce fungicide use. Previous work has targeted a chitin synthase in S. sclerotiorum and reduced disease severity through HIGS in transgenic tobacco plants (Andrade et al., 2016). However, it is important to test this approach on specific crops of interest, identify additional useful targets that can be stacked to ensure durable resistance in a particular host, and importantly, identify and target unique sequences such as virulence factors to avoid sufficiently conserved genes across taxa to prevent off-target effects, including within beneficial fungal communities.
Herein, we target a S. sclerotiorum oxaloacetate acetylhydrolase (Ssoah1) gene that encodes an enzyme that converts oxaloacetate to oxalate and acetate. The deletion of this gene in S. sclerotiorum abolished OA accumulation and resulted in restricted lesions in which the infectious hyphae gradually lose viability on tomato, soybean, and Arabidopsis (Liang et al., 2015 and Xu et al., 2015). Additionally, Ssoah1 expression appears to play a particularly critical role during infection of soybean, compared to other S. sclerotiorum hosts (Westrick et al., 2019 and Xu et al., 2015). Using a combination of confocal microscopy, gene expression studies, small RNA (sRNA) sequencing, and virus-induced gene silencing, we show that S. sclerotiorum mycelia take up environmental double-stranded and small interfering (ds/siRNAs) targeting Ssoah1, leading to the downregulation of Ssoah1 gene expression both in vitro and in planta. In planta, silencing assays were conducted using bean pod mottle virus (BPMV)-mediated HIGS. Importantly, our HIGS assays significantly suppressed S. sclerotiorum development and disease onset. This work demonstrates that RNAi strategies provide an additional tool for improving soybean resistance to S. sclerotiorum and managing SSR.
Materials and Methods
In vitro Synthesis of dsRNA and siRNA
RNA was extracted from actively growing S. sclerotiorum mycelia using a Maxwell RSC Plant RNA Kit (Promega) and cDNA was synthesized using a ProtoScript First Strand cDNA Synthesis Kit (New England Biolabs) following manufacturers' instructions. A 365-base pair fragment corresponding to Ssoah1 was amplified and added with a T7 promoter sequence by PCR using specific primers, T7_OA_F and T7_OA_R (Supplementary File 1). PCR reactions (50 μl) were set up using a KAPA HiFi PCR Kit (Roche Life Science) following the manufacturer's instructions and placed in a thermal cycler using the following conditions: 95°C for 3 min, 35 cycles of: 98°C for 30 s, 62°C for 15 s, and 72°C for 30 s, and last-step extension at 72°C for 30 s. The resulting PCR product was purified using a Zymoclean Gel DNA Recovery Kit (Zymo Research) and further used for in vitro transcription using a MEGAscript RNAi Kit (Life Technologies) according to the manufacturer's instructions. Ssoah1-siRNA was processed from Ssoah1-dsRNA using ShortCut RNase III (New England Biolabs) and subsequently purified by ethanol precipitation protocol according to the manufacturer's specifications. To perform fluorescent microscopy, Ssoah1-dsRNA and siRNA that were processed from dsRNA were independently labeled with Cy3 dye using a Silencer siRNA Labeling Kit (Life Technologies). Both non-labeled and labeled Ssoah1-dsRNA and siRNA were subsequently quantified by spectrophotometer (NanoDrop 2000, Thermo Scientific) following the protocols provided in the user guides for the MEGAscript RNAi Kit and Silencer siRNA Labeling Kit. GFP-dsRNA was used as a negative control, and a 353 bp-fragment of the GFP gene was also selected for in vitro transcription of GFP-dsRNA using procedures mentioned above and the primer pair, T7_GFP_F and T7-GFP-R (Supplementary File 1).
Detection of Fluorescently Labeled ds/siRNAs in Sclerotinia sclerotiorum
A 5-mm hyphal plug of S. sclerotiorum from the leading edge of a 2-day old fungal culture was placed on a cellophane layer on the surface of a Petri plate containing potato dextrose agar (PDA). After 2 days, 1 ml of potato dextrose broth (PDB) was added to the Petri plate, and fungal mycelia were collected from cellophane and transferred into a 1.5 ml Eppendorf tube. After centrifugation at 4,000 rpm for 5 min, the supernatant was removed, and mycelia were resuspended in 200 μl of PDB. Next, 4 μg of Cy3-labeled Ssoah1-dsRNA/siRNA were mixed into the PDB containing S. sclerotiorum mycelia and incubated on a shaker at room temperature at 100 rpm for 13 h. Cy3-Ssoah1-dsRNA/siRNA treated S. sclerotiorum mycelia were washed by adding 800 μl of sterile H2O, and mycelia were then collected by centrifugation, at 4,000 rpm for 5 min, for microscopy. To produce protoplasts from Cy3-Ssoah1-dsRNA/siRNA treated fungal material, mycelia collected from cellophane on PDA were suspended in 1 ml PDB and transferred into a 25 ml flask. After 18 h of incubation with 20 μg dsRNA or siRNA, mycelia were collected by centrifugation, and protoplasts were processed by adding 2.4 ml protoplasting buffer containing 10 mg/ml lysing enzyme (modified from Rollins, 2003). To determine if S. sclerotiorum takes up dsRNAs non-selectively, 4 μg of Cy3-labeled siRNA with no targets in S. sclerotiorum (Silencer Cy3-labeled Negative Control No. 1 siRNA, Thermo Fisher Scientific) was also applied to the fungus exogenously. Water, applied to S. sclerotiorum, was used as negative control for microscopy. Mycelia and protoplasts were visualized by confocal microscopy using a Zeiss LSM 710.
RNA Extraction, cDNA Synthesis, and RT-qPCR
A 5-mm S. sclerotiorum hyphal plug was placed into a 25-ml flask with 7 ml PDB at pH 8. Non-labeled Ssoah1-dsRNA were applied to PDB to make a final dsRNA concentration of 500 ng/ml, and the flasks were kept at room temperature on a shaker at 100 rpm. The application of GFP-dsRNA in PDB with S. sclerotiorum hyphal plugs, maintained in the same conditions as above, was considered a negative control for this experiment. Mycelia were collected at 96 h after placement in the PDB. RNA extraction and cDNA synthesis were performed as described previously. The transcriptional level of Ssoah1 was determined by RT-qPCR using synthesized cDNA as a template. Reactions (20 μl) were set up using SsoFast EvaGreen Supermix (Bio-Rad) according to the manufacturer's instructions. Quantitative PCR (qPCR) reactions were conducted using the Bio-Rad CFX96 Real-Time System and following program: 95°C for 30 s, followed by 40 cycles of: 95°C for 5 s, and 58°C for 5 s. Melt curves were generated by heating from 65°C to 95°C with 0.5°C increments. Relative expression of Ssoah1 to S. sclerotiorum H3 were calculated as 2−ΔCt (Livak and Schmittgen, 2001) prior to calculating expression relative to the GFP negative control. Experiments were repeated three times with at least three biological replicates for each repetition. Primers used for qPCR are listed as follows: Ssoah1_F1, Ssoah1_R1, H3_F1, and H3_R1 (Supplementary File 1). In conjunction with expression assays, PDB was collected from 7 ml cultures containing 500 ng/ml of dsRNA or GFP at 24, 48, and 72 h to analyze oxalate production using the method described in Willbur et al. (2017) (Supplementary File 2). Plates were read at 590 nm using a Bio-Rad iMark (Bio-Rad Laboratories) microplate reader.
Construction of BPMV Silencing Vectors
In order to silence the fungal gene Ssoah1, a 365 base pair sequence (Supplementary File 3) was selected within the mRNA of S. sclerotiorum strain 1980 (GenBank Accession XM_001590428). Total RNA was extracted from S. sclerotiorum using Trizol, and cDNA was synthesized via reverse transcription-PCR (RT-PCR) using an AMV first strand cDNA synthesis kit (New England Biolabs, Catalog # E6550). Through PCR with a high-fidelity DNA polymerase (KAPA HiFi®, Kapa Biosystems), BAMH1 restriction sites (to form “sticky ends”) were incorporated onto the double-stranded cDNA using specific primers (BAMH1_OA_F and BAMH1_OA_R, Supplementary Table 1), and the following PCR conditions: 95°C for 2 min followed by 35x at 98°C for 20 s, 68°C for 30 s, 72°C for 15 s, and a final extension at 72°C for 4 min. After digestion and gel purification (QIAquick Gel Extraction Kit®, QIAGEN), the construct was ligated into the viral vector RNA2 plasmid, pBPMV-IA-V1 to form pBPMV-OA (Zhang et al., 2010). The vector plasmids were then transformed into DH5α competent cells using 5 μl of the purified ligation product per 60 μl of competent cells, a 30 min ice incubation, 45 s heat shock in a 42°C water bath, and incubation in 500 μl of Luria broth (LB) for 1 h at 37°C. Colony PCR with MyTaqTM Red Mix (Bioline) was used to confirm successfully transformed clones. After an initial step for cell breakage at 95°C for 6 min, 35 cycles (95°C for 15 s, 68°C for 15 s, and 72°C for 10 s) were performed to confirm successful transformations, and colonies from a mini preparation were sent for sequencing to confirm antisense insertion. Using glycerol stocks, midi preparations were conducted (Fast Ion Plasmid Midi Kit®, IBI Scientific) for subsequent biolistic inoculations.
Plant Inoculation With BPMV Vectors
The modified system, used in this work and developed by Zhang et al. (2010), uses bean pod mottle virus (BPMV) which is bipartite. The two viral segments consist of a 6 kb RNA1 segment (pBPMV-IA-R1M) encoding viral proteins needed for genome replication and a 3.6 kb RNA2 segment (pBPMV-IA-V1) that encodes viral proteins needed for capsid assembly and movement and contains the cloning site where the silencing construct was inserted. As described in the protocol by Whitham et al. (2016), the RNA1 and RNA2 segments were co-bombarded into 10-day old Williams 82 (PI 518671) seedlings using gold particles coated with plasmid DNA. Two inoculations were made per plant, with one on each unifoliate leaf. Prior to inoculations, unifoliate seedlings were forced to etiolate in the dark for 24 h. After inoculations, plants were sprayed with water and kept in plastic bags for 24 h to maintain humidity. The inoculated plants were placed in a growth chamber at 22°C during the day and 20°C at night on a 16 h photoperiod. Soil was checked daily to determine if water was needed, and the plants were fertilized once weekly with Miracle-Gro® (Scotts Miracle-Gro Co.).
Successful infection of Williams 82 soybean plants was confirmed visually, as phenotypes resembled those of plants infected with BPMV, including rugosity and mottling (Supplementary File 4). Further confirmation using RT-PCR was also performed (Supplementary File 5), using the primers for silencing construct development, and one corresponding to the RNA2, VS_R2_F, and VS_R2_R (Supplementary File 1). Williams 82 soybean plants confirmed to contain the viral construct were lyophilized and stored at −80°C to be used in subsequent rub inoculations of more soybean plants. These inoculations were performed by grinding 50 mg of lyophilized tissue in 50 mM potassium phosphate buffer and rubbing the plant serum onto the soybean variety Traff (PI 470930), as described in Whitham et al. (2016). The Traff-BPMV HIGS system was previously validated by Ranjan et al. (2018). Traff has more tolerance to BPMV than Williams 82 and Sclerotinia and BPMV are able to co-infect in this system. Observations suggest that BPMV-infected Williams 82 inhibits S. sclerotiorum infection (Grau, data not shown). Viral symptoms of mottling and rugosity were apparent but milder in Traff than in Williams 82 (Supplementary File 6). RNA was extracted from the third trifoliates of plants with viral symptoms, from one experimental replicate, and plants were confirmed to contain the silencing constructs via RT-PCR and sequencing.
Characterization of Processed sRNAs
To confirm the processing of siRNA from pBPMV-OA, a high concentration of total RNA was extracted from three leaves in three biological replicates from the pBPMV-OA and pBPMV-EV- containing plants. Leaves for RNA extractions were collected from the variety Traff, prior to inoculating soybeans symptomatic for BPMV with S. sclerotiorum in one experimental replicate of the expression experiment. The third trifoliate leaves were immediately frozen in liquid nitrogen. Two grams of leaf tissue per biological replicate was ground into fine powder using liquid nitrogen in a mortar and pestle and used in RNA extractions. A phenol chloroform extraction with a LiCl purification was performed as described in Aragão et al. (2013) to yield high RNA concentrations of 1.2–2.8 μg/μl. A cDNA library was constructed by Novogene (CA, USA) using a Small RNA Sample Pre Kit, and Illumina sequencing was conducted according to company workflow, using 20 million reads. Raw data were filtered for quality as determined by reads with a quality score >5, reads containing N <10%, no 5′ primer contaminants, and reads with a 3′ primer and insert tag. The 3′ primer sequence was trimmed and reads with a poly A/T/G/C were removed.
Pathogen Aggressiveness Assays
In order to assess the aggressiveness of S. sclerotiorum on plants containing pBPMV-OA and plants containing pBPMV-EV constructs, BPMV-symptomatic Traff plants were inoculated using the cut petiole method as detailed by Peltier and Grau (2008), and lesions were measured 24–120 h post inoculation (HPI). Five pots, completely randomized per biological replicate and treatment, were thinned to one to three symptomatic plants per pot. Viral symptoms were apparent by the V3-V5 (vegetative stage at which the plant had three to five trifoliates) growth stage.
At the V4-V5 growth stage, the third trifoliate was excised at 2 cm and inoculated with S. sclerotiorum strain 1980. All cultures were grown in the same manner prior to inoculations. The isolate was obtained from sclerotia harvested from plants grown in a growth chamber, surface disinfested for 1 min in 95% ethanol and 1 min in 10% concentrated household bleach (8.25% NaClO, prior to dilution), and grown on standard, 15-mm deep Petri plates containing PDA, before sub-culturing onto 20-mm Petri plates with a thicker depth of PDA prior to inoculations. After 3 days of growth on these latter plates, a plug from the leading edge of mycelia was cut using an inverted one-ml pipette tip. The inverted pipette tip with agar plug, was then slid onto the excised petiole, which was cut to 2 cm prior to inoculation. Lesions were measured with digital calipers 24–120 HPI. These data were used to compare lesion size, across time points, and used to calculate the area under the disease progress curve transformation (Shaner and Finney, 1977). This experiment was repeated three times.
Plant Expression Assays
Similarly, for expression assays, Traff soybean plants were challenged with S. sclerotiorum strain 1980 using the cut petiole technique at the third petiole when plants were at the V4 or V5 growth stage. Expression assays were conducted in experiments independent of the pathogen aggressiveness assay. However, the 120 h time point used in expression assays included plant samples from the pathogen aggressiveness assay, since lesions in the aggressiveness assay were last measured at 120 h. In order to understand whether expression decreased in plants containing the silencing construct compared to plants with pBPMV-EV, 6 cm of stem tissue (3 cm above and 3 cm below the inoculation site) was collected 48–120 HPI and immediately frozen in liquid nitrogen. One-to-three stems per pot were combined for RNA extraction and three biological replicates were used for extractions. Expression experiments were repeated three times. Tissues were ground in liquid nitrogen using a mortar and pestle, and RNA was extracted using a Maxwell® RSC Plant RNA Kit. Two primer pairs were used for RT-qPCR to determine whether expression was reduced in plants containing pBPMV-OA (as described previously for expression analyses with exogenous applications of RNA). One primer pair (Ssoah1_F1 and Ssoah1_R1) corresponded to Ssoah1, outside of the region used for the silencing precursor, as to not amplify sequences from the silencing vector, and another corresponded to the S. sclerotiorum endogenous control, Histone H3 (H3_F1 and H3_R1). Relative expression of Ssoah1 to S. sclerotiorum H3 were calculated as 2−ΔCt (Livak and Schmittgen, 2001) prior to calculating expression relative to plants containing the empty vector negative control, pBPMV-EV.
Statistical Analyses
Expression, AUDPC, and lesion size differences were evaluated using a mixed-model analysis of variance (ANOVA) using PROC GLIMMIX in the SAS statistical software package (v 9.4, SAS Institute, Inc. Cary, NC, United States). Significance was reported at the α = 0.05 significance level. For the lesion size differences in aggressiveness assays, a lognormal distribution and compound symmetry covariance structure were used in the model for data analyses.
Results
Uptake of External RNAs by Sclerotinia sclerotiorum
Effective HIGS strategies are reliant on the ability of fungal organisms to uptake external ds/siRNAs. Accordingly, we tested whether such molecules can translocate into S. sclerotiorum hyphae. Double-stranded RNAs (dsRNAs) corresponding to a S. sclerotiorum Ssoah1 were transcribed in vitro, processed into small interfering RNAs (siRNAs) and labeled with Cy3 fluorescent dye. Fluorescent signals were detected in S. sclerotiorum hyphae, after a two-day growth period in potato dextrose broth (PDB) supplemented with Cy3-labeled Ssoah1-dsRNA/siRNAs. Clear, fluorescent signals were detected in Cy3-Ssoah1-dsRNA/siRNA treated samples using fluorescent confocal microscopy, but not in the negative control (Figure 1). To further confirm that RNAs were within the fungal cellular space, rather than on the surface of the hyphae, protoplasts were generated from Cy3-Ssoah1-dsRNA/siRNA treated mycelia. Similarly, fluorescent signals were clearly visible within protoplasts, indicating their translocation across cell walls and plasma membranes (Supplementary File 7). The uptake of environmental RNAs appears to occur readily in S. sclerotiorum, as Cy3-labeled siRNAs with no specific targets in this fungus were also observed in fungal mycelia and protoplasts (Figure 1 and Supplementary File 7). No fluorescence was observed in water-treated S. sclerotiorum mycelia or protoplasts (Figure 1 and Supplementary File 7). Taken together, these results indicate that environmental RNAs can freely translocate into S. sclerotiorum cells, and can thus potentially be used to target gene expression within the pathogen.
Figure 1. External RNAs are taken up by Sclerotinia sclerotiorum. Treatments were added to potato dextrose broth (PDB) containing 2 day old S. sclerotiorum mycelia originally grown on potato dextrose agar (PDA). Fluorescent signals were detected using confocal fluorescent microscopy. (A) An application of water, as a negative control, yielded no signal. (B) The positive control, application of Cy3-labeled-nocoding siRNA (Cy3-NC-siRNA, with no specific targets in S. sclerotiorum), (C) 4 μg application Cy3-labeled-Ssoah1-dsRNA, and (D) application of 4 μg of Cy3-labeled-Ssoah1-siRNA (targeting Ssoah1 in S. sclerotiorum) produced fluorescent signals that were detected 13 h post inoculation. Scale bars = 20 μm.
External Application of RNAs Targeting Ssoah1 Reduced Transcript Levels in vitro
We next tested whether the external application of dsRNAs leads to the silencing of the target gene. In three independent experiments, actively growing S. sclerotiorum hyphal plugs were placed in potato dextrose broth (PDB) containing in vitro transcribed Ssoah1-dsRNAs; GFP-dsRNA treatment served as a negative control. Transcript levels of Ssoah1 were assessed by RT-qPCR at 96 h post-inoculation. Compared to the GFP-dsRNA controls, Ssoah1 gene expression was markedly reduced (by >60%) in mycelia exposed to Ssoah1-dsRNAs (Figure 2) (p < 0.01), suggesting the uptake of Ssoah1-dsRNAs resulted in the downregulation of Ssoah1 transcripts in S. sclerotiorum and thus potentially depriving the fungus from a key virulence factor during infection. As expected, the downregulation of Ssoah1 also corresponds to a marked decrease in OA production by the fungus (Supplementary File 2).
Figure 2. External uptake of Ssoah1-dsRNA decreased Ssoah1 transcript levels in Sclerotinia sclerotiorum. A concentration of 500 ng/mL of Ssoah1-dsRNA were applied to potato dextrose broth (pH 8) and a potato dextrose agar plug containing mycelia of S. sclerotiorum was placed in the broth. Treatment with a same concentration of GFP-dsRNA served as negative control. Mycelia were collected 96 h after addition to PDB and the transcript level of Ssoah1 were analyzed. Relative expression of Ssoah1 was normalized to S. sclerotiorum H3 and relative to the GFP negative control. Ssoah1 transcript levels were significantly reduced across three experimental replicates (*p < 0.01). Error bars represent ± standard errors of three (for R1 and R2) or for (for R3) biological replicates across three experimental repetitions.
Overall, our in vitro assays show that S. sclerotiorum uptakes environmental RNAs targeting Ssoah1, and this uptake leads to reduced expression of the target gene and OA synthesis. Accordingly, we reasoned that dsRNA/siRNAs generated as part of a HIGS strategy targeting a vital virulence factor might hinder S. sclerotiorum pathogenic development in soybean.
Accumulation of siRNAs Corresponding to Ssoah1 in Soybean Mediated by BPMV
With evidence that S. sclerotiorum could uptake exogenous siRNAs, we evaluated whether plants could express silencing constructs for uptake by S. sclerotiorum through HIGS. To express siRNAs targeting Ssoah1, a 365 bp Ssoah1 fragment was cloned into the viral vector pBPMV-IA-V1 (Zhang et al., 2010) after the viral translation stop codon to create pBPMV-OA. Soybeans were inoculated with pBPMV-OA and an empty vector control (pBPMV-EV, negative control) using a combination of particle bombardment and rub inoculation. Small RNAs (sRNAs), sequenced from plant total RNA, were mapped to the 365 bp region carried by the pBPMV viral vector (Figure 3). Of the sRNAs profiled, 21 nt sequences were particularly abundant, with over 10,000 copies (reads per millions X.0235) estimated to be contained in leaf extracts tested (Figure 3A). This size corresponds to the 20–25 nt size characteristic of siRNAs, knowing that 21 nt siRNAs are the most efficient at mRNA degradation (Elbashir et al., 2001). These sRNAs mapped to various regions across the 365 bp cDNA silencing construct (Figures 3B,C and Supplementary File 3). The location distribution was similar in all three biological replicates from which RNA was sequenced (Figure 3 and Supplementary Files 8, 9). As expected, reads corresponding to the silencing target were only present in the plants containing the pBPMV-OA vector, and not in plants containing the empty vector, except for five, negligible reads in one of the three samples sequenced for sRNAs. Therefore, HIGS plants were confirmed to generate abundant siRNAs corresponding to the S. sclerotiorum gene, Ssoah1.
Figure 3. Total sRNAs were profiled from Traff soybean plants bombarded with pBPMV -OA and pBPMV-EV. The size distribution of sRNA reads <27 nt and corresponding to the 365 nt target within Ssoah1 were normalized to a count, using the conversion factor of 0.0235 X reads per million (RPM). (A) The size distribution of Ssoah1-dsRNA-derived sRNAs corresponding to Ssoah1 from a pBPMV-OA single plant are presented for sRNAs ranging from 18 to 27 nt. The 5′ end of 18–27 nt sRNAs from plants containing (B) pBPMV-OA and (C) pBPMV-EV vectors were aligned to the 365 nt sequence cDNA sequence used for silencing construct development. Antisense reads are presented in orange and sense reads are presented in blue.
Targeted Silencing of Ssoah1 in planta Using BPMV-Mediated HIGS
We tested whether expression of the target gene, Ssoah1, could be reduced in planta using a modified BPMV vector. Plant containing pBPMV-OA and empty vector controls were challenged with S. sclerotiorum strain 1980 (ATCC18683) in a time-course experiment. Stem tissues were collected at 48, 72, 96, and 120 HPI for RNA extractions and subsequent RT-qPCR. Across three independent experimental repetitions, the expression of Ssoah1 was significantly reduced in plants inoculated with the pBPMV-OA silencing vector compared to the empty vector control (Figure 4). The most significant reduction in expression was observed at 48 HPI (60% reduction, p < 0.01) and 96 HPI (42% reduction, p = 0.02; Figure 4) and all-time points trended toward a reduction in Ssoah1 expression. Although differences were not statistically significant at 72 and 120 HPI. These results suggest that BPMV-mediated HIGS targeting Ssoah1, effectively reduced the expression of this S. sclerotiorum gene in our time course.
Figure 4. Transcript levels of the target gene, Ssoah1, were decreased in expression assays. Plants containing virus-induced gene silencing constructs, pBPMV-OA and empty vectors, pBPMV-EV were challenge with S. sclerotiorum strain 1980, using a cut petiole technique. RNA was extracted from stem tissue from three biological replicates and three experimental repetitions. Relative expression of Ssoah1 was normalized to S. sclerotiorum H3 and relative to the pBPMV-EV. Expression in pBPMV-OA treated plants tended to decrease for all time points, however, the difference was significant at 48 h (*p < 0.01) and 96 h (p = 0.02).
Disease Development Was Impeded in BPMV-Mediated HIGS Plants
To test whether the reduced expression of Ssoah1 resulted in reduced pathogen aggressiveness, pBPMV-OA was introduced into the soybean variety, Traff. Symptomatic Traff seedlings were similarly challenged with S. sclerotiorum strain 1980 (ATCC18683) using a cut petiole method (Peltier and Grau, 2008; Ranjan et al., 2018), and lesions were measured 24-120 HPI. Visual differences in lesions were apparent, with plants containing pBPMV-EV having large, girdling lesions 96 HPI compared to markedly restricted lesions in pPBMV-OA plants (Figure 5A). Significant differences in lesion size were detected at 72 HPI and onward (Figure 5B, p < 0.01). Additionally, the area under the disease progress curve (AUDPC) in plants containing the silencing construct was reduced by 73% compared to pBPMV-EV plants (Figure 5C, p < 0.01), indicating a delay and reduction in SSR disease development. Reduced pathogen development using BPMV-mediated HIGS in soybean challenged with S. sclerotiorum suggest that such a strategy is effective in controlling SSR.
Figure 5. Pathogen aggressiveness was reduced in virus induced gene silencing (VIGS) assays. Plants containing VIGS silencing constructs (pBPMV-OA) and empty vectors (pBPMV-EV) were challenged with S. sclerotiorum strain 1980, using a cut petiole technique. Lesions were measured 1–5 days post inoculation (DPI). Data are the result of five biological replicates and three experimental repetitions. (A) Lesion development was delayed, and lesions were smaller in plants containing pBPMV-OA, while plants containing pBPMV-EV often showed girdling lesions at 96 h post inoculation (HPI). (B) Lesion size was reduced after two DPI, with an overall reduced lesion size (*p < 0.01). (C) The overall area under the disease progress curve was lower in plants containing pBPMV-OA (*p < 0.01).
Discussion
The first application of in planta HIGS as a disease control strategy was patented in 2008 for control of pathogenic fungi and oomycetes in Glycine max and Arabidopsis (Roberts et al., 2008). Since, various strides have been made in the application of HIGS (Peltier and Grau, 2008; Nowara et al., 2010; Tinoco et al., 2010; Yin et al., 2010; Panwar et al., 2013; Ghag et al., 2014; Andrade et al., 2016; Koch et al., 2016; Dou et al., 2020; Mahto et al., 2020) and SIGS (Wang et al., 2016; McLoughlin et al., 2018; Koch et al., 2019) in the control of a range of phytopathogens. The aforementioned works demonstrate that RNAi can be an effective gene silencing strategy in many plant-fungal interactions. Extending such methods to other fungal pathosystems could constitute an exciting alternative to chemical management, particularly in crops where breeding efforts have failed to produce highly resistant commercial cultivars. This study targets an essential pathogenicity factor of S. sclerotiorum, oxalic acid (OA) to enhance the resistance of soybean to Sclerotinia stem rot (SSR). Oxalic acid is an important virulence factor in S. sclerotiorum, where OA deficient mutants show reduced colonization and disease development (Liang et al., 2015; Xu et al., 2015). Accordingly, we targeted Ssoah1, a gene encoding an oxaloacetate acetylhydrolase enzyme that catalyzes the last step in OA biosynthesis, and the deletion of which abolishes OA production (Liang et al., 2015; Xu et al., 2015). Several lines of evidence presented within this study are consistent with the following conclusions: (i) environmental dsRNAs/siRNAs translocate into S. sclerotiorum hyphae; (ii) transcripts encoding Ssoah1 are markedly reduced when S. sclerotiorum is exposed to the corresponding dsRNAs in vitro and in planta; and (iii) RNAi of Ssoah1 using HIGS can reduce pathogen aggressiveness and SSR disease development.
We first determined whether S. sclerotiorum can uptake environmental RNAs, as uptake of dsRNA is not a ubiquitous feature of fungi. While exogenous dsRNA uptake occurs in the closely related Botrytis cinerea, recent findings indicate that the wheat pathogen Zymoseptoria tritici does not uptake environmental dsRNA despite encoding key components of the RNAi pathway (Wang et al., 2016; Kettles et al., 2019). After applying fluorescently labeled dsRNA/siRNA to fungal cultures, we detected fluorescent signals within fungal hyphae treated with both Ssoah1 derived-dsRNA/siRNAs and non-targeting siRNA (Figure 1), indicating that S. sclerotiorum readily uptakes these molecules in a non-specific manner.
While the mechanism by which RNA molecules are able to translocate across fungal cell walls and membranes is yet to be determined, it is clear that the non-selective movement of these molecules occurs in the predominantly necrotrophic pathogen S. sclerotiorum. Within fungi, these translocated dsRNAs are processed to siRNA after uptake and both dsRNA and siRNA efficiently silence gene expression (Dang et al., 2011 and Koch et al., 2016). Though exceptions have been noted, ds/siRNA uptake is likely a common occurrence in the fungal kingdom considering the successful utilization of HIGS against several fungal species (Nowara et al., 2010; Tinoco et al., 2010; Yin et al., 2010; Panwar et al., 2013; Pliego et al., 2013; Ghag et al., 2014; Andrade et al., 2016; Koch et al., 2016; Cooper and Campbell, 2017; Dou et al., 2020; Hu et al., 2020; Mahto et al., 2020). Thus, RNA spray regimes or HIGS has the potential to be utilized across many fungal pathogens including S. sclerotiorum. Interestingly, these fungal targeted RNAs also translocate, and are amplified within plant cells as well. Long dsRNA and dsRNA-derived siRNA efficiently translocate in barley (Koch et al., 2016). Song et al. (2018) also reported the successful in planta amplification of siRNA targeting Fusarium in wheat through plant RdRp, providing evidence for the potential of sustained dsRNA/siRNA-based SIGS in a field setting. This may also alleviate concerns about the stability of such molecules in nature and may confer prolonged protection under field conditions. Spray-induced gene-silencing studies are needed to evaluate the effectiveness of exogenously-applied dsRNAs for silencing Ssoah1 at local and distal locations, relative to the application site, in S. sclerotiorum hosts and to evaluate the effectiveness of SIGS targeting Ssoah1 to manage disease in a field setting.
To confirm that these dsRNA/siRNAs were not only absorbed but actively utilized in gene silencing, Ssoah1 transcripts in S. sclerotiorum were analyzed following treatment with corresponding dsRNAs in vitro. Oxalic acid production in S. sclerotiorum is a pH-responsive process and Ssoah1 transcripts are concordantly induced in alkaline conditions (Rollins, 2003; Kim et al., 2008). To ensure that Ssoah1 was highly expressed and a reduction of transcripts in S. sclerotiorum would be observable after taking up Ssoah1-dsRNA, the treated liquid media was alkalized to a pH of 8. As anticipated, uptake of Ssoah1-dsRNA did occur and resulted in significant reduction of Ssoah1 mRNA level at 96 h after dsRNA treatment (Figure 2). These results were encouraging, and clearly show that the translocation of environmental RNAs in S. sclerotiorum significantly affect transcript abundance of corresponding sequences.
In planta assays were performed to evaluate the use of HIGS targeting Ssoah1 to reduce the aggressiveness of S. sclerotiorum. We took advantage of our working virus-based HIGS system (Zhang et al., 2010; Ranjan et al., 2018) that uses BPMV as a vehicle to produce dsRNAs corresponding to our fungal target. Soybeans were inoculated with a BPMV vector containing an Ssoah1 fragment (pBPMV-OA), and the presence of corresponding sRNAs was determined via RNA sequencing. Indeed, sRNAs aligning to Ssoah1 were found within the expected range of siRNAs (21-24 nt, Figure 3A) in pBPMV-OA inoculated plants indicating successful processing by the soybean RNAi machinery. The sRNAs corresponding to Ssoah1 were not detected in the empty vector (pBPMV-EV) control plants. Importantly, a marked reduction in disease symptoms and reduced AUDPC was observed in pBPMV-OA soybeans compared to control empty vector plants (Figure 5). Small lesions did develop on plants expressing Ssoah1-dsRNAs, but in most replicates, these lesions were not girdling the stems at 96 HPI, as occurred in plants inoculated with pBPMV-EV. This reduction in symptom development was accompanied by reduced transcript levels of Ssoah1, indicating that HIGS is an effective strategy for reducing expression of S. sclerotiorum genes in soybean and conferring resistance to this pathogen. Future work will focus on the generation and evaluation of stable, transgenic soybean plants that use HIGS to target virulence factors, such as Ssoah1, in S. sclerotiorum.
A similar, previous approach, aimed to reduce the aggressiveness of S. sclerotiorum by targeting the structural gene chitin synthase in tobacco (Andrade et al., 2016). This study, however, is the first case of HIGS being used to target a pathogenicity rather than developmental factors of S. sclerotiorum. Because Ssoah1 orthologs are primarily present in other plant pathogenic fungi, the risk of off-target effects from this silencing construct on beneficial and non-pathogenic fungi is low. Targeting pathogenicity-related genes with RNAi in cereal crops has provided effective disease resistance to Puccinia graminis and the targeting of Ssoah1 and other S. sclerotiorum genes may provide similar resistance in soybean (Panwar et al., 2018). Ongoing research efforts in our lab aim to identify S. sclerotiorum specific virulence determinants that can be utilized in this manner. This is particularly important, since targeting structural components or housekeeping genes did not yield positive results (data not shown).
Our findings pave the way for RNAi-based approaches that target virulence factors to achieve enhanced resistance to S. sclerotiorum in soybean and potentially other valuable field crops. Virus-based HIGS using pBPMV, as developed by Zhang et al. (2010), can also be deployed in functional studies to both understand the role of S. sclerotiorum genes in infecting soybean hosts and to evaluate host factors involved in S. sclerotiorum pathogenic development (Ranjan et al., 2018). As a novel tool in SSR management, RNAi has the potential to reduce the use of chemical control in soybean systems and to close the gap of incomplete resistance observed in commercial and public soybean cultivars (Kim and Diers, 2000; Hoffman et al., 2002; McCaghey et al., 2017; Conley et al., 2018; Huzar-Novakowiski and Dorrance, 2018).
Data Availability Statement
Illumina sequencing data of small RNA from plants containing empty vectors and silencing vectors are publicly available through NCBI as accession GSE180131, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE180131.
Author Contributions
MM, DSh, SC, DSm, and MK designed all the experiments. MM and AR identified the silencing target. MM designed the in planta and exogenously applied silencing constructs and related primers. MM, JK, and AL generated pBPMV-OA silencing constructs, conducted the in planta assays, and analyzed data with assistance from DSm. DSh synthesized the exogenously applied dsRNA, conducted the in vitro experiments, and conducted microscopy and imaging. SW shared pBPMV vectors. MM, DSh, BW, DSm, and MK wrote the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was generously supported by the Wisconsin Soybean Marketing Board and the United States Department of Agriculture's National Sclerotinia Initiative (ARS Agreement Number: 58-3060-8-023). Graduate student, MM was also supported by the Department of Plant Pathology at the University of Wisconsin-Madison and The Storkan Hanes McCaslin Foundation.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.677631/full#supplementary-material
References
Allen, T. W., Bradley, C. A., Sisson, A. J., Byamukama, E., Chilvers, M. I., Coker, C. M., et al. (2017). Soybean yield loss estimates due to diseases in the United States and Ontario, Canada, from 2010 to 2014. Plant Health Prog. 18, 19–27. doi: 10.1094/PHP-RS-16-0066
Andrade, C. M., Tinoco, M. L. P., Rieth, A. F., Maia, F. C. O., and Aragão, F. J. L. (2016). Host-induced gene silencing in the necrotrophic fungal pathogen Sclerotinia sclerotiorum. Plant Pathol. 65, 626–632. doi: 10.1111/ppa.12447
Aragão, F. J. L., Nogueira, E. O. P. L., Tinoco, M. L. P., and Faria, J. C. (2013). Molecular characterization of the first commercial transgenic common bean immune to the Bean golden mosaic virus. J. Biotechnol. 166, 42–50. doi: 10.1016/j.jbiotec.2013.04.009
Boland, G. J., and Hall, R. (1986). Growthroom evaluation of soybean cultivars for resistance to Sclerotinia sclerotiorum. Can. J. Plant Sci. 66, 559–564. doi: 10.4141/cjps86-075
Carbonell, A., Ángel-Emilio Martínez de Alba, A. E., Flores, R., and Gago, S. (2008). Double-stranded RNA interferes in a sequence-specific manner with the infection of representative members of the two viroid families. Virology 371, 44–53. doi: 10.1016/j.virol.2007.09.031
Chun, D., Kao, L. B., Lockwood, J. L., and Isleib, T. G. (1987). Laboratory and field assessment of resistance in soybean to stem rot caused by Sclerotinia sclerotiorum. Plant Dis. 71, 811–815. doi: 10.1094/PD-71-0811
Conley, S., Roth, A., Gaska, J., and Smith, D. (2018). 2018 Wisconsin Soybean Performance Trials (No. A3654). University of Wisconsin-Madison.
Cooper, B., and Campbell, K. B. (2017). Protection against common bean rust conferred by a gene-silencing method. Phytopathology 107, 920–927. doi: 10.1094/PHYTO-03-17-0095-R
Dang, Y., Yang, Q., Xue, Z., and Liu, Y. (2011). RNA interference in fungi: pathways, functions, and applications. Eukaryotic Cell. 10, 1148–1155. doi: 10.1128/EC.05109-11
Dou, T., Shao, X., Hu, C., Liu, S., Sheng, O., Bi, F., et al. (2020). Host-induced gene silencing of Foc TR 4 ERG 6/11 genes exhibits superior resistance to Fusarium wilt of banana. Plant Biotechnol. J. 18, 11–13. doi: 10.1111/pbi.13204
Dubrovina, A. S., and Kiselev, K. V. (2019). Exogenous RNAs for gene regulation and plant resistance. Int. J. Mol. Sci. 20:2282. doi: 10.3390/ijms20092282
Dutton, M. V., and Evans, C. S. (1996). Oxalate production by fungi: its role in pathogenicity and ecology in the soil environment. Can. J. Microbiol. 42, 881–895. doi: 10.1139/m96-114
Elbashir, S. M., Martinez, J., Patkaniowska, A., Lendeckel, W., and Tuschl, T. (2001). Functional anatomy of siRNAs for mediating efficient RNAi in Drosophila melanogaster embryo lysate. EMBO J. 20, 6877–6888. doi: 10.1093/emboj/20.23.6877
Ghag, S. B., Shekhawat, U. K., and Ganapathi, T. R. (2014). Host-induced post-transcriptional hairpin RNA-mediated gene silencing of vital fungal genes confers efficient resistance against Fusarium wilt in banana. Plant Biotechnol. J. 12, 541–553. doi: 10.1111/pbi.12158
Gilbert, M. K., Majumdar, R., Rajasekaran, K., Chen, Z. Y., Wei, Q., Sickler, C. M., et al. (2018). RNA interference-based silencing of the alpha-amylase (amy1) gene in Aspergillus flavus decreases fungal growth and aflatoxin production in maize kernels. Planta 247, 1465–1473. doi: 10.1007/s00425-018-2875-0
Gogoi, A., Sarmah, N., Kaldis, A., Perdikis, D., and Voloudakis, A. (2017). Plant insects and mites uptake double-stranded RNA upon its exogenous application on tomato leaves. Planta 246, 1233–1241. doi: 10.1007/s00425-017-2776-7
Hoffman, D. D., Diers, B. W., Hartman, G. L., Nickell, C. D., Nelson, R. L., Pedersen, W. L., et al. (2002). Selected soybean plant introductions with partial resistance to Sclerotinia sclerotiorum. Plant Dis. 86, 971–980. doi: 10.1094/PDIS.2002.86.9.971
Hu, D., Chen, Z. Y., Zhang, C., and Ganiger, M. (2020). Reduction of Phakopsora pachyrhizi infection on soybean through host-and spray-induced gene silencing. Mol. Plant Pathol. 21, 794–807. doi: 10.1111/mpp.12931
Huzar-Novakowiski, J., and Dorrance, A. E. (2018). Ascospore inoculum density and characterization of components of partial resistance to Sclerotinia sclerotiorum in soybean. Plant Dis. 102, 1326–1333. doi: 10.1094/PDIS-11-17-1786-RE
Iquira, E., Humira, S., and François, B. (2015). Association mapping of QTLs for sclerotinia stem rot resistance in a collection of soybean plant introductions using a genotyping by sequencing (GBS) approach. BMC Plant Biol. 15:5. doi: 10.1186/s12870-014-0408-y
Kabbage, M., Williams, B., and Dickman, M. B. (2013). Cell death control: the interplay of apoptosis and autophagy in the pathogenicity of Sclerotinia sclerotiorum. PLoS Pathog. 9:e1003287. doi: 10.1371/journal.ppat.1003287
Kettles, G. J., Hofinger, B. J., Hu, P., Bayon, C., Rudd, J. J., Balmer, D. A., et al. (2019). SRNA profiling combined with gene function analysis reveals a lack of evidence for cross-kingdom RNAi in the wheat-Zymoseptoria tritici pathosystem. Front. Plant Sci. 10:892. doi: 10.3389/fpls.2019.00892
Kim, H. S., and Diers, B. W. (2000). Inheritance of partial resistance to sclerotinia stem rot in soybean. Crop Sci. 40, 55–61. doi: 10.2135/cropsci2000.40155x
Kim, K. S., Min, J. Y., and Dickman, M. B. (2008). Oxalic acid is an elicitor of plant programmed cell death during Sclerotinia sclerotiorum disease development. Mol. Plant-Microbe Interact. 21, 605–612. doi: 10.1094/MPMI-21-5-0605
Koch, A., Biedenkopf, D., Furch, A., Weber, L., Rossbach, O., Abdellatef, E., et al. (2016). An RNAi-based control of Fusarium graminearum infections through spraying of long dsRNAs involves a plant passage and is controlled by the fungal silencing machinery. PLoS Pathog. 12:e1005901. doi: 10.1371/journal.ppat.1005901
Koch, A., Höfle, L., Werner, B. T., Imani, J., Schmidt, A., Jelonek, L., et al. (2019). SIGS vs HIGS: a study on the efficacy of two dsRNA delivery strategies to silence Fusarium FgCYP51 genes in infected host and non-host plants. Mol. Plant Pathol. 20, 1636–1644. doi: 10.1111/mpp.12866
Li, H., Guan, R., Guo, H., and Miao, X. (2015). New insights into an RNAi approach for plant defense against piercing-sucking and stem-borer insect pests. Plant Cell Environ. 38, 2277–2285. doi: 10.1111/pce.12546
Liang, X., Liberti, D., Li, M., Kim, Y.-T., Hutchens, A., Wilson, R., et al. (2015). Oxaloacetate acetylhydrolase gene mutants of Sclerotinia sclerotiorum do not accumulate oxalic acid but do produce limited lesions on host plants. Mol. Plant Pathol. 16, 559–571. doi: 10.1111/mpp.12211
Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real- time quantitative PCR and the 2– ΔΔCT method. Methods 25, 402–408. doi: 10.1006/meth.2001.1262
Mahto, B. K., Singh, A., Pareek, M., Rajam, M. V., Dhar-Ray, S., and Reddy, P. M. (2020). Host-induced silencing of the Colletotrichum gloeosporioides conidial morphology 1 gene (CgCOM1) confers resistance against Anthracnose disease in chili and tomato. Plant Mol. Biol. 104, 381–395. doi: 10.1007/s11103-020-01046-3
McCaghey, M., Willbur, J., Ranjan, A., Grau, C. R., Chapman, S., Diers, B., et al. (2017). Development and evaluation of Glycine max germplasm lines with quantitative resistance to Sclerotinia sclerotiorum. Front. Plant Sci. 8:1495. doi: 10.3389/fpls.2017.01495
McLoughlin, A. G., Wytinck, N., Walker, P. L., Girard, I. J., Rashid, K. Y., de Kievit, T., et al. (2018). Identification and application of exogenous dsRNA confers plant protection against Sclerotinia sclerotiorum and Botrytis cinerea. Sci. Rep. 8, 1–14. doi: 10.1038/s41598-018-25434-4
Mitter, N., Worrall, E. A., Robinson, K. E., Li, P., Jain, R. G., Taochy, C., et al. (2017). Clay nanosheets for topical delivery of RNAi for sustained protection against plant viruses. Nat. Plants 3, 1–10. doi: 10.1038/nplants.2016.207
Nowara, D., Gay, A., Lacomme, C., Shaw, J., Ridout, C., Douchkov, D., et al. (2010). HIGS: Host-induced gene silencing in the obligate biotrophic fungal pathogen, Blumeria graminis. Plant Cell 22, 3130–3141. doi: 10.1105/tpc.110.077040
Panwar, V., Jordan, M., McCallum, B., and Bakkeren, G. (2018). Host?induced silencing of essential genes in Puccinia triticina through transgenic expression of RNA i sequences reduces severity of leaf rust infection in wheat. Plant Biotechnol. J. 16, 1013–1023. doi: 10.1111/pbi.12845
Panwar, V., McCallum, B., and Bakkeren, G. (2013). Host-induced gene silencing of wheat leaf rust fungus Puccinia triticina pathogenicity genes mediated by the Barley stripe mosaic virus. Plant Mol. Biol. 81, 595–608. doi: 10.1007/s11103-013-0022-7
Peltier, A. J., Bradley, C. A., Chilvers, M. I., Malvick, D. K., Mueller, D. S., Wise, K. A., et al. (2012). Biology, yield loss and control of Sclerotinia stem rot of soybean. J. Integr. Pest Manag. 3, B1–B7. doi: 10.1603/IPM11033
Peltier, A. J., and Grau, C. R. (2008). The influence of light on relationships between sclerotinia stem rot of soybean in field and controlled environments. Plant Dis. 92, 1510–1514. doi: 10.1094/PDIS-92-11-1510
Pliego, C., Nowara, D., Bonciani, G., Gheorghe, D. M., Xu, R., Surana, P., et al. (2013). Host-induced gene silencing in barley powdery mildew reveals a class of ribonuclease-like effectors. Mol. Plant Microbe Interact. 26, 633–642. doi: 10.1094/MPMI-01-13-0005-R
Ranjan, A., Jayaraman, D., Grau, C., Hill, J. H., Whitham, S. A., Ané, J.-M., et al. (2018). The pathogenic development of Sclerotinia sclerotiorum in soybean requires specific host NADPH oxidases. Mol. Plant Pathol. 19, 700–714. doi: 10.1111/mpp.12555
Roberts, J., Pitkin, J., Adams, T., and Monsanto Technology, L. L. C. (2008). In Planta RNAi Control of Fungi. U.S. Patent Application 11/670409.
Rollins, J. A. (2003). The Sclerotinia sclerotiorum pac1 gene is required for sclerotial development and virulence. Mol. Plant Microbe Interact. 16, 785–795. doi: 10.1094/MPMI.2003.16.9.785
Shaner, G., and Finney, R. E. (1977). The effect of nitrogen fertilization on the expression of slow-mildewing resistance in Knox wheat. Phytopathology 67, 1051–1056. doi: 10.1094/Phyto-67-1051
Shanmugam, V., Sharma, V., Bharti, P., Jyoti, P., Yadav, S. K., Aggarwal, R., et al. (2017). RNAi induced silencing of pathogenicity genes of Fusarium spp. for vascular wilt management in tomato. Ann. Microbiol. 67, 359–369. doi: 10.1007/s13213-017-1265-3
Song, X. S., Gu, K. X., Duan, X. X., Xiao, X. M., Hou, Y. P., Duan, Y. B., et al. (2018). Secondary amplification of siRNA machinery limits the application of spray-induced gene silencing. Mol. Plant Pathol. 19, 2543–2560. doi: 10.1111/mpp.12728
Tinoco, M. L. P., Dias, B. B., Dall'Astta, R. C., Pamphile, J. A., and Aragão, F. J. (2010). In vivo trans-specific gene silencing in fungal cells by in planta expression of a double-stranded RNA. BMC Biol. 8:27. doi: 10.1186/1741-7007-8-27
Tiwari, I. M., Jesuraj, A., Kamboj, R., Devanna, B. N., Botella, J. R., and Sharma, T. R. (2017). Host delivered RNAi, an efficient approach to increase rice resistance to sheath blight pathogen (Rhizoctonia solani). Sci. Rep. 7, 1–14. doi: 10.1038/s41598-017-07749-w
Wagner, T. A., Cai, Y., Bell, A. A., Puckhaber, L. S., Magill, C., Duke, S. E., et al. (2020). RNAi suppression of CYP82D P450 hydroxylase, an enzyme involved in gossypol biosynthesis, enhances resistance to Fusarium wilt in cotton. J. Phytopathol. 168, 103–112. doi: 10.1111/jph.12873
Wang, M., Weiberg, A., Lin, F.-M., Thomma, B. P. H. J., Huang, H.-D., and Jin, H. (2016). Bidirectional cross-kingdom RNAi and fungal uptake of external RNAs confer plant protection. Nat. Plants 2, 1–10. doi: 10.1038/nplants.2016.151
Westrick, N. M., Ranjan, A., Jain, S., Grau, C. R., Smith, D. L., and Kabbage, M. (2019). Gene regulation of Sclerotinia Sclerotiorum during infection of glycine max: on the road to pathogenesis. BMC Genomics 20:157. doi: 10.1186/s12864-019-5517-4
Whitham, S. A., Lincoln, L. M., Chowda-Reddy, R. V., Dittman, J. D., O'Rourke, J. A., and Graham, M. A. (2016). Virus-induced gene silencing and transient gene expression in soybean (Glycine max) using Bean pod mottle virus infectious clones. Curr. Protoc. Plant Biol. 1, 263–283. doi: 10.1002/cppb.20012
Willbur, J., McCaghey, M., Kabbage, M., and Smith, D. L. (2019). An overview of the Sclerotinia sclerotiorum pathosystem in soybean: impact, fungal biology, and current management strategies. Trop. Plant Pathol. 44, 3–11. doi: 10.1007/s40858-018-0250-0
Willbur, J. F., Ding, S., Marks, M. E., Lucas, H., Grau, C. R., Groves, C. L., et al. (2017). Comprehensive Sclerotinia stem rot screening of soybean germplasm requires multiple isolates of Sclerotinia sclerotiorum. Plant Dis. 101, 344–353. doi: 10.1094/PDIS-07-16-1055-RE
Williams, B., Kabbage, M., Kim, H. J., Britt, R., and Dickman, M. B. (2011). Tipping the balance: Sclerotinia sclerotiorum secreted oxalic acid suppresses host defenses by manipulating the host redox environment. PLoS Pathog. 7:e1002107. doi: 10.1371/journal.ppat.1002107
Workneh, F., and Yang, X. B. (2000). Prevalence of Sclerotinia stem rot of soybeans in the north- central United States in relation to tillage, climate, and latitudinal positions. Phytopathology 90, 1375–1382. doi: 10.1094/PHYTO.2000.90.12.1375
Wrather, A., Shannon, G., Balardin, R., Carregal, L., Escobar, R., Gupta, G. K., et al. (2010). Effect of diseases on soybean yield in the top eight producing countries in 2006. Plant Health Progress 11:5. doi: 10.1094/PHP-2010-0102-01-RS
Xu, L., Xiang, M., White, D., and Chen, W. (2015). pH dependency of sclerotial development and pathogenicity revealed by using genetically defined oxalate-minus mutants of Sclerotinia sclerotiorum. Environ. Microbiol. 17, 2896–2909. doi: 10.1111/1462-2920.12818
Yin, C., Downey, S. I., Klages-Mundt, N. L., Ramachandran, S., Chen, X., Szabo, L. J., et al. (2015). Identification of promising host-induced silencing targets among genes preferentially transcribed in haustoria of Puccinia. BMC Genomics 16:579. doi: 10.1186/s12864-015-1791-y
Yin, C., Jurgenson, J. E., and Hulbert, S. H. (2010). Development of a host-induced RNAi system in the wheat stripe rust fungus Puccinia striiformis f. sp. tritici. Mol. Plant. Microbe Interact. 24, 554–561. doi: 10.1094/MPMI-10-10-0229
Yin, C., Park, J. J., Gang, D. R., and Hulbert, S. H. (2014). Characterization of a tryptophan 2 monooxygenase gene from Puccinia graminis f. sp. tritici involved in auxin biosynthesis and rust pathogenicity. Mol. Plant. Microbe Interact. 27, 227–235. doi: 10.1094/MPMI-09-13-0289-FI
Keywords: host induced gene silencing, RNAi, Sclerotinia sclerotiorum, soybean, oxalic acid, Ssoah1, bean pod mottle virus, VIGS
Citation: McCaghey M, Shao D, Kurcezewski J, Lindstrom A, Ranjan A, Whitham SA, Conley SP, Williams B, Smith DL and Kabbage M (2021) Host-Induced Gene Silencing of a Sclerotinia sclerotiorum oxaloacetate acetylhydrolase Using Bean Pod Mottle Virus as a Vehicle Reduces Disease on Soybean. Front. Plant Sci. 12:677631. doi: 10.3389/fpls.2021.677631
Received: 08 March 2021; Accepted: 26 May 2021;
Published: 20 July 2021.
Edited by:
Maria Del Mar Jimenez-Gasco, Pennsylvania State University (PSU), United StatesReviewed by:
Bo Li, Huazhong Agricultural University, ChinaJia-He Wu, Institute of Microbiology, Chinese Academy of Sciences, China
Copyright © 2021 McCaghey, Shao, Kurcezewski, Lindstrom, Ranjan, Whitham, Conley, Williams, Smith and Kabbage. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Mehdi Kabbage, kabbage@wisc.edu; Damon L. Smith, damon.smith@wisc.edu
†These authors have contributed equally to this work