A putative new ampelovirus associated with grapevine leafroll disease.

A putative new ampelovirus was detected in Vitis vinifera cv. Carnelian showing mild leafroll symptoms and molecularly characterized. The complete genome consisted of 13,625 nt and had a structure similar to that of members of subgroup I in the genus Ampelovirus (fam. Closteroviridae). In-depth analyses showed that the virus from cv. Carnelian is the most distinct member of the "GLRaV-4 lineage" of ampeloviruses, which comprises GLRaV-4, -5, -6, -9, and the recently characterized GLRaV-Pr, and GLRaV-De. This virus appears to be a new member of the family Closteroviridae, for which the provisional name grapevine leafroll-associated Carnelian virus is proposed.

Leafroll disease (LRD) is one of the most economically important and widespread viral diseases of cultivated grapevines. Several filamentous, phloem-restricted viruses referred to as grapevine leafroll-associated viruses (GLRaVs) have been identified and characterized from leafroll-affected grapevines [21]. Some of these viruses have been recognized as belonging to definitive species in the genera Closterovirus (GLRaV-2) and Ampelovirus (GLRaV -1, -3 and -5) or unassigned species in the family Closteroviridae (GLRaV-7) [22], whereas several viruses/variants (i.e. GLRaV-4, -6, -9, -Pr, -De) still await official taxonomic assignment. The ampeloviruses associated with LRD are transmitted by grafting and by several species of mealybugs and soft scales [reviewed in ref. 23].

Mild downward rolling and premature reddening of the leaves, resembling leafroll disease, were observed in 2005 and 2006 on a grapevine (Vitis vinifera) cv. Carnelian [(Carignane × Cabernet Sauvignon) × Grenache] in a grape virus collection at the University of California, Davis. The presence of leafroll infection was confirmed by grafting onto the leafroll-specific indicator V. vinifera cv. Cabernet Franc (Fig. 1a). The original plant repeatedly tested negative in ELISA and RT-PCR for known grapevine leafroll-associated viruses, thus justifying further investigation. The emphasis in this paper is given to a detailed comparison of this virus to other “GLRaV-4-like” viruses/variants.

Fig. 1

a Mild leafroll symptoms induced by GLRaCV on indicator Vitis vinifera cv. Cabernet Franc. b Electrophoretic patterns of dsRNAs extracted from cv. Carnelian (lane 1) compared to dsRNAs extracted from grapevine infected with grapevine rupestris stem pitting-associated virus (lane 2) and from GLRaV-2-infected Nicotiana benthamiana (lane 3)

Double-stranded RNA (dsRNA) was isolated from diseased tissues using double phenol–chloroform extractions and CF-11 column chromatography [6], purified by selective enzymatic digestions [27] and analyzed by polyacrylamide gel electrophoresis (PAGE). The multiple high-molecular-weight dsRNA molecules in cortical tissue of the diseased cv. Carnelian were interpreted to be replicative forms of a clostero-like virus. The largest dsRNA molecule, interpreted as a full genomic-size replicon, migrated slightly faster than the corresponding molecule of grapevine leafroll-associated virus 2, indicating a smaller genome size. Based on differences in migration rates compared to standards (GLRaV-2 and grapevine rupestris stem pitting-associated virus, GRSPaV), the full genome-size molecule was estimated to be ca. 13.5 kbp (Fig. 1b).

Initial molecular characterization of this virus was carried out on a portion of the HSP70-homologue (HSP70h) cistron between motifs P1 and P2 that was amplified using universal degenerate primers for closterovirids [30]. Initial sequence data showed that this virus was fairly different from other LRD-associated viruses (ca. 70% identity at the amino acid level) and other members of the genus Ampelovirus.

