De novo reconstruction of consensus master genomes of plant RNA and DNA viruses from siRNAs

Diagnosis and Characterization of Plant Viruses Using HTS to Support Virus Management and Tomato Breeding

Viral diseases pose a significant threat to tomato crops (Solanum lycopersicum L.), one of the world's most economically important vegetable crops. The limited genetic diversity of cultivated tomatoes contributes to their high susceptibility to viral infections. To address this challenge, tomato breeding programs must harness the genetic resources found in native populations and wild relatives. Breeding efforts may aim to develop broad-spectrum resistance against the virome. To identify the viruses naturally infecting 19 advanced lines, derived from native tomatoes, high-throughput sequencing (HTS) of small RNAs and confirmation with PCR and RT-PCR were used. Single and mixed infections with tomato mosaic virus (ToMV), tomato golden mosaic virus (ToGMoV), and pepper huasteco yellow vein virus (PHYVV) were detected. The complete consensus genomes of three variants of Mexican ToMV isolates were reconstructed, potentially forming a new ToMV clade with a distinct 3' UTR. The absence of reported mutations associated with resistance-breaking to ToMV suggests that the Tm-1, Tm-2, and Tm-22 genes could theoretically be used to confer resistance. However, the high mutation rates and a 63 nucleotide insertion in the 3' UTR, as well as amino acid mutations in the ORFs encoding 126 KDa, 183 KDa, and MP of Mexican ToMV isolates, suggest that it is necessary to evaluate the capacity of these variants to overcome Tm-1, Tm-2, and Tm-22 resistance genes. This evaluation, along with the characterization of advanced lines using molecular markers linked to these resistant genes, will be addressed in future studies as part of the breeding strategy. This study emphasizes the importance of using HTS for accurate identification and characterization of plant viruses that naturally infect tomato germplasm based on the consensus genome sequences. This study provides crucial insights to select appropriate disease management strategies and resistance genes and guide breeding efforts toward the development of virus-resistant tomato varieties.

Application of High-Throughput Sequencing for Comprehensive Virome Profiling in Grapevines Shows Yellows in Iran

A comprehensive study on the whole spectrum of viruses and viroids in five Iranian grapevine cultivars was carried out using sRNA libraries prepared from phloem tissue. A comparison of two approaches to virus detection from sRNAome data indicated a significant difference in the results and performance of the aligners in viral genome reconstruction. The results showed a complex virome in terms of viral composition, abundance, and richness. Thirteen viruses and viroids were identified in five Iranian grapevine cultivars, among which the grapevine red blotch virus and grapevine satellite virus were detected for the first time in Iranian vineyards. Grapevine leafroll-associated virus 1 (GLRaV1) and grapevine fanleaf virus (GFLV) were highly dominant in the virome. However, their frequency and abundance were somewhat different among grapevine cultivars. The results revealed a mixed infection of GLRaV1/grapevine yellow speckle viroid 1 (GYSVd1) and GFLV/GYSVd1 in grapevines that exhibited yellows and vein banding. We also propose a threshold of 14% of complete reconstruction as an appropriate threshold for detection of grapevine viruses that can be used as indicators for reliable grapevine virome profiling or in quarantine stations and certification programs.

Molecular Characterization of Three Apple Geminivirus Isolates in Crabapples Detected in Inner Mongolia, China

Apple geminivirus 1 (AGV) in the genus Maldovirus of the family Geminiviridae was first identified infecting apple trees in the year 2015 in China. In this work, we characterized three isolates of the AGV in the Chinese pearleaf crabapple (Malus asiatica) in Inner Mongolia Autonomous Region. The viruses were detected by Illumina sequencing and its existence was confirmed by reverse transcription-polymerase chain reaction (RT-PCR) amplification of an AGV fragment. Between the three AGV isolates and the initially characterized AGV isolate PL2015, the nucleotide sequence identities of the complete genome ranged from 91.2 to 91.7%, of the coat protein gene (V1) ranged from 95.4% to 97.3%, and of the replicase gene (C1) ranged from 87.3% to 88.0%. Phylogenetic analysis indicated that the three isolates formed a monophyletic group together with the AGV, separated from the current genera in the family Geminiviridae. This is the first description of the AGV infecting crabapples.

Virus Elimination from Naturally Infected Field Cultivars of Potato (Solanum tuberosum) by Transgenic RNA Interference

Tissue culture methods enable virus elimination from vegetatively propagated crop plants but cannot prevent new infections. Here we used a tissue culture transgenic approach for curing field cultivars of Solanum tuberosum through the stimulation of RNA interference (RNAi)-based antiviral defenses. Expression cassettes carrying inverted repeats of potato virus S (PVS, genus Carlavirus) movement or coat protein sequences were used for the transformation of potato cultivars naturally infected with PVS and/or a related carlavirus potato virus M (PVM), without or with potato virus Y (PVY, genus Potyvirus). A high proportion of transformants PCR-positive for transgenes were cured from both carlaviruses and PVY. After 3-year field trials, 22 transgenic lines representing seven cultivars remained free of any virus or became infected only with PVY. Vegetative progenies of the transgenic lines of cultivar Zeren (initially coinfected with PVS, PVM, and PVY), sampled after in vitro propagation or field trials, and other field cultivars accumulated transgene-derived 21, 22, and 24 nt small interfering (si)RNAs almost exclusively from the PVS inverted repeats. Additionally, some field progenies accumulated 21–22 nt siRNAs from the entire PVY genome, confirming PVY infection. Taken together, transgenic RNAi is effective for virus elimination from naturally infected potato cultivars and their sequence-specific immunization against new infections.

In Silico Methods for the Identification of Viral-Derived Small Interfering RNAs (vsiRNAs) and Their Application in Plant Genomics

The current era of high-throughput sequencing (HTS) technology has expedited the detection and diagnosis of viruses and viroids in the living system including plants. HTS data has become vital to study the etiology of the infection caused by both known as well as novel viral elements in planta, and their impact on overall crop health and productivity. Viral-derived small interfering RNAs are generated as a result of defence response by the host via RNAi machinery. They are immensely exploited for performing exhaustive viral investigations in plants using bioinformatics as well as experimental approaches.This chapter briefly presents the basics of virus-derived small interfering RNAs (vsiRNAs ) biology in plants and their applications in plant genomics and highlights in silico strategies exploited for virus/viroid detection. It gives a systematic pipeline for vsiRNAs identification using currently available bioinformatics tools and databases. This will surely work as a quick beginner's recipe for the in silico revelation of plant vsiRNAs as well as virus/viroid diagnosis using high-throughput sequencing data.

Grapevine Red Blotch Disease Etiology and Its Impact on Grapevine Physiology and Berry and Wine Composition

Grapevine red blotch virus (GRBV) has become widespread in the United States since its identification in 2012. GRBV is the causative agent of grapevine red blotch disease (GRBD), which has caused detrimental economic impacts to the grape and wine industry. Understanding viral function, plant–pathogen interactions, and the effects of GRBV on grapevine performance remains essential to developing potential mitigation strategies. This comprehensive review examines the current body of knowledge regarding GRBV, to highlight gaps in the knowledge and potential mitigation strategies for grape growers and winemakers.

Integrating CRISPR-Cas and Next Generation Sequencing in Plant Virology

