Next Article in Journal
Immune Reconstitution and Safe Metabolic Profile after the Switch to Bictegravir/Emtricitabine/Tenofovir Alafenamide Fumarate among Virologically Controlled PLWH: A 96 Week Update from the BICTEL Cohort
Previous Article in Journal
African Swine Fever Virus Interaction with Host Innate Immune Factors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Global Population Structure of Apple Mosaic Virus (ApMV, Genus Ilarvirus)

1
Department of Plant Protection, Faculty of Agriculture, Bolu Abant İzzet Baysal University, 14030 Bolu, Türkiye
2
Scientifical Industrial and Technological Application and Research Center, Bolu Abant İzzet Baysal University, 14030 Bolu, Türkiye
3
Directorate of Central Plant Protection Research Institute, Gayret Mah. Fatih Sultan Mehmet Bulv., Yenimahalle, 06172 Ankara, Türkiye
4
Department of Plant Protection, Faculty of Agriculture, Universitas Gadjah Mada, Jl. Flora No. 1, Sleman, Yogyakarta 55281, Indonesia
*
Authors to whom correspondence should be addressed.
Viruses 2023, 15(6), 1221; https://doi.org/10.3390/v15061221
Submission received: 12 March 2023 / Revised: 2 May 2023 / Accepted: 9 May 2023 / Published: 23 May 2023
(This article belongs to the Section Viruses of Plants, Fungi and Protozoa)

Abstract

:
The gene sequence data for apple mosaic virus (ApMV) in NCBI GenBank were analyzed to determine the phylogeny and population structure of the virus at a global level. The phylogenies of the movement protein (MP) and coat protein (CP) genes, encoded by RNA3, were shown to be identical and consisted of three lineages but did not closely correlate with those of P1 and P2, suggesting the presence of recombinant isolates. Recombination Detection Program (RDP v.4.56) detected significant recombination signal in the P1 region of K75R1 (KY883318) and Apple (HE574162) and the P2 region of Apple (HE574163) and CITH GD (MN822138). Observation on several diversity parameters suggested that the isolates in group 3 had higher divergence among them, compared to isolates in groups 1 and 2. The neutrality tests assigned positive values to P1, indicating that only this region experiencing balanced or contracting selection. Comparisons of the three phylogroups demonstrated high Fixation index (FST) values and confirmed genetic separation and the lack of gene flow among them. Additionally, ±500 bp of partial MP + ‘intergenic region’ + partial CP coding regions of two Turkish isolates from apple and seven from hazelnut were sequenced and determined that their phylogenetic positions fell within group 1 and 3, respectively.

1. Introduction

Various herbaceous plants such as hop (Humulus lupulus), rose (Rosa spp.), and strawberry (Fragaria × ananassa) and woody plants such as apple (Malus spp.), pear (Pyrus communis L.), plum (Prunus domestica), hazelnut (Corylus avellana L.), cherry (Prunus avium), peach (Prunus persica), apricot (Prunus armeniaca), and lichens worldwide were naturally infected by apple mosaic virus (ApMV, genus Ilarvirus, Family Bromoviridae) [1,2,3,4,5,6,7]. There is no known vector of ApMV [8], but it is probably transmitted through seed and pollen in hazelnut [9]. The virus provides epidemiological challenge as it also persists in the vegetatively propagated material and is readily graft transmissible [10]. While the virus possibly remains latent in the host, the plants associated with ApMV infection exhibit different types of symptoms including vein banding of leaves, yellow discolorations, and oak-leaf patterns [11].
The virus has a tripartite, positive sense, single stranded RNA genome consisting of RNA-1 (3.4 kb in size), RNA-2 (2.9 kb), and RNA-3 (2.0 kb) [12,13]. Non-structural P1 and P2 proteins that include a methyltransferase domain (MT) and polymerase domain (RdRP) involved in viral replication are encoded by RNA1 and RNA2, respectively. Coat protein (CP), expressed via sub-genomic RNA4, and movement protein (MP) are encoded by the bicistronic RNA3 [14,15].
The genetic variation of ApMV has mostly been investigated by the analysis of the CP gene of isolates from various hosts and geographic locations. Initial analysis of the CP gene of ApMV isolates from Korea revealed three evolutionary groupings distinguished by plant hosts [16]. However, later studies on ApMV isolate sequences from diverse geographical origins throughout the world suggested the existence of two to five evolutionary groups. Seven Indian isolates formed a distinct group separated from the other four groups [2]. Apple and rose isolates from Poland and Belarus were positioned together in one of three constructed phylogroups [5]. These older data maintained some correlation of phylogroups and isolate origins, but perhaps due to frequent global trading, recent reports found no clear association of phylogroups with hosts or locations [6,17,18,19]. Despite this knowledge, detailed diversity and evolutionary analyses involving all regions of the ApMV genome has never been performed, making it difficult to have a general agreement on the exact nature of the virus phylogeny and evolution.
Clear understanding on the population structure of a plant virus could be the key for its accurate detection and effective control [20,21]. Therefore, genome sequences of all global ApMV isolates in the National Center for Biotechnology Information (NCBI) GenBank were analyzed in this study to advance our understanding of the virus at the molecular level. Additionally, Turkey’s 684,000 tons of hazelnut production [22] accounts for 69% of the total world production. The country is also the leader of hazelnut export with a world market share of 61%, followed by Italy, the USA, and Azerbaijan. The presence of ApMV in different host species planted in different regions of the country needs to be molecularly identified as a first step for disease control since considerable yield reduction, up to 42%, had been observed in hazelnut [23,24]. Thus, this study could also clarify the position of known Turkish isolates among global populations.

2. Materials and Methods

2.1. Alignment Datasets and Recombination Analysis