Further sequence data were generated by random-primed reverse transcription of purified dsRNA preparations and cloning into pGEM-T Easy plasmid (Promega Corporation, USA). Selected plasmids were sequenced at the UC Davis and/or Mississippi State University DNA Sequencing Facility. Sequence data were assembled and analyzed using DNAStar (Lasergene) software. The sequences that were generated were used to design virus-specific primers in order to progress towards sequencing the whole genome. Each genomic region was verified by sequencing multiple independent clones in order to achieve at least 5× coverage/nucleotide. Uniform sequence data among multiple sequenced clones was further proof of the absence of related GLRaVs in the accession of cv. Carnelian. Sequences of the viral termini were determined by artificial polyadenylation (3′end) and/or RACE (5′end) methodology using a GeneRacer kit (Invitrogen, USA). Comparisons with sequences available in the NCBI/GenBank database were performed using the BLAST [2] and CDD [20] resources. Phylogenetic analysis with bootstraping consisting of 1,000 pseudoreplicates was performed with ClustalX software [17] using the neighbor-joining algorithm, and trees were visualized with the Dendroscope program [15].

The genome of this virus, provisionally denoted grapevine leafroll-associated Carnelian virus (GLRaCV), consists of 13,625 nt (Fig. 2) and closely resembles that of several recently sequenced “short” ampeloviruses: plum bark necrosis stem pitting-associated virus (PBNSPaV) [3], GLRaV-Pr and -De [18, 19], and pineapple mealybug wilt-associated viruses 1 and 3 (PMWaV-1 and PMWaV-3) [24, 29].

Fig. 2

Diagrammatic representation of the GLRaCV genome to scale. Boxes depict putative ORFs. PRO protease; MTR methyltransferase; HEL helicase; RdRp RNA-dependent RNA polymerase; p5 hydrophobic 5 K protein; HSP70h heat shock 70 protein homologue; p60 60 K protein; CP coat protein; p23 23 K protein with unknown function

The GLRaCV genome starts with an AU-rich (56%) 214-nt-long non-coding region (5′NCR). The first ORF (ORF1a) extends for 2,288 codons and codes for a replication-associated polyprotein with an estimated molecular mass of 254.7 kDa (p255) containing hallmark domains of (in the 5′-to-3′ direction) papain-like protease [12], viral methyltransferase (MTR, pfam 01660) [26], and viral helicase superfamily 1 (HEL, pfam 01443) [11]. Additionally, p255 contains a conserved AlkB domain (2OG-Fe(II) oxygenase superfamily; pfam 03171) [5], which is located between MTR and HEL and was reported previously to be present in some other ampeloviruses, i.e. little cherry virus 2, GLRaV-3 [5] and GLRaV-Pr [19]. Functional domains identified in p255 showed considerable conservation of amino acids with related products encoded by GLRaV-Pr and -9. Helicase was the most conserved domain, sharing ca. 70% aa identity with GLRaV-Pr and -9, whereas about 60% of the residues of MTR and AlkB are identical to those of these closely related viruses (Table 1). As in other closterovirids, conserved functional domains were separated by regions with variable aa content and length, reducing the overall identity between GLRaCV and GLRaV-Pr (the closest fully sequenced virus) to ca. 45%. ORF1a terminates with the sequence guuUAAca (stop codon in capital letters), which is identical to the corresponding region in PMWaV-1 [24] and is presumably involved in a +1 ribosomal frameshift phenomenon. ORF1b overlaps ORF1a for 8 nucleotides and codes for a polypeptide with a deduced molecular mass of 59.1 kDa containing the conserved sequence motifs of a viral RNA-dependent RNA polymerase (RdRp) of the “alphavirus-like” viruses [16]. This protein shares 76–77% identity with orthologs of the most closely related ampeloviruses (GLRaV-9 and –Pr; Table 1). Identities with other short ampeloviruses were much lower (i.e. 37% with PBNSPaV and 53% with PMWaV-1).

Table 1

Percent nucleotide (nt) or amino acid (aa) identity between genomes/genome products of GLRaCV and closely related grapevine leafroll-associated viruses/isolates of the subgroup I in the genus Ampelovirus

Virus	% nt identity			% aa identity
	Genome	5′UTR	3′UTR	MTR	HEL	ORF1b	p5	HSP70h	p60	CP	p23