Plant pathology has been revolutionized by the emergence and intervention of next-generation sequencing technologies (NGS) which provide a fast, cost-effective, and reliable diagnostic for any class of pathogens. NGS has made tremendous advancements in the area of research and diagnostics of plant infecting viromes and has bridged plant virology with other advanced research fields like genome editing technologies. NGS in a broader perspective holds the potential for plant health improvement by diagnosing and mitigating the new or unusual symptoms caused by novel/unidentified viruses. CRISPR-based genome editing technologies can enable rapid engineering of efficient viral/viroid resistance by directly targeting specific nucleotide sites of plant viruses and viroids. Critical genes such as eIf (iso) 4E or eIF4E have been targeted via the CRISPR platform to produce plants resistant to single-stranded RNA (ssRNA) viruses. CRISPR/Cas-based multi-target DNA or RNA tests can be used for rapid and accurate diagnostic assays for plant viruses and viroids. Integrating NGS with CRISPR-based genome editing technologies may lead to a paradigm shift in combating deadly disease-causing plant viruses/viroids at the genomic level. Furthermore, the newly discovered CRISPR/Cas13 system has unprecedented potential in plant viroid diagnostics and interference. In this review, we have highlighted the application and importance of sequencing technologies on covering the viral genomes for precise modulations. This review also provides a snapshot vision of emerging developments in NGS technologies for the characterization of plant viruses and their potential utilities, advantages, and limitations in plant viral diagnostics. Furthermore, some of the notable advances like novel virus-inducible CRISPR/Cas9 system that confers virus resistance with no off-target effects have been discussed.

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.

High-Throughput Sequencing of Small RNAs for the Sanitary Certification of Viruses in Grapevine

Biological indexing is the method generally recognized for the certification of propagative grapevines in many countries, and it is mandatory in the European Union. It consists of the evaluation of the plant material after grafting on indicators that are inspected for symptom development. This is a lengthy process that requires well-trained workers, testing field, etc. Alternative diagnostic methods such as serology and RT-qPCR have been discarded for certification because of their intrinsic drawbacks. In turn, high-throughput sequencing (HTS) of plant RNA has been proposed as a plausible alternative to bioassay, but before it is accepted, different aspects of this process must be evaluated. We have compared the HTS of small RNAs with bioassays and other diagnostic methods from a set of 40 grapevine plants submitted for certification. The results allowed the authors the identification of numerous grapevine viruses in the samples, as well as different variants. Besides, relationships between symptom expression and viromes were investigated, in particular leafroll-associated viruses. We compared HTS results using analytical and bioinformatics approaches in order to define minimum acceptable quality standards for certification schemes, resulting in a pipeline proposal. Finally, the comparison between HTS and bioassay resulted favorable for the former in terms of reliability, cost, and timing.

Integrating Viral Metagenomics into an Ecological Framework

Viral metagenomics has expanded our knowledge of the ecology of uncultured viruses, within both environmental (e.g., terrestrial and aquatic) and host-associated (e.g., plants and animals, including humans) contexts. Here, we emphasize the implementation of an ecological framework in viral metagenomic studies to address questions in virology rarely considered ecological, which can change our perception of viruses and how they interact with their surroundings. An ecological framework explicitly considers diverse variants of viruses in populations that make up communities of interacting viruses, with ecosystem-level effects. It provides a structure for the study of the diversity, distributions, dynamics, and interactions of viruses with one another, hosts, and the ecosystem, including interactions with abiotic factors. An ecological framework in viral metagenomics stands poised to broadly expand our knowledge in basic and applied virology. We highlight specific fundamental research needs to capitalize on its potential and advance the field. Expected final online publication date for the Annual Review of Virology, Volume 8 is September 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

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.

Endogenous activated small interfering RNAs in virus-infected Brassicaceae crops show a common host gene-silencing pattern affecting photosynthesis and stress response

Viral infections are accompanied by a massive production of small interfering RNAs (siRNAs) of plant origin, such as virus-activated (va)siRNAs, which drive the widespread silencing of host gene expression, and whose effects in plant pathogen interactions remain unknown. By combining phenotyping and molecular analyses, we characterized vasiRNAs that are associated with typical mosaic symptoms of cauliflower mosaic virus infection in two crops, turnip (Brassica rapa) and oilseed rape (Brassica napus), and the reference plant Arabidopsis thaliana. We identified 15 loci in the three infected plant species, whose transcripts originate vasiRNAs. These loci appear to be generally affected by virus infections in Brassicaceae and encode factors that are centrally involved in photosynthesis and stress response, such as Rubisco activase (RCA), senescence-associated protein, heat shock protein HSP70, light harvesting complex, and membrane-related protein CP5. During infection, the expression of these factors is significantly downregulated, suggesting that their silencing is a central component of the plant's response to virus infections. Further findings indicate an important role for 22 nt long vasiRNAs in the plant's endogenous RNA silencing response. Our study considerably enhances knowledge about the new class of vasiRNAs that are triggered in virus-infected plants and will help to advance strategies for the engineering of gene clusters involved in the development of crop diseases.

siRNA biogenesis and advances in topically applied dsRNA for controlling virus infections in tomato plants

A non-transgenic approach based on RNA interference was employed to induce protection against tomato mosaic virus (ToMV) infection in tomato plants. dsRNA molecules targeting the cp gene of ToMV were topically applied on plants prior to virus inoculation. Protection was dose-dependent and sequence-specific. While no protection was achieved when 0-16 µg dsRNA were used, maximum rates of resistance (60 and 63%) were observed in doses of 200 and 400 µg/plant, respectively. Similar rates were also obtained against potato virus Y when targeting its cp gene. The protection was quickly activated upon dsRNA application and lasted for up to 4 days. In contrast, no detectable antiviral response was triggered by the dsRNA from a begomovirus genome, suggesting the method is not effective against phloem-limited DNA viruses. Deep sequencing was performed to analyze the biogenesis of siRNA populations. Although long-dsRNA remained in the treated leaves for at least 10 days, its systemic movement was not observed. Conversely, dsRNA-derived siRNA populations (mainly 21- and 22-nt) were detected in non-treated leaves, which indicates endogenous processing and transport through the plant. Altogether, this study provides critical information for the development of novel tools against plant viruses; strengths and limitations inherent to the systems are discussed.

Quasi-species nature and differential gene expression of severe acute respiratory syndrome coronavirus 2 and phylogenetic analysis of a novel Iranian strain

A novel coronavirus related to severe acute respiratory syndrome virus, (SARS-CoV-2) is the causal agent of the COVID-19 pandemic. Despite the genetic mutations across the SARS-CoV-2 genome being recently investigated, its transcriptomic genetic polymorphisms at inter-host level and the viral gene expression level based on each Open Reading Frame (ORF) remains unclear. Using available High Throughput Sequencing (HTS) data and based on SARS-CoV-2 infected human transcriptomic data, this study presents a high-resolution map of SARS-CoV-2 single nucleotide polymorphism (SNP) hotspots in a viral population at inter-host level. Four throat swab samples from COVID-19 infected patients were pooled, with RNA-Seq read retrieved from SRA NCBI to detect 21 SNPs and a replacement across the SARS-CoV-2 genomic population. Twenty-two RNA modification sites on viral transcripts were identified that may cause inter-host genetic diversity of this virus. In addition, the canonical genomic RNAs of N ORF showed higher expression in transcriptomic data and reverse transcriptase quantitative PCR compared to other SARS-CoV-2 ORFs, indicating the importance of this ORF in virus replication or other major functions in virus cycle. Phylogenetic and ancestral sequence analyses based on the entire genome revealed that SARS-CoV-2 is possibly derived from a recombination event between SARS-CoV and Bat SARS-like CoV. Ancestor analysis of the isolates from different locations including Iran suggest shared Chinese ancestry. These results propose the importance of potential inter-host level genetic variations to the evolution of SARS-COV-2, and the formation of viral quasi-species. The RNA modifications discovered in this study may cause amino acid sequence changes in polyprotein, spike protein, product of ORF8 and nucleocapsid (N) protein, suggesting further insights to understanding the functional impacts of mutations in the life cycle and pathogenicity of SARS-CoV-2.

