Next Article in Journal
Side-By-Side Evaluation of Three Commercial ELISAs for the Quantification of SARS-CoV-2 IgG Antibodies
Next Article in Special Issue
Interspecies Recombination-Led Speciation of a Novel Geminivirus in Pakistan
Previous Article in Journal
I226R Protein of African Swine Fever Virus Is a Suppressor of Innate Antiviral Responses
Previous Article in Special Issue
Identification of a New Badnavirus in the Chinaberry (Melia azedarach) Tree and Establishment of a LAMP-LFD Assay for Its Rapid and Visual Detection
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 14
Issue 3
10.3390/v14030576
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Characteristics and Incidence of Apple Rubbery Wood Virus 2 and Citrus Virus A Infecting Pear Trees in China

by
Yanxiang Wang
1,2,†,
Ying Wang
1,†,
Guoping Wang
1,
Qingyu Li
3,
Zhe Zhang
1,
Liu Li
1,
Yuzhuo Lv
1,
Zuokun Yang
1,
Jiashu Guo
1 and
Ni Hong
1,4,*
1
Key Lab of Plant Pathology of Hubei Province, College of Plant Science and Technology, Huazhong Agricultural University, Wuhan 430070, China
2
Hubei Key Laboratory of Edible Wild Plants Conservation and Utilization, College of Life Sciences, Hubei Normal University, Huangshi 435002, China
3
Yantai Academy of Agricultural Science, Yantai 264000, China
4
Key Laboratory of Horticultural Crop (Fruit Trees) Biology and Germplasm Creation of the Ministry of Agriculture, Huazhong Agricultural University, Wuhan 430070, China
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Viruses 2022, 14(3), 576; https://doi.org/10.3390/v14030576
Submission received: 19 January 2022 / Revised: 28 February 2022 / Accepted: 5 March 2022 / Published: 11 March 2022
(This article belongs to the Special Issue Evolution, Ecology and Diversity of Plant Virus)
Browse Figures
Versions Notes

Abstract

:
Apple rubbery wood virus 2 (ARWV-2) and citrus virus A (CiVA) belong to a recently approved family Phenuiviridae in the order Bunyavirales and possess negative-sense single-stranded RNA genomes. In this study, the genome sequence of three ARWV-2 isolates (S17E2, LYC2, and LYXS) and a CiVA isolate (CiVA-P) infecting pear trees grown in China were characterized using high-throughput sequencing combined with conventional reverse-transcription PCR (RT-PCR) assays. The genome-wide nt sequence identities were above 93.6% among the ARWV-2 isolates and above 93% among CiVA isolates. Sequence comparisons showed that sequence diversity occurred in the 5′ untranslated region of the ARWV-2 genome and the intergenic region of the CiVA genome. For the first time, this study revealed that ARWV-2 proteins Ma and Mb displayed a plasmodesma subcellular localization, and the MP of CiVA locates in cell periphery and can interact with the viral NP in bimolecular fluorescence complementation assays. RT-PCR tests disclosed that ARWV-2 widely occurs, while CiVA has a low incidence in pear trees grown in China. This study presents the first complete genome sequences and incidences of ARWV-2 and CiVA from pear trees and the obtained results extend our knowledge of the viral pathogens of pear grown in China.

1. Introduction

Apple rubbery wood virus 1 (ARWV-1) together with apple rubbery wood virus 2 (ARWV-2) were discovered in apple trees showing rubbery wood disease symptoms in Germany and USA [1]. Both ARWV-1 and ARWV-2 possess a tripartite negative-sense single-stranded RNA (-ssRNA) genome and are associated with apple rubbery wood disease [1]. The complementary RNA (cRNA) of each genomic RNA contains one open reading frame (ORF), with RNAs1–3 (large (L), medium (M), and small (S)) encoding an RNA-dependent RNA polymerase (RdRp), a movement protein (MP), and a nucleocapsid protein (NP), respectively. The major difference between ARWV-1 and ARWV-2 is that their encoded proteins have identities ranging from 59% (RdRp proteins) to ~66–68% (NP proteins). In addition, a few isolates of ARWV-2 have two distinct M segments and two distinct S RNA segments, referred as ARWV 2 Ma, 2 Mb, 2 Sa, and 2 Sb. The 2 Ma and 2 Mb segments have nt sequence identity of ~66%, and 2 Sa and 2 Sb segments have nt sequence identity of ~56% [1]. According to the genomic structures, phylogenetic relationships and coding protein features with those of viruses in the Bunyavirales order, the two viruses with species names Apple rubodvirus 1 and Apple rubodvirus 2 are assigned to the genus Rubodvirus in the recently approved family Phenuiviridae in the Bunyavirales order. In addition, the genus Rubodvirus contains two other virus species Grapevine rubodvirus 1 and Grapevine rubodvirus 2 identified from grapevine [2]. Recently, ARWV-1 and ARWV-2 have been identified in apple grown in China [3,4], and ARWV-2 was reported to be able to infect pear in China [5] and apple in Italy [6]. ARWV-1 and ARWV-2 together with citrus concave gum-associated virus (CCGaV) were also identified in apple-decline-afflicted samples in the United States and apple trees from a global apple collection [7].
Citrus virus A (CiVA) was described in a non-symptomatic, field-grown sweet orange tree in southern Italy [8]. The CiVA genome consists of two -ssRNA segments, with a large RNA (RNA1) encoding a RdRp in its complementary RNA and an ambisense RNA (RNA2) encoding MP and NP. The species name Coguvirus eburi has been established for CiVA, and the virus belongs to genus Coguvirus in the recently approved family Phenuiviridae [9]. The genus Coguvirus contains three species, including the type species Citrus coguvirus (the species name of CCGaV), Coguvirus eburi, and Grapevine coguvirus (the species name of grapevine associated cogu-like virus 1, GaCLV-1) [2,10]. Now, CiVA has been identified in pear in France [11], and is also found to infect pear and citrus in South Africa [12,13]. A recent study showed that the virus was related to the impietratura symptoms of citrus by using high-throughput sequencing (HTS)for identification of the impietratura infectious agent along with a comparative analysis among symptomatic and asymptomatic orange fruits in Greece [14].
Pear (Pyrus spp.) is widely grown in China. During field investigations, viral-disease-like symptoms are frequently observed on pear trees grown in China. In a recent study, we molecularly characterized a pathogenic emaravirus named pear chlorotic leaf spot associated virus (PCLSaV), which is associated with pear chlorotic leaf spot (PCLS) disease [15]. Meanwhile, leaf mosaic, chlorotic rings, and vein yellowing diseases also occur on some pear trees. To investigate the virome in the viral-disease-affected pear trees, leaves showing mosaic or chlorotic symptoms were collected from four pear trees, and subjected to RNA sequencing (RNA-Seq) and bioinformatic analyses. In addition to PCLSaV, two other -ssRNA viruses ARWV-2 and CiVA were identified. This is the first finding of these two viruses infecting pear grown in China. Here, we report the complete genome sequences of ARWV-2 and CiVA infecting pear by using high throughput sequencing combined with conventional Sanger sequencing and their infection status in pear trees grown in China. In addition, for the first time the subcellular localization of the predicated movement proteins encoded by the two viruses was determined.

2. Materials and Methods

2.1. RNA Sequencing

