PCR amplicon next generation sequencing (NGS) analysis offers a broadly applicable and targeted approach to detect populations of both high- or low-frequency virus variants in one or more plant samples. In this study, amplicon NGS was used to explore the diversity of the tripartite genome virus, Prunus necrotic ringspot virus (PNRSV) from 53 PNRSV-infected trees using amplicons from conserved gene regions of each of PNRSV RNA1, RNA2 and RNA3. Sequencing of the amplicons from 53 PNRSV-infected trees revealed differing levels of polymorphism across the three different components of the PNRSV genome with a total number of 5040, 2083 and 5486 sequence variants observed for RNA1, RNA2 and RNA3 respectively. The RNA2 had the lowest diversity of sequences compared to RNA1 and RNA3, reflecting the lack of flexibility tolerated by the replicase gene that is encoded by this RNA component. Distinct PNRSV phylo-groups, consisting of closely related clusters of sequence variants, were observed in each of PNRSV RNA1, RNA2 and RNA3. Most plant samples had a single phylo-group for each RNA component. Haplotype network analysis showed that smaller clusters of PNRSV sequence variants were genetically connected to the largest sequence variant cluster within a phylo-group of each RNA component. Some plant samples had sequence variants occurring in multiple PNRSV phylo-groups in at least one of each RNA and these phylo-groups formed distinct clades that represent PNRSV genetic strains. Variants within the same phylo-group of each Prunus plant sample had ≥97% similarity and phylo-groups within a Prunus plant sample and between samples had less ≤97% similarity. Based on the analysis of diversity, a definition of a PNRSV genetic strain was proposed. The proposed definition was applied to determine the number of PNRSV genetic strains in each of the plant samples and the complexity in defining genetic strains in multipartite genome viruses was explored.

The distribution of Ilarvirus species populations amongst 61 Australian Prunus trees was determined by next generation sequencing (NGS) of amplicons generated using a genus-based generic RT-PCR targeting a conserved region of the Ilarvirus RNA2 component that encodes the RNA dependent RNA polymerase (RdRp) gene. Presence of Ilarvirus sequences in each positive sample was further validated by Sanger sequencing of cloned amplicons of regions of each of RNA1, RNA2 and/or RNA3 that were generated by species specific PCRs and by metagenomic NGS. Prunus necrotic ringspot virus (PNRSV) was the most frequently detected Ilarvirus, occurring in 48 of the 61 Ilarvirus-positive trees and Prune dwarf virus (PDV) and Apple mosaic virus (ApMV) were detected in three trees and one tree, respectively. American plum line pattern virus (APLPV) was detected in three trees and represents the first report of APLPV detection in Australia. Two novel and distinct groups of Ilarvirus-like RNA2 amplicon sequences were also identified in several trees by the generic amplicon NGS approach. The high read depth from the amplicon NGS of the generic PCR products allowed the detection of distinct RNA2 RdRp sequence variant populations of PNRSV, PDV, ApMV, APLPV and the two novel Ilarvirus-like sequences. Mixed infections of ilarviruses were also detected in seven Prunus trees. Sanger sequencing of specific RNA1, RNA2, and/or RNA3 genome segments of each virus and total nucleic acid metagenomics NGS confirmed the presence of PNRSV, PDV, ApMV and APLPV detected by RNA2 generic amplicon NGS. However, the two novel groups of Ilarvirus-like RNA2 amplicon sequences detected by the generic amplicon NGS could not be associated to the presence of sequence from RNA1 or RNA3 genome segments or full Ilarvirus genomes, and their origin is unclear. This work highlights the sensitivity of genus-specific amplicon NGS in detection of virus sequences and their distinct populations in multiple samples, and the need for a standardized approach to accurately determine what constitutes an active, viable virus infection after detection by molecular based methods.

In Australia, Stylosanthes little leaf (StLL) phytoplasma has been detected in Stylosanthes scabra Vogel, Arachis pintoi Krapov, Saccharum officinarum L., Carica papaya L., Medicago sativa L., and Solanum tuberosum L. The 16S rRNA gene sequence of StLL phytoplasma strains from S. scabra, C. papaya, S. officinarum and S. tuberosum were compared and share 99.93–100 % nucleotide sequence identity. Phylogenetic comparisons between the 16S rRNA genes of StLL phytoplasma and other ‘Candidatus Phytoplasma’ species indicate that StLL represents a distinct phytoplasma lineage. It shares its most recent known ancestry with ‘Ca. Phytoplasma luffae’ (16SrVIII-A), with which it has 97.17–97.25 % nucleotide identity. In silico RFLP analysis of the 16S rRNA amplicon using iPhyClassifier indicate that StLL phytoplasmas have a unique pattern (similarity coefficient below 0.85) that is most similar to that of ‘Ca. Phytoplasma luffae’. The unique in silico RFLP patterns were confirmed in vitro. Nucleotide sequences of genes that are more variable than the 16S rRNA gene, namely tuf (tu-elongation factor), secA (partial translocation gene), and the partial ribosomal protein (rp) gene operon (rps19-rpl22-rps3), produced phylogenetic trees with similar branching patterns to the 16S rRNA gene tree. Sequence comparisons between the StLL 16S rRNA spacer region confirmed previous reports of rrn interoperon sequence heterogeneity for StLL, where the spacer region of rrnB encodes a complete tRNA-Isoleucine gene and the rrnA spacer region does not. Together these results suggest that the Australian phytoplasma, StLL, is unique according to the International Organization for Mycoplasmology (IRPCM) recommendations. The novel taxon ‘Ca. Phytoplasma stylosanthis’ is proposed, with the most recent strain from a potato crop in Victoria, Australia, serving as the reference strain (deposited in the Victorian Plant Pathology Herbarium as VPRI 43683).