# Evidence of differential spreading events of grapevine pinot Gris virus in Italy using datamining as a tool

Jean-Michel Hily, Véronique Komar, Nils Poulicard, Amandine Velt, Lauriane

Renault, Pierre Mustin, Emmanuelle Vigne, Anne-Sophie Spilmont, Olivier

Lemaire

► To cite this version:

Jean-Michel Hily, Véronique Komar, Nils Poulicard, Amandine Velt, Lauriane Renault, et al. Evidence of differential spreading events of grapevine pinot Gris virus in Italy using datamining as a tool. European Journal of Plant Pathology, 2021, 161 (3), pp.735-742. 10.1007/s10658-021-02343-3. hal-03341480

HAL Id: hal-03341480

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

Submitted on 10 Sep 2021

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_1.png" /></figure>

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

Check for
updates

# Evidence of differential spreading events of grapevine pinot Gris virus in Italy using datamining as a tool

Jean-Michel Hily ID • Véronique Komar • Nils Poulicard ID • Amandine Velt •

Lauriane Renault • Pierre Mustin • Emmanuelle Vigne ID • Anne-Sophie Spilmont •

Olivier Lemaire

Accepted: 24 July 2021

© The Author(s) 2021

**Abstract** Since its identification in 2003, grapevine Pinot gris virus (GPGV, *Trichovirus*) has now been detected in most grape-growing countries. So far, little is known about the epidemiology of this newly emerging virus. In this work, we used datamining as a tool to monitor *in-silico* the sanitary status of three vineyards in Italy. All data used in the study were recovered from a work that was already published and for which data were publicly available as SRA (Sequence Read Archive, NCBI) files. While incomplete, knowledge gathered from this work was still important, with evidence of differential accumulation of the virus in grapevine according to year, location, and variety-rootstock association. Additional data regarding GPGV genetic diversity were collected. Some advantages and pitfalls of datamining are discussed.

**Keywords** Grapevine · GPGV · Detection · Datamining

J.-M. Hily (✉) · A.-S. Spilmont
IFV, Le Grau-Du-Roi, France
e-mail: jean-michel.hily@vignevin.com

V. Komar · A. Velt · L. Renault · P. Mustin · E. Vigne ·
O. Lemaire
Université de Strasbourg, INRAE, SVQV UMR-A 1131,
F-68000 Colmar, France

N. Poulicard
PHIM Plant Health Institute, Univ Montpellier, IRD, CIRAD,
INRAE, Institut Agro, Montpellier, France

Since its characterization in Italy (Giampetruzzi et al., 2012), grapevine Pinot gris virus (GPGV, *Trichovirus*, *Betaflexiviridae*) has been detected in most grapevine growing regions around the world. Generally, the virus is detected using serological and/or molecular tools. In this work, we describe datamining as a potential additional method to identify grapevine infected with this virus, better estimating its distribution worldwide. While this specific work cannot be considered as an epidemiological study per se, it still unquestionably offers valuable information on the virus (i.e., its geographic distribution and genetic composition), providing a snapshot of the situation in three different vineyards in Italy at a specific time, giving new insight on GPGV accumulation, introduction and transmission.

This particular work is based on the data provided by a study on the contribution of genotype, the environment and their interaction to the berry transcriptome that was previously published (Dal Santo et al., 2018). Two cultivars, Cabernet Sauvignon and Sangiovese, were planted in three different locations: Montalcino, Bolgheri and Riccione. The former two Italian cities are located in the Tuscany hills and Tuscany coast respectively, while the latter is positioned on the Adriatic coast (Fig. 1). To minimize genetic variation, researchers used the same clonal material for each cultivar, with clones R5 and VCR23 of Cabernet Sauvignon and Sangiovese, respectively. In addition, three different rootstocks were tested in the study: Kober-5BB, 420A and 161.49 C. After uploading the 72 SRA files generated from this work, all samples were analyzed for the presence of GPGV using Workbench 12.0 software (CLC Genomics Workbench,
<figure><img src="image_3.png" /><figcaption>**Fig. 1** Maximum-likelihood tree inferred from sequences (7206 nt) of grapevine Pinot gris virus genome isolated from two cultivars, Cabernet Sauvignon clone R5 (star) and Sangiovese clone VCR23 (circle). Rootstocks are also indicated with 161.49 C (square), Kober 5BB (triangle) and 420A (diamond). Only bootstraps above 0.5 are shown. Colors correspond to the location</figcaption></figure>

in Italy where samples were recovered, Bolgheri (blue) and Riccione (red), see map on the upper right corner. Identity percentages between sequences are indicated on the right of the ML-tree. Measurements of population's differentiation (fixation index, F<sub>ST</sub>) and associated statistics (*P* value) are on the upper left corner

Aarhus, Denmark) as previously described (Hily et al., 2018). This was first assessed by mapping reads to a collection of curated GPGV reference sequences. For those displaying reads corresponding to GPGV, de novo assembly steps were performed and further extended by multiple rounds of residual reads mapping as previously described (Nourinejad Zarghani et al., 2018). Genome sequences being produced were ascertained using very stringent mapping parameters (length of 0.95/similarity of 0.95).