Mutations in virus-derived small RNAs

RNA viruses exist as populations of genome variants. Virus-infected plants accumulate 21–24 nucleotide small interfering RNAs (siRNAs) derived from viral RNA (virus-derived siRNAs) through gene silencing. This paper describes the profile of mutations in virus-derived siRNAs for three members of the family Potyviridae: Turnip mosaic virus (TuMV), Papaya ringspot virus (PRSV) and Wheat streak mosaic virus (WSMV). For TuMV in Arabidopsis thaliana, profiles were obtained for mechanically inoculated rosette leaves and systemically infected cauline leaves and inflorescence. Results are consistent with selection pressure on the viral genome imposed by local and systemic movement. By genetically removing gene silencing in the plant and silencing suppression in the virus, our results showed that antiviral gene silencing imposes selection in viral populations. Mutations in siRNAs derived from a PRSV coat protein transgene in the absence of virus replication showed the contribution of cellular RNA-dependent RNA polymerases to the generation of mutations in virus-derived siRNAs. Collectively, results are consistent with two sources of mutations in virus-derived siRNAs: viral RNA-dependent RNA polymerases responsible for virus replication and cellular RNA-dependent RNA polymerases responsible for gene silencing amplification.

Robust Virome Profiling and Whole Genome Reconstruction of Viruses and Viroids Enabled by Use of Available mRNA and sRNA-Seq Datasets in Grapevine (Vitis vinifera L.)

Next-generation sequencing (NGS) based virome analyses of mRNA and sRNA have recently become a routine approach for reliable detection of plant viruses and viroids. In the present study we identified the viral/viroidal spectrum of several Indian grapevine cultivars and reconstructed their whole genomes using the publically available mRNAome and sRNAome datasets. Twenty three viruses and viroids (including two variants of grapevine leafroll associated virus 4) were identified from two tissues (fruit peels and young leaves) of three cultivars among which nine unique grapevine viruses and viroids were identified for the first time in India. Irrespective of the assemblers and tissues used, the mRNA based approach identified more acellular pathogens than the sRNA based approach across cultivars. Further, the mRNAome was on par with the whole transcriptome in viral identification. Through de novo assembly of transcriptomes followed by mapping against reference genome, we reconstructed 19 complete/near complete genomes of identified viruses and viroids. The reconstructed viral genomes included four larger RNA genomes (>13 kb), a DNA genome (RG grapevine geminivirus A), a divergent genome (RG grapevine virus B) and a genome for which no reference is available (RG grapevine virus L). A large number of SNPs detected in this study ascertained the quasispecies nature of viruses. Detection of three recombination events and phylogenetic analyses using reconstructed genomes suggested the possible introduction of viruses and viroids into India from several continents through the planting material. The whole genome sequences generated in this study can serve as a resource for reliable indexing of grapevine viruses and viroids in quarantine stations and certification programs.

Whole-Genome Deep Sequencing Reveals Host-Driven in-planta Evolution of Columnea Latent Viroid (CLVd) Quasi-Species Populations

Columnea latent viroid (CLVd) is one of the most serious tomato diseases. In general, viroids have high mutation rates. This generates a population of variants (so-called quasi-species) that co-exist in their host and exhibit a huge level of genetic diversity. To study the population of CLVd in individual host plants, we used amplicon sequencing using specific CLVd primers linked with a sample-specific index sequence to amplify libraries. An infectious clone of a CLVd isolate Chaipayon-1 was inoculated on different solanaceous host plants. Six replicates of the amplicon sequencing results showed very high reproducibility. On average, we obtained 133,449 CLVd reads per PCR-replicate and 79 to 561 viroid sequence variants, depending on the plant species. We identified 19 major variants (>1.0% mean relative abundance) in which a total of 16 single-nucleotide polymorphisms (SNPs) and two single nucleotide insertions were observed. All major variants contained a combination of 4 to 6 SNPs. Secondary structure prediction clustered all major variants into a tomato/bolo maka group with four loops (I, II, IV and V), and a chili pepper group with four loops (I, III, IV and V) at the terminal right domain, compared to the CLVd Chaipayon-1 which consists of five loops (I, II, III, IV and V).

PVsiRNAPred: Prediction of plant exclusive virus-derived small interfering RNAs by deep convolutional neural network

Plant exclusive virus-derived small interfering RNAs (vsiRNAs) regulate various biological processes, especially important in antiviral immunity. The identification of plant vsiRNAs is important for understanding the biogenesis and function mechanisms of vsiRNAs and further developing anti-viral plants. In this study, we extracted plant vsiRNA sequences from the PVsiRNAdb database. We then utilized deep convolutional neural network (CNN) to develop a deep learning algorithm for predicting plant vsiRNAs based on vsiRNA sequence composition, known as PVsiRNAPred. The key part of PVsiRNAPred is the CNN module, which automatically learns hierarchical representations of vsiRNA sequences related to vsiRNA profiles in plants. When evaluated using an independent testing dataset, the accuracy of the model was 65.70%, which was higher than those of five conventional machine learning method-based classifiers. In addition, PVsiRNAPred obtained a sensitivity of 67.11%, specificity of 64.26% and Matthews correlation coefficient (MCC) of 0.31, and the area under the receiver operating characteristic (ROC) curve (AUC) of PVsiRNAPred was 0.71 in the independent test. The permutation test with 1000 shuffles resulted in a p value of<0.001. The above results reveal that PVsiRNAPred has favorable generalization capabilities. We hope PVsiRNAPred, the first bioinformatics algorithm for predicting plant vsiRNAs, will allow efficient discovery of new vsiRNAs.

Current and Emerging Methods for the Synthesis of Single-Stranded DNA

Methods for synthesizing arbitrary single-strand DNA (ssDNA) fragments are rapidly becoming fundamental tools for gene editing, DNA origami, DNA storage, and other applications. To meet the rising application requirements, numerous methods have been developed to produce ssDNA. Some approaches allow the synthesis of freely chosen user-defined ssDNA sequences to overcome the restrictions and limitations of different length, purity, and yield. In this perspective, we provide an overview of the representative ssDNA production strategies and their most significant challenges to enable the readers to make informed choices of synthesis methods and enhance the availability of increasingly inexpensive synthetic ssDNA. We also aim to stimulate a broader interest in the continued development of efficient ssDNA synthesis techniques and improve their applications in future research.

Plant virome reconstruction and antiviral RNAi characterization by deep sequencing of small RNAs from dried leaves

