Next Article in Journal
Klebsiella pneumoniae Phage M198 and Its Therapeutic Potential
Previous Article in Journal
Expanding the Diversity of Actinobacterial Tectiviridae: A Novel Genus from Microbacterium
Previous Article in Special Issue
A Survey of Wild Indigenous Cryptostylis ovata Orchid Populations in Western Australia Reveals Spillover of Exotic Viruses
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genomic, Evolutionary, and Pathogenic Characterization of a New Polerovirus in Traditional Chinese Medicine Viola philippica

1
Research Center for Life Sciences Computing, Zhejiang Lab, Hangzhou 311100, China
2
Hainan Institute, Zhejiang University, Sanya 572025, China
3
State Key Laboratory of Rice Biology, Key Laboratory of Biology of Crop Pathogens and Insects of Zhejiang Province, Institute of Biotechnology, Zhejiang University, Hangzhou 310058, China
*
Authors to whom correspondence should be addressed.
Viruses 2025, 17(1), 114; https://doi.org/10.3390/v17010114
Submission received: 9 December 2024 / Revised: 10 January 2025 / Accepted: 14 January 2025 / Published: 15 January 2025
(This article belongs to the Special Issue Emerging and Reemerging Plant Viruses in a Changing World)

Abstract

:
Viola philippica, a medicinal herbaceous plant documented in the Chinese Pharmacopoeia, is a promising candidate for research into plant-derived pharmaceuticals. However, the study of newly emerging viruses that threaten the cultivation of V. philippica remains limited. In this study, V. philippica plants exhibiting symptoms such as leaf yellowing, mottled leaves, and vein chlorosis were collected and subjected to RNA sequencing to identify potential viral pathogens. A novel polerovirus, named Viola Philippica Polerovirus (VPPV), was identified in V. philippica. VPPV possesses a linear, positive-sense, single-stranded RNA genome consisting of 5535 nucleotides (nt) and encodes seven highly overlapping open reading frames (ORFs). Two potential recombination events were identified within ORF2, ORF3a, and ORF3, providing insights into the genetic diversity and evolution history of this novel polerovirus. An infectious cDNA clone of VPPV was successfully constructed and shown to infect Nicotiana benthamiana. Using a PVX-based heterologous expression system, the VPPV P0 protein was shown to trigger a systemic hypersensitive response (HR)-like reaction in N. benthamiana, indicating that P0 functions as the main pathogenicity determinant. These findings contributed to the detection and understanding of pathogenic mechanisms and control strategies for VPPV in V. philippica.

1. Introduction

Viola philippica Cav., an herbaceous perennial plant primarily indigenous to East Asia, has garnered significant attention in botany and pharmacology due to its well-established medicinal properties in traditional Chinese medicine [1,2]. Its efficacy in treating respiratory infections, skin inflammations, and fevers is attributed to its high content of bioactive compounds, such as flavonoids, alkaloids, and saponins, which exhibit notable antimicrobial, anti-inflammatory, and antiviral activities [3,4,5]. These characteristics make V. philippica a promising candidate for further research into plant-derived pharmaceuticals. However, while extensive research has focused on the therapeutic potential of V. philippica, its role as a host for plant viruses has not been explored.
The identification of novel viruses in plants has been revolutionized by advances in high-throughput sequencing technologies, particularly RNA sequencing (RNA-seq) [6]. As obligate intracellular parasites, viruses depend on host cellular machinery for replication and often exhibit significant genetic variability. This diversity presents major challenges for the detection and identification of previously unknown or highly divergent viral species using traditional methods, such as electron microscopy, serological assays, or polymerase chain reaction (PCR)/reverse transcription-PCR (RT-PCR) [7]. However, high-throughput sequencing has emerged as a powerful tool to overcome these limitations, providing comprehensive insights into plant viromes [8]. In recent years, RNA-seq has played a pivotal role in discovering new viral species in a wide range of plants, including grain crops [9,10,11], medicinal plants [12], and ornamental plants [13]. These discoveries have contributed to a deeper understanding of plant–virus interactions, viral evolution, and the potential threats posed by these novel pathogens to global agriculture.
Poleroviruses belong to the genus Polerovirus and have genomes comprising polycistronic, positive-sense, single-stranded RNA, typically ranging from 5.5 to 6.2 kb in size [14]. Their virions are non-enveloped icosahedral particles, 26–34 nm in diameter [14]. In 2021, the Polerovirus genus was reclassified from the family Luteoviridae to the family Solemoviridae (order Sobelivirales) [15]. A characteristic Polerovirus genome features a covalently attached virus protein genome-linked (VPg) at the 5′-end and lacks a poly(A) tail at the 3′-end. The genome typically encodes five to seven overlapping open reading frames (ORFs), some of which are translated from subgenomic RNAs. ORF0, ORF1, and ORF2 are translated from genomic RNA, while ORF3a to ORF7 are translated from subgenomic RNA [16]. These viral proteins are primarily involved in RNA interference suppression, genome replication, virion assembly, in-plant movement, and vector transmission [17,18,19].
Poleroviruses are primarily transmitted by aphids in a circulative, non-propagative manner, causing significant economic losses in agricultural crops worldwide. The phloem-limited nature of poleroviruses restricts their movement within plants, resulting in systemic symptoms such as leaf vein yellowing, stunting, and reduced vigor, which profoundly impact crop productivity and quality [20,21]. While extensive research was conducted on poleroviruses affecting major crops, such as cereals, vegetables, and fruits, their presence and impact in traditional Chinese medicinal plants, such as V. philippica remain unexplored [22,23]. In this study, we utilized RNA-seq to identify viral pathogens in V. philippica plants exhibiting virus-like symptoms and discovered a novel virus belonging to the genus Polerovirus, named Viola philippica polerovirus (VPPV). Furthermore, an infectious cDNA clone of VPPV was constructed, and the virulence determinants of VPPV were investigated.

2. Materials and Methods

2.1. Plant Materials and Plant Growth

In April 2021, a V. philippica plant showing symptoms of yellowing, mottling, and vein chlorosis was sampled in the Summer Palace in Beijing, China. Nicotiana benthamiana plants were cultivated in a growth chamber maintained at 25 °C (day) and 22 °C (night), with a relative humidity of 60% and a photoperiod of 16 h (light) and 8 h (dark).

2.2. RNA-Seq and Data Analyses

