Next Article in Journal
Morphological and Molecular Identification of Dactylogyrus gobiocypris (Monogenea: Dactylogyridae) on Gills of a Model Fish, Gobiocypris rarus (Cypriniformes: Gobionidae)
Next Article in Special Issue
Discovery and Genome Characterization of a Closterovirus from Wheat Plants with Yellowing Leaf Symptoms in Japan
Previous Article in Journal
Taenia crassiceps Cysticercosis in a Wild Muskrat and a Domestic Dog in the Northeastern United States
Previous Article in Special Issue
Genetic Diversity of Viral Populations Associated with Ananas Germplasm and Improvement of Virus Diagnostic Protocols
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Pathogens
Volume 12
Issue 2
10.3390/pathogens12020205
pathogens-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Analysis of Hop Stunt Viroid Diversity in Grapevine (Vitis vinifera L.) in Slovakia: Coexistence of Two Particular Genetic Groups

by
Peter Alaxin
1,2,
Lukáš Predajňa
1,
Adam Achs
1,3,
Zdeno Šubr
1,
Michaela Mrkvová
2 and
Miroslav Glasa
1,2,*
1
Institute of Virology, Biomedical Research Center of Slovak Academy of Sciences, Dúbravská cesta 9, 845 05 Bratislava, Slovakia
2
Department of Biology, Faculty of Natural Sciences, University of Ss. Cyril and Methodius in Trnava, Námestie J. Herdu 2, 917 01 Trnava, Slovakia
3
Department of Molecular Biology, Faculty of Natural Sciences, Comenius University in Bratislava, Ilkovičova 6, 842 15 Bratislava, Slovakia
*
Author to whom correspondence should be addressed.
Pathogens 2023, 12(2), 205; https://doi.org/10.3390/pathogens12020205
Submission received: 19 December 2022 / Revised: 25 January 2023 / Accepted: 26 January 2023 / Published: 28 January 2023
(This article belongs to the Special Issue Plant Virus Genome Diversity in Plant Hosts and Insect Vectors)
Browse Figures
Versions Notes

Abstract

:
The hop stunt viroid (HSVd) is a widespread subviral pathogen infecting a broad spectrum of plant hosts including grapevine (Vitis vinifera L.). Despite its omnipresence in virtually all grapevine growing areas around the world, molecular data characterizing HSVd populations are missing from Slovakia. Analysis of the complete nucleotide sequences of 19 grapevine variants revealed the existence of two genetic HSVd groups in Slovakia (internally named the “6A” and “7A” groups based on the particular stretch of adenines at nucleotide positions 39–44/45, respectively). Despite their sampling at different times in various unrelated vineyards, the 6A and 7A groups are characterized by low intra-group divergence (~0.3 and 0.2%, respectively). On the other hand, inter-group divergence reached 2.2% due to several mutations, seven of which were found to be group-specific and mainly (except for one) located in the region of the pathogenic domain. Interestingly, in addition to their frequent co-existence within the same geographical location, the mixed infection of the 6A and 7A type sequence variants was also unequivocally and repeatedly proven within single grapevine plants. The RNA secondary structure analysis of representative isolates from each of these two genetic groups indicated a potential compensatory explanation of such mutations. These group-specific sites could be pointing towards the evolutionary selection linked to the necessity of the viroid to retain its structural conformational integrity, crucial for its functional biochemical ability to interact with specific grapevine cellular host factors required for HSVd propagation.

1. Introduction

Viroids represent a small (regarding the overall number of species [1]) but very widespread and successful group of phytopathogens, known only from higher plants, characterized by a very simple molecular composition. Categorized as subviral infectious agents, their whole structure is constituted of only one single-stranded RNA molecule, with covalently bound ends forming a circle [2,3]. Viroid “genomes” do not code for any gene [4], and thus are lacking any of their own protein or lipid components. Despite this fact, these infectious RNA molecules can still have a significant impact on overall plant physiology. Compared to their structurally more complex viral counterparts, viroids can induce a wide range of symptoms on plants, often resulting in important agronomical losses reported from different countries around the world [5,6,7,8,9]. Viroids owe their remarkable biochemical activity to the secondary structure of their RNA molecules (forming even three-dimensional functional motifs [10]), which are folded because of intramolecular base-pairing (conventional Watson–Crick as well as wobble). The resulting conformation thus basically functions as double-stranded RNA, which is a hypothesized model for viroids’ ability to hijack eucaryotic cellular proteins (some of which even normally do not recognize RNA as their substrate) in order to exploit them for their own replication and cell-to-cell transport [11,12,13,14].
The hop stunt viroid (HSVd, a type species of the genus Hostuviroid) belongs to the Pospiviroidae family, the larger one of the two viroid families recognized to date [15]. A typical characteristic of Pospiviroidae members is the localization in the cellular nucleus regarding their replication strategy. This is in contrast with the chloroplast-targeting and even own-ribozymal activity possessed in the family Avsunviroidae [1,3,15]. The HSVd has a broad host range, thus posing as a potential threat to agriculture in the world. Its first description was triggered by the study of the causal agent of shoot growth dwarfism (stunt symptom) in hop (Humulus lupulus L.; Cannabaceae) in Japan [16,17]. The HSVd also includes four other sequence variants from mutually different hosts, named after their host-specific symptomatology, i.e., the cucumber pale fruit viroid (CPFVd; previously regarded as a separate species [18,19]); the citrus cachexia viroid (CCaVd, the causal agent of the so-called cachexia xyloporosis disease of citrus trees [20,21,22,23]); the dapple peach (fruit disease) viroid, and the dapple plum (fruit disease) viroid in Prunus spp. [24,25,26,27]. Moreover, several other agriculturally and economically important plants can also be infected with HSVds, e.g., apricot, almond (Rosacea), fig, mulberry (Moraceae), pistachio (Anacardiaceae), and kumquat (Rutaceae) [27,28,29,30]. However, in many cases, the HSVd infection can stay asymptomatic [22,27].
Grapevine (Vitis vinifera L.; Vitaceae) is another well-known, widespread, and often symptomless natural HSVd host [31,32]. Nonetheless, the research on viroids in grapevines in Slovakia has previously received little attention [33]. As for other pathogens, the diversity of HSVd populations plays a significant role in the epidemiology and etiology of the viroid disease, and its understanding constitutes an important prerequisite to develop accurate diagnostics and efficient control strategies. Therefore, in this work, standard and high-throughput sequencing (HTS) tools were used to analyze HSVd molecular variability in Slovakia, leading to the identification of two grapevine-specific molecular groups of variants, acting often in mixed infections.