In plants, RNA interference (RNAi) generates small interfering (si)RNAs from entire genomes of viruses, satellites and viroids. Therefore, deep small (s)RNA sequencing is a universal approach for virome reconstruction and RNAi characterization. We tested this approach on dried barley leaves from field surveys. Illumina sequencing of sRNAs from 2 plant samples identified in both plants Hordeum vulgare endornavirus (HvEV) and barley yellow mosaic bymovirus (BaYMV) and, additionally in one plant, a novel strain of Japanese soil-borne wheat mosaic furovirus (JSBWMV). De novo and reference-based sRNA assembly yielded complete or near-complete genomic RNAs of these viruses. While plant sRNAs showed broad size distribution, viral sRNAs were predominantly 21 and 22 nucleotides long with 5′-terminal uridine or adenine, and were derived from both genomic strands. These bona fide siRNAs are presumably processed from double-stranded RNA precursors by Dicer-like (DCL) 4 and DCL2, respectively, and associated with Argonaute 1 and 2 proteins. For BaYMV (but not HvEV, or JSBWMV), 24-nucleotide sRNAs represented the third most abundant class, suggesting DCL3 contribution to anti-bymovirus defence. Thus, viral siRNAs are well preserved in dried leaf tissues and not contaminated by non-RNAi degradation products, enabling both complete virome reconstruction and inference of RNAi components mediating antiviral defense.

From current knowledge to best practice: A primer on viral diagnostics using deep sequencing of virus-derived small interfering RNAs (vsiRNAs) in infected plants

Plants have evolved many defense strategies for combating viral infections. One major surveillance strategy adopted by them is manipulating viral sequences to generate distinct small RNA products via Dicer-like enzymes (DCL), and thereby restricting virus multiplication through the RNA interference (RNAi) mechanism. The power of high-throughput sequencing technologies, with diverse computational tools to handle small RNA sequencing (sRNA-Seq) data, bestows unprecedented opportunities to answer fundamental questions in plant virology. Here, we present some basic concepts of virus-derived, small interfering RNA (vsiRNA) biogenesis in plants, optimization strategies, caveats, and best practices for efficient discovery and diagnosis of known as well as novel plant viruses/viroids using deep sequencing of small RNA (sRNA) pools.

Analyses of virus/viroid communities in nectarine trees by next-generation sequencing and insight into viral synergisms implication in host disease symptoms

We analyzed virus and viroid communities in five individual trees of two nectarine cultivars with different disease phenotypes using next-generation sequencing technology. Different viral communities were found in different cultivars and individual trees. A total of eight viruses and one viroid in five families were identified in a single tree. To our knowledge, this is the first report showing that the most-frequently identified viral and viroid species co-infect a single individual peach tree, and is also the first report of peach virus D infecting Prunus in China. Combining analyses of genetic variation and sRNA data for co-infecting viruses/viroid in individual trees revealed for the first time that viral synergisms involving a few virus genera in the Betaflexiviridae, Closteroviridae, and Luteoviridae families play a role in determining disease symptoms. Evolutionary analysis of one of the most dominant peach pathogens, peach latent mosaic viroid (PLMVd), shows that the PLMVd sequences recovered from symptomatic and asymptomatic nectarine leaves did not all cluster together, and intra-isolate divergent sequence variants co-infected individual trees. Our study provides insight into the role that mixed viral/viroid communities infecting nectarine play in host symptom development, and will be important in further studies of epidemiological features of host-pathogen interactions.

High-Throughput Sequencing Analysis of Small RNAs Derived from Coleus Blumei Viroids

Characterization of viroid-derived small RNAs (vd-sRNAs) is important to understand viroid–host interactions; however, vd-sRNAs belonging to the genus Coleviroid are yet to be identified and characterized. Herein, we used coleus plants singly infected with coleus blumei viroid (CbVd)-1, -5, or -6 and doubly infected with CbVd-1 and -5 to identify and analyze their vd-sRNAs. We found sense and antisense vd-sRNAs for CbVd-1, -5 and -6, and 22-nt vd-sRNAs were the most abundant; moreover, the 5′-terminal nucleotides (nts) of CbVd-1, -5, and -6 were biased toward U and C, and sRNAs derived from these three viroids were unevenly distributed along their genomes. We also noted that CbVd-5 and -6 share a fragment that forms the right half of the rod-like secondary structure of these viroids, which implied that they generated almost the same type of vd-sRNAs. This finding indicated that vd-sRNA biogenesis is mainly determined by the primary sequence of their substrates. More importantly, we found two complementary vd-sRNAs (22 nt) that were generated from the central conserved region (CCR) of these three viroids, suggesting an important role of CCR in vd-sRNA biogenesis. In conclusion, our results provide novel insight into the biogenesis of vd-sRNAs and the biological roles of CCR.

Spatial Associations of Vines Infected With Grapevine Red Blotch Virus in Oregon Vineyards

Spread and in-field spatial patterns of vines infected with grapevine red blotch virus (GRBV) were documented in Oregon vineyards using field sampling, molecular diagnostics, and spatial analysis. Grapevine petiole tissue collected from 2013 to 2016 was tested using quantitative polymerase chain reaction for GRBV. At Jacksonville in southern Oregon, 3.1% of vines were infected with GRBV in 2014, and GRBV incidence reached 58.5% of study vines by 2016. GRBV-infected plants and GRBV-uninfected plants were spatially aggregated at this site in 2015, and infected plants were spatially associated between years 2015 and 2016. In a southern Oregon vineyard near Talent, 10.4% of vines were infected with GRBV in 2014, and infection increased annually to 21.5% in 2016. At Talent, distribution of the infected vines was spatially associated across all years. GRBV infection was highest at Yamhill, in the Willamette Valley, where 31.7% of the tested vines had GRBV infection in 2014. By 2016, 59.2% of the vines tested positive for GRBV. Areas of aggregation increased and were spatially associated across all years. From 2013 to 2015, GRBV was not detected at Milton-Freewater in eastern Oregon. Spatial patterns of GRBV infection support evidence of spread by a mobile insect vector. GRBV is a significant threat to Oregon wine grape production because of its drastic year-over-year spread in affected vineyards.

High Throughput Sequencing For Plant Virus Detection and Discovery

Over the last decade, virologists have discovered an unprecedented number of viruses using high throughput sequencing (HTS), which led to the advancement of our knowledge on the diversity of viruses in nature, particularly unraveling the virome of many agricultural crops. However, these new virus discoveries have often widened the gaps in our understanding of virus biology; the forefront of which is the actual role of a new virus in disease, if any. Yet, when used critically in etiological studies, HTS is a powerful tool to establish disease causality between the virus and its host. Conversely, with globalization, movement of plant material is increasingly more common and often a point of dispute between countries. HTS could potentially resolve these issues given its capacity to detect and discover. Although many pipelines are available for plant virus discovery, all share a common backbone. A description of the process of plant virus detection and discovery from HTS data are presented, providing a summary of the different pipelines available for scientists' utility in their research.

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.

Shannon Entropy to Evaluate Substitution Rate Variation Among Viral Nucleotide Positions in Datasets of Viral siRNAs