RNA-seq and data analyses were performed as previously described [10]. Briefly, the total RNA was extracted from the leaf tissue of the collected sample using the MiniBEST Plant RNA Extraction Kit (TaKaRa, Tokyo, Japan) according to the manufacturer’s instructions. The extracted total RNA was treated with the TruSeq RNA Sample Prep Kit (Illumina, San Diego, CA, USA) to remove ribosomal RNAs. Subsequently, libraries were prepared and sequenced on an Illumina HiSeq X-ten platform (Biomarker Technologies, Beijing, China). Raw sequencing data underwent quality filtering and adapter trimming using CLC Genomics Workbench 9.5 software (QIAGEN, Germantown, MD, USA), resulting in high-quality clean reads that were de novo assembled into contigs using Trinity v2.3.2 program (Broad Institute and the Hebrew University of Jerusalem, Cambridge, MA, USA and Jerusalem, Israel). These contigs were then subjected to BLASTx and BLASTp searches in NCBI databases as described by Wang et al. [10].

2.3. Reverse Transcription-Polymerase Chain Reaction (RT-PCR) and 5′/3′ Rapid Amplification of cDNA Ends (5′/3′ RACE)

To obtain the full-length genome sequences of VPPV, total RNA was extracted from the collected plant samples. Reverse transcription (RT) was conducted using the ReverTra Ace qPCR RT Master Mix with gDNA Remover (Toyobo, Osaka, Japan) to generate cDNAs, followed by PCR to detect the virus in the collected plant samples using the Green Taq Mix (Vazyme, Nanjing, China) and specific primer pairs (Table S1) according to the assembled contig sequences. For the accurate amplification of viral genomic fragments, RT was performed using M-MLV reverse transcriptase (TaKaRa, Kusatsu, Japan) and fragment-specific primers, followed by super-fidelity PCR using KOD One PCR Master Mix (TOYOBO, Osaka, Japan) and specific primers (Table S2). To amplify the 3′ and 5′ ends of viral genomic RNA, RACE was performed using the HiScript-TS 5′/3′ RACE Kit (Vazyme, Nanjing, China) according to the manufacturer’s instructions. Notably, Poly(U) polymerase (NEB) was used to add a poly(A) tail to 3′ end of viral RNA genome, followed by RT using M-MLV and an oligo-dT-adapter primer to generate cDNAs. Subsequently, nested-PCR amplification was performed using gene-specific primers (GSPs) and adapter primers (Table S3). The amplified products of the expected size were gel-purified and cloned into the pCE2 TA/Blunt-Zero vector (Vazyme) for Sanger sequencing. The resulting sequences were checked and used to assemble the complete viral RNAs.

2.4. Organization Analysis of Viral Genome and Prediction of Viral Proteins

