Strawberry mottle virus (SMoV), classified in the family
Secoviridae, order
Picornavirales (
Sanfacon et al., 2009;
Thompson et al., 2002) is one of the most common viruses on strawberry plantations in Europe. SMoV is transmitted by strawberry aphids
Chaetosiphon sp. and
Aphis gossypii in a semi-persistent manner (
Frazier and Sylvester, 1960). Despite the fact that infected plants usually do not show any characteristic symptoms, severe strains of SMoV may reduce vigor and yield by 20 to 30% (
Mellor and Krczal, 1987). The incidence of SMoV on strawberry plantations has been reported in many countries of Europe, and in Argentina, USA, Canada, New Zealand, and Australia. The studies conducted in Poland so far have been aimed at occurrence and detection of the virus, its impact on the yield of strawberries and reaction of some herbaceous plant species on SMoV infection (
Babini et al., 2004;
Cieślińska, 2004;
Cieślińska and Stankiene, 2004;
Cieślińska and Zawadzka, 1991,
1996). Time-consuming and labor-intensive biological indexing used to be applied as the standard for the detection of SMoV until diagnostics techniques based on the detection of pathogen’s nucleic acids became an alternative to these tests. The development of new detection methods was possible due to the knowledge of the SMoV genome (
Thompson and Jelkmann, 2003;
Thompson et al., 2002,
2003). The positive-sense RNA genome is divided into two segments: RNA1 which encodes polyprotein P1 with domains for a putative RNA helicase, a 3C-like protease and an RNA-dependent RNA polymerase and RNA2 encodes a polyprotein P2 consisting of, among others, movement protein and coat protein (
Thompson et al., 2002). The complete sequence of RNA2 is currently available for nine SMoV isolates found in the Netherlands (
Thompson et al., 2002), Canada (
Bhagwat et al., 2016), and the Czech Republic (
Koloniuk et al., 2018).
The aim of this study was to characterize the molecular properties of the seven SMoV isolates found in strawberry plants in Poland and to examine their genetic variability based on the sequence of RNA2 fragment.
Seven symptomatic Fragaria vesca var. semperflorens ‘Alpine’ indicator plants grafted previously with leaflets of SMoV infected strawberry: Grandarosa’, ‘Pink Rosa’, ‘Pegat’, ‘Markat’, ‘Granat’ (isolates: Granda-1108, Pink-1108, Pegat-1108, Markat-1108, and Granat-1108 isolates, respectively) and two unknown cultivars (isolates: karkas-1108 and Unkn-1108 isolates) were used as the plant material in this study. Mild to severe mottling and chlorotic spots symptoms were observed on the indicators’ leaves. All of the infected plants were maintained in greenhouse in virus collection of the Horticulture Institute in Skierniewice, Poland.
The silica capture method described by
Boom et al. (1990) with adaptations by
Malinowski (1997) was used for extraction of the nucleic acids from the symptomatic leaves of
F. vesca.
One-tube reverse transcription-polymerase chain reaction (RT-PCR) was conducted to amplify the RNA2 fragment of the virus using primers 2SMoV1108-F (5′ GCGGACGCTGTTGCCACTT 3′) and 2SMoV2386-R (5′ TCCGCCGCTCCCACTGATT 3′). This primer pair was designed during this study based on RNA2 sequence of Ontario isolate from Canada (GenBank accession number: KU177219). RT-PCR assay was performed with a thermocycler PTC-200 (MJ Research, Inc., USA), and 7 μl of the amplification products were resolved by electrophoresis in 1.2% agarose gel in 0.5 X TBE buffer, stained in ethidium bromide and visualized on UV transilluminator (Syngen, Cambridge, England). The molecular weight of the PCR products was estimated by comparison with 100 bp DNA ladder (ThermoFisher Scientific). cDNA from a leaf sample of an healthy F. vesca ‘Alpine’ and from SMoV infected strawberry were included as the negative and positive controls, respectively.
Five μl of the RT-PCR products were separately digested with by BfaI, FauI, HaeIII, HincI and TaqI enzymes (Fermentas Life Science) to conduct the restriction fragment length polymorphism analysis (RFLP) of RNA2 fragment of the SMoV isolates. The generated restriction patterns were analyzed by electrophoresis in 8% polyacrylamide gels in 1% TBE buffer.
The RT-PCR products were purified from the gels using QIAquick
® Gel Extraction Kit (Qiagen, Poland) and sequenced in both strains using the same primers employed for amplification. The obtained sequences were aligned using ClustalW of the DNASTAR’s Lasergene 7.1 software (DNASTAR Inc). The consensus sequences were compared with sequences available in GenBank using the BLAST algorithm (
http://ncbi.nlm.nih.gov/BLAST/). The genetic relationships analysis of SMoV strains and the parallel RNA2 nucleotide sequences of the analogous SMoV sequences deposited at GenBank was carried out with the maximum likelihood method and subsequent bootstrap analysis in the MEGA software v. 5.2 (
Tamura et al., 2007).
The presence of SMoV was confirmed in seven
F. vesca var. semperflorens ‘Alpine’ indicator plants by RT-PCR using Smdetncr4a 5 and Sm2ncr1b 5 primers (
Thompson and Jelkmann, 2003). For all the seven samples ~1.3 kb products specific for partial RNA2 were also amplified when using primer pair 2SMoV1108F-2SMoV2386R (data not shown). The further study was aimed to evaluate of the possible genetic variability within this region of the seven detected SMoV isolates.
The RFLP analysis showed a restriction fragment length polymorphism of the amplified RNA2 fragment of the seven virus isolates. Six different profiles were obtained after digestion of RT-PCR products with the enzymes
BfaI,
FauI,
HaeIII,
HincI and
TaqI. Only two isolates - Pink-1108 and Granat-1108 showed the same restriction patterns for each of the separately used enzymes (
Fig. 1,
Table 1). This result indicated that when selecting the suitable restriction enzymes for digesting of the RNA2 fragment, the RFLP technique can be useful and reliable method for the study on genetic variability of SMoV strains.
The nucleotide sequences analysis of the partial RNA2, amplified using the primers designed during this study, confirmed the genetic diversity of the seven Polish SMoV isolates. The similarity of their nucleotide sequence within this region (1213 nt) was 94.7-100%, with the greatest identity between Granat-1108 and Pink-1108 isolates. In turn, lowest nucleotide similarity (94.7%) was found between Granda-1108 and Pegat-1108 isolates. Nucleotide sequences of the RNA2 fragment of the Pegat-1108, Markat-1108, karkas-1108 Granda-1108, Unkn-1108 and Granat-1108 isolates of SMoV were deposited in GenBank under accession numbers MH644000 - MH644005, respectively. The sequences of partial RNA2 of Pink-1108 and Granat-1108 isolates were identical thus only this one for Granat-1108 was submitted to the GenBank database.
A comparative analysis of the nucleotide sequences showed that the isolates found in Poland shared 91.3-92.1% identity with corresponding sequences of the SMoV strain ‘1134’ from the Netherlands (GeneBank accession number: AJ311876), 75.7-99.2% with the strains ‘SMoV_RNA2_A’, ‘SMoV_RNA2_B’, and ‘SMoV_RNA2_C’ found in the Czech Republic (GeneBank accession numbers: MH013325, MH01326, and MH013327, respectively), and 82.1-93.2% with the sequences of Canadian strains ‘Ontario’, ‘NB926’, ‘NSper3’, ‘Nsper17’, and ‘Nsper51’ (GeneBank accession numbers: KU177219, KU200454, KU200457, KU200459, and KU200461, respectively). Partial RNA2 sequence of the Czech strain ‘SMoV_RNA2_C’ showed the highest similarity with Polish isolates ranged from 96.1% for Pegat-1108 to 99.2% for Unkn-1108. In turn, surprisingly, the all seven Polish isolates showed the lowest similarity with the Czech strain ‘SMoV_RNA2_A’ ranging from 75.7% for Granat-1108 isolate to 76.3% for Granda-1108 isolate.
Variability of the RNA2 sequence of the SMoV isolates from different geographic origin allowed to classify them into two main clusters (
Fig. 2). One of them grouped the virus strains from Canada (KU177219, KU200454, KU200457, KU200459, KU200461), the Netherlands (AJ311876), and two strains from the Czech Republic (MH013325, MH013326). The other cluster included seven Polish isolates found in course of this study and the ‘SMoV_RNA2_C’ strain from the Czech Republic (GeneBank accession number MH013327). However, only Unkn-1108 isolate formed the distinct subgroup together with the Czech strain separately from the other Polish isolates. Common origin or spreading of the Czech strain ‘SMoV_RNA2_C’ by aphids could explain its high genetic similarity to seven Polish isolates.
Chaetosiphon sp. have not been found in Poland so far, but probably, due to climate changes these species of aphids have migrated from the warmer regions to our country spreading SMoV. Results of the phylogenetic analysis of amino acids sequences of the SMoV isolates were similar to these obtained in RNA2 nucleotide sequence analysis (data not shown).
This study has shown the presence of genetic diversity among seven Polish isolates of SMoV, and also has confirmed the variability within the strains identified in other countries. The high variability of the partial RNA2 among the Czech strains indicate that the occurrence of different SMoV isolates is not associated with geographical region. Although the newly designed primers allowed amplification of the RNA2 fragment of the all seven SMoV isolates, it is not sure that they would be suitable for the study of genetic variability of the all known and not identified so far virus isolates.