On 7 February 2023, 9, 5, 20, and 114 isolates (Supplementary Table S1) with complete P1 gene sequence (spans nts 79–3219 in RNA1; reference isolate: accession no. NC_003464), P2 (spans nts 80–2707 in RNA2; reference isolate: accession no. NC_003465), MP (spans nts 169–1029 in RNA3; reference isolate: accession no. NC_003480), and CP (spans nts 1126–1794 in RNA3; reference isolate: accession no. NC_003480), respectively, were retrieved from NCBI GenBank and then aligned according to the respective genome region using Clustal W algorithm suited in MEGA X software v.10.2.4 [25]. Some divergent isolates have different genome lengths with their respective reference isolate due to indel mutations.
Possible recombination signals on the P1, P2, MP, and CP sequences were searched using the RDP, GENECONV, BootScan, MaxChi, ChiMaera, SiScan, and 3Seq algorithms with a Bonferroni-corrected p value of <0.05 implemented in Recombination Detection Program (RDP v.4.56) [26]. Only those supported by at least five of the algorithms were included as credible signals.

2.2. Phylogenetic and Comparative Nucleotide and Amino Acid Analyses

Four phylogenetic trees based on the P1, P2, MP, and CP genes nucleotide (nt) sequences were built using Neighbor-joining (NJ) algorithms in MEGA X software v.10.2.4. The Kimura 2-parameter models [27] with uniform rate among sites and complete deletion of missing data treatment was determined by the lowest Bayesian information criterion score to be the best for nt substitution of all four alignments. The statistical significance of isolate clusters was determined using 1000 bootstrap replicates. Estimation of the percentage pairwise identities of the P1, P2, MP, and CP coding regions at the nucleotide (nt) and amino acid (aa) levels was calculated using Sequence Demarcation Tool software (SDT v1.2) [28].

2.3. Genetic Diversity, Polymorphism, and Neutral Selection Analyses

The diversities on the nt sequences of the P1, P2, MP, and CP genes were evaluated using DnaSP software v.6.12.03 [29] based on the confidence intervals of population genetics-related parameters: the number of haplotypes = h, haplotype diversity = Hd, average nt differences between sequences = k, nt diversity (per site) = π, total number of mutations = η, and the number of variable sites = S. Transcriptional selection was determined using a ratio of nonsynonymous to synonymous sites (ω = dN/dS). Neutral selection tests implemented in DnaSP v.6.12.03: Tajima’s D [30], and Fu and Li’s D* and F* [31] were performed with a window length of 100 sites and step size of 25 sites to measure genetic divergence on individual coding sequences of the ApMV genome.

2.4. Gene Flow and Genetic Differentiation among Populations

The KS*, KST*, Z*, Snn, and FST values [32,33] were estimated using DnaSP software v.6.12.03. These values determined the genetic differentiation of MP and CP genes nt sequences among compared phylogroups. When genetic divergence among the compared populations is non-existent, the KST* value nears zero [34]. A low genetic isolation is translated into a small Z* value [32]. The Snn value ranges from a minimum of 0.5 (to show that the compared population was identical) to a maximum of 1 (distinctly separated populations) [33]. The FST value ranges between 0, which is describing strictly identical populations, to 1, which is describing totally distinct populations [32,34]. A high gene flow and large genetic isolation between tested populations are usually indicated by an FST value of at least 0.33 [35,36].

2.5. Acquisition of Partial Sequences of Turkish Isolates and Their Phylogeny

Fifteen hazelnut leaf samples from Ordu province and five apple leaf samples from Ankara province, Turkey were collected during the summer of 2022 due to recent reports of viral symptomatic plants in several plantations. Total RNA was extracted from samples using NucleoZOL, one phase RNA purification, following the manufacturer’s procedures (Macherey-Nagel GmbH & Co. KG, Düren, Germany).
Nucleotide sequences of RNA3 of ApMV apple isolates (acc. no. HE574164 (India), KY883320 (Australia), KY971021 (Canada), and MT303103 (Ethiopia)) and hazelnut isolates from Poland (acc. no. HG328282-HG328285) acquired from NCBI GenBank database were aligned using the Clustal W [37] implemented in the MEGA X v.10.2.4 software [25]. The alignment was subjected to Primer-BLAST [38] to generate ApMV-1F (5′-AAGACCCGAAGCCGTAGTTG-3′)/ApMV-1R (5’-GCAAGATCCAGGGTGAGTGT-3′) primers applied in the amplification of a 527 bp of partial MP + ‘intergenic region’ + partial CP coding regions using single step RT-PCR. The 25 μL reaction mixture was composed as follows: 5 µL 5× GoTaq Flexi Buffer, 25 mM MgCl2, 10 mM dNTPs, 10 μM each of forward and reverse primer, 200 U Reverse-transcriptase enzyme, 0.2 μL RNase inhibitor (5000 U/mL), 0.25 µL GoTaq polymerase enzyme (5 U/µL), 0.4 µL DMSO, and 2 μL total RNA. The thermal cycler program was as follows: 50 °C for 45 min, 95 °C for 10 min, 35 cycles of 94 °C for 60 s, 62.5 °C for 60 s, and 72 °C for 45 s, followed by a final extension at 72 °C for 7 min. RT-PCR products were loaded onto a 1.5% 1× TAE agarose gel stained with Pronosafe DNA fluorescent marking (Condalab, Madrid, Spain). Electrophoresis was run at 80 Volt for 1 h in a machine (ThermoFisher Scientific, Waltham, MA, USA). The presence of any band was visualized under a UV transilluminator (Vilber, Marne-la-Vall’ee, France). Next, RT-PCR products with the expected band size were sequenced bi-directionally using the Sanger method in BM Laboratory Systems, Turkey.
The recovered nt sequences of Turkish isolates were aligned with those of 14 isolates in NCBI GenBank that have full sequence of RNA3. The alignment was trimmed according to Turkish isolates length. A phylogenetic tree based on the alignment was constructed to investigate their positions in global phylogrouping, using the Neighbor-joining (NJ) algorithms in MEGA X software v.10.2.4 with Kimura 2-parameter as the best supported nt substitution model [27], uniform rate among sites, and complete deletion of missing data treatment. The statistical significance of isolate branching was determined using 1000 bootstrap replicates. Differences in the partial aa sequence of the aligned isolates were examined using BioEdit v.7.2.5 [39].

