Quality Assessment and Validation of High-Throughput Sequencing for Grapevine Virus Diagnostics

Validation of High-Throughput Sequencing (HTS) for Routine Detection of Citrus Viruses and Viroids

The management of plant diseases relies on the accurate identification of pathogens that requires a robust and validated tool in terms of specificity, sensitivity, repeatability, and reproducibility. High-throughput sequencing (HTS) has become the method of choice for virus detection when either a complete viral status of a plant is required in a single assay or if an unknown viral agent is expected. To ensure that the most accurate diagnosis is made from an HTS data analysis, a standardized protocol per pathosystem is required. This chapter presents a detailed protocol for the detection of viruses and viroids infecting citrus using HTS. The protocol describes all the steps from sample processing, nucleic acid extraction, and bioinformatic analyses validated to be an efficient method for detection in this pathosystem. The protocol also includes a section on citrus tristeza virus (CTV) genotype differentiation using HTS data.

Cont-ID: detection of sample cross-contamination in viral metagenomic data

BACKGROUND

High-throughput sequencing (HTS) technologies completed by the bioinformatic analysis of the generated data are becoming an important detection technique for virus diagnostics. They have the potential to replace or complement the current PCR-based methods thanks to their improved inclusivity and analytical sensitivity, as well as their overall good repeatability and reproducibility. Cross-contamination is a well-known phenomenon in molecular diagnostics and corresponds to the exchange of genetic material between samples. Cross-contamination management was a key drawback during the development of PCR-based detection and is now adequately monitored in routine diagnostics. HTS technologies are facing similar difficulties due to their very high analytical sensitivity. As a single viral read could be detected in millions of sequencing reads, it is mandatory to fix a detection threshold that will be informed by estimated cross-contamination. Cross-contamination monitoring should therefore be a priority when detecting viruses by HTS technologies.

RESULTS

We present Cont-ID, a bioinformatic tool designed to check for cross-contamination by analysing the relative abundance of virus sequencing reads identified in sequence metagenomic datasets and their duplication between samples. It can be applied when the samples in a sequencing batch have been processed in parallel in the laboratory and with at least one specific external control called Alien control. Using 273 real datasets, including 68 virus species from different hosts (fruit tree, plant, human) and several library preparation protocols (Ribodepleted total RNA, small RNA and double-stranded RNA), we demonstrated that Cont-ID classifies with high accuracy (91%) viral species detection into (true) infection or (cross) contamination. This classification raises confidence in the detection and facilitates the downstream interpretation and confirmation of the results by prioritising the virus detections that should be confirmed.

CONCLUSIONS

Cross-contamination between samples when detecting viruses using HTS (Illumina technology) can be monitored and highlighted by Cont-ID (provided an alien control is present). Cont-ID is based on a flexible methodology relying on the output of bioinformatics analyses of the sequencing reads and considering the contamination pattern specific to each batch of samples. The Cont-ID method is adaptable so that each laboratory can optimise it before its validation and routine use.

G, Rivera Y, Mavrodieva VA, Nakhla MK, Roy A. High-throughput sequencing application in the detection and discovery of viruses associated with the regulated citrus leprosis disease complex

