Molecular Characterization of Two Isolates of Wild Tomato Mosaic Virus and Chilli Veinal Mottle Virus Co-infecting Chilli Pepper in China

The present study reports observation of a eld chilli pepper disease consisting of a co-infection with two potyviruses: Wild tomato mosaic virus Dehong isolate (WTMV-Dh) and Chili veinal mottle virus Dehong isolate (ChiVMV-Dh). We obtained the complete genome sequences of these two viruses by NGS sequencing. The WTMV-Dh is 9,598 nucleotides (nt) in length and encodes a complete polyprotein of 3,075 amino acids (aa). The polyprotein of WTMV-Dh shares 76.1–82.6% nt and 85.3–89.5% aa identities with the other three WTMV isolates reported previously. The ChiVMV-Dh is 9688 nt in length and encodes a complete polyprotein with 3, 089 aa. The polyprotein of ChiVMV-Dh shares 80.8–92.2% nt and 85.3–95.6% aa identities with the other ChiVMV isolates reported previously. Following phylogenetic analysis based on the polyprotein sequences of other potyviruses, WTMV-Dh clustered with the Vietnam strain WTMV-Laichau while ChiVMV-Dh clustered with several ChiVMV Sichuan isolates. Evaluation of the recombination events within the WTMV and ChiVMV subgroups indicated that some putative recombination events occurred in critical regions. These regions include the N-terminal of HC-Pro and P1 region of WTMV-Dh, CP and the P3 to CI region of ChiVMV-Dh, which may be new evidence of adaptive evolution of potyviruses.

