Impact Factor 4.076

The 3rd most cited journal in Microbiology

Original Research ARTICLE

Front. Microbiol., 09 August 2017 | https://doi.org/10.3389/fmicb.2017.01514

Two Negative-Strand RNA Viruses Identified in Watermelon Represent a Novel Clade in the Order Bunyavirales

  • 1State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing, China
  • 2National Citrus Engineering Research Center, Citrus Research Institute, Southwest University, Chongqing, China
  • 3Institute of Plant Protection, Henan Academy of Agricultural Sciences, Zhengzhou, China

Two novel negative-sense, single-stranded (ss) RNA viruses were identified in watermelon plants and named watermelon crinkle leaf-associated virus 1 and 2 (WCLaV-1 and -2), respectively. The multipartite genomes consist of three RNA molecules of ~6.8, 1.4, and 1.3 kb. The genomes and the deduced proteins of RNA1 and RNA3 show features resembling those of members in the genus Phlebovirus and Tenuivirus; however, the predicted proteins encoded by RNA2 are related to the movement protein (MP) in the genus Ophiovirus and Emaravirus. Furthermore, these two viruses define a novel clade in the family Phenuiviridae, order Bunyavirales, which is phylogenetically related to the viruses in the above four genera. Moreover, after mechanical inoculation with WCLaV-1 seedlings of the natural host watermelon plants develop crinkling similar to those observed in the field. These findings enhance our understanding of the evolution and the classification of ssRNA viruses.

Introduction

Watermelon (Citrullus lanatus) is an economically important fruit crop in many countries and is nutritionally beneficial to humans (Weng and Sun, 2012). In open fields, watermelon plants are frequently exposed to viruses such as watermelon mosaic virus (WMV), zucchini yellow mosaic virus (ZYMV), cucumber mosaic virus (CMV), cucumber green mottle mosaic virus (CGMMV), tobacco mosaic virus (TMV), and squash mosaic virus (SqMV) that can cause great losses in production (Liu et al., 2009). In recent years, novel viruses that infect watermelon have been discovered through the use of next-generation sequencing (NGS). For example, Citrullus lanatus cryptic virus (CiLCV), a putative member of the genus Deltapartitivirus in the family Partitiviridae and two orthotospoviruses, watermelon bud necrosis virus (WBNV), and melon severe mosaic virus (MeSMV) in the family Tospoviridae have been detected in watermelon plants (Ciuffo et al., 2009; Li et al., 2011; Sela et al., 2013; Xin et al., 2017a).

