Search

Citrus virus A (CiVA), a novel negative-sense single-stranded RNA virus assigned to the species Coguvirus eburi in the genus Coguvirus, was detected in South Africa with the use of high-throughput sequencing after its initial discovery in Italy. CiVA is closely related to citrus concave gum-associated virus (CCGaV), recently assigned to the species Citrus coguvirus. Disease association with CiVA is, however, incomplete. CiVA was detected in grapefruit (C. paradisi Macf.), sweet orange [C. sinensis (L.) Osb.], and clementine (C. reticulata Blanco) in South Africa, and a survey to determine the distribution, symptom association, and genetic diversity was conducted in three provinces and seven citrus production regions. The virus was detected in ‘Delta’ Valencia trees in six citrus production regions, and a fruit rind symptom was often observed on CiVA-positive trees. Additionally, grapefruit showing symptoms of citrus impietratura disease were positive for CiVA. This virus was primarily detected in older orchards that were established prior to the application of shoot tip grafting for virus elimination in the South African Citrus Improvement Scheme. The three viral-encoded genes of CiVA isolates from each cultivar and region were sequenced to investigate sequence diversity. Genetic differences were detected between the Delta Valencia, grapefruit, and clementine samples, with greater sequence variation observed with the nucleocapsid protein (NP) compared with the RNA-dependent RNA polymerase (RdRp) and the movement protein (MP). A real-time detection assay, targeting the RdRp, was developed to simultaneously detect citrus-infecting coguviruses, CiVA and CCGaV, using a dual priming reverse primer to improve PCR specificity.

Virus-like symptoms on fig tree leaves are a common occurrence worldwide and has mostly been attributed to fig mosaic disease (FMD). Even though only fig mosaic virus (FMV) has been shown to cause FMD, many other viruses have been identified in diseased fig trees. In 2021, total RNA was extracted from fig leaf samples displaying symptoms of mosaic and chlorotic mottling and was subjected to high-throughput sequencing (HTS) to construct the first virome profile of a South African fig tree. Bioinformatic analyses identified FMV, fig leaf mottle-associated virus 1 (FLMaV1), fig leaf mottle-associated virus 2 (FLMaV2), fig latent virus 1 (FLV1), fig badnavirus 1 (FBV1) and grapevine badnavirus 1 (GBV1) in the data. Reverse transcription polymerase chain reaction (RT-PCR) was conducted, for these viruses, on 24 additional fig leaf samples collected in the Western Cape. FBV1, GBV1, FMV, FLMaV1, FLV1, FLMaV2, and fig fleck-associated virus (FFkaV) were detected in 100%, 96%, 92%, 54%, 46%, 21%, and 12.5% of the samples, respectively. This is the first report on the presence of FMV, FLMaV2, FLV1, FFkaV, FBV1 and GBV1 in South Africa and offers a preliminary insight into the virus status of fig trees in the country.

It has been nearly 100 years since citrus growers in two distinct regions in the northern provinces of South Africa noticed unusual symptoms in their citrus trees, causing significant crop losses. They had no idea that these symptoms would later become part of an almost global pandemic of a disease called greening or huanglongbing (HLB). The rapid spread of the disease indicated that it might be caused by a transmissible pathogen, but it took >50 years to identify the causative agent as ‘Candidatus Liberibacter africanus’. Recently, the disease appeared in more African countries, spreading by both infected planting material and Trioza erytreae. To date, five ‘Ca. L. africanus’ subspecies have been identified in various rutaceous species, with ‘Ca. L. africanus subsp. clausenae’ the only subspecies for which a biovar was detected in citrus. Efforts to detect and differentiate HLB-causing Liberibacter species are ongoing, and recent developments are discussed here. This review focuses on aspects of the African form of HLB, including its specific bacterial species and subspecies, its main insect vector, its geographic distribution, and current management strategies.

