Correction: Haegeman et al. Looking beyond Virus Detection in RNA Sequencing Data: Lessons Learned from a Community-Based Effort to Detect Cellular Plant Pathogens and Pests. Plants 2023, 12, 2139.

In the original publication [...].

Looking beyond Virus Detection in RNA Sequencing Data: Lessons Learned from a Community-Based Effort to Detect Cellular Plant Pathogens and Pests.

High-throughput sequencing (HTS), more specifically RNA sequencing of plant tissues, has become an indispensable tool for plant virologists to detect and identify plant viruses. During the data analysis step, plant virologists typically compare the obtained sequences to reference virus databases. In this way, they are neglecting sequences without homologies to viruses, which usually represent the majority of sequencing reads. We hypothesized that traces of other pathogens might be detected in this unused sequence data. In the present study, our goal was to investigate whether total RNA-seq data, as generated for plant virus detection, is also suitable for the detection of other plant pathogens and pests. As proof of concept, we first analyzed RNA-seq datasets of plant materials with confirmed infections by cellular pathogens in order to check whether these non-viral pathogens could be easily detected in the data. Next, we set up a community effort to re-analyze existing Illumina RNA-seq datasets used for virus detection to check for the potential presence of non-viral pathogens or pests. In total, 101 datasets from 15 participants derived from 51 different plant species were re-analyzed, of which 37 were selected for subsequent in-depth analyses. In 29 of the 37 selected samples (78%), we found convincing traces of non-viral plant pathogens or pests. The organisms most frequently detected in this way were fungi (15/37 datasets), followed by insects (13/37) and mites (9/37). The presence of some of the detected pathogens was confirmed by independent (q)PCRs analyses. After communicating the results, 6 out of the 15 participants indicated that they were unaware of the possible presence of these pathogens in their sample(s). All participants indicated that they would broaden the scope of their bioinformatic analyses in future studies and thus check for the presence of non-viral pathogens. In conclusion, we show that it is possible to detect non-viral pathogens or pests from total RNA-seq datasets, in this case primarily fungi, insects, and mites. With this study, we hope to raise awareness among plant virologists that their data might be useful for fellow plant pathologists in other disciplines (mycology, entomology, bacteriology) as well.

Nuances of Responses to Two Sources of Grapevine Leafroll Disease on Pinot Noir Grown in the Field for 17 Years.

Grapevine leafroll disease (GLD) is one of the most economically damaging virus diseases in grapevine, with grapevine leafroll-associated virus 1 (GLRaV-1) and grapevine leafroll-associated virus 3 (GLRaV-3) as the main contributors. This study complements a previously published transcriptomic analysis and compared the impact of two different forms of GLD to a symptomless control treatment: a mildly symptomatic form infected with GLRaV-1 and a severe form with exceptionally early leafroll symptoms (up to six weeks before veraison) infected with GLRaV-1 and GLRaV-3. Vine physiology and fruit composition in 17-year-old Pinot noir vines were measured and a gradient of vigor, yield, and berry quality (sugar content and berry weight) was observed between treatments. Virome composition, confirmed by individual RT-PCR, was compared with biological indexing. Three divergent viromes were recovered, containing between four to seven viruses and two viroids. They included the first detection of grapevine asteroid mosaic-associated virus in Switzerland. This virus did not cause obvious symptoms on the indicators used in biological indexing. Moreover, the presence of grapevine virus B (GVB) did not cause the expected corky bark symptoms on the indicators, thus underlining the important limitations of the biological indexing. Transmission of GLRaV-3 alone or in combination with GVB by Planococcus comstocki mealybug did not reproduce the strong symptoms observed on the donor plant infected with a severe form of GLD. This result raises questions about the contribution of each virus to the symptomatology of the plant.

An Advanced One-Step RT-LAMP for Rapid Detection of Little cherry virus 2 Combined with High-Throughput Sequence-Based Phylogenomics Reveal Divergent Flowering Cherry Isolates.

