High Throughput Sequencing For Plant Virus Detection and Discovery

Identification, molecular characterization and prevalence of a novel cytorhabdovirus infecting zucchini crops in Greece

A new cytorhabdovirus was identified in zucchini (Cucurbita pepo) in Greece with the aid of high-throughput sequencing technology. The negative-sense, single-stranded genomic RNA of the new virus was determined and includes seven open reading frames in the order 3'-N-P-P3-P4-M-G-L-5' in the antigenomic orientation. Typical rhabdovirus-like particles were observed in infected leaf material. Comparative sequence analysis and phylogenetic reconstructions suggested that the described virus is a new member of the genus Cytorhabdovirus, and it was tentatively named cucurbit cytorhabdovirus 1 (CuCV1). To our knowledge CuCV1 is the first cytorhabdovirus infecting cucurbits in nature. Our surveys indicated that it occurs in a percentage of 36.7 % in zucchini crops in Greece.

High-Throughput Sequencing Facilitates Discovery of New Plant Viruses in Poland

Viruses cause epidemics on all major crops of agronomic importance, and a timely and accurate identification is essential for control. High throughput sequencing (HTS) is a technology that allows the identification of all viruses without prior knowledge on the targeted pathogens. In this paper, we used HTS technique for the detection and identification of different viral species occurring in single and mixed infections in plants in Poland. We analysed various host plants representing different families. Within the 20 tested samples, we identified a total of 13 different virus species, including those whose presence has not been reported in Poland before: clover yellow mosaic virus (ClYMV) and melandrium yellow fleck virus (MYFV). Due to this new finding, the obtained sequences were compared with others retrieved from GenBank. In addition, cucurbit aphid-borne yellows virus (CABYV) was also detected, and due to the recent occurrence of this virus in Poland, a phylogenetic analysis of these new isolates was performed. The analysis revealed that CABYV population is highly diverse and the Polish isolates of CABYV belong to two different phylogenetic groups. Our results showed that HTS-based technology is a valuable diagnostic tool for the identification of different virus species originating from variable hosts, and can provide rapid information about the spectrum of plant viruses previously not detected in a region.

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.

RNAseq Analysis of Rhizomania-Infected Sugar Beet Provides the First Genome Sequence of Beet Necrotic Yellow Vein Virus from the USA and Identifies a Novel Alphanecrovirus and Putative Satellite Viruses

“Rhizomania” of sugar beet is a soilborne disease complex comprised of beet necrotic yellow vein virus (BNYVV) and its plasmodiophorid vector, Polymyxa betae. Although BNYVV is considered the causal agent of rhizomania, additional viruses frequently accompany BNYVV in diseased roots. In an effort to better understand the virus cohort present in sugar beet roots exhibiting rhizomania disease symptoms, five independent RNA samples prepared from diseased beet seedlings reared in a greenhouse or from field-grown adult sugar beet plants and enriched for virus particles were subjected to RNAseq. In all but a healthy control sample, the technique was successful at identifying BNYVV and provided sequence reads of sufficient quantity and overlap to assemble > 98% of the published genome of the virus. Utilizing the derived consensus sequence of BNYVV, infectious RNA was produced from cDNA clones of RNAs 1 and 2. The approach also enabled the detection of beet soilborne mosaic virus (BSBMV), beet soilborne virus (BSBV), beet black scorch virus (BBSV), and beet virus Q (BVQ), with near-complete genome assembly afforded to BSBMV and BBSV. In one field sample, a novel virus sequence of 3682 nt was assembled with significant sequence similarity and open reading frame (ORF) organization to members within the subgenus Alphanecrovirus (genus Necrovirus; family Tombusviridae). Construction of a DNA clone based on this sequence led to the production of the novel RNA genome in vitro that was capable of inducing local lesion formation on leaves of Chenopodium quinoa. Additionally, two previously unreported satellite viruses were revealed in the study; one possessing weak similarity to satellite maize white line mosaic virus and a second possessing moderate similarity to satellite tobacco necrosis virus C. Taken together, the approach provides an efficient pipeline to characterize variation in the BNYVV genome and to document the presence of other viruses potentially associated with disease severity or the ability to overcome resistance genes used for sugar beet rhizomania disease management.