A leaf sample (ID: S17E2) from a P. pyrifolia cv. Cuiguan tree, from which ARWV-2 was identified, as reported previously [5], was included in this analysis. Samples LYC2 and LYXS were individually collected from a P. bretschneideri cv. Chili tree and a P. bretschneideri cv. Xiangshui tree grown in the Laiyang area, Shandong Province, China. The leaf samples FJCG and FJHH were individually collected from a P. pyrifolia cv. Cuiguan and a P. pyrifolia cv. Huanghua grown in Fujian Province, China. The three trees S17E2, LYC2, and LYXS exhibited brilliant yellow mosaic symptoms and yellow rings on some leaves of LYXS, the FJCG plant exhibited symptoms of chlorotic mottle and vein clearing, and the FJHH plant exhibited no obvious symptoms (Figure 1). These five samples were used for RNA-Seq analyses as described previously [15]. Briefly, ribosomal RNA (rRNA) in total RNA extracts was removed using an Epicentre Ribo-ZeroTM rRNA removal kit (Epicentre, Madison, WI, USA). The prepared rRNA-depleted RNA sample was used to construct a cDNA library with a TruSeq RNA Sample Prep Kit v2 (Illumina, San Diego, CA, USA) and sequenced on an Illumina HiSeq XTen sequencing machine (Illumina, San Diego, CA, USA) with a paired-end 150 bp setup (Biomarker Biology Technology Ltd. Company, Beijing, China). Analysis of the sequence data was performed with CLC Genomics Workbench V.10.1.1 (QIAGEN, Hilden, Germany).
The raw RNA reads from the Illumina platform were trimmed of adaptor sequences and filtered for low-quality reads using FASTP version 1.5.6 to remove adapter sequences and reads with more than 5% Ns or with 20% base quality values (Q20) less than 20. Then, the obtained RNA reads from each sample were de novo assembled into larger contigs using Velvet version 1.2.08 [16] with a k-mer of 15–17 and IDBA-UD version 1.1.1 [17] with k-mer values of 80, 90, and 110. Contigs were subsequently screened for sequence identities against the NCBI databases (http://www.ncbi.nlm.nih.gov, accessed on 16 October 2019) using BlastX and BlastN programs.

2.2. Amplification of the Viral Genomic RNAs

For the amplification of the CiVA genomic RNAs, cDNA templates were generated using random primer pd(N)6. Specific primers (Table S1) were designed based on the sequences of assembled contigs and their positions mapped to reported CiVA genomic RNAs (accession numbers. MG764565 and MG764566). For the amplification of the ARWV-2 genomic RNAs, the gaps in the RNA-Seq-derived contig sequences were filled by nested RT-PCR amplifications using primers designed based on contig sequences and reported genomic RNA sequences of ARWV-2 (MF062139–MF062143). For RNA1 amplification, total RNAs from ARWV-2-infected samples were reverse-transcribed using the random primer pd(N)6. For the amplification of M and S RNAs, cDNA templates were generated using primers targeting the 5′ termini of the viral RNAs (Table S1).
The terminal regions of the genomic RNAs of CiVA and ARWV-2 were determined by the RACE strategy using a commercial kit (GeneRacer, Invitrogen Carlsbad, CA USA) according to the manufacturer’s instructions. For 3′ RACE, poly (A) tails were added to the 3′ ends of the total RNAs by using the poly (A) polymerase kit (TaKaRa, Dalian, China), and cDNA was generated using the oligo (dT) primer provided in the 3′ RACE kit. For 5′ RACE, specific reverse primers were designed based on the nucleotide sequences conserved at the 5′ termini of the genomic RNAs of each virus (Table S1).
RT-PCR solutions and conditions were similar to those reported previously [15], except that annealing temperature and extension time varied depending on the primer sets used in each reaction and the size of the expected PCR products. PCR products were gel-purified and ligated into the pMD18-T vector (TaKaRa, Dalian, China). At least three positive clones of each PCR product were sequenced at Shanghai Sangon Biological Engineering & Technology and Service Co. Ltd., Shanghai, China. The obtained sequences were assembled into contiguous sequences by overlapping common regions of the amplicons.

2.3. Sequence Analyses

Prediction of ORFs was performed using the NCBI ORFFinder program (https://www.ncbi.nlm.nih.gov/orffinder, accessed on 12 July 2021) with minimal ORF length of 75 nt. Multiple sequence alignments and identity analyses were performed using the Muscle method in MEGA 7.0 [18]. Phylogenetic trees were constructed using the neighbor-joining method with 1000 bootstrap replicates. The corresponding sequences of an ARWV-1 isolate (accession numbers: MF062130, MF062131, and MF062132) and a CCGaV isolate (accession numbers. KX960112 and KX960111) were used as outgroups in the ARWV-2 and CiVA sequence-based trees. The numbers of mapped reads in virus genomic RNAs were measured by Samtools version 1.5 [15] with the default parameters. The read coverage was calculated using software Fastv version 1.0 [19] with the default parameters. The conserved protein domains were identified using the Pfam database [20] with the default parameters.

2.4. RT-PCR for Virus Detection

To understand the infection status of ARWV-2, leaf samples were collected from 173 pear trees grown in five provinces in China. Of these samples, 99 samples showed chlorotic leaf spots or yellow mosaic symptoms and 74 samples were asymptomatic. Total RNA extraction and reverse transcription (RT) were done as described above. For the efficient detection of ARWV-2, a nested RT-PCR (RT-nPCR) method was developed by using SaF1/SaR1 as an outer primer set and NP1-F/NP1-R as an inner primer set. The two sets of primers were designed in the sequence of ARWV-2 Sa segment (Table S1). The RT-nPCR solutions and conditions were similar to those reported previously [21]. For the detection of CiVA, a primer set R2-F/R2-R was designed from the sequence of CiVA ORF2b (Table S1). In these tests, a leaf sample of a virus-free seedling of P. betuleafolia was used as a negative control. PCR products were separated by electrophoresis on 2% agarose gels, stained with ethidium bromide, and visualized under UV light.

2.5. Subcellular Localization and Bimolecular Fluorescence Complementation Analysis

Two ORFs (without stop codons) coding for Ma and Mb of ARWV-2 isolate LYXS, and ORFs (without stop codons) coding for NP and MP of CiVA were amplified using gene-specific primers flanked with an attB recombination sequence (Table S2). The construction of vectors used for protein localization and bimolecular fluorescence complementation (BiFC) assays were done as described previously [22].
Recombinant plasmids were individually transformed into Agrobacterium tumefaciens strain GV3101 (Weidi Bio, Shanghai, China), and agro-infiltrated into leaves of Nicotiana benthamiana plants (5 weeks old). The CMV3a-mCherry and H2B-mCherry vectors were used as plasmodesma (PD) and nuclear markers, respectively. Fluorescence signals of fusion proteins transiently expressed in infiltrated leaves of N. benthamiana plants were viewed at 2 days post infiltration (dpi) using confocal laser scanning microscopy (CLSM; TCS-SP8, Leica Microsystems, Heidelberg, Germany) with an HC PL APO CS2 63×/1.20 WATER objective.

3. Results

3.1. RNA-Seq of ARWV-2 and CiVA

Totally, 87,218,776, 69,655,848, 76,050,182, 83,312,860, and 124,323,162 clean reads were obtained from samples S17E2, LYC2, LYXS, FJCG, and FJHH, respectively. BlastN and BlastX searches against the NCBI database using assembled contigs revealed plant viruses presenting in these samples. From samples S17E2 and LYC2, three contigs with lengths of 381–1499 bp and six contigs with lengths of 261–475 bp matched sequences of ARWV-2 genomic RNAs named L, Ma, and Sa. From sample LYXS, six contigs with lengths of 282–7367 bp matched sequences of genomic RNA segments L, Ma, Mb, Sa, and Sb of ARWV-2 isolates 982-11, R7, and R12. From sample FJCG, two contigs with lengths of 6703 bp and 2765 bp covering near 95% of CiVA genomic RNA1 and RNA2 were identified. From sample FJHH, five contigs with lengths of 406–4597 bp matched sequences of CiVA genomic RNA1 and RNA2. Meanwhile, two well-documented viruses, apple stem pitting virus (ASPV) and apple chlorotic leafspot virus (ACLSV), were identified in the sample S17E2, and apple stem grooving virus (ASGV) and ASPV were identified in the samples LYC2 and LYXS. The specimens FJCG and FJHH were infected with ASPV, ASGV, and ACLSV. Here, we considered only the two -ssRNA viruses ARWV-2 and CiVA for further analyzes. The identified contigs of ARWV-2 and CiVA are presented in Table 1.
Complete and near-complete sequences of ARWV-2 genomic RNAs from the three samples S17E2, LYC2, and LYXS, and complete genomic sequences (GenBank accession numbers: MZ819702–MZ819703) of CiVA from the sample FJCG were reconstructed by Sanger sequencing of RT-PCR products. The reconstructed genomic sequences showed over 99% identities with corresponding contig sequences generated from RNA assembly, indicating that the assembled contig sequences from RNA-Seq data were reliable.

3.2. Genomic Characteristics and Sequence Diversity of ARWV-2 Isolates Infecting Pear

The genome structures of ARWV-2 determined from samples S17E2 and LYC2, here named isolates S17E2 and LYC2, were identical to those of reported isolates BR-Gala, 355-1, and R7, and consisted of three RNA segments L, Ma, and Sa. The genome structure of ARWV-2 from sample LYXS, here named isolate LYXS, was identical to that of reported isolate R12, and consisted of five RNA segments L, Ma, Mb, Sa, and Sb. Complete sequences of L, Ma, and Sa segments of isolate S17E2 were 7370 nt, 1594 nt, and 1486 nt (GenBank accession numbers: MN163133–MN163135). For isolate LYC2, a near full-length L sequence (6736 nt, 314–7049 nt), a full-length Ma sequence (1607 nt, excluding 19 nt at 5′terminal), and a near full-length Sa sequence (1387 nt, excluding 127 nt at 3′terminal) (GenBank accession numbers: MZ819704–MZ819706) were determined. For isolate LYXS, a near full-length L sequence (7366 nt, excluding 127 nt at 3′terminal) and Ma sequence (1547 nt, excluding 61 nt at 3′terminal), a full-length Mb sequence (1602 nt), a full-length Sa sequence (1510 nt), and a near full-length Sb sequence (1243 nt, excluding 75 nt at 3′terminal) (GenBank accession numbers: MZ819707–MZ819711) were determined. The obtained sequences, together with the corresponding sequences of reported ARWV-2 isolates, were comparatively analyzed (Table 2). The genome-wide nt sequence identities among the ARWV-2 isolates were about 93.6–99.5%, 97.7–99.3%, 97.7–98.1%, 96.4–99.5%, and 97.1% for their ORFs in the RNA segments L, Ma, Mb, Sa, and Sb, respectively. Although the size of each segment was variable among isolates, except for the size variation occurring in Sa-ORF, other ORF sizes were the same among ARWV-2 isolates. Each ORF of isolate LYXS had the same size as the corresponding ORF of isolate R12. The size difference of each segment among these isolates was a result of length variability of the 5′untranslated region (UTR). The L segment of isolates LYXS, 355-1, R7, 982-11, and BR-Gala had a 30-nt insert as compared with isolate R12 (Figure S1). There were also different inserts or deletions in the 5′UTRs of other RNA segments of these isolates (Figure S1). Although both isolates LYXS and R12 had five segments, the inserts or deletions in the 5′UTRs of each RNA segment of the two isolates were not conserved. Most of the inserts occurred at the “A” enriched regions (Figure S1).
Neighbor-joining (NJ) trees constructed using the nucleotide sequences of the viral RNA segment L showed that the three ARWV-2 isolates from Chinese pear samples clustered into a large clade with the previously reported seven ARWV-2 isolates, and one isolate H2803 (host: apple) from the USA formed a separated clade. In the trees of segments M and S, the three Chinese ARWV-2 isolates clustered into two clades, with Mb and Sb segments of LYXS separated from M and S segments of S17E2 and LYC2, respectively (Figure 2).

3.3. Genomic Characteristics of CiVA

By RT-PCR amplifications and 3′ RACE reactions, the full-length sequences of the viral RNA1 and RNA2 of CiVA from pear sample FJCG (here named as isolate CiVA-P) were determined. Moreover, an expected amplicon with a size 1973 bp was obtained by using a primer pair N155F/M240R designed based on the sequences of the covering partial sequences of the two ORFs in the viral RNA2. The result confirmed that the ORF2a and ORF2b were encoded in the same RNA strand. The nt sequences of CiVA-P RNA1 and RNA2 obtained by Sanger sequencing were 99.9% and 99.6% identical to the corresponding contig sequences generated by de novo assembly, respectively. The overall nt identities between the isolate CiVA-P from a Chinese pear sample and each of the nine isolates available in GenBank were about 95% for their RNA1 and 93–94% for their RNA2. Among these CiVA isolates, the intergenic region (IGR) between ORF2a and ORF2b ranged from 301 nt to 313 nt, and CiVA-P determined here showed relatively low nt identities with other isolates reported from other countries, ranging from 83.9% to 88.0% (Table 3). The IGR had a high AU content of 76% and could form a secondary structure with a stop signal CUCUGCU conserved in phleboviruses [23].
The CiVA-FJCG RNA1 was 6690 nt in length with one ORF (ORF1, 6640–86 nt) in its complementary strand. The ORF1 encoded a putative RdRp consisting of 2184 aa with a molecular mass of 251 kDa, which contained six motifs conserved in the RdRp of viruses in the Bunyaviriles [24]. ORF1 of CiVA-P isolate shared 94.6–95.3% nt and 95.9–96.9% aa identities with the corresponding sequences of other CiVA isolates (Table 3). The RNA2 of CiVA-P was 2739 nt in length and contained two ORFs in opposite directions. ORF2a (57–1224 nt) encodes a MP of 395 aa with a molecular mass of 44.4 kDa. ORF2b (2667–1555 nt) encodes a NP of 370 aa with a molecular mass of 41.8 kDa. Its ORF2a shared 94.3–95.6% nt and 93.5–98.2% aa identities with the corresponding sequences of other CiVA isolates (Table 3). Its OR2b showed 94.0–95.0% nt and 93.0–95.7% aa identities with the corresponding sequences of other CiVA isolates.
In the NJ trees constructed using the nucleotide sequences of the viral RNA1, RNA2, and IGR, two Chinses isolates CiVA-P and CiVA-HH were always separated from other-CiVA isolates reported from other countries, and formed a distinct subclade, indicating genetic diverse of CiVA isolates CiVA-P and CiVA-HH, identified from pear grown in China, from other reported CiVA isolates (Figure 3).

3.4. RNA Profiles of Viruses ARWV-2 and CiVA

In total, 251 reads, 77 reads, and 12,911 reads matching ARWV-2 genomic RNAs were derived from samples S17E2, LYC2, and LYXS and 90,617 CiVA RNA reads were derived from sample FJCG. The high RNA read productivities of ARWV-2 from LYXS and CiVA from FJCG made it possible to evaluate the distribution of RNA reads on each segment of the two viral isolates. It was found that ARWV-2 RNA reads deriving from sample LYXS covered the full lengths of the five RNA segments (L, Ma, Mb, Sa, and Sb) of the viral genome (Figure S2). On each RNA segment, the RNA reads distributed unevenly, with several peaks at different positions. In particular, a prevalent read peak appeared in the 5′UTR of the Ma segment. The RNA read depths also differed among the RNA segments, with relatively high depths occurring in Sb and Mb segments (Figure S2).
CiVA RNA reads derived from sample FJCG also covered the full lengths of the viral RNAs 1 and 2. The overall read abundance on the viral genomic RNA1 was much lower than that of the viral genomic RNA2 (Figure S2). RNA reads of CiVA separated almost evenly in the viral genomic RNA1 except for several hotspots, whilst, RNA reads on the viral genomic RNA2 had two prevalent read peaks separated by an IGR.

3.5. Subcellular Localization of ARWV-2 and CiVA Proteins in Planta

Some ARWV-2 isolates have two additional RNAs coding proteins Sb and Mb, which show over 60% aa sequence identity with Sa and Ma. It was postulated that both Ma and Mb might function in the virus movement [1]. Pfam analysis showed that both protein Ma (interval: 123–236, E-value: 3.8 × 10−12) and Mb (interval: 130–242, E-value: 8.11 × 10−7) of ARWV-2 isolate LYXS contained the plant viral MP family (pfam01107) domains (Figure S3). To have a primary look at the subcellular localizations of the two proteins, the Ma and Mb were transiently expressed in epidermal cells of N. benthamiana leaves by using an agrobacterium infiltration method. It was found that proteins Ma and Mb showed a similar localization profile with fluorescence signals in periphery and punctate spots along cell membrane, which were co-localized with the plasmodesma (PD) marker CMV-3a-mCherry (Figure 4A,B), indicating that the viral Ma and Mb had the typical subcellular localization features of the movement proteins of plant viruses.
Although a core 30K viral MP domain containing a LxD/N50-70G motif was recognized in the MP of CiVA [8], whether the MP has the location features of viral movement proteins is unknown. This study showed that the fusion protein MP-eYFP of CiVA distributed smoothly along cell membrane and periphery, while the PD marker CMV-3a-mCherry distributed as punctate spots (Figure 4C). Considering that the cell to cell movement of some plant viruses needs MP and structure protein(s), we further localized the viral NP and tested the interaction between the viral NP and MP by BiFC assays. Results showed that NP-eYFP located in periphery and formed aggregates (Figure 4D), which looked like viral replication complexes (VRC) [25]. BiFC results showed that proteins NP and MP interacted with each other in both orientations. The interaction signal was observed along the cell periphery (Figure 4E,F), which was similar to the distribution of the viral MP, but different from PD marker CMV-3a-mCherry distribution.

3.6. RT-PCR Assays for ARWV-2 and CiVA in Pear Trees

RT-nPCR assays showed that ARWV-2 presented in 34 symptomatic and 24 asymptomatic samples, accounting for 34.3% and 32.4% of the 173 tested samples, respectively. The virus was detected in P. communis, P. bretschneideri, P. pyrifolia, and some hybrids (Table 4). In particular, 6 of 10 tested P. communis samples were positive to ARWV-2, with an incidence of 60%. Leaf samples from 57 pear trees (including the sample FJCG and FJHH used for the RNA-Seq analysis) grown in Hubei, Sichuan, and Fujian provinces in China were subjected to RT-PCR assay for CiVA using primer set R2-F/R2-R (Table S1). The virus was detected in four samples, accounting for 7.0% of the tested samples (Table 4). Of the CiVA-positive samples, one showed leaf mottle and two showed chlorotic leaf spot symptoms.
The 486-bp amplicons of ARWV-2 from 19 samples were cloned and sequenced. It was found that these amplified sequences shared about 99% nt and aa similarities and clones from each sample had over 99% nt similarity. The result indicated low sequence variation of the NP gene among ARWV-2 isolates. These sequences and the corresponding sequences of the three HTS detected isolates LYC2, LYXS, and S17E2 clustered into three distinct phylogenetic clades, with clade I represented with S17E2, clade II represented with LYC2 and LYXS, and the clade III consisting three isolates from Hubei Province (Figure 5).

4. Discussion

In this study, HTS was used to identify viruses in four diseased pear samples and revealed the presence of two recently described -ssRNA viruses, ARWV-2 and CiVA, and also three common pear-infecting viruses in the samples. This is the first report of the complete genome sequences of the two viruses from pear trees and the first detection of CiVA in China. The presented results extend our knowledge of the molecular characteristics and geographical distribution of the two viruses. The coverages of RNA-Seq reads on ARWV-2 genomic RNA segments varied among samples. The near full-length sequences of five RNA segments of ARWV-2 were assembled from the LYXS sample, which showed very strong mosaic disease symptom, while from sample LYC2, the ARWV-2 RNA read coverage was very low, which resulted in short assembled contigs (less than 500 nt). From sample S17E2, near full-length sequences of Sa and Ma segments were assembled, but only one short contig of 148 bp of the viral L segment was identified. CiVA reads had high cover depth and near full-length sequences of CiVA genomic RNA1 and RNA2 were assembled from FJCG library. HTS technologies have been proved to be a powerful tool for the identification of known and novel viruses [26,27]. Since the quality of these libraries were reliable with QC20 and QC30 values of 94%–99%, the differences in viral RNA read coverages among samples could be caused by the virus titer or the quality of RNA-Seq libraries [28,29]. Interestingly, the RNA read profiles of Mb and Sb segments of ARWV-2 isolate LYXS were significantly higher than those of their homolog Ma and Sa segments. The biological functions of the redundancy proteins encoded in Ma and Mb, and Sa and Sb segments are still unknown [1]. The high distribution profiles of RNA reads on the viral Mb and Sb segments might be related to the genome expression strategy and involved in specific functions. Similarly, RNA segments encoding for redundancy proteins have been reported in several emaraviruses [15,22,30]. The RNA read peaks on the genomic RNA2 of CiVA are separated by an IGR, which is similar to the characteristics for some other viruses previously reported [31].
The overall genomic structures of ARWV-2 and CiVA determined from pear plants grown in China were the same as those of the corresponding viruses infecting apple [1] and citrus [11], as reported previously. The three pear isolates of ARWV-2 identified in this study showed over 97% genomic identity with each other and with reported apple isolates. One isolate LYXS had five RNA segments, which were identical to those of an apple isolate R12 [1]. Our results revealed that the size of each RNA segment was variable among ARWV-2 isolates as a result of their 5′-UTR variability. Except for sample LYXS, we did not detect the Mb and Sb segments of ARWV-2 from other pear samples, indicating a low occurrence frequency of the ARWV-2 Mb and Sb segments in pear grown in China. The complete genomic sequences of three CiVA isolates from citrus [8,12], two CiVA isolates from apple [14], and a near-completed genomic sequence of two CiVA isolates from pear have been reported recently [11,13]. Here, we reported the complete genomic sequence of a CiVA-P isolate from a pear tree grown in China. The overall nt identities between the CiVA-P and each of CiVA isolates available in GenBank were about 95% for their RNA1 and 93–94% for their RNA2. Among these isolates, IGR between ORF2a and ORF2b had a size of 309 nt or 310 nt, and showed nt sequence diversity ranging from 11.4% to 13.8%. The sequence variation might be involved in the adaption of the virus in different host plants [32,33].
ARWV-2 was originally described in apple trees with rubbery wood disease [1]. In recent studies, the virus was detected in apple trees showing decline disease in the United States, asymptomatic trees from a global apple collection [7], and apple trees grown in China [4]. The latent infection of an apple rubbery-wood-disease-associated agent in pear has been recognized for some time [34]. In the present study, RT-nPCR assays showed that ARWV-2 commonly occurred in the investigated P. communis, P. bretschneideri, P. pyrifolia, and pear hybrids grown in several Chinese provinces. Although ARWV-2 was first detected by HTS in three pear trees showing leaf mosaic symptom, the virus also presented in some asymptomatic trees as detected by using RT-nPCR. Our previous bioassays demonstrated the infectious nature of ARWV-2 in pear plants, but did not find the associated symptoms on pear A20 [5]. The leaf mosaic symptom of the three RNA-Seq-analyzed samples was different from the chlorotic leafspot disease caused by PCLSaV [15], but like the mosaic symptom in apple trees infected by apple necrotic mosaic virus (ApNMV) [35]. However, we did not identify either PCLSaV or ApNMV by HTS in the samples. Further studies are needed to evaluate viruses potentially associated with the observed leaf mosaic disease. CiVA was originally characterized from an asymptomatic orange leaf sample [8]. Recently, the virus was found to be associated with impietratura disease of citrus in Greece [14]. In this study, we found that CiVA occurred only in four P. pyrifolia trees grown in Fujian and Sichuan provinces, indicating a lower incidence. Although CiVA was identified in pear trees showing viral-disease-like symptoms (one leaf mottle and two chlorotic leafspot), we couldn’t conclude the association of the virus with pear disease since the most tested samples showing similar symptoms were negative to the virus. In addition, we did not find the co-infection of the two viruses.
Considering the close taxon position of apple and pear belonging to the family Rosaceae, and the common infection of previously characterized viruses and viroids, including ASGV, ASPV, ACLSV, and apple scar skin viroid (ASSVd) [36], together with the two -ssRNA viruses characterized here, it is reasonable to speculate that similar virus origins might have occurred in apple and pear. Vegetative propagation has been proposed to play prevalent roles for the transmission of viruses and other graft-transmissible agents of perennial fruit trees. The wide infection of ARWV-2 as well as the three common viruses ASGV, ASPV, and ACLSV and one common viroid ASSVd indicate a relatively long occurring history of these viruses. The low occurrence frequency of CiVA might be a result of its recent introduction.
The 30K family MPs can form MP–RNA complexes to facilitate cell to cell movement of viral genome RNA via the PD [1]. Here, for the first time, we showed that two proteins Ma and Mb of ARWV-2 isolate LYXS displayed PD location, further supporting that the proteins might co-function on the virus movement [37,38]. However, the predicated MP of CiVA showed a cytoplasm location, which was different from the typical PD location of plant virus movement proteins, but similar to the locations of the movement protein NSm of two orthotospoviruses, impatiens necrotic spot (INSV), and iris yellow spot (IYSV) [39,40]. Likely, the movement functions of some plant virus proteins might not always be correlated with PD location [41,42]. Similar to nucleocapsid proteins of several orthotospoviruses, including capsicum chlorosis virus (CaCV), tomato spotted wilt virus (TSWV), INSV, and IYSV [25,39,40,42], the CiVA protein NP formed VRC-like aggregates in the cytoplasm of agroinfiltrated cells. The NP–MP complexes of orthotospoviruses are involved in ribonucleoprotein (RNP) movement to adjacent cells [43]. Thus, the interaction of CiVA proteins NP and MP might be necessary for the virus movement.
In conclusion, ARWV-2 and CiVA were identified in pear trees grown in China using HTS. The genome of ARWV-2 pear isolates is divergent, with three RNA segments or five RNA segments. ARWV-2 occurs widely, but CiVA has low incidence in pear trees grown in China. ARWV-2 isolate LYXS codes two movement proteins Ma and Mb, which display PD localization. The MP of CiVA is located in cell periphery and can interact with the viral NP.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/v14030576/s1, Figure S1. Multiple alignment of nucleotide sequences at the 5′un-translation region (5′UTR) of RNA segments of ARWV-2 amplified from pear samples and the corresponding sequences of reported isolates. Each sequence was labeled with a sample ID. Figure S2. Profile distribution of RNA reads from LYXS and FJCG along genomic RNAs of apple rubbery wood virus (ARWV-2) and citrus virus A (CiVA), respectively. The open reading frames were marked with green bold boxes under each RNA segment of ARWV-2 and CiVA. Figure S3. Pfam analysis the conserved domains of the proteins Ma (a) and Mb (b) encoded by ARWV-2. Table S1: Primers used for the genomic sequence amplification and RT-PCR detection of ARWV-2 and CiVA. Table S2: Primers designed for amplification the apple rubbery woody ARWV-2 ORFs of Ma and Mb and CiVA ORFs of NP and MP.

Author Contributions

Data curation, Y.W. (Yanxiang Wang) and Z.Y.; funding acquisition, G.W., and N.H.; investigation, Y.W. (Ying Wang), G.W., J.G., Z.Z., Q.L., L.L. and N.H.; methodology, J.G., Y.W. (Ying Wang) and Y.L.; project administration, G.W. and N.H.; resources, Q.L.; supervision, N.H. and G.W.; writing—original draft, N.H. and Y.W. (Yanxiang Wang); writing—review and editing, G.W. and N.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the Key National Project (grant number 2019YFD10 01800); the earmarked fund for Pear Modern Agro-industry Technology Research System (grant number CARS-28-16).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in article and supplementary materials.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Rott, M.E.; Kesanakurti, P.; Berwarth, C.; Rast, H.; Boyes, I.; Phelan, J.; Jelkmann, W. Discovery of negative-sense RNA viruses in trees infected with apple rubbery wood disease by next-generation sequencing. Plant Dis. 2018, 102, 1254–1263. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Kuhn, J.H.; Adkins, S.; Agwanda, B.R.; Al Kubrusli, R.; Alkhovsky, S.V.; Amarasinghe, G.K.; Avsic-Zupanc, T.; Ayllon, M.A.; Bahl, J.; Balkema-Buschmann, A.; et al. 2021 Taxonomic update of phylum Negarnaviricota (Riboviria: Orthornavirae), including the large orders Bunyavirales and Mononegavirales. Arch. Virol. 2021, 166, 3513–3566. [Google Scholar] [CrossRef] [PubMed]
  3. Hu, G.; Dong, Y.; Zhang, Z.; Fan, X.; Ren, F.; Lu, X. First report of apple rubbery wood virus 1 in apple in China. Plant Dis. 2021, 105, 3770. [Google Scholar] [CrossRef] [PubMed]
  4. Hu, G.; Dong, Y.; Zhang, Z.; Fan, X.; Ren, F.; Lu, X. First report of apple rubbery wood virus 2 infection of apples in China. Plant Dis. 2020, 105, 519. [Google Scholar] [CrossRef] [PubMed]
  5. Wang, Y.; Wang, G.P.; Hong, N.; Wang, Y.X.; Yang, Z.K.; Guo, J.S.; Zhang, Z.; Li, L.; Li, Y.J.; Li, Q.Y.; et al. First report of apple rubbery wood virus 2 infecting pear (Pyrus spp.) in China. Plant Dis. 2019, 103, 3293. [Google Scholar] [CrossRef]
  6. Minutolo, M.; Cinque, M.; Chiumenti, M.; Di Serio, F.; Alioto, D.; Navarro, B. Identification and characterization of citrus concave gum-associated virus infecting citrus and apple trees by serological, molecular and high-throughput sequencing approaches. Plants 2021, 10, 2039. [Google Scholar] [CrossRef]
  7. Wright, A.A.; Szostek, S.A.; Beaver-Kanuya, E.; Harper, S.J. Diversity of three bunya-like viruses infecting apple. Arch. Virol. 2018, 163, 3339–3343. [Google Scholar] [CrossRef]
  8. Navarro, B.; Zicca, S.; Minutolo, M.; Saponari, M.; Alioto, D.; Di Serio, F. A negative-stranded RNA virus infecting citrus trees: The second member of a new genus within the order Bunyavirales. Front. Microbiol. 2018, 9, 2340. [Google Scholar] [CrossRef]
  9. Kuhn, J.H.; Adkins, S.; Alioto, D.; Alkhovsky, S.V.; Amarasinghe, G.K.; Anthony, S.J.; Avsic-Zupanc, T.; Ayllon, M.A.; Bahl, J.; Balkema-Buschmann, A.; et al. 2020 taxonomic update for phylum Negarnaviricota (Riboviria: Orthornavirae), including the large orders Bunyavirales and Mononegavirales. Arch. Virol. 2020, 165, 3023–3072. [Google Scholar] [CrossRef]
  10. Diaz-Lara, A.; Navarro, B.; Di Serio, F.; Stevens, K.; Hwang, M.S.; Kohl, J.; Vu, S.T.; Falk, B.W.; Golino, D.; Al Rwahnih, M. Two novel negative-sense RNA viruses infecting grapevine are members of a newly proposed genus within the family Phenuiviridae. Viruses 2019, 11, 685. [Google Scholar] [CrossRef] [Green Version]
  11. Svanella-Dumas, L.; Marais, A.; Depasse, F.; Faure, C.; Lefebre, M.; Brans, Y.; Castaing, J.; Latour, F.; Bechti, A.; Candresse, T. First report of citrus virus A (CiVA) infecting pear (Pyrus communis) in France. Plant Dis. 2019, 103, 2703. [Google Scholar] [CrossRef]
  12. Bester, R.; Karaan, M.; Cook, G.; Maree, H.J. First report of citrus virus A in citrus in South Africa. J. Citrus Pathol. 2021, 8, 8. [Google Scholar] [CrossRef]
  13. Bougard, K.; Maree, H.J.; Pietersen, G.; Meitz-Hopkins, J.C.; Bester, R. First report of Coguvirus eburi infecting pear (Pyrus communis) in South Africa. Plant Dis. 2021, 106, 772. [Google Scholar] [CrossRef] [PubMed]
  14. Beris, D.; Ioanna, M.; Vassilakos, N.; Theologidis, I.; Rampou, A.; Kektsidou, O.; Massart, S.; Varveri, C. Association of citrus virus A to citrus impietratura disease symptoms. Phytopathology 2021, 111, 1782–1789. [Google Scholar] [CrossRef] [PubMed]
  15. Liu, H.; Wang, G.; Yang, Z.; Wang, Y.; Zhang, Z.; Li, L.; Waqas, M.; Hong, N.; Liu, H.; Wang, G.; et al. Identification and characterization of a pear chlorotic leaf spot-associated virus, a novel emaravirus associated with a severe disease of pear trees in China. Plant Dis. 2020, 104, 2786–2798. [Google Scholar] [CrossRef] [PubMed]
  16. Zerbino, D.R.; Birney, E. Velvet: Algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008, 18, 821–829. [Google Scholar] [CrossRef] [Green Version]
  17. Peng, Y.; Leung, H.C.M.; Yiu, S.M.; Chin, F.Y.L. IDBA-UD: A de novo assembler for single-cell and metagenomic sequencing data with highly uneven depth. Bioinformatics 2012, 28, 1420–1428. [Google Scholar] [CrossRef] [Green Version]
  18. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [Green Version]
  19. Chen, S.; He, C.; Li, Y.; Li, Z.; Melancon, C.E. A computational toolset for rapid identification of SARS-CoV-2, other viruses and microorganisms from sequencing data. Brief Bioinform. 2021, 22, 924–935. [Google Scholar] [CrossRef]
  20. Mistry, J.; Chuguransky, S.; Williams, L.; Qureshi, M.; Salazar, G.A.; Sonnhammer, E.L.L.; Tosatto, S.C.E.; Paladin, L.; Raj, S.; Richardson, L.J.; et al. Pfam: The protein families database in 2021. Nucleic Acids Res. 2021, 49, D412–D419. [Google Scholar] [CrossRef]
  21. Yao, B.; Lu, Y.; Li, L.; Wang, Y.; Li, Q.; Wang, G.; Hong, N. A sensitive nested multiplex RT-PCR assay for the simultaneous detection of three common viruses infecting pear plants. J. Virol. Methods 2019, 263, 105–110. [Google Scholar] [CrossRef] [PubMed]
  22. Guo, J.S.; Wang, Y.X.; Wang, G.P.; Hong, J.; Yang, Z.K.; Bai, J.Y.; Hong, N. Molecular characteristics of jujube yellow mottle-associated virus infecting jujube (Ziziphus jujuba Mill.) grown at Aksu in Xinjiang of China. Viruses 2021, 13, 25. [Google Scholar] [CrossRef] [PubMed]
  23. Albarino, C.G.; Bird, B.H.; Nichol, S.T. A shared transcription termination signal on negative and ambisense RNA genome segments of Rift Valley fever, sandfly fever Sicilian, and Toscana viruses. J. Virol. 2007, 81, 5246–5256. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Reguera, J.; Weber, F.; Cusack, S. Bunyaviridae RNA polymerases (L-protein) have an N-terminal, influenza-like endonuclease domain, essential for viral cap-dependent transcription. PLoS Pathog. 2010, 6, e1001101. [Google Scholar] [CrossRef] [PubMed]
  25. Feng, Z.; Chen, X.; Bao, Y.; Dong, J.; Zhang, Z.; Tao, X. Nucleocapsid of tomato spotted wilt tospovirus forms mobile particles that traffic on an actin/endoplasmic reticulum network driven by myosin XI-K. New Phytol. 2013, 200, 1212–1224. [Google Scholar] [CrossRef]
  26. 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]
  27. Rott, M.; Xiang, Y.; Boyes, I.; Belton, M.; Saeed, H.; Kesanakurti, P.; Hayes, S.; Lawrence, T.; Birch, C.; Bhagwat, B.; et al. Application of next generation sequencing for diagnostic testing of tree fruit viruses and viroids. Plant Dis. 2017, 101, 1489–1499. [Google Scholar] [CrossRef] [Green Version]
  28. Maliogka, V.I.; Minafra, A.; Saldarelli, P.; Ruiz-Garcia, A.B.; Glasa, M.; Katis, N.; Olmos, A. Recent advances on detection and characterization of fruit tree viruses using high-throughput sequencing technologies. Viruses 2018, 10, 436. [Google Scholar] [CrossRef] [Green Version]
  29. Pecman, A.; Kutnjak, D.; Gutiérrez-Aguirre, I.; Adams, I.; Fox, A.; Boonham, N.; Ravnikar, M. Next generation sequencing for detection and discovery of plant viruses and viroids: Comparison of two approaches. Front. Microbiol. 2017, 8, 1998. [Google Scholar] [CrossRef] [Green Version]
  30. Zheng, Y.; Navarro, B.; Wang, G.; Wang, Y.; Yang, Z.; Xu, W.; Zhu, C.; Wang, L.; Serio, F.D.; Hong, N. Actinidia chlorotic ringspot-associated virus: A novel emaravirus infecting kiwifruit plants. Mol. Plant Pathol. 2017, 18, 569–581. [Google Scholar] [CrossRef]
  31. Cao, M.; Zhang, S.; Li, M.; Liu, Y.; Dong, P.; Li, S.; Kuang, M.; Li, R.; Zhou, Y. Discovery of four novel viruses associated with flower yellowing disease of green sichuan pepper (Zanthoxylum Armatum) by virome analysis. Viruses 2019, 11, 696. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Hassani-Mehraban, A.; Saaijer, J.; Peters, D.; Goldbach, R.; Kormelink, R. Molecular and biological comparison of two tomato yellow ring virus (TYRV) isolates: Challenging the Tospovirus species concept. Arch. Virol. 2007, 152, 85–96. [Google Scholar] [CrossRef] [PubMed]
  33. Heinze, C.; Letschert, B.; Hristova, D.; Yankulova, M.; Kauadjouor, O.; Willingmann, P.; Atanassov, A.; Adam, G. Variability of the N-protein and the intergenic region of the S RNA of tomato spotted wilt tospovirus (TSWV). New Microbiol. 2001, 24, 175–187. [Google Scholar] [CrossRef] [PubMed]
  34. Campbell, A.I. The Growth of Young Pear Trees after Elimination of Some Viruses by Heat Treatment; Long Ashton Research Station: England, UK, 1966; pp. 111–116. [Google Scholar]
  35. Noda, H.; Yamagishi, N.; Yaegashi, H.; Xing, F.; Xie, J.P.; Li, S.F.; Zhou, T.; Ito, T.; Yoshikawa, N. Apple necrotic mosaic virus, a novel ilarvirus from mosaic-diseased apple trees in Japan and China. J. Gen. Plant Pathol. 2017, 83, 83–90. [Google Scholar] [CrossRef]
  36. Hadidi, A.; Barba, M. Economic Impact of Pome and Stone Fruit Viruses and Viroids. In Virus and Virus-like Diseases of Pome and Stone Fruits; Hadidi, A., Barba, M., Candresse, T., Jelkmann, W., Eds.; APS Press: St. Paul, MN, USA, 2011; pp. 1–7. ISBN 978-0-89054-396-2. [Google Scholar]
  37. Yu, C.; Karlin, D.G.; Lu, Y.; Wright, K.; Chen, J.; MacFarlane, S. Experimental and bioinformatic evidence that raspberry leaf blotch emaravirus P4 is a movement protein of the 30K superfamily. J. Gen. Virol. 2013, 94, 2117–2128. [Google Scholar] [CrossRef]
  38. Tilsner, J.; Cowan, G.H.; Roberts, A.G.; Chapman, S.N.; Ziegler, A.; Savenkov, E.; Torrance, L. Plasmodesmal targeting and intercellular movement of potato mop-top pomovirus is mediated by a membrane anchored tyrosine-based motif on the lumenal side of the endoplasmic reticulum and the C-terminal transmembrane domain in the TGB3 movement protein. Virology 2010, 402, 41–51. [Google Scholar] [CrossRef] [Green Version]
  39. Dietzgen, R.G.; Martin, K.M.; Anderson, G.; Goodin, M.M. In planta localization and interactions of impatiens necrotic spot tospovirus proteins. J. Gen. Virol. 2012, 93, 2490–2495. [Google Scholar] [CrossRef] [Green Version]
  40. Tripathi, D.; Raikhy, G.; Goodin, M.M.; Dietzgen, R.G.; Pappu, H.R. In vivo localization of iris yellow spot tospovirus (Bunyaviridae)-encoded proteins and identification of interacting regions of nucleocapsid and movement proteins. PLoS ONE 2015, 10, e0118973. [Google Scholar] [CrossRef]
  41. Mann, K.S.; Bejerman, N.; Johnson, K.N.; Dietzgen, R.G. Cytorhabdovirus P3 genes encode 30K-like cell-to-cell movement proteins. Virology 2016, 489, 20–33. [Google Scholar] [CrossRef] [Green Version]
  42. Widana Gamage, S.M.K.; Dietzgen, R.G. Intracellular localization, interactions and functions of capsicum chlorosis virus proteins. Front. Microbiol. 2017, 8, 612. [Google Scholar] [CrossRef] [Green Version]
  43. Soellick, T.R.; Uhrig, J.F.; Bucher, G.L.; Kellmann, J.W.; Schreier, P.H. The movement protein NSm of tomato spotted wilt tospovirus (TSWV): RNA binding, interaction with the TSWV N protein, and identification of interacting plant proteins. Proc. Natl. Acad. Sci. USA 2000, 97, 2373–2378. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Viruses 14 00576 g001 550
Figure 1. Disease symptoms on pear leaves used for high-throughput sequencing. (A) Cuiyu (ID: S17E2), (B) Chili (ID: LYC2), (C) Xiangshui (ID: LYXS), (D) Cuiguan (ID: FJCG), and (E) Huanghua (ID: FJHH).
Figure 1. Disease symptoms on pear leaves used for high-throughput sequencing. (A) Cuiyu (ID: S17E2), (B) Chili (ID: LYC2), (C) Xiangshui (ID: LYXS), (D) Cuiguan (ID: FJCG), and (E) Huanghua (ID: FJHH).
Viruses 14 00576 g001
Viruses 14 00576 g002 550
Figure 2. Unrooted NJ phylogenetic trees generated from the nucleotide sequences of L (A), M (B) and S (C) RNA segments of ARWV-2. In the trees based on M (B) and S (C) segment sequences, Ma and Mb, and Sa and Sb sequences from pear LYXS and the two reported isolates were included in the assay. The corresponding sequence of a reported ARWV-1 isolate was used as an outgroup in each tree. The referred sequences are marked by their isolate names followed by GenBank accession numbers. Three ARWV-2 isolates determined by high throughput sequencing are highlighted by black dots.
Figure 2. Unrooted NJ phylogenetic trees generated from the nucleotide sequences of L (A), M (B) and S (C) RNA segments of ARWV-2. In the trees based on M (B) and S (C) segment sequences, Ma and Mb, and Sa and Sb sequences from pear LYXS and the two reported isolates were included in the assay. The corresponding sequence of a reported ARWV-1 isolate was used as an outgroup in each tree. The referred sequences are marked by their isolate names followed by GenBank accession numbers. Three ARWV-2 isolates determined by high throughput sequencing are highlighted by black dots.
Viruses 14 00576 g002
Viruses 14 00576 g003 550
Figure 3. Unrooted NJ phylogenetic trees generated from the nucleotide sequences of RNA1 (A), RNA2 (B) segments of CiVA. The intergenic region (IGR) (C) sequences of CiVA RNA2 were also used for phylogenetic analysis. The corresponding sequence of a reported CCGaVisolate was used as an outgroup in each tree. The referred sequences are marked by their isolate names followed by GenBank accession numbers. Two Chinese CiVA isolates determined by high throughput sequencing are highlighted by black dots.
Figure 3. Unrooted NJ phylogenetic trees generated from the nucleotide sequences of RNA1 (A), RNA2 (B) segments of CiVA. The intergenic region (IGR) (C) sequences of CiVA RNA2 were also used for phylogenetic analysis. The corresponding sequence of a reported CCGaVisolate was used as an outgroup in each tree. The referred sequences are marked by their isolate names followed by GenBank accession numbers. Two Chinese CiVA isolates determined by high throughput sequencing are highlighted by black dots.
Viruses 14 00576 g003
Viruses 14 00576 g004 550
Figure 4. Subcellular localization (A–D) and bimolecular fluorescence (BiFC) (E,F) assays of two movement proteins Ma and Mb of ARWV-2 and MP and NP of CiVA in epidermal cells of wild-type Nicotiana benthamiana leaves. The fusion proteins CMV3a-mCherry and H2B-mCherry were used as plasmodesmata (PD) and nuclear markers, respectively. Colocalization dots of Ma and Mb with CMV3a-mCherry at PD are indicated by arrow heads (A,B). The aggregated bodies formed by protein NP of CiVA are denoted by arrows. The images were acquired 2 days after agroinfiltration under a confocal microscope at 63x/1.20 WATER objective. Scale bar = 20 μm.
Figure 4. Subcellular localization (A–D) and bimolecular fluorescence (BiFC) (E,F) assays of two movement proteins Ma and Mb of ARWV-2 and MP and NP of CiVA in epidermal cells of wild-type Nicotiana benthamiana leaves. The fusion proteins CMV3a-mCherry and H2B-mCherry were used as plasmodesmata (PD) and nuclear markers, respectively. Colocalization dots of Ma and Mb with CMV3a-mCherry at PD are indicated by arrow heads (A,B). The aggregated bodies formed by protein NP of CiVA are denoted by arrows. The images were acquired 2 days after agroinfiltration under a confocal microscope at 63x/1.20 WATER objective. Scale bar = 20 μm.
Viruses 14 00576 g004
Viruses 14 00576 g005 550
Figure 5. Unrooted NJ phylogenetic trees generated from the nucleotide sequences of partial S segment of ARWV-2. The corresponding sequences of three isolates S17E2, LYC2, and LYXS were included in the assay.
Figure 5. Unrooted NJ phylogenetic trees generated from the nucleotide sequences of partial S segment of ARWV-2. The corresponding sequences of three isolates S17E2, LYC2, and LYXS were included in the assay.
Viruses 14 00576 g005
Table
Table 1. Blast analysis of contigs derived by RNA-seq from five pear samples.
Table 1. Blast analysis of contigs derived by RNA-seq from five pear samples.
SampleRNAContig (bp)Matched Sequencentaa
Site%Site%
S17E2L381 aNC_055534.12224–2373, 2457–260497.0721–771, 791–847100
Ma1430MF062145.114–143197.11–37099.7
Sa1499MF062146.15–123796.91–37098.4
LYC2L261MW322828.17056–731597.32331–2376100
321MF062133.15732–601299.11884–198999.1
326MW322828.15366–569199.41768–186599.1
Sa326MF062146.1444–76998.5126–23398.2
288MK936228.1156–44399.334–104100
Ma475NC_055536.1817–129196.4253–37094.9
LYXSL7367MF062128.11–736698.31–237698.2
282MF062128.13466–37461001135–1228100
Sb1316NC_055533.116–131698.01–28698.6
Sa1511MF062146.11–148896.81–28898.6
Mb1611NC_055537.11–160897.91–39398.5
Ma1461MF062145.11–143796.01–37098.4
FJCGRNA16703MT922052.114–667995.31–218496.8
RNA22765MT720884.14–273694.41–395 (MP)97.9
1–370 (NP)94.0
FJHHRNA1661MT922052.11–56394.1202–218494.9
1722MT922052.1456–218695.21486–205896.3
4597MK273077.12088–668495.645–459597.6
RNA22562MT720884.19–253094.11–395 (MP)97.6
47–370 (NP)94.1
406MK273078.11888–219496.01–7895.0
a In the contig, 84 nucleotides between 2373 and 2457 nt were absent.
Table
Table 2. RNA sequence comparison of ARWV-2 isolates with that of an isolate LYXS determined in this study.
Table 2. RNA sequence comparison of ARWV-2 isolates with that of an isolate LYXS determined in this study.
HostIsolate aORF1 (L)ORF2a (Ma)ORF2b (Mb)ORF3a (Sa)ORF3b (Sb)
ntnt%ntnt%ntnt%ntnt%ntnt%
PearLYXS7367/1134 1182/867/855/
S17E2737097.6113498//86796.4//
LYC2673699.5113499.3//810 (partial)99.5//
AppleR12734097.4113497.8118298.186798.686197.1
BR-Gala734993.6113497.7//86798.7//
355-1738197.4//118297.786798.3//
R7736997.3113497.4//86798.7//
a Isolates determined in this study are in bold.
Table
Table 3. RNA sequence comparison of CiVA from two Chinese pear samples with that of reported CiVA isolates.
Table 3. RNA sequence comparison of CiVA from two Chinese pear samples with that of reported CiVA isolates.
HostIsolate aORF1ORF2aIGRORF2b
ntNt %aa %ntnt %aa %ntnt %ntnt %aa %
PearFJCG6555--1188--310 1113--
HH///111895.698.231398.7///
P1492655595.396.8/94.393.531085.3111395.095.4
CitrusCG301655595.196.7118895.597.030984.4111394.494.3
TXNO655594.996.5118895.095.730988.0111394.193.5
W4655594.795.9118894.695.431086.6111394.093.8
1.8655595.396.7118894.895.630184.1111993.593.0
1.5655595.096.3118894.994.930987.4111394.794.6
IMP-GR655594.696.3118895.296.730383.9111394.295.4
A7-Gr655595.096.9118895.296.730387.3111394.395.7
A3-Gr655594.996.8118895.096.530487.3111394.395.4
a The isolate determined in this study is in bold.
Table
Table 4. Incidence of ARWV-2 and CiVA in the pear leaf samples of five provinces of China as detected by RT-PCR.
Table 4. Incidence of ARWV-2 and CiVA in the pear leaf samples of five provinces of China as detected by RT-PCR.
VirusSpeciesNo. of SamplesInfected/Tested
AsymptomaticSymptomatic
ARWV-2P. communis105/91/1
P. bretschneideri218/191/2
P. pyrifolia11611/4525/71
hybrids220/17/21
unknown400/4
total17324/7434/99
CiVAP. pyrifolia2503/25
hybrids000
unknown2200
total5700
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, Y.; Wang, Y.; Wang, G.; Li, Q.; Zhang, Z.; Li, L.; Lv, Y.; Yang, Z.; Guo, J.; Hong, N. Molecular Characteristics and Incidence of Apple Rubbery Wood Virus 2 and Citrus Virus A Infecting Pear Trees in China. Viruses 2022, 14, 576. https://doi.org/10.3390/v14030576

AMA Style

Wang Y, Wang Y, Wang G, Li Q, Zhang Z, Li L, Lv Y, Yang Z, Guo J, Hong N. Molecular Characteristics and Incidence of Apple Rubbery Wood Virus 2 and Citrus Virus A Infecting Pear Trees in China. Viruses. 2022; 14(3):576. https://doi.org/10.3390/v14030576

Chicago/Turabian Style

Wang, Yanxiang, Ying Wang, Guoping Wang, Qingyu Li, Zhe Zhang, Liu Li, Yuzhuo Lv, Zuokun Yang, Jiashu Guo, and Ni Hong. 2022. "Molecular Characteristics and Incidence of Apple Rubbery Wood Virus 2 and Citrus Virus A Infecting Pear Trees in China" Viruses 14, no. 3: 576. https://doi.org/10.3390/v14030576

APA Style

Wang, Y., Wang, Y., Wang, G., Li, Q., Zhang, Z., Li, L., Lv, Y., Yang, Z., Guo, J., & Hong, N. (2022). Molecular Characteristics and Incidence of Apple Rubbery Wood Virus 2 and Citrus Virus A Infecting Pear Trees in China. Viruses, 14(3), 576. https://doi.org/10.3390/v14030576

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