comparison (MEGA 6.06) of WTMV-Dh to other potyviruses based on polyprotein amino acid sequence showed that the WTMV subgroup is divided into two branches. In the rst branch, WTMV-Dh clustered with the Vietnam strain Laichau together, while the second branch was composed of the other two Sichuan strains (Fig. 2).
Sequence analysis showed that the complete sequence of ChiVMV-Dh contains a complete polyprotein ORF but lacks 39 nt on the 3'UTR. The polyprotein is composed of 3,089 amino acids and is expected to cleave into 10 mature proteins, similar to other ChiVMV isolates [5] (S1 Fig). The coding region of the ChiVMV-Dh polyprotein is more closely related to Sichuan strains than Yunnan strains. The polyprotein shares 95.6%, 95.3%, and 95.3% aa sequence identity with the Sichuan strains Yp8 (KC711055), Pp4 (KC711056), and LZ (MK405594), respectively, and 90.4% aa sequence identity with the Yunnan strain YN-t (JX088636). But ChiVMV-Dh shared the highest identity with YN-t in the 5′UTR (80.2% nt) and P1 (84.4% nt, 82.7% aa) regions compared with other ChiVMV isolates ( Table 3). The 6K1 of ChiVMV-Dh shared a low nucleotide homology but the highest amino acid identity with YN-t (100% aa). Among all the ChiVMV isolates, the P3 of ChiVMV-Dh shared the least homology with the YN-t isolate (76.6% nt, 76.5% aa). The recombination events were examined by RDP4 based on 22 potyviruses sequences downloaded from the NCBI (the sequence information is shown in Fig. 2) that are closely correlated with ChiVMV-Dh. As shown in Table 4, a total of four putative recombination events related to ChiVMV-Dh were found. One of the putative recombination events was detected in CP and the 3′UTR, from nucleotide 8828 to 9781 of ChiVMV-Dh. The major parent was predicted as Dzh-Qyg, while YN-t was shown to bethe minor parent. The other three putative recombination events were detected in the same region of ChiVMV that covered parts of P3, CI, and 6K1. All the recombinants were from Indian and Pakistan isolates, and ChiVMV-Dh was assumed to be the minor parent, while the Korean isolate AM90971 was predicted as the major parent.
Phylogenetic analysis was conducted based on the polyprotein amino acid sequence to determine the phylogenetic relationship of ChiVMV-Dh with other potyviruses. The resultant neighbor-joining tree divided all ChiVMV isolates into two distinct clades, and ChiVMV-Dh clustered in the clade2 group. The ChiVMV clade2 has two branches, and ChiVMV-Dh clustered with Sichuan strains, while YN-t clustered in another branch (Fig. 2).
To evaluate the biological characteristics of the two viruses, Chenopodium quinoa was mechanically inoculated with the viruses using co-infected chilli pepper samples. Subsequently, single local lesions were separated three times to obtain puri ed virus isolates [19]. Puri ed virus isolates were identi ed as either ChiVMV-Dh or WTMV-Dh by RT-PCR ampli cation using ChiVMV and WTMV speci c primers (S2 Table) and sequencing, respectively. The ChiVMV-Dh and WTMV-Dh strains were then single or mix inoculated to chilli pepper using the methods above. Fourteen days post-inoculation (dpi), different symptoms were observed on chilli pepper, and RT-PCR was conducted to detect virus species (Fig. 3). The result showed WTMV-Dh induced mosaicism and deformity in developing leaves and yellowing of inoculated leaves. Meanwhile, ChiVMV-Dh induced shrinkage and yellowing of developing leaves, and veinal yellowing of inoculated leaves. Co-inoculated chilli pepper showed yellowing and mosaicism on developing leaves and yellowing and necrotic spotting on inoculated leaves.
In conclusion, we reported the co-infection of chilli pepper with WTMV and ChiVMV under natural conditions for the rst time. We also report the complete genome sequence of WTMV and ChiVMV from chilli pepper. By co-inoculation experiment, the symptoms of these two viruses on chilli pepper under singleinfection and co-infection were determined. It was found that the co-infection symptoms of WTMV and ChiVMV on chilli pepper are a combination of the single-infection symptoms and show more severe symptoms than single-infection. This nding indicates that co-infection of these two viruses may play a synergistic effect and increase virulence.
Recombination events are considered to be a signi cant source of genetic diversity and enhance host adaptability of plant viruses. Co-infection facilitates genetic exchange and recombination between two viruses. Although we did not nd direct recombination events between WTMV-Dh and ChiVMV-Dh, some recombination events were found within the WTMV or ChiVMV subgroups. The recombination events include P1, P3 to CP, and HC-Pro region.
Previous reports have indicated the P1 region is related to host adaptation and de nesthe host range [20,21]. Another study found recombination events in the P1 region of WTMV-Sn, but the minor parents of WTMV-Sn are unknown [14]. Herein, analysis of a WTMV-Dh recombination event found that WTMV-Dh may be a putative minor parent of WTMV-Sn in the P1 region.
Also found were recombination events of ChiVMV-Dh from P3 to the CI region, as the minor parent of isolates originated from India and Pakistan. Previous studies have found that Indian and Pakistan isolates of ChiVMV are closely related to Chinese isolates. Thus, a transboundary movement of infected chilli seedlings or other host plants is speculated [22,23]. However, our study indicated that the ChiVMV isolates from India and Pakistan are more likely to be recombinant viruses and that ChiVMV-Dh may serve as an intermediary between East Asian and South Asian ChiVMV isolates.
HC-Pro is a multifunctional protein composed of three functional regions, and the N-terminal region is necessary for aphid transmission [24,25]. We blasted the complete aminoacid sequence of HC-Pro in WTMV-Dh and ChiVMV-Dh. The results showed that the HC-Pro N-terminal regions of these two viruses differ substantially, but the central and C-terminal regions are highly conserved. Overall, the HC-Pro regions of the two viruses share 59%, 79.5% and 93.63% sequence identity in the N-terminal (1-100 aa), central (101-300 aa) and C-terminal (301-457 aa) regions, respectively. Therefore, we speculated the differences in the HC-Pro N-terminal sequences of WTMV-Dh and ChiVMV-Dh might be involved in the a nity with aphid species.
We also found a recombination event in the HC-Pro N-terminal region of WTMV-Dh. The major parent was Laichau, while the minor parent was the Sichuan isolate XC-1. This result further indicates that the N-terminal of HC-Pro is a hot spot for recombination and variability, which may participate in the adaptive evolution of aphid transmission. Table 1 Analyses for read coverage of reads mapping to WTMV-Dh and ChiVMV-Dh in co-infected samples.

Coding region
ChiVMV-Dh reads number