Citrus leprosis (CiL) is one of the destructive emerging viral diseases of citrus in the Americas. Leprosis syndrome is associated with two taxonomically distinct groups of Brevipalpus-transmitted viruses (BTVs), that consist of positive-sense Cilevirus, Higrevirus, and negative-sense Dichorhavirus. The localized CiL symptoms observed in multiple citrus species and other alternate hosts indicates that these viruses might have originated from the mites and eventually adopted citrus as a secondary host. Genetic diversity in the genomes of viruses associated with the CiL disease complex have complicated current detection and diagnostic measures that prompted the application of High-Throughput Sequencing (HTS) protocols for improved detection and diagnosis. Two cileviruses are known to infect citrus, and among them only citrus leprosis virus C2 (CiLV-C2) hibiscus strain (CiLV-C2H) has been reported in hibiscus and passion fruit in the US. Based on our current CiL disease complex hypothesis, there is a high probability that CiL disease is associated with more viruses/strains that have not yet been identified but exist in nature. To protect the citrus industry, a Ribo-Zero HTS protocol was utilized for detection of cileviruses infecting three different hosts: Citrus spp., Swinglea glutinosa, and Hibiscus rosa-sinensis. Real-time RT-PCR assays were used to identify plants infected with CiLV-C2 or CiLV-C2H or both in mixed infection in all the above-mentioned plant genera. These results were further confirmed by bioinformatic analysis using HTS generated data. In this study, we utilized HTS assay in confirmatory diagnostics to screen BTVs infecting Dieffenbachia sp. (family: Araceae), Passiflora edulis (Passifloraceae), and Smilax auriculata (Smilacaceae). Through the implementation of HTS and downstream data analysis, we detected not only the known cileviruses in the studied hosts but also discovered a new strain of CiLV-C2 in hibiscus from Colombia. Phylogenetically, the new hibiscus strain is more closely related to CiLV-C2 than the known hibiscus strain, CiLV-C2H. We propose this strain to be named as CiLV-C2 hibiscus strain 2 (CiLV-C2H2). The findings from the study are critical for citrus growers, industry, regulators, and researchers. The possible movement of CiLV-C2H2 from hibiscus to citrus by the Brevipalpus spp. warrants further investigation.

Virion-Associated Nucleic Acid-Based Metagenomics: A Decade of Advances in Molecular Characterization of Plant Viruses

Over the last decade, viral metagenomic studies have resulted in the discovery of thousands of previously unknown viruses. These studies are likely to play a pivotal role in obtaining an accurate and robust understanding of how viruses affect the stability and productivity of ecosystems. Among the metagenomics-based approaches that have been developed since the beginning of the 21st century, shotgun metagenomics applied specifically to virion-associated nucleic acids (VANA) has been used to disentangle the diversity of the viral world. We summarize herein the results of 24 VANA-based studies, focusing on plant and insect samples conducted over the last decade (2010 to 2020). Collectively, viruses from 85 different families were reliably detected in these studies, including capsidless RNA viruses that replicate in fungi, oomycetes, and plants. Finally, strengths and weaknesses of the VANA approach are summarized and perspectives of applications in detection, epidemiological surveillance, environmental monitoring, and ecology of plant viruses are provided. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

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.

First Report of Watermelon Mosaic Virus Naturally Infecting Coriander (Coriandrum sativum) and causing a leaf mottling disease in California