Raspberry leaf blotch emaravirus in Bosnia and Herzegovina: population structure and systemic movement

Raspberry leaf blotch virus (RLBV) is the putative agent of the homonymous disease and even though Bosnia and Herzegovina is a major producer worldwide there is no report of the virus presence in the country. We studied the virus population structure and assessed its ability to move systemically. RLBV is widespread in production areas and has a homogeneous population structure; leading to the hypothesis that the primary mode of dissemination is propagation material. The ability of the virus to move systemically eliminates propagation of root cuttings as a viable option to obtain RLBV-free plants, leaving RT-PCR screening as the better option to propagate RLBV- free plants in the absence of clean-up facilities or certification programs in the country.

Plant Viruses Infecting Solanaceae Family Members in the Cultivated and Wild Environments: A Review

Plant viruses infecting crop species are causing long-lasting economic losses and are endangering food security worldwide. Ongoing events, such as climate change, changes in agricultural practices, globalization of markets or changes in plant virus vector populations, are affecting plant virus life cycles. Because farmer’s fields are part of the larger environment, the role of wild plant species in plant virus life cycles can provide information about underlying processes during virus transmission and spread. This review focuses on the Solanaceae family, which contains thousands of species growing all around the world, including crop species, wild flora and model plants for genetic research. In a first part, we analyze various viruses infecting Solanaceae plants across the agro-ecological interface, emphasizing the important role of virus interactions between the cultivated and wild zones as global changes affect these environments on both local and global scales. To cope with these changes, it is necessary to adjust prophylactic protection measures and diagnostic methods. As illustrated in the second part, a complex virus research at the landscape level is necessary to obtain relevant data, which could be overwhelming. Based on evidence from previous studies we conclude that Solanaceae plant communities can be targeted to address complete life cycles of viruses with different life strategies within the agro-ecological interface. Data obtained from such research could then be used to improve plant protection methods by taking into consideration environmental factors that are impacting the life cycles of plant viruses.

Smash++: an alignment-free and memory-efficient tool to find genomic rearrangements

BACKGROUND

The development of high-throughput sequencing technologies and, as its result, the production of huge volumes of genomic data, has accelerated biological and medical research and discovery. Study on genomic rearrangements is crucial owing to their role in chromosomal evolution, genetic disorders, and cancer.

RESULTS

We present Smash++, an alignment-free and memory-efficient tool to find and visualize small- and large-scale genomic rearrangements between 2 DNA sequences. This computational solution extracts information contents of the 2 sequences, exploiting a data compression technique to find rearrangements. We also present Smash++ visualizer, a tool that allows the visualization of the detected rearrangements along with their self- and relative complexity, by generating an SVG (Scalable Vector Graphics) image.

CONCLUSIONS

Tested on several synthetic and real DNA sequences from bacteria, fungi, Aves, and Mammalia, the proposed tool was able to accurately find genomic rearrangements. The detected regions were in accordance with previous studies, which took alignment-based approaches or performed FISH (fluorescence in situ hybridization) analysis. The maximum peak memory usage among all experiments was ∼1 GB, which makes Smash++ feasible to run on present-day standard computers.

The population structure of Rose rosette virus in the USA

Rose rosette virus (RRV) (genus Emaravirus) is the causal agent of the homonymous disease, the most destructive malady of roses in the USA. Although the importance of the disease is recognized, little sequence information and no full genomes are available for RRV, a multi-segmented RNA virus. To better understand the population structure of the virus we implemented a Hi-Plex PCR amplicon high-throughput sequencing approach to sequence all 7 segments and to quantify polymorphisms in 91 RRV isolates collected from 16 states in the USA. Analysis revealed insertion/deletion (indel) polymorphisms primarily in the 5' and 3' non-coding, but also within coding regions, including some resulting in changes of protein length. Phylogenetic analysis showed little geographical structuring, suggesting that topography does not have a strong influence on virus evolution. Overall, the virus populations were homogeneous, possibly because of regular movement of plants, the recent emergence of RRV and/or because the virus is under strong purification selection to preserve its integrity and biological functions.