Agapanthus praecox Willd. is an ornamental flowering plant that is indigenous to southern Africa and was reported to be a host of tomato spotted wilt orthotospovirus (TSWV) in Australia in 2000 (Wilson et al. 2000). Tomato spotted wilt orthotospovirus (TSWV) belonging to the genus Orthotospovirus of the family Tospoviridae is a single-stranded negative sense RNA virus known to cause disease symptoms in many crops and ornamental plant species. This virus is in the top 10 of most economically important plant viruses worldwide (Rybicki 2015; Scholthof et al. 2011). In May 2021, leaf material from three agapanthus (Agapanthus praecox) plants displaying chlorotic mottling, and yellow lesions (Supplementary material 1A) was collected in Mbombela, South Africa. One gram of symptomatic leaf material was used for total RNA extraction from each of the three samples using a CTAB extraction protocol (Ruiz-García et al. 2019). The three RNA extracts were pooled, and a sequencing library was constructed using the Ion Total RNA-Seq Kit v2.0 and RiboMinus™ Plant Kit for RNA-Seq (ThermoFisher Scientific) (Central Analytical Facility (CAF), Stellenbosch University). The RNA library was sequenced on an Ion Torrent Proton Instrument (CAF). A total of 34,392,939 single-end reads were obtained. Data was trimmed for quality with Trimmomatic (CROP:250, MINLEN:50). De novo assembly was performed on the remaining 32,281,645 trimmed reads (average readlength: 100 nt, range: 50-250 nt) using SPAdes 3.13.0 and resulted in 4,788 contigs. BLASTn analysis identified viral contigs longer than 1,000 nucleotides (nts) with high nucleotide (nt) identity to TSWV (6 contigs), as well as to the newly discovered viruses, agapanthus tungro virus (AgTV) (1 contig), and agapanthus velarivirus (AgVV) (4 contigs) (Read et al 2021). Read mapping was performed against the relevant reference sequence with the highest nt identity to the contigs. For TSWV, 4995, 21221 and 14574 reads mapped to segment L (KY250488), M (KY250489) and S (KY250490) of isolate LK-1, respectively resulting in 99.97%, 100.00% and 99.97% genome coverage of the reference accessions. The nt identity between the reference accessions and the consensus sequences generated (OP921761-OP921763) were 97.26%, 97.64% and 97.82% for segment L, M and S. The presence of TSWV was confirmed in the HTS sample using an RT-PCR assay (primers L1 and L2) targeting the L segment of TSWV (Mumford et al. 1994). In July 2022, additional leaf samples displaying symptoms of chlorotic mottling, streaking, and ringspots were collected from 31 symptomatic and 3 asymptomatic agapanthus plants in public gardens in Stellenbosch, South Africa. Using the above-mentioned RT-PCR assay, 13 of the symptomatic samples tested positive for TSWV. All six plants displaying ring spot symptoms (Supplementary material 1B) were infected with TSWV. However, plants that displayed yellow streaking (five samples) and chlorotic mottling (two samples) (Supplementary material 1C-D) were also positive for TSWV which could be due to the presence of other viruses, plant growth stage, infection time or just variable symptom expression in a single host species as reported previously (Sherwood et al. 2003). The 275 bp RT-PCR amplicons of the HTS sample and three additional positive samples were validated with bidirectional Sanger sequencing (CAF) and had 96% identity to accession KY250488. The pairwise nt identity between amplicons was 98.55-100%. This is the first report of TSWV infecting agapanthus in South Africa. This study contributes information towards the distribution and incidence of TSWV and highlights the need for nurseries to screen plant material before propagation.