3. Results

3.1. Phylogenetic and Recombination Analyses

The constructed MP and CP phylogenetic trees showed exactly the same topology, in which global isolates were separated into three main lineages. Fourteen isolates with complete RNA3 sequences were consistently clustered in the same groups in the respective trees (Figure 1; Supplementary Table S1). However, the phylogenies of the P1 and P2 genes could not be resolved in the same way, probably due to recombination events that were found in the P1 region of K75R1 (KY883318) and Apple (HE574162) and the P2 region of Apple (HE574163) and CITH GD (MN822138) (Table 1). Recombinant isolates were removed from subsequent analyses to improve the accuracy of the results.

3.2. Comparative Nucleotide and Amino Acid Analyses

According to SDT analysis, the MP and CP genes of ApMV isolates maintained very high nt sequence identities of 88–100% and 87–100%, respectively (Table 2). However, the relatively few nt changes in some isolates (Negret 2 for MP and CP, and TrapezK5, TrapezK, and almond for CP) resulted in aa changes, causing low aa identities of them with those of other tested isolates.

3.3. Genetic Diversity and Polymorphism Analyses

Results of the analyses showed that there is a large genomic variation in the sequences of the P1 and P2 genes of few currently known isolates, as shown by the very high Hd, S, η, and k values obtained by both regions. The movement protein gene recorded lower Hd, S, η, and π values than CP, suggesting that the region is highly conserved when compared to the CP (Table 3). However, there were much less MP sequences available to be compared than the CP sequences.
For both MP and CP, isolates in group 3 exhibited greater diversity among themselves than isolates in group 1 and 2, respectively (Table 3). Therefore, the obtained data could also be an indication that group 3 consists of divergent isolates that might be evolutionary adapted to infect certain hosts, especially hazelnut, since no hazelnut isolate was positioned in the other two groups.

3.4. Neutral Selection Analysis

The three parameters of the neutrality test invariably assigned negative numbers, some with statistically significant p values, to three phylogroups of the MP and CP genes, which are in RNA3, and gave positive numbers to the P1 gene, which is in RNA1 (Table 4).

3.5. Gene Flow and Genetic Differentiation among Populations

Analysis was performed only on the MP and CP genes since there was no clear phylogrouping of the P1 and P2 genes due to a lack of available recombinant free sequences currently listed in NCBI GenBank. The comparisons among groups 1, 2, and 3 of both MP and CP raised high and statistically significant KS*, KST*, Z*, and Snn metrics and the Fixation index (FST) > 0.33, which thus supported grouping of the constructed phylogenetic trees into three lineages. However, the values presented by FST analysis on the sequences of the MP and CP genes suggested that there were more divergences between group 1 vs. 2 than group 1 vs. 3 and 2 vs. 3 (Table 5).

3.6. Phylogeny of Novel Turkish Isolates

All of 20 samples were positive for ApMV infection according to results of a molecular test using RT-PCR (Figure 2). Infected hazelnut showed distinctive ‘oak-leaf patterns’ symptoms, whereas infected apple exhibited white discoloration symptoms (Figure 3). Sequencing of nine positive RT-PCR products gave around 500 nts of the ApMV genome, covering partial MP + ‘intergenic region’ + partial CP gene (spans nts 637–1156 in RNA3; reference isolate: accession no. NC_003480). Some Turkish isolates showed different nt and aa lengths due to indel mutations (Figure 4). Newly obtained sequences were registered in NCBI GenBank with accession no. OP374166-OP374174.
The constructed phylogenetic trees positioned two Turkish isolates from apple and six other apple isolates in group 1, whereas seven Turkish isolates from hazelnut were clustered with other four hazelnut isolates identified in Poland in group 3 (Figure 5). Therefore, there was some degree of correlation between ApMV phylogeny with host species. The results were in line with those of polymorphism analysis, which indicated that group 3 belongs to evolutionary successful isolates capable to acquire hazelnut as an additional host. Results of the molecular study also verified the new primer pair capability to amplify partial genome of genetically diverse isolates belong to different hosts and phylogroups.

4. Discussion

