Whole Genome Characterization of Prunus Necrotic Ringspot Virus and Prune Dwarf Virus Infecting Stone Fruits in Russia

: We conducted a survey of the phytosanitary status of the Prunus germplasm collection in the Nikita Botanical Gardens, Yalta, Russia. The virome of plants displaying virus-like symptoms was studied using Illumina MiSeq high-throughput sequencing. Reads related to prunus necrotic ringspot virus (PNRSV), prune dwarf virus (PDV), and ourmia-like virus 1 (OuLV1) were generated in a number of samples. Near complete genomes of two divergent PNRSV isolates, PDV isolate, and a contig partly covered OuLV1 genome were assembled de novo using the metaSPAdes program. The structure of the genomic RNA1, RNA2, and RNA3 of the new ilarvirus isolates was shown to be typical of PNRSV and PDV. This is the ﬁrst report and characterization of the PNRSV and PDV full-length genomes from Russia, expanding the information on their geographical distribution and genetic diversity. An open reading frames (ORF)-based phylogeny of all full-length PNRSV and PDV genomes available in GenBank divided each ORF into two or three main clusters. A number of isolates migrated from one cluster to another cluster, depending on the analyzed genome segment. The different branching order may indicate reassortment in the evolutionary history of some PDV and PNRSV isolates.

PNRSV and PDV belong to the genus Ilarvirus from the family Bromoviridae.Their segmented genome consists of three single-stranded positive-sense RNA molecules, each of which is packed into separate particles.The 5 -ends of genome RNAs are capped and the 3 -terminal sequences are organized into an array of stem-loops interspersed with nonpaired sequences.Five open reading frames (ORFs) encode replicase (ORF1), RNAdependent RNA polymerase (RdRp, ORF2a), 2b protein (ORF2b), movement protein (MP, ORF3a), and coat protein (CP, ORF3b).However, RNA2 of PNRSV and PDV are monocistronic and encode only RdRp.Replicase and RdRp ensure replication of the viral genome.MP provides cell-to-cell transport of the virus, and CP is responsible for encapsidation of the viral genome, its spread across the whole plant, and activation of the infection by binding to the 3 -ends of the genome RNAs [8].
Among plants, PNRSV and PDV are transmitted through seeds, pollen, and vegetative propagation.Insect pollinators play an important role in the transmission of these viruses [9].In nature, PDV infects only Prunus spp.and is detected mainly in sweet and sour cherries, while PNRSV has a wider host range, including not only Prunus but also hop, apple, and rose.Symptoms are dependent on the virus/host combination, but chlorotic ringspots and line patterns are common on the leaves of PNRSV-and PDV-infected plants.In addition, PNRSV may cause shot holes in leaves and PDV is often associated with sour cherry yellows disease.Both viruses, especially in mixed infections, can be highly pathogenic for stone fruits and able to cause significant yield losses (20% to 60%, depending on the host) by reducing fruit growth, delaying fruit maturity, and making fruits unmarketable.They can also decrease pollen germination and bud-take after grafting (40% to 50%) and lead to stunting of infected peach and plum trees ( [10,11] and the references therein).
PNRSV and PDV occur worldwide.Their global distribution seems to stem from a wide international exchange of infected plant material, since the potential damage caused by these viruses has long been underestimated.Although movement of fruit tree cultivars is regulated, some plants apparently have not been virus-indexed or were symptomless but infected [11].On the other hand, PNRSV and PDV are quite well studied.In particular, hundreds of partial and complete sequences of their genomes from different hosts and regions of the world are deposited in GenBank, enabling the evaluation of the genetic diversity of these viruses and the improvement of measures for their control.
In Russia, these two ilarviruses were detected on cultivated Prunus species by ELISA in the 2000s [12,13].However, no genome of the Russian PDV or PNRSV isolates has been sequenced or characterized.The detection and molecular characterization of new Russian PNRSV and PDV isolates would obviously extend the information on the geographical distribution and genetic diversity of these viruses.
One of the largest germplasm collections of stone fruits in Russia is located in the Nikita Botanical Gardens (NBG), Yalta.Over three thousand various Prunus species, cultivars, hybrids, and breeding forms, both local and introduced from North America, Southern Europe, and Central Asia, are maintained in the collection as a repository of genetic diversity.This gene pool is widely used in breeding and biotechnology programs.Many Prunus genotypes originating from this collection are sent in hybrid and cultivar test plots to other regions of Russia [14].Meanwhile, a visual examination of the collection revealed trees displaying virus-like symptoms on the leaves.In this connection, we conducted a survey of the phytosanitary status of the NBG Prunus germplasm collection to identify assumed viruses and prevent their further spread to nighboring plants and new regions.The virome of some symptomatic trees was studied using high-throughput sequencing (HTS).Reads related to the PNRSV and PDV genomes were generated in a number of samples.
The objectives of this work were (i) confirming PNRSV and PDV in these samples using enzyme-linked immunosorbent assay (ELISA) and reverse transcription polymerase chain reaction (RT-PCR), (ii) sequencing, assembling, and characterizing full-length genomes of Russian PNRSV and PDV isolates from HTS-generated reads, (iii) evaluating the genetic diversity of these viruses based on the phylogeny of full-length genome sequences available in GenBank, and (iv) comparing the Russian PNRSV and PDV isolates with those from other sources.The results of this work expand information on the geographical distribution and genetic diversity of these viruses and can be implemented in breeding programs to improve virus control.