Little cherry virus 2 (LChV-2, genus Ampelovirus) is considered to be the main causal agent of the economically damaging little cherry disease, which can only be controlled by removal of infected trees. The widespread viral disease of sweet cherry (Prunus avium L.) is affecting the survival of long-standing orchards in North America and Europe, hence the dire need for an early and accurate diagnosis to establish a sound disease control strategy. The endemic presence of LChV-2 is mainly confirmed using laborious time-consuming reverse-transcription (RT-PCR). A rapid reverse transcription loop-mediated isothermal amplification (RT-LAMP) assay targeting a conserved region of the coat protein was developed and compared with conventional RT-PCR for the specific detection of LChV-2. This affordable assay, combined with a simple RNA extraction, deploys desirable characteristics such as higher ability for faster (<15 min), more analytically sensitive (100-fold), and robust broad-range diagnosis of LChV-2 isolates from sweet cherry, ornamental flowering cherry displaying heterogenous viral etiology and, for the first time, newly identified potential insect vectors. Moreover, use of Sanger and total RNA high-throughput sequencing as complementary metaviromics approaches confirmed the LChV-2 RT-LAMP detection of divergent LChV-2 isolates in new hosts and the relationship of their whole-genome was exhaustively inferred using maximum-likelihood phylogenomics. This entails unprecedented critical understanding of a novel evolutionary clade further expanding LChV-2 viral diversity. In conclusion, this highly effective diagnostic platform facilitates strategical support for early in-field testing to reliably prevent dissemination of new LChV-2 outbreaks from propagative plant stocks or newly postulated insect vectors. Validated results and major advantages are herein thoroughly discussed, in light of the knowledge required to increase the potential accuracy of future diagnostics and the essential epidemiological considerations to proactively safeguard cherries and Prunus horticultural crop systems from little cherry disease.

Identification and Molecular Characterization of a Novel Hordeivirus Associated With Yellow Mosaic Disease of Privet (Ligustrum vulgare) in Europe.

Wild plants serve as a large reservoir of known and yet-unknown viruses and as a source of viral pathogens of cultivated plants. Yellow mosaic disease of forest shrub Ligustrum vulgare (privet) was recurrently observed in Europe for more than 100 years. Using a universal virus identification approach based on deep sequencing and de novo assembly of viral small interfering (si)RNAs we identified a causative agent of this disease in Switzerland and reconstructed its complete 3-segmented RNA genome. Notably, a short 3'-terminal common region (CR) attached to each segment via a ∼53-71 nucleotide poly(A) tract, as determined by RT-PCR sequencing, was initially identified as an orphan siRNA contig with conserved tRNA-like secondary structure. Phylogenomic analysis classified this virus as a novel member in the genus Hordeivirus of family Virgaviridae, which we named ligustrum mosaic virus (LigMV). Similar to other hordeiviruses, LigMV formed rod-shape virions (visualized by electron microscopy), was transmitted through seeds and could also be mechanically transmitted to herbaceous hosts Chenopodium quinoa and Nicotiana benthamiana. Blot hybridization analysis identified genomic and subgenomic RNAs, sharing the 3'-CR and likely serving as monocistronic mRNAs for seven evolutionarily-conserved viral proteins including two subunits of viral RNA-dependent RNA polymerase, coat protein, triple gene block proteins mediating viral movement and cysteine-rich suppressor of RNA silencing. Analysis of size, polarity, and hotspot profiles of viral siRNAs suggested that they are produced by the plant antiviral Dicer-like (DCL) proteins DCL2 and DCL4 processing double-stranded intermediates of genomic RNA replication. Whole genome sequencing of French and Austrian isolates of LigMV revealed its genetic stability over a wide geographic range (>99% nucleotide identity to Swiss isolates and each other), suggesting its persistence and spread in Europe via seed dispersal.

A Primer on the Analysis of High-Throughput Sequencing Data for Detection of Plant Viruses.