ApMV is a major concern in stone fruits and other high value crops cultivated around the world, but its complete phylogeny is still unclear. In this current work, all complete sequences of the P1, P2, MP, and CP regions of the ApMV genome currently listed in NCBI GenBank were analyzed to better understand the population structure and evolution course of the virus. ApMV has also been studied quite well in Turkey, and some isolates have been sequenced. However, sequences of hazelnut isolate, which was reported to be genetically divergent [6], are not yet available. Therefore, hazelnut and apple isolates from Ordu and Ankara provinces obtained in this study were partially sequenced and then studied to know their molecular characters.
RDP analysis of this study did not only confirm the recombination in the P2 region of MN822138 recently reported by [11], but also found strong signals in the P2 region of HE574163 and in the P1 region of HE574162 and KY883318. However, a recombination event was not found in the sequences of the MP and CP genes. Phylogenetic trees further showed that while MP and CP shared an exact phylogeny, P1 and P2 did not have the same phylogeny. Given the frequently occurred signals on these few known sequences of P1 and P2, it is likely that recombination plays more important roles in shaping ApMV evolution, particularly in the two regions than reassortment. However, more sequences of P1 and P2 are clearly needed to draw a strong conclusion on this issue.
There was no strong agreement on ApMV phylogeny as most of the phylogenetic trees of the virus were constructed using the Neighbor-joining (NJ) method [2,5,6,18], while at least one applied Maximum-likelihood [17]. Thus, these proposed trees sometimes produced rather distinctive branches. To resolve this, phylogenetic trees constructed in the current study employed the most commonly used NJ method. Furthermore, unrooting the trees from any out-group seemed to increase the ease and accuracy of phylogrouping. As a result, the trees based on the complete sequences of MP and CP regions (both are in RNA3) showed exactly the same topology in which both trees were divided into three main phylogroups (groups 1, 2, and 3) (Figure 1). Prominently, the trees showed no correlation between the groups and hosts and origin, indicating a wide host range and distribution of any ApMV strain. Phylogrouping proposed in this study could be a strong base for future molecular study of ApMV since all GenBank isolates from different hosts and countries were involved in the analysis.
Both MP and CP regions actually maintained relatively high sequence identity among the isolates according to the SDT study. However, changes in the nt sequences of several isolates that were mostly recovered from hazelnut and belong to a subcluster in group 3 (Negret 2 for MP and CP) and TrapezK5, TrapezK, and almond for CP were translated into many non-synonymous substitutions; thus, their aa sequences shared very low identity with other isolates. Out of curiosity, separate polymorphism and neutral selection analyses were performed on the CP sequences of the 10 isolates in this subgroup within group 3, named ‘Hazelnut’ (all the available hazelnut isolates already belong to group 3 in MP comparison). In the CP comparison, the only 10 isolates in the Hazelnut subcluster were shown to have higher variance among them, according to all diversity parameters, and experienced weaker selection pressure (ω = 0.971) than those of the whole members of groups 1 and 2 (Table 3). The Hazelnut subgroup also received negative values from three neutrality tests (Table 4). Therefore, the results further distinguished the divergence of these isolates that adapted to infect hazelnut.
In this study, seven hazelnut isolates from Turkey were partially sequenced to complement our knowledge on the only known hazelnut isolates obtained in Poland [6]. Results of this phylogenetic study showed that even though new hazelnut and apple isolates were collected from nearby provinces within the same country, the Turkish hazelnut isolates were related much closer to other Polish isolates from hazelnut than to Turkish apple isolates. As a matter of fact, the Turkish apple isolates were clustered with other apple isolates from different countries (Figure 5). Similar to the Polish isolates, there were nt deletions in the RNA3 sequence of Turkish hazelnut isolates at positions nts 1003–1005 and 1089–1091 (reference isolate: accession no. NC_003480), which were reflected into deletion of one aa residue in the respective positions at translation of this segment. However, all seven Turkish isolates also experienced unique additional ten nt deletions at positions nts 1098–1100 and 1106–112 that resulted in the further removal of a total of four aa residues (Figure 4). Most of these aa changes at RNA3 between the apple and hazelnut isolates happened at the non-coding intergenic region (Figure 4), thus unlikely to affect isolates fitness.
Overall, all ApMV phylogroups were under weak evolutionary constraint, which could be one of the reasons for high genetic diversity among isolates in all three groups. In line with this finding, it can also be suggested that different ApMV populations, which all were constantly assigned negative neutrality values, are experiencing further expansion (population growth due to lack of subdivision) or bottleneck selections. Interestingly, these results were rather contrary to those experienced by prunus necrotic ringspot virus (PNRSV), another member of Ilarvirus that also infects stone fruits worldwide, which was found under very strong purifying pressures and received positive neutrality values for some of its populations [40].
High and statistically significant KS*, KST*, Z*, and Snn metrics obtained by all group comparisons confirmed that the clustering of isolates into three phylogroups in both MP and CP trees has been done correctly. Furthermore, estimation of genetic separateness based on FST values > 0.33 showed that members of these groups were genetically distinct, and there was a lack of gene flow among them.

5. Conclusions

The P1, P2, MP, and CP regions in the genome of all the isolates registered in NCBI Genbank were analyzed to study ApMV phylogeny and also to understand the global population structure of the virus. Both MP and CP trees shared the same topology in which isolates were clustered into three major groups (1, 2, and 3) without correlation with host species and countries. Two new Turkish isolates from apple and seven from hazelnut were positioned within group 1 and 3, respectively. However, the clustering cannot confirm P1 and P2 phylogenies yet due to a lack of available data. For MP and CP, calculation on different parameters estimated relatively higher genetic diversity among isolates in group 3 than those in groups 1 and 2. P1 was the only genome region sustaining a balanced or contracting selection. It is clear that additional molecular data are still needed to complete the analysis on the P1 and P2 regions.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v15061221/s1, Supplementary Table S1: ApMV isolates tested in this study.

Author Contributions

Conceptualization, A.Ç. and A.I.S.; methodology, A.Ç., A.I.S. and A.F.M.; formal analysis, A.Ç. and A.I.S.; investigation, A.Ç., A.F.M. and S.C.; resources, S.C.; data curation, S.C.; writing—original draft preparation, A.Ç. and A.I.S.; writing—review and editing, A.F.M. and S.C.; visualization, A.I.S.; supervision, A.Ç.; project administration, A.F.M. and S.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The new data presented in this study are openly available in the NCBI GenBank database, with accession numbers OP374166-OP374174.

Acknowledgments