A Novel Divergent Geminivirus Identified in Asymptomatic New World Cactaceae Plants

Cactaceae comprise a diverse and iconic group of flowering plants which are almost exclusively indigenous to the New World. The wide variety of growth forms found amongst the cacti have led to the trafficking of many species throughout the world as ornamentals. Despite the evolution and physiological properties of these plants having been extensively studied, little research has focused on cactus-associated viral communities. While only single-stranded RNA viruses had ever been reported in cacti, here we report the discovery of cactus-infecting single-stranded DNA viruses. These viruses all apparently belong to a single divergent species of the family Geminiviridae and have been tentatively named Opuntia virus 1 (OpV1). A total of 79 apparently complete OpV1 genomes were recovered from 31 different cactus plants (belonging to 20 different cactus species from both the Cactoideae and Opuntioideae clades) and from nine cactus-feeding cochineal insects (Dactylopius sp.) sampled in the USA and Mexico. These 79 OpV1 genomes all share > 78.4% nucleotide identity with one another and < 64.9% identity with previously characterized geminiviruses. Collectively, the OpV1 genomes display evidence of frequent recombination, with some genomes displaying up to five recombinant regions. In one case, recombinant regions span ~40% of the genome. We demonstrate that an infectious clone of an OpV1 genome can replicate in Nicotiana benthamiana and Opuntia microdasys. In addition to expanding the inventory of viruses that are known to infect cacti, the OpV1 group is so distantly related to other known geminiviruses that it likely represents a new geminivirus genus. It remains to be determined whether, like its cactus hosts, its geographical distribution spans the globe.

Viromes of Ten Alfalfa Plants in Australia Reveal Diverse Known Viruses and a Novel RNA Virus

Alfalfa plants in the field can display a range of virus-like symptoms, especially when grown over many years for seed production. Most known alfalfa viruses have RNA genomes, some of which can be detected using diagnostic assays, but many viruses of alfalfa are not well characterized. This study aims to identify the RNA and DNA virus complexes associated with alfalfa plants in Australia. To maximize the detection of RNA viruses, we purified double-stranded RNA (dsRNA) for high throughput sequencing and characterized the viromes of ten alfalfa samples that showed diverse virus-like symptoms. Using Illumina sequencing of tagged cDNA libraries from immune-captured dsRNA, we identified sequences of the single-stranded RNA viruses, alfalfa mosaic virus (AMV), bean leafroll virus, a new emaravirus tentatively named alfalfa ringspot-associated virus, and persistent dsRNA viruses belonging to the families Amalgaviridae and Partitiviridae. Furthermore, rolling circle amplification and restriction enzyme digestion revealed the complete genome of chickpea chlorosis Australia virus, a mastrevirus (family Geminiviridae) previously reported only from chickpea and French bean that was 97% identical to the chickpea isolate. The sequence data also enabled the assembly of the first complete genome (RNAs 1–3) of an Australian AMV isolate from alfalfa.

Sixty Years After the First Description: Genome Sequence and Biological Characterization of European Wheat Striate Mosaic Virus Infecting Cereal Crops

High-throughput sequencing technologies were used to identify plant viruses in cereal samples surveyed from 2012 to 2017. Fifteen genome sequences of a tenuivirus infecting wheat, oats, and spelt in Estonia, Norway, and Sweden were identified and characterized by their distances to other tenuivirus sequences. Like most tenuiviruses, the genome of this tenuivirus contains four genomic segments. The isolates found from different countries shared at least 92% nucleotide sequence identity at the genome level. The planthopper Javesella pellucida was identified as a vector of the virus. Laboratory transmission tests using this vector indicated that wheat, oats, barley, rye, and triticale, but none of the tested pasture grass species (Alopecurus pratensis, Dactylis glomerata, Festuca rubra, Lolium multiflorum, Phleum pratense, and Poa pratensis), are susceptible. Taking into account the vector and host range data, the tenuivirus we have found most probably represents European wheat striate mosaic virus first identified about 60 years ago. Interestingly, whereas we were not able to infect any of the tested cereal species mechanically, Nicotiana benthamiana was infected via mechanical inoculation in laboratory conditions, displaying symptoms of yellow spots and vein clearing evolving into necrosis, eventually leading to plant death. Surprisingly, one of the virus genome segments (RNA2) encoding both a putative host systemic movement enhancer protein and a putative vector transmission factor was not detected in N. benthamiana after several passages even though systemic infection was observed, raising fundamental questions about the role of this segment in the systemic spread in several hosts.