GLRaV-4								66*	65	78	52*
GLRaV-5*	63*							66	63	77	56
GLRaV-6								66*
GLRaV-9*	60*	53*	87*	55	69	77	73	67	64	78	56
GLRaV-Pr	60	58	76	58	71	76	69	69	64	77	50
GLRaV-De	63*							66	63	76	62*

* Partial sequences

After a 4-nt-long intergenic region, the coding part of the genome continues with a small ORF (ORF2) encoding a 46-aa-long hydrophobic protein (p5) with membrane-binding properties. It shared the highest level of amino acid sequence identity (73%) with the hydrophobic protein encoded by the GLRaV-9 genome.

The 58-kDa polypeptide encoded by ORF3 contains the conserved motifs of cellular chaperones and shares <70% identical amino acids with HSP70 homologues encoded by the GLRaV-4, -5, -6, -9, -Pr and -De genomes (Table 1). Identities with other related viruses were in the range of 43 (PBNSPaV) to 59% (PMWaV-1). ORF4 partially overlaps with the previous cistron and extends for 1,620 nt. It encodes a 539-aa-long protein related to the ca. 60-kDa proteins of other closterovirids with a role in virion/tail assembly/virus movement [8]. This ORF has fewer than 65% of its amino acids in common with the corresponding products of closely related “GLRaV-4-like” viruses/variants, and only 48–49% with PMWaV-1 and -3.

ORF5 is separated by 69 nt from the previous cistron and encodes a viral coat protein with a predicted molecular mass of 29.2 kDa. Serine, arginine and aspartic acid residues reported to be conserved in coat proteins (CP) of filamentous and rod-shaped phytoviruses [7] were identified at positions 125, 172 and 213, respectively. This protein shares 76–78% aa identity with orthologs of GLRaV-5 and related viruses, including the recently reported GLRaV-Pr and -De [18, 19]. Values of amino acid content conservation with CPs of PBNSPaV, PMWaV-1 and -3 are 27, 55 and 62%, respectively.

The 3′-end-proximal ORF6 codes for a putative p23 protein with 50–56% identity with the corresponding genome products of related viruses. Direct comparison of p23 with the viral CP showed low overall identity (16%) scattered randomly throughout the protein. The absence of significant amino acid conservation between the two proteins, as well as the lack of the “closter_coat motif” (pfam 01785) suggests that p23 may not be a paralog of the viral coat protein and may have a completely different function. The genome terminates with a 131-nt-long 3′ NCR that is 87 and 76% identical to the corresponding regions of GLRaV-9 and GLRaV-Pr, respectively. This protein shares 30 and 37% aa identity with PMVaV-1 and -3, and <15% with PBNSPaV.

In phylogenetic analysis, performed on several phylogenetically relevant gene products (i.e. RdRp, HSP70h, CP), GLRaCV always grouped with members of subgroup I in the genus Ampelovirus. Nevertheless, independent of the gene used in analysis, GLRaCV appeared to be the most distinct member of the “GLRaV-4 clade”, as it was consistently placed slightly apart from the rest of these viruses (Fig. 3). This topology was also supported by Bayesian inference and maximum-likelihood methods applied in MrBayes [14] and PhyML [13] software, respectively (not shown).

Fig. 3