High-throughput sequencing (HTS) technologies have become indispensable tools assisting plant virus diagnostics and research thanks to their ability to detect any plant virus in a sample without prior knowledge. As HTS technologies are heavily relying on bioinformatics analysis of the huge amount of generated sequences, it is of utmost importance that researchers can rely on efficient and reliable bioinformatic tools and can understand the principles, advantages, and disadvantages of the tools used. Here, we present a critical overview of the steps involved in HTS as employed for plant virus detection and virome characterization. We start from sample preparation and nucleic acid extraction as appropriate to the chosen HTS strategy, which is followed by basic data analysis requirements, an extensive overview of the in-depth data processing options, and taxonomic classification of viral sequences detected. By presenting the bioinformatic tools and a detailed overview of the consecutive steps that can be used to implement a well-structured HTS data analysis in an easy and accessible way, this paper is targeted at both beginners and expert scientists engaging in HTS plant virome projects.

A novel foveavirus identified in wild grapevine (Vitis vinifera subsp. sylvestris).

We report the genome sequence of a putative new foveavirus infecting non-cultivated Vitis vinifera, tentatively named "grapevine foveavirus A" (GFVA). This virus was identified by high-throughput sequencing analysis of a European wild Vitis collected in Switzerland. Phylogenetic analysis revealed that this virus clustered with known grapevine virus T (GVT) isolates but was clearly distinct from any of them. If considering the International Committee of Taxonomy of Viruses (ICTV)-suggested foveavirus species demarcation criterion based on sequence similarity in the replicase gene/protein, this virus should be considered a member of a new species closely related to GVT. On the other hand, comparison of capsid gene/protein sequences using the same criteria indicates that GFVA is at the border of species demarcation. Whether this virus represents a highly divergent GVT isolate or a member of a distinct but closely related species is discussed.

Single berry reconstitution prior to RNA-sequencing reveals novel insights into transcriptomic remodeling by leafroll virus infections in grapevines.

Leafroll viruses are among the most devastating pathogens in viticulture and are responsible for major economic losses in the wine industry. However, the molecular interactions underlying the effects on fruit quality deterioration are not well understood. The few molecular studies conducted on berries from infected vines, associated quality decreases with the repression of key genes in sugar transport and anthocyanin biosynthesis. Sampling protocols in these studies did however not account for berry heterogeneity and potential virus induced phenological shifts, which could have biased the molecular information. In the present study, we adopted an innovative individual berry sampling protocol to produce homogeneous batches for RNA extraction, thereby circumventing berry heterogeneity and compensating for virus induced phenological shifts. This way a characterization of the transcriptomic modulation by viral infections was possible and explain why our results differ significantly from previously reported repression of anthocyanin biosynthesis and sugar metabolism. The present study provides new insights into the berry transcriptome modulation by leafroll infection, highlighting the virus induced upregulation of plant innate immunity as well as an increased responsiveness of the early ripening berry to biotic stressors. The study furthermore emphasizes the importance of sampling protocols in physiological studies on grapevine berry metabolism.

Identification and Characterization of a Novel Robigovirus Species from Sweet Cherry in Turkey.

High throughput sequencing of total RNA isolated from symptomatic leaves of a sweet cherry tree (Prunus avium cv. 0900 Ziraat) from Turkey identified a new member of the genus Robigovirus designated cherry virus Turkey (CVTR). The presence of the virus was confirmed by electron microscopy and overlapping RT-PCR for sequencing its whole-genome. The virus has a ssRNA genome of 8464 nucleotides which encodes five open reading frames (ORFs) and comprises two non-coding regions, 5' UTR and 3' UTR of 97 and 296 nt, respectively. Compared to the five most closely related robigoviruses, RdRp, TGB1, TGB2, TGB3 and CP share amino acid identities ranging from 43-53%, 44-60%, 39-43%, 38-44% and 45-50%, respectively. Unlike the four cherry robigoviruses, CVTR lacks ORFs 2a and 5a. Its genome organization is therefore more similar to African oil palm ringspot virus (AOPRV). Using specific primers, the presence of CVTR was confirmed in 15 sweet cherries and two sour cherries out of 156 tested samples collected from three regions in Turkey. Among them, five samples were showing slight chlorotic symptoms on the leaves. It seems that CVTR infects cherry trees with or without eliciting obvious symptoms, but these data should be confirmed by bioassays in woody and possible herbaceous hosts in future studies.