Next-generation sequencing has opened the door to the reconstruction of viral populations and examination of the composition of mutant spectra in infected cells, tissues, and host organisms. In this chapter we present details on the use of the Shannon entropy method to estimate the site-specific nucleotide relative variability of turnip crinkle virus, a positive (+) stranded RNA plant virus, in a large dataset of short RNAs of Cicer arietinum L., a natural reservoir of the virus. We propose this method as a viral metagenomics tool to provide a more detailed description of the viral quasispecies in infected plant tissue. Viral replicative fitness relates to an optimal composition of variants that provide the molecular basis of virus behavior in the complex environment of natural infections. A complete description of viral quasispecies may have implications in determining fitness landscapes for host-virus coexistence and help to design specific diagnostic protocols and antiviral strategies.

Next-Generation Sequencing and CRISPR/Cas13 Editing in Viroid Research and Molecular Diagnostics

Viroid discovery as well as the economic significance of viroids and biological properties are presented. Next-generation sequencing (NGS) technologies combined with informatics have been applied to viroid research and diagnostics for almost a decade. NGS provides highly efficient, rapid, low-cost high-throughput sequencing of viroid genomes and of the 21⁻24 nt vd-sRNAs generated by the RNA silencing defense of the host. NGS has been utilized in various viroid studies which are presented. The discovery during the last few years that prokaryotes have heritable adaptive immunity mediated through clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated Cas proteins, have led to transformative advances in molecular biology, notably genome engineering and most recently molecular diagnostics. The potential application of the CRISPR-Cas13a system for engineering viroid interference in plants is suggested by targeting specific motifs of three economically important viroids. The CRISPR-Cas13 system has been utilized recently for the accurate detection of human RNA viruses by visual read out in 90 min or less and by paper-based assay. Multitarget RNA tests by this technology have a good potential for application as a rapid and accurate diagnostic assay for known viroids. The CRISPR/Cas system will work only for known viroids in contrast to NGS, but it should be much faster.

Small RNA-Omics for Plant Virus Identification, Virome Reconstruction, and Antiviral Defense Characterization

RNA interference (RNAi)-based antiviral defense generates small interfering RNAs that represent the entire genome sequences of both RNA and DNA viruses as well as viroids and viral satellites. Therefore, deep sequencing and bioinformatics analysis of small RNA population (small RNA-ome) allows not only for universal virus detection and genome reconstruction but also for complete virome reconstruction in mixed infections. Viral infections (like other stress factors) can also perturb the RNAi and gene silencing pathways regulating endogenous gene expression and repressing transposons and host genome-integrated endogenous viral elements which can potentially be released from the genome and contribute to disease. This review describes the application of small RNA-omics for virus detection, virome reconstruction and antiviral defense characterization in cultivated and non-cultivated plants. Reviewing available evidence from a large and ever growing number of studies of naturally or experimentally infected hosts revealed that all families of land plant viruses, their satellites and viroids spawn characteristic small RNAs which can be assembled into contigs of sufficient length for virus, satellite or viroid identification and for exhaustive reconstruction of complex viromes. Moreover, the small RNA size, polarity and hotspot profiles reflect virome interactions with the plant RNAi machinery and allow to distinguish between silent endogenous viral elements and their replicating episomal counterparts. Models for the biogenesis and functions of small interfering RNAs derived from all types of RNA and DNA viruses, satellites and viroids as well as endogenous viral elements are presented and discussed.

Foxtail mosaic virus: A Viral Vector for Protein Expression in Cereals

Rapid and cost-effective virus-derived transient expression systems for plants are invaluable in elucidating gene function and are particularly useful in plant species for which transformation-based methods are unavailable or are too time and labor demanding, such as wheat (Triticum aestivum) and maize (Zea mays). The virus-mediated overexpression (VOX) vectors based on Barley stripe mosaic virus and Wheat streak mosaic virus described previously for these species are incapable of expressing free recombinant proteins of more than 150 to 250 amino acids, are not suited for high-throughput screens, and have other limitations. In this study, we report the development of a VOX vector based on a monopartite single-stranded positive sense RNA virus, Foxtail mosaic virus (genus Potexvirus). In this vector, PV101, the gene of interest was inserted downstream of the duplicated subgenomic promoter of the viral coat protein gene, and the corresponding protein was expressed in its free form. The vector allowed the expression of a 239-amino acid-long GFP in both virus-inoculated and upper uninoculated (systemic) leaves of wheat and maize and directed the systemic expression of a larger approximately 600-amino acid protein, GUSPlus, in maize. Moreover, we demonstrated that PV101 can be used for in planta expression and functional analysis of apoplastic pathogen effector proteins such as the host-specific toxin ToxA of Parastagonospora nodorum Therefore, this VOX vector opens possibilities for functional genomics studies in two important cereal crops.

Causative Role of Grapevine Red Blotch Virus in Red Blotch Disease

Grapevine red blotch virus (GRBV) has a monopartite single-stranded DNA genome and is the type species of the genus Grablovirus in the family Geminiviridae. To address the etiological role of GRBV in the recently recognized red blotch disease of grapevine, infectious GRBV clones were engineered from the genome of each of the two previously identified phylogenetic clades for Agrobacterium tumefaciens-mediated inoculations of tissue culture-grown Vitis spp. plants. Following agroinoculation and one or two dormancy cycles, systemic GRBV infection was detected by multiplex polymerase chain reaction (PCR) in Vitis vinifera exhibiting foliar disease symptoms but not in asymptomatic vines. Infected rootstock genotype SO4 (V. berlandieri × V. riparia) exhibited leaf chlorosis and cupping, while infection was asymptomatic in agroinoculated 110R (V. berlandieri × V. rupestris), 3309C (V. riparia × V. rupestris), and V. rupestris. Spliced GRBV transcripts of the replicase-associated protein coding region accumulated in leaves of agroinfected vines, as shown by reverse-transcription PCR; this was consistent with systemic infection resulting from virus replication. Additionally, a virus progeny identical in nucleotide sequence to the infectious GRBV clones was recovered from agroinfected vines by rolling circle amplification, cloning, and sequencing. Concomitantly, subjecting naturally infected grapevines to microshoot tip culture resulted in an asymptomatic plant progeny that tested negative for GRBV in multiplex PCR. Altogether, our agroinoculation and therapeutic experiments fulfilled Koch's postulates and revealed the causative role of GRBV in red blotch disease.

Small RNA-Omics for Virome Reconstruction and Antiviral Defense Characterization in Mixed Infections of Cultivated Solanum Plants

In plants, RNA silencing-based antiviral defense generates viral small RNAs (sRNAs) faithfully representing the viral genomes. We employed sRNA sequencing and bioinformatics (sRNA-omics) to characterize antiviral defense and to reconstruct the full genomic sequences and their variants in the evolving viral quasispecies in cultivated solanaceous plants carrying mixed infections. In naturally infected Solanum tuberosum (potato), one case study revealed a virome comprising Potato virus Y (genus Potyvirus) and Potato virus X (genus Potexvirus), which was reconstructed by de novo-assembling separate genome-size sRNA contigs. Another case study revealed a virome comprising NTN and O strains of Potato virus Y, whose sRNAs assembled in chimeric contigs, which could be disentangled on the basis of reference genome sequences. Both viromes were stable in vegetative potato progeny. In a cross-protection trial of Solanum lycopersicum (tomato), the supposedly protective mild strain CH2 of Pepino mosaic virus (genus Potexvirus) was tested for protection against strain LP of the same virus. Reciprocal mechanical inoculations eventually resulted in co-infection of all individual plants with CH2 and LP strains, reconstructed as separate sRNA contigs. LP invasions into CH2-preinfected plants and vice versa were accompanied by alterations of consensus genome sequences in viral quasispecies, indicating a potential risk of cross-protection measures. Additionally, the study also revealed, by reconstruction from sRNAs, the presence of the mechanically nontransmissible Southern tomato virus (genus Amalgavirus) in some plants. Our in-depth analysis of sRNA sizes, 5'-nucleotide frequencies and hotspot maps revealed similarities in sRNA-generating mechanisms in potato and tomato, differential silencing responses to virome components and potential for sRNA-directed cross-targeting between viral strains which could not, however, prevent the formation of stable viromes.