Sampling
Trees displaying virus-like symptoms were found in the Prunus germplasm collection of the NBG situated in the steppe zone of the Crimean peninsula (N45 • 06 50 ; E34 • 00 50 ) in June of 2020.Attention was mainly paid to the most prospective hybrids and cultivars.Individual samples composed of five to eight symptomatic leaves were taken from each selected tree.The samples were bagged, labeled, and delivered to the Virology Department of Lomonosov Moscow State University, where they were stored at 4 • C until use.

Sequence Analyses
To analyze the whole genomes of the Russian PNRSV and PDV isolates, the available sequences of these viruses were retrieved from GenBank.Multiple nucleotide (nt) sequences alignments using ClustalW, phylogenetic analysis, and visualization of trees were performed in MEGA7 [18].Phylogenetic trees were reconstructed using the neighborjoining method, the Kimura 2-parameter model, and bootstrap testing from 1000 replicates.Codon positions included were 1st + 2nd + 3rd + noncoding.All positions containing gaps and missing data were eliminated.Sequence identities were calculated using BLASTn (https://blast.ncbi.nlm.nih.gov/Blast.cgi,accessed on 19 September 2021).ORFs in the complete virus genomes were identified using an ORF finder (https://ncbi.nlm.nih.gov/orffinder, accessed on 19 September 2021) with default parameters.Conserved domains in virus proteins were mapped using the Conserved Domain Database (CDD, https://ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi, accessed on 19 September 2021).The Recombination Detection Program (RDP4.101)[19] was used to search recombination events in the genomes of PNRSV and PDV isolates using the default setting, except that the options "sequences are linear" and "list events detected by >4 methods" were chosen.

Enzyme-Linked Immunosorbent Assay (ELISA)
The samples were analyzed by double antibody sandwich (DAS) ELISA using LOEWE Biochemica GmbH (Sauerlach, Germany) reagent sets 07051 and 07052 for the detection of PDV and PNRSV, respectively.The leaf tissue was ground in phosphate-buffered saline (PBS) supplemented with 0.02% (v/v) Tween 20, 2% (w/v) polyvinylpyrrolidone (MW of about 40,000), 0.5% (v/v) Triton X-100, and 0.02% (w/v) sodium azide using a 1/20 ratio (w/v).The extracts were clarified by low-speed centrifugation on a MiniSpin centrifuge (Eppendorf, Hamburg, Germany) and placed in duplicate in MaxiSorp microplate wells (Nunc, New York, NY, USA), pre-coated with the polyclonal virus-specific antibodies, and incubated for 2 h at 37 • C. The subsequent steps were performed according to the manufacturer's instructions.The positive and negative controls were from the mentioned kits.The optical densities were measured 40 min after substrate addition at the wavelength of 405 nm.

Plant Material and Virus Detection
Three ten-year-old trees with symptoms of suspected viral infection were selected for metatranscriptomic analysis.The samples taken from these trees were as follows: P10 from P. domestica cultivar " Čačanska beste" grafted on P. cerasifera; P11 from P. domestica seedling of unknown cultivar; and P13 from ((P. cerasifera Ehrh.x P. armeniaca L.) x P. brigantiaca Vill.) hybrid grafted on P. cerasifera.The symptoms observed on the leaves of these plants are presented in Figure 1.

HTS Results
The reads related to PNRSV and PDV were generated from samples P10, P11, and P13 (Table 1).They covered the assembled ilarvirus-specific contigs with an average depth of coverage 145× to 900×.In addition, reads related to the fungal ourmia-like virus 1 (OuLV1) genome were found in sample P10.

Number of Clean Reads a
Number of Virus- Both P10 and P13, as well as the P11 samples, tested positive in ELISA for PNRSV and PDV, respectively (Table S1).The amplicons of the expected sizes, 380 bps, were obtained by RT-PCR assay (Figure S1).Their sequencing by the Sanger method confirmed PNRSV in samples P10 and P13 and PDV in sample P11.Thus, the HTS results were corroborated by ELISA and RT-PCR.

HTS Results
The reads related to PNRSV and PDV were generated from samples P10, P11, and P13 (Table 1).They covered the assembled ilarvirus-specific contigs with an average depth of coverage 145× to 900×.In addition, reads related to the fungal ourmia-like virus 1 (OuLV1) genome were found in sample P10.

Characterization of the PDV Genome
Three contigs related to PDV genomic RNAs were assembled in the P11 sample.According to BLASTn, the contig of 3341 nt was most closely related (95.2-95.5%) to RNA1 of Canadian PDV isolates 13C233 (MZ220979), 13C259 (MZ220983), and Salmo BC cherry (U57648) from sweet cherry and Australian isolate NS9R1 (KY883326) from Prunus.This contig covered genomic RNA1 of the related isolates nearly completely, except for seven to 26 3 -terminal and eight 5 -terminal nt.ORF1 3168 nt long encoding replicase of 1055 amino acid (aa) residues was identified in the contig.The methyltransferase (MET) and helicase (HEL) domains were predicted at positions 95-394 and 766-1025, respectively.ORF1 was flanked with 5 -and 3 -non-coding regions (NCR) of 29 and 144 nt, respectively.
Another contig, 2585 nt long, was most closely related (92.1-92.3%identity) to RNA2 of Canadian isolates Niagara D5 (MK522388), 13C226 (MZ220866), and 13C233 (MZ220870) from sweet cherry, encompassing RNA2 of these isolates nearly completely, except for four to 49 nt at their 5 or 3 ends.ORF2 2364 nt long encoding a protein of 787 aa was identified in the contig.Bromoviridae RdRp motif was found at positions 298-669 of this protein.ORF2 was flanked with 5 -and 3 -NCRs of 17 and 204 nt, respectively.
One more contig, of 2124 nt, was most closely related (95-95.9%identity) to RNA3 of PDV-SWM1 (GU181406) and PDV-PA63 (GU181405) Poland isolates from sweet cherry and PDV-SOF15P11 (GU181404) and PDV-SOF17P17 (GU181399) Italian isolates from sour cherry.ORF3a 882 nt and ORF3b 657 nt long were identified, which encode the MP and CP of 293 and 218 aa, respectively.ORF3a and ORF3b were separated by the intergenic region (IR) of 71 nt.The 5 -and 3 -NCRs consisted of 254 and 260 nt, respectively.The lysine (K) and arginine (R)-rich motif KxxxKxxRxxxxRxxxKxK and the crucial arginine residue, both of which are responsible for binding to the 3 -NCRs of viral RNAs [8], were detected in the N-termini of the MP and CP, respectively.
Thus, the nearly complete genome of the Russian PDV isolate from plum was sequenced.Three genomic RNAs of the virus were assembled and four ORFs typical of PDV were identified.The genome sequences of the P11 isolate were deposited in GenBank under accession numbers OQ995020 to OQ995022.

Phylogenetic Analysis of PDV
For phylogenetic analysis, all complete sequences of three PDV RNAs available in GenBank on May 2023 (Table S2) were used.However, their alignments showed that the 5and 3 -NCRs of different isolates were heterogeneous and varied significantly in length.Therefore, an ORF sequence-based phylogeny was generated.In the case of RNA3, the analyzed genome region also included the IR separating ORF3a and ORF3b.
ORF1, ORF2, and ORF3a-IR-ORF3b were divided into three main clusters (Figure 2a-c) that were in agreement with the previous phylogenies based on the CP, MP, or RNA3 sequences [21][22][23][24][25].Most isolates were affiliated with the same cluster, irrespective of the analyzed genome segment.At the same time, a number of isolates (13C206, 13C226, 13C227, 13C233, 13C257, 13C258, 13C259, 13C277, 13TF105, Niagara D5, and 19SP003, all from Canada) relocated between clusters, depending on the genome segment studied.The genome segments of the isolate P11 were all assigned to cluster II (Figure 2a-c).However, the composition of this group differed from one tree to another, due to both the mentioned migration of some isolates among clusters and the presence of many new isolates in the tree reconstructed from RNA3-derived sequences.Thus, the immediate surroundings of P11 in the cluster varied from tree to tree.Nevertheless, the clustering of the P11 ORF1 and ORF3a-IR-ORF3b was consistent with the sequence identity analysis, while ORF2 clustering did not match the BLASTn results.

Characterization of the PNRSV Genomes
Three PNRSV-related contigs were assembled from the reads generated in both the P10 and P13 samples.In P10, the contig of 3314 nt was most closely related (over 99% identity) to RNA1 of the sour cherry isolates Cigany (MF069044) and PV-0962 (OP357942) from the Czech Republic and Germany, respectively.This contig covered RNA1 of these isolates nearly completely, except for several nt at the 5 -and 3 -termini.ORF1 of 3138 nt encoding replicase of 1045 aa was identified and MET and HEL domains were found at positions 51-387 and 753-1010 of the protein.ORF1 was flanked with the 5 -and 3 -NCRs of 16 and 160 nt, respectively.The 2568 nt long contig was most closely related (over 99% identity) to RNA2 of the Prunus isolates PV-0962 (OP357943), TNpeach5 (OL800570), Cigany (MF069045), and CNU1 (LC752755) from Germany, the USA, the Czech Republic, and South Korea, as well as from the Australian isolate SA from the honey bee Apis mellifera (MH427284).ORF2 of 2400 nt encoding a protein of 787 aa was identified in the contig.RdRp motif was found at positions 298-669 of this protein.ORF2 was flanked with 5and 3 -NCRs of 11 and 157 nt, respectively.The third contig of 1944 nt was most closely related (99% identity) to RNA3 of the Prunus isolates PV-0962 (OP357944) and Cigany (MF069046) from Germany and the Czech Republic.ORF3a 852 nt and ORF3b 675 nt long were identified, which encode the MP and CP of 283 and 224 aa, respectively.ORF3a and ORF3b were separated by the IR of 74 nt.The 5 -and 3 -NCRs consisted of 174 and 169 nt, respectively.The RNA-binding motif (KxxxRxxKKxxRxxKxKxKxR) and two RNA-binding arginine residues [8] were detected in the N-termini of the MP and CP, respectively.
Nearly complete genomes of PNRSV isolates P10 and P13 were deposited in GenBank under accession numbers OR090897 to OR090902.

Phylogenetic Analysis of PNRSV
As for PDV, the phylogeny was based on ORFs and the analyzed genome region of RNA3 included the IR between ORF3a and ORF3b.All complete sequences of three PNRSV RNAs available in GenBank on May 2023 (Table S3) were used for the phylogenetic analysis.
Most PNRSV isolates grouped into two or three main clusters, while others occupied an intermediate position between them (Figure 3a-c).This result showed high PNRSV genetic diversity, in accordance with previous data [22,[26][27][28][29]. ORF1, ORF2, and ORF3a-IR-ORF3b of P10 were assigned to cluster I, which consisted of a large number of isolates found on different stone fruits in many regions of the world.In contrast, there was no clear-cut position for P13 at the phylogenetic trees.At the same time, positions of the genomic segments of two Russian isolates on the phylogenetic trees corresponded to the sequence identity analysis.Like PDV, a few isolates (13TF133, CH57, Acot, P. cerasifera 1, Pch12, 13C227, 13C260, and Rannaja 46) migrated between clusters, depending on the analyzed genomic segment.

Recombination in the PDV and PNRSV Complete Genomes
Three genomic RNAs of some PNRSV and PDV isolates were assigned to different phylogenetic clusters (Figures 2a-c and 3a-c).Recombination, which can be the reason for the incongruence of phylogenetic trees [30], was detected in RNA3 of a few Slovakian PDV isolates [21].The ORF-based recombination analysis was employed using RDP4.101program.Several recombination events were inferred in PDV and PNRSV isolates (Table 2).They were detected in each genomic RNAs of these viruses, except for PDV RNA2.For example, the PDV isolates 13C257 and 13C258, as well as PNRSV isolates 13C257 and 13C258, were recombinants, in which the MP gene was affected.Recombination could explain the relocation of these isolates to group III at the RNA3-based phylogenetic trees closer to potential major parents (Figures 2c and 3c, respectively).However, no recombination was detected in the rest of the relocated PDV and PNRSV isolates or in the Russian isolates P10, P11, and P13.

Partial Characterization of Ourmia-like Virus Genome
Reads related to OuLV1 were assembled into a contig of 806 nt.This contig covered approximately one-third (positions 585 to 1392) of the genomic RNA of the OuLV1 isolates CREA-VE-6AS1 (NC_076380, 97.4% identity) and DMS10_23903 (MN532617, 96.4% identity) from the plant pathogenic fungi Cladosporum cladosporioides and Plasmopara viticola, respectively.The partial genome sequence of the fungal virus isolate named P10-Fun was deposited in GenBank under accession number OR183588.The detection of a mycovirus indicated that the tree was affected by a fungal infection, but the fungus species was not determined.Possibly, the symptoms observed on the plum leaves (Figure 1, sample P10) were of fungal origin.

Discussion
Virus diseases of stone fruits are of great concern due to their wide distribution and high economic significance.The timely and sensitive diagnosis of viruses is a clue in controlling virus diseases, tracing sources of infection, and studying the diseases' impacts on infected trees.However, virus diagnosis in fruit trees can be difficult, due to the high genetic diversity of viruses, their presence in latent form, or uneven distribution in the host [31].The HTS approach potentially allows the finding of all viruses infecting a plant and characterizing their genomes.
In this work, an attempt was made to identify viruses in the field Prunus germplasm collection using HTS.The virome of symptomatic Prunus trees was studied.Based on HTS results, PNRSV, PDV, and fungal ourmia-like virus were detected.Although PDV is observed mainly on cherries (Table S2), in this work it was found on plum.Both ilarviruses were confirmed by RT-PCR and ELISA tests.The foliar symptoms looked not quite common for these ilarviruses.However, ilarvirus symptoms on stone fruits are known to be variable and depend on the isolate, the host, and the environmental conditions, ranging from typical symptoms to the absence of visible manifestations of the disease.In addition, PNRSV symptoms in the P10 sample may be masked by a putative fungal infection.As the studied trees bloom and bear fruit every year, the viruses could infect them through pollen.Samples P10 and P13 could be also infected through myrobalan (P.cerasifera) rootstock, which is no longer available and cannot be tested for viruses.Considering the economic impact of PNRSV and PDV on stone fruit crops, a broader survey is needed to determine the prevalence and genetic diversity of these viruses in the collection and limit their further spread.
Complete genomes of two divergent PNRSV isolates and one PDV isolate were sequenced.Four ORFs and all motifs in the virus proteins, typical for these viruses, were identified in the new isolates.This is the first report and characterization of the PNRSV and PDV full-length genomes from Russia, expanding the information on their geographical distribution and genetic diversity.Multiple alignments of all genome segments of the Russian PNRSV and PDV isolates, together with other isolates of these viruses, available in GenBank, revealed no specific features in the Russian isolates, apart from point mutations randomly distributed along the genome.
According to BLASTn and phylogenetic analysis, the Russian isolate P10 was very close (over 99% identity, irrespective to the analyzed genome segment) to the PNRSV isolates Cigany, PV-0962, TNpeach5, and CNU1 found on various Prunus species in the Czech Republic, Germany, the USA, and South Korea.P13 was most closely related to other PNRSV isolates from Prunus spp., such as 13C241 (RNA1, 97.5% identity), Ruzyne and Pea1 (RNA2, 95.5% identity), and Ruzyne (RNA3, 96.9% identity), that originated from the Czech Republic, China, and Canada.Thus, BLASTn and phylogenetic analysis showed that RNA1, RNA2, and RNA3 of P10 were most similar to the same PNRSV isolates available in GenBank.Three genome segments of P13 were also most closely related to PNRSV isolates available in GenBank., although unlike those of the P10, they were known isolates of this virus.In contrast, each of three genome RNA of the P11 was closest to quite different PDV isolates, sharing 92 to 96% identity (Section 3.3, Figure 2a-c).In addition, the CP gene of the P11 was shown to be closest to the Turkish PDV isolates Sa12 (KF718675), PD8 (EF524271), and M27 (KF718668) from sweet cherry (98% identity).(They were not included in the alignments arranged from complete ORFs.)Thus, the position of the Russian PDV isolate, as well as others, remains uncertain.As the number of full-length PDV genomes increases, the position of the Russian PDV isolate on phylogenetic trees will possibly become clearer.
Taken together, this shows that the new Russian isolates from plums and ((P.cerasifera x P. armeniaca) x P. brigantiaca) trihybrids were the closest to the PNRSV and PDV isolates found on various stone fruit crops in different countries and continents.This suggests that PNRSV and PDV isolates from different hosts, and of different geographical origin, can be very similar.
Phylogenetic analysis of the ORF1, ORF2, and ORF3a-IR-ORF3b genome segments of PDV and PNRSV showed substantial genetic diversity of these viruses, consistent with previous phylogenies of their RNA3 derived sequences [21][22][23][24][25][26][27][28][29].Despite this, no geographical clustering was detected for the full-length genomes of the PDV and PNRSV isolates compared in this work.One of the explanations may be that these viruses are currently continuing to spread, constantly mixing the existing populations.
In addition, phylogenetic analysis of genome-wide sequences revealed no host clustering among PNRSV and PDV isolates.This was perhaps due to easy hybridization between Prunus species [3], which might potentially contribute to the cross-species transmission of the viruses.On the other hand, the vast majority of full-length genome sequences deposited in GenBank are cherry isolates that can distort the results of the host clustering analysis.It cannot be ruled out that further study of PNRSV and PDV isolates from other Prunus species will enable us to explain the lack of the host clustering.
Comparison of the phylogenetic trees reconstructed from three genomic segments showed that some isolates moved from one phylogroup to another, depending on the exploring genome segment.For example, the PDV isolates Niagara D5 and 13C226 were assigned to group III (RNA1), I (RNA 2) or group II (RNA 3).The PDV isolates 13C259 and NS9R relocated from group II (RNA1) to group I (RNA2 and RNA3) (Figure 2a-c).The PNRSV isolate P. cerasifera-1 moved from group I (RNA1) to group II (RNA2 and RNA3), but the immediate surroundings of RNA2 and RNA3 of this isolate were quite different (Figure 3a-c).It should be noted that the mentioned clustering was supported by high bootstrap values (100%).Such relocations of the isolates change the branching order and lead to the incongruence of the trees reconstructed for each of three genome RNAs.Incongruence of the phylogenetic trees may stem from different evolutionary histories of the genome segments of some PDV and PNRSV isolates, due to recombination or reassortment in their ancestors.
Several recombination events were detected in PDV and PNRSV isolates (Table 2).However, these findings cannot explain the different branching order for those isolates that are unaffected by recombination.Reassortment may be another possible explanation.Since members of the genus Ilarvirus have a segmented genome, and each of their three genomic RNAs is encapsidated in a separate viral particle, mixed infection of a plant with different virus isolates can lead to the exchange of genome segments between them and the emergence of reassortants in offspring.In the genus Cucumovirus (which also belongs to the family Bromoviridae), viable reassortants of cucumber mosaic virus, peanut stunt virus, and tomato aspermy virus were obtained experimentally [32] and found in nature [33].Intra-specific reassortment was not studied in ilarviruses, although it has been suggested in the genesis of this genus [33,34].Further study of the putative reassortment in PNRSV and PDV isolates may enhance our understanding of the evolution and genetic diversity of ilarviruses.

Supplementary Materials:
The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/horticulturae9080941/s1,Text file S1: CTAB-based protocol for total RNA isolation; Figure S1: Agarose gel electrophoresis of PCR products generated by RT-PCR assay of ilarviruses using degenerate ilarvirus-specific primers [20].Sample names are shown above the picture.NC-negative control (see Section 2.4 in the text).The arrow to the right to the picture indicates PCR products of the corresponding size.M-GeneRuler 100 bp Plus DNA Ladder (Thermo Scientific).;Table S1: Detection of prunus necrotic ringspot virus (PNRSV) and prune dwarf virus (PDV) by ELISA (A 405 ); Table S2: Prune dwarf virus isolates employed in phylogenetic analysis; Table S3: Prunus necrotic ringspot virus isolates employed in phylogenetic analysis.

18 Figure 1 .
Figure 1.Virus-like symptoms on the leaves of Prunus plants infected with prunus necrotic ringspot virus (P10 and P13) and prune dwarf virus (P11).

Figure 1 .
Figure 1.Virus-like symptoms on the leaves of Prunus plants infected with prunus necrotic ringspot virus (P10 and P13) and prune dwarf virus (P11).

Figure 2 .
Figure 2. Phylogenetic analysis of ORF1 (a), ORF2 (b), and ORF3a-IR-ORF3b genome segment (c) of prune dwarf virus conducted in MEGA7 [18].The names of isolates are shown next to the end of branches.The trees were reconstructed using the neighbor-joining method.The evolutionary distances were computed using the Kimura 2-parameter model.Bootstrap values (>75%) from 1000 replicates are shown next to the corresponding nodes.Curly brackets combine isolates belonging to

Figure 2 .
Figure 2. Phylogenetic analysis of ORF1 (a), ORF2 (b), and ORF3a-IR-ORF3b genome segment (c) of prune dwarf virus conducted in MEGA7 [18].The names of isolates are shown next to the end of branches.The trees were reconstructed using the neighbor-joining method.The evolutionary distances were computed using the Kimura 2-parameter model.Bootstrap values (>75%) from 1000 replicates are shown next to the corresponding nodes.Curly brackets combine isolates belonging to the same cluster.The scale bar indicates the number of substitutions per nucleotide.The Russian isolate P11 is highlighted with a black circle (•).The numbers (I, II, III) belong to the same cluster.

Figure 3 .
Figure 3. Phylogenetic analysis of ORF1 (a), ORF2 (b), and ORF3a-IR-ORF3b genome region (c) of prunus necrotic ringspot virus conducted in MEGA7 [18].The names of isolates are shown next to the end of branches.The trees were reconstructed using the neighbor-joining method and the Kimura 2-parameter model.Bootstrap values (>75%) from 1000 replicates are shown next to the corresponding nodes.Curly brackets combine isolates belonging to the same cluster.The scale bar in-

Figure 3 .
Figure 3. Phylogenetic analysis of ORF1 (a), ORF2 (b), and ORF3a-IR-ORF3b genome region (c) of prunus necrotic ringspot virus conducted in MEGA7 [18].The names of isolates are shown next to the end of branches.The trees were reconstructed using the neighbor-joining method and the Kimura 2-parameter model.Bootstrap values (>75%) from 1000 replicates are shown next to the corresponding nodes.Curly brackets combine isolates belonging to the same cluster.The scale bar indicates the number of substitutions per nucleotide.The Russian isolates P10 and P13 are highlighted with a black circle (•).The numbers (I, II, III) belong to the same cluster.

Table 1 .
Results of high-throughput sequencing of Prunus samples.

Table 1 .
Results of high-throughput sequencing of Prunus samples.

Table 2 .
Recombination events inferred in complete genomes of prune dwarf virus (PDV) and prunus necrotic ringspot virus (PNRSV) by the Recombination Detection Program v.4.101.
* Multiple comparison (MC) corrected probability of false-positive detection of recombination event.