Pod pepper vein yellows virus, a new recombinant polerovirus infecting Capsicum frutescens in Yunnan province, China

Pepper vein yellows viruses (PeVYV) are phloem-restricted viruses in the genus Polerovirus, family Luteoviridae. Typical viral symptoms of PeVYV including interveinal yellowing of leaves and upward leaf curling were observed in pod pepper plants (Capsicum frutescens) growing in Wenshan city, Yunnan province, China. The complete genome sequence of a virus from a sample of these plants was determined by next-generation sequencing and RT-PCR. Pod pepper vein yellows virus (PoPeVYV) (MT188667) has a genome of 6015 nucleotides, and the characteristic genome organization of a member of the genus Polerovirus. In the 5′ half of its genome (encoding P0 to P4), PoPeVYV is most similar (93.1% nt identity) to PeVYV-3 (Pepper vein yellows virus 3) (KP326573) but diverges greatly in the 3′-part encoding P5, where it is most similar (91.7% nt identity) to tobacco vein distorting virus (TVDV, EF529624) suggesting a recombinant origin. Recombination analysis predicted a single recombination event affecting nucleotide positions 4126 to 5192 nt, with PeVYV-3 as the major parent but with the region 4126–5192 nt derived from TVDV as the minor parent. A full-length clone of PoPeVYV was constructed and shown to be infectious in C. frutescens by RT-PCR and the presence of icosahedral viral particles.