Our datamining study revealed that only samples from Bolgheri and Riccione were positive for GPGV. The virus was hardly detected in a few samples from Montalcino (Table 1); however, no complete sequence could be recovered. These ‘Low Read Count’ samples were probably the result of ‘intra-lane contamination’, as previously described in other studies (Vigne et al., 2018). When using RPKM (Reads per kilo base per million) data as a proxy for virus accumulation in the samples, our analyses revealed differential accumulation of GPGV according to many variables (Fig. 2). Indeed, GPGV seems to accumulate more in berries in 2011 than in 2012 (P < 10<sup>-5</sup>) and at a later stage of fruit development, at mid-ripening rather than pre-veraison (P < 10<sup>-4</sup>). Also, the association cultivar-rootstock seems to have its importance in virus

accumulation. Indeed, GPGV seems to accumulate more in Cabernet Sauvignon cultivar grafted onto either 161-49 or Kober-5BB rootstocks, rather than in Sangiovese grafted onto 420A at either location (P ≤ 10<sup>-4</sup>). In addition, differential accumulation of GPGV was also observed according to location where grapevines were grown (P < 10<sup>-5</sup>), with GPGV accumulating more in Riccione than in Bolgheri.

When delving into the genetics of the virus, other information was revealed. Overall, 47 complete genome GPGV sequences (or near complete, covering at least all open reading frames) were assembled (Table 1), all submitted to GenBank (BK011089-BK011101, and the other sequences are available upon request). After a phylogenetic analysis (Fig. 1), three major clades of GPGV were found to infect these grapevines, each displaying a high intra-clade nucleic acid identity percentage ≥ 98.10%. Interestingly, GPGV sequences seemed to cluster together very well by location (Fig. 1, colors) however independently from cultivar. Fixation index (F<sub>ST</sub>) analyses (Fig. 1) confirmed the genetic differentiation of the viral population according to location, showing a statistically significant high F<sub>ST</sub> value (F<sub>ST</sub> = 0.608, P ≤ 10<sup>-5</sup>). Such segregation by location was also highlighted for grapevines grafted
<table><caption>**Table 1** All information regarding the datamining analyses performed from the study from Dal Santo et al., 2018</caption><thead><tr><th>SEA #</th><th>hybridization #</th><th>Sample ID</th><th>Cultivar</th><th>Rootstock</th><th>Developmental Stage</th><th>Location</th><th>Vintage</th><th>GPGV</th><th>RPKM</th><th>Mapped read counts*</th><th>Total read counts</th><th>Genome length (nt)</th></tr></thead><tbody><tr><td>SRR5457593</td><td>4</td><td>CS_MO_PV_11_A</td><td>Cabernet Sauvignon</td><td>S04</td><td>Pre-veraison</td><td>Montalcino</td><td>2011</td><td></td><td></td><td></td><td>39,659,627</td><td></td></tr><tr><td>SRR5457594</td><td>5</td><td>CS_MO_PV_11_B</td><td>Cabernet Sauvignon</td><td>S04</td><td>Pre-veraison</td><td>Montalcino</td><td>2011</td><td></td><td></td><td></td><td>37,953,191</td><td></td></tr><tr><td>SRR5457595</td><td>6</td><td>CS_MO_PV_11_C</td><td>Cabernet Sauvignon</td><td>S04</td><td>Pre-veraison</td><td>Montalcino</td><td>2011</td><td></td><td></td><td></td><td>45,920,500</td><td></td></tr><tr><td>SRR5457596</td><td>7</td><td>CS_MO_MR_11_A</td><td>Cabernet Sauvignon</td><td>S04</td><td>Mid-ripening</td><td>Montalcino</td><td>2011</td><td></td><td></td><td></td><td>30,131,817</td><td></td></tr><tr><td>SRR5457597</td><td>8</td><td>CS_MO_MR_11_B</td><td>Cabernet Sauvignon</td><td>S04</td><td>Mid-ripening</td><td>Montalcino</td><td>2011</td><td></td><td></td><td></td><td>25,466,144</td><td></td></tr><tr><td>SRR5457598</td><td>9</td><td>CS_MO_MR_11_C</td><td>Cabernet Sauvignon</td><td>S04</td><td>Mid-ripening</td><td>Montalcino</td><td>2011</td><td></td><td></td><td></td><td>29,627,432</td><td></td></tr><tr><td>SRR5457599</td><td>16</td><td>SG_MO_PV_11_A</td><td>Sangiovese</td><td>420A</td><td>Pre-veraison</td><td>Montalcino</td><td>2011</td><td></td><td></td><td></td><td>30,253,594</td><td></td></tr><tr><td>SRR5457600</td><td>17</td><td>SG_MO_PV_11_B</td><td>Sangiovese</td><td>420A</td><td>Pre-veraison</td><td>Montalcino</td><td>2011</td><td></td><td></td><td></td><td>27,619,510</td><td></td></tr><tr><td>SRR5457601</td><td>18</td><td>SG_MO_PV_11_C</td><td>Sangiovese</td><td>420A</td><td>Pre-veraison</td><td>Montalcino</td><td>2011</td><td></td><td></td><td></td><td>24,825,638</td><td></td></tr><tr><td>SRR5457602</td><td>19</td><td>SG</td><td>Sangiovese</td><td>420A</td><td>Mid-ripening</td><td>Montalcino</td><td>2011</td><td></td><td></td><td></td><td>31,261,949</td><td></td></tr><tr><td>SRR5457603</td><td>20</td><td>SG</td><td>Sangiovese</td><td>420A</td><td>Mid-ripening</td><td>Montalcino</td><td>2011</td><td></td><td></td><td></td><td>37,850,541</td><td></td></tr><tr><td>SRR5457604</td><td>21</td><td>SG</td><td>Sangiovese</td><td>420A</td><td>Mid-ripening</td><td>Montalcino</td><td>2011</td><td></td><td></td><td></td><td>33,319,419</td><td></td></tr><tr><td>SRR5457605</td><td>28</td><td>CS_BO_PV_11_A</td><td>Cabernet Sauvignon</td><td>161-49</td><td>Pre-veraison</td><td>Polgheri</td><td>2011</td><td>2</td><td>91,43</td><td>19,934</td><td>30,211,399</td><td>7287, 7287</td></tr><tr><td>SRR5457606</td><td>29</td><td>CS_BO_PV_11_B</td><td>Cabernet Sauvignon</td><td>161-49</td><td>Pre-veraison</td><td>Polgheri</td><td>2011</td><td>2</td><td>44,32</td><td>10,113</td><td>31,519,652</td><td>7254, 7254</td></tr><tr><td>SRR5457607</td><td>30</td><td>CS_BO_PV_11_C</td><td>Cabernet Sauvignon</td><td>161-49</td><td>Pre-veraison</td><td>Polgheri</td><td>2011</td><td>1</td><td>100,96</td><td>25,001</td><td>34,310,824</td><td>7247</td></tr><tr><td>SRR5457605</td><td>31</td><td>CS_BO_MR_11_A</td><td>Cabernet Sauvignon</td><td>161-49</td><td>Mid-ripening</td><td>Polgheri</td><td>2011</td><td>1</td><td>173,45</td><td>42,993</td><td>34,345,114</td><td>7247</td></tr><tr><td>SRR5457609</td><td>32</td><td>CS_BO_MR_11_B</td><td>Cabernet Sauvignon</td><td>161-49</td><td>Mid-ripening</td><td>Polgheri</td><td>2011</td><td></td><td></td><td></td><td>32,004,939</td><td></td></tr><tr><td>SRR5457610</td><td>33</td><td>CS_BO_MR_11_C</td><td>Cabernet Sauvignon</td><td>161-49</td><td>Mid-ripening</td><td>Polgheri</td><td>2011</td><td></td><td></td><td></td><td>32,253,343</td><td></td></tr></tbody></table>
<table><caption>Table 1 (continued)</caption><thead><tr><th>SEA #</th><th>hybridization #</th><th>Sample ID</th><th>Cultivar</th><th>Rootstock</th><th>Developmental Stage</th><th>Location</th><th>Vintage</th><th>GPGV</th><th>RPKM</th><th>Mapped read counts*</th><th>Total read counts</th><th>Genome length (nt)</th></tr></thead><tbody><tr><td>SRR5457611</td><td>40</td><td>SG_BO_PV_11_A</td><td>Sangiovese</td><td>420A</td><td>Pre-veraison</td><td>Polgheri</td><td>2011</td><td>1</td><td>10.43</td><td>2425</td><td>32,216,454</td><td>7243</td></tr><tr><td>SRR5457612</td><td>41</td><td>SG_BO_PV_11_B</td><td>Sangiovese</td><td>420A</td><td>Pre-veraison</td><td>Polgheri</td><td>2011</td><td>1</td><td>9.64</td><td>2092</td><td>30,065,198</td><td>7240</td></tr><tr><td>SRR5457613</td><td>42</td><td>SG_BO_PV_11_C</td><td>Sangiovese</td><td>420A</td><td>Pre-veraison</td><td>Polgheri</td><td>2011</td><td>1</td><td>4.52</td><td>922</td><td>25,270,284</td><td>7213</td></tr><tr><td>SRR5457614</td><td>43</td><td>SG_BO_MR_11_A</td><td>Sangiovese</td><td>420A</td><td>Mid-ripening</td><td>Polgheri</td><td>2011</td><td>1</td><td>24.92</td><td>6360</td><td>35,361,602</td><td>7307</td></tr><tr><td>SRR5457615</td><td>44</td><td>SG_BO_MR_11_B</td><td>Sangiovese</td><td>420A</td><td>Mid-ripening</td><td>Polgheri</td><td>2011</td><td></td><td></td><td></td><td>30,185,292</td><td></td></tr><tr><td>SRR5457616</td><td>45</td><td>SG_BO_MR_11_C</td><td>Sangiovese</td><td>420A</td><td>Mid-ripening</td><td>Polgheri</td><td>2011</td><td>1</td><td>68.12</td><td>15.589</td><td>31,708,932</td><td>7290</td></tr><tr><td>SRR5457617</td><td>52</td><td>CS_RI_PV_11_A</td><td>Cabernet Sauvignon</td><td>Kober-5BB</td><td>Pre-veraison</td><td>Riccione</td><td>2011</td><td>1</td><td>128.03</td><td>28.440</td><td>30,778,512</td><td>7254</td></tr><tr><td>SRR5457618</td><td>53</td><td>CS_RI_PV_11_B</td><td>Cabernet Sauvignon</td><td>Kober-5BB</td><td>Pre-veraison</td><td>Riccione</td><td>2011</td><td>1</td><td>128.00</td><td>27.080</td><td>29,314,935</td><td>72.54</td></tr><tr><td>SRR5457619</td><td>54</td><td>CS_RI_PV_11_C</td><td>Cabernet Sauvignon</td><td>Kober-5BB</td><td>Pre-veraison</td><td>Riccione</td><td>2011</td><td>1</td><td>111.44</td><td>28.416</td><td>35,330,755</td><td>7254</td></tr><tr><td>SRR5457620</td><td>55</td><td>CS_RI_MR_11_A</td><td>Cabernet Sauvignon</td><td>Kober-5BB</td><td>Mid-ripening</td><td>Riccione</td><td>2011</td><td>1</td><td>2258.35</td><td>48.5449</td><td>29,784,834</td><td>7254</td></tr><tr><td>SRR5457621</td><td>56</td><td>CS_RI_MR_11_B</td><td>Cabernet Sauvignon</td><td>Kober-5BB</td><td>Mid-ripening</td><td>Riccione</td><td>2011</td><td>1</td><td>2225.20</td><td>455.935</td><td>28,390,737</td><td>72.54</td></tr><tr><td>SRR5457622</td><td>57</td><td>CS_RI_MR_11_C</td><td>Cabernet Sauvignon</td><td>Kober-5BB</td><td>Mid-ripening</td><td>Riccione</td><td>2011</td><td>1</td><td>1565.76</td><td>284.493</td><td>25,176,180</td><td>72.54</td></tr><tr><td>SRR5457623</td><td>64</td><td>SG_RI_PV_11_A</td><td>Sangiovese</td><td>420A</td><td>Pre-veraison</td><td>Riccione</td><td>2011</td><td>1</td><td>89.48</td><td>17.871</td><td>27,673,291</td><td>7258</td></tr><tr><td>SRR5457624</td><td>65</td><td>SG_RI_PV_11_B</td><td>Sangiovese</td><td>420A</td><td>Pre-veraison</td><td>Riccione</td><td>2011</td><td>1</td><td>58.23</td><td>11.621</td><td>27,651,896</td><td>7254</td></tr><tr><td>SRR5457625</td><td>66</td><td>SG_RI_PV_11_C</td><td>Sangiovese</td><td>420A</td><td>Pre-veraison</td><td>Riccione</td><td>2011</td><td>1</td><td>135.88</td><td>27.403</td><td>27,943,588</td><td>7254</td></tr><tr><td>SRR5457626</td><td>67</td><td>SG_RI_MR_11_A</td><td>Sangiovese</td><td>420A</td><td>Mid-ripening</td><td>Riccione</td><td>2011</td><td>1</td><td>299.59</td><td>48.006</td><td>22,202,853</td><td>7289</td></tr><tr><td>SRR5457627</td><td>68</td><td>SG_RI_MR_11_B</td><td>Sangiovese</td><td>420A</td><td>Mid-ripening</td><td>Riccione</td><td>2011</td><td>1</td><td>414.33</td><td>89.289</td><td>29,860,068</td><td>7254</td></tr><tr><td>SRR5457628</td><td>59</td><td>SG_RI_MR_11_C</td><td>Sangiovese</td><td>420A</td><td>Mid-ripening</td><td>Riccione</td><td>2011</td><td>1</td><td>300.74</td><td>61.377</td><td>28,278,938</td><td>7254</td></tr><tr><td>SRR5457629</td><td>91</td><td>SG_BO_PV_12_A</td><td>Sangiovese</td><td>420A</td><td>Pre-veraison</td><td>Bolgheri</td><td>2012</td><td>1</td><td>3.99</td><td>884</td><td>30,685,737</td><td>7214</td></tr><tr><td>SRR5457630</td><td>92</td><td>SG_BO_PV_12_B</td><td>Sangiovese</td><td>420A</td><td>Pre-veraison</td><td>Bolgheri</td><td>2012</td><td>1</td><td>3.29</td><td>684</td><td>28,765,541</td><td>7250</td></tr><tr><td>SRR5457631</td><td>93</td><td>SG_BO_PV_12_C</td><td>Sangiovese</td><td>420A</td><td>Pre-veraison</td><td>Bolgheri</td><td>2012</td><td>✓</td><td>1.87</td><td>455</td><td>33,797,617</td><td>7131</td></tr><tr><td>SRR5457632</td><td>94</td><td>SG_MO_PV_12_A</td><td>Sangiovese</td><td>420A</td><td>Pre-veraison</td><td>Montalcino</td><td>2012</td><td>✓</td><td>0.69</td><td>143</td><td>28,565,019</td><td>4800</td></tr><tr><td>SRR5457633</td><td>95</td><td>SG_MO_PV_12_B</td><td>Sangiovese</td><td>420A</td><td>Pre-veraison</td><td>Montalcino</td><td>2012</td><td>✓</td><td>2.99</td><td>675</td><td>31,322,839</td><td>7201</td></tr></tbody></table>