Authors are grateful to the Bolu Abant Izzet Baysal University, and Directorate of Central Plant Protection Research Institute, Ankara for their permission to use laboratories and necessary equipment.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Petrzik, K.; Lenz, O. Remarkable variability of apple mosaic virus capsid protein gene after nucleotide position 141. Arch. Virol. 2002, 147, 1275–1285. [Google Scholar] [CrossRef]
  2. Lakshmi, V.; Hallan, V.; Ram, R.; Ahmed, N.; Zaidi, A.A.; Varma, A. Diversity of Apple mosaic virus isolates in India based on coat protein and movement protein genes. Indian J. Virol. 2011, 22, 44–49. [Google Scholar] [CrossRef]
  3. Cho, I.S.; Cho, J.D.; Choi, S.K.; Choi, G.S. Occurrence of stone fruit viruses on peach trees (Prunus persica L. Batsch) in Korea. Res. Plant. Dis. 2012, 18, 391–395. [Google Scholar] [CrossRef]
  4. Petrzik, K.; Vondrák, J.; Barták, M.; Peksa, O.; Kubešová, O. Lichens—A new source or yet unknown host of herbaceous plant viruses? Eur. J. Plant Pathol. 2014, 138, 549–559. [Google Scholar] [CrossRef]
  5. Valasevich, N.; Cieślińska, M.; Kolbanova, E. Molecular characterization of Apple mosaic virus isolates from apple and rose. Eur. J. Plant Pathol. 2015, 141, 839–845. [Google Scholar] [CrossRef]
  6. Cieślińska, M.; Valasevich, N. Characterization of Apple mosaic virus isolates detected in hazelnut in Poland. J. Plant Dis. Prot. 2016, 123, 187–192. [Google Scholar] [CrossRef]
  7. Grimová, L.; Winkowska, L.; Zíka, L.; Ryšánek, P. Distribution of viruses in old commercial and abandoned orchards and wild apple trees. J. Plant Pathol. 2016, 98, 549–554. [Google Scholar] [CrossRef]
  8. Postman, J.D.; Mahlenbacher, S.A. Apple mosaic virus in hazelnut germplasm. Acta Hortic. 1994, 351, 601–605. [Google Scholar] [CrossRef]
  9. Aramburu, J.; Rovira, M. Incidence and natural spread of Apple mosaic ilarvirus in hazel in north-east Spain. Plant Pathol. 2000, 49, 423–427. [Google Scholar] [CrossRef]
  10. Digiaro, M.; Savino, V.; Terlizzi, B.D. Ilarviruses in apricot and plum pollen. Acta Hortic. 1992, 309, 93–98. [Google Scholar] [CrossRef]
  11. Nabi, S.U.; Baranwal, V.K.; Rao, G.P.; Mansoor, S.; Vladulescu, C.; Raja, W.H.; Jan, B.L.; Alansi, S. High-throughput RNA sequencing of mosaic infected and non-infected apple (Malus x domestica Borkh.) cultivars: From detection to the reconstruction of whole genome of viruses and viroid. Plants 2022, 11, 675. [Google Scholar] [CrossRef] [PubMed]
  12. Shiel, P.J.; Alrefai, R.H.; Domier, L.L.; Korban, S.S.; Berger, P.H. The complete nucleotide sequence of apple mosaic virus RNA-3. Arch. Virol. 1995, 140, 1247–1256. [Google Scholar] [CrossRef]
  13. Shiel, P.J.; Berger, P.H. The complete nucleotide sequence of apple mosaic virus (ApMV) RNA 1 and RNA 2: ApMV is more closely related to alfalfa mosaic virus than to other ilarviruses. J. Gen. Virol. 2000, 81, 273–278. [Google Scholar] [CrossRef]
  14. Codoner, F.M.; Elena, S.F. The promiscuous evolutionary history of the family Bromoviridae. J. Gen. Virol. 2008, 89, 1739–1747. [Google Scholar] [CrossRef] [PubMed]
  15. Pallas, V.; Aparicio, F.; Herranz, M.C.; Amari, K.; Sanchez-Pina, M.A.; Myrta, A.; Sanchez-Navarro, J.A. Ilarviruses of Prunus spp.: A continued concern for fruit trees. Phytopathology 2012, 102, 1108–1120. [Google Scholar] [CrossRef]
  16. Lee, G.P.; Ryu, K.H.; Kim, H.R.; Kim, C.S.; Lee, D.W.; Kim, J.S.; Park, M.H.; Noh, Y.M.; Choi, S.H.; Han, D.H.; et al. Cloning and phylogenetic characterization of coat protein genes of two isolates of Apple mosaic virus from ‘Fuji’ apple. Plant Pathol. J. 2002, 18, 259–265. [Google Scholar] [CrossRef]
  17. Grimová, L.; Winkowska, L.; Ryšánek, P.; Svoboda, P.; Petrzik, K. Reflects the coat protein variability of apple mosaic virus host preference? Virus Genes. 2013, 47, 119–125. [Google Scholar] [CrossRef]
  18. Korkmaz, G.; Sipahioglu, H.M.; Usta, M. Survey of Apple mosaic virus in apple-growing provinces of East Anatolia (Malatya and Van) by RNA probe hybridization assay and RT-PCR. Turk. J. Agric. For. 2013, 37, 711–718. [Google Scholar] [CrossRef]
  19. Ertunc, F.; Topkaya, S.; Sezer, A. Distribution and molecular detection of apple mosaic virus in apple and hazelnut in Turkey. Afr. J. Biotechnol. 2014, 13, 3144–3149. [Google Scholar] [CrossRef]
  20. Topkaya, Ş.; Çelik, A.; Santosa, A.I.; Jones, R.A.C. Molecular analysis of the global population of potato virus S redefines its phylogeny, and has crop biosecurity implications. Viruses 2023, 15, 1104. [Google Scholar] [CrossRef]
  21. Çelik, A.; Coşkan, S.; Morca, A.F.; Santosa, A.I.; Koolivand, D. Insight into population structure and evolutionary analysis of the emerging tomato brown rugose fruit virus. Plants 2022, 11, 3279. [Google Scholar] [CrossRef]
  22. Turkish Statistical Institute (TurkStat). Crop Production Statistics. 2021. Available online: https://data.tuik.gov.tr/ (accessed on 2 August 2022).
  23. Aramburu, J.M.; Rovira, M. Effects of apple mosaic virus (ApMV) on the growth and yield of ‘Negret’ hazelnut. Acta Hortic. 1995, 386, 565–568. [Google Scholar] [CrossRef]
  24. Akbaş, B.; Değirmenci, K. Incidence and natural spread of Apple mosaic virus on hazelnut in the West Black Sea coast of Turkey and its effect on yield. J. Plant Pathol. 2009, 91, 767–771. [Google Scholar]
  25. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular Evolutionary Genetics Analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef] [PubMed]
  26. Martin, D.P.; Murrell, B.; Golden, M.; Khoosal, A.; Muhire, B. RDP4: Detection and analysis of recombination patterns in virus genomes. Virus Evol. 2015, 1, 1–5. [Google Scholar] [CrossRef]
  27. Kimura, M. A simple method for estimating evolutionary rate of base substitutions through comparative studies of nucleotide sequences. J. Mol. Evol. 1980, 16, 111–120. [Google Scholar] [CrossRef]
  28. Muhire, B.M.; Varsani, A.; Martin, D.P. SDT: A virus classification tool based on pairwise sequence alignment and identity calculation. PLoS ONE 2014, 9, e108277. [Google Scholar] [CrossRef]
  29. Rozas, J.; Ferrer-Mata, A.; Sánchez-DelBarrio, J.C.; Guirao-Rico, S.; Librado, P.; Ramos-Onsins, S.E.; Sanchez-Gracia, A. DnaSP 6: DNA sequence polymorphism analysis of large data sets. Mol. Biol. Evol. 2017, 34, 3299–3302. [Google Scholar] [CrossRef]
  30. Tajima, F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 1989, 123, 585–595. [Google Scholar] [CrossRef]
  31. Fu, Y.X.; Li, W.H. Statistical tests of neutrality of mutations. Genetics 1993, 133, 693–709. [Google Scholar] [CrossRef]
  32. Hudson, R.R.; Boos, D.D.; Kaplan, N.L. A statistical test for detecting geographic subdivision. Mol. Biol. Evol. 1992, 9, 138–151. [Google Scholar] [CrossRef]
  33. Hudson, R.R. A new statistic for detecting genetic differentiation. Genetics 2000, 155, 2011–2014. [Google Scholar] [CrossRef] [PubMed]
  34. Tsompana, M.; Abad, J.; Purugganan, M.; Moyer, J. The molecular population genetics of the Tomato spotted wilt virus (TSWV) genome. Mol. Ecol. 2005, 14, 53–66. [Google Scholar] [CrossRef] [PubMed]
  35. Hudson, R.R.; Slatkin, M.; Maddison, W.P. Estimation of levels of gene flow from DNA sequence data. Genetics 1992, 132, 583–589. [Google Scholar] [CrossRef] [PubMed]
  36. Coşkan, S.; Morca, A.F.; Akbaş, B.; Çelik, A.; Santosa, A.I. Comprehensive surveillance and population study on plum pox virus in Ankara Province of Turkey. J. Plant Dis. Prot. 2022, 129, 981–991. [Google Scholar] [CrossRef]
  37. Thompson, J.D.; Higgins, D.G.; Gibson, T.J. CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22, 4673–4680. [Google Scholar] [CrossRef]
  38. National Center for Biotechnology Information. Primer-BLAST. 2023. Available online: https://www.ncbi.nlm.nih.gov/tools/primer-blast/index.cgi?GROUP_TARGET=on (accessed on 7 February 2023).
  39. Hall, T.A. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl. Acids Symp. Ser. 1999, 41, 95–98. [Google Scholar]
  40. Çelik, A.; Santosa, A.I.; Gibbs, A.J.; Ertunç, F. Prunus necrotic ringspot virus in Turkey: An immigrant population. Arch. Virol. 2022, 167, 553–562. [Google Scholar] [CrossRef]
