Abstract
“Candidatus Liberibacter asiaticus” (CLas) is a phloem-restricted Gram-negative bacterium that is the causal agent of citrus huanglongbing (HLB). In this study, we identified a CLas-encoded Sec-dependent secretory protein CLIBASIA_04405 that could contribute to the pathogenicity of this bacterium. The gene expression level of CLIBASIA_04405 was significantly higher in citrus than in psyllids. Transient overexpression of the mature CLIBASIA_04405 protein (m4405) in Nicotiana benthamiana leaves could suppress hypersensitive response (HR)-based cell death and H2O2 accumulation triggered by the mouse BAX and the Phytophthora infestans INF1. An alanine-substitution mutagenesis assay revealed the essential of amino acid clusters EKR45–47 and DE64–65 in cell death suppression. Challenge inoculation of the transgenic N. benthamiana-expressing m4405 with Pseudomonas syringae DC3000ΔhopQ1-1 demonstrated the greatly reduced bacterial proliferation. Remarkably, transcriptome profiling and RT-qPCR analysis disclosed that the gene expression of six small heat shock proteins (sHSPs), a set of plant defense regulators, were significantly elevated in the transgenic m4405 lines compared with those in wild-type N. benthamiana. In addition, the transgenic m4405 lines displayed phenotypes of dwarfism and leaf deformation. Altogether, these data indicated that m4405 was a virulence factor of CLas.
Introduction
Citrus huanglongbing (HLB), also known as citrus greening, is the most devastating worldwide disease of citrus at the moment and greatly damages the citrus industry by shortening lifespan of trees and reducing fruit yield and quality (Bové, 2006; Gottwald, 2010). The disease is associated with “Candidatus Liberibacter asiaticus” species that are phloem-limited, fastidious alpha-proteobacteria (Bové and Ayres, 2007; Wang and Trivedi, 2013). To date, three Ca. Liberibacter asiaticus species have been identified to be the causal agent of HLB, including “Candidatus Liberibacter asiaticus” (CLas), “Ca. L. africanus” (CLaf), and “Ca. L. americanus” (CLam) (Jagoueix et al., 1996; Wang et al., 2017). Among them, CLas is transmitted among citrus by Diaphorina citri Kuwayama [Asian citrus psyllid (ACP)] and is the most widespread and virulent species (Gottwald, 2010). So far, no citrus variety has been identified that possesses immunity to CLas (Folimonova et al., 2009; da Graça et al., 2016; Martinelli and Dandekar, 2017; Zheng et al., 2018). Given its serious impact on citrus production, CLas attracts increasing attention. However, the pathogenicity of CLas remains largely unknown mostly owing to its fastidiousness.
Bacterial pathogens often secrete proteins (also called effectors) that contribute to their pathogenicity (Hayes et al., 2010; Dou and Zhou, 2012). It is known that many Gram-negative bacteria depend on a type III secretion system (T3SS), T4SS, or T6SS to deliver effectors into the host cells and thus allow the manipulation of host cellular pathways to benefit the pathogen (Yamane et al., 2004; Wang et al., 2011). As an intracellular bacterium, CLas does not possess this conserved secretion mechanism (Duan et al., 2009). However, like phytoplasmas, another group of phloem-dwelling intracellular bacteria (MacLean et al., 2011; Sugio et al., 2011), CLas was inferred to secrete effector proteins to host cells via a Sec-dependent secretion system (Duan et al., 2009; Prasad et al., 2016).
A significant number of CLas-encoded proteins have been experimentally validated to be Sec-dependent secretory proteins (Prasad et al., 2016; Zhang et al., 2019), and their potential roles in CLas pathogenesis are beginning to be revealed. For example, CLIBASIA_03875 and CLIBASIA_04025 suppress programmed cell death (PCD) (Zhang et al., 2019; Pang et al., 2020), CLIBASIA_05315 targets the chloroplasts and induces cell death in Nicotiana benthamiana (Pitino et al., 2016, 2018) as well as physically interacting with the papain-like cysteine proteases (PL) and inhibiting their protease activity in citrus (Clark et al., 2018). In addition, a subset of Sec-dependent secretory proteins were deduced to be candidate effectors that play roles during CLas infection using a Tobacco mosaic virus-based expression system or a strategy for the analysis of temporal and spatial gene expression (Shi et al., 2019; Ying et al., 2019; Li et al., 2020). Recently, two non-classically secreted proteins, CLIBASIA_RS00445 and SC2_gp095, were shown to dramatically downregulate the gene expression of RbohB, a gatekeeper of H2O2-mediated defense signaling in Nicotiana tabacum (Jain et al., 2015, 2018), and the former was also identified to suppress oxylipin-mediated defense in citrus plants (Jain et al., 2019). These findings indicate that, apart from the Sec-dependent secretory proteins, CLas also employs non-classically secreted proteins to interfere with the plant defense mechanisms.
CLIBASIA_04405 was first annotated in CLas strain psy62 (Duan et al., 2009). However, no biological role has been linked to it yet. Here we characterized CLIBASIA_04405 as a Sec-dependent secretory protein of CLas and showed that its mature form (designated m4405) is able to suppress hypersensitive cell death triggered by both the mouse BAX (Lacomme and Santa Cruz, 1999) and the Phytophthora infestans INF1 (Kamoun et al., 1998), two well-known hypersensitive response (HR)-based PCD inducers. In addition, we showed that transgenic N. benthamiana expressing m4405 displayed improved resistance to Pseudomonas syringae DC3000ΔhopQ1-1 (Wei et al., 2007) as well as defects in growth and development compared with wild-type (WT) N. benthamiana. The study revealed an effector of CLas that might contribute to the bacterial pathogenicity.
Materials and Methods
In silico Analysis
The putative signal peptide (SP) of CLIBASIA_04405, namely, 4405SP, was predicted using SignalP (version 5.0) with default settings for Gram-negative bacteria.
Alkaline Phosphatase (PhoA) Assay
The coding sequence of 4405SP was amplified with the primers 4405F/4405R (Supplementary Table 1) and cloned into pET-mphoA (Liu et al., 2019) that was double digested with Nde I and Hind III. The resulting construct pET-4405SP-mphoA that harbored an in-frame gene fusion between 4405SP and mphoA was transformed into Escherichia coli BL21. To detect the alkaline phosphatase (PhoA) activities, the transformants were randomly selected and incubated on the indicator Luria–Bertani agar (90 μg/ml BCIP, 100 mM IPTG, and 75 mM Na2HPO4) at 37°C, as previously described (Liu et al., 2019). After 24 h of incubation, the transformants with a color change to blue indicate PhoA activity.
Detection of the Transcripts for CLIBASIA_04405 in CLas-Infected Citrus and Psyllids
Total RNA was extracted from CLas-infected Valencia sweet orange (Citrus sinensis) seedlings and CLas-bearing psyllids (Liu et al., 2019). Using the specific primers of CLIBASIA_04405 (Supplementary Table 1), reverse transcription-quantitative polymerase chain reaction (RT-qPCR) was conducted as previously described (Liu et al., 2019). Three biological replicates were carried out, using the CLas LasgyrA gene (GenBank No. CP001677.5) as an internal reference. The relative gene expression was determined using the 2–ΔΔCt method (Livak and Schmittgen, 2001), and statistical analyses were performed using the Student’s t-test.
Subcellular Localization of m4405 in N. benthamiana
The coding sequences of CLIBASIA_04405 and m4405 amplified with the specific primers (Supplementary Table 1) were inserted into pCAMBIA1300-35S-GFP digested with Kpn I and Xho I. The resulting p4405-GFP and pm4405-GFP were individually introduced into Agrobacterium tumefaciens EHA105, followed by co-infiltration with agrobacterial cells carrying a plasma membrane marker pm-rk CD3-1007 (Nelson et al., 2007) into leaves of the 4-week-old N. benthamiana, as previously described (Van der Hoorn et al., 2000). At 60 h post inoculation (hpi), the infiltrated patches were visualized for GFP fluorescence under a LSM700 confocal microscope (Zeiss, Germany).
Agrobacterium-Mediated PVX Infection Assay
The coding sequence of m4405 amplified with the primers m4405-F/m4405-R (Supplementary Table 1) was cloned into a Potato virus X (PVX)-based expression vector pGR107 (Jones et al., 1999) to generate pPVX-m4405. Additionally, nine alanine-substituted mutants of m4405 (Supplementary Figure 1) were produced by PCR using the primers listed in Supplementary Table 1 and individually inserted into pGR107, resulting in pPVX-A9-10, pPVX-A24-25, pPVX-A39-40, pPVX-A45-47, pPVX-A52-54, pPVX-A59-60, pPVX-A64-65, and pPVX-A75-78, pPVX-A93-94.
The PCD suppression assay was carried out as described previously (Wang et al., 2011). Briefly, the A. tumefaciens GV3101 that harbored pPVX-m4405, pPVX-A9-10, pPVX-A24-25, pPVX-A39-40, pPVX-A45-47, pPVX-A52-54, pPVX-A59-60, pPVX-A64-65, pPVX-A75-78, or pPVX-A93-94 were first infiltrated into the expanded leaves of the 6-week old N. benthamiana plants. The A. tumefaciens cells that harbored pPVX-GFP were included as a negative control. At 24 hpi, the infiltrated patches were further inoculated with the Agrobacterium cells that harbored pGR-BAX or pGR-INF1. At 2 days post second inoculation, one third of the infiltrated leaves were detached to perform histochemical staining with 3,3′-diaminobenzidine (DAB) as described (Vanacker et al., 2000). To record the cell death development, the rest leaves were continuously observed up to 5 days post second inoculation. The experiment was repeated three times, and each replicate included six to eight plants.
Generation of Transgenic N. benthamiana
The coding sequence of m4405 was inserted into a binary plasmid vector pBI121 (Chen et al., 2003), resulting in pBI-m4405 that was subsequently introduced into A. tumefaciens EHA105. The transgenic N. benthamiana lines expressing m4405 were generated using a leaf disk method (Gallois and Marinho, 1995). T1 generation seeds were screened on media containing 50 μg/ml of kanamycin (Kan), and the Kan-resistant seedlings were further confirmed by RT-PCR with the primers listed in Supplementary Table 1.
Pseudomonas syringae Infection Assays
The infection assay was performed as previously described (Wei et al., 2018). DC3000ΔhopQ1-1 (Wei et al., 2007), a knockout mutant of P. syringae pv. tomato strain DC3000 was cultured in King’s B (KB) broth with rifampicin (50 μg/ml) at 28°C overnight, and collected by centrifugation (2,500 × g, 5 min), followed by washing and resuspension in 10 mM MgCl2 to a final concentration of 3 × 104 CFU/ml. The leaves of 4-week-old N. benthamiana and transgenic m4405 lines were syringe infiltrated with the bacterial suspensions. At 6 days post inoculation (dpi), the infiltrated leaf tissues were collected and ground in distilled water. Serial 10-fold dilutions were spread on KB plates with rifampicin (50 μg/ml) and incubated at 28°C for 2 days. Bacterial counts were expressed as CFU/cm2 leaf disc and represented the mean ± SD (standard deviation) of three biological replicates. Statistical analysis was determined using Student’s t−test.
Transcriptome Sequencing and Analysis of Differentially Expressed Genes
Total RNA was extracted from leaves of the 4-week-old N. benthamiana and transgenic m4405 lines by using an RNeasy Mini kit (Qiagen, United States), followed by cDNA library preparation and in-depth sequencing using an Illumina HiSeq2000. The resulting clean reads were assembled into transcripts using the Trinity method (Grabherr et al., 2011). After removing redundancy, the unigenes were generated and subjected to a blast search against the public protein databases of NR (NCBI non-redundant protein sequences), Swiss-Prot, NT (NCBI nucleotide sequences), and Protein family (Pfam) by using NCBI blastX with a cutoff E-value of ≤1e-5. The clean reads from N. benthamiana and the transgenic m4405 lines were individually mapped to the assembled transcriptome using RSEM with no more than two nucleotide mismatches (Li and Dewey, 2011), and the resulting read count for each matched unigene was further normalized using the reads per kilo bases per million read (RPKM) method (Mortazavi et al., 2008). The differentially expressed genes (DEGs) were determined using DEGseq (Wang et al., 2010) with threshold as “q-value < 0.005 and | log2 (foldchange)| > 1.” Heatmap was created using the function heatmap.2 of gplots R package (Warnes et al., 2015). The nucleotide sequences of raw reads and the assembled unigene were submitted to NCBI Sequence Read Archive (SRA) database under the accession number PRJNA663964.
Detection of the Transcripts for sHSPs in the Transgenic m4405 N. benthamiana
Total RNA from leaves of the 4-week-old N. benthamiana and transgenic m4405 lines was subjected to RT-qPCR analysis as previously described (Zhang et al., 2019). The gene-specific primers used are listed in Supplementary Table 1. The experiment was performed with three biological replicates, using the N. benthamiana actin (GenBank No. JQ256516.1) as a reference gene. The relative gene expression was determined using the 2–ΔΔCt method (Livak and Schmittgen, 2001), and statistical analysis was conducted using Student’s t-test.
Results
CLIBASIA_04405 Was a Sec-Dependent Secretory Protein That Exhibited Higher Transcriptional Level in Citrus Than in Psyllids
The CLIBASIA_04405 of CLas is composed of 121 amino acids (aa) (Duan et al., 2009). A BLAST search of the GenBank database revealed that this protein was 100% conserved among the CLas strains with the full genome available (Supplementary Figure 2), but was absent from both CLaf and CLam. An in silico analysis with SignalP indicated that CLIBASIA_04405 had 4405SP, a 20-aa SP at its N-terminus (Figure 1A). To functionally validate 4405SP, the E. coli-based phoA fusion assay was performed based on the chromogenic reaction of the secreted PhoA protein on LB medium containing BCIP substrate (Liu et al., 2019). After 24 h of incubation, the BL21 cells that expressed mphoA lacking its native SP, which served as a negative control, remained white, whereas the bacteria that expressed 4405SP-mPhoA turned dark blue (Figures 1B,C) after 6 h of incubation, indicating that 4405SP successfully directed the extracellular secretion of the mPhoA moiety. The in silico prediction combined with the experimental validation supported the concept that CLIBASIA_04405 had an SP that could direct the protein secretion through the Sec-dependent pathway.
FIGURE 1
We next detected the transcripts for CLIBASIA_04405 in host plants and insect vectors. RT-qPCR analysis was performed on total RNA extracted from the CLas-infected sweet orange and ACPs, and the result showed that the gene expression level of CLIBASIA_04405 in citrus was significantly higher (∼50-fold) than that in psyllids (Figures 2A,B).
FIGURE 2
m4405 Was Localized in Plasma Membrane of N. benthamiana Cells
Since CLIBASIA_04405 was a Sec-dependent secretory protein, we investigated the subcellular localization of its mature protein m4405 in plants. The m4405 coding sequence was fused with the GFP gene (Figure 3A) and expressed in N. benthamiana via Agrobacterium-mediated transient expression. The N-terminal GFP-tagged m4405 was detected under confocal fluorescence microscopy. As shown in Figure 3B, m4405-GFP was predominately colocalized with the plasma membranal marker pm-rk CD3-1007 (Nelson et al., 2007), indicating the plasma membrane localization of m4405 in plants. Using the same strategy, we observed that CLIBASIA_04405 shared a similar subcellular localization in N. benthamiana cells with m4405 (Figure 3B). This suggested that the SP of CLIBASIA_04405 had little or no effect on plasma membrane localization of the protein in plants.
FIGURE 3
m4405 Suppressed BAX- and INF1-Triggered Hypersensitive Cell Death in N. benthamiana
On the basis of a plant virus expression vector pGR107 (Jones et al., 1999), the mouse BAX and the P. infestans INF1 have been widely utilized to identify the PCD suppressor of phytopathogens (Bos et al., 2006; Li et al., 2015; Xiang et al., 2016; Fang et al., 2019). We employed this strategy to evaluate the ability of m4405 in PCD suppression. As shown in Figure 4A, m4405 completely suppressed both BAX- and INF1-triggered HR-based cell death in N. benthamiana, whereas cell death suppression had not been observed for GFP and the infiltration buffer.
FIGURE 4
To further understand the role of m4405 that leads to cell death suppression, the accumulation of H2O2, a reactive oxygen species (ROS) critical in plant cell death (Petrov and Van Breusegem, 2012), was analyzed in the N. benthamiana leaves co-expressing m4405 plus either BAX or INF1. Using a DAB staining method (Vanacker et al., 2000), we detected that m4405 completely inhibited both BAX- or INF1-stimulated H2O2 burst (Figure 4B). By contrast, neither GFP nor the infiltration buffer interfered with H2O2 accumulation.
To identify the key domain(s) of the m4405 involved in PCD suppression, nine clusters of charged amino acids (ED9–10, KE24–25, ER39–40, EKR45–47, EER52–54, RK59–60, DE64–65, DEDK75–78, and KK93–94) of m4405 were individually selected for alanine-substitution mutagenesis (Supplementary Figure 1). The resulting mutants were expressed through pGR107 to test their abilities of PCD suppression. As results, there were two mutants (A45-47 and A64-65), whose abilities to suppress PCD induced by INF1 and BAX were both greatly reduced or even lost completely (Table 1), indicating that EKR45–47 and DE64–65 are necessary for m4405 ability to suppress cell death.
TABLE 1
| Mutanta | PCD Inducerb | |
| BAX | INF1 | |
| m4405 | 0/135 | 0/99 |
| GFP | 135/135 | 99/99 |
| A9-10 | 1/45 | 9/33 |
| A24-25 | 1/45 | 21/33 |
| A39-40 | 1/45 | 22/33 |
| A45-47 | 45/45 | 32/33 |
| A52-54 | 0/45 | 30/33 |
| A59-60 | 0/45 | 12/33 |
| A64-65 | 22/45 | 31/33 |
| A75-78 | 6/45 | 16/33 |
| A93-94 | 10/45 | 28/33 |
Identification of the key amino acids of m4405 involved in PCD suppression.
aThe m4405 and the alanine-substituted m4405 mutants were transiently expressed in Nicotiana benthamiana via Agroinfiltration. bThe mouse BAX and the Phytophthora infestans INF1 were used as PCD inducers. The development of cell death was recorded at 5 days post infiltration of BAX or INF1. The PCD suppression of m4405 and its related mutants are shown as “lesion number/infiltration number.”
The Transgenic N. benthamiana Expressing m4405 Exhibited Improved Resistance to P. syringae
To gain further insight into the impact of m4405 on plant immunity, we generated transgenic N. benthamiana that stably expressed m4405 under the control of Cauliflower mosaic virus (CaMV) 35S promoter. A total of six independent transgenic lines were produced, and each of them exhibited phenotypes of dwarfism and leaf malformation. The lines m4405-18 and m4405-25 (Figure 5A) were then randomly selected for the following experiments. In the first experiment, RT-qPCR analysis was performed on total RNA of transgenic N. benthamiana and disclosed that the expression level of m4405 positively correlated with the phenotypic changes of the transgenic lines (Figures 5A,B). Second, using N. benthamiana as a control, the transgenic plants were challenged with DC3000ΔhopQ1-1, a P. syringae DC3000 mutant that lacks the type III effector HopQ1-1 and thus is able to grow in N. benthamiana instead of inducing an HR reaction (Wei et al., 2007, 2013). Measurement of growth of DC3000ΔhopQ1-1 in the inoculated leaf tissues at 6 dpi showed that the transgenic m4405 plants supported significantly less bacterial proliferation than WT (Figure 5C), indicating that the transgenic plants exhibited improved resistance to DC3000ΔhopQ1-1.
FIGURE 5
To address the underlying mechanism of resistance to DC3000ΔhopQ1-1 in the transgenic m4405 lines, we compared the transcriptome profiles of m4405-18 and m4405-25 with that of WT and identified an array of DEGs (Supplementary Table 2). Mining the DEGs revealed that a subset of genes encoding small heat shock proteins (sHSPs), including sHSP17.7, sHSP17.3, sHSP18.1a, sHSP18.1b, sHSP21, and sHSP22, were significantly upregulated in the transgenic lines compared with WT (Figure 6A). RT-qPCR analysis confirmed that the mRNAs for the six sHSPs were significantly accumulated in the transgenic lines (Figure 6B). The sHSPs are known to be a group of proteins with a molecular mass of 15 to 42 kDa that act as molecular chaperones to protect the plants against various biotic and abiotic stress conditions (Haslbeck and Vierling, 2015). Taking these results together, we inferred that the upregulated sHSPs in the transgenic m4405 N. benthamiana could contribute to the increased resistance to DC3000ΔhopQ1-1.
FIGURE 6
Discussion
Although the molecular basis of the CLas pathogenicity remains largely unknown, much progress has been made toward elucidating its mechanisms of pathogenesis, particularly through the identification and characterization of CLas-encoded effectors, such as the Sec-dependent secretory proteins (Pitino et al., 2016, 2018; Clark et al., 2018; Liu et al., 2019; Shi et al., 2019; Ying et al., 2019; Zhang et al., 2019; Li et al., 2020). In a previous study, a genome-wide in silico analysis on CLas psy62 strain predicted a total of 166 proteins containing SPs that may target them out of the cytoplasm via the Sec-dependent pathway (Prasad et al., 2016). However, 51 of these proteins, including CLIBASIA_04405, were determined to be false-positive predictions (Prasad et al., 2016) when experimentally validated using an early developed phoA fusion system (Mdluli et al., 1995). In this study, we revisited the putative SP of CLIBASIA_04405 with an updated phoA system (Liu et al., 2019) and provided evidence that CLIBASIA_04405 was instead a Sec-dependent secretory protein of CLas. Furthermore, we found that the expression level of CLIBASIA_04405 was significantly higher in a citrus host than in insect vectors. It is worthy to emphasize that a similar expression pattern was also observed in several CLas genes, such as CLIBASIA_00255, SC2_gp095, CLIBASIA_RS00445, CLIBASIA_00460, and CLIBASIA_03875, which have been suggested to function in the context of the CLas infection on the host plants (Jain et al., 2015, 2018; Li et al., 2017; Liu et al., 2019; Zhang et al., 2019). In addition, we show that m4405, the mature form of CLIBASIA_4405, completely suppresses PCD triggered by both BAX and INF1 in N. benthamiana. The data collectively indicates that CLIBASIA_4405 served as a virulence factor and was implicated in CLas pathogenesis in plant.
Programmed cell death plays roles in multiple cellular processes of plants, including immunity (Huysmans et al., 2017; Locato and De Gara, 2018). Upon pathogen infection, the plants may prime a form of PCD termed HR, a strong immune response that is accompanied by localized cell death to prevent the spread of the pathogen (Coll et al., 2011; Stael et al., 2015). In contrast, the pathogens, including bacteria, fungi, oomycetes, and nematodes, often secrete an arsenal of effectors that contribute to cell death suppression to enable their successful infection (da Cunha et al., 2007; Dou and Zhou, 2012; Goverse and Smant, 2014; Lo Presti et al., 2015). In CLas, the mature forms of two Sec-dependent secretory proteins CLIBASIA_03875 and CLIBASIA_04025, termed m3875 and SDE15, have been determined as PCD suppressors (Zhang et al., 2019; Pang et al., 2020). Herein, we demonstrated that m4405, the mature form of the Sec-dependent secretory protein CLIBASIA_04405, also inhibited hypersensitive cell death. Together, the findings highlight that, despite lacking the conserved secretion machineries such as T3SS (Duan et al., 2009), CLas employs the Sec-dependent secretory proteins to overcome HR of the citrus host.
To inhibit the development of hypersensitive cell death, pathogen-encoded PCD suppressors may interact with the host immune receptors (Bos et al., 2010; Wilton et al., 2010; Chen et al., 2012), or interfere with ROS accumulation that accompanies cell death (Dong et al., 2011; Hemetsberger et al., 2012; Zhang et al., 2015a). It has been shown that the CLas PCD suppressor m3875 regulates the gene expression of cell death regulators, cyclic nucleotide-gated channels (CNGC), Bax-inhibitor 1 (BI-1), and WRKY9 (Zhang et al., 2019), while another CLas PCD suppressor SDE15 directly interacts with a citrus protein ACCELERATED CELL DEATH 2 (ACD2), a negative regulator of PCD (Pang et al., 2020). In this study, how m4405 suppressed hypersensitive cell death remains elusive. However, alanine-substitution mutagenesis of m4405 showed that two amino acid clusters (EKR45–47 and DE64–65) were crucial for suppressing both BAX- and INF1-triggered cell death. Further study on the two loss-of-function mutants (A45-47 and A64-65) as well as identification of the potential target protein(s) in plants could be helpful to elucidate the mechanism by which m4405 suppressed PCD.
A number of pathogen effectors have been shown to mediate disease resistance when expressed transgenically (Kawamura et al., 2009; Rajput et al., 2015; Zhang et al., 2015b; Wang et al., 2019). Likewise, the study here displayed that the transgenic N. benthamiana expressing m4405 exhibited improved resistance to the bacterial pathogen DC3000ΔhopQ1-1. Transcriptome profiling and RT-qPCR revealed that the gene expression of six sHSPs (sHSP17.7, sHSP17.3, sHSP18.1a, sHSP18.1b, sHSP21, and sHSP22) were significantly upregulated in the transgenic lines. The sHSPs are evolutionarily conserved proteins in almost all organisms, ranging from bacteria and fungi to plants and animals (Waters and Vierling, 1999; Haslbeck and Vierling, 2015). The plant sHSPs are usually below detection limits in vegetative tissues under normal growth conditions but instead are induced by environmental stresses and developmental stimuli (Jacob et al., 2017; Carra et al., 2019). So far, there is emerging evidence that sHSPs serve as defense regulators and function in diverse defense mechanisms (Ul Haq et al., 2019). For example, the N. tabacum sHSP17 is required for basic immune responses other than HR-based resistance in plants (Maimbo et al., 2007). In contrast, sHSP20s plays a positive role in the hypersensitive defense response, and silencing of HSP20s in N. benthamiana can compromise the HR induced by the tomato resistance (R) genes I-2 and Mi-1 (Van Ooijen et al., 2010). In addition, a subset of the Arabidopsis sHSPs, including sHSP22, has been identified to be downregulated upon Pst DC3000 infection (Bricchi et al., 2012), implying their functions in response to the P. syringae infection. Remarkably, sHSP22 was one of the six sHSPs that were transcriptionally activated in the transgenic m4405 N. benthamiana. The potential of sHSP22, as well as those of the other five sHSPs, in P. syringae resistance merits further investigation.
In conclusion, the study identified a Sec-dependent secretory protein of CLas and showed that its mature form m4405 could suppress hypersensitive cell death in N. benthamiana. Transgenic overexpression of m4405 in N. benthamiana interfered with the plant growth and development and, in particular, upregulated the transcription of a group of sHSP genes. The data collectively indicated the potential role of m4405 in CLas infection and might shed new light on the pathogenicity of CLas.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
WL, XW, and CZ designed the experiments. PD, HY, CZ, and ZZ performed the experiments. WL and CZ analyzed the data. WL, CZ, and XW wrote the manuscript. All authors contributed to the article and approved the submitted version.
Funding
The project has been supported by the National Key Research and Development Program (2018YFD0201500) and National Natural Science Foundation of China (31871925 and 31970126).
Acknowledgments
We thank Dr. Hailei Wei at the Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences, for providing P. syringae DC3000ΔhopQ1-1, and the American Editor (www.mjeditor.com) for its linguistic assistance during the preparation of this manuscript.
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. The reviewer PY declared a shared affiliation with several of the authors, CZ, PD, XW, and WL, to the handling editor at the time of review.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.594669/full#supplementary-material
References
1
BosJ. I.ArmstrongM. R.GilroyE. M.BoevinkP. C.HeinI.TaylorR. M.et al (2010). Phytophthora infestans effector AVR3a is essential for virulence and manipulates plant immunity by stabilizing host E3 ligase CMPG1.Proc. Natl. Acad. Sci. U.S.A.1079909–9914. 10.1073/pnas.0914408107
2
BosJ. I.KannegantiT. D.YoungC.CakirC.HuitemaE.WinJ.et al (2006). The C-terminal half of Phytophthora infestans RXLR effector AVR3a is sufficient to trigger R3a-mediated hypersensitivity and suppress INF1-induced cell death in Nicotiana benthamiana.Plant J.48165–176. 10.1111/j.1365-313X.2006.02866.x
3
BovéJ. M. (2006). Huanglongbing: a destructive, newly-emerging, century-old disease of citrus.J. Plant Pathol.887–37. 10.4454/jpp.v88i1.828
4
BovéJ. M.AyresA. J. (2007). Etiology of three recent diseases of citrus in Sao Paulo state: sudden death, variegated chlorosis and Huanglongbing.IUBMB Life59346–354. 10.1080/15216540701299326
5
BricchiI.BerteaC. M.OcchipintiA.PaponovI. A.MaffeiM. E. (2012). Dynamics of membrane potential variation and gene expression induced by Spodoptera littoralis, Myzus persicae, and Pseudomonas syringae in Arabidopsis.PLoS One7:e46673. 10.1371/journal.pone.0046673
6
CarraS.AlbertiS.BeneschJ.BoelensW.BuchnerJ.CarverJ. A.et al (2019). Small heat shock proteins: multifaceted proteins with important implications for life.Cell Stress Chaperon.24295–308. 10.1007/s12192-019-00979-z
7
ChenP. Y.WangC. K.SoongS. C. (2003). Complete sequence of the binary vector pBI121 and its application in cloning T-DNA insertion from transgenic plants.Mol. Breed.11287–293. 10.1023/A:1023475710642
8
ChenY.LiuZ.HaltermanD. A. (2012). Molecular determinants of resistance activation and suppression by Phytophthora infestans effector IPI-O.PLoS Pathog.8:e1002595. 10.1371/journal.ppat.1002595
9
ClarkK.FrancoJ. Y.SchwizerS.PangZ.HawaraE.LiebrandT.et al (2018). An effector from the Huanglongbing-associated pathogen targets citrus proteases.Nat. Commun.9:1718. 10.1038/s41467-018-04140-9
10
CollN. S.EppleP.DanglJ. L. (2011). Programmed cell death in the plant immune system.Cell Death. Differ.181247–1256. 10.1038/cdd.2011.37
11
da CunhaL.SreerekhaM. V.MackeyD. (2007). Defense suppression by virulence effectors of bacterial phytopathogens.Curr. Opin. Plant Biol.10349–357. 10.1016/j.pbi.2007.04.018
12
da GraçaJ. V.DouhanG. W.HalbertS. E.KeremaneM. L.LeeR. F.VidalakisG.et al (2016). Huanglongbing: an overview of a complex pathosystem ravaging the world’s citrus.J. Integr. Plant Biol.58373–387. 10.1111/jipb.12437
13
DongS.YinW.KongG.YangX.QutobD.ChenQ.et al (2011). Phytophthora sojae avirulence effector Avr3b is a secreted NADH and ADP-ribose pyrophosphorylase that modulates plant immunity.PLoS Pathog.7:e1002353. 10.1371/journal.ppat.1002353
14
DouD.ZhouJ. M. (2012). Phytopathogen effectors subverting host immunity: different foes, similar battleground.Cell Host Microb.12484–495. 10.1016/j.chom.2012.09.003
15
DuanY.ZhouL.HallD. G.LiW.DoddapaneniH.LinH.et al (2009). Complete genome sequence of citrus Huanglongbing bacterium, “Candidatus Liberibacter asiaticus” obtained through metagenomics.Mol. Plant Microb. Interact.221011–1120. 10.1094/mpmi-22-8-1011
16
FangA.GaoH.ZhangN.ZhengX.QiuS.LiY.et al (2019). A novel effector gene SCRE2 contributes to full virulence of Ustilaginoidea virens to rice.Front. Microbiol.10:845. 10.3389/fmicb.2019.00845
17
FolimonovaS. Y.RobertsonC. J.GarnseyS. M.GowdaS.DawsonW. O. (2009). Examination of the responses of different genotypes of citrus to Huanglongbing (citrus greening) under different conditions.Phytopathology991346–1354. 10.1094/phyto-99-12-1346
18
GalloisP.MarinhoP. (1995). Leaf disk transformation using Agrobacterium tumefaciens-expression of heterologous genes in tobacco.Methods Mol. Biol.4939–48. 10.1385/0-89603-321-X:39
19
GottwaldT. R. (2010). Current epidemiological understanding of citrus Huanglongbing.Annu. Rev. Phytopathol.48119–139. 10.1146/annurev-phyto-073009-114418
20
GoverseA.SmantG. (2014). The activation and suppression of plant innate immunity by parasitic nematodes.Annu. Rev. Phytopathol.52243–265. 10.1146/annurev-phyto-102313-050118
21
GrabherrM. G.HaasB. J.YassourM.LevinJ. Z.ThompsonD. A.AmitI.et al (2011). Full-length transcriptome assembly from RNA-Seq data without a reference genome.Nat. Biotechnol.29644–652. 10.1038/nbt.1883
22
HaslbeckM.VierlingE. (2015). A first line of stress defense: small heat shock proteins and their function in protein homeostasis.J. Mol. Biol.4271537–1548. 10.1016/j.jmb.2015.02.002
23
HayesC. S.AokiS. K.LowD. A. (2010). Bacterial contact-dependent delivery systems.Annu. Rev. Genet.4471–90. 10.1146/annurev.genet.42.110807
24
HemetsbergerC.HerrbergerC.ZechmannB.HillmerM.DoehlemannG. (2012). The Ustilago maydis effector Pep1 suppresses plant immunity by inhibition of host peroxidase activity.PLoS Pathog.8:e1002684. 10.1371/journal.ppat.1002684
25
HuysmansM.LemaA. S.CollN. S.NowackM. K. (2017). Dying two deaths - programmed cell death regulation in development and disease.Curr. Opin. Plant Biol.3537–44. 10.1016/j.pbi.2016.11.005
26
JacobP.HirtH.BendahmaneA. (2017). The heat-shock protein/chaperone network and multiple stress resistance.Plant Biotechnol. J.15405–414. 10.1111/pbi.12659
27
JagoueixS.BovéJ. M.GarnierM. (1996). PCR detection of the two Candidatus Liberobacter species associated with greening disease of citrus.Mol. Cell Prob.1043–50. 10.1006/mcpr.1996.0006
28
JainM.FleitesL. A.GabrielD. W. (2015). Prophage-encoded peroxidase in ‘Candidatus Liberibacter asiaticus’ is a secreted effector that suppresses plant defenses.Mol. Plant Microb. Interact.281330–1337. 10.1094/mpmi-07-15-0145-r
29
JainM.Munoz-BodnarA.GabrielD. W. (2019). ‘Candidatus Liberibacter asiaticus’ peroxiredoxin (LasBCP) suppresses oxylipin-mediated defense signaling in citrus.J. Plant Physiol.23661–65. 10.1016/j.jplph.2019.03.001
30
JainM.Munoz-BodnarA.ZhangS.GabrielD. W. (2018). A secreted ‘Candidatus Liberibacter asiaticus’ peroxiredoxin simultaneously suppresses both localized and systemic innate immune responses in planta.Mol. Plant Microb. Interact.311312–1322. 10.1094/mpmi-03-18-0068-r
31
JonesL.HamiltonA. J.VoinnetO.ThomasC. L.MauleA. J.BaulcombeD. C. (1999). RNA-DNA interactions and DNA methylation in post-transcriptional gene silencing.Plant Cell112291–2301. 10.1105/tpc.11.12.2291
32
KamounS.van WestP.VleeshouwersV. G.de GrootK. E.GoversF. (1998). Resistance of Nicotiana benthamiana to Phytophthora infestans is mediated by the recognition of the elicitor protein INF1.Plant Cell101413–1426. 10.1105/tpc.10.9.1413
33
KawamuraY.HaseS.TakenakaS.KanayamaY.YoshiokaH.KamounS.et al (2009). INF1 elicitin activates jasmonic acid- and ethylene-mediated signalling pathways and induces resistance to bacterial wilt disease in tomato.J. Phytopathol.157287–297. 10.1111/j.1439-0434.2008.01489.x
34
LacommeC.Santa CruzS. (1999). Bax-induced cell death in tobacco is similar to the hypersensitive response.Proc. Natl. Acad. Sci U.S.A.967956–7961. 10.1073/pnas.96.14.7956
35
LiB.DeweyC. (2011). RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome.BMC Bioinform.12:323. 10.1186/1471-2105-12-323
36
LiH.YingX.ShangL.RedfernB.KypraiosN.XieX.et al (2020). Heterologous expression of CLIBASIA_03915/CLIBASIA_04250 by Tobacco mosaicvirus resulted in phloem necrosis in the senescent leaves of Nicotiana benthamiana.Int. J. Mol. Sci.21:1414. 10.3390/ijms21041414
37
LiJ.PangZ.TrivediP.ZhouX.YingX.JiaH.et al (2017). ‘Candidatus Liberibacter asiaticus’ encodes a functional salicylic acid (SA) hydroxylase that degrades SA to suppress plant defenses.Mol. Plant Microb. Interact.30620–630. 10.1094/MPMI-12-16-0257-R
38
LiZ.YinZ.FanY.XuM.KangZ.HuangL. (2015). Candidate effector proteins of the necrotrophic apple canker pathogen Valsa mali can suppress BAX-induced PCD.Front. Plant Sci.6:579. 10.3389/fpls.2015.00579
39
LiuX.FanY.ZhangC.DaiM.WangX.LiW. (2019). Nuclear import of a secreted “Candidatus Liberibacter asiaticus” protein is temperature dependent and contributes to pathogenicity in Nicotiana benthamiana.Front. Microbiol.10:1684. 10.3389/fmicb.2019.01684
40
LivakK. J.SchmittgenT. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2 (-Delta Delta C (T)) method.Methods25402–408. 10.1006/meth.2001.1262
41
Lo PrestiL.LanverD.SchweizerG.TanakaS.LiangL.TollotM.et al (2015). Fungal effectors and plant susceptibility.Annu. Rev. Plant Biol.66513–545. 10.1146/annurev-arplant-043014-114623
42
LocatoV.De GaraL. (2018). Programmed cell death in plants: an overview.Methods Mol. Biol.17431–8. 10.1007/978-1-4939-7668-3_1
43
MacLeanA. M.SugioA.MakarovaO. V.FindlayK. C.GrieveV. M.TóthR.et al (2011). Phytoplasma effector SAP54 induces indeterminate leaf-like flower development in Arabidopsis plants.Plant Physiol.157831–841. 10.1104/pp.111.181586
44
MaimboM.OhnishiK.HikichiY.YoshiokaH.KibaA. (2007). Induction of a small heat shock protein and its functional roles in Nicotiana plants in the defense response against Ralstonia solanacearum.Plant Physiol.1451588–1599. 10.1104/pp.107.105353
45
MartinelliF.DandekarA. M. (2017). Genetic mechanisms of the devious intruder Candidatus Liberibacter in Citrus.Front. Plant Sci.8:904. 10.3389/fpls.2017.00904
46
MdluliK. E.TreitJ. D.KerrV. J.NanoF. E. (1995). New vectors for the in vitro generation of alkaline phosphatase fusions to proteins encoded by G+C-rich DNA.Gene155133–134. 10.1016/0378-1119(94)00909-c
47
MortazaviA.WilliamsB. A.McCueK.SchaefferL.WoldB. (2008). Mapping and quantifying mammalian transcriptomes by RNA-Seq.Nat. Methods5621–628. 10.1038/nmeth.1226
48
NelsonB. K.CaiX.NebenführA. (2007). A multicolored set of in vivo organelle markers for co-localization studies in Arabidopsis and other plants.Plant J.511126–1136. 10.1111/j.1365-313X.2007.03212.x
49
PangZ.ZhangL.CoakerG.MaW.HeS.WangN. (2020). Citrus CsACD2 is a target of Candidatus Liberibacter asiaticus in Huanglongbing disease.Plant Physiol.184792–805. 10.1104/pp.20.00348
50
PetrovV. D.Van BreusegemF. (2012). Hydrogen peroxide-a central hub for information flow in plant cells.AoB Plants2012:pls014. 10.1093/aobpla/pls014
51
PitinoM.AllenV.DuanY. (2018). LasΔ5315 Effector induces extreme starch accumulation and chlorosis as Ca. Liberibacter asiaticus infection in Nicotiana benthamiana.Front. Plant Sci.9:113. 10.3389/fpls.2018.00113
52
PitinoM.ArmstrongC. M.CanoL. M.DuanY. (2016). Transient expression of Candidatus Liberibacter asiaticus effector induces cell death in Nicotiana benthamiana.Front. Plant Sci.7:982. 10.3389/fpls.2016.00982
53
PrasadS.XuJ.ZhangY.WangN. (2016). SEC-translocon dependent extracytoplasmic proteins of Candidatus Liberibacter asiaticus.Front. Microbiol.7:1989. 10.3389/fmicb.2016.01989
54
RajputN. A.ZhangM.ShenD.LiuT.ZhangQ. (2015). Overexpression of a Phytophthora cytoplasmic CRN effector confers resistance to disease, salinity and drought in Nicotiana benthamiana.Plant Cell Physiol.562423–2435. 10.1093/pcp/pcv164
55
ShiQ.PitinoM.ZhangS.KrystelJ.CanoL. M.ShattersR. G.et al (2019). Temporal and spatial detection of Candidatus Liberibacter asiaticus putative effector transcripts during interaction with Huanglongbing-susceptible, -tolerant, and -resistant citrus hosts.BMC Plant Biol.19:122. 10.1186/s12870-019-1703-4
56
StaelS.KmiecikP.WillemsP.Van Der KelenK.CollN. S.TeigeM.et al (2015). Plant innate immunity–sunny side up?Trends Plant Sci.203–11. 10.1016/j.tplants.2014.10.002
57
SugioA.KingdomH. N.MacLeanA. M.GrieveV. M.HogenhoutS. A. (2011). Phytoplasma protein effector SAP11 enhances insect vector reproduction by manipulating plant development and defense hormone biosynthesis.Proc. Natl. Acad. Sci. U.S.A.108E1254–E1263. 10.1073/pnas.1105664108
58
Ul HaqS.KhanA.AliM.KhattakA. M.GaiW. X.ZhangH. X.et al (2019). Heat shock proteins: dynamic biomolecules to counter plant biotic and abiotic stresses.Int. J. Mol. Sci.20:5321. 10.3390/ijms20215321
59
Van der HoornR. A.LaurentF.RothR.De WitP. J. (2000). Agroinfiltration is a versatile tool that facilitates comparative analyses of Avr9/Cf-9-induced and Avr4/Cf-4-induced necrosis.Mol. Plant Microb. Interact.13439–446. 10.1094/mpmi.2000.13.4.439
60
Van OoijenG.LukasikE.Van Den BurgH. A.VossenJ. H.CornelissenB. J.TakkenF. L. (2010). The small heat shock protein 20 RSI2 interacts with and is required for stability and function of tomato resistance protein I-2.Plant J.63563–572. 10.1111/j.1365-313X.2010.04260.x
61
VanackerH.CarverT. L.FoyerC. H. (2000). Early H2O2 accumulation in mesophyll cells leads to induction of glutathione during the hyper-sensitive response in the barley-powdery mildew interaction.Plant Physiol.1231289–1300. 10.1104/pp.123.4.1289
62
WangC.LiC.DuanG.WangY.ZhangY.YangJ. (2019). Overexpression of Magnaporthe oryzae systemic defense trigger 1 (MoSDT1) confers improved rice blast resistance in Rice.Int. J. Mol. Sci.20:4762. 10.3390/ijms20194762
63
WangL.FengZ.WangX.WangX.ZhangX. (2010). DEGseq: an R package for identifying differentially expressed genes from RNA-seq data.Bioinformatics26136–138. 10.1093/bioinformatics/btp612
64
WangN.TrivediP. (2013). Citrus Huanglongbing: a newly relevant disease presents unprecedented challenges.Phytopathology103652–665. 10.1094/phyto-12-12-0331-rvw
65
WangN.PiersonE. A.SetubalJ. C.XuJ.LevyJ. G.ZhangY.et al (2017). The Candidatus Liberibacter-host interface: insights into pathogenesis mechanisms and disease control.Annu. Rev. Phytopathol.55451–482. 10.1146/annurev-phyto-080516-35513
66
WangQ.HanC.FerreiraA. O.YuX.YeW.TripathyS.et al (2011). Transcriptional programming and functional interactions within the Phytophthora sojae RXLR effector repertoire.Plant Cell232064–2086. 10.1105/tpc.111.086082
67
WarnesG. R.BolkerB.BonebakkerL.GentlemanR.HuberW.LiawA.et al (2015). gplots: Various R Programming Tools for Plotting Data. Seattle, WA, USA. Available online: https://www.scienceopen.com/document?vid=0e5d8e31-1fe4-492f-a3d8-8cd71b2b8ad9(accessed September 16, 2020).
68
WatersE. R.VierlingE. (1999). The diversification of plant cytosolic small heat shock proteins preceded the divergence of mosses.Mol. Biol. Evol.16127–139. 10.1093/oxfordjournals.molbev.a0
69
WeiC. F.KvitkoB. H.ShimizuR.CrabillE.AlfanoJ. R.LinN. C.et al (2007). A Pseudomonas syringae pv. tomato DC3000 mutant lacking the type III effector HopQ1-1 is able to cause disease in the model plant Nicotiana benthamiana.Plant J.5132–46. 10.1111/j.1365-313X.2007.03126.x
70
WeiH. L.ChakravarthyS.WorleyJ. N.CollmerA. (2013). Consequences of flagellin export through the type III secretion system of Pseudomonas syringae reveal a major difference in the innate immune systems of mammals and the model plant Nicotiana benthamiana.Cell Microbiol.15601–618. 10.1111/cmi.12059
71
WeiH. L.ZhangW.CollmerA. (2018). Modular study of the Type III effector repertoire in Pseudomonas syringae pv. tomato DC3000 reveals a matrix of effector interplay in pathogenesis.Cell Rep.231630–1638. 10.1016/j.celrep.2018.04.037
72
WiltonM.SubramaniamR.ElmoreJ.FelsensteinerC.CoakerG.DesveauxD. (2010). The type III effector HopF2Pto targets Arabidopsis RIN4 protein to promote Pseudomonas syringae virulence.Proc. Natl. Acad. Sci. U.S.A.1072349–2354. 10.1073/pnas.0904739107
73
XiangJ.LiX.WuJ.YinL.ZhangY.LuJ. (2016). Studying the mechanism of Plasmopara viticola RxLR effectors on suppressing plant immunity.Front. Microbiol.7:709. 10.3389/fmicb.2016.00709
74
YamaneK.BunaiK.KakeshitaH. (2004). Protein traffic for secretion and related machinery of Bacillus subtilis.Biosci. Biotechnol. Biochem.682007–2023. 10.1271/bbb.68.2007
75
YingX.WanM.HuL.ZhangJ.LiH.LvD. (2019). Identification of the virulence factors of Candidatus Liberibacter asiaticus via heterologous expression in Nicotiana benthamiana using tobacco mosaic virus.Int. J. Mol. Sci.20:5575. 10.3390/ijms20225575
76
ZhangC.WangX.LiuX.FanY.ZhangY.ZhouX.et al (2019). A novel ‘Candidatus Liberibacter asiaticus’-encoded Sec-dependent secretory protein suppresses programmed cell death in Nicotiana benthamiana.Int. J. Mol. Sci.20:5802. 10.3390/ijms20225802
77
ZhangM.LiQ.LiuT.ShenD.ZhuY.LiuP.et al (2015a). Two cytoplasmic effectors of Phytophthora sojae regulate plant cell death via interactions with plant catalases.Plant Physiol.167164–175. 10.1104/pp.114.252437
78
ZhangM.RajputN. A.ShenD.ZengW.LiuT. (2015b). A Phytophthora sojae cytoplasmic effector mediates disease resistance and abiotic stress tolerance in Nicotiana benthamiana.Sci. Rep.5:10837. 10.1038/srep10837
79
ZhengZ.ChenJ.DengX. (2018). Historical perspectives, management, and current research of citrus HLB in Guangdong province of China, where the disease has been endemic for over a hundred years.Phytopathology1081224–1236. 10.1094/phyto-07-18-0255-ia
Summary
Keywords
Candidatus Liberibacter asiaticus, sec-dependent secretory protein, hypersensitive cell death, small heat shock protein, virulence factor
Citation
Zhang C, Du P, Yan H, Zhu Z, Wang X and Li W (2020) A Sec-Dependent Secretory Protein of the Huanglongbing-Associated Pathogen Suppresses Hypersensitive Cell Death in Nicotiana benthamiana. Front. Microbiol. 11:594669. doi: 10.3389/fmicb.2020.594669
Received
14 August 2020
Accepted
04 November 2020
Published
30 November 2020
Volume
11 - 2020
Edited by
Carmen R. Beuzón, University of Malaga, Spain
Reviewed by
Norma Elena Leyva-López, National Polytechnic Institute (IPN), Mexico; Ping Yang, Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, China
Updates
Copyright
© 2020 Zhang, Du, Yan, Zhu, Wang and Li.
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: Xuefeng Wang, wangxuefeng@cric.cnWeimin Li, liweimin01@caas.cn
This article was submitted to Microbe and Virus Interactions with Plants, a section of the journal Frontiers in Microbiology
Disclaimer
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.