Background
Pepper vein yellows viruses (PeVYV) are phloemrestricted viruses in the genus Polerovirus, family Luteoviridae and are currently classified into six species (International Committee on Taxonomy of Viruses [ICTV] 2019 release), named Pepper vein yellows virus 1 to 6 [1][2][3][4][5][6]. They have 86.2-94.6% nucleotide identity between them. Polerovirus genomes have seven open reading frames (ORF0 to ORF5 and ORF3a), putatively encoding proteins P0 to P5 and P3a [6]. Recombination is an important source of genetic variability in viruses, particularly for viruses possessing an RNA genome. PeVYVs have quite close nucleotide similarities to tobacco vein distorting virus (TVDV) in the 5′ half of their genomes, and are considered to have arisen from recombination between TVDV and other poleroviruses [4,7,8]. We here report a new recombinant of PeVYV with high identity to PeVYV-3 in the 5′ half of its genome and to TVDV in the 3′ part.
Longshu village (a, 45 samples), Bai Shiyan village (b, 28 samples) and Da Shuduan village (c, 16 samples). All had typical viral symptoms of interveinal leaf yellowing and fruit discoloration (Additional file 2: Fig. S1). These leaves (0.4 g per sample) were mixed into three pooled samples by origin, and then sent to the Central Laboratory of Zhejiang Academy of Agricultural Sciences (Hangzhou, China) for Next-generation RNA-Seq sequencing (NGS). A total amount of 1 μg RNA per sample was used as input material for the RNA  In order to verify the virus sequences, total RNA was isolated from each sample using TRIzol ™ Reagent (Invitrogen) in compliance with the manufacturer's instructions. Reverse transcription (RT) polymerase chain reaction (PCR) analyses were performed using the ReverTra AceTM qPCR RT Master Mix (Toyobo) and KOD-plus-Neo (Toyobo) following the manufacturer's protocol. RT was performed at 42 °C for 60 min with random primers followed by 72 °C for 10 min. The cycling conditions for the subsequent PCR were: 98 °C 3 min, and then 35 cycles of 98 °C for 30 s, 55 °C for 30 s, 68 °C for 1 kb/min; and 68 °C for 10 min. Primers designed from Contig_62 were used to amplify a coat protein (CP) fragment of the virus (PeVYV-CP) (Additional file 1: Table S1). Fragments of the expected size (621 bp) were obtained in 58 of the 89 symptomatic samples but not from healthy plants (raised from seed in the laboratory). After adding dATP at the 3′-terminal of the PCR products using1 μL Taq ™ (TaKaRa) at 72 °C for 10 min, they were cloned into pGEM-T Easy vector (Promega) and sequenced commercially (Sangon Biotech). The sequences were aligned using MUSCLE, and pairwise nucleotide sequence comparisons were done using the SDT (Species Demarcation Tool) v1.2 program [9,10].
The 116 amplicons were 93.0-98.4% identical to PeVYV-3, and so we tentatively designated this virus in pod pepper as Pod Pepper vein yellows virus (PoPeVYV).
In order to acquire the full-length sequence of PoP-eVYV in pod pepper, 5′ and 3′ RACE reactions were then performed to obtain the complete terminal sequences. In brief, the 3′ end of the viral RNA was polyadenylated using Poly(A) Tailing Kit (Invitrogen) and first-strand cDNA was synthesized using M4T primer (Additional file 1: Table S1). The 3′ terminal end was PCR amplified from the cDNA using M4 and 3′-RACE-F (Additional file 1: Table S1). Similarly, the 5′ end of the cDNA, after purification by treatment with Gel extraction Kit (Omega) was polyadenylated using the Poly(A) Tailing Kit. The purified cDNA/RNA heterocomplex was ligated with ZHM1 using T4 RNA ligase (Thermo Scientific) and PCR amplified using ZHM2 and 5′ RACE-R (Additional file 1: Table S1; Fig. 1a, b). To avoid errors in sequence assembly, the whole viral sequence was then amplified using two overlapping sections with the primer pairs PoPeVYV-1 and PoPeVYV-2 (Additional file 1: Table S1; Fig. 1c). RT-PCR products (expected sizes 4614 bp and 2824 bp) were obtained using the methods of KOD-plus-Neo (Toyobo) and the two products (which overlapped by 1423 bp) were cloned into the pGEM-T Easy Vector (Promega) and sequenced. We obtained 14 amplicons of PoPeVYV-1 and 4 amplicons of PoPeVYV-2, and found that 12 amplicons of PoPeVYV-1 exactly matched the overlapping region of PoPeVYV-2. The complete sequence was 6015 nt long (GenBank accession number: MT188667).
Poleroviruses are prone to recombination among themselves or with viruses belonging to other genera and the relationships between PoPeVYV, the previously described PeVYV isolates and TVDV suggests that there has been a recombination event affecting the 3′-end of the genome. This was confirmed using a variety of methods on the RDP4 recombinant platform [11]. A single recombination event affecting nucleotide positions 4126-5952, with TVDV and PeVYV-3 as the respective minor and major parents was consistently identified using GENECONV (P value of 1.035 E-93), RDP (8.410 E-85), BootScan (7.726 E-80), MaxChi (4.294 E-35), Chimera (7.719 E-05), SiScan (3.341 E-68), and Phylpro (2.331 E−15) (Fig. 2). Alignment of amino acid sequences of the PoPeVYV proteins with those of the six species of PeVYVs and TVDV also indicated a recombination event (Table 1). Phylogenetic trees were constructed from the complete genome sequences, 5′-half or 3′-half of other poleroviruses and an enamovirus as an outgroup (Additional file 4: Fig. S3). PoPeVYV, PeVYVs and TVDV formed an independent group separate from the other poleroviruses. The analysis confirmed that the genome and 3′ half of PoPeVYV was similar to TVDV, but that the 5′ half of PoPeVYV was similar to the PeVYVs (Additional file 4: Fig. S3). This is reflected in the phylogenetic analysis of the amino acid sequences of the separate gene products [12]: the P0 of PoPeVYV is most similar to PeVYV-1/4, P1/P2/P3/P4 are closely related to other PeVYVs but the RTD is most similar to TVDV (Fig. 3). However, proteins encoded by the 5′ half of the genome (P0 to P4) had the highest identity to those of PeVYV-3 and PeVYV-6, whereas proteins translated from the 3′ half (P5/RTD) were more closely related to those of TVDV (Table 1, Fig. 3). The various PeVYVs recognized usually have distinct P0 protein sequences and the topology of the P0 phylogenetic tree may depict the ongoing selection for a protein exhibiting better RNA silencing suppression capacity after a recent host jump [4,7,13]. The P0 of PoPeVYV has 77.5-87.1% amino acid identity and 85.2-90.5% nucleotide identity to that of the other PeVYVs, and 76.3% amino acid identity (81.3% nucleotide identity) to that of TVDV (Table 1). In  the P0 phylogenetic tree PoPeVYV was most similar to PeVYV-1 and PeVYV-4 and was clearly more similar to the PeVYVs than to TVDV (Fig. 3). These results suggest that PoPeVYV is a new recombinant polerovirus. An infectious clone of the virus was constructed for further investigation. RT-PCR was performed using KOD-plus-Neo (Toyobo) following the manufacturer's protocol. To generate infectious clones, the CloneExpress MultiS One Step Cloning Kit (Vazyme) was used for homologous recombination. Two overlapped PCR products were amplified with primers (Inf-PoPeVYV-1/2), and recombined with the linearized binary vector pCB301-MD (a modified version of pCB301 [14]), which includes the double 35S promoter and nopaline synthase terminator (NOS) (Additional file 5: Fig. S4) [15,16]. In this way, the full-length PoPeVYV cDNA was inserted between an upstream 35S promoter and a downstream hepatitis delta virus (HDV) ribozyme and NOS terminator in the binary vector to construct pCB-PoPeVYV. The PoPeVYV cDNA was ligated with 3′-terminal of 35S promoter without any extra nucleotide (Additional file 6: Fig. S5) [15,16]. This clone was transformed into Agrobacterium tumefaciens which was then delivered to C. frutescens plantlets by infiltration. There was mild upward leaf curling 45 days after inoculation (Fig. 4a), and RT-PCR using primers to detect the coat protein gene in the newly-emerged non-inoculated leaves showed that viral RNA was present and had spread systemically in all the inoculated plants (12/12) but not in the controls (Fig. 4b). Virions were purified from C. frutescens leaves using the method described previously [17]. Isometric particles about 25 nm in diameter were observed in the purified preparation from the inoculated plants (Fig. 4c) but not from the controls. RT-PCR and subsequent sequencing confirmed that the virus in the symptomatic (systemic) leaves had the same sequence as the PoPeVYV cDNA clone inoculated (data not show). These results demonstrate the infectivity of the full-length PoPeVYV to C. frutescens.
In fields of cultivated pod pepper in Wenshan, we had observed severe yellowing and curling symptoms with small leaves (Additional file 2: Fig. S1). Those symptoms were more severe than we observed in plants inoculated with the infectious clone of PoPeVYV by Agrobacterium tumefaciens which also did not have the obvious vein  yellowing usually caused by PeVYVs. Previous studies have shown that the P0 proteins of poleroviruses display variable RNA silencing suppression activity, and induce distinct symptoms [18,19]. The P0 of PoPeVYV has only 77.5-87.1% amino acid identity to the other PeVYVs, which could account for the distinct symptoms but the field plants were also infected with other viruses (PeVYVs, ChiVMV, ChiRSV, CMV etc.) as frequently happens in the field and which has also been reported in Taiwan and Thailand [20,21]. The severe viral symptoms in the field may therefore be a synergistic effect of mixed infection.
Recombination is an important source of genetic variability in viruses, particularly for viruses possessing an RNA genome. PeVYVs have higher identities to TVDV in the 5′ half of the genome and are considered to be recombinants of TVDV and other poleroviruses [4,7,8]. The new recombinant identified here has higher identity to TVDV at the 3′ part of genome, indicating a different sort of recombinant event. The 5′ half of the PoPeVYV genome shares high identity with PeVYV-3, which was reported from Hunan province in China [3]. The 3′ half is much more distant from PeVYV-3 but is more homologous with that of TVDV. TVDV cause a devastating tobacco disease in many tobacco producing areas include Hunan, Guizhou and Yunnan provinces in China [22,23]. Additionally, PeYVY has been reported to infect tobacco in Guizhou province suggesting that co-infection of tobacco might have provided an opportunity for recombination [24]. Recombination poses a A 621-bp fragment of the p3 was amplified, with a plasmid containing the p3 as a positive control. c Virions purified from leaves infected by the PoPeVYV infectious clone and observed by TEM problem for classification. The currently-recommended species demarcation criteria in the family Luteoviridae suggest that different species should have > 10% difference in amino acid sequence identity in any gene product from their closest relative. The P0 and P5 proteins of PoPeVYV have respectively 12.9-23.7% and 10.3-70.5% difference in amino acid identity to those of PeVYV1-6 ( Table 1). If applied here, those criteria suggest that PoPeVYV could be representative of a distinct species.

Conclusions
A full-length sequence of Pod pepper vein yellows virus (MT188667) was determined. Alignment and recombination analysis predicted a single recombination event with PeVYV-3 as the major parent but with the region 4126-5192 nt derived from TVDV as the minor parent. PoPeVYV is a new recombinant polerovirus infecting C. frutescens in Yunnan province, China.