Watermelon mosaic virus (WMV, genus Potyvirus, family Potyviridae) is a species of considerable economic importance to cucurbit crops worldwide (Keinath et al. 2017). This virus has a wide host range that includes more than 170 plant species from 27 families (Dong et al. 2017; Lecoq et al. 2011). In 2018, leaves of coriander (Coriandrum sativum) plants in a student garden (C-SG) at UC Davis, and in a home garden in Davis, CA (C-Pet) (~1.1 miles apart) showed symptoms of light green mottling and crumpling. Symptomatic leaves from each location were weakly positive with the general potyvirus immunostrip test (Agdia, Elkhart, IN). In RT-PCR tests with total RNA extracts (RNeasy Plant Mini Kit Qiagen, Germantown, MD) of these leaves and the potyvirus degenerate primer pair CIFor/CIRev (Ha et al. 2008), the expected-size ~0.7 kb fragment was amplified. These fragments were gel-purified and sequenced, and a BLASTn search revealed highest identities of 91.6% (C-SG) and 97.9% (C-Pet) with the sequence of an isolate of WMV from watermelon in the U.S. (TX29, KU246036). Thus, these isolates are designated WMV-C-SG-18 and WMV-C-Pet-18. Mechanical inoculation experiments were next performed with sap prepared with symptomatic coriander leaf tissue in ice-cold 0.01 M phosphate buffer (pH 7.0) in a 1:4 wt/vol ratio. First, to obtain pure isolates, sap was inoculated onto celite-dusted leaves of Chenopodium quinoa plants (3-4 leaf stage). As expected for WMV, leaves inoculated with sap of each isolate developed chlorotic local lesions ~9 d post-inoculation (dpi) (Moreno et al. 2004). One lesion for each isolate was excised, ground in phosphate buffer, and the sap was mechanically inoculated onto leaves of Nicotiana benthamiana plants. By ~14 dpi, newly emerged leaves showed mild mottling and crumpling, and were weakly positive with the potyvirus immunostrip test. To confirm that these plants were only infected with WMV, total RNA was extracted from symptomatic leaves and used for high throughput sequencing (HTS) (Soltani et al. 2021) at the Foundation Plant Services at UC Davis. The HTS analyses revealed infection with only WMV, i.e., no other viral contigs were identified, and allowed for determination of the complete sequences (~10,000 nt) of WMV [US-CA-C-SG-18] and WMV [US-CA-C-Pet-18] with GenBank accession numbers: OM746964 and OM746965, respectively. Whole genome sequence comparisons revealed that the sequences are 99.0% identical, and 97.3% identical to the sequence of WMV TX29. Sap from symptomatic N. benthamiana leaves infected with each isolate was mechanical inoculated onto leaves of coriander plants (30-35 d old). Newly emerged leaves developed epinasty, crumpling and light green mottling by 14 dpi, and WMV infection was confirmed by RT-PCR with the WMV-specific primer pair WMV-UNI-1F and WMV-UNI-1R (Kim et al. 2019). Thus, Koch's postulates were fulfilled for this leaf mottling disease of coriander. Furthermore, the isolates from coriander induced stunting and distortion and mosaic in leaves of melon, pumpkin and squash plants by 7 dpi, whereas watermelon plants developed stunting and small leaves with mild mottling by 20 dpi. Similar results were obtained with sap prepared from infected coriander leaves. Thus, infected coriander plants are a potential inoculum source for cucurbits via several aphid vectors (Keinath et al. 2017). This is the first report of a mottle disease of coriander caused by WMV, and adds to the wide host range of the virus.

Reproducibility and Sensitivity of High-Throughput Sequencing (HTS)-Based Detection of Citrus Tristeza Virus and Three Citrus Viroids

The credibility of a pathogen detection assay is measured using specific parameters including repeatability, specificity, sensitivity, and reproducibility. The use of high-throughput sequencing (HTS) as a routine detection assay for viruses and viroids in citrus was previously evaluated and, in this study, the reproducibility and sensitivity of the HTS assay were assessed. To evaluate the reproducibility of HTS, the same plants assayed in a previous study were sampled again, one year later, and assessed in triplicate using the same analyses to construct the virome profile. The sensitivity of the HTS assay was compared to routinely used RT-PCR assays in a time course experiment, to compensate for natural pathogen accumulation in plants over time. The HTS pipeline applied in this study produced reproducible and comparable results to standard RT-PCR assays for the detection of CTV and three viroid species in citrus. Even though the limit of detection of HTS can be influenced by pathogen concentration, sample processing method and sequencing depth, detection with HTS was found to be either equivalent or more sensitive than RT-PCR in this study.

Grapevine virome and production of healthy plants by somatic embryogenesis

