Comparison of High Throughput Sequencing to Standard Protocols for Virus Detection in Berry Crops

Comparative Performance Evaluation of Double-Stranded RNA High-Throughput Sequencing for the Detection of Viral Infection in Temperate Fruit Crops

There is limited information on the compared performances of biological, serological. and molecular assays with high-throughput sequencing (HTS) for viral indexing in temperate fruit crops. Here, using a range of samples of predetermined virological status, we compared two performance criteria (inclusivity and analytical sensitivity) of enzyme-linked immunosorbent assay (ELISA), molecular hybridization, reverse transcription (RT)-PCR, and double-stranded RNA (dsRNA) HTS for the detection of a total of 14 viruses (10 genera) and four viroids (three genera). When undiluted samples from individual plants were used, ELISA had the lowest performance, with an overall detection rate of 68.7%, followed by RT-PCR (82.5%) and HTS (90.7%; 100% if considering only viruses). The lower performance of RT-PCR reflected the inability to amplify some isolates as a consequence of point mutations affecting primer-binding sites. In addition, HTS identified viruses that had not been identified by other assays in nearly two-thirds of the samples. Analysis of serial dilutions of fruit tree samples allowed comparison of analytical sensitivities for various viruses. ELISA showed the lowest analytical sensitivity, but RT-PCR showed higher analytical sensitivity than HTS for most of the samples. Overall, these results confirm the superiority of HTS over biological indexing in terms of speed and inclusivity and show that while the absolute analytical sensitivity of RT-PCR tends to be higher than that of HTS, PCR inclusivity is affected by viral genetic diversity. Taken together, these results make a strong case for the implementation of HTS-based approaches in fruit tree viral testing protocols supporting quarantine and certification programs.

Complete genome sequence of vaccinium-associated virus C, a new member of the family Totiviridae from Vaccinium floribundum

Blueberries (Vaccinium sp.) are a major crop grown in the Pacific Northwest region. Currently, there are at least 17 known viruses that infect blueberry plants, and some of them cause a wide range of symptoms and economic losses. A new virus, vaccinium-associated virus C (VaVC) (family Totiviridae, genus Totivirus) was identified in an imported blueberry accession from the USDA-ARS National Clonal Germplasm Repository in Corvallis, Oregon. The complete genomic sequence of VaVC was determined, but the biological significance of VaVC is unknown and requires further study. Additional Vaccinium sp. accessions should be screened to investigate the incidence of this new virus.

Integration of Rubus yellow net virus in the raspberry genome: A story centuries in the making

Rubus yellow net virus (RYNV) belongs to genus Badnavirus. Badnaviruses are found in plants as endogenous, inactive or activatable sequences, and/or in episomal (infectious and active) forms. To assess the state of RYNV in Rubus germplasm, we sequenced the genomes of various cultivars and mined eight raspberry whole genome datasets. Bioinformatics analysis revealed the presence of a diverse array of endogenous RYNV (endoRYNV) sequences that differ significantly in their structure; some lineages have nearly complete, yet non-functional genomes whereas others have rudimentary, short sequence fragments. We developed assays to genotype the main lineages as well as the only known episomal lineage present in the United States. This study discloses the widespread presence of endoRYNVs in commercial raspberries, likely because breeding efforts have focused on a limited pool of germplasm that harbored endoRYNVs.

Phyllocoptes parviflori is a distinct species and a vector of the pervasive blackberry leaf mottle associated virus

Several viruses are transmitted by eriophyid mites (Acariformes: Eriophyoidea) including blackberry leaf mottle-associated emaravirus (BLMaV) (Emaravirus rubi). BLMaV is transmitted by an unidentified eriophyid species and is involved in blackberry yellow vein, a devastating disease in the southeastern United States. In this study, we assessed the eriophyid mite Phylocoptes parviflori as a vector of BLMaV and clarified its taxonomic status as it was previously synonymized with Phyllocoptes gracilis. P. parviflori can efficiently transmit BLMaV. The virus was found to cause yellow vein disease symptoms on 'Ouachita' blackberry marking a paradigm shift as disease symptoms have always been associated with multiple virus infections. Therefore, we propose renaming the virus to blackberry leaf mottle virus. The occurrence of P. parviflori on wild and cultivated blackberries, as well as its ability to colonize other Rubus species, enhances its importance as a major contributor to the spread of yellow vein disease.

Development of Rapid and Affordable Virus-Mimicking Artificial Positive Controls

A major bottleneck in the development of detection assays is the availability of positive controls. Their acquisition can be problematic, their maintenance is expensive, and without them, assays cannot be validated. Herein, we present a novel strategy for the development of virus-mimicking artificial positive controls (ViMAPCs). The time between design and application is less than 5 days, unlike alternatives which normally take several weeks to obtain and implement. The ViMAPCs provide a realistic representation of natural infection unlike alternatives and allow for an effortless recognition of laboratory-based contamination. The feasibility and adaptability of the strategy was evaluated using several RNA and DNA plant viruses. ViMAPCs can be used in diagnostics laboratories but also in the monitoring of pathogen outbreaks where rapid response is of utmost importance.