2. Materials and Methods

2.1. Grapevine Samples

The HSVd variants characterized in this work were obtained from grapevine plants (V. vinifera L.) sampled in different unrelated vineyards and from different cultivars from 2012 to 2021 as a part of the effort to monitor the health status of grapevines in Slovakia (Table 1). Targeted plants were selected randomly with no specific relation to plant symptomatology. Total RNAs were then extracted either from cortical scrapings of one-year-old canes (when collected in the dormant period) or leaf petioles (from leaves collected during the vegetation period) using a SpectrumTM Plant Total RNA Kit (Sigma-Aldrich, St. Louis, MO, USA) and stored at −80 °C prior to further analysis.

2.2. High-Throughput Sequencing and Analysis

Ribosomal-RNA-depleted preparations (Ribo-Zero rRNA Removal Kit, Illumina, San Diego, CA, USA) from three grapevine samples (SK30, SK503, and SK933; Table 1) were used for double-stranded cDNA synthesis using a SuperScript II Kit (Thermo Fisher Scientific, Waltham, MA, USA). The samples were further processed with a Transposon-based Chemistry Library Preparation Kit (Nextera XT, Illumina, San Diego, CA, USA) and the cDNA library was sequenced on the Illumina MiSeq platform (Illumina, San Diego, CA, USA; 300-bp paired-end sequencing) as previously described [34,35]. From 1.1 to 2.4 million high-quality reads per sample were used for de novo assembling using CLC Genomics Workbench 7.5 and Geneious v. 8.1.9 software. The obtained contigs were blasted to the viral/viroid genome database (ftp://ftp.ncbi.nih.gov/genomes/Viruses/all.fna.tar.gz). The complete nucleotide sequences of the HSVd genomes obtained from HTS data reported herein were deposited into the GenBank database under accession numbers MN548397, OP918897, and OP918898.

2.3. RT-PCR Detection, Sanger Sequencing, and Phylogenetic Analysis

Additional total RNAs isolated from grapevine samples were used as templates for a two-step RT-PCR. The first-strand cDNA was synthesized with the reverse transcription of total RNA using hexameric random DNA primers and an AMV Reverse Transcriptase (both from Promega, Madison, WI, USA). An aliquot of cDNA was then used as a template in a standard PCR for the complete viroid sequence amplification using TaKaRa Ex Taq polymerase (Takara Bio Inc., Shiga, Japan) and the HSVd-specific primer pair HSV-78P (5′ aaCCCGGGGCAACTCTTCTC-3′; sense) and HSV-83M (5′ aaCCCGGGGCTCCTTTCTCA-3′; antisense, lowercase letters represent the non-viroid sequence) partially overlapping in the central conserved region (CCR) of the HSVd [31]. The PCR cycling conditions were as follows: initial denaturation 95 °C for 2 min; 40 cycles (denaturation at 95 °C for 30 s; annealing at 57 °C for 30 s; elongation at 72 °C for 20 s); and final elongation at 72 °C for 5 min. Amplified PCR fragments were then analyzed by the standard agarose gel electrophoresis on 5% agarose gel (UltraPureTM Agarose-1000; Thermo Fisher Scientific, Waltham, MA, USA).
Selected amplicons were subsequently purified from the gel using Wizard® SV Gel and PCR Clean-Up System kits (Promega, Madison, WI, USA). Purified PCR products were subjected to Sanger sequencing (Eurofins Genomics, Vienna, Austria) primed by the same HSVd-specific primers used for RT-PCR. To study the intra-isolate variability, two purified PCR products were cloned into pGEM®-T Easy Vector (Promega, Madison, WI, USA) and the individual clones (6 clones per amplicon) were bidirectionally Sanger sequenced using pUC/M13 primers (Promega, Madison, WI, USA). The obtained HSVd complete sequences were deposited into the GenBank database under accession numbers OP918899–OP918914 (Table 1).
The genomic sequences obtained in this work, together with the representative HSVd sequences retrieved from the public database (https://www.ncbi.nlm.nih.gov/nucleotide/) were analyzed using the Molecular Evolution Genetics Analysis (MEGA7) program [36], using the p-distance model.
The prediction of RNA secondary structures of the selected variants was performed using the online web tool Mfold (available at http://www.unafold.org/mfold/applications/rna-folding-form.php) [37], with the default parameters for circular RNAs.

3. Results and Discussion

3.1. Identification of HSVd in the Grapevine Virome

The previous HTS analyses of three grapevine samples (SK30, SK503, and SK933) revealed a complex infection of plants by several viruses, i.e., grapevine leafroll-associated virus-2; grapevine leafroll-associated virus-3 (GLRaV-3); grapevine rupestris stem pitting-associated virus (GRSPaV); grapevine rupestris vein feathering virus (GRVFV); and grapevine Syrah virus-1 (GSyV-1) [38,39,40]. In addition to these several viral pathogens, the HSVd was also identified to be present in the samples [33,40]. In all three samples, the obtained HTS data enabled the reconstruction of the complete genomic sequences of HSVd with high average coverage (17.1–395.4x). Interestingly, although originating from unrelated grapevines sampled over a five-year period, the three HSVd genomes were strictly collinear (297 nts long) and 100% identical. High sequence similarity (identical for some variants even across the time and geographical distribution) can be similarly observed among database grapevine HSVd sequences [32]. Such sequence rigidity could be a consequence of the antiquity of the HSVd pathogen in grapevines as an infectious RNA [31,32] and the vegetative propagation of grapevines and their distribution around the world.
To further assess the occurrence of HSVd in Slovak grapevines (not studied previously), an additional 52 grapevine samples from different vineyards were analyzed with RT-PCR [31], from which 38 in total tested PCR-positive (the HSVd prevalence in tested samples reached 73.08%). This suggested the ubiquitous presence of this viroid species in Slovak vineyards.

3.2. Identification of Two Molecularly Different HSVd Groups

Out of 38 HSVd-positive samples, we further analyzed 14 samples, selected to ensure their diversity based on the vineyard location or cultivar differences (Table 1). From 12 samples, the complete HSVd genomes were obtained with the direct dideoxy Sanger sequencing of amplicons. Contrary to previous HTS-based reconstructed HSVd genomes, the sequences of 12 of these HSVd genomes (OP918899–OP918910) showed a clear length dichotomy with the occurrence of 297 nucleotides, as well as 298-nucleotide-long variants. The multiple alignments of Slovak HSVd genomes confirmed their splitting into two genetically different molecular groups, characterized by a low mean nucleotide intragroup divergence (reaching 0.3% and 0.2%), while the inter-group divergence reached 2.2%.
The first group (consisting of 3 isolates) was characterized by a genome 297 nt in length and a typical stretch of 6 adenosines at positions 39–44 (referred to here as the 6A group, Figure 1 and Figure 2). On the other hand, 9 isolates formed the second group, with a genome consisting of 298 nt and characterized by 7 adenosines at positions 39–45 (referred to here as the 7A group). Interestingly, one of the adenosines at positions 39–45/46 was responsible for the 1 nt difference in length between the two genetic groups.
In addition to the different adenosine stretches, the alignment revealed a total of seven mutations located in the positions spanning the pathogenic and terminal left domains, having a clear group-specific character (Figure 1). This grouping was also confirmed by phylogenetic analyses (Figure 2). However, due to the mixed infection of grapevines involving HSVd and other viral pathogens mentioned above, the contribution of the two HSVd groups to symptomatology (if any) in grapevine could not be determined. The fact that HSVds are reported mostly asymptomatically in grapevine (in contrast to in hop) could be linked to mutual long-term adaptation and the assumption that the HSVd pathogen originates from grapevine as its natural source [31].

3.3. Mixed Infection of HSVd Variants

Although in most of the cases, direct dideoxy Sanger sequencing resulted in unequivocal sequence reads, we were not able to obtain clean whole-genome HSVd sequences from 2 of 14 samples. The visual inspection of chromatograms from directly sequenced amplicons of the SK08-21 and SK11-21 samples showed an occurrence of several mixed signals, interestingly, mainly in the same 7 positions responsible for HSVd grouping in the previously characterized 12 samples (Supplementary Figure S1). To elucidate this problem, we cloned the PCR amplicons and sequenced six individual clones per amplicon. As expected, two distinct sequence variants belonging to the 6A and 7A groups were obtained from each of these two samples (OP918911 and OP918913 for SK08-21, and OP918912 and OP918914 for SK11-21) (Figure 1). For both samples, the ratio between the 7A and 6A group variants was 2:1.
The identification of the HSVd-variant mixed infection in single grapevine plants prompted us to re-analyze the HTS datasets of the SK30, SK503, and SK933 samples, from which the 6A group variants were determined by default Geneious mapping (MN548397, OP918897, and OP918898). Interestingly, by remapping the datasets to represent the 6A and 7A sequence variants (SK11-21_clone1 and SK11-21_clone2, respectively) using strict parameters (Geneious mapper, identity 99%, overlap 20–30 bp), we found a mix of both types of groups in the SK503 and SK933 samples. For example, the coverage of the two variants identified in the SK503 dataset was 116.0x for the 6A group variant and 63.0x for the 7A group variant. Moreover, the visual inspection of mapping confirmed the presence of reads representing both group variants, containing the same seven group-dividing positions. These results confirmed and further pinpointed the complex and heterologous nature of viroid populations in some grapevine samples, already noted for grapevine viruses [39].
Because of the sole-6A-variant character of SK30 and the “prevailing” 6A group variant in the SK503 and SK933 samples, only those sequences were added to the final alignment. Moreover, due to the mutual shared nucleotide identity across the extensive parts of the HSVd genome (Figure 1A,B), the major quantity of 7A-like reads would have been unified and indistinguishably mixed with the identical 6A-like reads. The overall actual ratio between the 6A and 7A groups was thus very difficult to determine from the HTS datasets. Nonetheless, it would be interesting to also study the prevalence of these two groups in grapevine plants, especially in conjunction with the observations of whether one of the groups is more abundant (more copies of genomic RNA) in this host or if their numbers undergo some sort of oscillation across time.
Our molecular data fit with previously published results of molecular diversity in conjunction with different host adaptations performed in Japan [32]. In their study, natural HSVd-grapevine variations resulted in the original observation of the exact same seven mutated positions. Interestingly, in this study, the most abundant mutations were those belonging to the 6A group, contrary to the larger prevalence of the 7A group variants in Slovakia.
To further evaluate the significance of the 6A- and 7A-group-belonging mutations, we also analyzed all available grapevine HSVd genomes retrieved from the GenBank database ((n = 228), accessed on 1 October 2022). The robust alignment data again pinpointed the dichotomic nature of the 7A-group-specific mutations (i.e., occurring only as 6A- or 7A-specific ones, without random or uncorrelated mutations, Table 2). Similar to Kawaguchi-Ito et al. [32], the 6A group variants were found to be more abundant within the database. It should be noted, however, that some grapevine HSVd isolates displayed particular patterns of mutations in these parts of the genome, which did not support the clear-cut division [41] observed in the Slovak samples. This indicated the long-term diversification of the HSVd in grapevine hosts and the possible role of non-host selection pressures.

3.4. Group-Defining Mutations vs. RNA Secondary Structure

Previous sequence analyses led us to define the seven specific mutations located in the positions spanning the pathogenic and terminal left domains as the common divergence factor within the Slovak HSVd variants (Figure 1). Taking also into consideration the overall low molecular variability observed in the remaining parts of the HSVd genome in both groups, we analyzed the relationship between these mutations and predicted an RNA secondary structure. We based our approach on the necessity of the viroid molecule to maintain its secondary RNA structure in order to preserve its biological functionality [11,12,42].
Therefore, two representative genomes from each genetic group (SK805 (OP918903) from the 6A group and SK723 (OP918901) from the 7A group) were subjected to the secondary structure analysis.
As previously predicted for all members of the Pospiviroidae family [1,3,43], a typical rod-like shape secondary structure was obtained after folding [37] for both genome types (Figure 3). Nevertheless, the secondary structure analysis of the two HSVd variants provided a potential explanation for the seven-group-forming mutated positions. Hypothetically, those could have started as compensations for the nucleotide-base-length shift, probably triggered by the one specific adenosine insertion in positions 44/45 (Figure 1) due to the eucaryotic-DNA-dependent RNA II polymerase “slip” and base addition when replicating poly(N) sequences [44,45,46]. One can hypothesize that this could be the first candidate for the mutation having caused the splitting of these two sequence groups and triggered subsequent compensatory mutations in the HSVd-grapevine population in order to maintain the infectious functionality of their RNA secondary structures. The hypothesis of a six-adenosine stretch being the original variant is seemingly also supported by the previously mentioned study of HSVd-grapevine variants in Japan [32]. Finally, our samples SK30 (OP918897); SK503 (OP918898); SK933 (MN548397); SK803 (OP918902); SK805 (OP918903); and SK809 (OP918904) from the 6A group were also 100% identical to the AB219944 genome [31] obtained from 100-year-old grapevine with the corresponding length of 297 nucleotides.
Indeed, when observing the positions 256–257 (257–258 in the 7A group), we noted the exchange of order from U–G (6A group) to G–U (7A group), probably due to the necessity of base-pairing with the additional A in position 45 for the 7A group (Figure 3). The inserted A would have caused the next-in-line C (position 46 in SK723) to be opposite to U instead of G (positions 256 and 257, respectively, as it is in SK805). Such an exchange from U–G (256–257; 6A group) to G–U (257–258; 7A group), observable in our data, would be a very simple and effective way to deal with such a length-changing mutation. Additionally, the exchange of A–A in positions 46–47 (6A group) for U–U in positions 47–48 in the 7A group could also be an explanation for the maintenance of structural integrity. This exchange seemed to stabilize the loop consisting of two uracils in the 6A group, which was formed only on the lower strand. In the 7A group, this loop was formed on both strands and consisted of five uracils (Figure 3). Yet again, one adenosine insertion would have caused these two uracils (positions 253–254; 6A group) to be in opposition to the two adenosines (A–A in mutated positions 46–47; 6A group), which would have led to intra-molecular base-pairing resulting in potential loop closing. The presence of such loop motifs in the secondary structure of RNA is crucial for viroid replication and trafficking in host cells, described mainly in Potato spindle tuber viroid (PSTVd) species [11,14,47]. Interestingly, in the case of the PSTVd, the periphery of its pathogenic domain contains motifs responsible for supposed interactions with cellular proteins involved in viroid trafficking [11,47], which could also be the case in HSVd populations. This discussed loop is, although larger, still preserved in the 7A group due to the substitutions exchanging A–A to U–U (in positions 47–48). Therefore, we denoted all these four positions (plus the above-mentioned A insertion) as they contained the potentially structurally compensational mutations due to the one-adenosine insertion. Similar compensation was also observed in the case of the two unique insertions of the variant SK08-21_clone2 (OP918912) (additional C in position 31 was stabilized via base pairing with an additional G in positions 274/275) (Supplementary Figure S2), again supporting such a type of mechanism for viroid RNA adaptability when dealing with random nucleotide insertions.
The remaining two group-specific mutated positions were apparently without any structural consequences, as both loops they were mapped into had in both groups the same structural integrity (Figure 3). They included an exchange of C (6A group) to A (7A group) in position 26 and the exchange of A (6A group) to U (7A group) in position 32. Because of the lack of structural interpretation, they could potentially have some functionally important compensatory meaning. For example, in the PSTVd, the very beginning of the RNA sequence (the terminal left domain) is crucial for replication, as it contains just those motifs recognized by RNA Pol II as well as by the unique splicing form of transcription factor IIIA (TFIIIA-7ZF) [11,13,14,47]. Experimental data regarding HSVd infectivity and, more importantly, symptom severity in cucumber (Cucumis sativus L.) [48], provided us with another potential clue for their interpretation.
Moreover, Xia et al. [48] used two HSVd variants (hop and grapevine) to study the involvement of mutations to symptomatology in cucumber. Interestingly, from 4 mutually differing positions, one mutation was located at position 26. The grapevine variant, causing severe symptoms in cucumber, had C in this position, similar to the 6A group. On the other hand, the hop variant, causing only mild symptoms, had A in position 26, similar to the 7A group.

4. Conclusions

The importance of studying the exact points and changes in the nucleotide sequence of the viroids for epidemiology purposes was supported by previous ecological as well as experimental studies. Those range from new host adaptations (e.g., the PSTVd on Nicotiana tabacum via a single-nucleotide substitution in the lower strand of the CCR [49]) to symptom expression (e.g., in the case of the HSVd, only a few, even if single nucleotide changes in the variable domain, can alter the severity of citrus cachexia symptoms [50]). This simple mutation dependence reflected two intertwining molecular principles: (i) the functional simplicity of viroids (RNA sequence) interacting with functionally similar cellular proteins of (ii) slightly different substrate specificities in different plant hosts due to their mutual structural and, more importantly, conformational (subtle) changes [1].
In this paper, we report the analysis of the complete genomic sequences of 19 HSVd-grapevine variants, which led to the finding of two specific-mutation-dependent molecular groups spread throughout Slovakia. It is noteworthy that these same specific mutations were previously reported by Kawaguchi-Ito et al. [32]. This fact ought to be explained by two factors: first, the antiquity of the HSVd pathogen as an infectious RNA molecule and its necessity to maintain its secondary RNA structure. To preserve its biological functionality, only a small number of mutations can be tolerated, and, at the same time, only in the specific locations of the RNA molecule. The second is the vegetative propagation of V. vinifera cultivars around the world, mainly for agricultural and breeding purposes. Moreover, grapevine, as the hypothesized original source of the HSVd [31], can create the specific host intracellular evolutionary pressure via its specific viroid-exploited enzymes.
Although the importance of group-specific mutations highlighted in this work was strongly supported by in silico data, their biological significance and correlation to HSVd infectivity and/or symptom severity must be experimentally tested. Such an understanding of the diversity of grapevine HSVd populations would have a significant role in the epidemiology of the viroid and constitutes an important prerequisite for the development of accurate diagnostics and efficient control strategies.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/pathogens12020205/s1, Figure S1: chromatogram of the PCR amplicon of the SK11-21 sample obtained by Sanger sequencing indicative of mixed infection by different HSVd variants; Figure S2: secondary structure of the SK08-21_clone2 (6A group).

Author Contributions

Conceptualization, P.A. and M.G.; methodology, A.A., L.P. and M.G.; investigation, P.A., A.A., Z.Š., M.M. and M.G.; resources, Z.Š., M.M. and M.G.; data curation, P.A. and M.G.; writing—original draft preparation, P.A.; writing—review and editing, P.A., L.P., A.A., Z.Š. and M.G.; funding acquisition, M.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by grant VEGA 2/0030/20 from the Scientific Grant Agency of the Ministry of Education and Slovak Academy of Sciences and by the European Union’s Horizon 2020 Research and Innovation Program (EVA-GLOBAL project, grant agreement number 871029). P.A. and M.G. acknowledge support from the Operational Programme Integrated Infrastructure (313011ASN4), co-financed by the European Regional Development Fund. This publication was created thanks to support under the Operational Programme Integrated Infrastructure for the project: Strengthening of Research, Development, and Innovation Capacities of Translational Biomedical Research of Human Diseases, IMTS: 313021BZC9, co-financed by the European Regional Development Fund.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The nucleotide sequences reported in this paper are deposited into the GenBank database (www.ncbi.nlm.nih.gov) under the accession numbers listed in the text.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Flores, R.; Daròs, J.-A.; Hernández, C.; Navarro, B.; Di Serio, F. Viroids. eLS 2020, 1, 192–203. [Google Scholar] [CrossRef]
  2. Diener, T.O. Potato spindle tuber ‘‘virus’’. IV. A replicating, low-molecular weight RNA. Virology 1971, 45, 411–428. [Google Scholar] [CrossRef] [PubMed]
  3. Diener, T.O. Discovering viroids—A personal perspective. Nat. Rev. Microbiol. 2003, 1, 75–79. [Google Scholar] [CrossRef] [PubMed]
  4. Katsarou, K.; Adkar-Purushothama, C.R.; Tassios, E.; Samiotaki, M.; Andronis, C.; Lisón, P.; Nikolaou, C.; Perreault, J.-P.; Kalantidis, K. Revisiting the Non-Coding Nature of Pospiviroids. Cells 2022, 11, 265. [Google Scholar] [CrossRef]
  5. Hammond, R.W.; Owens, R.A. Viroids: New and Continuing Risks for Horticultural and Agricultural Crops; APS: St. Paul. MN, USA, 2006. [Google Scholar] [CrossRef]
  6. Silletti, M.R.; Navarro, B.; Bozzano, G.; Trisciuzzi, V.N.; Di Serio, F. The PSTVd viroid, a threat to tomatoes and potatoes. Informatore Agrario 2009, 65, 89–90. [Google Scholar]
  7. Hadidi, A.; Vidalakis, G.; Sano, T. Chapter 2—Economic Significance of Fruit Tree and Grapevine Viroids. In Viroids and Satellites; Hadidi, A., Flores, R., Randles, J.W., Palukaitis, P., Eds.; Academic Press-Elsevier: Cambridge, MA, USA, 2017; pp. 15–25. ISBN 978-0-12-801498-1. [Google Scholar]
  8. Hammond, R.W. Economic significance of viroids in vegetable and field crops. In Viroids and Satellites; Hadidi, A., Flores, R., Randles, J.W., Palukaitis, P., Eds.; Academic Press-Elsevier: Cambridge, MA, USA, 2017; pp. 5–13. ISBN 978-0-12-801498-1. [Google Scholar]
  9. Mackie, A.E.; Barbetti, M.J.; Rodoni, B.; McKirdy, S.J.; Jones, R.A.C. Effects of a Potato Spindle Tuber Viroid Tomato Strain on the Symptoms, Biomass, and Yields of Classical Indicator and Currently Grown Potato and Tomato Cultivars. Plant Dis. 2019, 103, 3009–3017. [Google Scholar] [CrossRef]
  10. Takeda, R.; Petrov, A.I.; Leontis, N.B.; Ding, B. A Three-Dimensional RNA Motif in Potato spindle tuber viroid Mediates Trafficking from Palisade Mesophyll to Spongy Mesophyll in Nicotiana benthamiana. Plant Cell 2011, 23, 258–272. [Google Scholar] [CrossRef] [Green Version]
  11. Zhong, X.; Archual, A.J.; Amin, A.A.; Ding, B. A Genomic Map of Viroid RNA Motifs Critical for Replication and Systemic Trafficking. Plant Cell 2008, 20, 35–47. [Google Scholar] [CrossRef] [Green Version]
  12. Flores, R.; Gas, M.-E.; Molina-Serrano, D.; Nohales, M.-Á.; Carbonell, A.; Gago, S.; De la Peña, M.; Daròs, J.-A. Viroid Replication: Rolling-Circles, Enzymes and Ribozymes. Viruses 2009, 1, 317–334. [Google Scholar] [CrossRef] [Green Version]
  13. Wang, Y.; Qu, J.; Ji, S.; Wallace, A.J.; Wu, J.; Li, Y.; Gopalan, V.; Ding, B. A Land Plant-Specific Transcription Factor Directly Enhances Transcription of a Pathogenic Noncoding RNA Template by DNA-Dependent RNA Polymerase II. Plant Cell 2016, 28, 1094–1107. [Google Scholar] [CrossRef] [PubMed]
  14. Mudiyanselage, S.D.D.; Qu, J.; Tian, N.; Jiang, J.; Wang, Y. Potato Spindle Tuber Viroid RNA-Templated Transcription: Factors and Regulation. Viruses 2018, 10, 503. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Flores, R.; Daròs, J.A.; Hernández, C. Avsunviroidae family: Viroids containing hammerhead ribozymes. Adv. Virus Res. 2000, 55, 271–323. [Google Scholar] [CrossRef] [PubMed]
  16. Yamamoto, H.; Kagami, Y.; Kurokawa, M.; Nishimura, S.; Ukawa, S.; Kubo, S. Studies on hop stunt disease in Japan. Rep. Res. Lab. Kirin Brew. 1973, 16, 49–62. [Google Scholar]
  17. Sasaki, M.; Shikata, E. On Some Properties of Hop Stunt Disease Agent, a Viroid. Proc. Jpn. Acad. Ser. B 1977, 53, 109–112. [Google Scholar] [CrossRef] [Green Version]
  18. van Dorst, H.J.M.; Peters, D. Some biological observations on pale fruit, a viroid-incited disease of cucumber. Neth. J. Plant Pathol. 1974, 80, 85–96. [Google Scholar] [CrossRef]
  19. Sano, T.; Uyeda, I.; Shikata, E.; Ohno, T.; Okada, Y. Nucleotide sequence of cucumber pale fruit viroid: Homology to hop stunt viroid. Nucleic Acids Res. 1984, 12, 3427–3434. [Google Scholar] [CrossRef] [Green Version]
  20. Semancik, J.S.; Roistacher, C.N.; Rivera-Bustamante, R.; Duran-Vila, N. Citrus Cachexia Viroid, a New Viroid of Citrus: Relationship to Viroids of the Exocortis Disease Complex. J. Gen. Virol. 1988, 69, 3059–3068. [Google Scholar] [CrossRef]
  21. Reanwarakorn, K.; Semancik, J.S. Correlation of Hop Stunt Viroid Variants to Cachexia and Xyloporosis Disease of Citrus. Phytopathology 1999, 89, 568–574. [Google Scholar] [CrossRef] [Green Version]
  22. Hataya, T.; Tsushima, T.; Sano, T. Hop stunt viroid. In Viroids and Satellites; Hadidi, A., Flores, R., Randles, J., Palukaitis, P., Eds.; Academic Press-Elsevier: Cambridge, MA, USA, 2017; pp. 199–210. ISBN 978-0-12-801498-1. [Google Scholar]
  23. Vamenani, R.; Rahimian, H.; Alavi, S.M.; Parizi, A.P.; Razzaz, T.M. Genetic diversity of hops stunt viroid from symptomatic and asymptomatic citrus trees in Iran. J. Phytopathol. 2019, 167, 484–489. [Google Scholar] [CrossRef]
  24. Sano, T.; Hataya, T.; Terai, Y.; Shikata, E. Hop stunt viroid strains from dapple fruit disease of plum and peach in Japan. J. Gen. Virol. 1989, 70, 1311–1319. [Google Scholar] [CrossRef]
  25. Astruc, N.; Marcos, J.F.; Macquaire, G.; Candresse, T.; Pallás, V. Studies on the diagnosis of hop stunt viroid in fruit trees: Identification of new hosts and application of a nucleic acid extraction procedure based on non-organic solvents. Eur. J. Plant Pathol. 1996, 102, 837–846. [Google Scholar] [CrossRef]
  26. Xu, L.; Wang, J.W.; Zhu, D.Z.; Zong, X.J.; Wei, H.R.; Chen, X.; Hammond, R.W.; Liu, Q.Z. First Report of Hop stunt viroid From Sweet Cherry With Dapple Fruit Symptoms in China. Plant Dis. 2017, 101, 394. [Google Scholar] [CrossRef]
  27. Marquez-Molins, J.; Gomez, G.; Pallas, V. Hop stunt viroid: A polyphagous pathogenic RNA that has shed light on viroid-host interactions. Mol. Plant Pathol. 2021, 22, 153–162. [Google Scholar] [CrossRef]
  28. Lin, C.-Y.; Wu, M.-L.; Shen, T.-L.; Yeh, H.-H.; Hung, T.-H. Multiplex detection, distribution, and genetic diversity of Hop stunt viroid and Citrus exocortis viroid infecting citrus in Taiwan. Virology J. 2015, 12, 11. [Google Scholar] [CrossRef] [Green Version]
  29. Di Serio, F.; Ambrós, S.; Sano, T.; Flores, R.; Navarro, B. Viroid Diseases in Pome and Stone Fruit Trees and Koch’s Postulates: A Critical Assessment. Viruses 2018, 10, 612. [Google Scholar] [CrossRef] [Green Version]
  30. Maddahian, M.; Massumi, H.; Heydarnejad, J.; Hosseinipour, A.; Khezri, A.; Sano, T. Biological and molecular characterization of hop stunt viroid variants from pistachio trees in Iran. J. Phytopathol. 2019, 167, 163–173. [Google Scholar] [CrossRef]
  31. Sano, T.; Mimura, R.; Ohshima, K. Phylogenetic analysis of hop and grapevine isolates of hop stunt viroid supports a grapevine origin for hop stunt disease. Virus Genes 2001, 22, 53–59. [Google Scholar] [CrossRef]
  32. Kawaguchi-Ito, Y.; Li, S.-F.; Tagawa, M.; Araki, H.; Goshono, M.; Yamamoto, S.; Tanaka, M.; Narita, M.; Tanaka, K.; Liu, S.-X.; et al. Cultivated Grapevines Represent a Symptomless Reservoir for the Transmission of Hop Stunt Viroid to Hop Crops: 15 Years of Evolutionary Analysis. PLoS ONE 2009, 4, e8386. [Google Scholar] [CrossRef]
  33. Vozárová, Z.; Sihelská, N.; Predajňa, L.; Šoltys, K.; Glasa, M. First report of Grapevine yellow speckle viroid-1 infecting grapevine in Slovakia. J. Plant Pathol. 2016, 98, 697. [Google Scholar] [CrossRef]
  34. Glasa, M.; Predajňa, L.; Šoltys, K.; Sabanadzovic, S.; Olmos, A. Detection and molecular characterisation of Grapevine Syrah virus-1 isolates from Central Europe. Virus Genes 2015, 51, 112–121. [Google Scholar] [CrossRef]
  35. Glasa, M.; Predajňa, L.; Sihelská, N.; Šoltys, K.; Ruiz-García, A.B.; Olmos, A.; Wetzel, T.; Sabanadzovic, S. Grapevine virus T is relatively widespread in Slovakia and Czech Republic and genetically diverse. Virus Genes 2018, 54, 737–741. [Google Scholar] [CrossRef] [PubMed]
  36. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [Green Version]
  37. Zuker, M. Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31, 3406–3415. [Google Scholar] [CrossRef]
  38. Glasa, M.; Predajňa, L.; Komínek, P.; Nagyová, A.; Candresse, T.; Olmos, A. Molecular characterization of divergent grapevine Pinot gris virus isolates and their detection in Slovak and Czech grapevines. Arch. Virol. 2014, 159, 2103–2107. [Google Scholar] [CrossRef] [PubMed]
  39. Glasa, M.; Predajňa, L.; Šoltys, K.; Sihelská, N.; Nagyová, A.; Wetzel, T.; Sabanadzovic, S. Analysis of Grapevine rupestris stem pitting-associated virus (GRSPaV) in Slovakia reveals differences in intra-host population diversity and naturally occurring recombination events. Plant Pathol. J. 2017, 33, 34–42. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Glasa, M.; Predajňa, L.; Sihelská, N.; Šoltys, K.; Ruiz-García, A.B. Analysis of Virome by High-Throughput Sequencing Revealed Multiple Infection and Intra-Virus Diversity in a Single Grapevine Plant. Acta Hortic. Regiotect. 2020, 23, 35–39. [Google Scholar] [CrossRef]
  41. Vončina, D.; Almeida, R.P.P. Screening of some Croatian autochthonous grapevine varieties reveals a multitude of viruses, including novel ones. Arch. Virol. 2018, 163, 2239–2243. [Google Scholar] [CrossRef] [Green Version]
  42. Meshi, T.; Ishikawa, M.; Watanabe, Y.; Yamaya, J.; Okada, Y.; Sano, T.; Shikata, E. The sequence necessary for the infectivity of hop stunt viroid cDNA clones. Mol. Gen. Genet. 1985, 200, 199–206. [Google Scholar] [CrossRef]
  43. López-Carrasco, A.; Flores, R. Dissecting the secondary structure of the circular RNA of a nuclear viroid in vivo: A “naked” rod-like conformation similar but not identical to that observed in vitro. RNA Biol. 2017, 14, 1046–1054. [Google Scholar] [CrossRef] [Green Version]
  44. Mühlbach, H.-P.; Sänger, H.L. Viroid replication is inhibited by α-amanitin. Nature 1979, 278, 185–188. [Google Scholar] [CrossRef]
  45. Rackwitz, H.R.; Rohde, W.; Sänger, H.L. DNA-dependent RNA polymerase II of plant origin transcribes viroid RNA into full-length copies. Nature 1981, 291, 297–301. [Google Scholar] [CrossRef] [PubMed]
  46. Schindler, I.M.; Mühlbach, H.P. Involvement of nuclear DNA-dependent RNA polymerases in potato spindle tuber viroid replication: A reevaluation. Plant Sci. 1992, 84, 221–229. [Google Scholar] [CrossRef]
  47. Jiang, J.; Smith, H.N.; Ren, D.; Mudiyanselage, S.D.D.; Dawe, A.L.; Wang, L.; Wang, Y. Potato Spindle Tuber Viroid Modulates Its Replication through a Direct Interaction with a Splicing Regulator. J. Virol. 2018, 92, e01004-18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Xia, C.; Li, S.; Hou, W.; Fan, Z.; Xiao, H.; Lu, M.; Sano, T.; Zhang, Z. Global Transcriptomic Changes Induced by Infection of Cucumber (Cucumis sativus L.) with Mild and Severe Variants of Hop Stunt Viroid. Front. Microbiol. 2017, 8, 2427. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Wassenegger, M.; Spieker, R.L.; Thalmeir, S.; Gast, F.U.; Riedel, L.; Sänger, H.L. A single nucleotide substitution converts potato spindle tuber viroid (PSTVd) from a noninfectious to an infectious RNA for nicotiana tabacum. Virology 1996, 226, 191–197. [Google Scholar] [CrossRef] [PubMed]
  50. Serra, P.; Gago, S.; Duran-Vila, N. A single nucleotide change in Hop stunt viroid modulates citrus cachexia symptoms. Virus Res. 2008, 138, 130–134. [Google Scholar] [CrossRef]
Pathogens 12 00205 g001 550
Figure 1. The alignment of 19 whole-viroid sequences of HSVds from 17 tested grapevines from Slovakia. (A,C) the positions of the crucial and cluster-dividing mutations spread across the pathogenicity and terminal left domains (yellow arrows and rectangles with red outliners) with the one (red arrow) possibly initiating the emergence of all of the others. Other colored bases represent the positions of only one-variant-specific single-nucleotide mutations. (A,B) the region (representing the central, variability, and terminal right domains) that was uniformly conserved among all 19 Slovak grapevine isolates (from positions 100 to 200, approximately). Isolates belonging to the 6A cluster are in blue rectangles and isolates belonging to the 7A cluster are in yellow rectangles. The three sequences determined from the HTS data are marked by asterisks. The sequence of the forward primer HSV-78P [31] is marked by grey rectangles.
Figure 1. The alignment of 19 whole-viroid sequences of HSVds from 17 tested grapevines from Slovakia. (A,C) the positions of the crucial and cluster-dividing mutations spread across the pathogenicity and terminal left domains (yellow arrows and rectangles with red outliners) with the one (red arrow) possibly initiating the emergence of all of the others. Other colored bases represent the positions of only one-variant-specific single-nucleotide mutations. (A,B) the region (representing the central, variability, and terminal right domains) that was uniformly conserved among all 19 Slovak grapevine isolates (from positions 100 to 200, approximately). Isolates belonging to the 6A cluster are in blue rectangles and isolates belonging to the 7A cluster are in yellow rectangles. The three sequences determined from the HTS data are marked by asterisks. The sequence of the forward primer HSV-78P [31] is marked by grey rectangles.
Pathogens 12 00205 g001
Pathogens 12 00205 g002 550
Figure 2. Neighbor-joining phylogenetic tree constructed in MEGA7 from alignments of 19 complete HSVd sequences generated in this study and 23 database sequences. Slovak HSVd isolates are in blue or yellow rectangles, depending on their grouping. The database sequences are identified by their accession numbers, host, and geographical location. Only bootstrap values ≥ 70% (1000 bootstrap re-samplings) are indicated. The scale bar indicates a genetic distance of 0.01.
Figure 2. Neighbor-joining phylogenetic tree constructed in MEGA7 from alignments of 19 complete HSVd sequences generated in this study and 23 database sequences. Slovak HSVd isolates are in blue or yellow rectangles, depending on their grouping. The database sequences are identified by their accession numbers, host, and geographical location. Only bootstrap values ≥ 70% (1000 bootstrap re-samplings) are indicated. The scale bar indicates a genetic distance of 0.01.
Pathogens 12 00205 g002
Pathogens 12 00205 g003 550
Figure 3. The secondary structure of two representative complete HSVd sequences, i.e., SK805 (6A group, blue color) and SK723 (7A group, yellow color). The Gibbs free energy for the most thermodynamically stable conformation, calculated via Mfold software [37], was −114.0 and −114.10, respectively. The central conserved region (CCR), into which the diagnostic primers were targeted, is marked (red rectangle). The speculated compensational mutations were spread only across the terminal left and pathogenic domains. Each mutated base in both isolates are marked by arrows (purple for structurally (also indicated by double red arrow) and pink for functionally considered mutations).
Figure 3. The secondary structure of two representative complete HSVd sequences, i.e., SK805 (6A group, blue color) and SK723 (7A group, yellow color). The Gibbs free energy for the most thermodynamically stable conformation, calculated via Mfold software [37], was −114.0 and −114.10, respectively. The central conserved region (CCR), into which the diagnostic primers were targeted, is marked (red rectangle). The speculated compensational mutations were spread only across the terminal left and pathogenic domains. Each mutated base in both isolates are marked by arrows (purple for structurally (also indicated by double red arrow) and pink for functionally considered mutations).
Pathogens 12 00205 g003
Table
Table 1. List of Slovak HSVd isolates with some characteristics.
Table 1. List of Slovak HSVd isolates with some characteristics.
SampleHostLocalityYear of SamplingGenBank Accession Number
SK30 *V. vinifera, cv. VeltlinerPezinok2012OP918897
SK503 *V. vinifera, cv. DornfelderŠenkvice2015OP918898
SK933 *V. vinifera, unknown white berriesPezinok2017MN548397
SK717V. vinifera, unknown white berriesPezinok2016OP918899
SK721V. vinifera, unknown red berriesPezinok2016OP918900
SK723V. vinifera, cv. Muller-ThurgauPezinok2016OP918901
SK803V. vinifera, cv. PalavaPezinok2017OP918902
SK805V. vinifera, cv. VeltlinerPezinok2017OP918903
SK809V. vinifera, cv. WelshrieslingPezinok2017OP918904
SK928V. vinifera, cv. SilvanerLimbach2017OP918905
SK934V. vinifera, cv. DornfelderCífer2017OP918906
SK02-21V. vinifera, unknown red berriesPezinok2021OP918907
SK03-21V. vinifera, unknown red berriesPezinok2021OP918908
SK04-21V. vinifera, cv. TraminerPezinok2021OP918909
SK07-21V. vinifera, unknown white berriesPezinok2021OP918910
SK08-21V. vinifera, cv. VeltlinerPezinok2021OP918911 (clone 1)
OP918912 (clone 2)
SK11-21V. vinifera, cv. Muller-ThurgauPezinok2021OP918913 (clone 1)
OP918914 (clone 2)
* Whole viroid sequences generated from the HTS data.
Table
Table 2. Frequency of the occurrence of seven group-defining mutations in complete HSVd sequences from grapevine available in GenBank. The total number of analyzed sequences was 228; however, some contain degenerated bases in the analyzed positions; thus, the overall sum of the counted positive sequences is >228 in each row (for clarity, the data are also converted to percentage).
Table 2. Frequency of the occurrence of seven group-defining mutations in complete HSVd sequences from grapevine available in GenBank. The total number of analyzed sequences was 228; however, some contain degenerated bases in the analyzed positions; thus, the overall sum of the counted positive sequences is >228 in each row (for clarity, the data are also converted to percentage).
Nucleotide Position *6A Group
Nucleotide/Number of Sequences/% 		7A Group
Nucleotide/Number of Sequences/%
26C/151/62.66A/90/37.34
32A/156/64.20U/85/35.00 **
44/45Δ/158/65.56A/83/34.44
46/47A/162/66.94U/80/33.06
47/48A/185/77.40U/53/22.18 ***
256/257U/170/70.83G/70/29.17
257/258G/166/69.17U/74/30.83
* Listed according to SK805_OP918903/SK723_OP918901; ** 2 out of 228 references (0.80%) had a different base in this position (G); *** 1 out of 228 references (0.42%) had a different base in this position (C); Δ: nucleotide deletion.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Alaxin, P.; Predajňa, L.; Achs, A.; Šubr, Z.; Mrkvová, M.; Glasa, M. Analysis of Hop Stunt Viroid Diversity in Grapevine (Vitis vinifera L.) in Slovakia: Coexistence of Two Particular Genetic Groups. Pathogens 2023, 12, 205. https://doi.org/10.3390/pathogens12020205

AMA Style

Alaxin P, Predajňa L, Achs A, Šubr Z, Mrkvová M, Glasa M. Analysis of Hop Stunt Viroid Diversity in Grapevine (Vitis vinifera L.) in Slovakia: Coexistence of Two Particular Genetic Groups. Pathogens. 2023; 12(2):205. https://doi.org/10.3390/pathogens12020205

Chicago/Turabian Style

Alaxin, Peter, Lukáš Predajňa, Adam Achs, Zdeno Šubr, Michaela Mrkvová, and Miroslav Glasa. 2023. "Analysis of Hop Stunt Viroid Diversity in Grapevine (Vitis vinifera L.) in Slovakia: Coexistence of Two Particular Genetic Groups" Pathogens 12, no. 2: 205. https://doi.org/10.3390/pathogens12020205

APA Style

Alaxin, P., Predajňa, L., Achs, A., Šubr, Z., Mrkvová, M., & Glasa, M. (2023). Analysis of Hop Stunt Viroid Diversity in Grapevine (Vitis vinifera L.) in Slovakia: Coexistence of Two Particular Genetic Groups. Pathogens, 12(2), 205. https://doi.org/10.3390/pathogens12020205

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop