Yingru Xu1†Xiaotong Ji1†Zhuangzhuang Xu1Yanping Yuan1Xiling Chen1Derong Kong1
Yanlong Zhang1,2
Daoyang Sun1,2*
- 1College of Landscape Architecture and Arts, Northwest A&F University, Yangling, China
- 2National Engineering Technology Research Center for Oil Peony, Northwest A&F University, Yangling, China
RNA silencing is a common antiviral mechanism in eukaryotic organisms. However, the transcriptional regulatory mechanism that controls the RNA silencing process remains elusive. Here, we performed high-depth transcriptome analysis on petunia (Petunia hybrida) leaves infected with tobacco rattle virus (TRV) strain PPK20. A total of 7,402 differentially expressed genes (DEGs) were identified. Of them, some RNA silencing-related transcripts, such as RNA-dependent RNA polymerases (RDRs), Dicer-like RNase III enzymes (DCLs), and Argonautes (AGOs), were induced by viral attack. Furthermore, we performed TRV-based virus-induced gene silencing (VIGS) assay on 39 DEGs encoding putative transcription factors (TFs), using green fluorescent protein (GFP) and phytoene desaturase (PhPDS) as reporters. Results showed that the down-regulation of PhbHLH41, PhbHLH93, PhZPT4-3, PhCOL4, PhHSF-B3A, PhNAC90, and PhWRKY75 led to enhanced TRV accumulation and inhibited PhPDS-silenced photobleaching phenotype. In contrast, silencing of PhERF22 repressed virus accumulation and promoted photobleaching development. Thus, these TFs were identified as potential positive and negative regulators of antiviral RNA silencing, respectively. One positive regulator PhCOL4, belonging to the B-box zinc finger family, was selected for further functional characterization. Silencing and transient overexpression of PhCOL4 resulted in decreased and increased expression of several RNA silencing-related genes. DNA affinity purification sequencing analysis revealed that PhCOL4 targeted PhRDR6 and PhAGO4. Dual luciferase and yeast one-hybrid assays determined the binding of PhCOL4 to the PhRDR6 and PhAGO4 promoters. Our findings suggest that TRV-GFP-PhPDS-based VIGS could be helpful to identify transcriptional regulators of antiviral RNA silencing.
Introduction
RNA silencing is a highly conserved and sequence-specific RNA degradation mechanism in plants. Under biotic stresses such as virus infection, RNA silencing can function as a natural defense mechanism against foreign nucleic acids introduced by viruses (Wang and Metzlaff, 2005). RNA silencing is triggered by the specific RNAs with double-stranded features (Llave, 2010). The whole process can be divided into three phases: initiation, effector, and amplification (Csorba and Burgyán, 2016). RNA silencing initiates when double-stranded RNA (dsRNA) is recognized by Dicer-like RNase III enzyme (DCL) and is cleaved into small pieces of 21–24-nt small interfering RNAs (siRNAs). The 3′ terminal of siRNA is 2′-O-methylated by a methyltransferase HEN1 to maintain its structural stability (Ruiz-Ferrer and Voinnet, 2009). Then, the mature siRNA is incorporated into the RNA-induced silencing complex (RISC), which contains the core argonaute (AGO) protein for targeting and down-regulating homologous mRNA (Liu and Chen, 2016). Next, RNA-dependent RNA polymerase (RDR) uses primary siRNA as a primer to synthesize dsRNA, which is re-processed into secondary siRNAs by DCL enzyme to amplify the silencing signals (Pumplin and Voinnet, 2013). At the end of RNA silencing process, transcription of viral RNAs can be substantially reduced at the molecular level, thereby enhancing host resistance to viral attack.
A set of crucial proteins such as RDRs, DCLs, and AGOs have been known to be implicated in antiviral RNA silencing pathway in plants (Garcia-Ruiz et al., 2010). In the model plant Arabidopsis, a total of 6 RDRs, 4 DCLs, and 10 AGOs have been identified and characterized (Eamens et al., 2008). RDR1, RDR2, and RDR6 can resist the infections of various RNA viruses through the generation of secondary siRNAs (Gaffar and Koch, 2019). The RDR6 mutant displays increased susceptibility to cucumber mosaic virus (CMV) and potato virus X (PVX) in Nicotiana tabacum (Li et al., 2012). RDR1 mutation results in compromised resistance of tobacco plants to tobacco mosaic virus (TMV) (Vaucheret, 2006). Simultaneous suppression of RDR1, RDR2, and RDR6 causes a decreased resistance to turnip mosaic virus (TuMV) in Arabidopsis (Garcia-Ruiz et al., 2010).
In addition, DCL1 cleaves the imperfect hairpin RNA (hpRNA) into 21-nt miRNA (Henderson et al., 2006), while DCL2, DCL3, and DCL4 process long and nearly perfect dsRNA into 22-nt, 24-nt, and 21-nt siRNAs, respectively (Bologna and Voinnet, 2014). DCL1 has a minor contribution to the cleavage of dsRNA since only a small amount of virus-derived siRNAs are detected in PVX-infected Arabidopsis dcl2/dcl3/dcl4 triple mutant (Ruiz-Ferrer and Voinnet, 2009; Jaubert et al., 2011). Instead, DCL1 mainly functions by negatively regulating the transcription of DCL3 and DCL4 (Qu et al., 2008). DCL2 and DCL4 exhibit major, redundant, but hierarchical functions in antiviral defense. DCL2 carves out an important role when DCL4 is absent or suppressed. Reduced transcription of DCL2 and DCL4 leads to hyper-susceptibility to virus infection (Deleris et al., 2006; Pumplin and Voinnet, 2013). In contrast, DCL3 seems to play a minor role in RNA silencing-induced antiviral defense but a vital role in RNA-directed DNA methylation (Singh et al., 2019).
Most AGO proteins in Arabidopsis possess antiviral activities. Among them, AGO1 and AGO2 are the primary effectors against RNA viruses, while AGO4 functions predominately in defending against DNA viruses (Carbonell and Carrington, 2015; Kenesi et al., 2021). AGO1 has been revealed to be involved in antiviral response to turnip crinkle virus (Azevedo et al., 2010) and CMV (Zhang et al., 2006). AGO2 plays an essential role in resisting CMV (Harvey et al., 2011), PVX (Jaubert et al., 2011), and TuMV (Manacorda et al., 2021). AGO1 and AGO7 function cooperatively to ensure the removal of viral RNAs with different secondary structures (Qu et al., 2008). Furthermore, AGO proteins exhibit variable binding properties to siRNAs arising from cleaved dsRNA. AGO4 binds to DCL3-dependent 24-nt siRNA and participates in chromatin modification (Gao et al., 2010). AGO1 interacts with 21-nt and 22-nt small RNAs through different mechanisms. AGO4, AGO6, and AGO9 bind to 24-nt siRNA to act specifically in the DNA methylation pathway (Havecker et al., 2010; Zhu et al., 2011).
Apart from the key enzymes mentioned above, antiviral RNA silencing in plants also requires the function of transcription factors (TFs) which regulate the spatiotemporal expression of the above enzyme-encoding genes. Various TF families, including AP2/ERF (Fischer and Dröge-Laser, 2004), ARF (Zhang et al., 2020), bHLH (Aparicio and Pallás, 2017), bZIP (Herath and Verchot, 2021), HD-ZIP (Zou et al., 2016), HSF (Merkling et al., 2015), MADS (Yang et al., 2014), MYB (Li et al., 2018), NAC (Sun et al., 2019), WRKY (Freeborough et al., 2021), and zinc finger (Sera, 2017), have been demonstrated to participate in the defense response to virus infection. A broad range of biological pathways associated with innate immunity or defensive hormone signaling are potentially controlled by these TFs. However, few studies on TFs involved in RNA silencing-elicited antiviral defense have been reported.
Virus-induced gene silencing (VIGS) is a high-throughput reverse genetic technology, which was developed based on antiviral RNA silencing mechanism (Unver and Budak, 2009). The VIGS system requires the construction of an artificial viral vector, such as tobacco rattle virus (TRV), carrying endogenous gene fragments. The inoculation of viral constructs can cause the degradation of target genes through RNA silencing machinery in plants. In previous studies, we identified an ethylene-responsive element binding factor PhERF2 and an ocs element binding factor PhOBF1, which played critical roles in antiviral RNA silencing in petunia (Sun et al., 2016, 2017). Inclusion of PhERF2 or PhOBF1 fragment in the TRV vector containing reporter genes PhPDS and PhCHS led to impaired leaf photobleaching and white-corollas phenotypes. PhERF2 and PhOBF1 functioned as crucial participants in the defense against TRV, CMV, or TMV infection, and affected the transcription of several RDR, DCL, and AGO genes in the RNA silencing pathway. These results suggested that PhPDS or PhCHS can be used as an effective reporter gene for VIGS-based RNA silencing studies. A similar approach has been reported in Arabidopsis, in which a TRV-AtPDS construct was used to assess the antiviral function of AGO1-27 and DCL2/DCL3/DCL4 (Ma et al., 2015). These studies prompt us to identify novel regulators of antiviral RNA silencing using VIGS method. Here, we combined high-depth Illumina RNA sequencing (RNA-Seq) and TRV-GFP-PhPDS-based VIGS assays to screen key TFs potentially involved in antiviral RNA silencing in petunia. A number of candidate TFs that can modulate, either positively or negatively, viral resistance through RNA silencing were obtained. In particular, a zinc finger TF, annotated as PhCOL4, was validated for its essential role in antiviral RNA silencing.
Materials and Methods
Plant Materials and Growth Conditions
The seeds of petunia (Petunia hybrida) cultivar ‘Primetime Blue’ were purchased from Goldsmith Seeds Inc. (Gilroy, CA, United States). They were planted in a 72-well plastic tray filled with the soil mix of peat moss, perlite, and vermiculite at 2:1:1 by volume. The tray was placed into a growth chamber at 25/22°C temperature (day/night) with a photoperiod of 16/8 h (light/dark) and a 70% relative humidity. Young seedlings were transplanted to small plastic pots containing the same soil mix for inoculation with wild TRV strain and Agrobacterium culture bearing TRV-GFP-PhPDS or pCAMBIA2300 constructs. Wild TRV strain-inoculated petunia leaves at three stages were used for Illumina RNA-Seq. The leaves infected with various TRV constructs were used for VIGS assay. The leaves infiltrated with pCAMBIA2300 constructs were used for transient overexpression assay. To examine the abundances of transcripts, the inoculated or uppermost systemically infected leaves were collected for RNA extraction.
Inoculation of the Wild Type Tobacco Rattle Virus PPK20
Tobacco rattle virus strain PPK20 was preserved in tobacco (N. tabacum) cultivar ‘NC82’ plantlets. For infectious sap preparation, the young leaves infected with TRV (PPK20) were harvested and homogenized with a small amount of kieselguhr in 100 mM phosphate buffer at pH 7.0 (1:6, w/v) using a mortar and pestle. The viral preparations were mixed with 10% of 600-mesh carborundum powder by weight, and then mechanically inoculated onto the fully expanded healthy young leaves (Tian and Valkonen, 2013). To ensure sufficient viral propagation, the rub-inoculation process was performed once more within half an hour. The inoculum was flushed with sterile deionized water. The symptom occurrence was observed during the post-inoculation periods. Relative accumulation levels of TRV (RNA1 and RNA2) were determined by real-time quantitative PCR (RT-qPCR) analysis.
RNA Preparation and Illumina Sequencing
Petunia leaves infected with TRV (PPK20) were collected to extract total RNA using TRIzol reagent (Invitrogen, Carlsbad, CA, United States). Three independent biological replicates were used for each stage. RNA samples were further purified to eliminate contaminant DNA using RNase-free DNase I (Promega, Madison, WI, United States). RNA integrity was examined by observing the major RNA species (18S and 28S rRNA) through 1.5% agarose gel electrophoresis. RNA yield was quantified using a NanoDrop ND-2000c Spectrophotometer (NanoDrop Technologies, Wilmington, DE, United States). Illumina sequencing and generation of transcriptome data were conducted at Gene Denovo Bio-Technology Co., LTD. (Guangzhou, China). In brief, the general experimental procedure includes mRNA enrichment and fragmentation, cDNA synthesis, adaptor ligation, and ligated products sequencing as previously described (Sun et al., 2019).
De novo Assembly of Transcripts
High-quality clean reads were produced by removing adaptor, low-quality, and unknown nucleotides of raw reads through a Perl script tool. The resulting clean reads were assembled via Trinity program (v2.1.1) (Grabherr et al., 2011). The redundancy of assembled products was removed to obtain unique fragments using iAssembler software (Zheng et al., 2011). The fragments over 200 bp were enriched using TIGR Gene Indices clustering tools (TGICL, v2.1) (Pertea et al., 2003), and regarded as unigenes. The unigenes were mapped to the reference genome1 of P. axillaris, which is one of P. hybrida’s parents. Novel transcripts were determined using Cuffcompare software through a threshold for length ≥ 200 bp and exon number ≥ 2 (Trapnell et al., 2012), and annotated against the public protein databases including NR, Swiss-Prot, COG, and KEGG.
Differentially Expressed Gene Analysis
Read counts of unigenes were calculated and standardized to the fragments per kilobase of transcript per million mapped reads (FPKM), which were used to quantify the expression levels of unigenes through Cufflinks program (v2.1.1) (Trapnell et al., 2010). Differentially expressed genes (DEGs) were identified using edgeR package (v3.12) (Chen et al., 2014). Statistical significance of DEGs was evaluated with a cutoff of fold change (FC) ≥ 2.0 and false discovery rate (FDR) ≤ 0.05. Blast2GO tool (v2.3.5) and WEGO program2 were used to gain knowledge about Gene Ontology (GO) classification of DEGs (Conesa et al., 2005; Ye et al., 2006). To assign putative metabolic processes to DEGs, the pathway enrichment analysis was carried out using a BLAST search against the Kyoto Encyclopedia of Genes and Genomes (KEGG) database3 (Wixon and Kell, 2000). P-values were calculated using a previously reported method and subjected to FDR correction (Lu et al., 2012). Significantly enriched pathways were characterized using FDR ≤ 0.05 as a cutoff.
Measurement of Endogenous Hormones
Petunia leaves inoculated with TRV (PPK20) at different time points were harvested. For ethylene (ET) detection, the leaves were placed in a 50-ml plastic tube at 25°C for 4 h, with an airproof rubber plug embedded into the cap. The parafilm was used to seal the connection gaps. A 3-ml gas sample was withdrawn from the tube using a gas-tight syringe with a needle, and applied to a gas chromatograph (GC-8A; Shimadzu, Kyoto, Japan) for measurement of ET release. The examination of other hormones, including abscisic acid (ABA), jasmonic acid (JA), and salicylic acid (SA), was performed as previously described (Mornya and Cheng, 2011) with some minor modification. Briefly, petunia samples were extracted in 80% methanol and 1 mM butylated hydroxytoluence on a rotary shaker at 4°C for 1 h. After a complete homogenization, the mixture was centrifuged at 12,000 rpm for 20 min at 4°C. The supernatant was collected and purified through a 0.45 μm filter. It was concentrated to near dryness in N2 gas and then dissolved in 50% methanol. The resulting solution was analyzed by a high performance liquid chromatography-electrospray ionization tandem mass spectrometry (HPLC-ESI-MS/MS). All chemical substances purchased from Sigma-Aldrich (St Louis, MO, United States) were used as the reference standards.
Real-Time Quantitative PCR
Total RNA from petunia leaves was extracted as mentioned above. After purification and quantification, about 2–5 μg of RNA samples were used to synthesize first-strand cDNA with SuperScript III reverse transcriptase (Invitrogen, Carlsbad, CA, United States) according to the manufacturer’s recommendations. The resulting cDNA products were fivefold diluted by ddH2O and restored at –20°C until further use. The specific primer pairs were designed based on the known cDNA sequences of genes tested via Primer Premier54 (Supplementary Table S1). 26S ribosomal RNA served as a reference gene. RT-qPCR assay was performed using the SYBR Green PCR Master Mix (2×) (Applied Biosystems, Foster City, CA, United Stataes) in a LightCycler480 Real-Time PCR System (Roche Diagnostic, Basel, Switzerland). The experiment was repeated three times using independent RNA samples. Relative expression levels of genes were analyzed using a previously described protocol (Livak and Schmittgen, 2001).
Virus-Induced Gene Silencing Assay
The TRV-GFP construct was generated as previously described (Tian et al., 2014). A 138-bp fragment of PhPDS cDNA was amplified using the forward and reverse primers with a SacI and an XhoI restriction site, respectively. It was ligated into the corresponding sites of TRV-GFP plasmid to form the TRV-GFP-PhPDS construct. Thirty-nine selected candidate TF genes were cloned into the TRV-GFP-PhPDS construct between BamHI and KpnI sites using the Seamless Cloning and Assembly Kit (Novoprotein, Shanghai, China). Oligonucleotide primers used for gene fragment amplification and expression analysis are listed in Supplementary Table S1. The resulting TRV plasmids were electrotransformed into A. tumefaciens strain GV3101. The agro-inoculation was performed according to a previously described protocol (Sun et al., 2017). Three-week-old seedlings of petunia were used for the agro-inoculation. For expression assessment of PhPDS and TF genes, the sequences beyond the inserted fragments were used for primer design (Supplementary Table S1). For GFP fluorescence detection, fluorescent foci were visualized using a portable blue light-emitting diode flashlight as an excitation source at 450 nm (LUYOR-3260RB, LUYOR, Shanghai, China). The images were taken using a digital camera (EOS 20D, Canon, Tokyo, Japan) with a bypass emission filter (LUV-495, LUYOR, Shanghai, China). The intensity of fluorescence was quantified using ImageJ software (NIH, Bethesda, MD, United States).
Transient Overexpression of PhCOL4
Agrobacterium-mediated transient assay was performed as previously described (Mußmann et al., 2011) with some minor changes. The complete coding region of PhCOL4 (Supplementary Figure S1) was cloned into pCAMBIA2300 vector between EcoRI and SacI sites. The resulting plasmids were transformed into A. tumefaciens strain LBA4404. The transformed bacteria were cultured in liquid LB media supplemented with 40 mg⋅L–1 kanamycin for overnight at 28°C. The Agrobacterium cells were harvested and resuspended to an OD600 of 4.0 in the infiltration buffer containing 10 mM MgCl2, 10 mM MES, and 200 μM acetosyringone. The suspension was infiltrated into fully developed petunia leaves through a syringe. At 6 days post infiltration, the leaves were collected for gene expression analysis.
DNA Affinity Purification Sequencing Assay
DAP-Seq binding experiment was carried out according to a previously described method (Yao et al., 2020). In brief, genomic DNA (gDNA) was extracted from petunia leaves using CTAB method, and fragmented through a Covaris M220 (Woburn, MA, United States). Purification of fragmented gDNA was performed using AMPure XP beads (Beckman Coulter, Indianapolis, IN, United States). DNA library was prepared using the NEXTFLEX Rapid DNA-Seq Kit (PerkinElmer, Austin, TX, United States). The coding sequence of PhCOL4 was introduced into pFN19K HaloTag T7 SP6 Flexi expression vector. The TNT SP6 Coupled Wheat Germ Extract System (Promega, Madison, WI, United States) was used to express the Halo-PhCOL4 fusion proteins, which were then captured by Magne Halo Tag Beads (Promega, Madison, WI, United States). The proteins bound to beads were incubated with gDNA fragments at 25°C for 2 h. A high temperature of 98°C for 10 min was used to denature the proteins and release the bound gDNA fragments. The fragments were PCR-amplified, purified, and sequenced. For data processing, all reads were mapped to the reference genome of P. axillaris5. Peak analysis was conducted using Macs2 tool. Distribution of peaks in petunia genome sequence was analyzed using Homer (Heinz et al., 2010). Functional annotation of corresponding genes was determined through a blast against NR, NT, Swissprot, and Pfam databases. Protein-bound motifs were predicted using MEME-Chip suite 5.0.5 (Machanick and Bailey, 2011).
Dual Luciferase Assay
Dual luciferase assay was performed using a method as previously described (Sun et al., 2019). The coding sequence of PhCOL4 was ligated into pGreenII 62-SK vector, which was regarded as the effector. To construct the reporters, the promoter sequences of PhRDR6 and PhAGO4 (Supplementary Figure S1) were amplified using the primers listed in Supplementary Table S1, and cloned into pGreenII 0800-LUC vector expressing firefly luciferase (LUC) and Renilla luciferase (REN). The resulting plasmids were transformed into A. tumefaciens strain GV3101. Co-infiltration with effector and reporter constructs was carried out on petunia young leaves. LUC and REN activities were measured using a Tecan Infinite M200 luminometer (Männedorf, Switzerland). The data were expressed as a ratio of LUC to REN.
Yeast One-Hybrid Assay
The 35-bp fragments containing predicted PhCOL4-bound motifs in the PhRDR6 and PhAGO4 promoters were amplified using the primers (Supplementary Table S1). The fragments were ligated into pHIS2 vector to generate the baits. The full-length open reading frame region of PhCOL4 was fused to the GAL4 activation domain in pGADT7-Rec vector as the pray. All constructs were co-transformed into yeast strain Y187, which was grown and spotted on SD-Ura-His-Leu medium supplemented with 100 mM 3-aminotriazole (3-AT) at different dilutions.
Results
Accumulation of Tobacco Rattle Virus PPK20 in Petunia
To illustrate transcriptional regulatory mechanism of antiviral RNA silencing in petunia, we used an integrative strategy of TRV (PPK20)-infected Illumina RNA-seq and TRV-GFP-PhPDS-based VIGS to identify the TFs regulating the RNA silencing machinery. An important candidate TF was selected for further functional characterization. VIGS and transient overexpression assays were performed to analyze the impact of the identified TF on the expression of RNA silencing components. The target genes downstream of the TF were determined through DAP-Seq, dual luciferase, and yeast one-hybrid assays (Figure 1). Young leaves of P. hybrida at six-leaf-stage were inoculated with TRV (PPK20). A clear symptom development was observed in the inoculated leaves, which exhibited slight leaf yellowing at 3 days post inoculation (dpi) and severe leaf mottling or chlorosis at 6 dpi (Figure 2A). Two regions from TRV (PPK20) genomic RNA1 and RNA2 were chosen for quantification of virus RNA levels (Figure 2B). Virus accumulation analysis showed that virus RNA levels of both TRV-RNA1 and TRV-RNA2 increased gradually through 3 dpi and escalated sharply from 4 to 5 dpi. Then, the levels of TRV RNA transcripts continued to increase but at a much lower speed from 5 to 6 dpi (Figure 2C). Based on the trend of symptom development and viral growth, the infection process of TRV (PPK20) in petunia leaves was divided into three stages: 0 dpi (S0), 3 dpi (S3), and 6 dpi (S6), which were selected for transcriptome analysis.
Figure 1. Schematic representation of the experimental strategy to identify important transcription factors involved in antiviral RNA silencing in petunia.
Figure 2. Symptom development and transcriptome data evaluation in TRV PPk20-infected petunia leaves. (A) Representative phenotypes of petunia leaves at 0 (S0), 3 (S3), and 6 (S6) days post inoculation with TRV (PPK20). Scale bars = 5.0 mm. (B) Schematic diagram of the regions used for detection of TRV-RNA1 and TRV-RNA2 levels. (C) Relative accumulation levels of TRV-RNA1 and TRV-RNA2 in TRV (PPK20)-inoculated petunia leaves at intervals. Error bars represent standard error of the mean from three biological replicates. Transcript abundances were standardized to 26S rRNA. (D) Number of up-regulated and down-regulated unigenes in TRV (PPK20)-infected petunia leaves for the pairwise comparisons of S0-vs.-S3, S0-vs.-S6, and S3-vs.-S6. (E) Venn diagram of unique and shared differentially expressed unigenes for different pairwise comparisons.
Illumina RNA Sequencing and de novo Assembly
After data filtration, a total of 141,740,216, 118,866,564, and 117,466,686 clean RNA sequencing reads were obtained for petunia leaves at S0, S3, and S6, respectively. At each stage, 74.2% (105,214,806), 73.1% (86,903,750), and 73.3% (86,084,852) of these reads could be mapped to petunia reference genome. De novo assembly led to the generation of 34,486 unigenes, including 32,928 mapped and 1,558 novel genes (Table 1). The functions of 24,462 mapped unigenes were annotated based on protein databases (Table 1 and Supplementary Table S2). The sequences of all novel transcripts are shown in Supplementary Figure S2.
Table 1. Summary of TRV PPK20-infected petunia transcriptome sequencing dataset.
Detection of Differentially Expressed Genes During Tobacco Rattle Virus PPK20 Infection
Compared to S0, S3 displayed 3,413 up-regulated and 1,275 down-regulated unigenes. A large number of DEGs were detected for S6, which has 4,153 up-regulated and 1,782 down-regulated unigenes. When compared to S3, S6 only has a moderate number of DEGs: 895 up-regulated and 784 down-regulated (Figure 2D and Supplementary Tables S3–S5). The DEGs (7,402 in total) from pairwise comparison among S0, S3, and S6 were shown in Figure 2E. The highest number (3,449) of overlapped DEGs was observed between S0-S3 and S0-S6. Consistently, less shared DEGs (745) were found between S0-S3 and S3-S6, and 1,127 overlapped DEGs between S0–S6 and S3–S6. These data suggest that, compared to the control S0, the infected S3 and S6 have relatively closer transcriptome profiles with each other.
Gene Ontology (GO) enrichment analysis was carried out to investigate the functions of identified DEGs. At the molecular function level, the most significantly enriched GO terms correspond to “metabolic process,” “binding,” and “membrane” (Supplementary Figure S3A). At the biochemical pathway level, 1,581 DEGs for S0-S3, 2,007 DEGs for S0-S6, and 428 DEGs for S3-S6 were categorized into 126, 127, and 101 biochemical pathways, respectively. Among these pathways, “protein processing in endoplasmic reticulum,” “Plant–pathogen interaction,” and “Plant hormone signal transduction” were significantly enriched (Supplementary Figure S3B and Supplementary Table S6).
Induction of Endogenous Hormones and Their Biosynthetic Genes by Tobacco Rattle Virus PPK20 Infection
Plant hormones are known to play an important role in defense against viruses. “Plant hormone signal transduction” was detected as significantly enriched pathway in our identified DEGs after TRV (PPK20) infection. Therefore, we measured the levels of common plant hormones ET, ABA, JA, and SA during the TRV (PPK20)-infection process. Results showed that the accumulation of these hormones was significantly boosted by TRV (PPK20) infection. The highest levels of ET and JA in petunia leaves were found at 5 dpi, while ABA and SA contents peaked at 6 dpi (Figure 3A). These results are consistent with our transcriptome data, which showed that ET biosynthetic genes, including PhACO1, PhACO2, PhACO3, PhACO4, PhASC2, PhACS6, and PhACS7, were significantly up-regulated during TRV (PPK20) infection. The transcription of PhZEP, PhNCED3, PhNCED5, PhAAO4, PhSDR2a, PhSDR2b, and PhSDR3b rather than PhNCED7 in the ABA biosynthesis pathway was induced by TRV (PPK20) infection. In addition, most genes involved in JA production exhibited elevated expression except PhLOX3. Several SA biosynthetic genes, such as PhCM1, PhADTs (1 and 6), and PhPALs (1, 2a, and 2b), were also up-regulated upon TRV (PPK20) infection (Figure 3B).
Figure 3. Increased accumulation of several endogenous hormones and transcription of hormone biosynthetic genes by TRV PPK20 infection. (A) Accumulation levels of ethylene (ET), abscisic acid (ABA), jasmonic acid (JA), and salicylic acid (SA) in petunia leaves at various days post inoculation with TRV (PPK20). Error bars represent standard error of the mean from three biological replicates. (B) Expression levels of ET, ABA, JA, and SA biosynthesis-related genes in petunia leaves at three stages of TRV (PPK20) infection based on transcriptome data. The FPKMs of unigenes at the stage of S0-1 are normalized to 1.00, and transcript levels of all unigenes are transformed with Log2.
Identification of Antiviral RNA Silencing-Related Structural Genes and Transcriptional Regulators
Considering the essential roles of RDR, DCL, and AGO in antiviral RNA silencing, their transcripts differentially expressed during TRV (PPK20) infection were selected from the transcriptome data. Results showed that PhRDR1, PhRDR2, PhRDR6, PhDCL1a, PhDCL1b, PhDCL2, PhDCL3, PhDCL4, PhAGO1a, PhAGO1b, PhAGO2, PhAGO4, PhAGO5, PhAGO8, and PhAGO9 were dramatically up-regulated in TRV (PPK20)-infected petunia leaves. The highest transcriptional changes correspond to PhDCL1b whose expression rose by 130–260 times, followed by a 34–62 times increase in PhDCL2 transcription. In addition, the transcription of PhDCL1a and PhRDR6 increased moderately by less than three times (Figure 4).
Figure 4. Expression of transcripts involved in the RNA silencing pathway in TRV PPK20-infected petunia leaves. RNA-Seq and RT-qPCR analysis of transcript levels of RNA silencing-associated structural unigenes, including some RDRs, DCLs, and AGOs, in petunia leaves at different phases of TRV (PPK20) infection. 26S rRNA was used as a reference gene. Error bars represent standard error of the mean from three biological replicates. Statistical significance was determined using Student’s t-test (P < 0.05) and shown as asterisks.
To dissect the transcriptional regulatory mechanism for TRV defense in petunia, 469 DEGs encoding putative TFs were identified (Supplementary Tables S7–S9). Compared to S0, there are 286 (204 up-regulated and 82 down-regulated) TFs belonging to 32 families that were differentially expressed at S3. When S6 and S0 were compared, 374 (225 up-regulated and 149 down-regulated) TF-encoding (33 families) DEGs were identified. For S3 and S6, 148 (58 up-regulated and 90 down-regulated) TF-encoding DEGs from 24 families were found. The highest number of TFs included those from bHLH, bZIP, ZFP, ERF, MYB, NAC, and WRKY families (Figure 5A).
Figure 5. Distribution of differentially expressed transcription factors in petunia leaves upon TRV PPK20 infection. (A) Number of up-regulated (up) and down-regulated (down) transcription factors in TRV (PPK20)-inoculated petunia leaves for the pairwise comparisons of S0-vs.-S3, S0-vs.-S6, and S3-vs.-S6. (B) Schematic diagram of the TRV-GFP-PhPDS construct carrying the fragment of candidate transcription factor.
Screening of Antiviral RNA Silencing-Associated Transcription Factors Through Virus-Induced Gene Silencing
To identify critical transcriptional regulators involved in antiviral RNA silencing, 39 TRV-induced TFs were selected for TRV-GFP-PhPDS-based VIGS assay (Table 2). The fragments of these TFs were cloned into the multiple cloning sites of TRV-GFP-PhPDS vector to generate the recombinant constructs (Figure 5B). GFP and PhPDS were used as reporter genes for virus accumulation and gene silencing. Normal green fluorescence and leaf photobleaching were observed due to GFP expression and PhPDS silencing. Compared to TRV-GFP-PhPDS control, the inoculation with TRV-GFP-PhPDS constructs containing the fragments of PhARF11, PhbHLH41, PhbHLH93, PhZPT4-3, PhCOL4, PhGATA11, PhHSF-B3A, PhLOB1, PhNAC90, and PhWRKY75 resulted in enhanced fluorescent signals at 4 dpi. In contrast, the fluorescent brightness was suppressed in the leaves infected with TRV-GFP-PhPDS constructs containing PhPIF3, PhZAT10, PhERF22, PhERF1B, PhHSF-B3B, PhHSF-A2, and PhWRKY40 inserts (Figure 6). In addition, TRV-GFP-PhPDS/TGA2.1-, ZFP6-, and ERF22-infected systemic leaves showed promoted leaf photobleaching at 14 dpi, in comparison with the control. However, the plants inoculated with TRV-GFP-PhPDS vectors containing PhbHLH41, PhbHLH93, PhZPT4-3, PhZAT10, PhCOL4, PhBBX20, PhERF3, PhKD1-3, PhHSF-B3A, PhHSF-B3B, PhNAC90, PhJUB1B, and PhWRKY75 fragments displayed impaired photobleaching phenotypes (Figure 6).
Table 2. Selected candidate transcription factors used for screening antiviral RNA silencing-related genes.
Figure 6. Development of GFP fluorescence and photobleaching phenotype in petunia leaves with VIGS silencing of various transcription factors. Representative GFP-expressed fluorescence signals and PhPDS-silenced photobleaching phenotypes in petunia leaves inoculated with Agrobacterium containing 39 transcription factors (TFs)-inserted TRV-GFP-PhPDS constructs. Three-week-old seedlings of petunia were used for the agro-inoculation. Photographs were taken at 4 and 14 days post inoculation for fluorescence and photobleaching phenotypes. The TRV-GFP-PhPDS construct containing no TF insert was used as the control. In comparison with the control, enhanced and suppressed GFP fluorescence signals are highlighted by solid and hollow triangles, while improved and impaired leaf photobleaching phenotypes are denoted by corresponding squares, respectively. The predicted TFs regulating antiviral RNA silencing are marked by red squares.
Assessment of Viral Multiplication and Silencing Efficiency in Virus-Induced Gene Silencing Experiments
Transcript levels of 39 candidate TFs were examined in the uppermost leaves of petunia plants infected with TRV constructs. Results showed that all TF transcripts were significantly down-regulated in the corresponding TRV construct-infected petunia leaves (Supplementary Figure S4). Compared to TRV-GFP-PhPDS control, higher levels of GFP fluorescence intensity were observed with the down-regulation of PhARF11, PhbHLH41, PhbHLH93, PhZPT4-3, PhCOL4, PhGATA11, PhHSF-B3A, PhLOB1, PhNAC90, and PhWRKY75. On the contrary, silencing of PhPIF3, PhZAT10, PhERF22, PhHSF-B3B, PhHSF-A2, and PhWRKY40 led to lower fluorescence levels. TRV-RNA1 and TRV-RNA2 accumulation levels exhibited highly similar patterns to GFP fluorescence intensities (Figure 7A).
Figure 7. Virus accumulation and silencing efficiency in various transcription factors-silenced petunia leaves. (A) Relative GFP fluorescence intensities and TRV (RNA1 and RNA2) accumulation levels in the inoculated leaves at 4 days post inoculation (dpi) with Agrobacterium containing 39 transcription factors (TFs)-inserted TRV-GFP-PhPDS plasmids. (B) Percentage of photobleached leaves and expression of PhPDS in the uppermost leaves at 14 dpi with Agrobacterium containing various constructs. 26S rRNA was used as an internal control. Error bars represent standard error of the mean from three biological replicates. Asterisks indicate statistical significance as determined by Student’s t-test at P < 0.05.
At the phenotype level, TRV-GFP-PhPDS/TGA2.1-, ZFP6-, and ERF22-infected plants displayed remarkably higher rates of photobleached leaves than the control (TRV-GFP-PhPDS). In contrast, TRV-GFP-PhPDS/bHLH41-, bHLH93-, ZPT4-3-, ZAT10-, COL4-, BBX20-, ERF3-, KD1-3-, HSF-B3A-, HSF-B3B-, NAC90-, JUB1B-, and WRKY75-infiltrated plants demonstrated much lower rates of photobleached leaves. Compared to the control, the transcription of PhPDS was lower in the leaves with enhanced photobleaching, but higher in the leaves with reduced photobleaching (Figure 7B), which is consistent with the leaf phenotype.
PhCOL4 Modulates the Expression of RNA Silencing-Related Genes
In view of the prominent impact of PhCOL4 (SGN accession no. Peaxi162Scf00045g01824) silencing on GFP fluorescence and leaf photobleaching in the VIGS assay, we selected it for further functional characterization. Expression levels of PhCOL4 were examined in its silenced and transiently overexpressed petunia leaves. At 4 dpi with TRV-GFP-PhPDS/COL4, petunia leaves exhibited decreased transcript levels of PhCOL4 and several RNA silencing-related genes, including PhRDR6, PhDCL2, PhDCL4, PhAGO2, and PhAGO4, compared to TRV-GFP-PhPDS control (Figure 8A). Transient overexpression of PhCOL4 in petunia leaves resulted in increased transcript abundances of PhRDR6, PhDCL2, PhDCL4, PhAGO2, and PhAGO4, compared to empty vector control (Figure 8B).
Figure 8. Impact of PhCOL4 silencing and transient overexpression on RNA silencing-related genes. (A) RT-qPCR analysis of transcript levels of PhCOL4 and RNA silencing-related genes in the inoculated leaves at 4 days post inoculation (dpi) with Agrobacterium containing TRV-GFP-PhPDS or TRV-GFP-PhPDS/COL4. (B) RT-qPCR analysis of transcript levels of PhCOL4 and RNA silencing-related genes in the infiltrated leaves at 6 dpi with Agrobacterium containing pCAMBIA2300-PhCOL4 (pCAM2300-PhCOL4). The empty vector pCAMBIA2300 (pCAM2300) was used as the control. Expression levels were standardized to 26S rRNA. Error bars represent standard error of the mean from three biological replicates. Significance of difference was calculated by Student’s t-test at P < 0.05 and indicated by asterisks.
DAP-Seq Analysis of PhCOL4 Target Genes
To screen downstream target genes of PhCOL4, a DAP-Seq assay was performed. A total of 94,863,846 high-quality reads were generated from PhCOL4-bound products, compared to input control (Figure 9A). A number of peaks were located before transcription start sites of genes (Figure 9B). Of the 1,626 peaks, 660 (40.59%) were located within the 2kb upstream promoter regions. The remaining peaks were located within the regions of 5′UTR, 3′UTR, exon, intron, and downstream (2 kb after transcription end sites) (Figure 9C). Analysis of MEME and DREME revealed two top PhCOL4-bound motifs with the core element TTCTT (or AAGAA in the opposite strand) (Figure 9D). Through the functional annotation of peaks, two important genes, PhRDR6 and PhAGO4, involved in the RNA silencing pathway were found (Figure 9E).
Figure 9. Binding peaks and motifs of PhCOL4 by DNA affinity purification sequencing in petunia. (A) Abundance of reads from PhCOL4-bound DNA affinity purification sequencing (DAP-Seq) database. (B) Distribution of summits in peaks around transcription start sites (TSS) and end sites (TES) of genes. Summit positions within transcribed regions are shaded. (C) Number of peaks in different genomic categories. The promoter indicates the 2 kb region upstream of TSS, and the downstream represents the 2 kb region downstream of TES. (D) Two top DNA motifs bound by PhCOL4 from MEME and DREME analyses of DAP-Seq database. (E) Two putative PhCOL4-bound peaks and corresponding RNA silencing-related genes.
PhCOL4 Binds to the PhRDR6 and PhAGO4 Promoters
To determine the transactivation of PhRDR6 and PhAGO4 promoters by PhCOL4, a dual luciferase assay based on the effector and reporter constructs was carried out (Figure 10A). Compared to empty vector control, the co-expression of 35S:PhCOL4 with pPhRDR6:LUC or pPhAGO4:LUC caused a dramatic increase in LUC activity (Figure 10B). To further confirm the interaction of PhCOL4 with predicted motifs by DAP-Seq, a 35-bp fragment containing motif sequence in the PhRDR6 or PhAGO4 promoter was used as a probe for yeast one-hybrid assay. The constructs of promoter and PhCOL4 were used as the bait and pray, respectively (Figure 10C). The yeast cells co-transformed with PhCOL4 and wild-type promoter fragments grew well on the selective medium in the presence of 100 mM 3-AT. However, the growth of yeast transformants with mutant promoter fragments were remarkably inhibited (Figure 10D). These results revealed the direct interaction between PhCOL4 and the promoters of PhRDR6 and PhAGO4.
Figure 10. Transactivation of PhRDR6 and PhAGO4 promoters by PhCOL4. (A) Schematic diagrams of the effector and reporter constructs for dual luciferase assay. REN, Renilla luciferase; LUC, firefly luciferase. (B) Dual luciferase assay of the pPhRDR6 and pPhAGO4. The activation was expressed as a ratio of LUC to REN. Error bars indicate standard error of the mean from three biological replicates. Statistical significance was determined using Student’s t-test (P < 0.05) and shown as asterisks. (C) Schematic diagrams of the bait and pray constructs for yeast one-hybrid assay. The fragments containing wild-type (wt) and mutant (m) PhCOL4-bound motifs were used for construction of bait vectors. (D) Growth of yeast cells transformed with the bait and pray plasmids on SD-Ura-His-Leu plates supplemented with 100 mM 3-aminotriazole (3-AT) at different dilutions.
Discussion
RNA silencing is an efficient antiviral mechanism in plants, which can suppress the infection of various viruses (Yang and Li, 2018; Guo et al., 2019). The lack of an effective screening system has limited our ability to identify the regulatory machinery that controls antiviral RNA silencing. In this study, we adopted an integrative approach of high-depth RNA-Seq and transient VIGS to screen important TFs regulating the RNA silencing process. This approach has the advantages of fast detection, high throughput, and excellent selectivity. It has been widely employed to characterize the biotic stress-defensive genes in Zea mays (Linde and Doehlemann, 2013), Solanum lycopersicum (Paz et al., 2011; Zuluaga et al., 2013), Gossypium hirsutum (Gao et al., 2016; Zhang et al., 2016), N. benthamiana (Luo et al., 2019), and Rosa chinensis (Gao et al., 2021). A recent publication reported that two functional genes, GmSEOB and GmPAP27, were identified to confer resistance to soybean mosaic virus through integrative method of transcriptome and VIGS (Sun et al., 2021).
Endogenous hormones have proven to be tuners of plant response to environmental stimuli and biotic stresses, such as virus infection. In this study, we found increased production of ET, ABA, JA, and SA in TRV (PPK20)-infected petunia leaves (Figure 3A), suggesting their potential roles in modulating antiviral defense. These hormones have been revealed to affect plant resistance to various viruses, including tobacco necrosis virus (Iriti and Faoro, 2008), PVX (Lee et al., 2011), tomato bushy stunt virus (TBSV) (Sansregret et al., 2013), TMV (Zhu et al., 2014), bamboo mosaic virus (BaMV) (Alazem et al., 2017), and plum pox virus (PPV) (Pasin et al., 2020). These studies suggest that ABA, JA, and SA may have positive contribution to viral defense, whereas the function of ET varies across species infected by different types of viruses. For example, the mutant of ACS1, a key ET biosynthetic gene, showed reduced susceptibility to TMVcg (Chen et al., 2013). However, a contrasting report showed that the accumulation of TBSV was dramatically disrupted in ET-sensitive but not in ET-insensitive plants (Sansregret et al., 2013). Our data showed that a rapid production of ABA and JA occurred at 1 and 2 dpi with TRV (PPK20) (Figure 3A). Previous studies have revealed that ABA- and JA-mediated plant defense against viral attack occurs at early stages of infection (Asselbergh et al., 2008; García-Marcos et al., 2013). Their defensive effects can be compromised at later phases, due to their antagonistic interactions with SA (Alazem and Lin, 2015). This may explain the divergent induction profiles between ABA/JA and other hormones during the initial stage of TRV (PPK20) infection.
It is worth mentioning that the plants respond differently to virus infection between local and systemic tissues, especially of endogenous hormones. For example, SA is crucial for the establishment of local and systemic resistance to viral attack. The infection with TMV leads to a remarkable elevation in SA levels in both inoculated and systemic leaves of tobacco (Vlot et al., 2009). ABA is active in the inoculated Arabidopsis leaves with BaMV, whereas this activeness is only found in the systemic leaves upon infection with TMVcg (Chen et al., 2013). Obuda pepper virus infection significantly increases the accumulation of some hormones, such as SA, ABA, and JA, in the inoculated leaves rather than the systemic leaves of pepper (Dziurka et al., 2016). These results suggest that the hormones may play variable roles in the defense against virus invasion in different plant tissues. This possibility requires further examination in the subsequent studies.
There is increasing evidence showing the participation of hormones in RNA silencing machinery. It has been reported that ABA modulates the transcription of AGOs, and ABA-elicited defense response against BaMV depends on AGO2 and AGO3. Mutation of the ABA biosynthesis pathway attenuates the expression of AGO2 and accelerates the replication and systemic movement of PVX (Alazem et al., 2017). Expression of NahG, a gene contributing to the degradation of SA, decreases the abundances of PPV-derived small RNAs (Alamillo et al., 2006). SA treatment results in elevated transcript levels of RDR1 in tobacco (Hunter et al., 2013) and RDR1, RDR2, DCL1, and DCL2 in tomato (Campos et al., 2014). Some TFs associated with ABA and SA production share similar expression patterns with a few RNA silencing-related genes. These TF-bound cis-acting elements could be detected in the promoter regions of some RDRs, DCLs, or AGOs (Alazem et al., 2019). In addition, a rice TF JAMYB participates in the JA signaling pathway and specifically targets the downstream gene AGO18, a core RNA silencing component imparting resistance to rice stripe virus (Yang et al., 2020). This implies that ABA, SA, and JA may regulate different components of the RNA silencing machinery. Although the direct link between ET and antiviral RNA silencing has not been described, our previous work suggested that PhERF2, an ethylene-responsive element binding factor, regulated gene silencing efficiency and resistance to TRV and CMV infections (Sun et al., 2016). Another example is that ET-inducible TF RAV2 is essential for HC-Pro and P38 (two plant viral suppressors of silencing) to inhibit the silencing process (Endres et al., 2010). These observations suggest that the involvement of ET in RNA silencing-based antiviral defense is a complicated process. An extensive investigation of the impact of hormone treatment on expression of RNA silencing-related genes will be required in future work.
In the present study, we verified a number of TRV (PPK20)-induced TFs for their effects on GFP-labeled virus accumulation and PhPDS-silenced photobleaching phenotype. These TFs may serve as putative transcriptional regulators of antiviral RNA silencing (Figure 6). We found that these TFs belong to a wide range of families. Silencing of PhARF11, PhbHLH41, PhbHLH93, PhZPT4-3, PhCOL4, PhGATA11, PhHSF-B3A, PhLOB1, PhNAC90, and PhWRKY75 led to increased fluorescence intensities and TRV accumulation. Consistently, PhbHLH41, PhbHLH93, PhZPT4-3, PhCOL4, PhHSF-B3A, PhNAC90, and PhWRKY75 silencing greatly reduced photobleaching phenotypes (Figures 6, 7). Therefore, these 7 TFs can be identified as positive regulators of antiviral RNA silencing. Moreover, the other TFs, whose down-regulation increased virus accumulation, may play their antiviral roles beyond the RNA silencing machinery. In contrast to the above TFs, petunia leaves with down-regulation of PhPIF3, PhZAT10, PhERF22, PhERF1B, PhHSF-B3B, PhHSF-A2, and PhWRKY40 displayed decreased GFP-monitored virus accumulation, implying their negative roles in modulating antiviral responses. Considering the enhanced photobleaching development in TRV-GFP-PhPDS/ERF22-infiltrated leaves, PhERF22 may function as a suppressor of antiviral RNA silencing. Numerous studies have indicated that viruses usually hijack the host cellular machinery or components to evade plant’s immunity system and facilitate their infection (Patarroyo et al., 2013; Fang et al., 2019). A recent example showed that beet severe curly top virus can activate the Arabidopsis gene VIM5, encoding an E3 ubiquitin ligase, to repress the RNA silencing process (Benoit, 2020). It remains to be investigated whether a cross-talk between TRV proteins and petunia PhERF22 exists.
It is noteworthy that the down-regulation of a few TFs only led to altered photobleaching phenotype but no change to GFP intensity, i.e., virus accumulation. These TFs included PhBBX20, PhERF3, PhKD1-3, PhTGA2.1, and PhZFP6 (Figure 6). We reasoned that this phenomenon may be attributed to the environment–plant interactions and intra-plant variation of RNA silencing-based traits. It has been reported that the VIGS efficiency probably changes under certain environmental factors. Low temperature and humidity can increase TRV-induced gene silencing efficiency in tomato (Fu et al., 2006). In previous studies, we have also observed variable silencing phenotypes in various TRV-PhPDS-infected petunia plants (Sun et al., 2017). These findings demonstrate that the spread of virus or virus-derived silencing signals can be influenced by external factors, leading to varied levels of virus infection. Further study is necessary to investigate the potential impact of these external factors on virus infection.
In this work, we found a zinc finger TF, PhCOL4, plays an essential role in antiviral RNA silencing by specifically activating PhRDR6 and PhAGO4, two important components in the RNA silencing pathway (Figure 10). Some reports have revealed the important functions of PhRDR6 (Li et al., 2012) and PhAGO4 (Carbonell and Carrington, 2015; Kenesi et al., 2021) in RNA silencing-mediated antiviral defense. In Arabidopsis, AtCOL4, a homolog of PhCOL4, has been suggested to function as a modulator of flowering time (Steinbach, 2019) and abiotic stress tolerance (Min et al., 2015). Our results provided the direct evidence for the involvement of PhCOL4 in the antiviral RNA silencing process, thus corroborating the effectiveness of integrative approach of transcriptome and VIGS in identifying transcriptional regulators. However, how PhCOL4, PhRDR6, and PhAGO4 synergistically modulate plant defense against virus is still unknown. Future studies should include an investigation of antiviral roles of PhCOL4 and its target genes in petunia.
Conclusion
We adopted an integrative approach of high-depth RNA-Seq and VIGS assay to investigate the regulatory mechanism of antiviral RNA silencing in petunia. Our results showed that TRV (PPK20) infection induced significant transcriptional changes in petunia, which involve various TFs and RNA silencing-related genes. A number of TFs possibly regulating antiviral RNA silencing were identified through the TRV-GFP-PhPDS-based VIGS screening system. One potential positive regulator, PhCOL4, was further validated for its crucial role in antiviral RNA silencing by specifically activating the PhRDR6 and PhAGO4 promoters. These transcriptional regulators may be used for future genetic engineering of virus resistance in this important ornamental plant.
Data Availability Statement
The RNA-Seq and DAP-Seq data are available at GenBank (https://www.ncbi.nlm.nih.gov/genbank/) with accession BioProject PRJNA693880 and PRJNA693899, respectively.
Author Contributions
YX and DS designed the research and wrote the manuscript. YX, XJ, ZX, YY, and XC carried out the experiments. DS and DK analyzed the RNA-Seq and DAP-Seq data. YX performed the graphic arrangement. DS and YZ revised the manuscript and improved the English. All the authors contributed to the article and approved the submitted version.
Funding
This work was funded by the National Natural Science Foundation of China (Grant no. 31801895), Shaanxi Natural Science Basic Research Plan Project (Grant no. 2019JQ-417), Postdoctoral Special Funding Project of China (Grant no. 2019T120958), and Shaanxi Postdoctoral Research Project (Grant no. 2018BSHYDZZ72).
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.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We are appreciative of Xiaoning Luo’s work for the detection of endogenous hormones. We are grateful to Lingling Zeng for technical assistance on the measurement of GFP fluorescence signals. We thank Mengchen Li for the advice on data statistical analysis. We appreciate Ayla Norris’s kind help for editing this manuscript. We also appreciate the helpful and invaluable comments from editor and reviewers for improving this manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.876428/full#supplementary-material
Footnotes
- ^ https://solgenomics.net/organism/Petunia_axillaris/genome
- ^ http://wego.genomics.org.cn/cgi-bin/wego/index.pl
- ^ http://www.genome.jp/kegg/
- ^ http://www.premierbiosoft.com/primerdesign/
- ^ https://solgenomics.net/organism/Petunia_axillaris/genome
References
Alamillo, J. M., Saénz, P., and García, J. A. (2006). Salicylic acid-mediated and RNA-silencing defense mechanisms cooperate in the restriction of systemic spread of plum pox virus in tobacco. Plant J. 48, 217–227. doi: 10.1111/j.1365-313X.2006.02861.x
Alazem, M., and Lin, N. S. (2015). Roles of plant hormones in the regulation of host-virus interactions. Mol. Plant Pathol. 16, 529–540. doi: 10.1111/mpp.12204
Alazem, M., He, M. H., Moffett, P., and Lin, N. S. (2017). Abscisic acid induces resistance against bamboo mosaic virus through argonaute 2 and 3. Plant Physiol. 174, 339–355. doi: 10.1104/pp.16.00015
Alazem, M., Kim, K. H., and Lin, N. S. (2019). Effects of abscisic acid and salicylic acid on gene expression in the antiviral RNA silencing pathway in Arabidopsis. Int. J. Mol. Sci. 20:2538. doi: 10.3390/ijms20102538
Aparicio, F., and Pallás, V. (2017). The coat protein of alfalfa mosaic virus interacts and interferes with the transcriptional activity of the bHLH transcription factor ILR3 promoting salicylic acid-dependent defence signalling response. Mol. Plant Pathol. 18, 173–186. doi: 10.1111/mpp.12388
Asselbergh, B., De Vleesschauwer, D., and Höfte, M. (2008). Global switches and fine-tuning-ABA modulates plant pathogen defense. Mol. Plant-Microbe 21, 709–719. doi: 10.1094/MPMI-21-6-0709
Azevedo, J., Garcia, D., Pontier, D., Ohnesorge, S., Yu, A., Garcia, S., et al. (2010). Argonaute quenching and global changes in dicer homeostasis caused by a pathogen-encoded GW repeat protein. Gene. Dev. 24, 904–915. doi: 10.1101/gad.1908710
Benoit, M. (2020). The great escape: how a plant DNA virus hijacks an imprinted host gene to avoid silencing. Plant Cell 32, 3051–3052. doi: 10.1105/tpc.20.00631
Bologna, N. G., and Voinnet, O. (2014). The diversity, biogenesis, and activities of endogenous silencing small RNAs in Arabidopsis. Annu. Rev. Plant Biol. 65, 473–503. doi: 10.1146/annurev-arplant-050213-035728
Campos, L., Granell, P., Tárraga, S., López-Gresa, P., Conejero, V., Bellés, J. M., et al. (2014). Salicylic acid and gentisic acid induce RNA silencing-related genes and plant resistance to RNA pathogens. Plant Physiol. Biochem. 77, 35–43. doi: 10.1016/j.plaphy.2014.01.016
Carbonell, A., and Carrington, J. C. (2015). Antiviral roles of plant argonautes. Curr. Opin. Plant Biol. 27, 111–117. doi: 10.1016/j.pbi.2015.06.013
Chen, L., Zhang, L., Li, D., Wang, F., and Yu, D. (2013). WRKY8 transcription factor functions in the TMV-cg defense response by mediating both abscisic acid and ethylene signaling in Arabidopsis. Proc. Natl. Acad. Sci. U S A. 110, 1963–1971. doi: 10.1073/pnas.1221347110
Chen, Y., Lun, A. T., and Smyth, G. K. (2014). Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21, 3674–3676. doi: 10.1093/bioinformatics/bti610
Conesa, A., Götz, S., García-Gómez, J., Terol, J., Talón, M., and Robles, M. (2005). Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21, 3674–3676. doi: 10.1093/bioinformatics/bti610
Csorba, T., and Burgyán, J. (2016). “Antiviral silencing and suppression of gene silencing in plants,” in Current Research Topics in Plant Virology, eds A. Wang and X. Zhou (Cham: Springer International Publishing).
Deleris, A., Gallego-Bartolome, J., Bao, J., Kasschau, K. D., Carrington, J. C., and Voinnet, O. (2006). Hierarchical action and inhibition of plant dicer-like proteins in antiviral defense. Science 313, 68–71. doi: 10.1126/science.1128214
Dziurka, M., Janeczko, A., Juhász, C., Gullner, G., Oklestková, J., Novák, O., et al. (2016). Local and systemic hormonal responses in pepper leaves during compatible and incompatible pepper-tobamovirus interactions. Plant Physiol. Bioch. 109, 355–364. doi: 10.1016/j.plaphy.2016.10.013
Eamens, A., Wang, M. B., Smith, N. A., and Waterhouse, P. M. (2008). RNA silencing in plants: yesterday, today, and tomorrow. Plant Physiol. 147, 456–468. doi: 10.1104/pp.108.117275
Endres, M. W., Gregory, B. D., Gao, Z., Foreman, A. W., Mlotshwa, S., Ge, X., et al. (2010). Two plant viral suppressors of silencing require the ethylene-inducible host transcription factor RAV2 to block RNA silencing. PLoS Pathog. 6:e1000729. doi: 10.1371/journal.ppat.1000729
Fang, X. D., Yan, T., Gao, Q., Cao, Q., Gao, D. M., Xu, W. Y., et al. (2019). A cytorhabdovirus phosphoprotein forms mobile inclusions trafficked on the actin/ER network for viral RNA synthesis. J. Exp. Bot. 70, 4049–4062. doi: 10.1093/jxb/erz195
Fischer, U., and Dröge-Laser, W. (2004). Overexpression of NtERF5, a new member of the tobacco ethylene response transcription factor family enhances resistance to tobacco mosaic virus. Mol. Plant-Microbe 17, 1162–1171. doi: 10.1094/MPMI.2004.17.10.1162
Freeborough, W., Gentle, N., and Rey, M. E. (2021). WRKY transcription factors in cassava contribute to regulation of tolerance and susceptibility to cassava mosaic disease through stress responses. Viruses 13:1820. doi: 10.3390/v13091820
Fu, D. Q., Zhu, B. Z., Zhu, H. L., Zhang, H. X., Xie, Y. H., Jiang, W. B., et al. (2006). Enhancement of virus-induced gene silencing in tomato by low temperature and low humidity. Mol. Cells 21, 153–160.
Gaffar, F. Y., and Koch, A. (2019). Catch me if you can! RNA silencing-based improvement of antiviral plant immunity. Viruses 11:673. doi: 10.3390/v11070673
Gao, P., Zhang, H., Yan, H., Wang, Q., Yan, B., Jian, H., et al. (2021). RcTGA1 and glucosinolate biosynthesis pathway involvement in the defence of rose against the necrotrophic fungus Botrytis cinerea. BMC Plant Biol. 21:223. doi: 10.1186/s12870-021-02973-z
Gao, W., Long, L., Xu, L., Lindsey, K., Zhang, X., and Zhu, L. (2016). Suppression of the homeobox gene HDTF1 enhances resistance to Verticillium dahliae and Botrytis cinerea in cotton. J. Integr. Plant Biol. 58, 503–513. doi: 10.1111/jipb.12432
Gao, Z., Liu, H. L., Daxinger, L., Pontes, O., He, X., Qian, W., et al. (2010). An RNA polymerase II-and AGO4-associated protein acts in RNA-directed DNA methylation. Nature 465, 106–109. doi: 10.1038/nature09025
García-Marcos, A., Pacheco, R., Manzano, A., Aguilar, E., and Tenllado, F. (2013). Oxylipin biosynthesis genes positively regulate programmed cell death during compatible infections with the synergistic pair potato virus X-potato virus Y and tomato spotted wilt virus. J. Virol. 87, 5769–5783. doi: 10.1128/JVI.03573-12
Garcia-Ruiz, H., Takeda, A., Chapman, E. J., Sullivan, C. M., Fahlgren, N., Brempelis, K. J., et al. (2010). Arabidopsis RNA-dependent RNA polymerases and dicer-like proteins in antiviral defense and small interfering RNA biogenesis during turnip mosaic virus infection. Plant Cell 22, 481–496. doi: 10.1105/tpc.109.073056
Grabherr, M. G., Haas, B. J., Yassour, M., Levin, J. Z., Thompson, D. A., Amit, I., et al. (2011). Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 29, 644–652. doi: 10.1038/nbt.1883
Guo, Z., Li, Y., and Ding, S. W. (2019). Small RNA-based antimicrobial immunity. Nat. Rev. Immunol. 19, 31–44. doi: 10.1038/s41577-018-0071-x
Harvey, J. J., Lewsey, M. G., Patel, K., Westwood, J., Heimstädt, S., Carr, J. P., et al. (2011). An antiviral defense role of AGO2 in plants. PLoS One 6:e14639. doi: 10.1371/journal.pone.0014639
Havecker, E. R., Wallbridge, L. M., Hardcastle, T. J., Bush, M. S., Kelly, K. A., Dunn, R. M., et al. (2010). The Arabidopsis RNA-directed DNA methylation argonautes functionally diverge based on their expression and interaction with target loci. Plant Cell 22, 321–334. doi: 10.1105/tpc.109.072199
Heinz, S., Benner, C., Spann, N., Bertolino, E., Lin, Y. C., Laslo, P., et al. (2010). Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol. Cell 38, 576–589. doi: 10.1016/j.molcel.2010.05.004
Henderson, I. R., Zhang, X., Lu, C., Johnson, L., Meyers, B. C., Green, P. J., et al. (2006). Dissecting Arabidopsis thaliana dicer function in small RNA processing, gene silencing and DNA methylation patterning. Nat. Genet. 38, 721–725. doi: 10.1038/ng1804
Herath, V., and Verchot, J. (2021). Insight into the bzip gene family in Solanum tuberosum: genome and transcriptome analysis to understand the roles of gene diversification in spatiotemporal gene expression and function. Int. J. Mol. Sci. 22:253. doi: 10.3390/ijms22010253
Hunter, L. J., Westwood, J. H., Heath, G., Macaulay, K., Smith, A. G., MacFarlane, S. A., et al. (2013). Regulation of RNA-dependent RNA polymerase 1 and isochorismate synthase gene expression in Arabidopsis. PLoS One 8:e66530. doi: 10.1371/journal.pone.0066530
Iriti, M., and Faoro, F. (2008). Abscisic acid is involved in chitosan-induced resistance to tobacco necrosis virus (TNV). Plant Physiol. Biochem. 46, 1106–1111. doi: 10.1016/j.plaphy.2008.08.002
Jaubert, M., Bhattacharjee, S., Mello, A. F., Perry, K. L., and Moffett, P. (2011). Argonaute2 mediates RNA-silencing antiviral defenses against potato virus X in Arabidopsis. Plant Physiol. 156, 1556–1564. doi: 10.1104/pp.111.178012
Kenesi, E., Lopez-Moya, J. J., Orosz, L., Burgyán, J., and Lakatos, L. (2021). Argonaute2 controls antiviral activity against sweet potato mild mottle virus in Nicotiana benthamiana. Plants 10:867. doi: 10.3390/plants10050867
Lee, W. S., Fu, S. F., Verchot-Lubicz, J., and Carr, J. P. (2011). Genetic modification of alternative respiration in Nicotiana benthamiana affects basal and salicylic acid-induced resistance to potato virus X. BMC Plant Biol. 11:41. doi: 10.1186/1471-2229-11-41
Li, S., Yu, F., Wang, M., Guo, X., and Li, H. (2012). Molecular characterization of a Nicotiana tabacum NtRDR6 gene. Plant Mol. Biol. Rep. 30, 1375–1384. doi: 10.1007/s11105-012-0455-1
Li, T., Zhang, X. Y., Huang, Y., Xu, Z. S., Wang, F., and Xiong, A. S. (2018). An R2R3-MYB transcription factor, SlMYB28, involved in the regulation of TYLCV infection in tomato. Sci. Hortic. 237, 192–200. doi: 10.1016/j.scienta.2018.04.020
Linde, K., and Doehlemann, G. (2013). Utilizing virus-induced gene silencing for the functional characterization of maize genes during infection with the fungal pathogen Ustilago maydis. Methods Mol. Biol. 975:47. doi: 10.1007/978-1-62703-278-0_4
Liu, L., and Chen, X. (2016). RNA quality control as a key to suppressing RNA silencing of endogenous genes in plants. Mol. Plant 9, 826–836. doi: 10.1016/j.molp.2016.03.011
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
Llave, C. (2010). Virus-derived small interfering RNAs at the core of plant-virus interactions. Trends Plant Sci. 15, 701–707. doi: 10.1016/j.tplants.2010.09.001
Lu, J., Du, Z. X., Kong, J., Chen, L. N., Qiu, Y. H., Li, G. F., et al. (2012). Transcriptome analysis of Nicotiana tabacum infected by cucumber mosaic virus during systemic symptom development. PLoS One 7:e43447. doi: 10.1371/journal.pone.0043447
Luo, C., Wang, Z. Q., Liu, X., Zhao, L., Zhou, X., and Xie, Y. (2019). Identification and analysis of potential genes regulated by an alphasatellite (TYLCCNA) that contribute to host resistance against tomato yellow leaf curl China virus and its betasatellite (TYLCCNV/TYLCCNB) infection in Nicotiana benthamiana. Viruses 11:442. doi: 10.3390/v11050442
Ma, X., Nicole, M. C., Meteignier, L. V., Hong, N., Wang, G., and Moffett, P. (2015). Different roles for RNA silencing and RNA processing components in virus recovery and virus-induced gene silencing in plants. J. Exp. Bot. 66, 919–932. doi: 10.1093/jxb/eru447
Machanick, P., and Bailey, T. L. (2011). MEME-ChIP: motif analysis of large DNA datasets. Bioinformatics 27, 1696–1697. doi: 10.1093/bioinformatics/btr189
Manacorda, C. A., Tasselli, S., Marano, M. R., and Asurmendi, S. (2021). TuMV infection alters miR168/AGO1 and miR403/AGO2 systems regulation in Arabidopsis. bioRxiv [prprint] doi: 10.1101/2021.02.17.431672
Merkling, S. H., Overheul, G. J., van Mierlo, J. T., Arends, D., Gilissen, C., and van Rij, R. P. (2015). The heat shock response restricts virus infection in Drosophila. Sci. Rep. 5:12758. doi: 10.1038/srep12758
Min, J. H., Chung, J. S., Lee, K. H., and Kim, C. S. (2015). The CONSTANS-like 4 transcription factor, AtCOL4, positively regulates abiotic stress tolerance through an abscisic acid-dependent manner in Arabidopsis. J. Integr. Plant Biol. 57, 313–324. doi: 10.1111/jipb.12246
Mornya, P., and Cheng, F. (2011). The levels of hormone and carbohydrate in autumn and non-autumn flowering tree peonies. Can. J. Plant Sci. 91, 991–998. doi: 10.4141/cjps2011-043
Mußmann, V., Serek, M., and Winkelmann, T. (2011). Selection of transgenic petunia plants using the green fluorescent protein (GFP). Plant Cell Tiss. Org. 107, 483–492. doi: 10.1007/s11103-005-7704-8
Pasin, F., Shan, H., García, B., Müller, M., San León, D., Ludman, M., et al. (2020). Abscisic acid connects phytohormone signaling with RNA metabolic pathways and promotes an antiviral response that is evaded by a self-controlled RNA virus. Plant Commun. 1:100099. doi: 10.1016/j.xplc.2020.100099
Patarroyo, C., Laliberté, J. F., and Zheng, H. (2013). Hijack it, change it: how do plant viruses utilize the host secretory pathway for efficient viral replication and spread? Front. Plant Sci. 3:308. doi: 10.3389/fpls.2012.00308
Paz, E. Y., Rehman, S., Bai, Y., Bermudez, R. S., Borrás-Hidalgo, O., and Thomma, B. (2011). “Transcriptome analysis on verticillium-infected tomato to identify genes involved in host defense,” in EPS PhD, Autumn School’Host-Microbe Interactomics’, (Wageningen).
Pertea, G., Huang, X., Liang, F., Antonescu, V., Sultana, R., Karamycheva, S., et al. (2003). TIGR Gene Indices clustering tools (TGICL): a software system for fast clustering of large EST datasets. Bioinformatics 19, 651–652. doi: 10.1093/bioinformatics/btg034
Pumplin, N., and Voinnet, O. (2013). RNA silencing suppression by plant pathogens: defence, counter-defence and counter-counter-defence. Nat. Rev. Microbiol. 11, 745–760. doi: 10.1038/nrmicro3120
Qu, F., Ye, X., and Morris, T. J. (2008). Arabidopsis DRB4, AGO1, AGO7, and RDR6 participate in a DCL4-initiated antiviral RNA silencing pathway negatively regulated by DCL1. Proc. Natl. Acad. Sci. U S A. 105, 14732–14737. doi: 10.1073/pnas.0805760105
Ruiz-Ferrer, V., and Voinnet, O. (2009). Roles of plant small RNAs in biotic stress responses. Annu. Rev. Plant Biol. 60, 485–510. doi: 10.1146/annurev.arplant.043008.092111
Sansregret, R., Dufour, V., Langlois, M., Daayf, F., Dunoyer, P., Voinnet, O., et al. (2013). Extreme resistance as a host counter-counter defense against viral suppression of RNA silencing. PLoS Pathog. 9:e1003435. doi: 10.1371/journal.ppat.1003435
Sera, T. (2017). Use of peptide aptamers, cationic peptides and artificial zinc finger proteins to generate resistance to plant viruses. Curr. Opin. Virol. 26, 120–124. doi: 10.1016/j.coviro.2017.07.023
Singh, J., Mishra, V., Wang, F., Huang, H. Y., and Pikaard, C. S. (2019). Reaction mechanisms of Pol IV, RDR2, and DCL3 drive RNA channeling in the siRNA-directed DNA methylation pathway. Mol. Cell 75, 576–589. doi: 10.1016/j.molcel.2019.07.008
Steinbach, Y. (2019). The Arabidopsis thaliana constans-LIKE 4 (COL4)-a modulator of flowering time. Front. Plant Sci. 10:651. doi: 10.3389/fpls.2019.00651
Sun, D., Li, S., Niu, L., Reid, M. S., Zhang, Y., and Jiang, C. Z. (2017). PhOBF1, a petunia ocs element binding factor, plays an important role in antiviral RNA silencing. J. Exp. Bot. 68, 915–930. doi: 10.1093/jxb/erw490
Sun, D., Nandety, R. S., Zhang, Y., Reid, M. S., Niu, L., and Jiang, C. Z. (2016). A petunia ethylene-responsive element binding factor, PhERF2, plays an important role in antiviral RNA silencing. J. Exp. Bot. 67, 3353–3365. doi: 10.1093/jxb/erw155
Sun, D., Zhang, X., Zhang, Q., Ji, X., Jia, Y., Wang, H., et al. (2019). Comparative transcriptome profiling uncovers a Lilium regale NAC transcription factor, LrNAC35, contributing to defence response against cucumber mosaic virus and tobacco mosaic virus. Mol. Plant Pathol. 20, 1662–1681. doi: 10.1111/mpp.12868
Sun, T., Sun, X., Li, F., Ma, N., Wang, M., Chen, Y., et al. (2021). H2O2 mediates transcriptome reprogramming during soybean mosaic virus-induced callose deposition in soybean. Crop J. 10, 262–272. doi: 10.1016/j.cj.2021.04.005
Tian, J., Pei, H., Zhang, S., Chen, J., Chen, W., Yang, R., et al. (2014). TRV-GFP: a modified tobacco rattle virus vector for efficient and visualizable analysis of gene function. J. Exp. Bot. 65, 311–322. doi: 10.1093/jxb/ert381
Tian, Y. P., and Valkonen, J. P. (2013). Genetic determinants of potato virus Y required to overcome or trigger hypersensitive resistance to PVY strain group O controlled by the gene Ny in potato. Mol. Plant-Microbe 26, 297–305. doi: 10.1094/MPMI-09-12-0219-R
Trapnell, C., Roberts, A., Goff, L., Pertea, G., Kim, D., Kelley, D. R., et al. (2012). Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc. 7, 562–578. doi: 10.1038/nprot.2012.016
Trapnell, C., Williams, B. A., Pertea, G., Mortazavi, A., Kwan, G., Van Baren, M. J., et al. (2010). Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 28, 511–515. doi: 10.1038/nbt.1621
Unver, T., and Budak, H. (2009). Virus-induced gene silencing, a post transcriptional gene silencing method. Int. J. Plant Genom. 2009:198680. doi: 10.1155/2009/198680
Vaucheret, H. (2006). Post-transcriptional small RNA pathways in plants: mechanisms and regulations. Gene. Dev. 20, 759–771. doi: 10.1101/gad.1410506
Vlot, A. C., Dempsey, D. A., and Klessig, D. F. (2009). Salicylic acid, a multifaceted hormone to combat disease. Annu. Rev. Phytopathol. 47, 177–206. doi: 10.1146/annurev.phyto.050908.135202
Wang, M. B., and Metzlaff, M. (2005). RNA silencing and antiviral defense in plants. Curr. Opin. Plant Biol. 8, 216–222. doi: 10.1016/j.pbi.2005.01.006
Wixon, J., and Kell, D. (2000). The Kyoto encyclopedia of genes and genomes-KEGG. Yeast (Chichester, England) 17, 48–55.
Yang, X., Caro, M., Hutton, S. F., Scott, J. W., Guo, Y., Wang, X., et al. (2014). Fine mapping of the tomato yellow leaf curl virus resistance gene Ty-2 on chromosome 11 of tomato. Mol. Breeding 34, 749–760. doi: 10.1007/s11032-014-0072-9
Yang, Z., and Li, Y. (2018). Dissection of RNAi-based antiviral immunity in plants. Curr. Opin. Virol. 32, 88–99. doi: 10.1016/j.coviro.2018.08.003
Yang, Z., Huang, Y., Yang, J., Yao, S., Zhao, K., Wang, D., et al. (2020). Jasmonate signaling enhances RNA silencing and antiviral defense in rice. Cell Host Microbe 28, 89–103. doi: 10.1016/j.chom.2020.05.001
Yao, J., Shen, Z., Zhang, Y., Wu, X., Wang, J., Sa, G., et al. (2020). Populus euphratica WRKY1 binds the promoter of H+-ATPase gene to enhance gene expression and salt tolerance. J. Exp. Bot. 71, 1527–1539. doi: 10.1093/jxb/erz493
Ye, J., Fang, L., Zheng, H., Zhang, Y., Chen, J., Zhang, Z., et al. (2006). WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 34, 293–297.
Zhang, H., Li, L., He, Y., Qin, Q., Chen, C., Wei, Z., et al. (2020). Distinct modes of manipulation of rice auxin response factor OsARF17 by different plant RNA viruses for infection. Proc. Natl. Acad. Sci. U S A. 117, 9112–9121. doi: 10.1073/pnas.1918254117
Zhang, W., Zhang, H., Qi, F., and Jian, G. (2016). Generation of transcriptome profiling and gene functional analysis in Gossypium hirsutum upon Verticillium dahliae infection. Biochem. Bioph. Res. Co. 473, 879–885. doi: 10.1016/j.bbrc.2016.03.143
Zhang, X., Yuan, Y. R., Pei, Y., Lin, S. S., Tuschl, T., Patel, D. J., et al. (2006). Cucumber mosaic virus-encoded 2b suppressor inhibits Arabidopsis argonaute1 cleavage activity to counter plant defense. Gene. Dev. 20, 3255–3268. doi: 10.1101/gad.1495506
Zheng, Y., Zhao, L., Gao, J., and Fei, Z. (2011). iAssembler: a package for de novo assembly of Roche-454/Sanger transcriptome sequences. BMC Bioinformatics 12:453. doi: 10.1186/1471-2105-12-453
Zhu, F., Xi, D. H., Yuan, S., Xu, F., Zhang, D. W., and Lin, H. H. (2014). Salicylic acid and jasmonic acid are essential for systemic resistance against tobacco mosaic virus in Nicotiana benthamiana. Mol. Plant-Microbe 27, 567–577. doi: 10.1094/MPMI-11-13-0349-R
Zhu, H., Hu, F., Wang, R., Zhou, X., Sze, S.-H., Liou, L. W., et al. (2011). Arabidopsis Argonaute10 specifically sequesters miR166/165 to regulate shoot apical meristem development. Cell 145, 242–256. doi: 10.1016/j.cell.2011.03.024
Zou, L. J., Deng, X. G., Han, X. Y., Tan, W. R., Zhu, L. J., Xi, D. H., et al. (2016). Role of transcription factor HAT1 in modulating Arabidopsis thaliana response to cucumber mosaic virus. Plant Cell Physiol. 57, 1879–1889. doi: 10.1093/pcp/pcw109
Keywords: petunia, transcriptome, virus-induced gene silencing, tobacco rattle virus, transcription factor
Citation: Xu Y, Ji X, Xu Z, Yuan Y, Chen X, Kong D, Zhang Y and Sun D (2022) Transcriptome Profiling Reveals a Petunia Transcription Factor, PhCOL4, Contributing to Antiviral RNA Silencing. Front. Plant Sci. 13:876428. doi: 10.3389/fpls.2022.876428
Received: 15 February 2022; Accepted: 21 March 2022;
Published: 14 April 2022.
Edited by:
Sachin Teotia, Sharda University, IndiaReviewed by:
Jinping Zhao, Texas A&M University, United StatesVed Prakash, Kansas State University, United States
Copyright © 2022 Xu, Ji, Xu, Yuan, Chen, Kong, Zhang and Sun. 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: Daoyang Sun, sundaoyang@nwsuaf.edu.cn
†These authors have contributed equally to this work