Molecular Characterization of a Novel Enamovirus Infecting Raspberry

Raspberry plants, valued for their fruits, are vulnerable to a range of viruses that adversely affect their yield and quality. Utilizing high-throughput sequencing (HTS), we identified a novel virus, tentatively named raspberry enamovirus 1 (RaEV1), in three distinct raspberry plants. This study provides a comprehensive characterization of RaEV1, focusing on its genomic structure, phylogeny, and possible transmission routes. Analysis of nearly complete genomes from 14 RaEV1 isolates highlighted regions of variance, particularly marked by indel events. The evidence from phylogenetic and sequence analyses supports the classification of RaEV1 as a distinct species within the Enamovirus genus. Among the 289 plant and 168 invertebrate samples analyzed, RaEV1 was detected in 10.4% and 0.4%, respectively. Most detections occurred in plants that were also infected with other common raspberry viruses. The virus was present in both commercial and wild raspberries, indicating the potential of wild plants to act as viral reservoirs. Experiments involving aphids as potential vectors demonstrated their ability to acquire RaEV1 but not to successfully transmit it to plants.

Going beyond consensus genome sequences: An innovative SNP-based methodology reconstructs different Ugandan cassava brown streak virus haplotypes at a nationwide scale in Rwanda

Cassava Brown Streak Disease (CBSD), which is caused by cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV), represents one of the most devastating threats to cassava production in Africa, including in Rwanda where a dramatic epidemic in 2014 dropped cassava yield from 3.3 million to 900,000 tonnes (1). Studying viral genetic diversity at the genome level is essential in disease management, as it can provide valuable information on the origin and dynamics of epidemic events. To fill the current lack of genome-based diversity studies of UCBSV, we performed a nationwide survey of cassava ipomovirus genomic sequences in Rwanda by high-throughput sequencing (HTS) of pools of plants sampled from 130 cassava fields in thirteen cassava-producing districts, spanning seven agro-ecological zones with contrasting climatic conditions and different cassava cultivars. HTS allowed the assembly of a nearly complete consensus genome of UCBSV in twelve districts. The phylogenetic analysis revealed high homology between UCBSV genome sequences, with a maximum of 0.8 per cent divergence between genomes at the nucleotide level. An in-depth investigation based on Single Nucleotide Polymorphisms (SNPs) was conducted to explore the genome diversity beyond the consensus sequences. First, to ensure the validity of the result, a panel of SNPs was confirmed by independent reverse transcription polymerase chain reaction (RT-PCR) and Sanger sequencing. Furthermore, the combination of fixation index (FST) calculation and Principal Component Analysis (PCA) based on SNP patterns identified three different UCBSV haplotypes geographically clustered. The haplotype 2 (H2) was restricted to the central regions, where the NAROCAS 1 cultivar is predominantly farmed. RT-PCR and Sanger sequencing of individual NAROCAS1 plants confirmed their association with H2. Haplotype 1 was widely spread, with a 100 per cent occurrence in the Eastern region, while Haplotype 3 was only found in the Western region. These haplotypes' associations with specific cultivars or regions would need further confirmation. Our results prove that a much more complex picture of genetic diversity can be deciphered beyond the consensus sequences, with practical implications on virus epidemiology, evolution, and disease management. Our methodology proposes a high-resolution analysis of genome diversity beyond the consensus between and within samples. It can be used at various scales, from individual plants to pooled samples of virus-infected plants. Our findings also showed how subtle genetic differences could be informative on the potential impact of agricultural practices, as the presence and frequency of a virus haplotype could be correlated with the dissemination and adoption of improved cultivars.

Analysis of Virus-Derived siRNAs in Strawberry Plants Co-Infected with Multiple Viruses and Their Genotypes

Plants can be infected with multiple viruses. High-throughput sequencing tools have enabled numerous discoveries of multi-strain infections, when more than one viral strain or divergent genomic variant infects a single plant. Here, we investigated small interfering RNAs (siRNAs) in a single strawberry plant co-infected with several strains of strawberry mottle virus (SMoV), strawberry crinkle virus (SCV) and strawberry virus 1 (StrV-1). A range of plants infected with subsets of the initial viral species and strains that were obtained by aphid-mediated transmission were also evaluated. Using high-throughput sequencing, we characterized the small RNA fractions associated with different genotypes of these three viruses and determined small RNA hotspot regions in viral genomes. A comparison of virus-specific siRNA (vsiRNA) abundance with relative viral concentrations did not reveal any consistent agreement. Strawberry mottle virus strains exhibiting considerable variations in concentrations were found to be associated with comparable quantities of vsiRNAs. Additionally, by estimating the specificity of siRNAs to different viral strains, we observed that a substantial pool of vsiRNAs could target all SMoV strains, while strain-specific vsiRNAs predominantly targeted rhabdoviruses, SCV and StrV-1. This highlights the intricate nature and potential interference of the antiviral response within a single infected plant when multiple viruses are present.