Figure 1. Unrooted phylogenetic trees on the basis of complete nt sequences of four regions of the ApMV genome that were generated by the Neighbor-joining method based on Kimura 2-parameter’s model (1000 bootstrap replicates, only >50% values were shown) with uniform rate among sites and complete deletion of missing data treatment implemented in MEGA X software. Phylogenetic analysis of the complete ORFs of (a) P1 gene (3145 nts), (b) P2 gene (2634 nts), (c) MP gene (862 nts), and (d) CP gene (675 nts) are shown. Isolates with a complete RNA3 sequence are printed in bold.
Figure 1. Unrooted phylogenetic trees on the basis of complete nt sequences of four regions of the ApMV genome that were generated by the Neighbor-joining method based on Kimura 2-parameter’s model (1000 bootstrap replicates, only >50% values were shown) with uniform rate among sites and complete deletion of missing data treatment implemented in MEGA X software. Phylogenetic analysis of the complete ORFs of (a) P1 gene (3145 nts), (b) P2 gene (2634 nts), (c) MP gene (862 nts), and (d) CP gene (675 nts) are shown. Isolates with a complete RNA3 sequence are printed in bold.
Viruses 15 01221 g001
Figure 2. Gel electrophoresis visualization of successful RT-PCR amplification of 527 bp of partial MP + ‘intergenic region’ + partial CP coding regions of new Turkish isolates. M ladder; 1 Isolate no. OP374173 from apple; 2 Negative control (uninfected sample); 3 OP374166 from hazelnut; 4 OP374167 from hazelnut; 5 OP374168 from hazelnut; 6 OP374169 from hazelnut; 7 OP374170 from hazelnut; 8 OP374171 from hazelnut.
Figure 2. Gel electrophoresis visualization of successful RT-PCR amplification of 527 bp of partial MP + ‘intergenic region’ + partial CP coding regions of new Turkish isolates. M ladder; 1 Isolate no. OP374173 from apple; 2 Negative control (uninfected sample); 3 OP374166 from hazelnut; 4 OP374167 from hazelnut; 5 OP374168 from hazelnut; 6 OP374169 from hazelnut; 7 OP374170 from hazelnut; 8 OP374171 from hazelnut.
Viruses 15 01221 g002
Figure 3. Symptom variation exhibited by ApMV infected plants in Turkey. (a) ‘Oak-leaf patterns’ on hazelnut leaves; (b) Severe form of ‘oak-leaf patterns’ on hazelnut leaves; (c,d) White patches on leaves of two different apple plants.
Figure 3. Symptom variation exhibited by ApMV infected plants in Turkey. (a) ‘Oak-leaf patterns’ on hazelnut leaves; (b) Severe form of ‘oak-leaf patterns’ on hazelnut leaves; (c,d) White patches on leaves of two different apple plants.
Viruses 15 01221 g003
Figure 4. Multiple alignment of the recovered amino acid sequences of partial MP (positions 1–130) + ‘intergenic region’ + partial CP (positions 166–176) genes of genome of Turkish ApMV isolates using different colors to show variation between apple and hazelnut isolates. ApMV RefSeq NC_003480 was included in the alignment.
Figure 4. Multiple alignment of the recovered amino acid sequences of partial MP (positions 1–130) + ‘intergenic region’ + partial CP (positions 166–176) genes of genome of Turkish ApMV isolates using different colors to show variation between apple and hazelnut isolates. ApMV RefSeq NC_003480 was included in the alignment.
Viruses 15 01221 g004
Figure 5. Unrooted phylogram based on ±500 nts of partial MP + ‘intergenic region’ + partial CP genes of genome of nine new ApMV Turkish isolates (highlighted with black squares) and 14 NCBI GenBank isolates with complete RNA3 sequence (printed in bold). The phylogram was generated by the Neighbor-joining method based on Kimura 2-parameter’s model (1000 bootstrap replicates, only >50% values were shown), with uniform rate among sites and complete deletion of missing data treatment implemented in MEGA X software.
Figure 5. Unrooted phylogram based on ±500 nts of partial MP + ‘intergenic region’ + partial CP genes of genome of nine new ApMV Turkish isolates (highlighted with black squares) and 14 NCBI GenBank isolates with complete RNA3 sequence (printed in bold). The phylogram was generated by the Neighbor-joining method based on Kimura 2-parameter’s model (1000 bootstrap replicates, only >50% values were shown), with uniform rate among sites and complete deletion of missing data treatment implemented in MEGA X software.
Viruses 15 01221 g005
Table 1. Putative recombination events detected in the P1 and P2 regions by RDP4 analysis of ApMV isolates.
Table 1. Putative recombination events detected in the P1 and P2 regions by RDP4 analysis of ApMV isolates.
No.RecombinantParents:
Major/Minor
Breakpoints 1
(Start/End)
RDP Implemented
Method 2 (p Value)
P1
1.KY883318 (Australia)NC_003464 (USA)/unknown2365/2790R (6.564 × 10−21)
G (8.039 × 10−20)
B (2.271 × 10−15)
M (1.136 × 10−15)
C (1.831 × 10−13)
S (3.149 × 10−10)
3S (5.593 × 10−28)
2.KY883318 (Australia)NC_003464 (USA)/KY965059 (Canada)360/3062R (9.084 × 10−12)
G (1.468 × 10−7)
B (8.045 × 10−12)
M (3.028 × 10−8)
C (5.966 × 10−7)
S (1.456 × 10−4)
3S (2.376 × 10−6)
3.HE574162 (India)KY965061 (Canada)/NC_003464 (USA)1288/2063R (6.740 × 10−16)
G (3.307 × 10−14)
B (3.122 × 10−16)
M (3.925 × 10−11)
C (1.939 × 10−11)
S (2.879 × 10−14)
3S (9.453 × 10−16)
4.HE574162 (India)KY965061 (Canada)/unknown2558/2879R (2.457 × 10−3)
G (1.522 × 10−3)
B (1.570 × 10−4)
M (6.996 × 10−4)
C (6.122 × 10−4)
S (4.158 × 10−4)
3S (5.398 × 10−4)
P2
1.HE574163 (India)KY965062 (Canada)/unknown493/1125R (2.001 × 10−8)
G (4.360 × 10−4)
B (4.893 × 10−5)
M (1.883 × 10−9)
C (2.128 × 10−3)
S (3.546 × 10−6)
3S (4.883 × 10−6)
2.MN822138 (India)KY883319 (Australia)/KY965062 (Canada)1132/2458R (3.834 × 10−11)
G (3.646 × 10−10)
B (3.876 × 10−9)
M (1.759 × 10−12)
C (4.075 × 10−6)
S (1.325 × 10−39)
3S (1.322 × 10−26)
1 Position in alignment; 2 R = RDP; G = GENECOV; B = BootScan; M = MaxChi; C = Chimaera; S = Siscan; 3S = 3Seq.
Table 2. Nucleotide (nt) and amino acid (aa) identities in the P1, P2, movement (MP), and coat protein (CP) genomic regions among ApMV phylogroups.
Table 2. Nucleotide (nt) and amino acid (aa) identities in the P1, P2, movement (MP), and coat protein (CP) genomic regions among ApMV phylogroups.
Genomic RegionIdentity (%)
ntaa
P191–9978–99
P287–9689–96
MP88–10069–100
group 1 × group 295–9783–92
group 1 × group 388–9069–75
group 2 × group 388–9070–77
CP87–10069–100
group 1 × group 293–9477–88
group 1 × group 387–9169–78
group 2 × group 387–9169–78
Table 3. Summary of genetic diversity and polymorphism analyses of four ApMV genomic regions from different populations.
Table 3. Summary of genetic diversity and polymorphism analyses of four ApMV genomic regions from different populations.
PhylogroupsNhHdSηkπdSdNω
P1771.000389407177.1430.0560.0920.0470.511
P2331.000381389256.6670.1120.2750.0470.171
MP20170.98417419446.6630.0550.0590.0530.898
Group 113110.974383811.5900.0130.0180.0120.667
Group 2320.667171711.3330.0130.0210.0110.524
Group 3 (=Hazenut)441.00010210555.3330.0640.0740.0620.838
CP114870.99425732137.8700.0630.0730.0590.808
Group 175540.99013614917.0810.0270.0310.0260.839
Group 2960.889636523.5830.0350.0380.0320.842
Group 330270.99119823549.0850.0760.0940.0690.734
Hazelnut10101.00016118358.8220.0890.1010.0980.971
N: number of isolates, h: number of haplotypes, Hd: haplotype diversity, S: number of variable sites, η: total number of mutations, k: average number of nucleotide differences between sequences, π: nucleotide diversity (per site), dN: non-synonymous nucleotide diversity, dS: synonymous nucleotide diversity, ω: dN/dS.
Table 4. Summary of demography test statistics of four ApMV genomic regions from different populations.
Table 4. Summary of demography test statistics of four ApMV genomic regions from different populations.
PhylogroupsFu and Li’s D *Fu and Li’s F *Tajima’s D
P10.78257 ns0.77054 ns0.39048 ns
P2ndndnd
MP−0.23355 ns−0.23355 ns−0.60877 ns
Group 1−0.94708 ns−0.86588 ns−0.23742 ns
Group 2ndndnd
Group 3 (=Hazelnut)−0.26300 ns−0.29969 ns−0.35436 ns
CP−2.54107 *−2.34983 *−1.24742 ns
Group 1−3.45373 **−3.1910 **−1.50761 ns
Group 20.29956 ns0.23542 ns−0.07131 ns
Group 3−0.98120 ns−1.03740 ns−0.66893 ns
Hazelnut−0.44376 ns−0.50406 ns−0.45212 ns
* 0.01 < p value < 0.05; ** 0.001 < p value < 0.02, ns = not significant, nd = not determined, at least four isolates are needed for analysis.
Table 5. Genetic differentiation estimates for lineages of ApMV, based on sequences of MP and CP genomic regions comparisons.
Table 5. Genetic differentiation estimates for lineages of ApMV, based on sequences of MP and CP genomic regions comparisons.
ComparisonαKS *αKST *p ValueαZ *p ValueSnnp ValueβFST
MP
Group 1 (n = 13)/Group 2 (n = 3)2.27530.16510.0020 **3.42110.0010 **1.00000.0020 **0.6874
Group 1 (n = 13)/Group 3 (n = 4)2.53390.21970.0000 ***3.54490.0000 ***1.00000.0000 ***0.6591
Group 2 (n = 3)/Group 3 (n = 4)3.15130.20260.0070 **1.42560.0330 *1.00000.0330 *0.6493
CP
Group 1 (n = 75)/Group 2 (n = 9)2.73410.06770.0000 ***6.94050.0000 ***1.00000.0000 ***0.5428
Group 1 (n = 75)/Group 3 (n = 30)2.90910.12890.0000 ***7.21770.0000 ***0.99050.0000 ***0.5093
Group 2 (n = 9)/Group 3 (n = 30)2.73410.06770.0000 ***6.94050.0000 ***1.00000.0000 ***0.5428
* 0.01 < p value > 0.05; ** 0.001 < p value > 0.01; *** p value < 0.001; αKS*, KST*, Z* and Snn are test statistics of genetic differentiation; βFST, coefficient of gene differentiation, which measures inter-population diversity.
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

Çelik, A.; Morca, A.F.; Coşkan, S.; Santosa, A.I. Global Population Structure of Apple Mosaic Virus (ApMV, Genus Ilarvirus). Viruses 2023, 15, 1221. https://doi.org/10.3390/v15061221

AMA Style

Çelik A, Morca AF, Coşkan S, Santosa AI. Global Population Structure of Apple Mosaic Virus (ApMV, Genus Ilarvirus). Viruses. 2023; 15(6):1221. https://doi.org/10.3390/v15061221

Chicago/Turabian Style

Çelik, Ali, Ali Ferhan Morca, Sevgi Coşkan, and Adyatma Irawan Santosa. 2023. "Global Population Structure of Apple Mosaic Virus (ApMV, Genus Ilarvirus)" Viruses 15, no. 6: 1221. https://doi.org/10.3390/v15061221

APA Style

Çelik, A., Morca, A. F., Coşkan, S., & Santosa, A. I. (2023). Global Population Structure of Apple Mosaic Virus (ApMV, Genus Ilarvirus). Viruses, 15(6), 1221. https://doi.org/10.3390/v15061221

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