<figure><img src="image_4.png" /></figure>
<table><caption>Table 1 (continued)</caption><thead><tr><th>SEA #</th><th>hybridization #</th><th>Sample ID</th><th>Cultivar</th><th>Rootstock</th><th>Developmental Stage</th><th>Location</th><th>Vintage</th><th>GPGV</th><th>RPKM</th><th>Mapped read counts*</th><th>Total read counts</th><th>Genome length (nt)</th></tr></thead><tbody><tr><td>SRR5457634</td><td>96</td><td>SG_MO_PV_12_C</td><td>Sangiovese</td><td>420A</td><td>Pre-veraison</td><td>Montalcino</td><td>2012</td><td>✓</td><td>1.43</td><td>312</td><td>30,193,456</td><td>6686</td></tr><tr><td>SRR5457635</td><td>97</td><td>SG_RI_PV_12_A</td><td>Sangiovese</td><td>420A</td><td>Pre-veraison</td><td>Riccione</td><td>2012</td><td>1</td><td>46.29</td><td>10,705</td><td>32,044,752</td><td>7257</td></tr><tr><td>SRR5457636</td><td>98</td><td>SG_RI_PV_12_B</td><td>Sangiovese</td><td>420A</td><td>Pre-veraison</td><td>Riccione</td><td>2012</td><td>1</td><td>41.83</td><td>7769</td><td>25,735,588</td><td>7253</td></tr><tr><td>SRR5457637</td><td>99</td><td>SG_RI_PV_12_C</td><td>Sangiovese</td><td>420A</td><td>Pre-veraison</td><td>Riccione</td><td>2012</td><td>1</td><td>38.53</td><td>8245</td><td>29,653,480</td><td>7271</td></tr><tr><td>SRR5457539</td><td>100</td><td>CS_MO_PV_12_A</td><td>Sangiovese</td><td>420A</td><td>Pre-veraison</td><td>Montalcino</td><td>2012</td><td></td><td></td><td></td><td>28,374,413</td><td></td></tr><tr><td>SRR5457639</td><td>101</td><td>CS_MO_PV_12_B</td><td>Cabernet Sauvignon</td><td>S04</td><td>Pre-veraison</td><td>Montalcino</td><td>2012</td><td></td><td></td><td></td><td>39,038,471</td><td></td></tr><tr><td>SRR5457640</td><td>102</td><td>CS_MO_PV_12_C</td><td>Cabernet Sauvignon</td><td>S04</td><td>Pre-veraison</td><td>Montalcino</td><td>2012</td><td></td><td></td><td></td><td>29,599,165</td><td></td></tr><tr><td>SRR5457641</td><td>103</td><td>CS_RI_PV_12_A</td><td>Cabernet Sauvignon</td><td>Kober-5BB</td><td>Pre-veraison</td><td>Riccione</td><td>2012</td><td>1</td><td>55.19</td><td>10,488</td><td>26,329,353</td><td>7253</td></tr><tr><td>SRR5457642</td><td>104</td><td>CS_RI_PV_12_B</td><td>Cabernet Sauvignon</td><td>Kober-5BB</td><td>Pre-veraison</td><td>Riccione</td><td>2012</td><td>1</td><td>50.26</td><td>11,045</td><td>30,452,556</td><td>7253</td></tr><tr><td>SRR5457643</td><td>105</td><td>CS_RI_PV_12_C</td><td>Cabernet Sauvignon</td><td>Kober-5BB</td><td>Pre-veraison</td><td>Riccione</td><td>2012</td><td>1</td><td>63.52</td><td>14,937</td><td>32,582,117</td><td>7253</td></tr><tr><td>SRR5457644</td><td>106</td><td>CS_BO_PV_12_A</td><td>Cabernet Sauvignon</td><td>161-49</td><td>Pre-veraison</td><td>Bolgheri</td><td>2012</td><td>1</td><td>20.46</td><td>2295</td><td>15,541,092</td><td>7277</td></tr><tr><td>SRR5457645</td><td>107</td><td>CS_BO_PV_12_B</td><td>Cabernet Sauvignon</td><td>161-49</td><td>Pre-veraison</td><td>Bolgheri</td><td>2012</td><td>1</td><td>14.68</td><td>2996</td><td>28,275,962</td><td>7223</td></tr><tr><td>SRR5457646</td><td>108</td><td>CS_BO_PV_12_C</td><td>Cabernet Sauvignon</td><td>161-49</td><td>Pre-veraison</td><td>Bolgheri</td><td>2012</td><td>1</td><td>26.84</td><td>12,934</td><td>66,769,968</td><td>7282</td></tr><tr><td>SRR5457647</td><td>109</td><td>SG_BO_MR_12_A</td><td>Sangiovese</td><td>420A</td><td>Mid-ripening</td><td>Bolgheri</td><td>2012</td><td>✓</td><td>0.97</td><td>216</td><td>30,804,911</td><td>6126</td></tr><tr><td>SRR5457648</td><td>110</td><td>SG_BO_MR_12_B</td><td>Sangiovese</td><td>420A</td><td>Mid-ripening</td><td>Bolgheri</td><td>2012</td><td>1</td><td>6.31</td><td>1463</td><td>32,129,314</td><td>7219</td></tr><tr><td>SRR5457649</td><td>111</td><td>SG_BO_MR_12_C</td><td>Sangiovese</td><td>420A</td><td>Mid-ripening</td><td>Bolgheri</td><td>2012</td><td>✓</td><td>2.15</td><td>388</td><td>25,018,444</td><td>6948</td></tr><tr><td>SRR5457650</td><td>112</td><td>SG_MO_MR_12_A</td><td>Sangiovese</td><td>420A</td><td>Mid-ripening</td><td>Montalcino</td><td>2012</td><td></td><td></td><td></td><td>24,003,382</td><td></td></tr><tr><td>SRR5457651</td><td>113</td><td>SG_MO_MR_12_B</td><td>Sangiovese</td><td>420A</td><td>Mid-ripening</td><td>Montalcino</td><td>2012</td><td></td><td></td><td></td><td>37,168,759</td><td></td></tr><tr><td>SRR5457652</td><td>114</td><td>SG_MO_MR_12_C</td><td>Sangiovese</td><td>420A</td><td>Mid-ripening</td><td>Montalcino</td><td>2012</td><td></td><td></td><td></td><td>29,938,586</td><td></td></tr><tr><td>SRR5457653</td><td>115</td><td>SG_RI_MR_12_A</td><td>Sangiovese</td><td>420A</td><td>Mid-ripening</td><td>Riccione</td><td>2012</td><td>1</td><td>28.48</td><td>7041</td><td>34,255,543</td><td>7250</td></tr><tr><td>SRR5457654</td><td>116</td><td>SG_RI_MR_12_B</td><td>Sangiovese</td><td>420A</td><td>Mid-ripening</td><td>Riccione</td><td>2012</td><td>1</td><td>31.09</td><td>6790</td><td>30,258,155</td><td>7250</td></tr><tr><td>SRR5457655</td><td>117</td><td>SG_RI_MR_12_C</td><td>Sangiovese</td><td>420A</td><td>Mid-ripening</td><td>Riccione</td><td>2012</td><td>1</td><td>7.48</td><td>1524</td><td>28,230,567</td><td>7255</td></tr></tbody></table>
<table><caption>**Table 1 (continued)**</caption><thead><tr><th>SEA #</th><th>hybridization #</th><th>Sample ID</th><th>Cultivar</th><th>Rootstock Stage</th><th>Developmental Stage</th><th>Location</th><th>Vintage</th><th>GPGV</th><th>RPKM</th><th>Mapped read counts*</th><th>Total read counts</th><th>Genome length (nt)</th></tr></thead><tbody><tr><td>SRR5457656</td><td>118</td><td>CS_MO_MR_12_A</td><td>Sangiovese</td><td>420A</td><td>Mid-ripening</td><td>Montalcino</td><td>2012</td><td>✓</td><td>1.05</td><td>224</td><td>29,549,033</td><td>6431</td></tr><tr><td>SRR5457657</td><td>119</td><td>CS_MO_MR_12_B</td><td>Cabernet Sauvignon</td><td>S04</td><td>Mid-ripening</td><td>Montalcino</td><td>2012</td><td></td><td></td><td></td><td>22,749,636</td><td></td></tr><tr><td>SRR5457658</td><td>120</td><td>CS_MO_MR_12_C</td><td>Cabernet Sauvignon</td><td>S04</td><td>Mid-ripening</td><td>Montalcino</td><td>2012</td><td></td><td></td><td></td><td>29,723,920</td><td></td></tr><tr><td>SRR5457659</td><td>121</td><td>CS_RI_MR_12_A</td><td>Cabernet Sauvignon</td><td>Kober-5BB</td><td>Mid-ripening</td><td>Riccione</td><td>2012</td><td>1</td><td>408,02</td><td>82,399</td><td>27,982,523</td><td>7250</td></tr><tr><td>SRR5457660</td><td>122</td><td>CS_RI_MR_12_B</td><td>Cabernet Sauvignon</td><td>Kober-5BB</td><td>Mid-ripening</td><td>Riccione</td><td>2012</td><td>1</td><td>923,80</td><td>235,636</td><td>35,343,145</td><td>7284</td></tr><tr><td>SRR5457661</td><td>123</td><td>CS_RI_MR_12_C</td><td>Cabernet Sauvignon</td><td>Kober-5BB</td><td>Mid-ripening</td><td>Riccione</td><td>2012</td><td>1</td><td>1336,91</td><td>318,331</td><td>32,992,910</td><td>7277</td></tr><tr><td>SRR5457662</td><td>124</td><td>CS_BO_MR_12_A</td><td>Cabernet Sauvignon</td><td>161-49</td><td>Mid-ripening</td><td>Bolgheri</td><td>2012</td><td>2</td><td>117,96</td><td>24,447</td><td>28,741,196</td><td>7217,7217</td></tr><tr><td>SRR5457663</td><td>125</td><td>CS_130_MR_12_B</td><td>Cabernet Sauvignon</td><td>161-49</td><td>Mid-ripening</td><td>Bolgheri</td><td>2012</td><td>2</td><td>74,49</td><td>209,41</td><td>38,951,034</td><td>7217,7217</td></tr><tr><td>SRR5457664</td><td>126</td><td>CS_130_MR_12_C</td><td>Cabernet Sauvignon</td><td>161-49</td><td>Mid-ripening</td><td>Bolgheri</td><td>2012</td><td>2</td><td>63,53</td><td>15,221</td><td>33,195,577</td><td>7217,7217</td></tr></tbody></table>

