# Comparison of srological and molecular methods with high-throughput sequencing for the detection and quantification of grapevine fanleaf virus in vineyard samples

Emmanuelle Vigne, Shahinez Garcia, Véronique Komar, Olivier Lemaire, Jean-Michel Hily

<figure><img src="image_1.png" /></figure>

To cite this version:

Emmanuelle Vigne, Shahinez Garcia, Véronique Komar, Olivier Lemaire, Jean-Michel Hily. Comparison of srological and molecular methods with high-throughput sequencing for the detection and quantification of grapevine fanleaf virus in vineyard samples. Frontiers in Microbiology, 2018, 9, 10.3389/fmicb.2018.02726. [hal-02623552](https://hal.inrae.fr/hal-02623552)

HAL Id: hal-02623552

https://hal.inrae.fr/hal-02623552

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Check for
updates

# Comparison of Serological and Molecular Methods With High-Throughput Sequencing for the Detection and Quantification of Grapevine Fanleaf Virus in Vineyard Samples

Emmanuelle Vigne, Shahinez Garcia, Véronique Komar, Olivier Lemaire and Jean-Michel Hily*

L'UMR Santé de la Vigne et Qualité du Vin, INRA-Université de Strasbourg, Colmar, France

## OPEN ACCESS

**Edited by:**

Heiko Ziebell,
Julius Kühn-Institut, Germany

**Reviewed by:**

Pasquale Saldarelli,
Istituto di Virologia Vegetale (CNR),
Italy
Feng Qu,
The Ohio State University,
United States

**\*Correspondence:**

Jean-Michel Hily
jean-michel.hily@inra.fr

**Specialty section:**

This article was submitted to
Virology,
a section of the journal
Frontiers in Microbiology

**Received:** 25 June 2018

**Accepted:** 24 October 2018

**Published:** 22 November 2018

**Citation:**

Vigne E, Garcia S, Komar V, Lemaire O and Hily J-M (2018) *Comparison of Serological and Molecular Methods With High-Throughput Sequencing for the Detection and Quantification of Grapevine Fanleaf Virus in Vineyard Samples*. Front. Microbiol. 9:2726. doi: 10.3389/fmicb.2018.02726

Grapevine fanleaf virus (GFLV) is the main causal agent of fanleaf degeneration, the most damaging viral disease of grapevine. GFLV is included in most grapevine certification programs that rely on robust diagnostic tools such as biological indexing, serological methods, and molecular techniques, for the identification of clean stocks. The emergence of high throughput sequencing (HTS) offers new opportunities for detecting GFLV and other viruses in grapevine accessions of interest. Here, two HTS-based methods, i.e., RNAseq and smallRNAseq (focusing on the 21 to 27 nt) were explored for their potential to characterize the virome of grapevine samples from two 30-year-old GFLV-infected vineyards in the Champagne region of France. smallrnaseq was optimal for the detection of a wide range of viral species within a sample and RNAseq was the method of choice for full-length viral genome assembly. The implementation of a protocol to discriminate between low GFLV titer and *in silico* contamination (intra-lane contamination due to index misassignment) during data processing was critical for data analyses. Furthermore, we compared the performance of semi-quantitative DAS-ELISA (double antibody enzyme-linked immunosorbent assay), RT-qPCR (Reverse transcription-quantitative polymerase chain reaction), Immuno capture (IC)-RT-PCR, northern blot for viral small interfering RNA (vsiRNA) detection and RNAseq for the detection and quantification of GFLV. While detection limits were variable among methods, as expected, GFLV diagnosis was consistently achieved with all of these diagnostic methods. Together, this work highlights the robustness of DAS-ELISA, the current method routinely used in the French grapevine certification program, for the detection of GFLV and offers perspectives on the potential of HTS as an approach of high interest for certification.

**Keywords:** grapevine, GFLV, detection, serological and molecular methods, high-throughput sequencing, contamination evaluation protocol
INTRODUCTION

More than 80 viruses and five viroids have been identified in *Vitis spp.*, making grapevine the most virus-infected agricultural commodity known to date (Martelli, 2017). While some viruses might not be directly linked to a particular disease, others are of great economic concern to the grapevine industry. The grapevine fanleaf virus (GFLV) from the genus *Nepovirus* in the family *Secoviridae* (Sanfaçon, 2015) causes fanleaf degeneration disease (Schmitt-Keichinger et al., 2017). This disease is of major economic concern (Andret-Link et al., 2004). The genome of GFLV is composed of two single-stranded positive-sense RNA molecules (RNA1 and RNA2). Both genomic RNAs are necessary for infection and can be associated with a satellite RNA (RNA3) (Schmitt-Keichinger et al., 2017). GFLV exhibits strong genetic diversity with genome sequences being divergent up to 20%, with many recombination events being identified on both genomic RNAs from field isolates (Vigne et al., 2004, 2008; Zhou et al., 2017; Hily et al., 2018c).

Preventive measures based on the use of clean grapevines derived from certified as free of GFLV stocks is the most reliable option to manage fanleaf disease. The identification of clean vines through extensive and robust diagnosis is foundational to certification schemes (Golino et al., 2017a,b). GFLV and other viruses need to be undetectable in order for vines to be certified. Several diagnostic methodologies are available for GFLV such as bioassays, serological assays and polymerase chain reaction (PCR)-based assays (Golino et al., 2017a). A combination of both ELISA and PCR-based techniques, known as immunocapture-PCR (IC-PCR), is also used for GFLV detection. Recently, high-throughput sequencing (HTS) is explored as a new approach for virus discovery and detection. HTS has the potential for inclusion in certification programs (Al Rwahnih et al., 2015; Maree et al., 2018). Different nucleic acid preparations and different sequencing platforms are routinely used for HTS (Adams et al., 2009; Al Rwahnih et al., 2009; Kreuze et al., 2009). The benefits and pitfalls of each of them have

been extensively reviewed (Kircher and Kelso, 2010; Barzon et al., 2013; Reuter et al., 2015; Roossinck et al., 2015). Briefly, a virus and its genetic make-up can be detected by (i) characterizing total RNA or DNA, with or without specific enrichment steps (Dayaram et al., 2012; Beuve et al., 2018), (ii) double-stranded (ds) RNA for viruses with dsRNA genomes and RNA viruses and viroids that form dsRNA during their replication (Coetzee et al., 2010; Blouin et al., 2016), (iii) products of RNA interference, an adaptive antiviral plant defense mechanism, such as vsiRNA (viral small-interfering RNA) (Donaire et al., 2009; Kreuze et al., 2009), and (iv) encapsidated nucleic acids using the VANA (Virion-Associated Nucleic Acid) approach (Filloux et al., 2015). HTS revolutionizes grapevine virus diagnostics with unbiased and massive sequencing, compared to serological and molecular assays which rely on specific antibodies or primers/probes (Al Rwahnih et al., 2015). However, prior to the adoption of HTS for routine and generic viral detection, many challenges have to be overcome and some factors need to be optimized such as: (i) the development of nucleic acid preparation protocols for consistent diagnostic results, (ii) the development of straightforward and easy-to-use bioinformatics tools, and (iii) the identification of contaminations such as 'between-run' contaminations or 'intra-lane' cross contamination, also known as 'index hopping.' The former type of contaminations derives from a carry-over of contaminants from previous runs, while the latter arises when multiple libraries are pooled in a single lane and sequenced simultaneously (Kircher et al., 2012; Maree et al., 2018).

In this study, we compared HTS, including RNAseq and smallRNAseq, with serological and molecular methods for the detection of GFLV. These methods were selected to maximize GFLV detection by targeting viral particles (DAS-ELISA); encapsidated viral RNAs (IC-RT-PCR), viral RNAs (RT-qPCR and RNAseq) and vsiRNAs (siRNA blots and small RNAseq (Figure 1A). We also developed a bioinformatics framework to identify intra-lane cross contamination for GFLV.

<figure><img src="image_4.png" /></figure>
# MATERIALS AND METHODS

## Plant and Viral Material

Samples of young apical leaves were collected from 20 selected vines in two 30-year-old vineyards of *Vitis vinifera* cv. Chardonnay highly affected by fanleaf degeneration in the Champagne region (France) (**Supplementary Table S1**). Samples were collected in June when GFLV titer is high (Walter et al., 1984; Walter and Etienne, 1987) and flash frozen. Tissues were ground and homogenized into a fine powder with a mortar and pestle in liquid nitrogen. Aliquots of ground tissue were stored at -80°C and used for further testing by DAS-ELISA, IC-RT-PCR, RT-qPCR, vsiRNA Northern blot, and RNAseq.

Another set of six GFLV-infected grapevine leaf samples was chosen to compare the performance of RNAseq and smallRNAseq. Four samples (IC-MaA8191 and IC-MaA8193 from *V. vinifera* cv. Grenache, and IC-P1a and IC-P2a from *V. vinifera* cv. Chardonnay) were from the INRA grapevine virus collection in Colmar, France; and two samples (Va1 and Va4 from *V. vinifera* cv. Chardonnay) were from a commercial vineyard in Bagneux-la-Fosse in the Champagne region of France. Two samples (ENTAV-E39 and ENTAV-E173 from *V. vinifera* cv. 110R and Grenache, respectively) were from the grapevine collection in the Institut Français de la Vigne et du Vin (ENTAV-IFV, Le Grau-du-Roi, France). Additionally, leaf samples of grapevines infected by GFLV isolates GHu, F13-Col, B844, and CO1 (A17b) and IC-MaA8193 (**Supplementary Table S2**) from the INRA grapevine virus collection in Colmar, France were selected to verify the robustness of the diagnosis methods tested in this study. Leaves collected from these GFLV-infected grapevines were prepared in ELISA grinding buffer and conserved at -20°C until further testing. Sequences of these GFLV isolates and other GFLV isolates that were used as reference for bioinformatics purposes in this study are listed in **Supplementary Table S2**.

## Total RNA Extraction From Grapevine Leaf Samples and High Throughput Sequencing

Total RNAs were extracted from 100 mg of grapevine leaf tissue using an adapted protocol for the recovery of small and high molecular weight RNAs with the RNeasy Plant mini kit (Qiagen, Venlo, Netherlands), as per the manufacturer's recommendations. Post-extraction, purity criteria (A260nm/A230nm and A260nm/A280nm both > 1.8) and quality levels (RNA integrity number > 8) were assessed via a Bioanalyzer (Agilent, Santa Clara, CA, United States). Total RNA was used to prepare cDNA libraries after a poly-A selection at the GeT-PlaGe Genotoul platform facility (INRA-Toulouse, France). The HTS approach of choice in this study was a paired-end 2 × 150 pb RNAseq completed on a Hiseq 3000 (Illumina, San Diego, CA, United States) following the manufacturer's instructions. Double indexing was used for some samples (Lane 1, **Supplementary Table S1**) while single indexing was performed for other samples (Lane 2, **Supplementary Table S1**). For siRNA

sequencing, samples were sequenced on a Hiseq Instrument at 1 × 50 pb and multiplexed at Fasteris (Plan-les-Ouates, Switzerland).

## HTS Data Analyses

Demultiplexing was performed by GeT-Genotoul and Fasteris. The GeT-Genotoul platform used Bcl2fastq version 2.20.0.422, allowing 0 barcode mismatches with a minimum trimmed read length of 35. Analyses of datasets were performed using Workbench 11.0 software (CLC Genomics Workbench, Aarhus, Denmark). Prior to any downstream analyses, a FastQC was performed. For RNAseq datasets, only reads above 70 nucleotides (nts) were kept after trimming and quality check (**Supplementary Tables S3, S4**). For smallRNAseq datasets, only reads between 20 and 27 nt were kept for downstream analyses after a quality check step (**Supplementary Table S4**).

The sanitary status of grapevine samples was assessed by mapping reads to a collection of curated reference sequences of viruses known to infect grapevine (Martelli, 2017), as previously described (Hily et al., 2018c). Briefly, relaxed mapping stringency (0.5/0.7 or 0.9/0.9, corresponding to read length/similarity parameters for RNAseq and smallRNAseq data, respectively) was used to span maximum genome diversity within virus species. RPKM (Read Per Kilobase Million) was used as normalization method. RPKM values express abundance of viral RNAs, and are obtained by accounting for the number of reads mapping a virus genome, the genome size of the virus and the total number of reads from the sample. *De novo* assemblies (parameters of word size of 17 and minimum contig length of 200) were performed after removal of reads that mapped to the *Vitis vinifera* genome<fnref n="1" />. Contigs were then tested against GenBank reference sequences using Blastn/Blastx<fnref n="2" />.

For the detection of possible intra-lane contamination, all reads that mapped to the GFLV consensus RNA1 and RNA2 sequences at relaxed stringency parameters (see above) were recovered from each sample. Reads were then mapped to GFLV sequences obtained from the same sequencing lane (e.g., IC-MaA8193, IC-P1a, IC-P2a, KX034936, and KX034888) or GFLV sequences that were not present in the same lane [F13-col, B844, CO1 (A17b), CO2 (A17d), and GHu] (**Supplementary Tables S1, S2**). Mapping was performed at a very high stringency in order to confirm the origin of the reads [0.95/0.98 and 0.99/0.99 of (length/similarity) for RNAseq and smallRNAseq datasets, respectively]. A schematic representation of our protocol for detecting possible contamination is shown as **Supplementary Data Sheet 1**.

## RT-qPCR for GFLV Diagnosis

For viral detection and absolute quantification with RT real-time PCR, 20 ng of total RNA, quantified by Nanodrop™, were reverse-transcribed into cDNA using Superscript III reverse transcriptase (Invitrogen, Carlsbad, CA, United States) with a mix of oligo (dT). Eva Green qPCR mix (Biorad™) was used with

<footnote n="1">¹http://www.plantgdb.org/XGDB/phplib/download.php?GDB=Vv (Genoscope 12x, last visited March 2018)</footnote>

<footnote n="2">²http://blast.ncbi.nlm.nih.gov/Blast.cgi (last visited March 2018)</footnote>
specific primers targeting GFLV-RNA1 and RNA2 for real time amplification (Biorad™ CFX96 real-time system). Degenerated primers (**Supplementary Table S5**) were designed within the most conserved regions after alignment of 32 complete GFLV genomes (from NCBI and our own sequence database) for the amplification of about 120 bp fragments. PCR was performed in 96-well optical reaction plates for 30 s at 95°C, 40 cycles of denaturation for 5 s at 95°C, annealing and elongation for 5 s at 65°C. Each sample/primers combination was carried out in triplicate. A melting curve analysis was performed to ascertain that a single product was produced in each well. Viral quantification was expressed as molecule of viral RNA per ng of total RNA. The detection threshold of GFLV was determined via a serial dilution of plasmids containing cDNA from GFLV RNA1 and RNA2 (50 to 5.10<sup>-7</sup> ng of DNA) (Vigne et al., 2013).

## Semi-Quantitative DAS-ELISA for GFLV Diagnosis

DAS-ELISA was performed using a commercial kit (Bioreba AG, Reinach, Switzerland), as previously described (Vigne et al., 2013). For each sample, 100 µL of a mixture of frozen tissue powder (0.5 g) re-suspended in 5 mL of extraction buffer were loaded into microtiter plates, as well as serial dilution 10<sup>-1</sup> and 10<sup>-2</sup> for each sample. For quantification purpose, each microplate contained a serial dilution of a purified GFLV preparation, ranging from 7,480 to 0.22 pg of virion per well. GFLV concentration for each sample was deduced from a linear regression.

## IC-RT-PCR for GFLV Diagnosis

For each DAS-ELISA performed, a second plate was run in parallel for IC-RT-PCR. After antibody/antigen incubation, 12 µL of sterile water were added for 10 min at 70°C and 5 min on ice to disrupt viral particles. The recovered viral RNAs were subjected to first-strand cDNA synthesis with degenerated primers (**Supplementary Table S5**) to amplify fragments of 387 pb for RNA1 and 567 pb for RNA2. The cycling PCR parameters were an initial denaturation at 95°C for 2 min, followed by 36 cycles of 30 s at 95°C, 30 s at 52°C and 45 s at 72°C, followed by a 5 min at 72°C. For IC-RT-qPCR, after recovering the encapsidated viral RNAs, the aforementioned RT-qPCR protocol was applied.

## Northern Blot for GFLV vsiRNAs Detection

Total RNA was extracted from 0.1 g of frozen tissue using 1 ml Concert Plant RNA Reagent following manufacturer's protocol (Invitrogen, Carlsbad, CA, United States). Low molecular weight RNA was analyzed using 30 µg of total RNA separated by denaturing polyacrylamide gel electrophoresis and transferred to Hybond N + membranes (GE Healthcare, Piscataway, NJ, United States). Following chemically crosslinking with EDC treatment (Pall et al., 2007), hybridization was done with a mix of five radiolabeled DNA probes amplified from five genetically distinct GFLV isolates (F13-col, B844, CO1 (A17b), GHu, and IC-MaA8193) in siRNA hot spot regions

of RNA1 and RNA2 according to our smallRNAseq datasets (**Supplementary Table S5**). DNA probes labeled with α-32P dCTP (Hartmann, Braunschweig, Germany) were obtained using the Prime-a-Gene® labeling system (Promega, Madison, WI, United States). Hybridization was conducted overnight at 42°C in PerfectHyb™ plus buffer (Sigma-Aldrich, St Louis, MO, United States). Washes were done three times for 10 min at 50°C in 2X SSC and 2% SDS and the membranes were exposed to Biorad FX phosphorimager (Biorad, Hercules, CA, United States). To determine the detection threshold of vsiRNA blot, we used 21 and 24 nt RNA oligonucleotides specific to RNA1 or RNA2 (**Supplementary Table S5**). We performed a serial dilution to test the detection level for each RNA oligonucleotide. The loading control was hybridized with a vvi-miR159c probe (**Supplementary Table S5**).

## RESULTS

### Comparison of RNAseq and smallRNAseq

We tested and compared two HTS techniques to detect viruses in grapevine samples (**Figure 1**) with RNAseq (Illumina 3000, 2 × 150 bp) detecting mostly poly-A-tailed virus sequences and smallRNAseq (1 × 50 bp) detecting vsiRNA as the final product of the adaptive antiviral plant defense mechanism (Hamilton et al., 2002). Samples from four grapevines from the INRA grapevine collection in Colmar, France (three being GFLV-infected and one appearing as non-infected by GFLV) and from two grapevines showing typical fanleaf symptoms (stunting and very short internodes) in a Champagne vineyard were selected for this study. Two grapevine accessions provided by ENTAV-IFV were considered as negative controls. Our main goal was to detect GFLV and other viruses. It is important to mention that the eight samples were multiplexed in the same sequencing lane in the two HTS techniques tested (**Supplementary Table S1**).

### GFLV Detection

The same strategy was used to analyze the RNAseq and smallRNAseq datasets with most detection parameters being relaxed for the recovery of maximal GFLV diversity. Although GFLV is genetically diverse (Vigne et al., 2004; Mekuria et al., 2009; Zhou et al., 2017), we mapped all reads to RNA1 and RNA2 consensus sequences obtained from seven GFLV genomes (Hily et al., 2018c). As a result, many reads mapping to GFLV consensus sequences were recovered in all samples tested (**Table 1**) with 75 reads mapping to the GFLV-RNA1 consensus sequence for sample ENTAV-E39 and 183,385 reads mapping to the GFLV-RNA2 consensus sequence for sample IC-P1a for the RNAseq dataset; and with 254 (sample IC-MaA8191) to 797,497 (sample IC-MaA8193) reads for the smallRNAseq dataset (**Table 1**). The eight samples used for this study had either a 'high read count' (HRC) with elevated reads (67,540 to 797,497) mapping to the GFLV consensus sequences (highlighted in gray in **Table 1**) or a 'low read count' (LRC) with less than 1,000 reads recovered from either HTS techniques. Such distinct contrast in read recovery, ranging from 1 to 10,000 fold, was specific to
<table><caption>TABLE 1 | Information about reads mapping to viral sequences.</caption><thead><tr><th rowspan="2">Sample name</th><th rowspan="2">DAS-ELISA</th><th rowspan="2">Technique</th><th rowspan="2">Total read #</th><th colspan="8">GFLV consensus</th><th colspan="2">GFLV de novo sequence</th><th rowspan="2">GVA</th><th rowspan="2">GRSPaV</th></tr><tr><th>#reads</th><th>RNA1</th><th>RNA2</th><th>Highest depth</th><th>RNA1</th><th>RNA2</th><th>Coverage %</th><th>read</th><th>RNA1</th><th>RNA2</th><th>Coverage %</th><th># reads</th></tr></thead><tbody><tr><td>IC-MaA8191</td><td>-</td><td>small RNAseq</td><td>18,343,185</td><td>374</td><td>254</td><td>29</td><td>27</td><td>38.2</td><td>38.7</td><td>-</td><td>-</td><td>-</td><td>-</td><td>2</td><td>6,631</td></tr><tr><td></td><td></td><td>RNAseq</td><td>165,823,267</td><td>482</td><td>783</td><td>26</td><td>42</td><td>98.1</td><td>96.4</td><td>-</td><td>-</td><td>-</td><td>-</td><td>0</td><td>33,799</td></tr><tr><td>IC-MaA8193</td><td>+</td><td>small RNAseq</td><td>20,366,738</td><td>797,497</td><td>518,542</td><td>46,686</td><td>22,620</td><td>56.5</td><td>48.4</td><td>2,259,514</td><td>1,593,089</td><td>99.9</td><td>99.8</td><td>0</td><td>6,528</td></tr><tr><td></td><td></td><td>RNAseq</td><td>23,461,101</td><td>93,017</td><td>167,810</td><td>5,009</td><td>5,195</td><td>100.0</td><td>97.9</td><td>96,756</td><td>182,407</td><td>100.0</td><td>100.0</td><td>4</td><td>4,092</td></tr><tr><td>ENTAV-E39</td><td>-</td><td>small RNAseq</td><td>18,161,215</td><td>835</td><td>500</td><td>68</td><td>43</td><td>41.4</td><td>45.8</td><td>-</td><td>-</td><td>-</td><td>-</td><td>0</td><td>671</td></tr><tr><td></td><td></td><td>RNAseq</td><td>22,164,625</td><td>75</td><td>97</td><td>8</td><td>21</td><td>58.7</td><td>82.1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>0</td><td>14</td></tr><tr><td>ENTAV-E173</td><td>-</td><td>small RNAseq</td><td>23,226,315</td><td>468</td><td>329</td><td>37</td><td>26</td><td>41.6</td><td>38.1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td><td>3,038</td></tr><tr><td></td><td></td><td>RNAseq</td><td>51,911,501</td><td>243</td><td>265</td><td>20</td><td>24</td><td>87.2</td><td>93.5</td><td>-</td><td>-</td><td>-</td><td>-</td><td>0</td><td>3,030</td></tr><tr><td>IC-P1a</td><td>+</td><td>small RNAseq</td><td>27,468,283</td><td>711,941</td><td>396,013</td><td>54,488</td><td>68,134</td><td>81.5</td><td>71.7</td><td>1,879,692</td><td>1,233,456</td><td>99.7</td><td>99.6</td><td>1</td><td>3,561</td></tr><tr><td></td><td></td><td>RNAseq</td><td>29,519,363</td><td>155,869</td><td>183,385</td><td>8,374</td><td>14,409</td><td>99.8</td><td>97.6</td><td>165,675</td><td>204,211</td><td>100.0</td><td>100.0</td><td>0</td><td>3,176</td></tr><tr><td>IC-P2a</td><td>+</td><td>small RNAseq</td><td>22,989,799</td><td>591,793</td><td>381,518</td><td>45,195</td><td>64,694</td><td>82.1</td><td>73.3</td><td>1,844,236</td><td>1,097,753</td><td>99.5</td><td>99.8</td><td>6</td><td>11,097</td></tr><tr><td></td><td></td><td>RNAseq</td><td>22,461,740</td><td>67,540</td><td>101,207</td><td>4,724</td><td>7,753</td><td>99.8</td><td>97.6</td><td>68,942</td><td>112,932</td><td>100.0</td><td>100.0</td><td>0</td><td>5,498</td></tr><tr><td>Va1</td><td>-</td><td>small RNAseq</td><td>24,068,532</td><td>388</td><td>255</td><td>27</td><td>34</td><td>36.4</td><td>34.8</td><td>-</td><td>-</td><td>-</td><td>-</td><td>0</td><td>5,971</td></tr><tr><td></td><td></td><td>RNAseq</td><td>11,633,150</td><td>130</td><td>191</td><td>10</td><td>58</td><td>78.8</td><td>94.7</td><td>-</td><td>-</td><td>-</td><td>-</td><td>0</td><td>3,833</td></tr><tr><td>Va4</td><td>-</td><td>small RNAseq</td><td>29,827,959</td><td>545</td><td>345</td><td>37</td><td>38</td><td>42.5</td><td>42.4</td><td>-</td><td>-</td><td>-</td><td>-</td><td>4</td><td>11,620</td></tr><tr><td></td><td></td><td>RNAseq</td><td>8,596,909</td><td>127</td><td>187</td><td>12</td><td>70</td><td>81.5</td><td>94.1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>0</td><td>10,805</td></tr></tbody></table>

The number of reads (# reads) mapping to the consensus grapevine fanleaf virus (GFLV) RNA1 and RNA2 sequences obtained from seven full-length GFLV genomes (Hily et al., 2018c) are indicated. The highest depth corresponds to the highest number of reads present at any position along the genome. Coverage % corresponds to the coverage of the reference. In bold are numbers of reads from 'low read count' samples and numbers of reads from 'high read count' samples are in italic and highlighted in gray. For GFLV positive samples, de novo assembled GFLV sequences were used as reference (–: not applicable since no genome could be assembled). Sequences of grapevine virus A (GVA) and grapevine rupestris stem pitting-associated virus (GRSPaV) were used as controls.
GFLV; it was not observed when mapping other virus genomes. For example, 0 to 6 reads mapped to the genome of grapevine virus A (GVA; whose genome size is close to that of GFLV RNA1) and 671 to 11,620 reads mapped to the genome of grapevine rupestris stem pitting-associated virus (GRSPaV), an ubiquitous virus of grapevine (Table 1).

Sequencing depth was proportional to the number of reads mapped to each GFLV sequence. By focusing on the genome coverage percentage, close to the full-length GFLV genome sequence was covered by reads from the RNAseq dataset for HRC samples (97.6 to 100% genome coverage) and, to a lesser extent, for LRC samples (58.7–98.1% genome coverage). The difference in genome coverage between the two categories of samples was also observed with the smallRNAseq dataset but with a lower coverage (48.4–82.1 and 34.8–45.8% for HRC and LRC samples, respectively). While the number of reads was antithetical between HRC and LRC samples, their presence was undeniable, covering a high percentage of the GFLV genome. These results begged the following questions: Is the GFLV titer in LRC samples extremely low? Or, do these reads correspond to some contamination due to multiplexing? These two questions are of particular interest for the two ENTAV-IFV samples: ENTAV-E39 and ENTAV-E173 which are resulting from a successful therapeutics treatment via microshoot tip culture (Golino et al., 2017b), never displayed any disease symptom and always tested negative for major viral diseases of grapevine.

## A Protocol to Identify GFLV Intra-Lane Contamination

To discriminate between low GFLV titer or intra-lane contamination due to the sample multiplexing step, we

designed a framework of basic bioinformatics steps for RNAseq datasets (Figure 2). First, all reads mapping to GFLV consensus sequences at low stringency (parameters 0.5/0.7, length/identity, see section “Materials and Methods”) were recovered. Reads were then mapped to a series of beforehand *de novo* assembled GFLV genomes that were either present within the same sequencing lane (e.g., IC-P1a, IC-P2a and IMA-GMO) (Supplementary Tables S1, S2) or to GFLV sequences that were not present on the same sequencing lane (e.g., F13-col and B844) (Supplementary Table S2). For this step, very stringent parameters (0.95/0.98) were used (Table 2). Even at extremely high stringency and using only a small portion of GFLV sequences recovered from the same sequencing lane, 69–80% reads for GFLV RNA1 and 60–96% reads for GFLV RNA2 corresponded to these sequences (Table 2). In contrast, from the same pool of reads, only 0–2% of them mapped to sequences from GFLV variants not present in the sequencing lane, with the exception of 24% for GFLV RNA1 of sample IC-MaA819. This dichotomy in percentages indicates a probable source of contamination from samples being multiplexed within a particular sequencing lane. This was true also for smallRNAseq datasets although a higher percentage of reads mapping to GFLV variant sequences from either the same (88–99%) or different (26–52%) sequencing lanes was obtained. The higher percentage obtained with smallRNAseq compared to RNAseq datasets could be explained by the fact that smaller sequences (20–27 nt vs. ≈ 146 nt, Supplementary Table S4) allowed for more reads to be mapped onto a genome by the former compared to the latter technology.

Our ‘intra-lane contamination’ protocol was applied to all LRC samples. Consequently, the intra-lane GFLV contaminant

<figure><img src="image_5.png" /><figcaption>**FIGURE 2** | Framework for Grapevine fanleaf virus (GFLV) intra-lane contamination detection from RNAseq libraries. Images below are taken from the CLC Workbench software and correspond to steps within the framework specifically looking at GFLV RNA2. A more detailed protocol is provided in **Supplementary Data Sheet 1**.</figcaption></figure>
<table><caption>**TABLE 2** | Assessment of GFLV intra-lane contamination in RNAseq and smallRNAseq datasets.</caption><thead><tr><th rowspan="3">Sample name</th><th rowspan="3">read #</th><th rowspan="3">Technique</th><th colspan="3">RNA1</th><th colspan="3">RNA2</th></tr><tr><th>GFLV</th><th>Same lane</th><th>Other lane</th><th>GFLV</th><th>Same lane</th><th>Other lane</th></tr><tr><th>Consensus</th><th>Variants</th><th>Variants</th><th>Consensus</th><th>Variants</th><th>Variants</th></tr></thead><tbody><tr><td rowspan="2">IC-MaA8191</td><td>18,343,185</td><td>smallRNAseq</td><td>374</td><td>340 (91%)</td><td>121 (32%)</td><td>254</td><td>223 (88%)</td><td>67 (26%)</td></tr><tr><td>165,823,267</td><td>RNAseq</td><td>482</td><td>395 (80%)</td><td>115 (24%)</td><td>783</td><td>472 (60%)</td><td>14 (2%)</td></tr><tr><td rowspan="2">ENTAV-E39</td><td>18,161,215</td><td>smallRNAseq</td><td>835</td><td>816 (98%)</td><td>382 (46%)</td><td>500</td><td>480 (96%)</td><td>186 (37%)</td></tr><tr><td>22,164,625</td><td>RNAseq</td><td>75</td><td>60 (80%)</td><td>1 (1%)</td><td>97</td><td>63 (65%)</td><td>0 (0%)</td></tr><tr><td rowspan="2">ENTAV-E173</td><td>23,226,315</td><td>smallRNAseq</td><td>468</td><td>461 (99%)</td><td>220 (47%)</td><td>329</td><td>318 (97%)</td><td>121 (37%)</td></tr><tr><td>51,911,501</td><td>RNAseq</td><td>243</td><td>178 (73%)</td><td>5 (2%)</td><td>265</td><td>223 (84%)</td><td>2 (1%)</td></tr><tr><td rowspan="2">Va1</td><td>24,068,532</td><td>smallRNAseq</td><td>388</td><td>377 (97%)</td><td>203 (52%)</td><td>255</td><td>247 (97%)</td><td>106 (42%)</td></tr><tr><td>11,633,150</td><td>RNAseq</td><td>130</td><td>99 (76%)</td><td>2 (2%)</td><td>191</td><td>126 (66%)</td><td>0 (0%)</td></tr><tr><td rowspan="2">Va4</td><td>29,827,959</td><td>smallRNAseq</td><td>545</td><td>524 (96%)</td><td>270 (50%)</td><td>345</td><td>330 (96%)</td><td>141 (41%)</td></tr><tr><td>8,596,909</td><td>RNAseq</td><td>127</td><td>87 (69%)</td><td>3 (2%)</td><td>187</td><td>180 (96%)</td><td>0 (0%)</td></tr></tbody></table>

Reads that first mapped to GFLV-consensus sequences are indicated (Table 1). They were recovered and then mapped to de novo assembled sequences (RNA1 and RNA2) obtained from positive samples present on the same sequencing lane (same lane variants) or to sequences absent from that lane (other lane variants, Supplementary Tables S1, S2). Number of reads and percentage (in parenthesis) that mapped are indicated. Mapping parameters were set very stringently, 0.95/0.98 and 0.99/0.99 (length/identity) for RNAseq and smallRNAseq datasets, respectively.

reads were removed from the datasets and all LRC samples were ultimately defined as negative (Figure 2 and Table 1). Even if a contamination might be observed for other viruses (see GRSPaV in sample ENTAV-E39, Table 1), this protocol could not be used in our datasets for viruses other than GFLV since their accumulation was not sufficient enough to trigger the implementation of the 'intra-lane contamination' method. For example, the number of reads mapping to GRSPaV sequences corresponded to only 1.1-4% of those mapping GFLV sequences in samples IC-MaA8193, IC-P1a and IC-P2a (Table 1).

## Analysis of the Sanitary Status of Grapevine Samples: Detection of Viruses and Identification of Viral Variants

The robustness of RNAseq and smallRNAseq for the detection of viruses other than GFLV and viroids was tested. Two methods were applied (Hily et al., 2018c): a direct mapping using a set of references and a Blastn/Blastx analysis following a *de novo* assembly. For most samples, the same viral species were detected by both HTS techniques (Table 3). GRSPaV, as well as hop stunt viroid (HSVd) and grapevine yellow speckle viroid-1 (GYSVd1), were readily detected and assembled from all samples, except for sample ENTAV-E39 for which only HSVd was detected. Also, it was not surprising to see that a grapevine leafroll-associated virus 3 (GLRaV-3), which does not have a poly-A-tailed genome, was detected by smallRNAseq but not by RNAseq in sample IC-MaA8193 (Table 3 and Supplementary Figure S3). Nonetheless, some viruses with a non-poly-A-tailed genome were detected via RNAseq probably because the viral titer was high or the viral sequence might contain a stretch of adenine that was selected by the poly-(A) selection process (see GLRaV2 and GLRaV3 in samples Va4 and IC-MaA8191, respectively, Table 3). What was unexpected was that grapevine Red Globe virus (GRGV), a virus with a poly-A-tailed genome from the

family *Tymoviridae*, was not detected by RNAseq in sample Va4 but identified using smallRNAseq (Table 3 and Supplementary Figure S3).

From this study, it seems that smallRNAseq libraries allow for the detection of more viral species than the RNAseq technique. However, complete viral genome sequences were rarely obtained via *de novo* assembly when using the smallRNAseq datasets because many genome regions were not covered by reads. In contrast, it was very straightforward to obtain complete (to near complete) viral genomes using RNAseq datasets. Furthermore, we were able to assemble and identify different variants of a virus species (GRSPaV and GRVFV) within the same sample from either libraries (Table 3), however, *de novo* assembly was much easier from RNAseq compared to smallRNAseq datasets. No new viral species were detected in any of the samples selected for this study using the RNAseq or smallRNAseq datasets.

## RNAseq Versus Traditional GFLV Detection Methods

We compared the performance of RNAseq, DAS-ELISA, RT-qPCR, IC-RT-PCR and vsiRNA blot for GFLV diagnosis using grapevine leaf samples from the INRA grapevine virus collection in Colmar (France) and from two vineyards in the Champagne region of France. Five GFLV isolates from the INRA grapevine virus collection were selected to identify the most robust broad-spectrum diagnostic methodology. These five GFLV isolates have 82-92% nucleotide sequence identity among their RNA1 and RNA2 (Supplementary Figure S1).

## Traditional Methods for GFLV Diagnosis

First, we compared the sensitivity of each technics by testing serial dilutions of GFLV-infected grapevine leaf crude sap from five GFLV isolates. All isolates were readily detected by three methods (Table 4). However, IC-RT-qPCR that was designed to detect GFLV RNA2 was the most sensitive technique (1E-04). It was 10
<table><caption>**TABLE 3** | Virus detection in eighth grapevine samples by RNAseq and smallRNAseq.</caption><thead><tr><th rowspan="2">Sample name</th><th rowspan="2">Technique</th><th colspan="10">GFLV</th></tr><tr><th>RNA1</th><th>RNA2</th><th>GRSPaV</th><th>GYSVd1</th><th>HSVd</th><th>GFkV</th><th>GRVFV</th><th>GRGV</th><th>GLRaV2</th><th>GLRaV3</th></tr></thead><tbody><tr><td rowspan="2">IC-MaA8191</td><td>smallRNAseq</td><td>-</td><td>-</td><td>✓</td><td>✓</td><td>✓</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td></tr><tr><td>RNAseq</td><td>-</td><td>-</td><td>2</td><td>1</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>✓</td></tr><tr><td rowspan="2">IC-MaA8193</td><td>smallRNAseq</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1**</td></tr><tr><td>RNAseq</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>**</td></tr><tr><td rowspan="2">ENTAV-E39</td><td>smallRNAseq</td><td>-</td><td>-</td><td>-</td><td>-</td><td>✓</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>RNAseq</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="2">ENTAV-E173</td><td>smallRNAseq</td><td>-</td><td>-</td><td>✓</td><td>✓</td><td>✓</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>RNAseq</td><td>-</td><td>-</td><td>1</td><td>1</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="2">IC-Pa1</td><td>smallRNAseq</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td><td>-</td><td>✓</td><td>✓</td><td>-</td><td>-</td></tr><tr><td>RNAseq</td><td>1</td><td>1</td><td>2</td><td>2</td><td>1</td><td>-</td><td>2</td><td>✓</td><td>-</td><td>-</td></tr><tr><td rowspan="2">IC-P2a</td><td>smallRNAseq</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td><td>✓</td><td>-</td><td>-</td><td>✓</td><td>-</td><td>-</td></tr><tr><td>RNAseq</td><td>1</td><td>1</td><td>4</td><td>1</td><td>1</td><td>-</td><td>-</td><td>✓</td><td>-</td><td>-</td></tr><tr><td rowspan="2">Va1</td><td>smallRNAseq</td><td>-</td><td>-</td><td>✓</td><td>✓</td><td>✓</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>RNAseq</td><td>-</td><td>-</td><td>2</td><td>2</td><td>1</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="2">Va4</td><td>smallRNAseq</td><td>-</td><td>-</td><td>✓</td><td>✓</td><td>✓</td><td>1</td><td>-</td><td>-</td><td>1*</td><td>1</td></tr><tr><td>RNAseq</td><td>-</td><td>-</td><td>3</td><td>1</td><td>1</td><td>1</td><td>-</td><td>-*</td><td>✓</td><td>-</td></tr></tbody></table>

Names of viruses detected (✓) and number of genomes assembled (# and highlighted in gray). *Sample Va4 showed 281,643 reads for grapevine Red Globe virus (GRGV) in the smallRNAseq datasets and 118 reads (23.6% genome covered) in RNAseq datasets. See **Supplementary Figure S3** for more information. **Sample IC-MaA8193 showed 60,875 reads for grapevine leafroll-associated virus 3 (GLRaV3) in the smallRNAseq datasets and 10 reads (4.71% genome covered) in the RNAseq datasets. See **Supplementary Figure S3** for more information.

times more sensitive than IC-RT-PCR which was also ten times more sensitive than DAS-ELISA, except for variant IC-MaA8193.

Then we assessed the detection limit for each method selected in this study. For the antibody-based detection techniques (DAS-ELISA and IC-RT-PCR), we started from a known quantity of a purified GFLV strain GHu preparation from which serial dilutions were performed from 7,480 to 0.22 pg. For DAS-ELISA, the detection limit was set at 117 pg of virus (Figure 3A, higher part of the panel), which corresponded to the dilution above the background set at 2.5× the optical density (OD405nm) of the negative control. From the same GFLV purification, 1.75 pg of virus was detected by IC-RT-PCR when using RNA1 specific primers and 0.44 pg of GFLV was detected when using RNA2 specific primers (Figure 3A, lower panel). This result indicated that IC-RT-PCR is 60–250x more sensitive than DAS-ELISA for GFLV diagnosis. In order to define the RT-qPCR-based detection limit, we used known amounts of plasmids carrying cDNAs of GFLV RNA1 or RNA2, from which serial dilutions were made.

The last dilution detected by qPCR was 5 fg for both plasmids which corresponds to 480 molecules of RNA1 and 780 molecules of RNA2 (Figure 3B). To estimate the vsiRNA blot detection limit, small RNAs of 21 and 24 nt corresponding to conserved sequences within the five probes used for hybridization were designed. As low as 0.1 pmol of siRNA RNA1 was detected by vsiRNA blot (Figure 3C).

## Field Samples Study

To compare the different techniques in a real detection setting, twenty samples from two GFLV-infected vineyard plots located in the Champagne region of France were evaluated. Thirteen out of the 20 samples were positive for GFLV by DAS-ELISA while seven samples were negative (Table 5). Viral quantity ranged from 1,497 to 24,712 ng of virion per gram of fresh leaf material in infected samples. An immuno-capture assay followed by RT-PCR confirmed the aforementioned results, with the same 13 samples testing positive for GFLV and the same seven samples

<table><caption>**TABLE 4** | Sensitivity of three diagnostic techniques for five GFLV isolates.</caption><thead><tr><th rowspan="2">GFLV isolates</th><th rowspan="2">DAS-ELISA</th><th colspan="2">IC-RT-PCR</th><th colspan="2">IC-RT-qPCR</th></tr><tr><th>RNA1</th><th>RNA2</th><th>RNA1</th><th>RNA2</th></tr></thead><tbody><tr><td>B844</td><td>1E-02</td><td>1E-03</td><td>1E-03</td><td>1E-04 (38.6)</td><td>1E-04 (36.5)</td></tr><tr><td>F13</td><td>1E-02</td><td>1E-03</td><td>1E-03</td><td>1E-03 (37.3)</td><td>1E-04 (36.2)</td></tr><tr><td>GHu</td><td>1E-02</td><td>1E-03</td><td>1E-04</td><td>1E-04 (39.5)</td><td>1E-04 (36.8)</td></tr><tr><td>C01 (A17b)</td><td>1E-02</td><td>1E-03</td><td>1E-03</td><td>1E-03 (37.9)</td><td>1E-04 (34.3)</td></tr><tr><td>MaA8193</td><td>1E-03</td><td>1E-01</td><td>1E-03</td><td>1E-04 (37.0)</td><td>1E-04 (35.2)</td></tr></tbody></table>

Sensitivity is expressed as the last dilution of crude leaf sap for which a positive reaction is obtained, with the first dilution being 0.5 g of frozen tissue in 5 mL of DAS-ELISA grinding buffer. For IC-RT-qPCR, C<sub>t</sub>-values at which the virus is detected are indicated in parenthesis.
<figure><img src="image_6.png" /></figure>

testing negative using either RNA1 or RNA2 specific primers (Table 5 and Supplementary Figure S2A). The same 13 samples were confirmed positive for GFLV by RT-qPCR, except sample 1 for which RNA1 was below the detection level (Table 5), while RNA2 was unquestionably detected in that particular sample. All seven samples found GFLV negative by the two antibody-based techniques were confirmed negative by RT-qPCR by using primers for both genomic RNAs. Northern blot for vsiRNA validated the results obtained by the other detection methods, with the same 13 samples displaying vsiRNA, while seven did not (Table 5 and Supplementary Figure S2B).

We finally performed RNAseq with total RNAs used in RT-qPCR. The twenty samples were multiplexed in a single RNAseq lane and an average of about 34 million raw reads (ranging from 31 to 44 M, Supplementary Table S3) was obtained per sample. Interestingly, all samples displayed reads

corresponding to the GFLV genome, ranging from 99 to 379,611 reads and 127 to 506,907 reads mapping to the GFLV consensus RNA1 and RNA2 sequences, respectively (Supplementary Table S6). However, treating the data with our intra-lane contamination detection protocol, the seven samples that were found negative using the previously ran methodologies were negative for GFLV, showing consistency among all the results (Table 5). However, for GFLV quantification, no correlation was found among diagnostic techniques, as expected (Supplementary Table S8).

From the RNAseq libraries, 26 GFLV RNA1, 18 GFLV RNA2, and three GFLV RNA3 complete to near complete (i.e., covering at the minimum the open reading frame) sequences were *de novo* assembled (Supplementary Table S7), underlying the fact that 10 grapevines were co-infected with more than one RNA1 molecule and five of them with two RNA2 molecules.
<table><caption>**TABLE 5** | Grapevine fanleaf virus detection in 20 vineyard samples by five techniques.</caption><thead><tr><th rowspan="2"></th><th colspan="14">Sample names</th></tr><tr><th>Pa1</th><th>Pa2</th><th>Pa3</th><th>Pa4</th><th>Pa5</th><th>Pa6</th><th>Pa7</th><th>Pa8</th><th>Pa9</th><th>Pa10</th><th>Pa11</th><th>Pa12</th><th>Pa13</th><th>Pa14</th><th>Pa15</th><th>Pa16</th><th>Py17</th><th>Py18</th><th>Py19</th><th>Py20</th></tr></thead><tbody><tr><td>DAS-ELISA (ng/g)</td><td>3,104</td><td>6,782</td><td>-</td><td>7,004</td><td>-</td><td>1,497</td><td>-</td><td>4,531</td><td>-</td><td>4,372</td><td>-</td><td>6,581</td><td>-</td><td>23,654</td><td>4,425</td><td>-</td><td>24,712</td><td>18,077</td><td>4,393</td><td>6,518</td></tr><tr><td>IC-RT-PCR</td><td>RNA1</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td></td><td>RNA2</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>RT-qPCR (mol/ng)</td><td>RNA1</td><td>-</td><td>94,258</td><td>-</td><td>115,494</td><td>-</td><td>68,575</td><td>-</td><td>12,321</td><td>-</td><td>10,944</td><td>-</td><td>60,934</td><td>-</td><td>59,227</td><td>52,485</td><td>-</td><td>13,132</td><td>11,489</td><td>4,776</td><td>16,898</td></tr><tr><td></td><td>RNA2</td><td>40,369</td><td>200,184</td><td>-</td><td>267,744</td><td>-</td><td>25,6424</td><td>-</td><td>26,717</td><td>-</td><td>91,501</td><td>-</td><td>94,271</td><td>-</td><td>193,214</td><td>174,371</td><td>-</td><td>219,388</td><td>159,640</td><td>166,384</td><td>124,855</td></tr><tr><td>vsiRNA blot</td><td>RNA1</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td></td><td>RNAseq (RPKM)</td><td>363</td><td>925</td><td>-</td><td>669</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td></td><td>RNA2</td><td>668</td><td>2,097</td><td>-</td><td>1,249</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td></td><td>RNA3</td><td>-</td><td>925</td><td>-</td><td>900</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr></tbody></table>

DAS-ELISA: ng of virions per gram of leaves (ng/g). RT-qPCR: number of RNA molecule per ng of total RNA (mol/ng) and RNAseq: reads per kilobase Million (RPKM). For IC-RT-PCR and vsiRNA, see Supplementary Figure S2. Negative samples (-) with a reaction below the detection limit (117 pg for ELISA, 480 RNA1 and 780 RNA2 molecules for RT-qPCR, no visible amplicons in IC-RT-PCR nor vsiRNA on the blots) and a negative outcome following our intra-lane contamination protocol for RNAseq are indicated.

## Additional Information on the Virome in Two Selected Vineyard Plots in the Champagne Region of France

From the RNAseq datasets, information on the presence of viroids and other poly-A-tailed virus genomes was obtained (Supplementary Table S7). GRSPaV, HSVd and GYSVd1 were readily detected and assembled from all samples. A majority of samples (70%) was also infected with grapevine fleck virus (GFkV). Grapevine rupestris vein feathering virus (GRVFV) and grapevine Syrah virus 1 (GSyV-1) were also easily detected. Sample Pa10 was infected with grapevine Pinot gris virus (GPGV). Additionally, the presence of grapevine leafroll-associated virus 2 (GLRaV-2) with a non-poly-A-tailed genome was unquestionably detected in samples Pa4, Pa8, Py14, and Py17, although the full-length viral genome was not obtained. The presence of these viruses was corroborated via RT-PCR: GPGV and GRVFV in sample Pa10, GSyV-1 in sample Pa5 and GLRaV2 in Py14 (data not shown). These results confirm that grapevine is rarely infected by a single virus species, or even a single variant of that species, but rather by multiple viruses in a vineyard. Only samples Pa3 and Py13 were only co-infected with GRSPaV and viroids that are present in most grapevines worldwide and are considered as part of the 'microbiome' of grapevine (Saldarelli et al., 2017). All the other samples were super-infected with at least one (samples Pa7, Pa9, Py16, and Py19) and up to five different viral species (sample Py17). Nonetheless, the virus count might not be final, knowing that RNAseq does not provide an exhaustive view of the virome.

## DISCUSSION

In this work we compared RNAseq with semi-quantitative DAS-ELISA, IC-RT-PCR, RT-qPCR and vsiRNA blot for the detection of GFLV. Our results showed that GFLV is readily detected by the five diagnostic methodologies (Figure 1). These techniques covered the viral diversity exemplified by five GFLV strains (Supplementary Figure S1 and Table 4). Diagnostic results were consistent across the five methods, from the very cheap and dependable DAS-ELISA to the massive parallel sequencing strategy offered by RNAseq. The reliability of DAS-ELISA was reported in previous studies with only rare cases of GFLV variants not being detected (Liebenberg et al., 2009; Bruisson et al., 2017). Our findings also confirmed that IC-RT-qPCR has a better detection sensitivity compared to IC-RT-PCR which was also better than that of DAS-ELISA (Table 4). Interestingly, detecting the presence of GFLV was more sensitive when targeting RNA2 compared to RNA1. Likewise, when dealing with RNAseq datasets, RPKMs (which reflect viral accumulation) were always higher for RNA2 compared to RNA1 sequences (Table 5). This higher accumulation of RNA2 was previously described from the estimation of the quantity of RNA species from purified GFLV particles (Pinck et al., 1988).

As expected, no correlation in GFLV quantification was observed between the different diagnostic techniques (Supplementary Table S8). What was more startling was
that no correlation in viral quantification was found between techniques that identify the same target, for example total RNA by RT-qPCR and RNAseq. This could be due to an inherent bias of each method such as the use of RPKM that normalize the coverage/quantification along the whole molecule for RNAseq compared to the involvement of anchored primers in the RT-qPCR method, amplifying only one small targeted portion of the viral genome. The use of specific primers might also prevent the amplification of all variants by RT-qPCR due to mismatches. However, correlations were found when comparing GFLV RNA1 and RNA2 quantification by the same detection technique, strongly supporting the fact that both RNAs are necessary for a successful infection to take place and probably at a steady equilibrium, with a RNA2/RNA1 ratio = 2.37 ± 0.16 according to our RNAseq dataset. This ratio was quite similar to the 2.33 ± 0.06 ratio from one of our previous metagenomics study (Hily et al., 2018c). This result reinforces the need to target RNA2 for a more efficient GFLV detection by molecular techniques.

High throughput sequencing offers unparalleled amount of sequence information from a single run. From this work using either RNAseq or smallRNAseq libraries, we learned that the number of reads or the genome coverage are not sufficient to proclaim a sample positive for a virus (Table 1). After a direct mapping analysis from the RNAseq datasets, both HRC and LRC samples displayed over 80% coverage of the GFLV genome. By performing a simple *de novo* assembly from reads mapping to GFLV-consensus sequences in the LRC samples, not a single complete GFLV sequence was obtained but rather a number of short contigs that matched complete GFLV genomes *de novo* assembled from HRC samples present in the same sequencing lane. From these observations, we developed a bioinformatics protocol to remove highly identical reads from different samples within the same sequencing lane (Figure 2 and Supplementary Data Sheet 1), allowing for a discrimination between potential low viral titer infection and intra-lane cross-contamination, also known as ‘index misassignment’ or ‘index hopping’ (Kircher et al., 2012). Such contamination arises when multiple libraries are pooled in a single lane and sequenced simultaneously. To address this issue, a unique sequence (called index) per sample is added to individual DNA fragments during a library preparation. Index hopping occurs post-sequencing and before final data analysis during a process called ‘demultiplexing.’ Levels of index hooping are different depending on the indexes being used<fnref n="3" />. This step consists in identifying all reads by sorting them computationally to its correct library. While index hopping is inevitable when multiplexing samples, few practices offer mitigation options, with the best of them consisting in the use of a unique dual indexing combinations (Kircher et al., 2012). However, we clearly show that the use of dual indexing is not sufficient to completely abolish such index hopping (see results from Lane 1). This is why other options to identify and evaluate contaminations are needed. Our contamination evaluation

protocol is quite robust, however, it will not discriminate contaminations generated prior to sequencing, for example, contamination produced during nucleic acid extraction, through the addition of indexes and PCR steps. Furthermore, one main drawback of our procedure is when a LRC sample is co-infected with a virus at low titer and a virus with the same (or very close) genome sequence to the HRC one(s) pooled in the same sequencing lane. In that particular case, it would be technically impossible to discriminate between contamination and low viral titer.

With this study, we confirmed that more than just a detection tool, these massively paralleled sequencing techniques can detect and discriminate, within a sample and from a single run, one or more variants of the same virus or of different viral species (Table 3, Supplementary Figure S3 and Supplementary Table S7; Beuve et al., 2018; Hily et al., 2018a,c). Such feat is not possible using traditional detection technics unless, for example, the PCR-based techniques have been beforehand specially designed for this purpose, requiring additional tune-up and development prior to the analysis. Other information can be mined when using HTS, focusing not only on a single virus but describing the full virome of a plant. This terminology, describing the collection of nucleic acid that make up the viral community within a sample, takes all its value when dealing with a perennial plant such as grapevine. In this work, we confirmed that the virome of a vineyard is often more complex than previously thought, with grapevines being simultaneously infected by many viral species (Table 3 and Supplementary Table S7; Coetzee et al., 2010; Eichmeier et al., 2016; Fajardo et al., 2017; Beuve et al., 2018). HTS studies previously revealed the presence of many mycoviruses (Al Rwahnih et al., 2011; Espach et al., 2012)<fnref n="4" />, as well as viruses with unknown host (Hily et al., 2018b) but associated with grapevine. All aforementioned reports are describing the complex interaction between microbial communities within a particular ecosystem, ultimately questioning the concept and the definition of a ‘healthy’ vine. Here, we confirm the presence of several virus species of the family *Tymoviridae* such as GFkV, GSyV-1, and GRVFV, in samples from two vineyards in the Champagne region of France. These might ultimately be considered as GRSPaV, HSVd and GYSVd to be part of the ‘background’ virome of grapevine (Saldarelli et al., 2017).

Our data highlighted some advantages and drawbacks of RNAseq and smallRNAseq libraries. For example, the latter offers a wider viral detection panel, not restricted to poly-A-tailed viruses (Table 3). On the other hand, *de novo* assembly and full-length genome sequence recovery was more difficult than from RNAseq libraries. Often some viral genomes were not fully recovered from smallRNAseq datasets, with small gaps along the genome (data not shown). This was also previously reported for other sequencing techniques (Beuve et al., 2018). Except for GFkV or other viruses with a cytosine-rich genome, RNAseq provided data that can yield complete viral genomes as a result

<footnote n="3">³https://www.illumina.com/content/dam/illumina-marketing/documents/products/whitepapers/index-hopping-white-paper-770-2017-004.pdf (last visited May 2018)</footnote>

<footnote n="4">⁴http://www.wineland.co.za/virus-diversity-associated-grapevine-leafroll-disease/ (last visited May 2018)</footnote>
of a straightforward *de novo* assembly protocol. Moreover, viral genetic diversity was better defined using RNAseq. For example, within our parameters (quality and sequencing depth), little to no difference in reads mapping were observed between GFLV-consensus and *de novo* assembled sequences from our RNAseq datasets. This indicated that using a consensus genome might be the best approach to recover a maximum viral diversity within a sample even when dealing with a highly divergent virus such as GFLV. On the other hand, read numbers nearly tripled when dealing with smallRNAseq data. More work in parameter setting is needed to easily detect the full genetic diversity within a sample when using smallRNAseq libraries. Another advantage that HTS offers compared to traditional methods is the possibility to go back to datasets which can be considered as ‘evergreen’ and can be analyzed endlessly. This is particularly advantageous because new viruses are detected almost monthly.

All detection techniques were in agreement for the characterization of the presence or absence of GFLV in samples from a grapevine virus collection or naturally infected vineyards. We can conclude that (i) DAS-ELISA is a robust and reliable method detecting as low as 117 pg of virion, (ii) Molecular approaches are, as expected, more sensitive than serological methods, regardless of the isolate tested with IC-RT-qPCR allowing the detection of GFLV as low as 0.44 pg of virions when targeting RNA2, (iii) SmallRNAseq libraries are useful for an exhaustive definition of the sanitary status of a sample, (iv) RNAseq is the method of choice for an easy detection of distinct GFLV variants, and (v) a suitable protocol discriminating between contamination and the potential low GFLV infection level needs to be implemented for RNAseq and small RNAseq datasets to confidently decide on the infectious status of grapevine samples.

## DATA AVAILABILITY

The raw data supporting the conclusions of this manuscript will be made available by the authors to any qualified researcher upon request.

## 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

Al Rwahnih, M., Daubert, S., Golino, D., Islas, C., and Rowhani, A. (2015). Comparison of next-generation sequencing versus biological indexing for the optimal detection of viral pathogens in grapevine. *Phytopathology* 105, 758–763. doi: 10.1094/PHYTO-06-14-0165-R

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

Al Rwahnih, M., Daubert, S., Urbez-Torres, J. R., Cordero, F., and Rowhani, A. (2011). Deep sequencing evidence from single grapevine plants reveals a virome dominated by mycoviruses. *Arch. Virol.* 156, 397–403. doi: 10.1007/s00705-010-0869-8

## AUTHOR CONTRIBUTIONS

EV, OL, and J-MH designed the experiments. EV, SG, VK, and J-MH performed the experiments. EV, OL, and J-MH wrote the paper. All authors contributed to the manuscript revision, read and approved the submitted version.

## FUNDING

This work was supported by the Institut National de la Recherche Agronomique (INRA) and the project VACCIVINE funded by the ‘Plan National Dépérissement du Vignoble’ (French Ministry of Agriculture, FranceAgrimer and CNIV), class of 2017. Additional funds were kindly provided by three French professional committees for viticulture (Interprofessional Committee of Champagne Wine, CIVC, Epernay; Interprofessional Office of Burgundy Wines, BIVB, Beaune; and Interprofessional Committee of Wines from Alsace, CIVA, Colmar), the Chamber of Agriculture of Vaucluse and Moët and Chandon.

## ACKNOWLEDGMENTS

We would like to particularly recognize Géraldine Uriel (CIVC), Véronique Bonnet (Moët and Chandon), Olivier Jacquet (the Chamber of Agriculture of Vaucluse), and Pascal Blois (IFV) for their help and support in the field and for providing samples used in this study. We thank Dr. Marc Fuchs for proofreading the manuscript. We acknowledge J. Misbach and the greenhouse team for technical support, L. Ley and the members of the experimental unit of INRA-Colmar for the production of plants.

## SUPPLEMENTARY MATERIAL

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

Andret-Link, P., Laporte, C., Valat, L., Ritzenthaler, C., Demangeat, G., Vigne, E., et al. (2004). Grapevine fanleaf virus: still a major threat to the grapevine industry. *J. Plant Pathol.* 86, 183–195.

Barzon, L., Lavezzo, E., Costanzi, G., Franchin, E., Toppo, S., and Palù, G. (2013). Next-generation sequencing technologies in diagnostic virology. *J. Clin. Virol.* 58, 346–350. doi: 10.1016/j.jcv.2013.03.003

Beuve, M., Hily, J. M., Alliaume, A., Reinbold, C., Le Maguet, J., Candresse, T., et al. (2018). A complex virome unveiled by deep sequencing analysis of RNAs from a French pinot noir grapevine exhibiting strong leafroll symptoms. *Arch. Virol.* 163, 2937–2946. doi: 10.1007/s00705-018-3949-9

Blouin, A. G., Ross, H. A., Hobson-Peters, J., O'Brien, C. A., Warren, B., and MacDiarmid, R. (2016.). A new virus discovered by immunocapture of double-stranded RNA, a rapid method for virus enrichment in metagenomic studies. *Mol. Ecol. Res.* 16, 1255–1263. doi: 10.1111/1755-0998.12525

Bruisson, S., Lebel, S., Walter, B., Prevotat, L., Seddas, S., and Schellenbaum, P. (2017). Comparative detection of a large population of grapevine viruses by TaqMan RT-qPCR® and ELISA. *J. Virol. Methods* 240, 73–77. doi: 10.1016/j.jviromet.2016.12.003
Coetzee, B., Freeborough, M.-J., Maree, H. J., Celton, J.-M., Rees, D. J. G., and Burger, J. T. (2010). Deep sequencing analysis of viruses infecting grapevines: virome of a vineyard. *Virology* 400, 157–163. doi: 10.1016/j.virol.2010.01.023

Dayaram, A., Opong, A., Jäschke, A., Hadfield, J., Baschiera, M., Dobson, R. C. J., et al. (2012). Molecular characterisation of a novel cassava associated circular ssDNA virus. *Virus Res.* 166, 130–135. doi: 10.1016/j.virusres.2012.03.009

Donaire, L., Wang, Y., Gonzalez-Ibeas, D., Mayer, K. F., Aranda, M. A., and Llave, C. (2009). Deep-sequencing of plant viral small RNAs reveals effective and widespread targeting of viral genomes. *Virology* 392, 203–214. doi: 10.1016/j.virol.2009.07.005

Eichmeier, A., Komínková, M., Komínek, P., and Baránek, M. (2016). Comprehensive virus detection using next generation sequencing in grapevine vascular tissues of plants obtained from the wine regions of Bohemia and Moravia (Czech Republic). *PLoS One* 11:e0167966. doi: 10.1371/journal.pone.0167966

Espach, Y., Maree, H. J., and Burger, J. T. (2012). "The use of next-generation sequencing to identify novel mycoviruses in single grapevine plants," in *Proceedings of the 17th International Council for the Study of Virus and Virus-like Diseases of the Grapevine*, (Davis, CA: ICVG).

Fajardo, T. V. M., Silva, F. N., Eiras, M., and Nickel, O. (2017). High-throughput sequencing applied for the identification of viruses infecting grapevines in Brazil and genetic variability analysis. *Trop. Plant Pathol.* 42, 250–260. doi: 10.1007/s40858-017-0142-8

Filloux, D., Dallot, S., Delaunay, A., Galzi, S., Jacquot, E., and Roumagnac, P. (2015). "Metagenomics approaches based on virion-associated nucleic acids (VANA): an innovative tool for assessing without a priori viral diversity of plants," in *Plant Pathology. Methods in Molecular Biology*, ed. C. Lacomme (New York, NY: Humana Press).

Golino, D. A., Fuchs, M., Al Rwahnih, M., Farrar, K., Schmidt, A., and Martelli, G. P. (2017a). "Regulatory aspects of grape viruses and virus diseases: certification, quarantine, and harmonization," in *Grapevine Viruses: Molecular Biology, Diagnostics and Management*, eds B. Meng, G. P. Martelli, D. A. Golino, and M. Fuchs (Cham: Springer International Publishing), 581–598. doi: 10.1007/978-3-319-57706-7_28

Golino, D. A., Fuchs, M., Sim, S., Farrar, K., and Martelli, G. P. (2017b). "Improvement of grapevine planting stock through sanitary selection and pathogen elimination," in *Grapevine Viruses: Molecular Biology, Diagnostics and Management*, eds B. Meng, G. P. Martelli, D. A. Golino, and M. Fuchs (Cham: Springer International Publishing), 561–579.

Hamilton, A., Voinnet, O., Chappell, L., and Baulcombe, D. (2002). Two classes of short interfering RNA in RNA silencing. *EMBO J.* 21, 4671–4679. doi: 10.1093/embj/cdf464

Hemmer, C., Djennane, S., Ackerer, L., Hleibieh, K., Marmonier, A., Gersch, S., et al. (2018). Nanobody-mediated resistance to *Grapevine fanleaf* virus in plants. *Plant Biotechnol. J.* 16, 660–671.

Hily, J. M., Beuve, M., Vigne, E., Demangeat, G., Candresse, T., and Lemaire, O. (2018a). A genome-wide study of grapevine rupestris stem pitting-associated virus. *Arch. Virol.* 163, 3105–3111. doi: 10.1007/s00705-018-3945-0

Hily, J. M., Candresse, T., Garcia, S., Vigne, E., Tannières, M., Komar, V., et al. (2018b). Next-Generation Sequencing and the viromic study of grapevine leaves: from detection of grapevine-infecting viruses to description a new environmental *Tymovirales. Front. Microbiol.* 9:1782. doi: 10.3389/fmicb.2018.01782

Hily, J. M., Demanèche, S., Poulicard, N., Tannières, M., Djennane, S., Beuve, M., et al. (2018c). Metagenomic-based impact study of transgenic grapevine rootstock on its associated virome and soil bacteriome. *Plant Biotechnol. J.* 16, 208–220. doi: 10.1111/pbi.12761

Kircher, M., and Kelso, J. (2010). High-throughput DNA sequencing – concepts and limitations. *Bioessays* 32, 524–536. doi: 10.1002/bies.200900181

Kircher, M., Sawyer, S., and Meyer, M. (2012). Double indexing overcomes inaccuracies in multiplex sequencing on the illumina platform. *Nucleic Acids Res.* 40:e3. doi: 10.1093/nar/gkr771

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

Liebenberg, A., Freeborough, M.-J., Visser, C. J., Bellstedt, D. U., and Burger, J. T. (2009). Genetic variability within the coat protein gene of *Grapevine*

*fanleaf* virus isolates from South Africa and the evaluation of RT-PCR, DAS-ELISA and immunostrips as virus diagnostic assays. *Virus Res.* 142, 28–35. doi: 10.1016/j.virusres.2009.01.016

Maree, H. J., Fox, A., Al Rwahnih, M., Boonham, N., and Candresse, T. (2018). Application of HTS for routine plant virus diagnostics: state of the art and challenges. *Front. Plant Sci.* 9:1082. doi: 10.3389/fplms.2018.01082

Martelli, G. P. (2017). "An overview on grapevine viruses, viroids, and the diseases they cause," in *Grapevine Viruses: Molecular Biology, Diagnostics and Management*, eds B. Meng, G. P. Martelli, D. A. Golino, and M. Fuchs (Cham: Springer International Publishing), 31–46. doi: 10.1007/978-3-319-57706-7_2

Mekuria, T. A., Gutha, L. R., Martin, R. R., and Naidu, R. A. (2009). Genome diversity and intra- and interspecies recombination events in *Grapevine fanleaf* virus. *Phytopathology* 99, 1394–1402. doi: 10.1094/PHYTO-99-12-1394

Pall, G. S., Codony-Servat, C., Byrne, J., Ritchie, L., and Hamilton, A. (2007). Carbodiimide-mediated cross-linking of RNA to nylon membranes improves the detection of siRNA, miRNA and piRNA by northern blot. *Nucleic Acids Res.* 35:e60. doi: 10.1093/nar/gkm112

Pinck, L., Fuchs, M., Pinck, M., Ravelonandro, M., and Walter, B. (1988). A satellite RNA in grapevine fanleaf virus strain F13. *J. Gen. Virol.* 69, 233–239. doi: 10.1099/0022-1317-69-1-233

Reuter, J. A., Spacek, D., and Snyder, M. P. (2015). High-throughput sequencing technologies. *Mol. Cell.* 58, 586–597. doi: 10.1016/j.molcel.2015.05.004

Roossinck, M. J., Martin, D. P., and Roumagnac, P. (2015). Plant virus metagenomics: advances in virus discovery. *Phytopathology* 105, 716–727. doi: 10.1094/PHYTO-12-14-0356-RVW

Saldarelli, P., Giampetruzzi, A., Maree, H. J., and Al Rwahnih, M. (2017). "High-throughput sequencing: advantages beyond virus identification," in *Grapevine Viruses: Molecular Biology, Diagnostics and Management*, eds B. Meng, G. P. Martelli, D. A. Golino, and M. Fuchs (Cham: Springer International Publishing), 625–642. doi: 10.1007/978-3-319-57706-7_30

Sanfaçon, H. (2015). *Secoviridae: A Family of Plant Picorna-Like Viruses with Monopartite or Bipartite Genomes*. Hoboken, NJ: John Wiley & Sons, Ltd.

Schmitt-Keichinger, C., Hemmer, C., Berthold, F., and Ritzenthaler, C. (2017). "Molecular, Cellular, and Structural Biology of Grapevine fanleaf virus," in *Grapevine Viruses: Molecular Biology, Diagnostics and Management*, eds B. Meng, G. P. Martelli, D. A. Golino, and M. Fuchs (Cham: Springer International Publishing), 83–107. doi: 10.1007/978-3-319-57706-7_4

Vigne, E., Bergdoll, M., Guyader, S., and Fuchs, M. (2004). Population structure and genetic variability within isolates of *Grapevine fanleaf* virus from a naturally infected vineyard in France: evidence for mixed infection and recombination. *J. Gen. Virol.* 85, 2435–2445. doi: 10.1099/vir.0.79904-0

Vigne, E., Gottula, J., Schmitt-Keichinger, C., Komar, V., Ackerer, L., Belval, L., et al. (2013). A strain-specific segment of the RNA-dependent RNA polymerase of grapevine fanleaf virus determines symptoms in *Nicotiana* species. *J. Gen. Virol.* 94, 2803–2813. doi: 10.1099/vir.0.057646-0

Vigne, E., Marmonier, A., and Fuchs, M. (2008). Multiple interspecies recombination events within RNA2 of *Grapevine fanleaf* virus and *Arabis mosaic* virus. *Arch. Virol.* 153, 1771–1776. doi: 10.1007/s00705-008-0182-y

Walter, B., and Etienne, L. (1987). Detection of the *Grapevine fanleaf* viruses away from the period of vegetation. *J. Phytopathol.* 120, 355–364. doi: 10.1111/j.1439-0434.1987.tb00499.x

Walter, B., Vuittenez, A., Kuszala, J., Stocky, G., Burckard, J., and van Regenmortel, M. H. V. (1984). Détection sérologique des virus du court-noué de la vigne par le test ELISA. *Agronomie* 4, 527–534. doi: 10.1051/agro:19840605

Zhou, J., Fan, X., Dong, Y., Zhang, Z., Ren, F., Hu, G., et al. (2017). Complete nucleotide sequence of a new variant of grapevine fanleaf virus from northeastern China. *Arch. Virol.* 162, 577–579. doi: 10.1007/s00705-016-3101-7

**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.

Copyright © 2018 Vigne, Garcia, Komar, Lemaire and Hily. 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) and the copyright owner(s) 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.