Plum viroid I (PlVd-I) was recently identified as a new viroid in 2020 present in Japanese plum (

High-throughput sequencing (HTS) was used to construct the virome profile of an old grapevine-leafroll-diseased grapevine (Vitis vinifera). De novo assembly of HTS data showed a complex infection, including a virus sequence with similarity to viruses of the genus Badnavirus, family Caulimoviridae. The complete genome sequence of this virus consists of 7090 nucleotides and has four open reading frames (ORFs). Genome organisation and phylogenetic analysis identify this virus as a divergent variant of grapevine Roditis leaf discoloration-associated virus (GRLDaV) with 90% nucleotide sequence identity to isolate w4 (NC_027131). This is the first genome sequence of a South African variant of GRLDaV.

BackgroundHigh-throughput sequencing (HTS) has been applied successfully for virus and viroid discovery in many agricultural crops leading to the current drive to apply this technology in routine pathogen detection. The validation of HTS-based pathogen detection is therefore paramount.MethodsPlant infections were established by graft inoculating a suite of viruses and viroids from established sources for further study. Four plants (one healthy plant and three infected) were sampled in triplicate and total RNA was extracted using two different methods (CTAB extraction protocol and the Zymo Research Quick-RNA Plant Miniprep Kit) and sent for Illumina HTS. One replicate sample of each plant for each RNA extraction method was also sent for HTS on an Ion Torrent platform. The data were evaluated for biological and technical variation focussing on RNA extraction method, platform used and bioinformatic analysis.ResultsThe study evaluated the influence of different HTS protocols on the sensitivity, specificity and repeatability of HTS as a detection tool. Both extraction methods and sequencing platforms resulted in significant differences between the data sets. Using a de novo assembly approach, complemented with read mapping, the Illumina data allowed a greater proportion of the expected pathogen scaffolds to be inferred, and an accurate virome profile was constructed. The complete virome profile was also constructed using the Ion Torrent data but analyses showed that more sequencing depth is required to be comparative to the Illumina protocol and produce consistent results. The CTAB extraction protocol lowered the proportion of viroid sequences recovered with HTS, and the Zymo Research kit resulted in more variation in the read counts obtained per pathogen sequence. The expression profiles of reference genes were also investigated to assess the suitability of these genes as internal controls to allow for the comparison between samples across different protocols.ConclusionsThis study highlights the need to measure the level of variation that can arise from the different variables of an HTS protocol, from sample preparation to data analysis. HTS is more comprehensive than any assay previously used, but with the necessary validations and standard operating procedures, the implementation of HTS as part of routine pathogen screening practices is possible.

Ornithogalum thyrsoides, a widely cultivated bulbous ornamental plant endemic to South Africa, has significant commercial value as a pot plant and for the production of cut flowers. However, infection by viruses threatens the success of commercial cultivation, as symptoms negatively affect the appearance of the plant and flowers. To date, four Ornithogalum-infecting viruses have been reported. Complete genome sequence data are available for three of these viruses, but the genome of the potyvirus ornithogalum virus 3 (OV3) has not been fully sequenced. In this study, the complete sequence of OV3 was determined by high-throughput sequencing (HTS) and validated by Sanger sequencing. Based on recognition of protease cleavage patterns and multiple sequence alignments with closely related viruses, the polyprotein of OV3 was predicted to be proteolytically cleaved to produce 10 mature peptides containing domains conserved in members of the genus Potyvirus. Phylogenetic analysis and species demarcation criteria confirm the previous classification of OV3 as a member of a separate species in this genus. This is the first report of a complete genome sequence of OV3.

Two isolates of the T68 genotype of citrus tristeza virus (CTV) were derived from a common source, GFMS12, by single aphid transmission. These isolates, named GFMS12-8 and GFMS12-1.3, induced stem pitting with differing severity in ‘Duncan’ grapefruit (Citrus × paradisi [Macfad.]). Full-genome sequencing of these isolates showed only minor nucleotide sequence differences totaling 45 polymorphisms. Numerous nucleotide changes, in relatively close proximity, were detected in the p33 open reading frame (ORF) and the leader protease domains of ORF1a. This is the first report of full-genome characterization of CTV isolates of a single genotype, derived from the same source, but showing differences in pathogenicity. The results demonstrate the development of intragenotype heterogeneity known to occur with single-stranded RNA viruses. Identification of genetic variability between isolates showing different pathogenicity will enable interrogation of specific genome regions for potential stem pitting determinants.

Over the past 2 decades, fruit symptoms resembling a marbling pattern on the fruit skin or corking of the fruit flesh were observed on Japanese plums in South Africa, resulting in unmarketable fruit. The ability of high-throughput sequencing (HTS) to detect known and unknown pathogens was exploited by assaying affected and unaffected fruit tree accessions to identify the potential aetiological agent of marbling and/or corky flesh disease. In this study, it is shown that the disease is associated with a previously undescribed small RNA with typical viroid structural features. The potential viroid was the only pathological agent consistently detected in all symptomatic trees by HTS, and the association with the symptoms was confirmed in field surveys over two seasons. To date, this RNA was not detectable by RT-PCR in seedlings raised from seeds collected from infected trees. Although the autonomous replication of this viroid-like RNA was not proven, it was shown to be transmissible by grafting and associated with a range of symptoms that include marbling on the fruit skin, corky flesh, reduced fruit size, irregular shape, and uneven fruit surface depending on the cultivar. Moreover, the circular RNA genome, consisting of 317 nucleotides, strongly supports that this viroid-like RNA is most likely a viroid for which the name plum viroid I (PVd-I) is proposed. The primary structure of this viroid showed a less than 90% nucleotide sequence identity to viroids of the genus Apscaviroid, with which it has close phylogenetic relationships and shares conserved structural motifs.

Citrus black spot (CBS), caused by

It has been nearly 100 years since citrus growers in two distinct regions in the northern provinces of South Africa noticed unusual symptoms in their citrus trees, causing significant crop losses. They had no idea that these symptoms would later become part of an almost global pandemic of a disease called greening or huanglongbing (HLB). The rapid spread of the disease indicated that it might be caused by a transmissible pathogen, but it took50 years to identify the causative agent as '

Although physical activity is a health-promoting, popular global pastime, regular engagement in strenuous exercises, such as long-distance endurance running races, has been associated with a variety of detrimental physiological and immunological health effects. The resulting altered physiological state has previously been associated with fluctuations in various key metabolite concentrations, however, limited literature exists pertaining to the global/holistic metabolic changes that are induced by such. This investigation subsequently aims at elucidating the metabolic changes induced by a marathon by employing an untargeted proton nuclear magnetic resonance (1H-NMR) spectrometry metabolomics approach. A principal component analysis (PCA) plot revealed a natural differentiation between pre- and post-marathon metabolic profiles of the 30-athlete cohort, where 17 metabolite fluctuations were deemed to be statistically significant. These included reduced concentrations of various amino acids (AA) along with elevated concentrations of ketone bodies, glycolysis, tricarboxylic acid (TCA) cycle, and AA catabolism intermediates. Moreover, elevated concentrations of creatinine and creatine in the post-marathon group supports previous findings of marathon-induced muscle damage. Collectively, the results of this investigation characterize the strenuous metabolic load induced by a marathon and the consequential regulation of main energy-producing pathways to accommodate this, and a better description of the cause of the physiological changes seen after the completion of a marathon.

Plum viroid I (PlVd-I) is found in marbling and corky flesh diseased plum trees in South Africa. In this study a reverse transcription loop-mediated isothermal amplification (RT-LAMP) assay for the high-throughput detection of PlVd-I was developed. This assay can be performed on crude extracts and detection can either be a pH dependent colorimetric reaction or a real-time fluorescent signal reaction. The false discovery rate was shown to be low and no decrease in sensitivity was detected compared to RT-PCR. The RT-LAMP assay allows for the fast and cost-effective detection of PlVd-I that will curtail the distribution of infected plant material.

Author(s): Aldrich, Dirk Jacobus; Bester, Rachelle; Cook, Glynnis; du Plessis, Anton; Burger, Johan Theodorus; Maree, Hans Jacob | Abstract: Citrus tristeza virus (CTV) is the most important viral pathogen of citrus. CTV-induced stem pitting negatively impacts grapefruit and sweet orange production. The mechanisms of stem pitting development in CTV-infected citrus remain unclear. This study evaluated the utility of high-resolution CT scanning as a tool to study stem pitting in live citrus material. CT scans were used to easily identify pits based on differences in tissue density. Stem pits were also mapped and modelled three-dimensionally along the length of the stem. Nano-CT scanning proved to be a potentially valuable, non-destructive method for stem pitting characterization in citrus.

Apple rubbery wood virus 2 (ARWV-2; Rott et al., 2018) belong to the species Apple rubodvirus 2, a member of the genus Rubodvirus (family Phenuiviridae; Kuhn et al 2020). ARWV-2 was first identified in apples and is associated with apple rubbery wood disease (ARWD) that is characterized by unusual flexibility of stems and branches, reduced growth, shortened internodes and increased cold sensitivity (Jakovljevic et al., 2017, Rott et al., 2018). ARWD was first reported in 1935 in England on apple and has since been found on quince and pear (Jakovljevic et al., 2017; Rott et al., 2018). In January 2021, leaves were collected from a pear tree (Pyrus communis cv. Forelle, F514) in a commercial orchard near Villiersdorp, South Africa. The tree displayed no foliar or tree branch symptoms, except for malformed fruits potentially due to insect feeding damage or pear stony pit disease previously associated with infection of apple stem pitting virus (ASPV) (Paunovic et al. 1999). Leaf petioles (one gram) were used for total RNA extraction, using a modified CTAB extraction protocol (Ruiz-Garcia et al. 2019). A sequencing library was constructed (Illumina TruSeq Stranded Total RNA with plant Ribo-Zero) and sequenced on an Illumina HiseqX instrument (Macrogen, South Korea). A total of 30,709,182 paired-end reads (100 nt) were obtained and trimmed for quality with Trimmomatic (SLIDINGWINDOW:3:20, MINLEN:20) (Bolger et al. 2014). De novo assembly, using default parameters of CLC Genomics Workbench 11.0.1 (Qiagen), resulted in 97,294 contigs. BLASTn analysis identified 17 viral contigs, with 14 contigs having high nucleotide identity to ASPV and three to ARWV-2. The latter contigs included all three segments of ARWV-2. The L contig was 7371 nts, M was 1289 nts and S was 1463 nts in length, generated with 7341, 626 and 9161 reads for segment L, M and S, respectively. Segment S had the highest read coverage (524.87x), followed by segment L (88.07x) and M (36.60x). The ARWV-2 GenBank accessions with the highest percentage identity to the contigs were MF062128.1 from United States of America (98.2% to segment L), MN163134.1 from China (97.5% to segment M) and NC_055535.1 from Germany (93.5% to segment S). The contigs spanned 100%, 80.92% and 100% of these accessions of segments L, M and S, respectively and were deposited in GenBank as accessions MZ593725- MZ593727. Reverse transcription polymerase chain reaction (RT-PCR) was used to validate the presence of ARWV-2 in sample F514, using primers directed at segments L (con708_178F/con708_666R), M (ARWaV-2S1_38F/ARWaV-2S1_682R) and S (ARWaV-2M567F/ARWaV-2M1342R) (Rott et al., 2018). Amplicon sequences (510 bp (L), 645 bp (M) and 799 bp (S)) were confirmed with bi-directional Sanger sequencing. Fifty-nine additional pear samples were surveyed in 2021 for ARWV-2 using the M segment assay as mentioned above. The survey included the Koue Bokkeveld and Elgin areas, and cultivars Bosc (22 samples), Abate (10 samples), Rosemarie (3 samples), Forelle (9 samples), Packham's Triumph (12 samples) and Early Bon Chretien (3 samples). A total of 27 samples (11 samples from the Koue Bokkeveld region and 16 samples from the Elgin region) tested positive for ARWV-2, demonstrating the common presence of this virus in pears in South Africa. This is the first report of ARWV-2 infecting pear in South Africa. Although no association with disease symptoms were observed, this study expands the data on the incidence and distribution of this virus in South Africa.

Huanglongbing (HLB, Asian Citrus Greening), the most devastating disease of citrus has not been detected in southern Africa (Gottwald, 2010). HLB is associated with 'Candidatus Liberibacter asiaticus' (CLas), a phloem-limited bacterium vectored by Diaphorina citri Kuwayama (Hemiptera: Liviidae), the Asian Citrus Psyllid (ACP). African Citrus Greening, associated with 'Candidatus Liberibacter africanus' (CLaf) and its vector the African Citrus Triozid, Trioza erytreae (Del Guercio) (Hemiptera: Triozidae), are endemic to Africa, although not previously reported from Angola. African Greening is less severe than HLB, largely due to heat sensitivity of CLaf and its vector. Introduction of HLB into southern Africa would be devastating to citrus production in commercial and informal sectors. Concern was raised that CLas or ACP might hae inadvertently been introduced into Angola. In July 2019, a survey was conducted in two citrus nurseries in Luanda and Caxito and in different orchards on 7 farms surrounding Calulo and Quibala. Yellow sticky traps for insects were placed at the various localities and collected after c. 3 weeks. Breeding signs of T. erytreae (pit galls) were observed on citrus in some locations, but no insect vectors were detected on traps. Trees were inspected for signs and symptoms of citrus pests and diseases, particularly those that resemble HLB (foliar blotchy mottle, shoot chlorosis, vein yellowing and corking, lopsided fruit with aborted seeds and colour inversion) and its vectors (pit galls on leaves or waxy exudates). Leaves and shoots with suspect symptoms were sampled for laboratory analysis (43 samples). DNA was extracted from petiole and midrib tissue of leaves using a modified CTAB extraction protocol of Doyle and Doyle (1990). Real-time PCR was done using universal Liberibacter primers of Roberts et al. (2015), CLaf specific primers of Li et al. (2006) and CLas specific primers of Bao et al. (2019). All real-time PCR protocols indicated the presence of CLaf in 6 samples (Tab. S1). CLas or other citrus Liberibacter species were not detected. The presence of CLaf in sample 37 was confirmed by constructing a library (NEXTFLEX® DNA Sequencing Kit, PerkinElmer) with extracted DNA and performing high-throughput sequencing on an Ion Torrent™ S5™ platform (Central Analytical Facility, Stellenbosch University). To improve the quality of the reads, all 233,617,700 obtained reads were trimmed from the 3' end to a maximum length of 240 nt using Trimmomatic (Bolger et al. 2014). The high quality reads were mapped to the Citrus sinensis reference genome (NC_023046.1) using Bowtie 2.3.4 (Langmead and Salzberg 2012) to subtract all the reads that had high identity to the host plant (number of mismatches allowed in the seed was set to 1). The 14,691,369 unmapped reads (6.2% of original data) were mapped to the CLaf reference genome NZ_CP004021.1 using CLC Genomics Workbench 10.1.1 (Qiagen) (Length fraction = 0.8; Similarity fraction = 0.9). A CLaf consensus genome was generated that spanned 99.7% of the reference genome and the 163001 mapped reads had a 22.9 mean read coverage. The consensus sequence was 99.7% identical to NZ_CP004021.1 and was submitted to Genbank as accession: CP054879. The positive CLaf detections were from trees with typical HLB or African Citrus Greening symptoms, viz. lopsided fruit with green stylar ends, aborted seed and stained columella at base of fruit button; yellow shoots with leaves showing symptoms of blotchy mottle and vein yellowing and corking (Fig. S1) in a commercial citrus farm outside Calulo and included 2 'Ponkan' mandarin (C. reticulata), 2 Valencia and 1 'Navelina' tree (C. sinensis), and a citrus nursery in Luanda (1 lime tree; C. aurantifolia) (Tab. S1). This first report of CLaf in Angola highlights the need to prevent spread by removing infected trees and managing the insect vector, as well as the need for further surveys to determine the occurrence of African Greening and its vectors in other provinces and to confirm the absence of exotic citrus pests and diseases. References Bao, M. et al. 2020. Plant Dis. 104:527 Bolger, A. M. et al. 2014. Bioinformatics. 30:2114-2120. Doyle, J.J. and Doyle, J.L. 1990. Focus 12:13 Gottwald, T.R. 2010. Annu. Rev. Phytopathol. 48:119 Langmead, B. and Salzberg, S. 2012. Nature Methods. 9:357-359. Li, W. et al. 2006. Jnl. Microbiol. Methods 66:104 Roberts, R. et al. 2015. Int. J. Syst. Evol. Micr. 65:723.

High-throughput sequencing (HTS) was used to construct the virome profile of an old grapevine-leafroll-diseased grapevine (Vitis vinifera). De novo assembly of HTS data showed a complex infection, including a virus sequence with similarity to viruses of the genus Badnavirus, family Caulimoviridae. The complete genome sequence of this virus consists of 7090 nucleotides and has four open reading frames (ORFs). Genome organisation and phylogenetic analysis identify this virus as a divergent variant of grapevine Roditis leaf discoloration-associated virus (GRLDaV) with 90% nucleotide sequence identity to isolate w4 (NC_027131). This is the first genome sequence of a South African variant of GRLDaV.

High-throughput sequencing (HTS) was used to investigate ringspots on ivy (Hedera helix) leaves. De novo assembly of HTS data generated from a total RNA extract from these leaves yielded a contig with sequence similarity to viruses of the genus Badnavirus, family Caulimoviridae. The complete genome sequence of this virus consists of 8,885 nucleotides and has three open reading frames (ORFs). Genome organisation and phylogenetic analysis identifies this newly identified virus as a new member of the genus Badnavirus for which we propose the name "ivy ringspot-associated virus" (IRSaV).

Coguvirus eburi is a member of the genus Coguvirus in the family Phenuviridae (Khun et al., 2020). The species Coguvirus eburi was established to include citrus virus A (CiVA), which is a negative-sense, single-stranded RNA virus that was first found infecting sweet orange in southern Italy via high-throughput sequencing (HTS) (Navarro et al., 2018). This virus was also found to infect pome fruits in France, such as pear (Svanella-Dumas et al., 2019). More recently CiVA infections have been associated with impietratura disease in citrus (Beris et al. 2021). In the summer of 2021, leaf samples were collected from a pear tree (Pyrus communis cv. Bosc, B175) in the Koue Bokkeveld, South Africa as part of a virus survey. Sample B175 displayed no visual disease symptoms. One gram of leaf petioles was used for total RNA extraction, using a modified CTAB extraction protocol (Ruiz-Garcia et al. 2019). Ribo-depleted RNA was prepared (Ribo-Zero Plant kit) and a sequencing library constructed (Illumina TruSeq Stranded Total RNA). The RNA library was paired-end (2 × 100 bp) sequenced on an Illumina HiSeqX instrument (Macrogen, South Korea). A total of 47,750,152 reads were obtained. Raw data was trimmed for quality with Trimmomatic (SLIDINGWINDOW:3:20, MINLEN:20) (Bolger et al. 2014). De novo assembly performed with CLC Genomics Workbench 11.0.1 (Qiagen) (Default parameters) using high quality reads yielded 75250 contigs. BLASTn analysis identified two viral contigs with high nucleotide (nt) identity to apple stem pitting virus (ASPV) and CiVA. The CiVA contig was 9400 nts and on closer examination, a concatemer of CiVA RNA1 and RNA2. The concatenation occurred due to the characteristic near-identical nucleotides shared at the 5' and 3' ends of RNA1 and RNA2 of these negative-stranded RNA viruses (Navarro et al., 2018). After splitting and curation, the RNA1 contig was 6664 nts and the RNA2 contig 2686 nts. A total of 51397 and 34820 reads were used to construct these contigs resulting in an average depth of coverage of 761 and 1281 for RNA1 and RNA2, respectively. The contigs had the highest nt identity to the complete CiVA GenBank accessions MT720885.1 (95.53%) and MW148460.1 (96.03%), spanning 99.6% and 98.1 % of the genomes of RNA1 and RNA2, respectively. These contigs were submitted as partial genomes to GenBank as accessions MZ463039 and MZ463040. Reverse transcription polymerase chain reaction (RT-PCR) was used to validate the presence of CiVA in sample B175. Two RT-PCR assays, directed at RNA1 and RNA2 respectively (Bester et al. (2021)) were used to generate amplicons. Amplicon sequences were confirmed with bi-directional Sanger sequencing. Twenty-one additional samples from the same orchard as B175 as well as other samples from the Koue Bokkeveld and Elgin areas, including cultivars Abate (10 samples), Forelle (10 samples), Early Bon Chretien (3 samples), Packham's Triumph (12 samples) and Rosemarie (3 samples), were all surveyed for CiVA using the same RT-PCR assays as mentioned above. Thirty-six of the 59 samples tested were positive for CiVA, which further confirms the presence and wide-spread distribution of this virus in the limited survey conducted in pears in South Africa. However, no association with any disease symptoms or specific cultivar were identified. This is the first report of CiVA infecting pear in South Africa. This study therefore contributed to investigating the distribution of this virus and will assist the South African plant material certification scheme to assess the incidence of CiVA in South Africa.

The molecular interactions between grapevine leafroll-associated virus 3 (GLRaV-3) and three grapevine cultivars were characterised by identifying small RNAs and genes with differential expression between healthy and infected grapevine using next-generation sequencing and RT-qPCR assays. Cultivar specificity was identified in the sRNA and gene expression analyses, and both approaches identified Chenin blanc-specific responses providing the first insight into the differential symptom expression observed between cultivars. Differentially expressed genes identified in all three cultivars, provide a potential pool of genes displaying a molecular response to GLRaV-3 infection, contributing to the characterisation of this virus-plant interaction.

The economic importance of the grapevine industry creates the need for the generation of genomic resources to ensure the sustainability of the industry. These resources will extend the knowledge of grapevine physiology, cultivar-specific characteristics and understanding how diseases and environmental conditions affect the plant. Small non-coding RNAs (sRNAs) are regulator molecules that can influence normal development and/or responses to environmental stimuli. In this study, next-generation sequencing (NGS) and bioinformatic analyses were used to identify sRNA species in grapevine phloem tissue. Forty-five novel mature Vitis vinifera miRNAs were identified and the presence of ten was validated using stemloop RT-qPCRs. Both non-coding and protein-coding gene regions were identified as putative phased loci of which the majority of the protein-coding gene loci were identified as disease resistance proteins. Potential phase-initiating miRNAs were also identified. A high number of sRNAs originating from repeat sequences were identified in all cultivar groups. This study also confirmed the non-random manner in which tRNA-derived sRNAs originate. The number of tRNA halves identified in Chenin blanc was lower compared to the other cultivar groups investigated. The identification of these different sRNAs extends the grapevine sRNA knowledge base significantly and can contribute to the unravelling of sRNA regulation in plants, potentially creating tools for grapevine functional studies.

MicroRNAs (miRNAs) are a class of endogenous small non-coding RNAs (sRNA) that play an essential role in the regulation of target mRNAs expressed during plant development and in response to stress. MicroRNA expression profiling has helped to identify miRNAs that regulate a range of processes, including the plant's defence response to pathogens. In this study, differential miRNA expression in own-rooted Vitis vinifera cv. Cabernet Sauvignon plants infected with grapevine leafroll-associated virus 3 was investigated with microarrays and next-generation sequencing (NGS) of sRNA and mRNA. These high-throughput approaches identified several differentially expressed miRNAs. Four miRNAs, identified by both approaches, were validated by stemloop RT-PCRs. Three of the predicted targets of the differentially expressed miRNAs were also differentially expressed in the transcriptome data of infected plants, and were validated by RT-qPCR. Identification of these miRNAs and their targets can lead to a better understanding of host-pathogen interactions involved in grapevine leafroll disease and the identification of possible targets for virus resistance.

High-throughput sequencing (HTS) is a powerful tool employed by plant virologists for the detection of viruses, the characterization of virus genomes and the study of host-pathogen interactions. Virus detection has been an important application of this technology, which has resulted in the discovery of novel viruses or viral strains as well as for the detection of known viruses in a plant sample. Here we describe the entire process that needs to be considered for the genome analysis of Citrus tristeza virus (CTV) by HTS, including the experimental design, sample preparation, nucleic acid purification, HTS library construction, and bioinformatic analysis.

High-throughput sequencing (HTS) is a powerful tool employed by plant virologists for the detection of viruses, the characterization of virus genomes and the study of host-pathogen interactions. Virus detection has been an important application of this technology, which has resulted in the discovery of novel viruses or viral strains as well as for the detection of known viruses in a plant sample. Here we describe the entire process that needs to be considered for the genome analysis of Citrus tristeza virus (CTV) by HTS, including the experimental design, sample preparation, nucleic acid purification, HTS library construction, and bioinformatic analysis.

A survey of viruses affecting grapevine in the wine regions of the Western Cape Province in South Africa was conducted. The survey determined the relative abundance of five different Grapevine leafroll-associated virus 3 (GLRaV-3) variants. Virus profiles were also determined for individual vines. A total of 315 plants were sampled and analysed over two growing seasons. Five GLRaV-3 variants were detected as either single or mixed infections, with GLRaV-3 variant groups II and VI being the most prominent as single infections and in combinations with each other and other variants. An analysis of the distribution of variants per region showed that single infections of variant groups II and VI occurred predominantly in certain regions, and were equally distributed in the red and white cultivars studied. The distribution of a recently identified, unclassified variant of GLRaV-3 (represented by isolates GH24 and GTG10) was included in the study. The overall analysis showed that infection with variant groups II and VI were the most abundant among the samples with 49.8 and 47.6 %, respectively, followed by variant group I, variants similar to isolate GH24 and variant group III with 16.2, 13.3 and 2.5 % infection, respectively. Mixed infections, representing 36 virus combinations, were found in 251 plants. The most abundant virus combination was GLRaV-3 with Grapevine virus E (GVE), found in 28 % of the plants. GLRaV-3 was the predominant virus detected in the samples with a frequency of 80 % detection, followed by GVE (57.4 %), Grapevine rupestris stem pitting-associated virus (GRSPaV) (36.8 %), Grapevine virus A (GVA) (19.3 %), Grapevine virus F (GVF) (16.25 %), Grapevine leafroll-associated virus 2 (GLRaV-2) (8.25 %), Grapevine leafroll-associated virus 1 (GLRaV-1) (1.58 % infection) and Grapevine leafroll-associated virus 4 (GLRaV-4 like) (0.6 %). Most of the plants tested were infected with multiple viruses. The complexity of virus populations detected in this study, highlights the need for detection methods able to identify all viruses and their variants in vineyards. The information generated in this study will assist in the development of reliable detection assays that will benefit the monitoring of disease spread and aid in the efficient management of Grapevine leafroll disease (GLD).

Accurate detection and quantitation of viruses can be beneficial to plant-virus interaction studies. In this study, three SYBR green real-time RT-PCR assays were developed to quantitate grapevine leafroll-associated virus 3 (GLRaV-3) in infected vines. Three genomic regions (ORF1a, coat protein and 3'UTR) were targeted to quantitate GLRaV-3 relative to three stably expressed reference genes (actin, GAPDH and α-tubulin). These assays were able to detect all known variant groups of GLRaV-3, including the divergent group VI, with equal efficiency. No link could be established between the concentration ratios of the different genomic regions and subgenomic RNA (sgRNA) expression. However, a significant lower virus concentration ratio for plants infected with variant group VI compared to variant group II was observed for the ORF1a, coat protein and the 3'UTR. Significant higher accumulation of the virus in the growth tip was also detected for both variant groups. The quantitation of viral genomic regions under different conditions can contribute to elucidating disease aetiology and enhance knowledge about virus ecology.

Objectives:\ud \ud To enable analysis and comparisons of different relative quantitation experiments, a web-browser application called Harbin was created that uses a quantile-based scoring system for the comparison of samples at different time points and between experiments.\ud \ud Results:\ud \ud Harbin uses the standard curve method for relative quantitation to calculate concentration ratios (CRs). To evaluate if different datasets can be combined the Harbin quantile bootstrap test is proposed. This test is more sensitive in detecting distributional differences between data sets than the Kolmogorov–Smirnov test. The utility of the test is demonstrated in a comparison of three grapevine leafroll associated virus 3 (GLRaV-3) RT-qPCR data sets.\ud \ud Conclusions:\ud \ud The quantile-based scoring system of CRs will enable the monitoring of virus titre or gene expression over different time points and be useful in other genomic applications where the combining of data sets are required.

The complete genome sequences of two newly identified genetic variants of grapevine leafroll-associated virus 3 (GLRaV-3) are described. Two isolates, GH11 (18671 nt) and GH30 (18576 nt) were sequenced and found to be less than 70 % similar to other South African GLRaV-3 variants. The genome organization of GH11 and GH30 is similar to that of previously described GLRaV-3 isolates. However, no corresponding open reading frame 2 (ORF2) could be identified. Phylogenetic analysis indicates that GH11 and GH30 cluster in a subgroup (group VI) that is separate from the previously described GLRaV-3 isolates and should be regarded as a different strain of GLRaV-3.

The evolutionary history of the exclusively grapevine (Vitis spp.) infecting, grapevine leafroll-associated virus 3 (GLRaV-3) has not been studied extensively, partly due to limited available sequence data. In this study we trace the evolutionary history of GLRaV-3, focussing on isolate GH24, a newly discovered variant. GH24 was discovered through the use of next-generation sequencing (NGS) and the whole genome sequence determined and validated with Sanger sequencing. We assembled an alignment of all 13 available whole genomes of GLRaV-3 isolates and all other publicly available GLRaV-3 sequence data. Using multiple recombination detection methods we identified a clear signal for recombination in one whole genome sequence and further evidence for recombination in two more, including GH24. We inferred phylogenetic trees and networks and estimated the ages of common ancestors of GLRaV-3 clades by means of relaxed clock models calibrated with asynchronous sampling dates. Our results generally confirm previously identified variant groups as well as two new groups (VII and VIII). Higher order groups were defined as supergroups designated A to D. Supergroup A includes variant groups I-V and supergroup B group VI and its related unclassified isolates. Supergroups C and D are less well known, including the newly identified groups VII (including isolate GH24) and VIII respectively. The inferred node ages suggest that the origins of the major groups of GLRaV-3, including isolate GH24, may have occurred prior to worldwide cultivation of grapevines, whilst the current diversity represents closely related isolates that diverged from common ancestors within the last century.

Background The use of next-generation sequencing has become an established method for virus detection. Efficient study design for accurate detection relies on the optimal amount of data representing a significant portion of a virus genome. Findings In this study, genome coverage at different sequencing depths was determined for a number of viruses, viroids, hosts and sequencing library types, using both read-mapping and de novo assembly-based approaches. The results highlighted the strength of ribo-depleted RNA and sRNA in obtaining saturated genome coverage with the least amount of data, while even though the poly(A)-selected RNA yielded virus-derived reads, it was insufficient to cover the complete genome of a non-polyadenylated virus. The ribo-depleted RNA data also outperformed the sRNA data in terms of the percentage of coverage that could be obtained particularly with the de novo assembled contigs. Conclusion Our results suggest the use of ribo-depleted RNA in a de novo assembly-based approach for the detection of single-stranded RNA viruses. Furthermore, we suggest that sequencing one million reads will provide sufficient genome coverage specifically for closterovirus detection. Electronic supplementary material The online version of this article (doi:10.1186/s12985-016-0539-x) contains supplementary material, which is available to authorized users.

Background Grapevine leafroll-associated virus 3 (GLRaV-3) is the main contributing agent of leafroll disease worldwide. Four of the six GLRaV-3 variant groups known have been found in South Africa, but their individual contribution to leafroll disease is unknown. In order to study the pathogenesis of leafroll disease, a sensitive and accurate diagnostic assay is required that can detect different variant groups of GLRaV-3. Methods In this study, a one-step real-time RT-PCR, followed by high-resolution melting (HRM) curve analysis for the simultaneous detection and identification of GLRaV-3 variants of groups I, II, III and VI, was developed. A melting point confidence interval for each variant group was calculated to include at least 90% of all melting points observed. A multiplex RT-PCR protocol was developed to these four variant groups in order to assess the efficacy of the real-time RT-PCR HRM assay. Results A universal primer set for GLRaV-3 targeting the heat shock protein 70 homologue (Hsp70h) gene of GLRaV-3 was designed that is able to detect GLRaV-3 variant groups I, II, III and VI and differentiate between them with high-resolution melting curve analysis. The real-time RT-PCR HRM and the multiplex RT-PCR were optimized using 121 GLRaV-3 positive samples. Due to a considerable variation in melting profile observed within each GLRaV-3 group, a confidence interval of above 90% was calculated for each variant group, based on the range and distribution of melting points. The intervals of groups I and II could not be distinguished and a 95% joint confidence interval was calculated for simultaneous detection of group I and II variants. An additional primer pair targeting GLRaV-3 ORF1a was developed that can be used in a subsequent real-time RT-PCR HRM to differentiate between variants of groups I and II. Additionally, the multiplex RT-PCR successfully validated 94.64% of the infections detected with the real-time RT-PCR HRM. Conclusion The real-time RT-PCR HRM provides a sensitive, automated and rapid tool to detect and differentiate different variant groups in order to study the epidemiology of leafroll disease.