Virus Detection by High-Throughput Sequencing of Small RNAs: Large-Scale Performance Testing of Sequence Analysis Strategies.

Recent developments in high-throughput sequencing (HTS), also called next-generation sequencing (NGS), technologies and bioinformatics have drastically changed research on viral pathogens and spurred growing interest in the field of virus diagnostics. However, the reliability of HTS-based virus detection protocols must be evaluated before adopting them for diagnostics. Many different bioinformatics algorithms aimed at detecting viruses in HTS data have been reported but little attention has been paid thus far to their sensitivity and reliability for diagnostic purposes. Therefore, we compared the ability of 21 plant virology laboratories, each employing a different bioinformatics pipeline, to detect 12 plant viruses through a double-blind large-scale performance test using 10 datasets of 21- to 24-nucleotide small RNA (sRNA) sequences from three different infected plants. The sensitivity of virus detection ranged between 35 and 100% among participants, with a marked negative effect when sequence depth decreased. The false-positive detection rate was very low and mainly related to the identification of host genome-integrated viral sequences or misinterpretation of the results. Reproducibility was high (91.6%). This work revealed the key influence of bioinformatics strategies for the sensitive detection of viruses in HTS sRNA datasets and, more specifically (i) the difficulty in detecting viral agents when they are novel or their sRNA abundance is low, (ii) the influence of key parameters at both assembly and annotation steps, (iii) the importance of completeness of reference sequence databases, and (iv) the significant level of scientific expertise needed when interpreting pipeline results. Overall, this work underlines key parameters and proposes recommendations for reliable sRNA-based detection of known and unknown viruses.

Absolute Quantification of Grapevine Red Blotch Virus in Grapevine Leaf and Petiole Tissues by Proteomics.

Grapevine red blotch is a recently identified viral disease that was first recognized in the Napa Valley of California. Infected plants showed foliar symptoms similar to leafroll, another grapevine viral disease, on vines testing negative for known grapevine leafroll-associated virus. Later, the Grapevine red blotch virus (GRBV) was independently discovered in the US states of California and New York and was demonstrated to be the causal agent of red blotch disease. Due to its wide occurrence in the United States, vector transmission, and impacts on grape industry, this virus has the potential to cause serious economic losses. Despite numerous attempts, it has yet not been possible to isolate or visualize viral particles from GRBV-infected plants, thereby hampering the development of a serological assay that would facilitate GRBV detection in grapevine. In this work, mass spectrometry approaches were applied in order to quantify GRBV in infected plants and identify potential biomarkers for viral infection. We present for the first time the physical detection on the protein level of the two GRBV genes V1 (coat protein) and V2 in grapevine tissue lysates. The GRBV coat protein load in petioles was determined to be in the range of 100-900 million copies per milligram wet weight by using three heavy isotope labeled reference peptides as internal standards. In leaves on the other hand, the V1 copy number per unit wet tissue weight appeared to be about six times lower than in petioles, and about 300 times lower in terms of protein concentration in the extractable protein mass, albeit these estimations could only be made with one reference peptide detectable in leaf extracts. Moreover, we found in leaf and petiole extracts of GRBV-infected plants a consistent upregulation of several enzymes involved in flavonoid biosynthesis by label-free shotgun proteomics, indicating the activation of a defense mechanism against GRBV, a plant response already described for Grapevine leafroll-associated virus infection on the transcriptome level. Finally and importantly, we identified some other microorganisms belonging to the grapevine leaf microbiota, two bacterial species (Novosphingobium sp. Rr 2-17 and Methylobacterium) and one virus, Grapevine rupestris stem pitting-associated virus.

Molecular Characterization of Divergent Closterovirus Isolates Infecting Ribes Species.