Fundamental Aspects of Plant Viruses-An Overview on Focus Issue Articles

Given the importance of and rapid research progress in plant virology in recent years, this Focus Issue broadly emphasizes advances in fundamental aspects of virus infection cycles and epidemiology. This Focus Issue comprises three review articles and 18 research articles. The research articles cover broad research areas on the identification of novel viruses, the development of detection methods, reverse genetics systems and functional genomics for plant viruses, vector and seed transmission studies, viral population studies, virus-virus interactions and their effect on vector transmission, and management strategies of viral diseases. The three review articles discuss recent developments in application of prokaryotic clustered regularly interspaced short palindromic repeats/CRISPR-associated genes (CRISPR/Cas) technology for plant virus resistance, mixed viral infections and their role in disease synergism and cross-protection, and viral transmission by whiteflies. The following briefly summarizes the articles appearing in this Focus Issue.

Detection of Plant Viruses and Disease Management: Relevance of Genetic Diversity and Evolution

Plant viruses cause considerable economic losses and are a threat for sustainable agriculture. The frequent emergence of new viral diseases is mainly due to international trade, climate change, and the ability of viruses for rapid evolution. Disease control is based on two strategies: i) immunization (genetic resistance obtained by plant breeding, plant transformation, cross-protection, or others), and ii) prophylaxis to restrain virus dispersion (using quarantine, certification, removal of infected plants, control of natural vectors, or other procedures). Disease management relies strongly on a fast and accurate identification of the causal agent. For known viruses, diagnosis consists in assigning a virus infecting a plant sample to a group of viruses sharing common characteristics, which is usually referred to as species. However, the specificity of diagnosis can also reach higher taxonomic levels, as genus or family, or lower levels, as strain or variant. Diagnostic procedures must be optimized for accuracy by detecting the maximum number of members within the group (sensitivity as the true positive rate) and distinguishing them from outgroup viruses (specificity as the true negative rate). This requires information on the genetic relationships within-group and with members of other groups. The influence of the genetic diversity of virus populations in diagnosis and disease management is well documented, but information on how to integrate the genetic diversity in the detection methods is still scarce. Here we review the techniques used for plant virus diagnosis and disease control, including characteristics such as accuracy, detection level, multiplexing, quantification, portability, and designability. The effect of genetic diversity and evolution of plant viruses in the design and performance of some detection and disease control techniques are also discussed. High-throughput or next-generation sequencing provides broad-spectrum and accurate identification of viruses enabling multiplex detection, quantification, and the discovery of new viruses. Likely, this technique will be the future standard in diagnostics as its cost will be dropping and becoming more affordable.

High-Throughput Sequencing Reveals Bell Pepper Endornavirus Infection in Pepper (Capsicum annum) in Slovakia and Enables Its Further Molecular Characterization