Grapevine (Vitis spp.) is a widespread fruit tree hosting many viral entities that interact with the plant modifying its responses to the environment. The production of virus‐free plants is becoming increasingly crucial for the use of grapevine as a model species in different studies. Using high‐throughput RNA sequencing, the viromes of seven mother plants grown in a germplasm collection vineyard were sequenced. In addition to the viruses and viroids already detected in grapevine, we identified 13 putative new mycoviruses. The different spread among grapevine tissues collected in vineyard, greenhouse and in vitro conditions suggested a clear distinction between viruses/viroids and mycoviruses that can successfully be exploited for their identification. Mycoviruses were absent in in vitro cultures, while plant viruses and viroids were particularly accumulated in these plantlets. Somatic embryogenesis applied to the seven mother plants was effective in the elimination of the complete virome, including mycoviruses. However, different sanitization efficiencies for viroids and grapevine pinot gris virus were observed among genotypes. The absence of mycoviruses in in vitro plantlets, associated with the absence of all viral entities in somaclones, suggested that this regeneration technique is also effective to eradicate endophytic/epiphytic fungi, resulting in gnotobiotic or pseudo‐gnotobiotic plants.

Detection of Banana Mild Mosaic Virus in Musa In Vitro Plants: High-Throughput Sequencing Presents Higher Diagnostic Sensitivity Than (IC)-RT-PCR and Identifies a New Betaflexiviridae Species

The banana mild mosaic virus (BanMMV) (Betaflexiviridae, Quinvirinae, unassigned species) is a filamentous virus that infects Musa spp. and has a very wide geographical distribution. The current BanMMV indexing process for an accession requires the testing of no less than four plants cultivated in a greenhouse for at least 6 months and causes a significant delay for the distribution of the germplasm. We evaluated the sensitivity of different protocols for BanMMV detection from in vitro plants to accelerate the testing process. We first used corm tissues from 137 in vitro plants and obtained a diagnostic sensitivity (DSE) of only 61% when testing four plants per accession. After thermotherapy was carried out to eliminate BanMMV infection, the meristem was recovered and further grown in vitro. The same protocol was evaluated in parallel on the corm tissue surrounding the meristem, as a rapid screening to evaluate virus therapy success, and was compared to the results obtained following the standard protocol. The obtained results showed 28% false negatives when conducting testing from corm tissues, making this protocol unsuitable in routine processes. Furthermore, RT-PCR and high-throughput sequencing (HTS) tests were applied on tissues from the base (n = 39) and the leaves (n = 36). For RT-PCR, the average DSE per sample reached 65% from either the base or leaves. HTS was applied on 36 samples and yielded 100% diagnostic specificity (DSP) and 100% DSE, whatever the sampled tissue, allowing the identification of a new Betaflexiviridae species infecting Musa. These results suggest that a reliable diagnostic of BanMMV from in vitro plants using RT-PCR or HTS technologies might represent an efficient alternative for testing after greenhouse cultivation.

Evaluation of sensitivity and specificity in RNA-Seq-based detection of grapevine viral pathogens

Virus detection is a crucial step for the implementation of clean stock programs that preserve healthy crop species. Viral infections in grapevine, a vegetatively propagated perennial crop, cannot be eradicated from the vineyards by the application of agrochemicals and must be curtailed at the stage of nursery production during the propagation of planting material. Viral detection is routinely performed using enzyme-linked immunosorbent assays (ELISA) or Reverse Transcription-quantitative Polymerase Chain Reactions (RT-qPCR). High throughput sequencing (HTS) approaches have the potential to detect all viral pathogens in a plant specimen. However, to date, no published HTS-based study has used threshold selection based on ROC curves for discriminating positive from negative samples. To fill this gap, we assessed the specificity and sensitivity of different sequencing and bioinformatics approaches for nine common viruses, which were tested in the same specimens using ELISA and/or RT-qPCR. The normalized detection thresholds giving the best results were 19.28 Fragments Per Kilobase of transcript per Million mapped reads (FPKM) for alignment-based total RNA-Seq approaches, 386 Reads Per Million mapped reads (RPM) for metagenomics-based total RNA-Seq, 1572 FPKM for alignment-based small RNA-Seq analysis and 0.97 % of contigs for de novo analysis of small RNA-Seq data. Validation of the proposed thresholds using independent specimens collected over time from the same stocks and other specimens collected from nearby stocks that had derived from the same propagating material showed that HTS approaches are accurate, with RNA-Seq approaches showing better performance than small RNA-Seq.