The 'number' in the GPGV column correspond to the number of complete genome assembled in de novo in each sample. ✓ indicates that reads have mapped onto GPGV genome, as shown in the Mapped read counts columns would indicate, however no complete genome from contiguous sequence could be obtained and RPKM (Read per Kilobase Million) were always below 3 when no genome were assembled. This work was performed using CLC-Workbench using very stringent mapping parameters * (0,95/0,95)

<figure><img src="image_5.png" /></figure>
<figure><img src="image_6.png" /><figcaption>**Fig. 2** Box plot diagrams of RPKM in function of different variables. From left to right: year, developmental stage (MR: mid-ripening, PV: pre-veraison), rootstock, overall location, Sangiovese grafted onto 420A, grapevine cultivated in Bolgheri and in Riccione (CS: Cabernet Sauvignon). On each box, the central line is the median, the edges of the boxes are the 25th and 75th percentiles, the whiskers extend to the most extreme data and the</figcaption></figure>

onto rootstock 161.49 C used exclusively in Bolgheri and grapevines onto Kober 5BB exclusively used in Riccione ($F_{ST} = 0.564, P \le 10^{-5}$). Comparison of sequences obtained from the 420A rootstock also displayed statistically significant $F_{ST}$ values; however, the values were lower than the ones mentioned above. This is most likely because 420A was used as a rootstock in both locations. Furthermore, the genetic background of the grapevine cultivar, which was also present in both locations, had no statistically significant impact on viral populations ($F_{ST} = 0.018, P = 0.185$).

