Viral Strain-Specific Activation of Pathogen-Associated Molecular Pattern-Triggered Immunity Enhances Symptom Severity in Broad Bean Wilt Virus 2 Infection

Broad bean wilt virus 2 (BBWV2) is an emerging virus in various economically important crops, especially pepper (Capsicum annuum L.), worldwide. Recently, the emergence of various BBWV2 strains that induce severe symptoms has increased damage to pepper crops. While the symptomatic variations among virus strains should be associated with differences in the transcriptomic reprogramming of host plants upon infection, underlying molecular mechanisms and associated genes are largely unknown. In the present study, we employed transcriptome analysis to identify responsible host factors for symptom enhancement in the BBWV2-pepper pathosystem using two distinct BBWV2 strains, PAP1 (a severe strain) and RP1 (a mild strain). Comparative analysis of the differentially expressed genes (DEGs) revealed that various genes associated with pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI) and ethylene signaling were significantly upregulated upon infection with the severe PAP1 strain, but not with the mild RP1 strain. Indeed, hormone analysis revealed that ethylene emission was significantly increased in pepper plants infected with PAP1. These observations imply that the activation of the PTI-associated defense responses reinforce symptom formation during BBWV2 infection in a virus strain-specific manner.


INTRODUCTION
The symptom severity of crop diseases is closely related to productivity (Gaunt, 1995). Plant viruses are obligate intracellular parasites that disturb the normal activity of host cells, resulting in a wide range of symptoms. Different strains of plant viruses often induce symptoms with different severity in the host plants. Such symptomatic variation can be triggered by molecular interactions between the host and viral strain-specific factors, leading to transcriptome reprogramming of numerous genes involved not only in defense responses, but also in plant development, physiology, and metabolism (Seo et al., 2018). In this regard, the identification of crucial host genes associated with symptomatic variation and their underlying molecular mechanisms can provide new insights into molecular plant physiological phenomena as well as symptom development.
Pepper (Capsicum annuum L.) is an economically important vegetable crop worldwide. However, commercial pepper production in many global areas continues to suffer from various viral diseases. While more than 100 viruses infect pepper plants worldwide (Kenyon et al., 2014), broad bean wilt virus 2 (BBWV2; genus Fabavirus, family Secoviridae) is an emerging virus in many global areas (Kwak et al., 2013a,b;Ferriol et al., 2014). The incidence of BBWV2 in pepper has increased rapidly over the past decade, because the virus can be easily transmitted by aphids and infect a wide range of plants, including various weed hosts that inhabit areas near pepper fields (Kwak et al., 2013a;Kim et al., 2014a). The BBWV2 genome is composed of two segments of positivesense single-stranded RNAs, RNA1 (~5.8 knt) and RNA2 (~3.3 knt; Ferrer et al., 2011). Each RNA segment encodes a single open reading frame (ORF) that translates into a large polyprotein precursor. Proteolysis of the polyprotein precursor encoded by RNA1 yields five mature proteins: protease cofactor (Co-Pro), NTP-binding motif (NTBM), VPg, protease (Pro), and RNA-dependent RNA polymerase (RdRp). After the same processing, RNA2 expresses the following three mature proteins: movement protein (MP), large coat protein (LCP), and small coat protein (SCP; Ferrer et al., 2011).
Broad bean wilt virus 2 causes diverse symptoms in pepper, depending on the compatibility between pepper varieties and virus strains, including mild mosaic, malformation, stunting, and chlorosis (Kwak et al., 2013b(Kwak et al., , 2016. In a previous study, we generated infectious cDNA clones of two pathogenically different BBWV2 strains, namely pBBWV2-RP1 (a mild strain) and pBBWV2-PAP1 (a severe strain; Kwak et al., 2016). Using the chimeric viruses and amino acid substitution mutant viruses derived from pBBWV2-RP1 and pBBWV2-PAP1, we demonstrated that MP is the viral strain-specific factor that determines symptom severity in BBWV2 . In the present study, we employed comparative transcriptome analysis to gain a better understanding of the molecular mechanisms associated with the development of distinct symptoms caused by the two BBWV2 strains in pepper.

Plant Growth and Virus Inoculation
Pepper (Capsicum annuum L. cv. Sinhong) and Nicotiana benthamiana plants were grown in an insect-free growth chamber in a cycle of 16 h light at 26°C and 8 h darkness at 24°C. Full-length infectious cDNA clones of two BBWV2 strains (RP1 and PAP1), as described in our previous study (Kwak et al., 2016), were used as the viral sources for each strain. Infectious cDNA clones of the virus strains were inoculated into the leaves of 2-week-old N. benthamiana plants by Agrobacterium-mediated infiltration (agroinfiltration), as previously described (Seo et al., 2009 To analyze virus accumulation levels in the systemically infected leaves, the same RNA preparations used for the RNA sequencing were subjected to quantitative realtime RT-PCR (qRT-PCR) using an iCycler iQ5 qRT-PCR detection system (Bio-Rad, United States) with the following specific primers: BBWV2-R1-RT-Fw (5'-TCACAGGTTATGCCG CTTGT-3') and BBWV2-R1-RT-Rv (5'-TCACTCGTCCCAAGC TGTTC-3') for BBWV2 RNA1 detection and UBI2-F (5'-TACCCTTCACCTTGTCCTCC-3') and UBI2-R (5'-GCCAT CCTCCAACTGTTTTC-3') for ubiquitin2 mRNA detection. The ubiquitin2 gene (CA.PGAv.1.6.scaffold337.91) was used as an internal reference to standardize the different samples. Three biological and technical replicates were used per sample.

Processing of mRNA Sequence Data
The Illumina pipeline filtrated raw sequence reads, which were then mapped to the reference transcripts of C. annuum cv. CM334 (Kim et al., 2014b) retrieved from the Pepper Genome Platform 1 using the RNA-seq mapping algorithm implemented in bowtie2 (v2.1.0) software (Langmead et al., 2009), allowing all aligning with a maximum of two mismatches. The raw data were deposited in the National Center for Biotechnology Information (NCBI) SRA database with BioProject accession number PRJNA751625. To eliminate bias caused by variations in sequencing depth, the number of mapped clean reads for each genes was measured and then normalized with the DESeq package in R software (Anders and Huber, 2010). Differentially expressed genes (DEGs) were identified by a ≤2-fold change in read coverage and a binomial test with a false discovery rate (FDR) ≤ 0.01. The FDR was applied to identify the threshold value of p for multiple tests and was calculated using DESeq. Correlation analysis and hierarchical clustering were performed to categorize the genes according to patterns of expression using the AMAP library in R (Lucas, 2014).

Gene Enrichment Analysis
Gene ontology (GO) analysis was performed to functionally annotate the DEGs based on the protein sequence similarity (e-value cutoff ≤ 1e-10) in the GO database (Ashburner et al., 2000). The number of DEGs assigned in each GO term was counted using the in-house scripts of SEEDERS Inc (Daejeon, South Korea). GO term enrichment was performed using the PANTHER overrepresentation test (Mi et al., 2013). Functional enrichment analysis was performed to assign biological relevance to the gene network modules using agriGO v2.0 (Tian et al., 2017). The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway was analyzed using the sequence similarity (e-value cut off ≤ 1e-10, identity ≥90) of pepper proteins in the KEGG database (Kanehisa and Goto, 2000).

DAB Staining
3,3'-diaminobenzidine (DAB) staining was performed to detect reactive oxygen species (ROS) as previously described . Briefly, upper uninoculated leaf samples were collected from the virus-infected and healthy plants at 14 days post-inoculation (dpi). The leaf samples were stained with 1 mg/ml DAB (Sigma, United States) by vacuum infiltration and incubated for 8 h at room temperature in the dark. The leaf samples were then destained with 90% ethanol. The destained samples were mounted in 50% glycerol and observed by light microscopy.

Ethylene Analysis
Ethylene production from virus-infected and healthy plants was measured at 14 dpi as described previously with minor modifications (Hou et al., 2018). Briefly, 5 g of the upper uninoculated leaves were sampled from three plants and placed in 400 ml glass jars for 5 h at 26°C in the dark. Next, 1 ml of the jar headspace gas sample was sampled and analyzed using a gas chromatograph equipped with a hydrogen flame ionization detector and an activated alumina column (Agilent Technologies, United States). Samples were compared to a standard of known concentration. The experiment was repeated three times, and data were averaged.

Transcriptome Analysis of Symptomatic Variation Caused by Two Distinct BBWV2 Strains in Pepper
Two distinct BBWV2 strains, PAP1 and RP1, cause symptoms with different severity in pepper and N. benthamiana: BBWV2-PAP1 induced severe symptoms of mosaic, leaf malformation, and stunting in pepper and N. benthamiana, whereas BBWV2-RP1 caused no visible symptoms in pepper and very mild symptoms in N. benthamiana (Figure 1).
We sought to examine host genes associated with viral strain-specific symptoms in the BBWV2-pepper pathosystem. To this end, Illumina RNA sequencing was performed to examine the transcriptomic changes in pepper upon infection with BBWV2-PAP1 and RP1. We first analyzed viral RNA accumulation of each strain in the systemic leaves using qRT-PCR. The results showed no significant difference in viral RNA accumulation between the two BBWV2 strains (Figure 2A), suggesting that the symptomatic difference between the two BBWV2 strains was not due to differences in virus accumulation levels. The same RNA preparation was subjected to cDNA library construction. Nine cDNA libraries (three for each sample) were sequenced by Illumina RNA sequencing. For mapping, we used a total of 35,884 reference transcripts of C. annuum cv. CM334 (Kim et al., 2014b). The Illumina pipeline filtered and trimmed the raw HiSeq reads, yielding approximately 81-86 million clean pair-end reads from each of the nine libraries ( Table 1). The acquired reads were mapped on the reference pepper transcripts; this resulted in the mapping of approximately 83.6-92.81% of the nucleotides ( Table 1). The nucleotide coverage for healthy, BBWV2-RP1, and BBWV2-PAP1 was 175, 164.68, and 168.45 times, respectively ( Table 1).

Identification of DEGs in Response to Infection With BBWV2-PAP1 or RP1
The fragments per kilobase of exon per million fragments mapped (FPKM) value was used to normalize the expression levels of the mapped genes. DEGs were identified by statistically comparing the FPKM values of two samples using a 2-fold change in expression with the FDR ≤ 0.01. By comparing virusinfected samples to healthy samples, we identified 160 and 22 DEGs in response to infection with BBWV2-PAP1 and RP1, respectively ( Figure 2B; Supplementary Table S2). When the DEGs were compared between the two infection conditions, a total of 147 genes were identified to be specifically expressed in response to BBWV2-PAP1 (upregulation of 141 genes and downregulation of six genes; Figure 2C; Table 2;  Supplementary Table S2). Comparatively few genes were specifically affected upon infection with BBWV2-RP1 (upregulation of four genes and downregulation of five genes; Figure 2C; Supplementary Table S2). Hierarchical clustering of the total DEGs revealed six distinct groups based on the expression patterns in each infection condition ( Figure 2D).

Gene Ontology Analysis of Identified DEGs
For a better understanding of the DEGs associated with severe symptom development caused by BBWV2-PAP1, GO analysis was performed using the DEGs upregulated upon BBWV2-PAP1 infection. A total of 37 GO terms were significantly enriched by infection with BBWV2-PAP1 (Supplementary Table S3). Most BBWV2-PAP1-specific GO terms were associated with defense response and response to stimulus, e.g., systemic acquired resistance (GO:0009627), response to ethylene (GO:0009723), defense response to fungus (GO:0050832), response to salicylic acid (GO:0009751), innate immune response (GO:0045087), response to reactive oxygen species (GO:0000302), and immune response (GO:0006955). Hierarchical GO enrichment analysis using agriGO v2.0 further showed that the identified GO terms were highly correlated in a network context (Figure 3). The analysis also highlighted a few downstream GO terms, including defense response to fungus (GO:0050832), response to hormone (GO:0009725), and response to ethylene (GO:0009723), which are likely directly correlated to symptom development in response to BBWV2-PAP1 (Figure 3).

Important KEGG Pathways Influenced by BBWV2-PAP1 Infection
Kyoto encyclopedia of genes and genomes pathway analysis was performed to understand the relationship between DEGs and the associated pathways. KEGG analysis using the DEGs upregulated in response to BBWV2-PAP1 infection identified some of the highly ranked pathways, including plant-pathogen interaction, metabolic pathways, plant hormone signal transduction, and pathways related to plant pigment metabolism (i.e., flavonoid biosynthesis, flavone and flavonol biosynthesis, and biosynthesis of secondary metabolites; Table 3). In particular, the most significant enrichment of upregulated DEGs was found in the KEGG pathway of plant-pathogen interaction (Pathway ID: cann04626; Table 3; Figure 4A). Various receptor-like protein kinase genes (Figure 4B), pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI) signaling genes (Figure 4C), and pathogenesis-related protein genes ( Figure 4D) involved in the KEGG pathway of plant-pathogen interaction were significantly upregulated in response to infection with BBWV2-PAP1. KEGG analysis also identified a few ethylene responsive transcription factor genes, including ERF5, in the plant hormone signal transduction pathway (Pathway ID: cann00943; Figure 4E). Our DEG analysis also revealed that key genes involved in ethylene biosynthesis (i.e., ACC oxidase genes) were significantly upregulated in response to infection  with BBWV2-PAP1 ( Figure 4E). Overall, our results revealed that the PTI and ethylene pathways were activated in BBWV2-PAP1-infectged pepper.

BBWV2-PAP1 Infection Increases ROS Accumulation and Ethylene Emission in Pepper
Pathogen-associated molecular pattern-triggered immunity response typically accompanies the accumulation of ROS and is sometimes associated with ethylene pathways (Peng et al., 2018). Because our transcriptomic analysis revealed that BBWV2-PAP1 infection activates PTI and ethylene pathways in pepper, we investigated whether ROS accumulation and ethylene production are enhanced in BBWV2-PAP1-infected leaf tissues.
To detect ROS accumulation, DAB staining was performed on systemic leaves of the pepper plants infected with BBWV2-PAP1 or RP1. The leaves infected with BBWV2-PAP1 showed more intense staining than BBWV2-RP1-infected leaves (Figure 5A), indicating that ROS accumulation was specifically elevated by infection with BBWV2-PAP1 but not RP1.
To verify whether activation of the ethylene pathway leads to an increase in ethylene production in BBWV2-PAP1-infected pepper plants, ethylene emission was assessed for detached leaves infected with BBWV2-PAP1 or RP1. There was no significant difference in ethylene production in healthy and BBWV2-RP1-infected leaves ( Figure 5B); however, ethylene production was considerably higher in BBWV2-PAP1-infected leaves than in BBWV2-RP1-infected leaves ( Figure 5B). This implies that, consistent with the transcriptional upregulation of the ethylene pathway genes, BBWV2-PAP1 infection causes an increase in ethylene production in pepper.

Overexpression of Ethylene Responsive Factor 5 Enhances Symptom Severity During BBWV2-RP1 Infection
Ethylene biosynthesis is strongly induced in PTI and plays an important role in defense responses (Yang et al., 2017;Peng et al., 2018). ERF5 functions as a positive regulators of ethylenemediated immunity in Arabidopsis (Moffat et al., 2012;Son et al., 2012). Thus, we hypothesized that if the activation of ethylene signaling pathways is correlated with symptom severity in pepper infected with BBWV2, the ectopic overexpression of ERF5 would enhance symptom severity in pepper infected with BBWV2-RP1 (a mild strain). Testing this hypothesis required the overexpression of ERF5 in the pepper cells infected with BBWV2-RP1. For this, we utilized a viral vector system generated by engineering BBWV2-RP1 to overexpress foreign genes in pepper (Choi et al., 2019). Specifically, we engineered the BBWV2-RP1-based vector to overexpress ERF5 in the BBWV2-RP1-infected cells, and the resulting construct was named pBBWV2-RP1-R2-ERF5 ( Figure 6A). To examine if the ERF5 overexpression affects symptom severity caused by BBWV2, pepper and N. benthamiana plants were infiltrated with a mixture of Agrobacterium cultures containing pBBWV2-RP1-R1 and pBBWV2-RP1-R2-ERF5. Interestingly, both plant species were systemically infected with the recombinant BBWV2 vector overexpressing ERF5 and developed severe symptoms of leaf malformation and stunting ( Figure 6B). In particular, necrotic cell death responses were also observed in the systemically infected pepper leaves ( Figure 6B). Therefore, our results revealed that the ethylene pathway-mediated responses specifically induced by BBWV2-PAP1 are associated with enhanced symptom severity in BBWV2 infection.

Quantitative Real-Time PCR Validation of RNA Sequencing Data
To validate the gene expression data from the RNA sequencing, the expression levels of 12 representative genes that significantly altered in response to virus infection were evaluated by qRT-PCR using the same RNA preparations used for RNA sequencing. The expression of target genes was normalized to that of the ubiquitin2 gene (CA.PGAv.1.6.scaffold337.91) as an internal control. For each gene, relative expression levels were calculated and compared to the RNA sequencing results. In general, the qRT-PCR results for all of the tested genes were consistent with the RNA sequencing results (Supplementary Table S4).   This demonstrates that the changes in gene expression determined by the RNA sequencing were accurate.

DISCUSSION
Broad bean wilt virus 2 is an emerging virus in various economically important crops worldwide. BBWV2 has a wide host range, and the viral strains isolated from different hosts are genetically and pathogenically differentiated (Ferrer et al., 2011;Kwak et al., 2013aKwak et al., , 2016. Symptomatic variation among different virus strains in a host plant is primarily due to strainspecific viral proteins. Strain-specific viral proteins can have different activities to their original functions or interact with other proteins. In BBWV2, MP was identified as a strain-specific elicitor that determine symptom severity in pepper and N. benthamiana (Kwak et al., 2016;Seo et al., 2017); however, the molecular mechanisms and host genes involved in symptom development upon BBWV2 infection remain unknown. Different symptoms induced by two different strains of BBWV2 (PAP1 and RP1) indicate that transcriptional responses must be different. Notably, BBWV2-RP1 induced no visible symptoms in pepper (Figure 1) while its replication level was similar to BBWV2-PAP1 (Figure 2A). Indeed, only a small number of DEGs were identified in pepper infected with BBWV2-RP1 (Figures 2B,C). On the other hand, BBWV2-PAP1 infection resulted in the upregulation of 150 DEGs (Figure 2B). Our GO term and KEGG pathway analyses revealed that many of the upregulated DEGs were PTI-associated genes, including receptor-like protein kinase genes, PTI signaling genes, and pathogenesis-related protein genes (Figures 3, 4).
Plant cells can recognize microbial invasions through the perception of PAMPs by pattern-recognition receptors (PRRs) localized on the plasma membrane, triggering the PTI activation (Peng et al., 2018). PTI involves a series of downstream defense responses, such as mitogen-activated protein kinase activation, ROS accumulation, plant hormone biosynthesis, and expression of pathogenesis-related proteins, resulting in the prevention of pathogen colonization (Peng et al., 2018). Our comparative transcriptome analysis revealed that PTI pathways can be activated in a virus strain-specific manner in pepper infected with BBWV2 (Figure 4). Many receptor-like protein kinase genes were specifically upregulated in response to BBWV2-PAP1 infection ( Figure 4B). In plants, most surface PRRs are receptorlike protein kinases in plants (Wang et al., 2016). The Arabidopsis leucine-rich repeat receptor kinases FLAGELLIN SENSING 2 and EF-TU RECEPTOR, which can recognize bacterial flagellin and EF-Tu, respectively, are well-characterized PRRs (Zipfel, 2009;Tang et al., 2017). PAMP perception upregulates a number of receptor-like protein kinases, indicating that PTI increases the capacity for pathogen recognition (Zipfel et al., 2004(Zipfel et al., , 2006. Broad bean wilt virus 2-PAP1 infection also resulted in the upregulation of various PTI signaling genes, including mitogenactivated protein kinase (CA.PGAv.1.6.scaffold200.11), WRKY transcription factor 70 (CA.PGAv.1.6.scaffold788.4), and WRKY transcription factor 70-like (CA.PGAv.1.6.scaffold1018.1; Figure 4C), which can lead to changes in the transcription of numerous genes involved in defense responses (Wang et al., 2016;Peng et al., 2018). Various pathogenesis-related protein genes were also upregulated in BBWV2-PAP1-infected pepper ( Figure 4D). The basic pathogenesis-related protein 1 gene (CA.PGAv.1.6.scaffold608.26) was highly upregulated in an incompatible interaction with Xanthomonas campestris pv. vesicatoria in pepper (Kim and Hwang, 2000). Interestingly, it was demonstrated that the upregulation of the basic pathogenesisrelated protein 1 gene was intrinsically associated with increased ethylene production (Kim and Hwang, 2000). In pepper, treatment with aminoethoxyvinylglycine (AVG), an inhibitor of ethylene biosynthesis, resulted in the suppression of basic pathogenesisrelated protein 1 gene expression, indicating that ethylene is an upstream regulator of this gene (Kim and Hwang, 2000). In addition, other ethylene-responsive pathogenesis-related protein genes, i.e., pathogenesis-related protein 2 and pathogenesis-related protein 3, which encode β-1,3-glucanases and chitinase, respectively (van Loon et al., 2006b;Mazarei et al., 2007), were specifically upregulated in response to infection with BBWV2-PAP1 ( Figure 4D). Therefore, our results suggest that PTI responses triggered by BBWV2-PAP1 infection are associated with ethylene signaling in pepper. Indeed, ethylene production and the pathway genes were significantly upregulated in pepper infected with BBWV2-PAP1 (Figures 4E, 5B).
Ethylene is one of the important regulating factors in disease responses as well as growth and development in plants (van Loon et al., 2006a;Huang et al., 2016). Accumulating evidence suggests that the involvement of ethylene signaling in defense responses is linked to the activation of ROS production (Mersmann et al., 2010;Yang et al., 2017). Similarly, BBWV2-PAP1 infection elevated ROS accumulation as well as ethylene production in pepper ( Figure 5). In defense responses, ROS can function as molecular signals to activate programmed cell death, thereby increasing symptom severity (Peng et al., 2018). Ethylene signaling has also been shown to positively regulate the FLS2 receptor accumulation in Arabidopsis (Mersmann et al., 2010), which is consistent with our observation that various receptor-like protein kinase genes were upregulated by BBWV2-PAP1 infection ( Figure 4A). Ethylene signaling involves plant-specific transcription factors, such as ERFs, which regulate various responses to environmental stimuli (Muller and Munne-Bosch, 2015). While various ERFs have been well characterized recently, ERF5 was shown to activate ethylene-dependent defense genes, and its overexpression resulted in increased resistance to the necrotrophic pathogen Botrytis cinerea (Moffat et al., 2012). However, no visible phenotypic alterations were observed in transgenic tobacco plants that overexpress ERF5 (Fischer and Dröge-Laser, 2004).
In this study, we demonstrated that overexpression of ERF5 can strongly enhance symptom severity in pepper even when infected with the mild BBWV2 RP1 strain (Figure 6). This further demonstrates that the ethylene signaling pathway is associated with the enhancement of symptom severity in BBWV2 infection.
A previous study showed that tobacco mosaic virus coat protein can be recognized as a PAMP to activate PTI and oxidative burst (Allan et al., 2001). Furthermore, for BBWV2, viral strain-specific MP was identified as a symptom severity determinant . In this study, we found that the viral strain-specific activation of PTI was associated with an increase in symptom severity in BBWV2 infection, that is, BBWV2 MP acts as a PAMP in a virus strainspecific manner. In this regard, it is interesting to consider how plant virus MP activates PTI because, in general, PTI is A B FIGURE 6 | Overexpression of ethylene responsive factor 5 (ERF5) enhances symptom severity in pepper infected with the BBWV2 mild strain RP1. (A) Schematic representation of the BBWV2-RP1 RNA2 recombinant constructs. The ERF5 gene was inserted between movement protein (MP) and large coat protein (LCP) to be expressed by proteolysis when the virus replicates. In vivo transcription of the BBWV2-RP1 RNA2 recombinant constructs is under the control of a double 35S promoter (P35SX2), a cis-cleaving ribozyme sequence (Rz), and a NOS terminator (NOSt). (B) Symptomatic responses in pepper and N. benthamiana plants upon infection with pBBWV2-RP1-ERF5. The inocula are indicated above each image. An enlargement of the boxed area is shown in the upper right panel.
activated upon the extracellular recognition of PAMPs by PRRs. BBWV2 MP is localized to plasmodesmata (PD), and forms tubule structures to facilitate virus cell-to-cell movement (Xie et al., 2016). PD are plasma membrane-lined channels that provide symplastic continuity between neighboring cells (Maule et al., 2011). PD membranes also contain various receptor proteins and receptor-like protein kinases (Faulkner, 2013). A previous study revealed that the interactions between viral MP and PD-located receptor-like proteins are crucial for tubule formation and virus cell-to-cell movement (Amari et al., 2010). Therefore, it is possible that BBWV2 MP may interact with PD-located PRRs, resulting in PTI signaling activation.
Disease symptom development is a complex physiological process involving large transcriptomic changes. Thus, analysis of the DEGs in plant-virus interactions is important to explore the molecular basis of various physiological phenomena in plants. The transcriptional responses in pepper were notably different after infection with the two different BBWV2 strains (PAP1 and RP1), which had different symptoms. Our comparative transcriptome analysis revealed molecular genetic clues for the increase in symptom severity upon infection with BBWV2-PAP1. Overall, our findings improve understanding of molecular mechanisms underlying disease symptom development in plants and provide a basis for the future exploration of the functions of pepper genes in fundamental plant physiology.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are publicly available. This data can be found here: National Center for Biotechnology Information (NCBI) BioProject database under accession number PRJNA751625.