In the tenth Report of the International Committee on Taxonomy of Viruses (ICTV), the Bunyaviridae Study Group of the International Committee on Taxonomy of Viruses (ICTV) established the new order Bunyavirales to accommodate nine families that comprise 13 genera (https://talk.ictvonline.org/taxonomy/). These viruses are a unique group that infect a wide range of hosts, including vertebrates, invertebrates, and plants (Webster et al., 2011). Three of the families (genera) contain members that have plants as their primary hosts: Fimoviridae (Emaravirus), Phenuiviridae (Tenuivirus), and Tospoviridae (Orthotospovirus). Members of the order Bunyavirales have similarities in genome organization and particle morphology, e.g., spherical in emaraviruses and orthotospoviruses and filamentous in tenuiviruses (Falk and Tsai, 1998; Mielke-Ehret and Mühlbach, 2012; Zhang et al., 2016). Their genomes consist of three to eight segments of linear negative-sense (or ambisense, depending on the genus) RNA, e.g., large (L), medium (M), and small (S), which encode the viral RNA polymerase, two glycoproteins as a single gene product that is usually co-translationally cleaved, and the nucleocapsid protein, respectively (Elliott and Blakqori, 2011; Léger and Lozach, 2015). The nucleotide sequences at the 3′-terminus and the 5′-terminus of each segment are complementary, resulting in RNAs that form panhandle structures (Kormelink et al., 2011; Jackson and Li, 2016).

In 2015 and 2016, during a field survey in Kāifēng, Hénán Province, China, virus-like symptoms on leaves of watermelon were observed. To identify potential causal agents, RNA sequencing (RNA-seq) and small RNA sequencing (sRNA-seq) were carried out. Here, we describe and discuss the molecular and biological properties of two new (-ss) RNA viruses that were discovered from the infected leaves, named watermelon crinkle leaf-associated viruses 1 and 2 (WCLaV-1 and -2), as possible representatives of a novel taxon in the family Phenuiviridae, order Bunyavirales.

Materials and Methods

Sample Collection, Small RNA Sequencing, and RNA Sequencing

In July 2015, leaf samples with mosaic and curling symptoms were collected from watermelon plants of cv. Hēixiùlóngjǔanfēng grown in the field in Hénán Province, China. Subsequently, in June 2016, watermelon samples with virus-like symptoms were randomly collected from Kāifēng County, Hénán Province. Samples were stored at −70°C.

Symptomatic leaves from cv. Hēixiùlóngjǔanfēng (sample KF-1) and healthy leaves (sample KF-0) were selected for small RNA-sequencing (sRNA-seq) in 2015; leaves coded KF-15 were used for RNA-sequencing (RNA-seq) in 2016. Total RNAs were extracted using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer's protocol but modified to repeat the chloroform extraction and subsequent centrifugation step (Liu et al., 2015). A BioAnalyzer 2100 (Agilent Technology, Santa Clara, CA, USA) was used to assess the quality of the purified RNA, then the purified RNA was quantified using a Nanodrop ND-2000 spectrophotometer (Infinigen Biotechnology, City of Industry, CA, USA). The Ribo-Zero Magnetic Kit (Plant Leaf) (Epicentre, Madison, WI, USA) was used to deplete the ribosomal RNA from transcriptomes for RNA-seq, then the library was constructed with a TruSeq RNA Sample Prep Kit (Illumina, San Diego, CA, USA), and a TruSeq Small RNA Sample Prep Kit (Illumina) was used to construct the sRNA libraries. Subsequently, the Illumina Hiseq2500 platform (1 × 50 bp read lengths) was used for sRNA-seq (SinoGenoMax, Beijing, China), and the Hiseq 4000 platform was used for RNA-seq with PE150 bp (Biomarker Biology Technology, Beijing, China). Raw reads from the Illumina platform were trimmed of adaptor sequences and low-quality reads by the CLC Genomics Workbench 9.5 (Qiagen, Valencia, CA, USA). The Velvet program was used for de novo assembly with a k-mer of 17 for small RNA data (Zerbino and Birney, 2008; Wu et al., 2010). The CLC Genomics Workbench 9.5 was used for de novo assembly for RNA-seq data. The assembled contigs were subsequently screened against the NCBI databases using a BLASTn and BLASTx search with standard parameters.

Amplification of Full Genomes

To verify the presence of viruses in field samples, we used RT-PCR and specific primers designed on the sequences of assembled contigs (Table S4). The sequences of the extreme ends of the genomic RNAs were determined employing the 5′- and 3′-RACE System for Rapid Amplification of cDNA Ends kits (Thermo Fisher Scientific, Waltham, MA, USA) (Than et al., 2016). The primers and other oligonucleotides used for RT-PCR and RACE analyses are listed in Table S3. The PCR products were separated by electrophoresis in 1.0% agarose gels and visualized by ethidium bromide staining, then purified with the Wizard SV Gel and PCR Clean-Up System (Promega, Madison, WI, USA). The final products were cloned into the PEASY-T5 cloning vector (TransGen Biotech, Beijing, China), and then inserted into Trans-T1 Chemically Competent Cells (TransGen). At least eight clones from each subclone were sequenced (Beijing Genomics Institute, Beijing, China).

Sequence and Phylogenetic Analyses

Sequences were analyzed and assembled using Vector NTI 11.5, then submitted to the GenBank database in NCBI using the software Sequin. The reconstructed genome was subjected to standard sequence analyses: (I) prediction of ORFs using the ORF Finder program in NCBI (http://www.ncbi.nlm.nih.gov/projects/gorf/); (II) identification of conserved and functional domains of the predicted proteins in WCLaV-1 and WCLaV-2 using the Conserved Domain Database (CDD) in NCBI (https://www.ncbi.nlm.nih.gov/Structure/cdd/cdd.shtml) (Marchler-Bauer et al., 2015) and the SMART tool (http://smart.embl-heidelberg.de/) (Letunic et al., 2015); analyses of the core RdRp and NP domains in other (-ss) RNA viruses using the CDD; (III) multiple sequence alignment of the core motifs in RdRp and the conserved nucleotide sequence at the 5′- and 3′-end by CLC Genomics Workbench 9.5; (IV) identity analyses using the needle program in WebLab (http://weblab.cbi.pku.edu.cn/); (V) multiple sequence alignments first using the Clustal W method and then MEGA 6.0 (Tamura et al., 2013) for phylogenetic tree construction using the neighbor-joining method with 1,000 bootstrap replicates, the Poisson model and pairwise gap deletion options. Abbreviations and accession numbers of viruses used in this study are listed in Table S5.

Mechanical Transmission in a Controlled Environment

The crude sap, extracted from leaves (about 1 g) infected by WCLaV-1 and homogenized in 10 mL 0.01 M phosphate buffer (PB) (pH = 7.2), was used for mechanical transmission tests on four-leaf-stage seedlings of watermelon cv. Hēixiùlóngjǔanfēng. Plants non-inoculated and inoculated only by PB were used as negative controls. After inoculation, the seedlings were maintained in a clean growth chamber at 25°C and 80% relative humidity with 16-h light/8-h dark. The mechanically inoculated plants were carefully observed for symptoms and subsequently tested for virus infection by RT-PCR assays with the virus-specific primer pairs listed in Table S4 and sequenced at 17, 42, and 67 days post inoculation (dpi). The mechanical transmission experiments were repeated three times with 20, 23, and 46 healthy watermelon plants inoculated, respectively. The cDNA was synthesized using Moloney murine leukemia virus (MoMLV) reverse transcriptase (Promega), and the PCR protocol employed 32 cycles after an initial denaturation at 95°C for 5 min. Each cycle consisted of denaturation at 95°C for 45 s, primer annealing for 45 s at 55°C and an extension for 1 min at 72°C, followed by an elongation step of 10 min at 72°C. The products were separated using 1.0% agarose gel electrophoresis.

Results

Identification of Two Unique Negative-Stranded RNA Viruses of Watermelon Using Next-Generation Sequencing

Our field investigation revealed diseased plants of watermelon with various virus-like symptoms, including leaf crinkling, mosaic, and bunchy top (Figures 1A,B); no symptoms were observed on the seemingly healthy plants (Figure 1C). To identify the viruses present in symptomatic samples, we subjected samples KF-1 and KF-15 to sRNA-seq and RNA-seq, respectively. Healthy plant KF-0 was also small RNA-sequenced. For sRNA-seq, a total of 22,235,704, and 20,667,689 reads with length between 18 and 27 nucleotides (nt) were obtained in sample KF-0 and KF-1, respectively. After de novo assembly and BlastX and BlastN searches against nr and nt databases, ZYMV (genus Potyvirus), WMV (genus Potyvirus), melon aphid-borne yellows virus (MABYV, genus Polerovirus), CiLCV (genus Deltapartitivirus), and four contigs homologous to viruses in the genera Phlebovirus and Ophiovirus were identified in sample KF-1 (Table S1). In contrast, contigs from KF-0 did not show any significant similarity with previously reported viruses. For RNA-seq, a total of 16,547,389 reads with a length of 150 nt were obtained after filtering the reads mapping to the watermelon genome (Guo et al., 2013). Interestingly, we identified CiLCV, CGMMV (genus Tobamovirus), watermelon virus A (WVA, family Betaflexiviridae), and three contigs homologous (high similarity with that identified in KF-1) to phleboviruses and ophioviruses in the sample KF-15 collected in 2016 (Table S1). Reverse transcription-polymerase chain reaction (RT-PCR) using specific primers (Table S4) and sequencing of the amplified products confirmed the presence of these contigs.

FIGURE 1
www.frontiersin.org

Figure 1. Symptoms on WCLaV-1- and WCLaV-2-infected watermelon leaves in the field (A) sample KF-1, crinkling, and mosaic, on the leaves of watermelon plants used for small RNA-sequencing for WCLaV-1 in 2015; (B) sample KF-15, cluster/crinkling on the top of leaves used for RNA-sequencing for WCLaV-2 in 2016; (C) symptomless leaves from uninfected watermelon plants. Red arrows in (A,B) indicate the crinkle leaf symptom.

Genomic Organization of WCLaV-1 and WCLaV-2

The genomes were obtained by RT-PCR and RACE-PCR (Figure S1) using virus-specific primers designed on the sequences of these contigs (Tables S3, S4). The assembled genomes are composed of three RNA segments, with complementary-sense RNAs (vcRNAs) of these two viruses predicted to contain a single ORF (Figures 2A,B). The nucleotide sequences reported here have been deposited in GenBank (respectively, accession KY781184, KY781185, KY781186 for RNA1, RNA2, and RNA3 of WCLaV-1; KY781187, KY781188, KY781189 for RNA1, RNA2, and RNA3 of WCLaV-2).

FIGURE 2
www.frontiersin.org

Figure 2. Schematic representation of the genome organizations of WCLaV-1 (A) and WCLaV-2 (B). Open reading frames (ORFs) and deduced products of each RNA are shown as gray boxes, with the amino acid (aa) length, estimated molecular weight (kDa), and function of the putative proteins. (C) Amino acid alignment between conserved RdRp premotif A and motifs A–E of WCLaV-1, WCLaV-2, and selected (-ss) RNA viruses. The conserved nucleotide stretches are present at the termini of WCLaV-1 (D) or WCLaV-2 (E) with tenuiviruses and phleboviruses. RC, reverse and complementary sequence; nt, nucleotides; v, viral RNA; EMARaV, European mountain ash ringspot-associated virus (AY563040); RRV, rose rosette virus (HQ871942); RLBV, raspberry leaf blotch virus (FR823299); TSWV, tomato spotted wilt virus (D10066); WSMoV, watermelon silver mottle virus (AF133128); RSV, rice stripe virus (D31879); RVFV, Rift Valley fever virus (DQ375403.1); SFTSV, severe fever with thrombocytopenia syndrome (SFTS) virus (HM745930); RGSV, rice grassy stunt virus (AB009656).

RNA1 for both viruses is the largest segment and has only one ORF (ORF1): 6,636 nt and 6,680 nt for WCLaV-1 and -2, which start, respectively, at position 6,584–6,586 and 6,634–6,636 and end at position 38–40 and 40–42. ORF1 encodes a putative RNA-dependent RNA polymerase (RdRp, p1) of 2,182 amino acids (aa) and 2,198 aa with a predicted molecular weight of 252.76 and 253.63 kDa, respectively, for WCLaV-1 and -2 (Figures 2A,B). A Blastp search with the algorithm PSI-BLAST using the deduced aa sequences of p1 show that they have about 25% identity to the L proteins of viruses in the genus Phlebovirus (Table S2). Additionally, SMART analyses recognized the conserved domain Bunya_RdRp (Pfame-value: 5.7e-51, 1.1e-47) in the p1 protein of these two viruses (Table S2). Interestingly, the conserved domain RRM_SF super family (RNA recognition motif superfamily, also known as RNA binding domain) was predicted by the Conserved Domain Database for p1 in WCLaV-2 (interval: 1,386–1,443; E-value: 4.22e-04), but was not found in WCLaV-1 (Marchler-Bauer et al., 2015). The alignment of RdRp of WCLaV-1, WCLaV-2, and other (-ss) RNA viruses indicates that the amino acid sequences of their RdRps contain the six conserved motifs (premotif A and motifs A–E), which represent highly conserved regions of the RdRp of the members in the order Bunyavirales (Figure 2C). Both WCLaV-1 and WCLaV-2 have the motif A (DATKWC), motif B (QGILHYTSS), and motif D (KS). Motif C includes the SDD sequence. The motif E tetrapeptide E(F/Y)xS, which is specifically conserved in the polymerases of segmented negative-sense RNA viruses (van Poelwijk et al., 1997; Kormelink et al., 2011), was also found in the proteins encoded by RNA1 of these two viruses (Figure 2C). In addition, we identified three basic residues (K, R, and R/K) in premotif A and a glutamic acid (E) downstream of premotif A (Figure 2C), which are also conserved in bunyavirus RdRps (Bruenn, 2003; Elbeaino et al., 2009). Moreover, the endonuclease domains H66-D78-PD95−96-E107-K126, similar to that in phlebo- and tenuiviruses and differing from the H-D-PD-DxK-T in orthotospoviruses, were also observed in the N-terminus of WCLaV-1 and WCLaV-2. These domains are probably involved in cap-snatching, a genome expression strategy to cap viral mRNAs proposed for members of the previous family Bunyaviridae and genus Tenuivirus (Biswas and Nayak, 1994; Reguera et al., 2010). In addition, the core RdRp domain of WCLaV-1 and WCLaV-2 is closely related to that of the tenuiviruses and phleboviruses, respectively; WCLaV-1 shares 27.6% aa sequence identity with Salehabad virus (SALV, genus Phlebovirus) and WCLaV-2 shares 27.2% aa sequence identity with rice stripe virus (RSV, genus Tenuivirus), both remarkably higher than with viruses in other genera, but far lower than the 75% identity between these two viruses.

RNA2 of WCLaV-1 (1,409 nt) and WCLaV-2 (1,423 nt) may, respectively, encode a protein (p2) of 44.29 and 46.11 kDa, spanning position 1,356–1,358 and 1,366–1,368 (AUG start codon) to position 196–198 and 155–157 (UGA stop codon) of vcRNA (Figures 2A,B). Interestingly, stretches highly rich in U were found in the 5′-terminus of RNA2 of both viruses. Blastp searches showed that the deduced aa sequences of RNA2 have 24% identity (28/116; E = 1.4) with the 54 kDa protein of citrus psorosis virus (CPsV, genus Ophiovirus), demonstrating a putative movement protein (MP) function, for WCLaV-1 and 19% identity (25/133; E = 1.4) to the 58 kDa MP of blueberry mosaic associated virus (BIMaV) in the genus Ophiovirus (Robles Luna et al., 2013; Borniego et al., 2016) for WCLaV-2 (Table S1). Therefore, the deduced protein p2 of WCLaV-1 and WCLaV-2 may be the putative MP. Nevertheless, no significant conserved domains in p2 of these two viruses have been identified using CDD and SMART analyses (Table S2). Pairwise comparisons showed that the aa sequence identities between the p2 of WCLaV-1, WCLaV-2 and other ssRNA viruses ranged, respectively, from 2.2 to 18.3% and 9.7 to 18.5%; the highest aa identities, 18.3 and 18.5%, respectively, were found with proteins encoded by RNA2 of groundnut bud necrosis virus (GBNV, the genus Orthotospovirus) (Table 1). Compared with the aa identity of 46.1% between these two viruses, the identities with known viruses are very low. Altogether, these data revealed that p2 encoded by RNA2 of these two viruses has unique features worth further study.

TABLE 1
www.frontiersin.org

Table 1. Amino acid sequences identities between WCLaV-1, WCLaV-2, and other (-ss) RNA viruses.

RNA3 of WCLaV-1 is 1,310 nt long and 1,260 nt for WCLaV-2. The ORFs start at position 1,251 and 1,198 nt, respectively, with an AUG start codon and terminate with an UAA stop codon at position 208 and 158 nt, respectively, putatively encoding a nucleoprotein (NP, p3) of 347 a and 346 aa with a molecular mass of 38.75 and 38.68 kDa, respectively (Figures 2A,B). It is worth noting that U-rich regions are also observed at the 5′-terminus in RNA3 of these two viruses, similar to the case in RNA2. According to the CDD and SMART analyses, RNA3 encodes proteins that contain the conserved domains Tenui_N super family (accession cl05345), the typical domains in nucleocapsid proteins of viruses in the genera Tenuivirus and Phlebovirus, with a Pfame-value of 5.4D-16, 3.2e-17, respectively, distinct from those in other genera. Moreover, the highest amino acid identities for the nucleocapsid protein, 17.2 and 17.4%, respectively, were found with SALV and sandfly fever Naples virus (SFNV, genus Phlebovirus) for WCLaV-1 and WCLaV-2, far lower than the 47.9% aa identity between these two viruses for nucleocapsid protein (Table 1). They share the highest aa identities, 21.2 and 23.2%, respectively, with the core NP domains of rice grassy stunt virus (RGSV, genus Tenuivirus) for WCLaV-1 and Punta Toro virus (PTV, genus Phlebovirus) for WCLaV-2 (Table 1).

In addition, the highly conserved nucleotides 5′-ACACAAAG-3′ at the 5′-terminus were found in each RNA1 segment, and the conserved nucleotides 3′-UGUGUUUC-5′ at the 3′-terminus were detected in each RNA2 and RNA3 segment (Figures 2D,E). The 5′- and 3′-terminal regions exhibited obvious complementarity, which is a common characteristic of genomic segments in tenuiviruses and phleboviruses (Elliott and Blakqori, 2011).

Phylogenetic Relationships of the WCLaV-1 and WCLaV-2 with Other (-ss) RNA Viruses

To reveal the relationships of WCLaV-1 and WCLaV-2 with known (-ss) RNA viruses, we constructed phylogenetic trees based on the amino acid sequences of the putative RdRp (Figure 3), MP (Figure S2A) and NP (Figure S2B). According to the generated trees based on amino acid sequences of RdRp, WCLaV-1, and WCLaV-2 are most closely related to viruses in the genera Tenuivirus, Phlebovirus, Goukovirus, Phasivirus, the family Phenuiviridae (Figure 3). However, phylogenetic analyses based on the amino acid sequences of RNA2-encoded proteins indicate these two viruses consistently cluster with emaraviruses and ophioviruses (Figure S2A). WCLaV-1 and WCLaV-2 form a distant and separate clade with the selected (-ss) viruses in all the phylograms, revealing that they may belong to a novel genus.

FIGURE 3
www.frontiersin.org

Figure 3. Phylogenetic analysis of WCLaV-1, WCLaV-2 based on the amino acid sequences of putative RdRp, and the corresponding proteins of representative negative-sense RNA viruses. The tree was generated by the neighbor-joining method using the Poisson model, and bootstrap values were estimated using 1,000 replicates and MEGA 6.0. The scale bar represents number of substitutions per amino acid site. Viruses in the same genera or family are labeled with the same colors. Viruses found in this study are marked with red squares; unclassified viruses are marked with black triangles.

Mechanical Transmission of WCLaV-1 in a Controlled Environment

The mechanical transmission experiment with WCLaV-1 was done three times with 20, 23, and 46 healthy watermelon plants inoculated, respectively. Each time, only one plant, respectively, of the inoculated samples was positive for WCLaV-1 infection by RT-PCR, thus conclusively showing infection rates of 5.0, 4.3, and 2.2%, respectively. For example, at 17 dpi, WCLaV-1 was detected only in watermelon plant No. 4 using RT-PCR with primer pair NP-1F/R, and the positive band was further identified by sequencing (Figure 4A, Table S4). Moreover, four uninoculated leaves of watermelon plant No. 4 were confirmed as infected by WCLaV-1 by RT-PCR detection with primer pairs NP-1F/R and MP-1F/R (Figure 4B), thus demonstrating that the watermelon plant was systemically infected by WCLaV-1. WCLaV-1 was not detected in the control plants non-inoculated and inoculated with phosphate buffer (PB). By 17 dpi, new leaves of watermelon seedlings had developed obvious mosaic and crinkling similar to the symptoms observed on plants in the field (Figure 4C). Subsequently, at 42 dpi, watermelon plant No. 4 displayed obvious crinkling and dwarf symptoms, which were lacking in the non-inoculated plant (Figure 4D). At 67 dpi, leaf crinkling was also observed (Figure 4E).

FIGURE 4
www.frontiersin.org

Figure 4. Test for mechanical transmission of watermelon crinkle leaf-associated virus 1 (WCLaV-1) to watermelon. (A,B) RT-PCR and agarose gel electrophoresis to detect infected watermelon leaves. (A) At 17 days post-inoculation (dpi), only sample No. 4 was infected; (B) detection of systematic infection in four leaves of sample No. 4 using two primer pairs in RT-PCR. Symptom development over time (C–E). (C) Crinkling at 17 dpi; (D) mosaic and stunt at 42 dpi compared with the non-inoculated plant and close ups of symptoms; (E) curling and mosaic at 67 dpi. M, DNA Marker DL2000 (Takara); PC, positive control; NC, uninfected watermelon control; WC, water control.

Discussion

The development of NGS and bioinformatics has accelerated the discovery of novel viruses (Adams et al., 2009; Hadidi et al., 2016; Massart et al., 2017). Unlike traditional techniques such as ELISA, PCR, or hybridization for diagnosing viruses, NGS technology requires no prior knowledge of the host or viral nucleotide sequences because virus-specific siRNAs are generated during an antiviral defense response in the plants infected by any virus, and the virus-specific siRNAs overlap in sequence, and thus can be assembled into long, contiguous fragments (contigs) of the invading viral genome (Al Rwahnih et al., 2009; Kreuze et al., 2009; Cao et al., 2014) With this technology, many viruses with negative-sense single-stranded (ss) RNA genomes, such as BIMaV and pigeonpea sterility mosaic virus 2 (PPSMV-2) have been characterized (Elbeaino et al., 2015; Li et al., 2015; Patil and Kumar, 2015). Indeed, our group previously used this technology to characterize the novel plant virus watermelon virus A (Xin et al., 2017b) and CiLCV, infecting watermelon (Xin et al., 2017a). In this study, we discovered and characterized two novel (-ss) RNA viruses from watermelon using NGS.

These two novel viruses share many common features with viruses in the newly established order Bunyavirales. First of all, the conserved sequences SDD in motif C, the signature of the core domain in RdRp for the segmented negative-stranded RNA viruses in the families Orthomyxoviridae, Arenaviridae, and the previous Bunyaviridae (now Peribunyaviridae), as opposed to the GDNQ characteristic of non-segmented negative-stranded RNA viruses, are found in the putative RdRp of these two viruses (Kormelink et al., 2011). Likewise, the sequence E(F/Y)xS in motif E, which is specifically conserved in the polymerases of segmented (-) RNA viruses (Müller et al., 1994; Reguera et al., 2010), were also discovered, demonstrating that these two viruses are segmented (-ss) RNA virus. Furthermore, the eight terminal conserved bases are identical to the corresponding regions in the genomes of phleboviruses and tenuiviruses; The deduced RdRp and NP are homologous to the corresponding proteins of bunyaviruses, sharing higher sequence identities and possessing the same conserved domains with viruses in the genera Phlebovirus and Tenuivirus (Kormelink et al., 2011). Phylogenetic analyses of viral RdRp and NP consistently place these two viruses in the bunyavirus supergroup, most closely related to, but distinct from phleboviruses and tenuiviruses (Elliott and Blakqori, 2011; Mielke-Ehret and Mühlbach, 2012). These observations suggest an evolutionary link between these two viruses with tenuiviruses and phleboviruses, distinguished from ophioviruses and emaraviruses.

Although WCLaV-1 and WCLaV-2 share sequence similarities with the members of the genera Tenuivirus, Phlebovirus, Ophiovirus, and Emaravirus, they have some unique properties and are markedly distant from all of the viruses in the above four genera. First, genomes of most (-ss) RNA viruses (except ophioviruses) generally contain a glycoprotein (GP) precursor, which has not been detected in WCLaV-1 and WCLaV-2 so far. Second, no conserved domains were found in the proteins encoded by RNA2, predicted as the MPs of WCLaV-1 and WCLaV-2, which are widely found in plant-infecting (-ss) RNA viruses (Kormelink et al., 2011). Third, all of the aa identities with the representative (-ss) RNA viruses in RNA1–3 are far lower than those between WCLaV-1 and WCLaV-2, e.g., the highest aa identity (22.7%) with known viruses was lower than the aa identity of 59.0% between these two viruses in RNA1. Fourth, they did not form clades that included any of the selected (-ss) RNA viruses; in contrast, these two viruses clustered on the same branch in the phylogenetic trees. Fifth, the genomes of all emaraviruses consist of four or more negative-sense RNAs, with the fully conserved 13-nt sequences in the 5′- and 3′-termini. Further, the longest RNA segments possess two ORFs in the same strand, a particular property of viruses in the genus Ophiovirus. The RNA3 segments in tenuiviruses and phleboviruses use an ambisense strategy, encoding one protein in the 5′-region of the vRNA and another protein in the vcRNA (Kormelink et al., 2011; Mielke-Ehret and Mühlbach, 2012). Because not all these features have been found in WCLaV-1 and WCLaV-2, they cannot now be assigned to any described genus based on molecular features and phylogenetic analyses. Taken together, these data indicate that WCLaV-1 and WCLaV-2 constitute part of a new virus cluster that bridges the genera Tenuivirus, Phlebovirus, Ophiovirus, and Emaravirus.

There is a highly conserved stretch 5′-ACACAAAG-3′ in the 5′-terminal of RNA1, but it is not found in RNA2 or RNA3. In contrast, the conserved nucleotides 5′-CUAUGUGU-3′ were only observed in the 3′-terminal of RNA2 and RNA3, not in RNA1, which may also imply a longer sequence. Although attempts to identify additional WCLaV-1 and WCLaV-2 genomic fragment(s) by RT-PCR and deep analyses of virus-derived sRNAs and contigs failed, we cannot exclude the potential existence of one or more additional genomic RNAs. Additionally, the sequences may include stretches due to cap snatching because the RACE was performed on total RNA.

Questions regarding the biology of WCLaV-1 and WCLaV-2 still need to be answered. Iris yellow spot virus (genus Orthotospovirus, family Tospoviridae) can be transmitted to healthy onions by mechanical inoculation in a controlled glasshouse at a relatively low transmission rate (Bag et al., 2015). Likewise, WCLaV-1 can also be mechanically transmitted to the natural host watermelon as confirmed in this study. The low transmission rates after mechanical inoculation indicate that other factors may be involved in the transmission of these viruses. Indeed, many bunyavirales are vectored by arthropods such as mosquitoes, flies, ticks, planthoppers, and thrips, and the associated arthropods also serve as virus reservoirs in many cases (Kormelink et al., 2011; Léger and Lozach, 2015). During the field survey, a large population of thrips was found in the watermelon plants, so we will investigate thrips species in the field where the samples were collected and test them in a transmission assay. We hope to find the major vector insects of these two viruses, then explore their roles in disease epidemics.

In conclusion, we are the first to report the genomic organization, phylogenetic relationship, and biological characterization of two new (-ss) RNA viruses, which are confirmed as members belonging to a previously undescribed novel clade in the family Phenuiviridae, order Bunyavirales.

Author Contributions

XW: conceived and designed the experiments. MX, WL, YR, and MC: performed the experiments. MX, WL, XZ, and MC: analyzed the data. XW, MX, and MC: wrote the manuscript. All authors read and approved the final manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This research was supported by the Inter-Governmental S&T Cooperation Proposal (2016YFE0131000) and the National Natural Science Foundation of China (31501611).

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fmicb.2017.01514/full#supplementary-material

References

Adams, I. P., Glover, R. H., Monger, W. A., Mumford, R., Jackeviciene, E., Navalinskiene, M., et al. (2009). Next-generation sequencing and metagenomic analysis: a universal diagnostic tool in plant virology. Mol. Plant Pathol. 10, 537–545. doi: 10.1111/j.1364-3703.2009.00545.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Al Rwahnih, M., Daubert, S., Golino, D., and Rowhani, A. (2009). Deep sequencing analysis of RNAs from a grapevine showing Syrah decline symptoms reveals a multiple virus infection that includes a novel virus. Virology 387, 395–401. doi: 10.1016/j.virol.2009.02.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Bag, S., Schwartz, H. F., Cramer, C. S., Havey, M. J., and Pappu, H. R. (2015). Iris yellow spot virus (Tospovirus: Bunyaviridae): from obscurity to research priority. Mol. Plant Pathol. 16, 224–237. doi: 10.1111/mpp.12177

PubMed Abstract | CrossRef Full Text | Google Scholar

Biswas, S. K., and Nayak, D. P. (1994). Mutational analysis of the conserved motif of Influenza A virus polymerase basic protein 1. J. Virol. 68, 1819–1826.

Google Scholar

Borniego, M. B., Karlin, D., Pe-a, E. J., Robles Luna, G., and García, M. L. (2016). Bioinformatic and mutational analysis of ophiovirus movement proteins, belonging to the 30K superfamily. Virology 498, 172–180. doi: 10.1016/j.virol.2016.08.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Bruenn, J. A. (2003). A structural and primary sequence comparison of the viral RNA dependent RNA polymerases. Nucleic Acids Res. 31, 1821–1829. doi: 10.1093/nar/gkg277

PubMed Abstract | CrossRef Full Text | Google Scholar

Cao, M. J., Du, P., Wang, X. B., Yu, Y. Q., Qiu, Y. H., Li, W. X., et al. (2014). Virus infection triggers widespread silencing of host genes by a distinct class of endogenous siRNAs in Arabidopsis. Proc. Natl. Acad. Sci. U.S.A. 111, 14613–14618. doi: 10.1073/pnas.1407131111

PubMed Abstract | CrossRef Full Text | Google Scholar

Ciuffo, M., Kurowski, C., Vivoda, E., Copes, B., Masenga, V., Falk, B. W., et al. (2009). A new Tospovirus sp. in cucurbit crops in Mexico. Plant Dis. 93, 467–474. doi: 10.1094/PDIS-93-5-0467

CrossRef Full Text | Google Scholar

Elbeaino, T., Digiaro, M., Alabdullah, A., De Stradis, A., Minafra, A., Mielke, N., et al. (2009). A multipartite single-stranded negative-sense RNA virus is the putative agent of fig mosaic disease. J. Gen. Virol. 90, 1281–1288. doi: 10.1099/vir.0.008649-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Elbeaino, T., Digiaro, M., Uppala, M., and Sudini, H. (2015). Deep-sequencing of dsRNAs recovered from mosaic-diseased pigeonpea reveals the presence of a novel emaravirus: pigeonpea sterility mosaic virus 2. Arch. Virol. 160, 2019–2029. doi: 10.1007/s00705-015-2479-y

CrossRef Full Text | Google Scholar

Elliott, R. M., and Blakqori, G. (2011). “Molecular biology of orthobunyaviruses,” in Bunyaviridae: Molecular and Cellular Biology, eds A. Plyusnin and R. M. Elliot (Norfolk, UK: Caister Academic Press), 1–39.

Google Scholar

Falk, B. W., and Tsai, J. H. (1998). Biology and molecular biology of viruses in the genus Tenuivirus. Annu. Rev. Phytopathol. 36, 139–163. doi: 10.1146/annurev.phyto.36.1.139

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, S. G., Zhang, J. G., Sun, H. H., Salse, J., Lucas, W. J., Zhang, H. Y., et al. (2013). The draft genome of watermelon (Citrullus lanatus) and resequencing of 20 diverse accessions. Nat. Genet. 45, 51–58. doi: 10.1038/ng.2470

PubMed Abstract | CrossRef Full Text | Google Scholar

Hadidi, A., Flores, R., Candresse, T., and Barba, M. (2016). Next-generation sequencing and genome editing in plant virology. Front. Microbiol. 7:1325. doi: 10.3389/fmicb.2016.01325

PubMed Abstract | CrossRef Full Text | Google Scholar

Jackson, A. O., and Li, Z. H. (2016). Developments in plant negative-strand RNA virus reverse genetics. Annu. Rev. Phytopathol. 54, 469–498. doi: 10.1146/annurev-phyto-080615-095909

PubMed Abstract | CrossRef Full Text | Google Scholar

Kormelink, R., Garcia, M. L., Goodin, M., Sasaya, T., and Haenni, A. L. (2011). Negative-strand RNA viruses: the plant-infecting counterparts. Virus Res. 162, 184–202. doi: 10.1016/j.virusres.2011.09.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Kreuze, J. F., Perez, A., Untiveros, M., Quispe, D., Fuentes, S., Barker, I., et al. (2009). Complete viral genome sequence and discovery of novel viruses by deep sequencing of small RNAs: a generic method for diagnosis, discovery and sequencing of viruses. Virology 388, 1–7. doi: 10.1016/j.virol.2009.03.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Léger, P., and Lozach, P. Y. (2015). Bunyaviruses: from transmission by arthropods to virus entry into the mammalian host first-target cells. Future Virol. 10, 859–881. doi: 10.2217/fvl.15.52

CrossRef Full Text | Google Scholar

Letunic, I., Doerks, T., and Bork, P. (2015). SMART: recent updates, new developments and status in 2015. Nucleic Acids Res. 43, 257–260. doi: 10.1093/nar/gku949

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, C. X., Shi, M., Tian, J. H., Lin, X. D., Kang, Y. J., Chen, L. J., et al. (2015). Unprecedented genomic diversity of RNA viruses in arthropods reveals the ancestry of negative-sense RNA viruses. Elife 4:e05378. doi: 10.7554/eLife.05378

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J. T., Yeh, Y. C., Yeh, S. D., Raja, J. A., Rajagopalan, P. A., Liu, L. Y., et al. (2011). Complete genomic sequence of watermelon bud necrosis virus. Arch. Virol. 156, 359–362. doi: 10.1007/s00705-010-0881-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, W. W., Gray, S., Huo, Y., Li, L., Wei, T. Y., and Wang, X. F. (2015). Proteomic analysis of interaction between a plant virus and its vector insect reveals new functions of hemipteran cuticular protein. Mol. Cell. Proteomics 14, 2229–2242. doi: 10.1074/mcp.M114.046763

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Y., Wang, Y. N., Wang, X. F., and Zhou, G. H. (2009). Molecular characterization and distribution of Cucumber green mottle mosaic virus in China. J. Phytopathol. 157, 393–399. doi: 10.1111/j.1439-0434.2008.01509.x

CrossRef Full Text | Google Scholar

Marchler-Bauer, A., Derbyshire, M. K., Gonzales, N. R., Lu, S., Chitsaz, F., Geer, L. Y., et al. (2015). CDD: NCBI's conserved domain database. Nucleic Acids Res. 43, 222–226. doi: 10.1093/nar/gku1221

PubMed Abstract | CrossRef Full Text | Google Scholar

Massart, S., Candresse, T., Gil, J., Lacomme, C., Predajna, L., Ravnikar, M., et al. (2017). A framework for the evaluation of biosecurity, commercial, regulatory, and scientific impacts of plant viruses and viroids identified by NGS technologies. Front. Microbiol. 8:45. doi: 10.3389/fmicb.2017.00045

PubMed Abstract | CrossRef Full Text | Google Scholar

Mielke-Ehret, N., and Mühlbach, H. P. (2012). Emaravirus: a novel genus of multipartite, negative strand RNA plant viruses. Viruses 4, 1515–1536. doi: 10.3390/v4091515

PubMed Abstract | CrossRef Full Text | Google Scholar

Müller, R., Poch, O., Delarue, M., Bishop, D. H. L., and Bouloy, M. (1994). Rift valley fever virus L segment: correction of the sequence and possible functional role of newly identified regions conserved in RNA-dependent polymerases. J. Gen. Virol. 75, 1345–1352. doi: 10.1099/0022-1317-75-6-1345

PubMed Abstract | CrossRef Full Text | Google Scholar

Patil, B. L., and Kumar, P. L. (2015). Pigeonpea sterility mosaic virus: a legume-infecting Emaravirus from South Asia. Mol. Plant Pathol. 16, 775–786. doi: 10.1111/mpp.12238

PubMed Abstract | CrossRef Full Text | Google Scholar

Reguera, J., Weber, F., and Cusack, S. (2010). Bunyaviridae RNA polymerases (L-protein) have an N-terminal, influenza-like endonuclease domain, essential for viral cap-dependent transcription. PLoS Pathog. 6, 8871–8887. doi: 10.1371/journal.ppat.1001101

PubMed Abstract | CrossRef Full Text | Google Scholar

Robles Luna, G., Pe-a, E. J., Borniego, M. B., Heinlein, M., and Garcia, M. L. (2013). Ophioviruses CPsV and MiLBVV movement protein is encoded in RNA 2 and interacts with the nucleocapsid protein. Virology 441, 152–161. doi: 10.1016/j.virol.2013.03.019

CrossRef Full Text | Google Scholar

Sela, N., Lachman, O., Reingold, V., and Dombrovsky, A. (2013). A new cryptic virus belonging to the family Partitiviridae was found in watermelon co-infected with Melon necrotic spot virus. Virus Genes 47, 382–384. doi: 10.1007/s11262-013-0937-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Tamura, K., Stecher, G., Peterson, D., Filipski, A., and Kumar, S. (2013). MEGA6: molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 30, 2725–2729. doi: 10.1093/molbev/mst197

PubMed Abstract | CrossRef Full Text | Google Scholar

Than, W., Qin, F. L., Liu, W. W., and Wang, X. F. (2016). Analysis of Sogatella furcifera proteome that interact with P10 protein of Southern rice black-streaked dwarf virus. Sci. Rep. 6:32445. doi: 10.1038/srep32445

PubMed Abstract | CrossRef Full Text | Google Scholar

van Poelwijk, F., Prins, M., and Goldbach, R. (1997). Completion of the impatiens necrotic spot virus genome sequence and genetic comparison of the L proteins within the family Bunyaviridae. J. Gen. Virol. 78, 543–546. doi: 10.1099/0022-1317-78-3-543

PubMed Abstract | CrossRef Full Text | Google Scholar

Webster, C. G., Reitz, S. R., Perry, K. L., and Adkins, S. (2011). A natural M RNA reassortant arising from two species of plant-and insect-infecting bunyaviruses and comparison of its sequence and biological properties to parental species. Virology 413, 216–225. doi: 10.1016/j.virol.2011.02.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Weng, Y. L., and Sun, S. J. (2012). Therapeutic hypothermia after cardiac arrest in adults: mechanism of neuroprotection, phases of hypothermia, and methods of cooling. Crit. Care Clin. 28, 231–243. doi: 10.1016/j.ccc.2011.10.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, Q. F., Luo, Y. J., Lu, R., Lau, N., Lai, E. C., Li, W. X., et al. (2010). Virus discovery by deep sequencing and assembly of virus-derived small silencing RNAs. Proc. Natl. Acad. Sci. U.S.A. 107, 1606–1611. doi: 10.1073/pnas.0911353107

PubMed Abstract | CrossRef Full Text | Google Scholar

Xin, M., Cao, M. J., Liu, W. W., Ren, Y. D., Lu, C. T., and Wang, X. F. (2017a). The genomic and biological characterization of Citrullus lanatus cryptic virus infecting watermelon in China. Virus Res. 232, 106–112. doi: 10.1016/j.virusres.2017.02.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Xin, M., Zhang, P. P., Liu, W. W., Ren, Y. D., Cao, M. J., and Wang, X. F. (2017b). The complete nucleotide sequence and genome organization of a novel betaflexivirus infecting Citrullus lanatus. Arch. Virol. doi: 10.1007/s00705-017-3461-7. [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

Zerbino, D. R., and Birney, E. (2008). Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 18, 821–829. doi: 10.1101/gr.074492.107

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z. K., Zheng, K. Y., Dong, J. H., Fang, Q., Hong, J., and Wang, X. F. (2016). Clustering and cellular distribution characteristics of virus particles of Tomato spotted wilt virus and Tomato zonate spot virus in different plant hosts. Virol. J. 13:11. doi: 10.1186/s12985-016-0466-x

CrossRef Full Text | Google Scholar

Keywords: watermelon (Citrullus lanatus), Next Generation Sequencing (NGS), negative-sense single-stranded RNA virus, novel virus clade, Bunyavirales

Citation: Xin M, Cao M, Liu W, Ren Y, Zhou X and Wang X (2017) Two Negative-Strand RNA Viruses Identified in Watermelon Represent a Novel Clade in the Order Bunyavirales. Front. Microbiol. 8:1514. doi: 10.3389/fmicb.2017.01514

Received: 28 April 2017; Accepted: 27 July 2017;
Published: 09 August 2017.

Edited by:

Sead Sabanadzovic, Mississippi State University, United States

Reviewed by:

Jens H. Kuhn, Rocky Mountain Laboratories, NIAID (NIH) Integrated Research Facility (IRF-Frederick), United States
Massimo Turina, Consiglio Nazionale Delle Ricerche (CNR), Italy

Copyright © 2017 Xin, Cao, Liu, Ren, Zhou and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Mengji Cao, caomengji@cric.cn
Xifeng Wang, xfwang@ippcaas.cn