Analysis of DNAs associated with coconut foliar decay disease implicates a unique single-stranded DNA virus representing a new taxon

The unique ecology, pathology and undefined taxonomy of coconut foliar decay virus (CFDV), found associated with coconut foliar decay disease (CFD) in 1986, prompted analyses of old virus samples by modern methods. Rolling circle amplification and deep sequencing applied to nucleic acid extracts from virion preparations and CFD-affected palms identified twelve distinct circular DNAs, eleven of which had a size of about 1.3 kb and one of 641 nt. Mass spectrometry-based protein identification proved that a 24 kDa protein encoded by two 1.3 kb DNAs is the virus capsid protein with highest sequence similarity to that of grabloviruses (family Geminiviridae), even though CFDV particles are not geminate. The nine other 1.3 kb DNAs represent alphasatellites coding for replication initiator proteins that differ clearly from those encoded by nanovirid DNA-R. The 641 nt DNA-gamma is unique and may encode a movement protein. Three DNAs, alphasatellite CFDAR, capsid protein encoding CFDV DNA-S.1 and DNA-gamma share sequence motifs near their replication origins and were consistently present in all samples analysed. These DNAs appear to be integral components of a possibly tripartite CFDV genome, different from those of any Geminiviridae or Nanoviridae family member, implicating CFDV as representative of a new genus and family.

Gene expression and population polymorphism of maize Iranian mosaic virus in Zea mays, and intracellular localization and interactions of viral N, P, and M proteins in Nicotiana benthamiana

Maize Iranian mosaic virus (MIMV; Mononegavirales, Rhabdoviridae, Nucleorhabdovirus) infects maize and several other poaceous plants. MIMV encodes six proteins, i.e., nucleocapsid protein (N), polymerase cofactor phosphoprotein (P), putative movement protein (P3), matrix protein (M), glycoprotein (G), and large RNA-dependent RNA polymerase (L). In the present study, MIMV gene expression and genetic polymorphism of an MIMV population in maize were determined. N, P, P3, and M protein genes were more highly expressed than the 5' terminal G and L genes. Twelve single nucleotide polymorphisms were identified across the genome within a MIMV population in maize from RNA-Seq read data pooled from three infected plants indicating genomic variations of potential importance to evolution of the virus. MIMV N, P, and M proteins that are known to be involved in rhabdovirus replication and transcription were characterized as to their intracellular localization and interactions. N protein accumulated exclusively in the nucleus and interacted with itself and with P protein. P protein accumulated in both the nucleus and cell periphery and interacted with itself, N and M proteins in the nucleus. M protein was localized in the cell periphery and on endomembranes, and interacted with P protein in the nucleus. MIMV proteins show a distinctive combination of intracellular localizations and interactions.

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.

Revisiting the Roles of Tobamovirus Replicase Complex Proteins in Viral Replication and Silencing Suppression

Tobamoviral replicase possesses an RNA-dependent RNA polymerase (RDR) domain and is translated from genomic (g)RNA via a stop codon readthrough mechanism at a one-to-ten ratio relative to a shorter protein lacking the RDR domain. The two proteins share methyltransferase and helicase domains and form a heterodimer implicated in gRNA replication. The shorter protein is also implicated in suppressing RNA silencing-based antiviral defenses. Using a stop codon mutant of Oilseed rape mosaic tobamovirus (ORMV), we demonstrate that the readthrough replicase (p182) is sufficient for gRNA replication and for subgenomic RNA transcription during systemic infection in Nicotiana benthamiana and Arabidopsis thaliana. However, the mutant virus displays milder symptoms and does not interfere with HEN1-mediated methylation of viral short interfering (si)RNAs or plant small (s)RNAs. The mutant virus tends to revert the stop codon, thereby restoring expression of the shorter protein (p125), even in the absence of plant Dicer-like activities that generate viral siRNAs. Plant RDR activities that generate endogenous siRNA precursors do not prevent replication or movement of the mutant virus, and double-stranded precursors of viral siRNAs representing the entire virus genome are likely synthesized by p182. Transgenic expression of p125 partially recapitulates the ORMV disease symptoms associated with overaccumulation of plant sRNAs. Taken together, the readthrough replicase p182 is sufficient for viral replication and transcription but not for silencing suppression. By contrast, the shorter p125 protein suppresses silencing, provokes severe disease symptoms, causes overaccumulation of unmethylated viral and plant sRNAs but it is not an essential component of the viral replicase complex.

Identification and characterization of privet leaf blotch-associated virus, a novel idaeovirus

A novel virus has been identified by next-generation sequencing (NGS) in privet (Ligustrum japonicum L.) affected by a graft-transmissible disease characterized by leaf blotch symptoms resembling infectious variegation, a virus-like privet disease with an unclear aetiology. This virus, which has been tentatively named 'privet leaf blotch-associated virus' (PrLBaV), was absent in non-symptomatic privet plants, as revealed by NGS and reverse transcription-polymerase chain reaction (RT-PCR). Molecular characterization of PrLBaV showed that it has a segmented genome composed of two positive single-stranded RNAs, one of which (RNA1) is monocistronic and codes for the viral replicase, whereas the other (RNA2) contains two open reading frames (ORFs), ORF2a and ORF2b, coding for the putative movement (p38) and coat (p30) proteins, respectively. ORF2b is very probably expressed through a subgenomic RNA starting with six nucleotides (AUAUCU) that closely resemble those found in the 5'-terminal end of genomic RNA1 and RNA2 (AUAUUU and AUAUAU, respectively). The molecular signatures identified in the PrLBaV RNAs and proteins resemble those of Raspberry bushy dwarf virus (RBDV), currently the only member of the genus Idaeovirus. These data, together with phylogenetic analyses, are consistent with the proposal of considering PrLBaV as a representative of the second species in the genus Idaeovirus. Transient expression of a recombinant PrLBaV p38 fused to green fluorescent protein in leaves of Nicotiana benthamiana, coupled with confocal laser scanning microscopy assays, showed that it localizes at cell plasmodesmata, strongly supporting its involvement in viral movement/trafficking and providing the first functional characterization of an idaeovirus encoded protein.

RNAi-mediated resistance to viruses: a critical assessment of methodologies

In plants, RNA interference (RNAi)-based antiviral defense is mediated by multigenic families of Dicer-like enzymes generating small interfering (si)RNAs from double-stranded RNA (dsRNA) produced during replication and/or transcription of RNA and DNA viruses, and Argonaute enzymes binding viral siRNAs and targeting viral RNA and DNA for siRNA-directed posttranscriptional and transcriptional silencing. Successful viruses are able to suppress or evade the production or action of viral siRNAs. In antiviral biotech approaches based on RNAi, transgenic expression or non-transgenic delivery of dsRNA cognate to a target virus pre-activates or boosts the natural plant antiviral defenses. Design of more effective antiviral RNAi strategies requires better understanding of viral siRNA biogenesis and viral anti-silencing strategies in virus-infected plants.