Phylogenetic tree constructed using HSP70h amino acid sequences of members of some definitive and tentative species in the family Closteroviridae. Bootstrap values are indicated on the nodes. Sequences used to construct the tree and their accession numbers are as follows: beet yellows virus (BYV, NC_001598), beet yellow stunt virus (BYSV, U51931), blackberry yellow vein associated virus (BYVaV, NC_006963), carnation necrotic fleck virus (CNFV, EU884443), carrot yellow leaf virus (CYLV, NC_013007), citrus tristeza virus (CTV, NC_001661), cucurbit yellow stunting disorder virus (CYSDV, NC_004810), grapevine leafroll-associated virus 1 (GLRaV-1, AF195822), 2 (GLRaV-2, NC_007448), 3 (GLRaV-3, NC_004667), 4 (GLRaV-4, AM162280), 5 (GLRaV-5, AF233934), 7 (GLRaV-7, Y15987), 9 (GLRaV-9, AY297819), De (GLRaV-De, AM494935), Pr (GLRaV-Pr, NC_011702), grapevine leafroll-associated Carnelian virus (GLRaCV, FJ907331), little cherry virus 1 (LChV-1, NC_001836), little cherry virus 2 (LChV-2, NC_005065), lettuce infectious yellows virus (LIYV, NC_003618), mint virus-1 (MV-1, NC_006944), pineapple mealybug wilt-associated virus 1 (PMWaV-1, NC_010178), pineapple mealybug wilt-associated virus 2 (PMWaV-2, AF283103), pineapple mealybug wilt-associated virus 3 (PMWaV-3, DQ399259), plum bark necrosis stem pitting associated virus (PBNSPaV, NC_009992), potato yellow vein virus (PYVV, NC_006063), raspberry leaf mottle virus (RLMoV, NC_008585), sweet potato chlorotic stunt virus (SPCSV, NC_004124) and tomato infectious chlorosis virus (TICV, NC_013259). Outgroup: putative heat shock protein 70 of Arabidopsis thaliana (AAN71949)

Based on our results, the virus identified in the LRD-affected grapevine plant of cv. Carnelian appears to be a member of a novel viral species in the genus Ampelovirus. It induced mild leafroll symptomatology in the original plant as well as in the leafroll-specific indicator host. Virus was not detected by RT–PCR using a panel of primers specific for GLRaV-1 to -7 and -9 that are routinely used for diagnostic purposes at the Foundation Plant Services laboratory at UC Davis.

Complete sequencing of this virus showed a genomic organization similar to those of GLRaV-Pr [19], GLRaV-9 [4], PMWaV-1 [24], PMWaV-3 [29] and PBNSPaV [3], the genomes of which are completely or almost completely known, and presumably to a number of other, partially characterized “short GLRaVs” (i.e. GLRaV-4, -5, -6, -De and -WC15) [1, 10, 18, 25]. However, in-depth analysis showed that the isolate from cv. Carnelian shares rather limited amino acid identity with sequenced members of this subgroup of ampeloviruses. With the exception of the viral CP, which shares 76–78% common residues with its closest relatives, the amino acid content of other putative products of the GLRaCV genome differs by more than 25% from orthologs in the closely related GLRaV-5, -9, -De and -Pr, as well as with GLRaV-4 and -6 (unpublished data). Independent of the genome product used for analysis, identities with other approved or putative “short” ampeloviruses (i.e. PMWaV-1, -3, PBNSPaV) were at least 10–15% lower than values with “GLRaV-4-like viruses”. The replication-associated polyprotein (p255) and p23 (encoded by ORF6) appear to be less conserved than their orthologs in GLRaV-9 and -Pr (Table 1). In the case of p255, the overall polyprotein shares only 44 and 48% aa identity with GLRaV-9 and GLRaV-Pr, respectively, due to the poorly conserved regions in the hallmark domains (MTR-AlkB-HEL). Indeed, when compared separately, the functional domains showed levels of amino acid identity comparable to the rest of the genome (see Table 1). The putative p23 protein shares 50–56% common residues with the corresponding products of GLRaV-4, -5, -9, -Pr and -De. These values are in the range of what has been reported for the corresponding products of PMWaV-1 and -3 (53%) [29], two closely related viruses with a similar organization that share 63–73% conservation in the rest of their genomes. All of these values are far below the current species demarcation thresholds for the family Closteroviridae [22], strongly suggesting that the virus isolate from cv. Carnelian could be a member of a novel species in the genus Ampelovirus. Considering that the current situation regarding the nomenclature and taxonomy of GLRaVs is rather complex due to the uncertain fate of GLRaV-4, -6 and a suite of recently reported similar viruses/strains [4, 9, 18, 19, 25, 28], the provisional name grapevine leafroll-associated Carnelian virus (GLRaCV) is proposed for this virus. Final assignment of GLRaCV will be discussed during the critical revision of the taxonomy of grapevine-infecting viruses of the genus Ampelovirus, which is currently underway [21, and G.P. Martelli personal communication].