Ribosomal RNA-depleted total RNAs from a sweet pepper plant (Capsicum annuum, labelled as N65) grown in western Slovakia and showing severe virus-like symptoms (chlorosis, mottling and deformation of leaf lamina) were subjected to high-throughput sequencing (HTS) on an Illumina MiSeq platform. The de novo assembly of ca. 5.5 million reads, followed by mapping to the reference sequences, revealed the coinfection of pepper by several viruses; i.e., cucumber mosaic virus (CMV), watermelon mosaic virus (WMV), pepper cryptic virus 2 (PCV2) and bell pepper endornavirus (BPEV). A complete polyprotein-coding genomic sequence (14.6 kb) of BPEV isolate N65 was determined. A comparison of BPEV-N65 sequences with BPEV genomes available in GenBank showed 86.1% to 98.6% identity at the nucleotide level. The close phylogenetic relationship with isolates from India and China resulted in their distinct grouping compared to the other BPEV isolates. Further analysis has revealed the presence of BPEV in sweet or chili peppers obtained from various sources and locations in Slovakia (plants grown in gardens, greenhouse or retail shop). Additionally, the partial sequencing of two genomic portions from 15 BPEV isolates revealed that the Slovak isolates segregated into two molecular clusters, indicating a genetically distinct population (mean inter-group nucleotide divergence reaching 12.7% and 14.5%, respectively, based on the genomic region targeted). Due to the mix infections of BPEV-positive peppers by potato virus Y (PVY) and/or CMV, the potential role of individual viruses in the observed symptomatology could not be determined. This is the first evidence and characterization of BPEV from the central European region.

Phytovirome Analysis of Wild Plant Populations: Comparison of Double-Stranded RNA and Virion-Associated Nucleic Acid Metagenomic Approaches

Metagenomic studies have indicated that the diversity of plant viruses was until recently far underestimated. As important components of ecosystems, there is a need to explore the diversity and richness of the viruses associated with plant populations and to understand the drivers shaping their diversity in space and time. Two viral sequence enrichment approaches, double-stranded RNA (dsRNA) and virion-associated nucleic acids (VANA), have been used and compared here for the description of the virome of complex plant pools representative of the most prevalent plant species in unmanaged and cultivated ecosystems. A novel bioinformatics strategy was used to assess viral richness not only at the family level but also by determining operational taxonomic units (OTU) following the clustering of conserved viral domains. A large viral diversity dominated by novel dsRNA viruses was detected in all sites, while a large between-site variability limited the ability to draw a clear conclusion on the impact of cultivation. A trend for a higher diversity of dsRNA viruses was nevertheless detected in unmanaged sites (118 versus 77 unique OTUs). The dsRNA-based approach consistently revealed a broader and more comprehensive diversity for RNA viruses than the VANA approach, whatever the assessment criterion. In addition, dissimilarity analyses indicated both approaches to be largely reproducible but not necessarily convergent. These findings illustrate features of phytoviromes in various ecosystems and a novel strategy for precise virus richness estimation. These results allow us to reason methodological choices in phytovirome studies and likely in other virome studies where RNA viruses are the focal taxa.IMPORTANCE There are today significant knowledge gaps on phytovirus populations and on the drivers impacting them but also on the comparative performance-methodological approaches for their study. We used and compared two viral sequence enrichment approaches, double-stranded RNAs (dsRNA) and virion-associated nucleic acids (VANA), for phytovirome description in complex pools representative of the most prevalent plant species in unmanaged and cultivated ecosystems. Viral richness was assessed by determining operational taxonomic units (OTU) following the clustering of conserved viral domains. There is some limited evidence of an impact of cultivation on viral populations. These results provide data allowing us to reason the methodological choices in virome studies. For researchers primarily interested in RNA viruses, the dsRNA approach is recommended because it consistently provided a more comprehensive description of the analyzed phytoviromes, but it understandably underrepresented DNA viruses and bacteriophages.

Metagenomic analysis of virome cross-talk between cultivated Solanum lycopersicum and wild Solanum nigrum

Wild plants and weeds growing close to crops constitute a potential reservoir for future epidemies or for the emergence of novel viruses but the frequency and directionality of viral flow between cultivated and wild plants remains poorly documented in many cases. Here, we studied the diversity of viral populations between tomato (Solanum lycopersicum) and neighboring european black nightshade (Solanum nigrum) using high throughput sequencing (HTS) based metagenomics. A large variability in virome richness with only 17.9% shared Operational Taxonomy Units between tomato and nightshade, but this richness could not be linked to a particular host or to local conditions. A detailed population analysis based on assembled contigs for potato virus Y (PVY), broad wilt bean virus 1 and a new ilarvirus tentatively named Solanum nigrum ilarvirus 1 provides information on the circulation of these viruses between these two Solanum species and enriches our knowledge of the tomato virome.