The pepper virome: natural co-infection of diverse viruses and their quasispecies

BACKGROUND

The co-infection of diverse viruses in a host plant is common; however, little is known about viral populations and their quasispecies in the host.

RESULTS

Here, we report the first pepper viromes that were co-infected by different types of viral genomes. The pepper viromes are dominated by geminivirus DNA-A followed by a novel carlavirus referred to as Pepper virus A. The two pepper cultivars share similar viral populations and replications. However, the quasispecies for double-stranded RNA virus and two satellite DNAs were heterogeneous and homogenous in susceptible and resistant cultivars, respectively, indicating the quasispecies of an individual virus depends on the host.

CONCLUSIONS

Taken together, we provide the first evidence that the host plant resistant to viruses has an unrevealed antiviral system, affecting viral quasispecies, not replication.

Viroid quasispecies revealed by deep sequencing

Viroids are non-coding single-stranded circular RNA molecules that replicate autonomously in infected host plants causing mild to lethal symptoms. Their genomes contain about 250-400 nucleotides, depending on viroid species. Members of the family Pospiviroidae, like the Potato spindle tuber viroid (PSTVd), replicate via an asymmetric rolling-circle mechanism using the host DNA-dependent RNA-Polymerase II in the nucleus, while members of Avsunviroidae are replicated in a symmetric rolling-circle mechanism probably by the nuclear-encoded polymerase in chloroplasts. Viroids induce the production of viroid-specific small RNAs (vsRNA) that can direct (post-)transcriptional gene silencing against host transcripts or genomic sequences. Here, we used deep-sequencing to analyze vsRNAs from plants infected with different PSTVd variants to elucidate the PSTVd quasipecies evolved during infection. We recovered several novel as well as previously known PSTVd variants that were obviously competent in replication and identified common strand-specific mutations. The calculated mean error rate per nucleotide position was less than [Formula: see text], quite comparable to the value of [Formula: see text] reported for a member of Avsunviroidae. The resulting error threshold allows the synthesis of longer-than-unit-length replication intermediates as required by the asymmetric rolling-circle mechanism of members of Pospiviroidae.

Viral small RNAs reveal the genomic variations of three grapevine vein clearing virus quasispecies populations

Viral small RNAs (vsRNAs) include viral small interfering RNAs (vsiRNAs) that are initiators and products of RNA silencing, and small RNAs that are derived from viral RNAs with function still unknown. Sequencing of vsRNAs allows assembling of viral genomes and revelation of viral population variations at genomic levels. Grapevine vein clearing virus (GVCV) is a new member of the family Caulimoviridae whose DNA genome is replicated by reverse transcription of pre-genomic RNA molecules. In this short report, three genomic sequences of GVCV were assembled from vsRNAs that were isolated and sequenced from three individual grapevines in commercial vineyards and compared to the GVCV-CHA reference genome. Profiles of single nucleotide polymorphism among three viral populations indicated a closer relatedness between two populations in different grape cultivars at the same location than those in the same grape cultivar at different locations, suggesting the spread of GVCV populations among vineyards of close proximity. Classic types of vsiRNAs (21-nt, 22-nt, and 24-nt) were found in the three GVCV vsiRNA populations, but these did not produce alignment hotspots on the GVCV-CHA reference genome. The number of 36-nt reads is the highest among vsRNAs, the role of these vsRNAs remains unclear. The analysis of vsRNAs provides a first holistic picture of genomic variations among GVCV viral quasispecies populations that help monitor epidemics and evolution of GVCV populations, an emerging virus that is becoming a threat to grape production in the Midwest region of the USA.

Red blotch disease alters grape berry development and metabolism by interfering with the transcriptional and hormonal regulation of ripening

Grapevine red blotch-associated virus (GRBaV) is a major threat to the wine industry in the USA. GRBaV infections (aka red blotch disease) compromise crop yield and berry chemical composition, affecting the flavor and aroma properties of must and wine. In this study, we combined genome-wide transcriptional profiling with targeted metabolite analyses and biochemical assays to characterize the impact of the disease on red-skinned berry ripening and metabolism. Using naturally infected berries collected from two vineyards, we were able to identify consistent berry responses to GRBaV across different environmental and cultural conditions. Specific alterations of both primary and secondary metabolism occurred in GRBaV-infected berries during ripening. Notably, GRBaV infections of post-véraison berries resulted in the induction of primary metabolic pathways normally associated with early berry development (e.g. thylakoid electron transfer and the Calvin cycle), while inhibiting ripening-associated pathways, such as a reduced metabolic flux in the central and peripheral phenylpropanoid pathways. We show that this metabolic reprogramming correlates with perturbations at multiple regulatory levels of berry development. Red blotch caused the abnormal expression of transcription factors (e.g. NACs, MYBs, and AP2-ERFs) and elements of the post-transcriptional machinery that function during red-skinned berry ripening. Abscisic acid, ethylene, and auxin pathways, which control both the initiation of ripening and stress responses, were also compromised. We conclude that GRBaV infections disrupt normal berry development and stress responses by altering transcription factors and hormone networks, which result in the inhibition of ripening pathways involved in the generation of color, flavor, and aroma compounds.

Metavisitor, a Suite of Galaxy Tools for Simple and Rapid Detection and Discovery of Viruses in Deep Sequence Data

Metavisitor is a software package that allows biologists and clinicians without specialized bioinformatics expertise to detect and assemble viral genomes from deep sequence datasets. The package is composed of a set of modular bioinformatic tools and workflows that are implemented in the Galaxy framework. Using the graphical Galaxy workflow editor, users with minimal computational skills can use existing Metavisitor workflows or adapt them to suit specific needs by adding or modifying analysis modules. Metavisitor works with DNA, RNA or small RNA sequencing data over a range of read lengths and can use a combination of de novo and guided approaches to assemble genomes from sequencing reads. We show that the software has the potential for quick diagnosis as well as discovery of viruses from a vast array of organisms. Importantly, we provide here executable Metavisitor use cases, which increase the accessibility and transparency of the software, ultimately enabling biologists or clinicians to focus on biological or medical questions.

Deep sequencing for discovery and evolutionary analysis of plant viruses

The advent of next generation sequencing (NGS), or deep sequencing, has allowed great advances to be made in discovery, diagnostics, and evolutionary studies in plant viruses. Various methods have been used for enrichment for virus-specific nucleic acids, each of which have some drawbacks. Many novel viruses have been discovered in plants by NGS technologies, and there is a good deal of promise for more comprehensive studies in virus evolution. However, each aspect of using NGS has its caveats that need to be considered, and there is still a need for better tools of analysis, as well as method for validation of sequence variation.

Synergies and antagonisms in virus interactions