For open reading frame (ORF) prediction, the viral genomic sequences were analyzed using the ORFfinder service (http://www.ncbi.nlm.nih.gov/orffinder, accessed on 25 January 2021). Amino acid sequences of the predicted viral proteins were then searched against the GenBank database using BLASTp. Conserved protein domains within the deduced viral proteins were identified utilizing the Conserved Domain Search Service (CD-Search) software available through NCBI (https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi, accessed on 9 December 2024). Genomic and protein sequences of closely related viruses were obtained from the GenBank database at NCBI (Table S4) and aligned with VPPV genomic and protein sequences using DNASTAR7.1 software (https://www.dnastar.com/, accessed on 9 December 2024).

2.5. Recombination and Phylogenetic Analysis

To assess viral recombination, the complete genomic sequences of VPPV and closely related viruses, including chickpea chlorotic stunt virus (CpCSV, AY956384.1), pepper vein yellows virus (PeVYV, KU999109.1), melon aphid-borne yellows virus (MABYV, EU000534.1), potato leafroll virus (PLRV, OQ446811.1), cotton leafroll dwarf virus (CLDV, GU167940.1), and cereal yellow dwarf virus (CYDV, AF235168.2), were aligned using clustalW, followed by recombination analysis conducted with the RDP5 program (http://web.cbio.uct.ac.za/~darren/rdp.html, accessed on 9 December 2024), as described by Martin et al. [24]. To elucidate, the evolutionary relationship of the virus, the amino acid sequences of RdRps or CP from VPPV, and closely related poleroviruses obtained from the GenBank database were aligned using clustalW. Subsequently, phylogenetic trees were subsequently constructed using MEGA X software (https://www.megasoftware.net/, accessed on 9 December 2024). Detailed sequence information for this study is provided in Table S4.

2.6. Plasmid Construction

The complete genomic sequence of VPPV was amplified by RT-PCR and cloned into the linear expression vector pCB301-RZ using MonCloneTM Hi-Fusion Cloning Mix V2 (Monad, Suzhou, China) following the manufacturer’s instructions. The viral protein genes were individually amplified from infected V. philippica plants via RT-PCR using specific primers listed in Table S5. The amplified products were subcloned into the linear expression vectors pGD or pGR106 utilizing MonCloneTM Single Assembly Cloning Mix (Monad, Suzhou, China) or MonCloneTM Hi-Fusion Cloning Mix V2 (Monad, Suzhou, China), according to the manufacturer’s instructions.

2.7. Agroinfiltration Assays

The recombinant vectors pGD or pCB301-RZ were introduced into Agrobacterium tumefaciens strain EHA105, while the pGR106-based vectors were individually introduced into A. tumefaciens strain GV3101 using the freeze-thawing procedure. The transformed A. tumefaciens cultures were diluted to OD600 = 0.6 with infiltration buffer (10 mM MgCl2, 10 mM MES, pH 5.6, and 100 μM acetosyringone) and then infiltrated into the N. benthamiana leaves.

2.8. Protein Extraction and Western Blot Assays

Plant leaves were ground to a fine powder in liquid nitrogen and homogenized in extraction buffer (50 mM Tris-HCl, pH 6.8, 9 M carbamide, 4.5% SDS, and 7.5% 2-mercaptoethanol). The resulting protein samples were subjected to electrophoresis on a 12% sodium dodecyl sulfate-polyacrylamide gel and subsequently transferred onto nitrocellulose membranes. The membranes were blocked with 5% skimmed milk in a 0.01 M phosphate-buffered saline (PBS) for 30 min at 37 °C, followed by incubation in anti-PVX-CP mouse monoclonal antibody for 60 min at 37 °C. After three rinses in PBS containing 0.05% Tween-20, the membranes were incubated in a goat anti-mouse secondary antibody conjugated horseradish peroxidase for another 60 min at 37 °C. The resulting nitrocellulose membranes were analyzed using ECL substrate and photographed by ImageQuant LAS 4000mini (GE HealthCare, New York, NY, USA).

3. Results

3.1. Identification of a Novel Polerovirus in Viola philippica

A V. philippica plant exhibiting symptoms of leaf vein yellowing, mottling, and vein chlorosis was sampled (Figure 1a,b). Total RNA was extracted and subjected to RNA-seq analysis to identify the viral pathogen(s) responsible for the observed symptoms. The de novo assembly of high-quality clean reads was performed, followed by the screening of contigs longer than 1000 bp by using BLASTx algorithms. The results identified two contigs with high sequence similarity to polerovirus: contig 1345 (4044 nt) and contig 21,591 (1465 nt). To validate the RNA-seq findings, RT-PCR assays using contig-specific primers (Table S1) successfully amplified segments from both contigs (586 bp and 596 bp) in the symptomatic sample (Figure 1a) but not from asymptomatic leaf (Figure 1c). Sanger sequencing confirmed the presence of contigs 1345 and 21,591 in the symptomatic V. philippica sample.
Based on the genome organization of known poleroviruses (Figure 1e), we speculated that contig 1345 represents the upstream portion of the novel viral genome, while contig 21,591 corresponds to the downstream portion, with an unidentified sequence gap between these two contigs. To obtain the complete viral genome, 5′ and 3′ RACE and super-fidelity RT-PCR amplification were separately conducted to determine the terminal sequences and fill the sequence gap between these two contigs. The results revealed that the novel virus has a positive-sense, single-stranded RNA genome of 5535 nt in length, including a 69-nt 5′ UTR that begins with the conserved sequence ACAAAAGA [25], and a 179-nt 3′ UTR, excluding the poly(A) tail (Figure 1e). Genome sequence alignment analysis demonstrated that this novel virus shares 36.7% to 56.1% genome identity with twelve known viruses in the genus Polerovirus, with the highest identity observed with chickpea chlorotic stunt virus (CpCSV) (Table 1). The complete genome sequence of the novel virus was deposited in Genbank under accession number PP770488 and was designated as viola philippica polerovirus (VPPV).

3.2. Genome Organization of VPPV

The VPPV genome contains seven putative ORFs, labeled ORF0 to ORF5 and ORF3a, which are predicted to encode proteins P0 to P5 and P3a. ORF0 (nt positions 70–795) is predicted to encode a 26.9 kDa RNA silencing suppressor protein (P0), which shares the highest amino acid sequence identity (24.2%) with P0 of pepper vein yellows virus 4 (PeVYV-4) (Table 1). The ORF1 (nt positions 197–2002) encodes a 66.1 kDa polyprotein (P1), exhibiting amino acid sequence identity ranging from 26.3% to 39.1% with P1 proteins from other known poleroviruses (Table 1). P1 is predicted to contain a conserved Peptidase_S39 superfamily domain (cl09540, 215–414 aa), which is involved in protein cleavage activity.
Given that ORF1 contains a conserved slippery heptanucleotide sequence 1564GGGAAAC1570, it is anticipated that a fusion protein, i.e., RdRp weighing 115.5 kDa, is produced via a −1 ribosomal frameshift between ORF1 and ORF2 (nt positions 197–1570 and 1570–3282). The RdRp features a conserved ps-ssRNAv_RdRp-like super family domain (cl40470, 552–1007 aa), which is essential for viral replication. VPPV RdRp shares amino acid sequence identity ranging from 44.2% to 56.2% with RdRps from twelve representative poleroviruses, with the highest identity observed with the RdRp of melon aphid-borne yellows virus (MABYV) (Table 1).
ORF3 (nt positions 3461–4021) encodes a putative 21.2 kDa coat protein (CP; P3), which shares amino acid sequence identity ranging from 30.9% to 37.2% with CPs from known poleroviruses (Table 1). P3 is predicted to contain a Luteo_coat superfamily domain (cl03007, 53–180 aa). Between the translation initiation sites of RdRp and CP, there is a 178-nt intergenic region. It is hypothesized that the non-AUG-initiated protein P3a begins at nt position 3305 and terminates at a UAA stop codon at nt position 3457.
ORF4 (nt positions 3561–3971), which overlaps with ORF3 but is in a different reading frame, is predicted to encode a movement protein (MP; P4) of 15.1 kDa. VPPV P4 shares the highest amino acid identity (26.6%) with P4 of potato leafroll virus (PLRV) (Table 1). ORF5 (nt positions 3461–5356) is translated via an in-frame readthrough of the ORF3 stop codon, encoding a minor capsid component of 71.0 kDa, also referred to as the read through protein (RTP). VPPV RTP shares amino acid identity ranging from 21.5% to 30.0% with P5 of other known poleroviruses (Table 1).

3.3. Recombination and Phylogenetic Analysis

Poleroviruses exhibit a propensity for recombination both within their own genus and with viruses from other genera. To assess the frequency and distribution of potential recombination events in VPPV, we aligned genome sequences of VPPV (PP770488.1), CpCSV (AY956384.1), PeVYV (KU999109.1), MABYV (EU000534.1), PLRV (OQ446811.1), CLDV (GU167940.1), and CYDV (AF235168.2) using the recombination program RDP5 [24]. A putative recombination event in VPPV, with CpCSV as the major parent and PeVYV as the minor parent, was strongly supported by the GENECONV method, with p-values of 1.043 × 10−24 (Figure 2a). The predicted beginning and ending breakpoints are located at nt positions 2020 and 3622 in the VPPV genome. This region is a known “hot spot” for recombination events in poleroviruses [26]. The recombinant fragment extends from the end of ORF1 to the non-coding intergenic region (IR) of the VPPV sub-genome and encodes the C-terminal part of RdRp (P1-P2) (Figure 2a). Another putative recombination event, spanning nt positions 3620–3874 within the VPPV genome, involves PeVYV as the major parent and CpCSV or MABYV as the minor parents. This event was strongly supported by the MaxChi (p-values: 1.676 × 10−23), Chimaera (p-values: 6.575 × 10−5), and 3Seq methods (p-values: 3.819 × 10−14) (Figure 2b). This recombinant fragment spans from ORF3a to ORF3 and encodes the near-complete P3a protein and the N-terminal part of the CP (P3a-P3) (Figure 2b).
To investigate the evolutionary relationship between VPPV and other members of the family Solemoviridae, we aligned the amino acid sequences of RdRp or CP proteins from 19 representative viruses across different genera within the family Solemoviridae, alongside with 3 viruses from the family Tombusviridae, to construct a phylogenetic tree (Table S4). The phylogenetic tree based on RdRp clearly indicated that VPPV clustered with twelve representative poleroviruses and exhibited the closest genetic relationship with MABYV (Figure 2c), suggesting its classification as a member of the Polerovirus genus. Additionally, the phylogenetic tree based on viral CP further supported this assignment of VPPV to the Polerovirus genus (Figure 2d).

3.4. Construction of the Infectious cDNA Clone of VPPV

To gain insights into virus biology, an infectious clone derived from a complete viral genome is an important tool. The full-length cDNA of VPPV was cloned into the binary vector pCB301 (pCB301-VPPV), which features the cauliflower mosaic virus 35S promoter to generate capped VPPV viral transcripts, and the autolytic hepatitis delta virus ribozyme (HDV Rz) to release the authentic VPPV viral 3′ end. N. benthamiana plants were infiltrated with A. tumefaciens harboring pCB301-VPPV or the empty pCB301 vector. At 12 days post infiltration (dpi), newly emerging leaves from the infiltrated N. benthamiana plants were collected and subjected to RT-PCR using VPPV-specific primer pairs (Table S1). RT-PCR products of the expected size (~586 bp) were detected in the systemic leaves of pCB301-VPPV infiltrated plants, while no products were detected in the negative control pCB301-infiltrated (CK-) plants (Figure 3a). This confirmed successful viral RNA amplification and VPPV infection in pCB301-VPPV-infiltrated N. benthamiana plants. However, despite the establishment of a stable viral infection, no discernible virus-like symptoms were observed in the pCB301-VPPV-infiltrated N. benthamiana plants (Figure 3b), suggesting that VPPV cause asymptomatic infection in its non-natural host, N. benthamiana.

3.5. Virulence Determinants of VPPV

Since the infectious cDNA clone of VPPV causes asymptomatic infection in N. benthamiana, a potato X virus (PVX)-based heterologous system was employed to better understand the pathology of this novel polerovirus and identify its viral virulence determinants. Seven PVX-based recombinant vectors were constructed: PVX-P0, PVX-P1, PVX-RdRp, PVX-P3a, PVX-P3, PVX-P4, and PVX-RTP. Notably, an AUG start codon was added before the putative start codon (AUA) of P3a to initiate non-canonical P3a expression. Additionally, the ORF3 stop codon (UGA) was mutated to GCG (Ala) to enable full-length ORF5 expression of the fusion protein RTP. A. tumefaciens carrying these expression vectors was individually infiltrated into N. benthamiana leaves, with the empty PVX vector serving as a negative control for comparison.
At 3 days post-infiltration (dpi), plants inoculated with PVX-P0 developed severe necrotic lesions on the infiltrated leaves. By 5 dpi, necrotic symptoms extended to the upper systemic leaves, ultimately leading to severe wilting and death of the plant (Figure 3c). These results indicate that VPPV P0 can elicit a hypersensitive response (HR) in N. benthamiana. At 10 dpi, PVX-P3a-infected plants developed punctate necrosis on the leaves, while plants infiltrated with PVX-RdRp, PVX-P3, PVX-P4, or PVX-RTP displayed only curling and mild mosaic symptoms, similar to those caused by negative control PVX infection (Figure 3c).
To confirm PVX infection and measure PVX accumulation in the infiltrated plants, we performed a Western blot assay using an anti-PVX coat protein (CP) antibody at 10 dpi (Figure 3d). The results confirmed successful infection by all PVX recombinant vectors in infiltrated N. benthamiana. However, due to the severe necrosis induced by PVX-P0, protein extraction from the PVX-P0 sample was unsuccessful (Figure 3d). Interestingly, PVX-P1 failed to induce any typical PVX-associated symptoms, such as leaf curling and mosaic patterns, throughout the entire observation period (Figure 3c). This lack of symptoms was consistent with our PVX accumulation assay, which showed significantly reduced PVX CP accumulation in PVX-P1 samples compared to the negative control (Figure 3d). Based on these results, we concluded that P0 and P3a are likely the pathogenicity determinants of VPPV.

4. Discussion

V. philippica has attracted considerable attention in the traditional Chinese medicine market due to its medicinal properties [1,2]. However, in natural ecosystems, V. philippica is not only susceptible to microbial pathogens but also exposed to a variety of plant-specific viruses, which pose a detrimental impact on the Chinese medicine planting industry. To date, no reports have documented viral infections in V. philippica plants. In this study, we employed RNA-seq analysis to investigate the presence of emerging viruses in herbaceous V. philippica exhibiting typical viral symptoms, such as leaf yellowing, mottling, and vein chlorosis. Our analysis led to the identification of a novel polerovirus, VPPV, and the determination of its complete genome sequence. To our knowledge, this is the first report of a novel polerovirus infecting V. philippica. According to the latest species demarcation criteria for the Polerovirus genus outlined by International Committee on Taxonomy of Viruses (ICTV), a difference in amino acid sequence identity greater than 10% in any gene product is required for species classification [14]. Our finding indicated that VPPV qualifies as a distinct species within the Polerovirus genus based on its viral genomic structure (Figure 1e), amino acid homology (Table 1), and phylogenetic tree analysis (Figure 2c,d).
Viruses in the economically important Polerovirus genus are known for their complex interactions with both their plant hosts and insect vectors [27,28]. The VPPV genomic RNA encodes three proteins, P0, P1, and RdRp, while its subgenomic RNA encodes four proteins: CP, MP, RTP, and a short non-AUG-initiated P3a. The polerovirus P0 protein is a viral suppressor of RNA silencing, targeting the ARGONAUTE1(AGO1) PAZ motif and its adjacent upstream sequence for degradation [29]. Furthermore, it was demonstrated that RNA silencing suppressors heterologous expressed by PVX can induce hypersensitive response [30,31]. Consistent with previous findings, the heterologous expression of the VPPV P0 protein by PVX resulted in plant necrosis (Figure 3c), suggesting that P0 is a key pathogenicity determinant, acting as a viral suppressor of RNA silencing and eliciting a systemic HR-like response in PVX-associated synergism.
The polyprotein (P1) in poleroviruses is translated into at least two forms (P1 and P1/P2 transframe protein)and contains a serine proteinase domain, serving as a precursor for VPg generation [32]. Further investigation is needed to clarify the self-cleavage activity and processed products of VPPV P1. The expression of VPPV P1 in heterologous PVX system mitigated typical PVX-associated symptoms (Figure 3c) and reduced PVX CP accumulation (Figure 3d), possibly due to its nucleic acid-binding properties, which negatively impact PVX replication. This is consistent with previous findings that PLRV P1/2 and P1 possess sequence-unrelated binding domains involved in PLRV RNA replication [33]. The RdRp (P1/2) of poleroviruses is a translational fusion protein generated via a -1 ribosomal frameshift mechanism [34,35]. Based on these findings, VPPV RdRp is speculated to be produced by ORF1 and ORF2 undergoing a −1 ribosomal frameshift at the slippery sequence 1564GGGAAAC1570.
The non-AUG-initiated P3a of polerovirus, encoded by ORF3a, is essential for viral long-distance movement and phloemrestriction [19]. We speculate that VPPV P3a may play a crucial role in intercellular viral movement and systemic viral infection in plants. The expression of VPPV P3a by PVX did not significantly alter the mild mosaic symptoms typically observed in PVX-infected N. benthamiana plants. However, occasional small scattered bleached spots were observed (Figure 3c), suggesting that this non-canonical-initiated viral protein may have subtle effects on the symptom phenotype.
In accordance with other poleroviruses, the coding region of the VPPV genomic ORF3 contains a P4, which is encoded by an alternative reading frame (Figure 1e) and is potentially associated with viral cell-to-cell movement [36]. ORF5 is expressed exclusively through an in-frame readthrough of the stop codon in ORF3 (Figure 1e), resulting in the production of RTP, a fusion protein with a C-terminus extension of CP [37,38]. VPPV RTP may facilitate viral transport through aphid salivary glands membranes, as poleroviruses and luteoviruses carrying mutations within RTP are not aphid transmissible [39,40,41].
Poleroviruses exhibit frequent intra- and interspecific recombination, particularly in regions associated with ORF2 and ORF3a [42,43]. Our findings reveal two potential recombination events concentrated at the boundaries of RdRp, P3a, and CP genes in VPPV (Figure 2a,b). One recombination event was observed in RdRp, with CpCSV identified as the major parent (Figure 2a). Another recombination event occurred in P3a and the N-terminal of CP, with PeVYV as the major parent (Figure 2b). These findings strongly suggest distinct evolutionary histories for these three genes in VPPV. This is further supported by significant incongruence between the phylogenetic trees of the RdRp and CP gene (Figure 2c,d). The clustering of recombination breakpoints in these three genes may be due to certain RdRp-P3a-CP combinations and confer a selective advantage in specific host/vector combinations [44]. These results provide valuable insights into the genetic diversity and evolution of poleroviruses.
Understanding the biology, epidemiology, and pathogenic mechanisms of poleroviruses is critical for developing effective management strategies. We successfully constructed a full-length infectious cDNA clone of VPPV and found that VPPV can asymptomatically infect N. benthamiana (Figure 3a,b). However, this infectivity assay should be further confirmed in its native host plant, V. philippica. The VPPV infectious cDNA clone will serve as a tool to investigate important viral biological traits, such as aphid vector interactions, susceptible hosts, the epidemiology, and impacts of this virus on the productivity of V. philippica. Furthermore, the role of V. philippica as a potential virus reservoir raises critical questions about the transmission pathways of poleroviruses and their impact on both wild and cultivated plants in shared environments. Understanding the interaction between V. philippica and poleroviruses will provide broader insights into viral epidemiology, evolution, and the development of plant virus resistance strategies.

5. Conclusions

This study identified and characterized a novel polerovirus, VPPV, which infects the medicinal herb V. philippica. A full-length infectious cDNA clone of VPPV was developed, and the P0 and P3a proteins were determined as two key pathogenic determinants. The discovery of VPPV infection in V. philippica advances our understanding of polerovirus diversity and evolutionary dynamics, laying the groundwork for future studies on its impact and management strategies in traditional Chinese medicinal plants.

Supplementary Materials

The following supporting information can be downloaded at: www.mdpi.com/article/10.3390/v17010114/s1, Table S1. Specific primers used to detect the VPPV contigs; Table S2. Specific primers used for RT-RCR amplification of VPPV genome sequence; Table S3. Specific primers used to perform RACE; Table S4. Virus information used to analyze the nucleotide and amino acid identities, recombination, and construct phylogenetic trees; Table S5. Specific primers used to construct vectors.

Author Contributions

Conceptualization, S.F. (Shuai Fu) and J.W.; writing—original draft preparation, S.F. (Shuai Fu) and Y.C.; writing—review and editing S.F. (Shuai Fu) and J.W.; funding acquisition, S.F. (Shuai Fu) and J.W.; methodology, and investigation, Y.C., J.Y., Y.Z. and S.F. (Shifang Fei); software, formal analysis, and data curation, G.C. and H.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the “Pioneer” and “Leading Goose” R&D Program of Zhejiang (2024SSYS0007; 2024C03004), the National Key Research and Development Program of China (2022YFD1401200), the Hainan Province Science and Technology Special Fund (ZDYF2022XDNY240), and the Research Fund of Zhejiang Lab (119003-PI2202).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data generated and analyzed in this study are available upon request.

Acknowledgments

We thank Xiaoxiao Feng (the Agricultural Experiment Station, Zhejiang University) for her assistance during the glasshouse plant planting.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Cao, D.L.; Zhang, X.J.; Xie, S.Q.; Fan, S.J.; Qu, X.J. Application of chloroplast genome in the identification of traditional Chinese medicine Viola philippica. BMC Genom. 2022, 23, 540. [Google Scholar] [CrossRef] [PubMed]
  2. Lu, G.; Qiao, J.; Wang, L.; Liu, H.; Wu, G.; Zhu, Y.; Zhao, Y.; Xie, G.; Qin, M. An integrated study of Violae Herba (Viola philippica) and five adulterants by morphology, chemical compositions and chloroplast genomes: Insights into its certified plant origin. Chin. Med. 2022, 17, 32. [Google Scholar] [CrossRef] [PubMed]
  3. Xie, C.; Kokubun, T.; Houghton, P.J.; Simmonds, M.S. Antibacterial activity of the Chinese traditional medicine, Zi Hua Di Ding. Phytother. Res. 2004, 18, 497–500. [Google Scholar] [CrossRef] [PubMed]
  4. Cao, P.; Yang, Y.; Uche, F.I.; Hart, S.R.; Li, W.W.; Yuan, C. Coupling Plant-Derived Cyclotides to Metal Surfaces: An Antibacterial and Antibiofilm Study. Int. J. Mol. Sci. 2018, 19, 793. [Google Scholar] [CrossRef] [PubMed]
  5. Wang, Y.L.; Zhang, L.; Li, M.Y.; Wang, L.W.; Ma, C.M. Lignans, flavonoids and coumarins from Viola philippica and their α-glucosidase and HCV protease inhibitory activities. Nat. Prod. Res. 2019, 33, 1550–1555. [Google Scholar] [CrossRef]
  6. Villamor, D.E.V.; Ho, T.; Al Rwahnih, M.; Martin, R.R.; Tzanetakis, I.E. High throughput sequencing for plant virus detection and discovery. Phytopathology 2019, 109, 716–725. [Google Scholar] [CrossRef]
  7. Mehetre, G.T.; Leo, V.V.; Singh, G.; Sorokan, A.; Maksimov, I.; Yadav, M.K.; Upadhyaya, K.; Hashem, A.; Alsaleh, A.N.; Dawoud, T.M.; et al. Current Developments and Challenges in Plant Viral Diagnostics: A Systematic Review. Viruses 2021, 13, 412. [Google Scholar] [CrossRef]
  8. Bejerman, N.; Debat, H.; Dietzgen, R.G. The Plant Negative-Sense RNA Virosphere: Virus Discovery Through New Eyes. Front. Microbiol. 2020, 11, 588427. [Google Scholar] [CrossRef]
  9. Yu, J.; Zeng, M.; Zhou, Y.; Wang, J.; Zhou, X.; Wu, J. Discovery of a novel whitefly- and aphid-transmitted polerovirus on rice plants with dwarfing and fewer tillering symptoms. Crop Health 2024, 2, 13. [Google Scholar] [CrossRef]
  10. Wang, D.; Fu, S.; Wu, H.; Cao, M.; Liu, L.; Zhou, X.; Wu, J. Discovery and genomic function of a novel rice dwarf-associated Bunya-like virus. Viruses 2022, 14, 1183. [Google Scholar] [CrossRef]
  11. Xie, Y.; Fu, S.; Xie, L.; Wang, Y.; Cao, M.; Zhou, X.; Wu, J. Identification and characterization of two novel Noda-like viruses from rice plants showing the dwarfing symptom. Viruses 2022, 14, 1159. [Google Scholar] [CrossRef] [PubMed]
  12. Zheng, L.; Fu, S.; Xie, Y.; Han, Y.; Zhou, X.; Wu, J. Discovery and characterization of a novel umbravirus from Paederia scandens plants showing leaf chlorosis and yellowing symptoms. Viruses 2022, 14, 1821. [Google Scholar] [CrossRef] [PubMed]
  13. Dong, J.; Chen, Y.; Xie, Y.; Cao, M.; Fu, S.; Wu, J. The identification of viral pathogens in a Physostegia virginiana plant using high-throughput RNA sequencing. Viruses 2023, 15, 1972. [Google Scholar] [CrossRef] [PubMed]
  14. Sõmera, M.; Fargette, D.; Hébrard, E.; Sarmiento, C.; Ictv Report Consortium. ICTV virus taxonomy profile: Solemoviridae. J. Gen. Virol. 2021, 102, 001707. [Google Scholar]
  15. Walker, P.J.; Siddell, S.G.; Lefkowitz, E.J.; Mushegian, A.R.; Adriaenssens, E.M.; Alfenas-Zerbini, P.; Harrach, B.; Harrison, R.L.; Junglen, S.; Knowles, N.J.; et al. Changes to virus taxonomy and to the International Code of Virus Classification and Nomenclature ratified by the International Committee on Taxonomy of Viruses. Arch. Virol. 2021, 166, 2633–2648. [Google Scholar] [CrossRef]
  16. Sõmera, M.; Sarmiento, C.; Truve, E. Overview on Sobemoviruses and a Proposal for the Creation of the Family Sobemoviridae. Viruses 2015, 7, 3076–3115. [Google Scholar] [CrossRef]
  17. Kaplan, I.B.; Lee, L.; Ripoll, D.R.; Palukaitis, P.; Gildow, F.; Gray, S.M. Point mutations in the potato leafroll virus major capsid protein alter virion stability and aphid transmission. J. Gen. Virol. 2007, 88 Pt 6, 1821–1830. [Google Scholar] [CrossRef]
  18. Csorba, T.; Lózsa, R.; Hutvágner, G.; Burgyán, J. Polerovirus protein P0 prevents the assembly of small RNA-containing RISC complexes and leads to degradation of ARGONAUTE1. Plant J. 2010, 62, 463–472. [Google Scholar] [CrossRef]
  19. Smirnova, E.; Firth, A.E.; Miller, W.A.; Scheidecker, D.; Brault, V.; Reinbold, C.; Rakotondrafara, A.M.; Chung, B.Y.-W.; Ziegler-Graff, V. Discovery of a small non-AUG-initiated ORF in poleroviruses and luteoviruses that is required for long-distance movement. PLoS Pathog. 2015, 11, e1004868. [Google Scholar] [CrossRef]
  20. Malmstrom, C.M.; Alexander, H.M. Effects of crop viruses on wild plants. Curr. Opin. Virol. 2016, 19, 30–36. [Google Scholar] [CrossRef]
  21. Edula, S.R.; Bag, S.; Milner, H.; Kumar, M.; Suassuna, N.D.; Chee, P.W.; Kemerait, R.C.; Hand, L.C.; Snider, J.L.; Srinivasan, R.; et al. Cotton leafroll dwarf disease: An enigmatic viral disease in cotton. Mol. Plant Pathol. 2023, 24, 513–526. [Google Scholar] [CrossRef] [PubMed]
  22. Aradottir, G.I.; Crespo-Herrera, L. Host plant resistance in wheat to barley yellow dwarf viruses and their aphid vectors: A review. Curr. Opin. Insect Sci. 2021, 45, 59–68. [Google Scholar] [CrossRef] [PubMed]
  23. Miller, W.A.; Lozier, Z. Yellow Dwarf Viruses of Cereals: Taxonomy and Molecular Mechanisms. Annu. Rev. Phytopathol. 2022, 60, 121–141. [Google Scholar] [CrossRef] [PubMed]
  24. Martin, D.P.; Varsani, A.; Roumagnac, P.; Botha GMaslamoney, S.; Schwab, T.; Kelz, Z.; Kumar, V.; Murrell, B. RDP5: A computer program for analyzing recombination in, and removing signals of recombination from, nucleotide sequence datasets. Virus Evol. 2021, 7, veaa087. [Google Scholar] [CrossRef]
  25. Gudeta, W.F.; Shin, A.Y.; Kim, S.E.; Jeong-A, K.; Seok-Yoon, K.; Moon, J.S. Complete genome sequence of stellaria aquatica virus B, a novel polerovirus that infects Stellaria aquatica. Arch. Virol. 2023, 168, 22. [Google Scholar] [CrossRef]
  26. Pagán, I.; Holmes, E.C. Long-term evolution of the Luteoviridae: Time scale and mode of virus speciation. J. Virol. 2010, 84, 6177–6187. [Google Scholar] [CrossRef]
  27. Delfosse, V.C.; Barrios Barón, M.P.; Distéfano, A.J. What we know about poleroviruses: Advances in understanding the functions of polerovirus proteins. Plant Pathol. 2021, 70, 1047–1061. [Google Scholar] [CrossRef]
  28. Patton, M.F.; Bak, A.; Sayre, J.M.; Heck, M.L.; Casteel, C.L. A polerovirus, Potato leafroll virus, alters plant-vector interactions using three viral proteins. Plant Cell Environ. 2020, 43, 387–399. [Google Scholar] [CrossRef]
  29. Bortolamiol, D.; Pazhouhandeh, M.; Marrocco, K.; Genschik, P.; Ziegler-Graff, V. The Polerovirus F box protein P0 targets ARGONAUTE1 to suppress RNA silencing. Curr. Biol. 2007, 17, 1615–1621. [Google Scholar] [CrossRef]
  30. Scholthof, H.B.; Scholthof, K.B.; Jackson, A.O. Identification of tomato bushy stunt virus host-specific symptom determinants by expression of individual genes from a potato virus X vector. Plant Cell 1995, 7, 1157–1172. [Google Scholar]
  31. Wang, X.; Luo, C.; Xu, Y.; Zhang, C.; Bao, M.; Dou, J.; Wang, Q.; Cheng, Y. Expression of the p24 silencing suppressor of Grapevine leafroll-associated virus 2 from Potato virus X or Barley stripe mosaic virus vector elicits hypersensitive responses in Nicotiana benthamiana. Plant Physiol. Biochem. 2019, 142, 34–42. [Google Scholar] [CrossRef] [PubMed]
  32. Li, X.; Ryan, M.D.; Lamb, J.W. Potato leafroll virus protein P1 contains a serine proteinase domain. J. Gen. Virol. 2000, 81 Pt 7, 1857–1864. [Google Scholar] [CrossRef] [PubMed]
  33. Prüfer, D.; Kawchuk, L.; Monecke2, M.; Nowok, S.; Fischer, R.; Rohde, W. Immunological analysis of potato leafroll luteovirus (PLRV) P1 expression identifies a 25 kDa RNA-binding protein derived via P1 processing. Nucleic Acids Res. 2019, 27, 421–425. [Google Scholar] [CrossRef]
  34. Jacks, T.; Madhani, H.D.; Masiarz, F.R.; Varmus, H.E. Signals for ribosomal frameshifting in the rous sarcoma virus gag-pol region. Cell 1988, 55, 447–458. [Google Scholar] [CrossRef]
  35. Nixon, P.L.; Cornish, P.V.; Suram, S.V.; Giedroc, D.P. Thermodynamic analysis of conserved loop-stem interactions in P1-P2 frameshifting RNA pseudoknots from plant Luteoviridae. Biochemistry 2002, 41, 10665–10674. [Google Scholar] [CrossRef]
  36. Shakir, S.; Boissinot, S.; Michon, T.; Lafarge, S.; Zaidi, S.S. Beyond movement: Expanding functional landscape of luteovirus movement proteins. Trends Plant Sci. 2024, 29, 1331–1341. [Google Scholar] [CrossRef]
  37. Bahner, I.; Lamb, J.; Mayo, M.A.; Hay, R.T. Expression of the genome of potato leafroll virus: Readthrough of the coat protein termination codon in vivo. J. Gen. Virol. 1990, 71 Pt 10, 2251–2256. [Google Scholar] [CrossRef]
  38. Boissinot, S.; Erdinger, M.; Monsion, B.; Ziegler-Graff, V.; Brault, V. Both structural and non-structural forms of the readthrough protein of cucurbit aphid-borne yellows virus are essential for efficient systemic infection of plants. PLoS ONE 2014, 9, e93448. [Google Scholar] [CrossRef]
  39. Brault, V.; van den Heuvel, J.F.; Verbeek, M.; Ziegler-Graff, V.; Reutenauer, A.; Herrbach, E.; Garaud, J.; Guilley, H.; Richards, K.; Jonard, G. Aphid Transmission of beet western yellows luteovirus requires the minor capsid read-through protein P74. EMBO J. 1995, 14, 650–659. [Google Scholar] [CrossRef]
  40. Chay, C.A.; Gunasinge, U.B.; Dinesh-Kumar, S.P.; Miller, W.A.; Gray, S.M. Aphid transmission and systemic plant infection determinants of barley yellow dwarf luteovirus-PAV are contained in the coat protein readthrough domain and 17-kDa protein, respectively. Virology 1996, 219, 57–65. [Google Scholar] [CrossRef]
  41. Schiltz, C.J.; Wilson, J.R.; Hosford, C.J.; Adams, M.C.; Preising, S.E.; DeBlasio, S.L.; MacLeod, H.J.; Van Eck, J.; Heck, M.L.; Chappie, J.S. Polerovirus N-terminal readthrough domain structures reveal molecular strategies for mitigating virus transmission by aphids. Nat. Commun. 2022, 13, 6368. [Google Scholar] [CrossRef] [PubMed]
  42. Moonan, F.; Molina, J.; Mirkov, T.E. Sugarcane yellow leaf virus: An emerging virus that has evolved by recombination between luteoviral and poleroviral ancestors. Virology 2000, 269, 156–171. [Google Scholar] [CrossRef]
  43. Zhao, K.; Yin, Y.; Hua, M.; Wang, S.; Mo, X.; Yuan, E.; Zheng, H.; Lin, L.; Chen, H.; Lu, Y.; et al. Pod pepper vein yellows virus, a new recombinant polerovirus infecting Capsicum frutescens in Yunnan province. China Virol. J. 2021, 18, 42. [Google Scholar] [CrossRef] [PubMed]
  44. LaTourrette, K.; Holste, N.M.; Garcia-Ruiz, H. Polerovirus genomic variation. Virus Evol. 2021, 7, veab102. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Symptoms of the diseased Viola philippica plants, RT-PCR detection of VPPV contigs and characteristics of the VPPV genome. (a,b) The symptomatic V. philippica plants show symptoms of leaf yellowing, mottling, and vein chlorosis. (c) Healthy V. philippica plants. (d) RT-PCR detection of contig 1345 and contig 21,591 in the symptomatic V. philippica plants and healthy plants (CK-). (e) Schematic diagram of the genome organization of VPPV. The positions of −1 ribosomal frameshift site in ORF1 and stop codon readthrough are marked in the schematic diagram.
Figure 1. Symptoms of the diseased Viola philippica plants, RT-PCR detection of VPPV contigs and characteristics of the VPPV genome. (a,b) The symptomatic V. philippica plants show symptoms of leaf yellowing, mottling, and vein chlorosis. (c) Healthy V. philippica plants. (d) RT-PCR detection of contig 1345 and contig 21,591 in the symptomatic V. philippica plants and healthy plants (CK-). (e) Schematic diagram of the genome organization of VPPV. The positions of −1 ribosomal frameshift site in ORF1 and stop codon readthrough are marked in the schematic diagram.
Viruses 17 00114 g001
Figure 2. Recombination and phylogenetic analysis of VPPV with other related poleroviruses in the family Solemoviridae and Tombusviridae. Results of the recombination analyzed by GENECONV method (a) and MaxChi method (b). Lines indicate the percentage of similarity per alignment, and the pink area indicates the recombinant part of the sequence. The phylogenetic trees were constructed based on RdRp (c) or CP protein (d) using the maximum likelihood method with 1000 bootstraps. VPPV is shown in red. The bootstrap values are indicated adjacent to the nodes. Detailed sequence information of these viruses is shown in Table S4.
Figure 2. Recombination and phylogenetic analysis of VPPV with other related poleroviruses in the family Solemoviridae and Tombusviridae. Results of the recombination analyzed by GENECONV method (a) and MaxChi method (b). Lines indicate the percentage of similarity per alignment, and the pink area indicates the recombinant part of the sequence. The phylogenetic trees were constructed based on RdRp (c) or CP protein (d) using the maximum likelihood method with 1000 bootstraps. VPPV is shown in red. The bootstrap values are indicated adjacent to the nodes. Detailed sequence information of these viruses is shown in Table S4.
Viruses 17 00114 g002
Figure 3. Agroinfiltration inoculation of the established VPPV infectious cDNA clone and identification of pathogenicity determinants of VPPV. (a) RT-PCR detection of VPPV in the systemic leaves of N. benthamiana plants infiltrated with A. tumefaciens culture carrying the pCB301-VPPV (line 1 and 2) or carrying the empty pCB301 vector (CK-). pCB301-VPPV plasmid was used as PCR positive control (CK+). (b) Symptoms of N. benthamiana plants infiltrated with A. tumefaciens culture carrying the pCB301-VPPV vector (left) or carrying the empty pCB301 vector (the CK- control, right). (c) Symptoms of N. benthamiana plants expressing VPPV P0, P1, RdRp, P3a, P3, P4, or RTP using PVX heterologous expression system. White arrows indicate necrotic patches. (d) Western blot analysis of PVX CP accumulation in systemic leaves of PVX inoculated N. benthamiana at 10 days post inoculation (dpi). The Rubisco stained by Ponceau S was used as the loading control.
Figure 3. Agroinfiltration inoculation of the established VPPV infectious cDNA clone and identification of pathogenicity determinants of VPPV. (a) RT-PCR detection of VPPV in the systemic leaves of N. benthamiana plants infiltrated with A. tumefaciens culture carrying the pCB301-VPPV (line 1 and 2) or carrying the empty pCB301 vector (CK-). pCB301-VPPV plasmid was used as PCR positive control (CK+). (b) Symptoms of N. benthamiana plants infiltrated with A. tumefaciens culture carrying the pCB301-VPPV vector (left) or carrying the empty pCB301 vector (the CK- control, right). (c) Symptoms of N. benthamiana plants expressing VPPV P0, P1, RdRp, P3a, P3, P4, or RTP using PVX heterologous expression system. White arrows indicate necrotic patches. (d) Western blot analysis of PVX CP accumulation in systemic leaves of PVX inoculated N. benthamiana at 10 days post inoculation (dpi). The Rubisco stained by Ponceau S was used as the loading control.
Viruses 17 00114 g003
Table 1. The nucleotide and amino acid identities between VPPV and the 12 representative viruses in the genus Polerovirus.
Table 1. The nucleotide and amino acid identities between VPPV and the 12 representative viruses in the genus Polerovirus.
Virus NameNucleotide Identities (%)Amino Acid Identities (%)
P0P1RdRpP3 (CP)P4 (MP)RTP
potato leafroll virus46.515.829.247.634.226.623.8
pepper vein yellows virus 445.324.235.153.532.315.021.5
pepper vein yellows virus 636.723.235.754.133.915.026.8
beet chlorosis virus47.919.829.649.535.924.627.2
carrot red leaf virus47.716.530.450.134.220.023.6
cereal yellow dwarf virus-RPS46.516.729.344.537.222.028.7
chickpea chlorotic stunt virus56.121.737.749.335.023.929.1
cotton leafroll dwarf virus50.519.236.852.335.523.130.0
maize yellow mosaic virus49.118.032.249.430.918.028.3
melon aphid-borne yellows virus51.823.339.156.235.722.629.6
pumpkin polerovirus48.923.226.344.234.121.629.1
turnip yellows virus50.822.137.554.236.425.427.8
Red numbers indicate the highest identity.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chen, Y.; Chen, G.; Yu, J.; Zhou, Y.; Fei, S.; Chen, H.; Wu, J.; Fu, S. Genomic, Evolutionary, and Pathogenic Characterization of a New Polerovirus in Traditional Chinese Medicine Viola philippica. Viruses 2025, 17, 114. https://doi.org/10.3390/v17010114

AMA Style

Chen Y, Chen G, Yu J, Zhou Y, Fei S, Chen H, Wu J, Fu S. Genomic, Evolutionary, and Pathogenic Characterization of a New Polerovirus in Traditional Chinese Medicine Viola philippica. Viruses. 2025; 17(1):114. https://doi.org/10.3390/v17010114

Chicago/Turabian Style

Chen, Yuanling, Gaoxiang Chen, Jiaping Yu, Yali Zhou, Shifang Fei, Haorong Chen, Jianxiang Wu, and Shuai Fu. 2025. "Genomic, Evolutionary, and Pathogenic Characterization of a New Polerovirus in Traditional Chinese Medicine Viola philippica" Viruses 17, no. 1: 114. https://doi.org/10.3390/v17010114

APA Style

Chen, Y., Chen, G., Yu, J., Zhou, Y., Fei, S., Chen, H., Wu, J., & Fu, S. (2025). Genomic, Evolutionary, and Pathogenic Characterization of a New Polerovirus in Traditional Chinese Medicine Viola philippica. Viruses, 17(1), 114. https://doi.org/10.3390/v17010114

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop