Vanja Miljanić1
Jernej Jakše
Urban Kunej
Denis Rusjan
Andreja Škvarč
Nataša Štajner- 1Department of Agronomy, Biotechnical Faculty, University of Ljubljana, Ljubljana, Slovenia
- 2Chamber of Agriculture and Forestry of Slovenia, Agriculture and Forestry Institute Nova Gorica, Nova Gorica, Slovenia
Diseases caused by viruses and virus-like organisms are one of the major problems in viticulture and grapevine marketing worldwide. Therefore, rapid and accurate diagnosis and identification is crucial. In this study, we used HTS of virus- and viroid-derived small RNAs to determine the virome status of Slovenian preclonal candidates of autochthonous and local grapevine varieties (Vitis vinifera L.). The method applied to the studied vines revealed the presence of nine viruses and two viroids. All viral entities were validated and more than 160 Sanger sequences were generated and deposited in NCBI. In addition, a complete description into the co-infections in each plant studied was obtained. No vine was found to be virus- and viroid-free, and no vine was found to be infected with only one virus or viroid, while the highest number of viral entities in a plant was eight.
Introduction
Grapevine is one of the most important fruit crops by acreage and economic importance (Torregrosa et al., 2015). According to the International Organization of Vine and Wine Intergovernmental Organization (OIV), 7.4 million hectares around the world were planted with grapevines in 2018, and vineyards in Spain, China, France, Italy, and Turkey represented 50% of the total world cultivated grapevine area. According to the OIV, the world production of grapes in 2018 was 77.8 million tons. In Slovenia in 2018, grapevines occupied an area of 15,630 hectares and annual production of grapes was 126,958 tons.1
Grapevine may harbor more than 86 viruses and viroids, belonging to different families and genera (Fuchs, 2020). Viruses and virus-like organisms cause severe damage to grapevine production worldwide. They cause leaf degeneration, malformation, puckering, leaf rolling, chlorosis, necrosis, ringspots, line patterns, mosaic patterns, vein-banding, vein-clearing, stunting, wilting, shortened internodes, fasciation, zigzag growth, grooving, cracking, and pitting of wood (Credi et al., 1981; Mavrič et al., 2003; Andret-Link et al., 2004; Bouyahia et al., 2005; Al Rwahnih et al., 2009; Lunden et al., 2010; Giampetruzzi et al., 2012; Mavrič Pleško et al., 2014; Sudarshana et al., 2015). They impact vine yield and wine quality, as viral entities delay ripening, affect grape quality, decrease sugar content, affect the content of pigments, various aromatic components and other metabolites, and increase the acidity of wines (Lee and Martin, 2009; Vega et al., 2011; Alabi et al., 2016; Girardello et al., 2020; Lee et al., 2021). Viral entities eventually lead to the death of chronically infected plants.
Therefore, rapid and accurate diagnosis and identification is very important. Most methods for detection and identification require prior knowledge of the potential pathogens (e.g., use of antibodies in serological methods or virus specific primers in PCR amplification), with the exception of the metagenomic approach called high-throughput sequencing technology (HTS). HTS is a powerful technology that enables rapid detection of viral entities in plant tissues, including unknown as well as known viruses and viroids in symptomatic and asymptomatic plants, without the need for prior knowledge (Al Rwahnih et al., 2009; Kreuze et al., 2009; Fajardo et al., 2020). HTS of small RNAs (small RNA sequencing; sRNA-seq) has been shown to be efficient in detecting plant viruses or viroids (Kreuze et al., 2009; Kashif et al., 2012; Vives et al., 2013; Jakse et al., 2015; Singh et al., 2020). This approach exploits a natural antiviral defense mechanism called RNA silencing or RNA interference (RNAi). The silencing mechanism is initiated by RNase III-like enzymes called Dicer-like enzymes (DCL) which cleave long double-stranded RNAs (dsRNAs) into short interfering (si)RNA and miRNA precursors with a hairpin or stem-loop structure into miRNA duplexes (miRNA/miRNA*) (Bernstein et al., 2001; Bartel, 2004; Baulcombe, 2004). During the process of viral infection small RNAs (sRNAs) accumulate abundantly in plants and can be detected by deep sequencing of infected plants. sRNA-seq provides a unique opportunity to easily detect and identify grapevine viruses and viroids due to the abundance of sRNAs (Navarro et al., 2009; Giampetruzzi et al., 2012; Eichmeier et al., 2016; Czotter et al., 2018; Demian et al., 2020; Li et al., 2021).
Slovenia is a traditional wine-growing country with many local and indigenous grapevine varieties revitalized in current clonal selection programs, according to the rules on the marketing of material for the vegetative propagation of vine (Official gazette, N°93/05 and 101/20) and OIV process for the clonal selection of vines (Resolution oiv-viti-564a-2017). In the past, propagation material was controlled just visually, which led to uncontrolled spread of viruses.
The aim of the presented work was to investigate the virome status of Slovenian preclonal grapevine candidates and to study their genetic diversity and co-infections using identification by sRNA-seq and confirmation by RT-PCR and Sanger sequencing.
Materials and Methods
Plant Material
A total of 82 dormant cuttings of 6 preclonal grapevine varieties (Vitis vinifera L.)—2 reds, “Refošk” (“Terrano”) and “Pokalca” (“Schioppettino”), and 4 whites, “Laški rizling” (“Welschriesling”), “Rebula” (“Ribolla Gialla”), “Malvazija” (“Malvasia d’Istria”), and “Zeleni Sauvignon” (“Sauvignon Vert”), were collected from the 3 vineyards [Pouzelce (P); Base (B) and Genebank (G), referenced to Table 1)] maintained by Centre of grapevine selection (STS Vrhpolje, Vipava) in Primorska wine-growing region in Slovenia in February 2019. After 3–4 weeks in water at room temperature one-bud cuttings started bud-bursting and the obtained leaves were collected and stored at –80°C for further analysis.
Table 1. Viruses and viroids detected in 12 libraries using the VirusDetect approach.
High-Throughput Sequencing of Virus- and Viroid-Derived Small RNAs
Eighty-two grapevine samples were pooled into 12 bulks/libraries, each bulk representing samples of the same variety; some varieties were represented with more than one bulk (Table 1). sRNAs were extracted by an enrichment procedure using a mirVana miRNA Isolation Kit (Ambion, Life Technologies) according to the manufacturer’s instructions. sRNA libraries were prepared using the Ion Total RNA-Seq kit and checked for quality using an Agilent 2100 Bioanalyzer (Agilent Technologies). Barcode-labeled cDNA libraries were sequenced on an Ion PI chip v3 using an Ion Proton Sequencer (Ion Torrent; Life Technologies) according to the manufacturer’s instructions. According to the Ion Torrent sequencing pipeline, raw reads had removed adapter sequences and they were deposited in the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) database under BioProject accession number PRJNA765925.
Analysis of High-Throughput Sequencing Data
Raw sequence reads were filtered based on quality score and read length using cutadapt tools (Martin, 2011). VirusDetect2 (Zheng et al., 2017), an automated bioinformatics pipeline, was employed for further analysis of obtained sequences. VirusDetect is a highly sensitive and enables efficient analysis of sRNA datasets for viral identification. The software package first align sRNA reads to viral GenBank references using Burrows-Wheeler Aligner (BWA). The mapped sRNA reads are then assembled into contigs according to the viral reference. This program also maps the sRNA reads to host reference sequences to discard host-derived sRNAs. VirusDetect performs de novo assembly of sRNAs using Velvet with automated parameter optimization. Contigs obtained from de novo assemblies were aligned to the grapevine genome and all contigs with nucleotide identity greater than 90% with the grapevine genome were discarded. De novo assembled contigs were concatenated with reference-guided generated contigs and then all redundancies were removed, according to the employed iAssembler pipeline. The obtained contigs were then compared with the viral GenBank references for their identification. This automated pipeline used the BLASTN algorithm to compare the contigs to the reference virus nucleotide sequences and BLASTX algorithm to compare the contigs to the reference virus protein sequences (Zheng et al., 2017).
Validation of Predicted Infections by RT-PCR, Direct Sanger Sequencing and Cloning, Sequence Analysis and Phylogenetic Studies
Validation of HTS results of predicted infections was performed by RT-PCR and Sanger sequencing. Total RNAs for all individual samples were extracted from 70 to 100 mg of frozen leaves using Monarch RNA Total Miniprep Kit (New England Biolabs) following recommended instructions. The RNAs concentration and purity were assessed with NanoVue Plus spectrophotometer (GE Healthcare Life Sciences). Due to the low RNA concentration and purity, three samples (Malvazija 20/2P, Malvazija 21/7P and Malvazija 23/4P) were excluded from further analysis. RT−PCRs were performed using a two−step protocol, where total RNAs were first reverse transcribed using High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems) according to the manufacturer’s instructions followed by PCR with specific primers (Supplementary Table 2). The PCR reaction mixture (20 μL total) contained 10.7 μL nuclease-free water, 4 μL 5 × PCR buffer (Promega), 1.6 μL MgCl2 (Kapa Biosystems), 1.6 μL dNTP mix (10 mM each of the 4 dNTPs) (Promega), 0.5 μL of each primer, 0.1 μL KAPA Taq DNA polymerase (Kapa Biosystems), and 1 μL of cDNA. RT-PCR products were analyzed by electrophoresis on a 1.4% agarose gel, stained with ethidium bromide, and visualized under a UV transilluminator and remaining reaction was cleaned by Exo-Sap treatment. RT-PCR products of a predicted sizes were sequenced directly in both directions for all viruses and viroids, except for Grapevine rupestris stem pitting-associated virus (GRSPaV), where RT-PCR products were ligated into the pGEM-Teasy Vector Systems cloning kit (Promega) and then transformed into Escherichia coli DH-5α competent cells. Blue/white screening was performed on the LB/carbenicillin/IPTG/X-gal/agar plates. The positive clones were randomly picked and then colony PCR was performed using specific primers (RSP 52/RSP 53) (Supplementary Table 2). After purification, direct RT-PCR or cloned products were sequenced using an Applied Biosystems 3130 Genetic Analyzer. After sequencing, the forward and reverse traces were trimmed and assembled using CodonCode Aligner 9.0.1 (CodonCode Corporation). All sequences generated in this study were deposited in the NCBI GenBank database.3 The generated virus and viroid sequences were compared using the ClustalW program (Thompson et al., 1994) implemented in MEGA X software (Kumar et al., 2018). A p-distance model was applied for nucleotide (nt) and deduced amino acid (aa) divergence sequence analysis. Phylogenetic trees were constructed using MEGA X software. The Modeltest implemented in MEGA X was applied to investigate the best-fitting model of nt substitution. The reliability of the obtained trees was evaluated using the bootstrap method based on 1,000 replicates and bootstrap values lower than 50% were omitted.
Results and Discussion
Viruses and Viroids Detected by High-Throughput Sequencing of Virus- and Viroid-Derived Small RNAs
sRNA-seq from the pooled grapevine samples yielded 4,206,135–17,668,261 reads. In VirusDetect pipeline 3,643,531–13,905,492 reads per pool were processed (Supplementary Table 1). Additional results of the detection pipeline are presented in Supplementary Table 1. Using the described approach, nine viruses: Raspberry bushy dwarf virus (RBDV), Grapevine Pinot gris virus (GPGV), Grapevine rupestris stem pitting-associated virus (GRSPaV), Grapevine fanleaf virus (GFLV), Grapevine leafroll-associated virus 3 (GLRaV-3), Grapevine fleck virus (GFkV), Grapevine Red Globe virus (GRGV), Grapevine rupestris vein feathering virus (GRVFV), Grapevine Syrah virus-1 (GSyV-1), and two viroids: Hop stunt viroid (HSVd) and Grapevine yellow speckle viroid- 1 (GYSVd-1) were identified. GRGV, GRVFV, and GSyV-1 were detected for the first time in Slovenia (paper in review).
The highest number of viral entities in a library was eight (in libraries 005 and 010), the lowest number was five (in libraries 007, 009, 014), while the remaining libraries contained six viral entities. The genomes of RBDV and GFLV are bipartite and consist of two single-stranded positive-sense RNAs (RNA1 and RNA2), therefore the consensus length, reference coverage, number of contigs and sequencing depth for both, RNA1 and RNA2 are shown (Table 1). GPGV, HSVd, and GYSVd-1 were detected in all analyzed libraries. RBDV and GLRaV-3 were detected in only one library, 005 (“Laški rizling” variety) and 010 (“Refošk” variety), respectively (Table 1).
Validation of Predicted Infections by RT-PCR, Direct Sanger Sequencing and Cloning, Sequence Analysis and Phylogenetic Studies
Raspberry Bushy Dwarf Virus
The natural occurrence of RBDV in grapevine was first confirmed in Slovenia in “Laški rizling” and “Štajerska belina” using DAS-ELISA and IC-RT-PCR (Mavrič et al., 2003). Reports that this virus infecting grapevines are rare, except in Slovenia it has also been reported in neighboring Sebria and Hungary (Jevremovic and Paunovic, 2011; Pleško et al., 2012; Czotter et al., 2018). In our study, RBDV was found only in “Laški rizling” variety (005 library). It was found in this variety in all Slovenian wine-growing regions (Mavrič Pleško et al., 2009). Complete or almost complete reference coverage of both RNA1 (100%) and RNA2 (98.7%) was obtained (Table 1). All four samples were confirmed positive with RT-PCR and partial RNAs2 were Sanger sequenced and deposited in NCBI (GenBank accession no. OK139039-OK139042). The MP sequences of our isolates shared 100% nt identity (100% aa identity), while in the CP gene region they shared 98.18–99.55% nt identities (97.95–100% aa identities). Considering phylogenetic analysis of partial sequences of the CP gene (438 bp), our isolates clustered among other isolates of Vitis sp. retrieved from NCBI and they were clearly separated from isolates of Rubus sp. (Supplementary Figure 1), which was also reported in other studies (Mavrič Pleško et al., 2009, 2020; Valasevich et al., 2011).
Grapevine Pinot Gris Virus
GPGV is an emerging virus associated with grapevine leaf mottling and deformation (GLMD) disease (Giampetruzzi et al., 2012), but has not yet been included in certification programs in Europe. In Slovenia, the first symptoms (shortened internodes, mottling, deformation, and poor leaf development) were observed in 2001, and samples were tested for eight viruses (ArMV, CLRV, GFLV, RBDV, SLRSV, TBRV, ToRSV, and TRSV), but none was confirmed by DAS-ELISA (Mavrič Pleško et al., 2014). In 2012, GPGV was discovered in Italy using sRNA-seq (Giampetruzzi et al., 2012), and in 2014 its occurrence was reported also from Slovenia (Mavrič Pleško et al., 2014). These authors reported that GPGV seems to be widespread in the Primorska wine-growing region, but they also observed it in different parts of Slovenia. In addition to Italy, the virus seems to be widespread in other neighboring countries, Hungary (17 out of 18 libraries) and Croatia (61.97%) (Czotter et al., 2018; Hančević et al., 2021). Based on HTS data, GPGV was the most prevalent virus in our study. It was found in all 12 libraries (95.3–100% coverage of complete reference sequences) (Table 1). Seventy-two out of 79 samples were positive (91.14%) (Figure 1). Forty sequences were generated and deposited in NCBI (GenBank accession no. OK139043-OK139082). A polymorphism showing C/T variation introducing a premature termination codon was observed in the MP sequence. The C/T polymorphism was observed in 13 samples making MP shorter by 18 nt (6 aa). This polymorphism was also observed for isolates analyzed in some other studies (Reynard, 2015; Saldarelli et al., 2015; Czotter et al., 2018; Abou Kubaa et al., 2020). The MP sequences of 40 Slovenian isolates shared nt identities of 93.94–100% (87.79–100% aa identities). 40 Slovenian CP sequences shared pairwise nt identities of 94.53–100% (97.86–100% aa identities). The phylogenetic tree was constructed using part of the sequences of the MP gene and CP gene (718 nt) and it showed partitioning of our isolates into two clades with isolates from geographically relatively close countries (Supplementary Figure 2).
Figure 1. Prevalence of nine viruses and two viroids in analyzed samples.
Grapevine Rupestris Stem Pitting-Associated Virus
Using the HTS approach, GRSPaV (member of the rugose wood complex) was detected in 10 libraries (Table 1). In eight libraries complete reference sequences were covered 48.4–99.1%, while in other two libraries, 012 (“Zeleni Sauvignon” variety) and 015 (“Malvazija” variety), only partial genome sequences were covered 62.3 and 81.4%, respectively. In all libraries GRSPaV had a low sequencing depth (5.3X–16.6X). It was confirmed in all libraries where it was predicted with RT-PCR, and even in two other libraries (007 and 014), where it was not detected by sRNA-seq/VirusDetect pipeline. In library 007 (“Rebula” variety) all eight samples were infected, while in library 014 (“Malvazija” variety) only two samples were infected. In Hungary, using sRNA-seq approach, the same contradictory results were reported by three independent studies (Czotter et al., 2018; Demian et al., 2020; Turcsan et al., 2020). The authors indicated that this could be due to technical issues, the possibility that concentrations were under detection threshold due to the bulk sequencing strategy or deeper biological aspects, as a possible long coexistence between grapevine and GRSPaV resulted in mutual adaptation (Gambino et al., 2012), and potentially the plant immune response was not activated. The reason why GRSPaV was not detected by sRNA-seq in some libraries in different studies conducted in different countries requires further studies. Overall, 70 out of 79 samples were tested positive (88.61%), making GRSPaV the second most abundant virus (Figure 1). The selected RT-PCR products were directly sequenced, but due to different genetic variants in the same sample, we were not able to generate high quality sequences, therefore cloning was performed. One RT-PCR product was selected from each variety. Three white colonies were randomly picked from each variety, and after colony PCR, three products were obtained for Laški rizling 3/45B, Rebula 16/1B, Zeleni Sauvignon 14/2P, and Malvazija 20/48P, and one product for Refošk 10/3B and Pokalca 9/27G. In total, 14 products were sequenced (GenBank accession no. OK138921–OK138934), and all were different from each other, even when originating from the same plant. The highest overall mean distance was revealed for three variants of Laški rizling 3/45B (17.14%), and the lowest overall mean distance for Malvazija 20/48P (6.32%). While the overall mean distance for all 14 sequenced variants was 14.06%. It can be concluded that at least three genetic variants exist in the selected samples which differ extensively in the analyzed genome region. The high genetic diversity could be due to the lack of proofreading activity of RdRp, errors in genome replication, frequent recombination, and grafting of individual plants onto differentially infected rootstocks and scions (Glasa et al., 2017). The phylogenetic tree also showed that the genetic variants, even if from the same plant, clustered in different clades (Supplementary Figure 3).
Grapevine Fanleaf Virus
GFLV, responsible for a fanleaf degeneration disease and one of the viruses causing the most significant damages on vines (Andret-Link et al., 2004), was detected in two of our libraries, 013 (“Zeleni Sauvignon” variety) and 016 (“Pokalca” variety) (Table 1). RNA1 was covered 92 and 71.7%, respectively, while RNA2 was covered 90.7 and 92.1%, respectively. Validation by RT-PCR using published primers resulted in one sample positive from bulk 013 and three samples from bulk 016. All positive samples were sequenced (GenBank accession no. OK139035–OK139038). Three isolates from “Pokalca” variety shared 99.67 or 99.84% nt identities (99.5 or 100% aa identities), whereas the isolate Zeleni Sauvignon 16/3P differed greatly from the Pokalca isolates with 87.27 or 87.44% nt identities (96.49 or 96.99% aa identities). In addition to the differences between our isolates they also differed from the isolates deposited in NCBI. Pokalca isolates shared the highest identity, 90.62 or 90.79%, with the isolate from France (MG731624), while Zeleni Sauvignon 16/3P shared the highest identity 91.29% with the isolate from Switzerland (MG731616). High sequence variation between GFLV isolates of partial or complete RNA2 (2AHP, 2BMP, and 2CCP) has been reported in several studies (Naraghi-Arani et al., 2001; Fattouch et al., 2005; Pompe-Novak et al., 2007; Elbeaino et al., 2014). This virus does not possess proofreading activity of RdRp and the large genetic variability indicates that the GFLV genome consists of quasispecies populations (Naraghi-Arani et al., 2001). Phylogenetic analysis showed that our isolates differed in the region of partial RNA2 from isolates retrieved from NCBI, even from the previously characterized isolates from Slovenia (Pompe-Novak et al., 2007), but they were the closest to the isolates from Italy and France (Supplementary Figure 4). Due to the differences observed among sequences, we designed new primers based on our sequencing data and repeated analysis. Positive amplifications were observed for additional two samples from bulk 013 and for four samples from bulk 016.
Grapevine Leafroll-Associated Virus 3
GLRaV-3 is the major causal agent of one of the most detrimental grapevine diseases named as grapevine leafroll disease (GLD) (Maree et al., 2013). GLRaV-3 was detected only in library 010 (“Refošk” variety). The reference sequence GQ352631 was 98.6% covered and a sequencing depth of 39.6X (Table 1). With the primer pair amplifying the CP gene region, only one sample (Refošk 11/4P) was positive, therefore we used the primer pair amplifying the HSP70h gene region and the same result was obtained. The sequencing of CP gene region of Refošk 11/4P (GenBank accession no. OK138920) showed the highest nt identity (99.76%) with 15 sequences originating from Greece, Portugal, United States, Canada and Pakistan and 3 sequences of unknown origin. Phylogenetic analysis showed that our isolate clustered together with the isolate from Portugal (Supplementary Figure 5). GLRaV-3 was the least prevalent virus in our sample set. Other grapevine leafroll-associated viruses, members of Ampelovirus or Closterovirus genus, were not detected. The main reason for the lack of detection of leafroll-associated viruses can be explained by the fact that after the mass selection, all selected vines (potential preclonal vines—ELITE) were screened with ELISA tests, which have a fairly good detection for viruses of GLD, therefore at that step, all infected vines were excluded by further selections and propagation. We are aware that serological tests can have quite large deviations in the detection of viruses, but in this case, it seems that we have been quite successful in the leafroll-associated viruses detection with the ELISA test.
Grapevine Fleck Virus
GFkV was detected in eight libraries. In five libraries (005, 007, 008, 010, and 012) complete genome reference sequences (AJ309022 or KT000362) were covered with 76.9–82.1%, while in two libraries of “Malvazija” variety (014 and 015), reference sequence AJ309022 was covered only with 30.4 and 32.2%, respectively, and in library 011 partial sequence KF417610 was covered with 33.3%. GFkV was validated in all predicted libraries and for 34 samples we got positive RT-PCR result. All samples were sequenced, but the results showed that seven products belonged to GRVFV, which is consistent with reports from Czotter et al. (2018) from Hungary, indicating high similarity between these two viruses and possible cross-amplification with primers. Sequences of two samples, Laški rizling 3/56B and 3/64B had lower quality, therefore they were excluded from further analysis. Twenty-five GFkV sequences were generated and deposited in NCBI (GenBank accession no. OK139010–OK139034). They shared nt identities of 91.6–100% (93.7–100% aa identities). GFkV is phloem-limited, not mechanically transmissible, and spreads by grafting and infected propagating material (Sabanadzovic et al., 2017). Our isolates shared the highest nt identities with isolates from Bosnia and Herzegovina, North Macedonia, Hungary and the United States. Phylogenetic studies showed that the sequenced isolates clustered in different clades depending on the variety, except “Laški rizling,” with isolates from neighboring countries, while all samples of “Refošk” variety cluster together with isolate from the United States (Supplementary Figure 6). A few decades ago the predecessors of our samples were grafted onto untested rootstocks imported from neighboring countries and from Davis University in California (Hrček, 1977). It seems that with the rootstocks GFkV was imported, but also a lot of grafts produced in Slovenia were exported in neighboring countries, especially in former Yugoslavia.
In addition to GFkV, three fleck-like viruses (GRGV, GRVFV, and GSyV-1) were detected for the first time in Slovenia (paper in review).
Hop Stunt Viroid
HSVd has a wide natural host range from different botanical families. Slovenia is one of the major hop producers, and the first report that HSVd infects hops in Slovenia was published in 2012 (Radisek et al., 2012), while it was recently confirmed on grapevines in co-infection with GV-Sat, GLRaV-1, GLRaV-2, GRSPaV, GPGV, and GYSVd-1 (Miljanić et al., 2021). According to the HTS results, HSVd was present in all libraries (Table 1). In ten libraries the reference sequence KJ810551 was covered 100%, while in the other two libraries (008 and 014) the reference sequence KY508372 was covered 99.4% (Table 1). It was validated by RT-PCR and all samples were positive. Forty complete genome sequences were generated and deposited in NCBI (GenBank accession no. OK138935–OK138974). Thirty-eight isolates were 100% identical to each other, while two isolates (Pokalca 3/4P and Pokalca 3/6P) were identical and shared 98% identities with other isolates. In the genome of Pokalca 3/4P and 3/6P isolates, insertions were observed at positions 123 and 257, while SNPs were observed at positions 171, 172, 238, 244, 259, and 260 (Supplementary Figure 9). Also, phylogenetic tree showed that this two isolates clustered completely different from other isolates (Supplementary Figure 7).
Grapevine Yellow Speckle Viroid-1
According to our analysis GYSVd-1 was found in all analyzed libraries with 100% references coverage (Table 1). According to RT-PCR, 71 samples were positive (89.87%) (Figure 1). To our knowledge, nine GYSVd-1 sequences of Slovenian autochthonous grapevine varieties have been deposited in NCBI so far (Štajner et al., 2019). In this study, 35 complete genome sequences were generated and deposited in NCBI (GenBank accession no. OK138975–OK139009). GYSVd-1 was less prevalent and showed higher genetic diversity than HSVd. Slovenian GYSVd-1 sequences had 95.35–100% nt identities. Multiple alignment with ClustalW revealed InDel mutations at four positions in the genome (63, 92, 163, and 287) (Supplementary Figure 10). Phylogenetic analysis showed that analyzed GYSVd-1 isolates clustered in different clades regardless of variety (Supplementary Figure 8).
Co-infections in Analyzed Samples
GPGV, GRSPaV, HSVd, and GYSVd-1 were the most prevalent in our sample set (Figure 1) and their co-infection were the most common (18.99%) (Figure 2). The second most prevalent co-infections were GPGV + GRSPaV + GRVFV + HSVd + GYSVd-1 and GPGV + GRSPaV + GFkV + GRVFV + HSVd + GYSVd-1 (16.46%) (Figure 2). There were no vines that were free of viruses or viroids, and there were no vines that were infected with only one viral entity. The highest number of tested plants were infected with five viral entities (37.97%), followed by six (24.05%) and four viral entities (22.78%), while one sample (Laški rizling 3/45B) was infected with eight viruses/viroids (Figure 3).
Figure 2. Percentages of co-infections in analyzed samples.
Figure 3. Percentages of mixed infections in analyzed samples.
Conclusion
The main advantage of using the HTS approach is the complete insight into virome of the analyzed samples (Czotter et al., 2018). When screening the virome status of selected plants the HTS approach is considered method of choice. The HTS approach used for virome screening is mainly based on bulked samples, which is cost effective, because in the analysis usually a lot of samples are included, and the main limitation is the possibility that due to the bulk sequencing strategy viral concentrations may be under detection threshold in some cases. In our study all individual samples were tested with RT-PCR for each HTS predicted infection, and all obtained results were consistent, except for GRSPaV, but due to the fact that similar results related to inconsistent detection of GRSPaV were obtained in different studies, it may have deeper biological aspect and required further analysis which are discussed.
The present study gives us a detailed insight into the virome status of preclonal candidates of autochthonous and local grapevine varieties in the Primorska wine-growing region of Slovenia. In this study significant number of sequences were generated for different viral pathogens and could further improve their routine diagnostics, which is especially important as they cannot be controlled by conventional plant protection methods.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
VM: HTS and bioinformatics analysis, validation of HTS data, data analysis, and writing the original draft. UK: HTS and bioinformatics analysis. DR: provided the plant material, review and editing of the original draft. AŠ: provided the plant material. JJ and NŠ: experimental design, review and editing of the original draft. All authors contributed to the article and approved the submitted version.
Funding
This research was funded by the Slovenian Research Agency (SRA–ARRS), grant no. P4-0077—research program: Genetics and Modern Technologies of Crops awarded to the research group at Chair of Genetics, Biotechnology, Statistics and Plant breeding at the University of Ljubljana, Biotechnical Faculty, and scholarship for nationals of Western Balkans to postgraduate study in Slovenia (Public Scholarship, Development, Disability, and Maintenance Fund of the Republic of Slovenia, 245. Public Announcement) awarded to VM.
Conflict of Interest
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.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We would like to acknowledge to the Slovenian Research Agency (SRA–ARRS), grant number P4-0077-research program: Genetics and Modern Technologies of Crops and Public Scholarship, Development, Disability and Maintenance Fund of the Republic of Slovenia.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.830866/full#supplementary-material
Footnotes
References
Abou Kubaa, R., Choueiri, E., Jreijiri, F., El Khoury, Y., and Saldarelli, P. (2020). First report of grapevine Pinot gris virus in Lebanon and the middle East. J. Plant Pathol. 102:565. doi: 10.1007/s42161-019-00453-w
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
Alabi, O. J., Casassa, L. F., Gutha, L. R., Larsen, R. C., Henick-Kling, T., Harbertson, J. F., et al. (2016). Impacts of grapevine leafroll disease on fruit yield and grape and wine chemistry in a wine grape (Vitis vinifera L.) cultivar. PLoS One 11:e0149666. doi: 10.1371/journal.pone.0149666
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. doi: 10.4454/jpp.v86i3.987
Bartel, D. P. (2004). MicroRNAs: genomics, biogenesis, mechanism, and function. Cell 116, 281–297. doi: 10.1016/S0092-8674(04)00045-5
Bernstein, E., Caudy, A. A., Hammond, S. M., and Hannon, G. J. (2001). Role for a bidentate ribonuclease in the initiation step of RNA interference. Nature 409, 363–366. doi: 10.1038/35053110
Bouyahia, H., Boscia, D., Savino, V., La Notte, P., Pirolo, C., Castellano, M. A., et al. (2005). Grapevine rupestris stem pitting-associated virus is linked with grapevine vein necrosis. Vitis J. Grapevine Res. 44, 133–137.
Credi, R., Babini, A. R., Betti, L., Bertaccini, A., and Gelli, C. (1981). A distinctive isolate of strawberry latent ringspot virus from grapevines in Italy. Phytopathologia 20, 56–63.
Czotter, N., Molnar, J., Szabó, E., Demian, E., Kontra, L., Baksa, I., et al. (2018). NGS of virus-derived small RNAs as a diagnostic method used to determine viromes of Hungarian Vineyards. Front. Microbiol. 9:122. doi: 10.3389/fmicb.2018.00122
Demian, E., Jaksa-Czotter, N., Molnar, J., Tusnady, G. E., Kocsis, L., and Varallyay, E. (2020). Grapevine rootstocks can be a source of infection with non-regulated viruses. Eur. J. Plant Pathol. 156, 897–912. doi: 10.1007/s10658-020-01942-w
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
Elbeaino, T., Kiyi, H., Boutarfa, R., Minafra, A., Martelli, G. P., and Digiaro, M. (2014). Phylogenetic and recombination analysis of the homing protein domain of grapevine fanleaf virus (GFLV) isolates associated with ‘yellow mosaic’ and ‘infectious malformation’ syndromes in grapevine. Arch. Virol. 159, 2757–2764. doi: 10.1007/s00705-014-2138-8
Fajardo, T. V. M., Bertocchi, A. A., and Nickel, O. (2020). Determination of the grapevine virome by high-throughput sequencing and grapevine viruses detection in Serra Gaucha, Brazil. Rev. Ceres 67, 156–163. doi: 10.1590/0034-737X202067020010
Fattouch, S., Acheche, H., M’Hirsi, S., Mellouli, L., Bejar, S., Marrakchi, M., et al. (2005). RT-PCR-RFLP for genetic diversity analysis of Tunisian virus isolates in their natural host plants. J. Virol. Methods 127, 126–132. doi: 10.1016/j.jviromet.2005.03.008
Fuchs, M. (2020). Grapevine viruses: a multitude of diverse species with simple but overall poorly adopted management solutions in the vineyard. J. Plant Pathol. 102, 643–653. doi: 10.1007/s42161-020-00579-2
Gambino, G., Cuozzo, D., Fasoli, M., Pagliarani, C., Vitali, M., Boccacci, P., et al. (2012). Co-evolution between Grapevine rupestris stem pitting-associated virus and Vitis vinifera L. leads to decreased defence responses and increased transcription of genes related to photosynthesis. J. Exp. Bot. 63, 5919–5933. doi: 10.1093/jxb/ers244
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 Res. 163, 262–268. doi: 10.1016/j.virusres.2011.10.010
Girardello, R. C., Cooper, M. L., Lerno, L. A., Brenneman, C., Eridon, S., Sokolowsky, M., et al. (2020). Impact of grapevine red blotch disease on cabernet sauvignon and merlot wine composition and sensory attributes. Molecules 25:3299. doi: 10.3390/molecules25143299
Glasa, M., Predajňa, L., Šoltys, K., Sihelská, N., Nagyová, A., Wetzel, T., et al. (2017). Analysis of grapevine rupestris stem pitting-associated virus in Slovakia reveals differences in intra-host population diversity and naturally occurring recombination events. Plant Pathol. J. 33, 34–42. doi: 10.5423/PPJ.OA.07.2016.0158
Hančević, K., Saldarelli, P., Čarija, M., Černi, S., Zdunić, G., Mucalo, A., et al. (2021). Predominance and diversity of GLRaV-3 in native vines of mediterranean Croatia. Plants 10, 1–14. doi: 10.3390/plants10010017
Jakse, J., Radisek, S., Pokorn, T., Matousek, J., and Javornik, B. (2015). Deep-sequencing revealed Citrus bark cracking viroid (CBCVd) as a highly aggressive pathogen on hop. Plant Pathol. 64, 831–842. doi: 10.1111/ppa.12325
Jevremovic, D., and Paunovic, S. (2011). Raspberry bushy dwarf virus: a grapevine pathogen in Serbia. Pestic. Fitomed. 26, 55–60. doi: 10.2298/pif1101055j
Kashif, M., Pietilä, S., Artola, K., Jones, R. A. C., Tugume, A. K., Mäkinen, V., et al. (2012). Detection of viruses in sweetpotato from Honduras and Guatemala augmented by deep-sequencing of small-RNAs. Plant Dis. 96, 1430–1437. doi: 10.1094/PDIS-03-12-0268-RE
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
Kumar, S., Stecher, G., Li, M., Knyaz, C., and Tamura, K. (2018). MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 35, 1547–1549. doi: 10.1093/molbev/msy096
Lee, J., and Martin, R. R. (2009). Influence of grapevine leafroll associated viruses (GLRaV-2 and -3) on the fruit composition of oregon Vitis vinifera L. cv. pinot noir: phenolics. Food Chem. 112, 889–896. doi: 10.1016/j.foodchem.2008.06.065
Lee, J., Rennaker, C. D., Thompson, B. D., and Karasev, A. V. (2021). Influence of grapevine red blotch virus (GRBV) on idaho ‘syrah’ grape composition. Sci. Hortic. 282:110055. doi: 10.1016/j.scienta.2021.110055
Li, H., Wei, L., Qin, C., Cheng, J., Zhang, X., Chen, W., et al. (2021). Characterization of viruses and viroids in Vitis vinifera ‘Kyoho’ in Hangzhou China by small RNA deep sequencing and molecular detection. J. Hortic. Sci. Biotechnol. 96, 400–406. doi: 10.1080/14620316.2020.1845985
Lunden, S., Meng, B., Avery, J., and Qiu, W. (2010). Association of Grapevine fanleaf virus, Tomato ringspot virus and Grapevine rupestris stem pitting-associated virus with a grapevine vein-clearing complex on var. Chardonnay. Eur. J. Plant Pathol. 126, 135–144. doi: 10.1007/s10658-009-9527-y
Maree, H. J., Almeida, R. P. P., Bester, R., Chooi, K. M., Cohen, D., Dolja, V. V., et al. (2013). Grapevine leafroll-associated virus 3. Front. Microbiol. 4:82. doi: 10.3389/fmicb.2013.00082
Martin, M. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet.J. 17, 10–12. doi: 10.14806/ej.17.1.200
Mavrič Pleško, I., Lamovšek, J., Lešnik, A., and Viršèek Marn, M. (2020). Raspberry bushy dwarf virus in Slovenia - geographic distribution, genetic diversity and population structure. Eur. J. Plant Pathol. 158, 1033–1042. doi: 10.1007/s10658-020-02115-5
Mavrič Pleško, I., Viršček Marn, M., Seljak, G., and Žežlina, I. (2014). First report of grapevine Pinot gris virus infecting grapevine in Slovenia. Plant Dis. 98:1014. doi: 10.1094/PDIS-11-13-1137-PDN
Mavrič Pleško, I., Viršček Marn, M., Širca, S., and Urek, G. (2009). Biological, serological and molecular characterisation of Raspberry bushy dwarf virus from grapevine and its detection in the nematode Longidorus juvenilis. Eur. J. Plant Pathol. 123, 261–268. doi: 10.1007/s10658-008-9362-6
Mavrič, I., Marn, M. V., Koron, D., and Žežlina, I. (2003). First Report of Raspberry bushy dwarf virus on Red Raspberry and Grapevine in Slovenia. Plant Dis. 87, 1148–1148. doi: 10.1094/pdis.2003.87.9.1148d
Miljanić, V., Jakše, J., Beber, A., Rusjan, D., Škvarč, A., and Štajner, N. (2021). First report of grapevine satellite virus in Slovenia. J. Plant Pathol. 98:1014. doi: 10.1007/s42161-021-00902-5
Naraghi-Arani, P., Daubert, S., and Rowhani, A. (2001). Quasispecies nature of the genome of Grapevine fanleaf virus. J. Gen. Virol. 82, 1791–1795. doi: 10.1099/0022-1317-82-7-1791
Navarro, B., Pantaleo, V., Gisel, A., Moxon, S., Dalmary, T., Bisztray, G., et al. (2009). Deep sequencing of viroid-derived small RNAs from grapevine provides new insights on the role of RNA silencing in plant-viroid interaction. PLoS One 4:e7686. doi: 10.1371/journal.pone.0007686
Pleško, I. M., Marn, M. V., Nyerges, K., and Lázár, J. (2012). First report of raspberry bushy dwarf virus Infecting Grapevine in Hungary. Plant Dis. 96, 1582–1582. doi: 10.1094/pdis-04-12-0383-pdn
Pompe-Novak, M., Gutiérrez-Aguirre, I., Vojvoda, J., Blas, M., Tomažiè, I., Vigne, E., et al. (2007). Genetic variability within RNA2 of Grapevine fanleaf virus. Eur. J. Plant Pathol. 117, 307–312. doi: 10.1007/s10658-006-9096-2
Radisek, S., Majer, A., Jakse, J., Javornik, B., and Matoušek, J. (2012). First report of hop stunt viroid infecting hop in Slovenia. Plant Dis. 96, 592–592. doi: 10.1094/pdis-08-11-0640-pdn
Reynard, J. (2015). “Survey of emerging viruses in Switzerland,” in Proceedings of the 18th Congress of the International Council for the Study of Virus and Virus-like Diseases of the Grapevine, Ankara, 223–224.
Sabanadzovic, S., Aboughanem-Sabanadzovic, N., and Martelli, G. P. (2017). “Grapevine fleck and similar viruses,” in Grapevine Viruses: Molecular Biology, Diagnostics and Management, eds B. Meng, G. Martelli, D. Golino, and M. Fuchs (Cham: Springer), 331–349. doi: 10.1007/978-3-319-57706-7_16
Saldarelli, P., Giampetruzzi, A., Morelli, M., Malossini, U., Pirolo, C., Bianchedi, P., et al. (2015). Genetic variability of Grapevine Pinot gris virus and its association with Grapevine leaf mottling and deformation. Phytopathology 105, 555–563. doi: 10.1094/PHYTO-09-14-0241-R
Singh, K., Jarošova, J., Fousek, J., Chen, H., and Kundu, J. K. (2020). Virome identification in wheat in the Czech Republic using small RNA deep sequencing. J. Integr. Agric. 19, 1825–1833. doi: 10.1016/S2095-3119(19)62805-4
Štajner, N., Beber, A., Rusjan, D., and Jakše, J. (2019). Occurrence of grapevine yellow speckle viroid 1 in autochthonous grapevines in Slovenia. J. Plant Pathol. 101:397. doi: 10.1007/s42161-018-0170-3
Sudarshana, M. R., Perry, K. L., and Fuchs, M. F. (2015). Grapevine red blotch-associated virus, an emerging threat to the grapevine industry. Phytopathology 105, 1026–1032. doi: 10.1094/PHYTO-12-14-0369-FI
Thompson, J. D., Higgins, D. G., and Gibson, T. J. (1994). CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 22, 4673–4680. doi: 10.1093/nar/22.22.4673
Torregrosa, L., Vialet, S., Adivèze, A., Iocco-Corena, P., and Thomas, M. R. (2015). Grapevine (Vitis vinifera L.). Methods Mol. Biol. 1224, 177–194. doi: 10.1007/978-1-4939-1658-0_15
Turcsan, M., Demian, E., Varga, T., Jaksa-Czotter, N., Szegedi, E., Olah, R., et al. (2020). Hts-based monitoring of the efficiency of somatic embryogenesis and meristem cultures used for virus elimination in grapevine. Plants 9, 1–10. doi: 10.3390/plants9121782
Valasevich, N., Kukharchyk, N., and Kvarnheden, A. (2011). Molecular characterisation of raspberry bushy dwarf virus isolates from Sweden and belarus. Arch. Virol. 156, 369–374. doi: 10.1007/s00705-010-0912-9
Vega, A., Gutiérrez, R. A., Peña-Neira, A., Cramer, G. R., and Arce-Johnson, P. (2011). Compatible GLRaV-3 viral infections affect berry ripening decreasing sugar accumulation and anthocyanin biosynthesis in Vitis vinifera. Plant Mol. Biol. 77, 261–274. doi: 10.1007/s11103-011-9807-8
Vives, M. C., Velázquez, K., Pina, J. A., Moreno, P., Guerri, J., and Navarro, L. (2013). Identification of a new enamovirus associated with citrus vein enation disease by deep sequencing of small RNAs. Phytopathology 103, 1077–1086. doi: 10.1094/PHYTO-03-13-0068-R
Keywords: Vitis vinifera L., preclonal candidates, HTS, viruses, viroids
Citation: Miljanić V, Jakše J, Kunej U, Rusjan D, Škvarč A and Štajner N (2022) Virome Status of Preclonal Candidates of Grapevine Varieties (Vitis vinifera L.) From the Slovenian Wine-Growing Region Primorska as Determined by High-Throughput Sequencing. Front. Microbiol. 13:830866. doi: 10.3389/fmicb.2022.830866
Received: 07 December 2021; Accepted: 10 January 2022;
Published: 21 February 2022.
Edited by:
Yi Xu, Nanjing Agricultural University, ChinaReviewed by:
Pasquale Saldarelli, Institute for Sustainable Plant Protection, National Research Council (CNR), ItalyKeith Perry, Cornell University, United States
Copyright © 2022 Miljanić, Jakše, Kunej, Rusjan, Škvarč and Štajner. 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.
*Correspondence: Nataša Štajner, bmF0YXNhLnN0YWpuZXJAYmYudW5pLWxqLnNp