Area Wide Monitoring of Plant and Honey Bee (Apis mellifera) Viruses in Blueberry (Vaccinium corymbosum) Agroecosystems Facilitated by Honey Bee Pollination

Healthy agroecosystems are dependent on a complex web of factors and inter-species interactions. Flowers are hubs for pathogen transmission, including the horizontal or vertical transmission of plant-viruses and the horizontal transmission of bee-viruses. Pollination by the European honey bee (Apis mellifera) is critical for industrial fruit production, but bees can also vector viruses and other pathogens between individuals. Here, we utilized commercial honey bee pollination services in blueberry (Vaccinium corymbosum) farms for a metagenomics-based bee and plant virus monitoring system. Following RNA sequencing, viruses were identified by mapping reads to a reference sequence database through the bioinformatics portal Virtool. In total, 29 unique plant viral species were found at two blueberry farms in British Columbia (BC). Nine viruses were identified at one site in Ontario (ON), five of which were not identified in BC. Ilarviruses blueberry shock virus (BlShV) and prune dwarf virus (PDV) were the most frequently detected viruses in BC but absent in ON, while nepoviruses tomato ringspot virus and tobacco ringspot virus were common in ON but absent in BC. BlShV coat protein (CP) nucleotide sequences were nearly identical in all samples, while PDV CP sequences were more diverse, suggesting multiple strains of PDV circulating at this site. Ten bee-infecting viruses were identified, with black queen cell virus frequently detected in ON and BC. Area-wide bee-mediated pathogen monitoring can provide new insights into the diversity of viruses present in, and the health of, bee-pollination ecosystems. This approach can be limited by a short sampling season, biased towards pollen-transmitted viruses, and the plant material collected by bees can be very diverse. This can obscure the origin of some viruses, but bee-mediated virus monitoring can be an effective preliminary monitoring approach.

A Novel, Divergent Member of the Rhabdoviridae Family Infects Strawberry

Strawberry (Fragaria × ananassa) is the most important berry crop worldwide and viruses pose a constant threat to the industry. In this communication, we describe a novel virus in the family Rhabdoviridae referred to as strawberry virus 3 (StrV-3). The virus does not show significant homology when compared with recognized rhabdoviruses and, therefore, the establishment of a new genus should be considered. A triplex reverse-transcription PCR test was developed and successfully employed in a survey of the National Clonal Germplasm Repository Fragaria collection. A CRISPR-Cas-based protocol was also developed and shown to detect the virus in as little as 1 fg of total RNA, a protocol to be used in the detection of the virus in candidate G1 plants. The strawberry aphid (Chaetosiphon fragaefolii) was evaluated-alas, unsuccessfully-as a potential vector of the virus. This work broadens our understanding of the family Rhabdoviridae and assists in the quest of releasing plant material free of viruses.

Genomic Analysis and Development of Infectious Clone of a Novel Carlavirus Infecting Blueberry

A new blueberry virus was discovered using high-throughput sequencing. Using sequence identity values, phylogenetics, and serological and biological properties, we propose the virus, putatively named blueberry virus S (BluVS), to be a distinct species within the genus Carlavirus (family Betaflexiviridae). The genome was analyzed in depth, and an infectious clone was developed to initiate studies on virus pathogenicity. Agroinfiltration of the binary vector construct produced severe systemic symptoms in Nicotiana occidentalis. Back-inoculation using sap from agroinfiltrated N. occidentalis produced identical symptoms to the recipient plants (N. occidentalis), and virus purification yielded flexuous carlavirus-like particles. However, unlike blueberry scorch virus (BlScV), BluVS caused symptomless infection in Chenopodium quinoa and reacted weakly to BlScV antibodies in an enzyme-linked immunosorbent assay. Collectively, the results provide evidence for the distinct speciation of BluVS. The availability of an infectious clone provides tools for future studies on the biology of the virus.

Viroscope: Plant viral diagnosis from high-throughput sequencing data using biologically-informed genome assembly coverage

High-throughput sequencing (HTS) methods are transforming our capacity to detect pathogens and perform disease diagnosis. Although sequencing advances have enabled accessible and point-of-care HTS, data analysis pipelines have yet to provide robust tools for precise and certain diagnosis, particularly in cases of low sequencing coverage. Lack of standardized metrics and harmonized detection thresholds confound the problem further, impeding the adoption and implementation of these solutions in real-world applications. In this work, we tackle these issues and propose biologically-informed viral genome assembly coverage as a method to improve diagnostic certainty. We use the identification of viral replicases, an essential function of viral life cycles, to define genome coverage thresholds in which biological functions can be described. We validate the analysis pipeline, Viroscope, using field samples, synthetic and published datasets, and demonstrate that it provides sensitive and specific viral detection. Furthermore, we developed Viroscope.io a web-service to provide on-demand HTS data viral diagnosis to facilitate adoption and implementation by phytosanitary agencies to enable precise viral diagnosis.