In addition to the presence/absence of GPGV in the samples, this work highlights two distinct situations at the viral genomic level. Indeed, one vineyard is infected by a single variant, identity percentage ≥ 99%, as previously defined for GPGV (Hily et al., 2020), represented here by samples from the Bolgheri region, while the other vineyard (Riccione) is infected by at least two (or more) variants. These results indirectly, but strongly, suggest probable independent introduction/transmission events of GPGV in two out of the three locations specifically looked at, in this transcriptomic study. These situations are probably the result of transmission events through grafting (Saldarelli et al., 2014) and movement of infected material as previously suggested (Al Rwahnih et al., 2016; Fajardo et al., 2017; Wu & Habili, 2017). They may also have occurred horizontally by vectors either in the nursery or in the vineyard,

dots refer to the outliers. Since RPKM values did not follow a normal distribution, a generalized linear model (GLM) with Poisson link function was used. The significance of the considered effect was tested using Wald chi2 test and the *p* values smaller than 0.05 threshold were considered statistically significant. All analyses and graphic representations were made with the R software version 4.0.2 (R core Team 2012).

with distinct variants of the virus being detected at each location, regardless of the clonal background of the grapevine. In addition, the detection of these different variants according to location, each displaying probable differences in fitness, may results in differential virus accumulation as observed above. Overall, this in silico work add onto the so-far limited knowledge on the natural transmission of GPGV in vineyards (Bertazzon et al., 2020; Hily et al., in press).