Metagenomic surveys and data from next generation sequencing revealed that mixed infections among plant viruses are probably a rule rather than an exception in natural pathosystems. The documented cases may range from synergism to antagonism, which may depend from the spatiotemporal order of arrival of the viruses on the host and upon the host itself. In synergistic interactions, the measurable differences in replication, phenotypic and cytopathological changes, cellular tropism, within host movement, and transmission rate of one of the two viruses or both are increased. Conversely, a decrease in replication, or inhibition of one or more of the above functions by one virus against the other, leads to an antagonistic interaction. Viruses may interact directly and by transcomplementation of defective functions or indirectly, through responses mediated by the host like the defense mechanism based on RNA silencing. Outcomes of these interactions can be applied to the risk assessment of transgenic crops expressing viral proteins, or cross-protected crops for the identification of potential hazards. Prior to experimental evidence, mathematical models may help in forecasting challenges deriving from the great variety of pathways of synergistic and antagonistic interactions. Actually, it seems that such predictions do not receive sufficient credit in the framework of agriculture.

Identification and characterization of an isolate of apple green crinkle associated virus involved in a severe disease of quince (Cydonia oblonga, Mill.)

A quince tree showing severe symptoms of a previously undescribed viral disease occurring in northern Apulia (Italy) was analysed using high-throughput sequencing of small RNA libraries, leading to the identification of a new strain of apple green crinkle associated virus (isolate AGCaV-CYD) showing peculiar traits. RT-PCR with specific primers detected AGCaV-CYD in consistent association with symptoms in the surveyed orchards. Molecular characterization of the reconstructed genome, together with phylogenetic analysis, showed it to be closely related to an AGCaV strain causing green crinkle disease in apple (AGCaV-AUR) and divergent from the type strain of apple stem pitting virus (ASPV-PA66).

Phylogeny of Geminivirus Coat Protein Sequences and Digital PCR Aid in Identifying Spissistilus festinus as a Vector of Grapevine red blotch-associated virus

Grapevine red blotch-associated virus (GRBaV) is a single-stranded DNA virus, proposed to be a member of the family Geminiviridae, and is associated with grapevines showing red blotch symptoms in North America. The existence of the virus was reported in 2012, and subsequently detected in grapevines in major grape production regions. We investigated if a vector exists that can transmit GRBaV in vineyards. Phylogenetic analysis of the predicted amino acid sequence of coat protein (CP) of GRBaV with the CP of 23 geminiviruses representing all seven genera of the family Geminiviridae revealed that GRBaV-CP was most similar to that of Tomato pseudo-curly top virus, a geminivirus known to be transmitted by a treehopper (Membracidae), a family that is closely related to leafhoppers (Cicadellidae). To identify vectors of GRBaV, hemipteran species within and nearby wine grape vineyards where virus spread was suspected were collected and transmission assays were conducted. Among the species tested, the three-cornered alfalfa hopper Spissistilus festinus (Hemiptera: Membracidae) was able to both acquire the virus from a grapevine infected with GRBaV and transmit the virus to healthy grapevines in the laboratory. In commercial vineyards, lateral shoots of grapevines girdled due to feeding injury by the adult three-cornered alfalfa hopper also tested positive for the virus using digital PCR. These findings represent an important step in understanding the biology of GRBaV and develop management guidelines.

Next-Generation Sequencing and Genome Editing in Plant Virology

Next-generation sequencing (NGS) has been applied to plant virology since 2009. NGS provides highly efficient, rapid, low cost DNA, or RNA high-throughput sequencing of the genomes of plant viruses and viroids and of the specific small RNAs generated during the infection process. These small RNAs, which cover frequently the whole genome of the infectious agent, are 21-24 nt long and are known as vsRNAs for viruses and vd-sRNAs for viroids. NGS has been used in a number of studies in plant virology including, but not limited to, discovery of novel viruses and viroids as well as detection and identification of those pathogens already known, analysis of genome diversity and evolution, and study of pathogen epidemiology. The genome engineering editing method, clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 system has been successfully used recently to engineer resistance to DNA geminiviruses (family, Geminiviridae) by targeting different viral genome sequences in infected Nicotiana benthamiana or Arabidopsis plants. The DNA viruses targeted include tomato yellow leaf curl virus and merremia mosaic virus (begomovirus); beet curly top virus and beet severe curly top virus (curtovirus); and bean yellow dwarf virus (mastrevirus). The technique has also been used against the RNA viruses zucchini yellow mosaic virus, papaya ringspot virus and turnip mosaic virus (potyvirus) and cucumber vein yellowing virus (ipomovirus, family, Potyviridae) by targeting the translation initiation genes eIF4E in cucumber or Arabidopsis plants. From these recent advances of major importance, it is expected that NGS and CRISPR-Cas technologies will play a significant role in the very near future in advancing the field of plant virology and connecting it with other related fields of biology.

Integrated analyses using RNA-Seq data reveal viral genomes, single nucleotide variations, the phylogenetic relationship, and recombination for Apple stem grooving virus

Background

Next-generation sequencing (NGS) provides many possibilities for plant virology research. In this study, we performed integrated analyses using plant transcriptome data for plant virus identification using Apple stem grooving virus (ASGV) as an exemplar virus. We used 15 publicly available transcriptome libraries from three different studies, two mRNA-Seq studies and a small RNA-Seq study.

Results

We de novo assembled nearly complete genomes of ASGV isolates Fuji and Cuiguan from apple and pear transcriptomes, respectively, and identified single nucleotide variations (SNVs) of ASGV within the transcriptomes. We demonstrated the application of NGS raw data to confirm viral infections in the plant transcriptomes. In addition, we compared the usability of two de novo assemblers, Trinity and Velvet, for virus identification and genome assembly. A phylogenetic tree revealed that ASGV and Citrus tatter leaf virus (CTLV) are the same virus, which was divided into two clades. Recombination analyses identified six recombination events from 21 viral genomes.

Conclusions

Taken together, our in silico analyses using NGS data provide a successful application of plant transcriptomes to reveal extensive information associated with viral genome assembly, SNVs, phylogenetic relationships, and genetic recombination.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2994-6) contains supplementary material, which is available to authorized users.

Characterization of a Novel Polerovirus Infecting Maize in China

A novel virus, tentatively named Maize Yellow Mosaic Virus (MaYMV), was identified from the field-grown maize plants showing yellow mosaic symptoms on the leaves collected from the Yunnan Province of China by the deep sequencing of small RNAs. The complete 5642 nucleotide (nt)-long genome of the MaYMV shared the highest nucleotide sequence identity (73%) to Maize Yellow Dwarf Virus-RMV. Sequence comparisons and phylogenetic analyses suggested that MaYMV represents a new member of the genus Polerovirus in the family Luteoviridae. Furthermore, the P0 protein encoded by MaYMV was demonstrated to inhibit both local and systemic RNA silencing by co-infiltration assays using transgenic Nicotiana benthamiana line 16c carrying the GFP reporter gene, which further supported the identification of a new polerovirus. The biologically-active cDNA clone of MaYMV was generated by inserting the full-length cDNA of MaYMV into the binary vector pCB301. RT-PCR and Northern blot analyses showed that this clone was systemically infectious upon agro-inoculation into N. benthamiana. Subsequently, 13 different isolates of MaYMV from field-grown maize plants in different geographical locations of Yunnan and Guizhou provinces of China were sequenced. Analyses of their molecular variation indicate that the 3' half of P3-P5 read-through protein coding region was the most variable, whereas the coat protein- (CP-) and movement protein- (MP-)coding regions were the most conserved.