Five isolates of a new member of the family Closteroviridae, tentatively named blackcurrant leafroll-associated virus 1 (BcLRaV-1), were identified in the currant. The 17-kb-long genome codes for 10 putative proteins. The replication-associated polyprotein has several functional domains, including papain-like proteases, methyltransferase, Zemlya, helicase, and RNA-dependent RNA polymerase. Additional open reading frames code for a small protein predicted to integrate into the host cell wall, a heat-shock protein 70 homolog, a heat-shock protein 90 homolog, two coat proteins, and three proteins of unknown functions. Phylogenetic analysis showed that BcLRaV-1 is related to members of the genus Closterovirus, whereas recombination analysis provided evidence of intraspecies recombination.

Grapevine red blotch virus: Absence in Swiss Vineyards and Analysis of Potential Detrimental Effect on Viticultural Performance.

Grapevine red blotch virus (GRBV) is a recently described virus that infects grapevine. Little information is available on the possible occurrence and distribution outside North America. Therefore, we surveyed commercial vineyards from the three major grape-growing regions in Switzerland to determine the presence or absence of GRBV. In total, 3,062 vines were analyzed by polymerase chain reaction. None of the vines tested positive for GRBV, suggesting the absence of GRBV from Swiss vineyards. We also investigated whether GRBV was present in 653 grapevine accessions in the Agroscope grapevine virus collection at Nyon, including dominantly Swiss (457) but also international accessions. Only six referential accessions were infected by GRBV, all originating from the United States, whereas all others from 10 European and 8 non-European origins tested negative. High-throughput sequencing analysis of Zinfandel A2V13, in the collection since 1985, confirmed close similarity of GRBV isolate Z_A2V13 to American isolates according to genomes deposited in GenBank. Because the Zinfandel A2V13 reference was also maintained grafted on the leafroll virus indicator Vitis vinifera 'Gamay', we evaluated the effect of GRBV on viticultural performance over a 3-year period. Our results showed clear detrimental effects of GRBV on grapevine physiology (vine vigor, leaf chlorophyll content, and gas exchange) and fruit quality. These findings underscore the importance of implementation of GRBV testing worldwide in certification and quarantine programs to prevent the dissemination of this virus.

A Framework for the Evaluation of Biosecurity, Commercial, Regulatory, and Scientific Impacts of Plant Viruses and Viroids Identified by NGS Technologies.

Recent advances in high-throughput sequencing technologies and bioinformatics have generated huge new opportunities for discovering and diagnosing plant viruses and viroids. Plant virology has undoubtedly benefited from these new methodologies, but at the same time, faces now substantial bottlenecks, namely the biological characterization of the newly discovered viruses and the analysis of their impact at the biosecurity, commercial, regulatory, and scientific levels. This paper proposes a scaled and progressive scientific framework for efficient biological characterization and risk assessment when a previously known or a new plant virus is detected by next generation sequencing (NGS) technologies. Four case studies are also presented to illustrate the need for such a framework, and to discuss the scenarios.

Biological, Serological, and Molecular Characterization of a Highly Divergent Strain of Grapevine leafroll-associated virus 4 Causing Grapevine Leafroll Disease.

The complete genome sequence of a highly divergent strain of Grapevine leafroll-associated virus 4 (GLRaV-4) was determined using 454 pyrosequencing technology. This virus, designated GLRaV-4 Ob, was detected in Vitis vinifera 'Otcha bala' from our grapevine virus collection at Agroscope. The GLRaV-4 Ob genome length and organization share similarities with members of subgroup II in the genus Ampelovirus (family Closteroviridae). Otcha bala was graft-inoculated onto indicator plants of cultivar Gamay to evaluate the biological properties of this new strain, and typical leafroll symptoms were induced. A monoclonal antibody for the rapid detection of GLRaV-4 Ob by enzyme-linked immunosorbent assay is available, thus facilitating large-scale diagnostics of this virus. Based on the relatively small size of the coat protein, the reduced amino acid identity and the distinct serological properties, our study clearly shows that GLRaV-4 Ob is a divergent strain of GLRaV-4. Furthermore, molecular and serological data revealed that the AA42 accession from which GLRaV-7 was originally reported is in fact co-infected with GLRaV-4 Ob and GLRaV-7. This finding challenges the idea that GLRaV-7 is a leafroll-causing agent.