Lately, datamining is becoming a very important and powerful tool to identify new pathogens, as well as new variants of known viruses, such as from the now well-known *Coronaviridae* family for example (https://virological.org/t/serratus-the-ultra-deep-search-to-discover-novel-coronaviruses/516) (last visited 04/2021). Datamining can be also utilized to increase the number of complete genome sequences for downstream studies on the evolutionary history of specific viruses for example (Hily et al., 2020). In this work, datamining can be considered as an *in-silico* tool to monitor *post facto* the sanitary status of any vineyards around the world from which data have already been collected, published and made publicly available. There are a few pitfalls regarding datamining as a tool. Indeed, we do not have always all the details regarding the samples (i.e. metadata about the samples such as the exact origin and location of collection). We do not have the choice of the technology with which data were obtained nor the quality of the sample. However, the
information being generated is still very valuable, it has already been paid for and therefore almost free (other than the time of analysis), it is available to anyone and most of all, it is ever growing.

**Funding** This work was supported by Institut National de Recherche pour l'Agriculture, l'Alimentation et l'Environnement (INRAE) and by Institut Français de la Vigne et du Vin (IFV) and by the projects 'VACCIVINE' and 'GPGV' funded through 'Plan National Dépérissément du vignoble' (French Ministry of Agriculture) in 2017 and 2019, respectively. A fellowship from Moët & Chandon, Comité Interprofessionnel du vin de Champagne (CIVC), Bureau Interprofessionnel des Vins de Bourgogne (BIVB) and Comité Interprofessionnel des Vins d'Alsace (CIVA) was awarded to JMH.

## Declarations

**Ethical approval** No human and/or animal participants were involved in the study.

**Consent to participate, submission and for publication** All authors have been personally and actively involved in the work leading to this manuscript and consent its submission and publication.

**Conflict of interest** The authors declare that they have no conflict of interest.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

## References

Al Rwahnih, M., Golino, D., & Rowhani, A. (2016). First report of grapevine pinot gris virus infecting grapevine in the United States. *Plant Disease, 100*(5), 1030–1030. https://doi.org/10.1094/pdis-10-15-1235-pdn

Bertazzon, N., Forte, V., & Angelini, E. (2020). Fast transmission of grapevine 'Pinot gris' virus (GPGV) in vineyard. *Vitis, 59*(1), 29–34.

Dal Santo, S., Zenoni, S., Sandri, M., De Lorenzis, G., Magris, G., De Paoli, E., et al. (2018). Grapevine field experiments reveal the contribution of genotype, the influence of environment and the effect of their interaction (G×E) on the berry transcriptome. *The Plant Journal, 93*(6), 1143–1159. https://doi.org/10.1111/tpj.13834

Fajardo, T. V. M., Eiras, M., & Nickel, O. (2017). First report of grapevine pinot gris virus infecting grapevine in Brazil. *Australasian Plant Disease Notes, 12*(1), 45. https://doi.org/10.1007/s13314-017-0270-5

Giampetruzzi, A., Roumi, V., Roberto, R., Malossini, U., Yoshikawa, N., La Notte, P., et al. (2012). A new grapevine virus discovered by deep sequencing of virus- and viroid-derived small RNAs in cv Pinot gris. *Virus Research, 163*(1), 262–268. https://doi.org/10.1016/j.virusres.2011.10.010

Hily, J.-M., Demanèche, S., Poulicard, N., Tannières, M., Djennane, S., Beuve, M., Vigne, E., Demangeat, G., Komar, V., Gertz, C., Marmonier, A., Hemmer, C., Vigneron, S., Marais, A., Candresse, T., Simonet, P., & Lemaire, O. (2018). Metagenomic-based impact study of transgenic grapevine rootstock on its associated virome and soil bacteriome. *Plant Biotechnology Journal, 16*(1), 208–220. https://doi.org/10.1111/pbi.12761

Hily, J.-M., Poulicard, N., Candresse, T., Vigne, E., Beuve, M., Renault, L., Velt, A., Spilmont, A. S., & Lemaire, O. (2020). Datamining, genetic diversity analyses and phylogeographic reconstructions redefine the worldwide evolutionary history of grapevine Pinot gris virus and grapevine berry inner necrosis virus. *Phytobiomes Journal, 4*(2), 165–177. https://doi.org/10.1094/PBIOMES-10-19-0061-R

Hily, J.-M., Komar, V., Poulicard, N., Vigne, E., Jacquet, O., Protet, N., Spilmont, A. S., & Lemaire, O. (in press). Biological evidence and molecular modeling of a grapevine Pinot gris virus outbreak in a vineyard. *Phytobiomes Journal*. https://doi.org/10.1094/pbiomes-11-20-0079-r

Nourinejad Zarghani, S., Hily, J. M., Glasa, M., Marais, A., Wetzel, T., Faure, C., Vigne, E., Velt, A., Lemaire, O., Boursiquot, J. M., Okic, A., Ruiz-Garcia, A. B., Olmos, A., Lacombe, T., & Candresse, T. (2018). Grapevine virus T diversity as revealed by full-length genome sequences assembled from high-throughput sequence data. *PLoS One, 13*(10), e0206010. https://doi.org/10.1371/journal.pone.0206010

Saldarelli, P., Giampetruzzi, A., Morelli, M., Malossini, U., Pirolo, C., Bianchedi, P., & Gualandri, V. (2014). Genetic variability of grapevine Pinot gris virus and its association with grapevine leaf mottling and deformation. *Phytopathology, 105*(4), 555–563. https://doi.org/10.1094/PHYTO-09-14-0241-R

Vigne, E., Garcia, S., Komar, V., Lemaire, O., & 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. *Frontiers in Microbiology, 22*(9), 2726. https://doi.org/10.3389/fmicb.2018.02726

Wu, Q., & Habili, N. (2017). The recent importation of grapevine Pinot gris virus into Australia. *Virus Genes, 53*(6), 935–938. https://doi.org/10.1007/s11262-017-1475-6