Full genome characterization of 12 citrus tatter leaf virus isolates for the development of a detection assay

Citrus tatter leaf virus (CTLV) threatens citrus production worldwide because it induces bud-union crease on the commercially important Citrange (Poncirus trifoliata × Citrus sinensis) rootstocks. However, little is known about its genomic diversity and how such diversity may influence virus detection. In this study, full-length genome sequences of 12 CTLV isolates from different geographical areas, intercepted and maintained for the past 60 years at the Citrus Clonal Protection Program (CCPP), University of California, Riverside, were characterized using next generation sequencing. Genome structure and sequence for all CTLV isolates were similar to Apple stem grooving virus (ASGV), the type species of Capillovirus genus of the Betaflexiviridae family. Phylogenetic analysis highlighted CTLV’s point of origin in Asia, the virus spillover to different plant species and the bottleneck event of its introduction in the United States of America (USA). A reverse transcription quantitative polymerase chain reaction assay was designed at the most conserved genome area between the coat protein and the 3’-untranslated region (UTR), as identified by the full genome analysis. The assay was validated with different parameters (e.g. specificity, sensitivity, transferability and robustness) using multiple CTLV isolates from various citrus growing regions and it was compared with other published assays. This study proposes that in the era of powerful affordable sequencing platforms the presented approach of systematic full-genome sequence analysis of multiple virus isolates, and not only a small genome area of a small number of isolates, becomes a guideline for the design and validation of molecular virus detection assays, especially for use in high value germplasm programs.

Next generation sequencing and proteomics in plant virology: how is Colombia doing?

Crop production and trade are two of the most economically important activities in Colombia, and viral diseases cause a high negative impact to agricultural sector. Therefore, the detection, diagnosis, control, and management of viral diseases are crucial. Currently, Next-Generation Sequencing (NGS) and ‘Omic’ technologies constitute a right-hand tool for the discovery of novel viruses and for studying virus-plant interactions. This knowledge allows the development of new viral diagnostic methods and the discovery of key components of infectious processes, which could be used to generate plants resistant to viral infections. Globally, crop sciences are advancing in this direction. In this review, advancements in ‘omic’ technologies and their different applications in plant virology in Colombia are discussed. In addition, bioinformatics pipelines and resources for omics data analyses are presented. Due to their decreasing prices, NGS technologies are becoming an affordable and promising means to explore many phytopathologies affecting a wide variety of Colombian crops so as to improve their trade potential.

Discovery of Four Novel Viruses Associated with Flower Yellowing Disease of Green Sichuan Pepper (Zanthoxylum Armatum) by Virome Analysis

An emerging virus-like flower yellowing disease (FYD) of green Sichuan pepper (Zanthoxylum armatum v. novemfolius) has been recently reported. Four new RNA viruses were discovered in the FYD-affected plant by the virome analysis using high-throughput sequencing of transcriptome and small RNAs. The complete genomes were determined, and based on the sequence and phylogenetic analysis, they are considered to be new members of the genera Nepovirus (Secoviridae), Idaeovirus (unassigned), Enamovirus (Luteoviridae), and Nucleorhabdovirus (Rhabdoviridae), respectively. Therefore, the tentative names corresponding to these viruses are green Sichuan pepper-nepovirus (GSPNeV), -idaeovirus (GSPIV), -enamovirus (GSPEV), and -nucleorhabdovirus (GSPNuV). The viral population analysis showed that GSPNeV and GSPIV were dominant in the virome. The small RNA profiles of these viruses are in accordance with the typical virus-plant interaction model for Arabidopsis thaliana. Rapid and sensitive RT-PCR assays were developed for viral detection, and used to access the geographical distributions. The results revealed a correlation between GSPNeV and the FYD. The viruses pose potential threats to the normal production of green Sichuan pepper in the affected areas due to their natural transmission and wide spread in fields. Collectively, our results provide useful information regarding taxonomy, transmission and pathogenicity of the viruses as well as management of the FYD.