Applications of next generation high throughput sequencing technologies in characterization, discovery and molecular interaction of plant viruses

Characterization of Hibiscus Chlorotic Ringspot Virus-Derived vsiRNAs from Infected Hibiscus rosa-sinensis in China

Lots of progress have been made about pathogen system of Hibiscus rosa-sinensis and hibiscus chlorotic ringspot virus (HCRSV), however, interactions between H. rosa-sinensis and HCRSV remain largely unknown. Hereon, firstly, HCRSV infection in H. rosa-sinensis from Zhangzhou city of China was confirmed by traditional electron microscopy, modern reverse transcription polymerase chain reaction and RNA-seq methods. Secondly, sequence feature analysis showed the full-length sequence of HCRSV-ZZ was 3,909 nucleotides (nt) in length and had a similar genomic structure with other carmovirus. It contains a 5' untranslated region (UTR), followed by seven open reading frames encoding for P28, P23, P81, P8, P9, P38, and P25, and the last a 3-terminal UTR. Thirdly, HCRSV- ZZ-derived vsiRNAs were identified and characterized for the first time from disease H. rosa-sinensis through sRNA-seq to reveal interactions between pathogen ant plant host. It was shown that the majority of HCRSV-ZZ-derived vsiRNAs were 21 nt, 22 nt, and 20 nt, with 21 nt being most abundant. The 5'-terminal nucleotide of HCRSV-ZZ vsiRNAs preferred U and C. HCRSV-ZZ vsiRNAs derived predominantly (72%) from the viral genome positive-strand RNA. The distribution of HCRSV-ZZ vsiRNAs along the viral genome is generally even, with some hot spots and cold spots forming in local regions. These hot spots and cold spots could be corresponded to the regions of stem loop secondary structures forming in HCRSV-ZZ genome by nucleotide paring. Taken together, our findings certify HCRSV infection in H. rosa-sinensis and provide an insight into interaction between HCRSV and H. rosa-sinensis and contribute to the prevention and treatment of this virus.

Advances in understanding plant-pathogen interactions: insights from tomato as a model system

The impact of plant diseases coupled with climate change on agriculture worldwide cannot be overemphasized from negative effects on crop yield as well as economy to food insecurity. The model plants are essential for understanding the intricacies of plant-pathogen interactions. One of such plants is the tomato (Solanum lycopersicum L.). Researchers hope to increase tomato productivity and boost its resilience to pathogen attacks by utilizing OMICS and biotechnological methods. With an emphasis on tomato viral infections, this review summarizes significant discoveries and developments from earlier research. The analysis elucidates ongoing efforts to advance plant pathology by exploring the implications for sustainability and tomato production.

Recent advances and challenges in plant viral diagnostics

Accurate and timely diagnosis of plant viral infections plays a key role in effective disease control and maintaining agricultural productivity. Recent advances in the diagnosis of plant viruses have significantly expanded our ability to detect and monitor viral pathogens in agricultural crops. This review discusses the latest advances in diagnostic technologies, including both traditional methods and the latest innovations. Conventional methods such as enzyme-linked immunosorbent assay and DNA amplification-based assays remain widely used due to their reliability and accuracy. However, diagnostics such as next-generation sequencing and CRISPR-based detection offer faster, more sensitive and specific virus detection. The review highlights the main advantages and limitations of detection systems used in plant viral diagnostics including conventional methods, biosensor technologies and advanced sequence-based techniques. In addition, it also discusses the effectiveness of commercially available diagnostic tools and challenges facing modern diagnostic techniques as well as future directions for improving informed disease management strategies. Understanding the main features of available diagnostic methodologies would enable stakeholders to choose optimal management strategies against viral threats and ensure global food security.

Deciphering the virome of Chunkung (Cnidium officinale) showing dwarfism-like symptoms via a high-throughput sequencing analysis

BACKGROUND

Viruses have notable effects on agroecosystems, wherein they can adversely affect plant health and cause problems (e.g., increased biosecurity risks and economic losses). However, our knowledge of their diversity and interactions with specific host plants in ecosystems remains limited. To enhance our understanding of the roles that viruses play in agroecosystems, comprehensive analyses of the viromes of a wide range of plants are essential. High-throughput sequencing (HTS) techniques are useful for conducting impartial and unbiased investigations of plant viromes, ultimately forming a basis for generating further biological and ecological insights. This study was conducted to thoroughly characterize the viral community dynamics in individual plants.

RESULTS

An HTS-based virome analysis in conjunction with proximity sampling and a tripartite network analysis were performed to investigate the viral diversity in chunkung (Cnidium officinale) plants. We identified 61 distinct chunkung plant-associated viruses (27 DNA and 34 RNA viruses) from 21 known genera and 6 unclassified genera in 14 known viral families. Notably, 12 persistent viruses (7 DNA and 5 RNA viruses) were exclusive to dwarfed chunkung plants. The detection of viruses from the families Partitiviridae, Picobirnaviridae, and Spinareoviridae only in the dwarfed plants suggested that they may contribute to the observed dwarfism. The co-infection of chunkung by multiple viruses is indicative of a dynamic and interactive viral ecosystem with significant sequence variability and evidence of recombination.

CONCLUSIONS

We revealed the viral community involved in chunkung. Our findings suggest that chunkung serves as a significant reservoir for a variety of plant viruses. Moreover, the co-infection rate of individual plants was unexpectedly high. Future research will need to elucidate the mechanisms enabling several dozen viruses to co-exist in chunkung. Nevertheless, the important insights into the chunkung virome generated in this study may be relevant to developing effective plant viral disease management and control strategies.

Viruses Infecting Trees and Herbs That Produce Edible Fleshy Fruits with a Prominent Value in the Global Market: An Evolutionary Perspective

Trees and herbs that produce fruits represent the most valuable agricultural food commodities in the world. However, the yield of these crops is not fully achieved due to biotic factors such as bacteria, fungi, and viruses. Viruses are capable of causing alterations in plant growth and development, thereby impacting the yield of their hosts significantly. In this work, we first compiled the world's most comprehensive list of known edible fruits that fits our definition. Then, plant viruses infecting those trees and herbs that produce fruits with commercial importance in the global market were identified. The identified plant viruses belong to 30 families, most of them containing single-stranded RNA genomes. Importantly, we show the overall picture of the host range for some virus families following an evolutionary approach. Further, the current knowledge about plant-virus interactions, focusing on the main disorders they cause, as well as yield losses, is summarized. Additionally, since accurate diagnosis methods are of pivotal importance for viral diseases control, the current and emerging technologies for the detection of these plant pathogens are described. Finally, the most promising strategies employed to control viral diseases in the field are presented, focusing on solutions that are long-lasting.

Whole-Genome Characterization of Alfalfa Mosaic Virus Obtained from Metagenomic Analysis of Vinca minor and Wisteria sinensis in Iran: with Implications for the Genetic Structure of the Virus

Alfalfa mosaic virus (AMV), an economically important pathogen, is present worldwide with a very wide host range. This work reports for the first time the infection of Vinca minor and Wisteria sinensis with AMV using RNA sequencing and reverse transcription polymerase chain reaction confirmation. De novo assembly and annotating of contigs revealed that RNA1, RNA2, and RNA3 genomic fragments consist of 3,690, 2,636, and 2,057 nucleotides (nt) for IR-VM and 3,690, 2,594, and 2,057 nt for IR-WS. RNA1 and RNA3 segments of IR-VM and IR-WS closely resembled those of the Chinese isolate HZ, with 99.23-99.26% and 98.04-98.09% nt identity, respectively. Their RNA2 resembled that of Canadian isolate CaM and American isolate OH-2-2017, with 97.96-98.07% nt identity. The P2 gene revealed more nucleotide diversity compared with other genes. Genes in the AMV genome were under dominant negative selection during evolution, and the P1 and coat protein (CP) proteins were subject to the strongest and weakest purifying selection, respectively. In the population genetic analysis based on the CP gene sequences, all 107 AMV isolates fell into two main clades (A, B) and isolates of clade A were further divided into three groups with significant subpopulation differentiation. The results indicated moderate genetic variation within and no clear geographic or genetic structure between the studied populations, implying moderate gene flow can play an important role in differentiation and distribution of genetic diversity among populations. Several factors have shaped the genetic structure and diversity of AMV: selection, recombination/reassortment, gene flow, and random processes such as founder effects.

Characterization of Hibiscus Latent Fort Pierce Virus-Derived siRNAs in Infected Hibiscus rosa-sinensis in China

Although limited progress have been made about pathogen system of Hibiscus rosa-sinensis and Hibiscus latent Fort Pierce virus (HLFPV), interaction between plant host and pathogen remain largely unknown, which led to deficiency of effective measures to control disease of hibiscus plants caused by HLFPV. In this study, infection of HLFPV in Hibiscus rosa-sinensis was firstly confirmed for the first time by traditional electron microscopy, modern reverse transcription polymerase chain reaction and RNA-seq methods in China (HLFPV-Ch). Sequence properties analyzing suggested that the full-length sequences (6,465 nt) of HLFPV-Ch had a high sequence identity and a similar genomic structure with other tobamoviruses. It includes a 5'-terminal untranslated region (UTR), followed by four open reading frames encoding for a 128.5-kDa replicase, a 186.5-kDa polymerase, a 31-kDa movement protein, 17.6-kDa coat protein, and the last a 3'-terminal UTR. Furthermore, HLFPV-Ch-derived virus-derived siRNAs (vsiRNAs) ant its putative target genes, reported also for the first time, were identified and characterized from disease Hibiscus rosa-sinensis through sRNA-seq and Patmatch server to investigate the interaction in this pathogen systems. HLFPV-Ch-derived vsiRNAs demonstrated several general and specific characteristics. Gene Ontology classification revealed predicted target genes by vsiRNAs are involved in abroad range of cellular component, molecular function and biological processes. Taken together, for first time, our results certified the HLFPV infection in China and provide an insight into interaction between HLFPV and Hibiscus rosa-sinensis.

Illuminating an Ecological Blackbox: Using High Throughput Sequencing to Characterize the Plant Virome Across Scales

The ecology of plant viruses began to be explored at the end of the 19th century. Since then, major advances have revealed mechanisms of virus-host-vector interactions in various environments. These advances have been accelerated by new technlogies for virus detection and characterization, most recently including high throughput sequencing (HTS). HTS allows investigators, for the first time, to characterize all or nearly all viruses in a sample without a priori information about which viruses might be present. This powerful approach has spurred new investigation of the viral metagenome (virome). The rich virome datasets accumulated illuminate important ecological phenomena such as virus spread among host reservoirs (wild and domestic), effects of ecosystem simplification caused by human activities (and agriculture) on the biodiversity and the emergence of new viruses in crops. To be effective, however, HTS-based virome studies must successfully navigate challenges and pitfalls at each procedural step, from plant sampling to library preparation and bioinformatic analyses. This review summarizes major advances in plant virus ecology associated with technological developments, and then presents important considerations and best practices for HTS use in virome studies.

Molecular Characterization of Hovenia Dulcis-Associated Virus 1 (HDaV1) and 2 (HDaV2): New Tentative Species within the Order Picornavirales

In a systematic field survey for plant-infecting viruses, leaf tissues were collected from trees showing virus-like symptoms in Brazil. After viral enrichment, total RNA was extracted and sequenced using the MiSeq platform (Illumina). Two nearly full-length picorna-like genomes of 9534 and 8158 nucleotides were found associated with Hovenia dulcis (Rhamnaceae family). Based upon their genomic information, specific primers were synthetized and used in RT-PCR assays to identify plants hosting the viral sequences. The larger contig was tentatively named as Hovenia dulcis-associated virus 1 (HDaV1), and it exhibited low nucleotide and amino acid identities with Picornavirales species. The smaller contig was related to insect-associated members of the Dicistroviridae family but exhibited a distinct genome organization with three non-overlapping open reading frames (ORFs), and it was tentatively named as Hovenia dulcis-associated virus 2 (HDaV2). Phylogenetic analysis using the amino acid sequence of RNA-dependent RNA polymerase (RdRp) revealed that HDaV1 and HDaV2 clustered in distinct groups, and both viruses were tentatively assigned as new members of the order Picornavirales. HDaV2 was assigned as a novel species in the Dicistroviridae family. The 5′ ends of both viruses are incomplete. In addition, a nucleotide composition analysis (NCA) revealed that HDaV1 and HDaV2 have similarities with invertebrate-infecting viruses, suggesting that the primary host(s) of these novel virus species remains to be discovered.

siRNAs Derived from Cymbidium Mosaic Virus and Odontoglossum Ringspot Virus Down-modulated the Expression Levels of Endogenous Genes in Phalaenopsis equestris

Interplay between Cymbidium mosaic virus (CymMV)/Odontoglossum ringspot virus (ORSV) and its host plant Phalaenopsis equestris remain largely unknown, which led to deficiency of effective measures to control disease of P. equestris caused by infecting viruses. In this study, for the first time, we characterized viral small interfering RNAs (vsiRNAs) profiles in P. equestris co-infected with CymMV and ORSV through small RNA sequencing technology. CymMV and ORSV small interfering RNAs (siRNAs) demonstrated several general and specific/new characteristics. vsiRNAs, with A/U bias at the first nucleotide, were predominantly 21-nt long and they were derived predominantly (90%) from viral positive-strand RNA. 21-nt siRNA duplexes with 0-nt overhangs were the most abundant 21-nt duplexes, followed by 2-nt overhangs and then 1-nt overhangs 21-nt duplexes in infected P. equestris. Continuous but heterogeneous distribution and secondary structures prediction implied that vsiRNAs originate predominantly by direct Dicer-like enzymes cleavage of imperfect duplexes in the most folded regions of the positive strand of both viruses RNA molecular. Furthermore, we totally predicted 54 target genes by vsiRNAs with psRNATarget server, including disease/stress response-related genes, RNA interference core components, cytoskeleton-related genes, photosynthesis or energy supply related genes. Gene Ontology classification showed that a majority of the predicted targets were related to cellular components and cellular processes and performed a certain function. All target genes were down-regulated with different degree by vsiRNAs as shown by real-time reverse transcription polymerase chain reaction. Taken together, CymMV and ORSV siRNAs played important roles in interplay with P. equestris by down modulating the expression levels of endogenous genes in host plant.

Tissue culture and next-generation sequencing: A combined approach for detecting yam (Dioscorea spp.) viruses

In vitro culture offers many advantages for yam germplasm conservation, propagation and international distribution. However, low virus titres in the generated tissues pose a challenge for reliable virus detection, which makes it difficult to ensure that planting material is virus-free. In this study, we evaluated next-generation sequencing (NGS) for virus detection following yam propagation using a robust tissue culture methodology. We detected and assembled the genomes of novel isolates of already characterised viral species of the genera Badnavirus and Potyvirus, confirming the utility of NGS in diagnosing yam viruses and contributing towards the safe distribution of germplasm.

Cassava brown streak disease: historical timeline, current knowledge and future prospects

Cassava is the second most important staple food crop in terms of per capita calories consumed in Africa and holds potential for climate change adaptation. Unfortunately, productivity in East and Central Africa is severely constrained by two viral diseases: cassava mosaic disease (CMD) and cassava brown streak disease (CBSD). CBSD was first reported in 1936 from northeast Tanzania. For approximately 70 years, CBSD was restricted to coastal East Africa and so had a relatively low impact on food security compared with CMD. However, at the turn of the 21st century, CBSD re-emerged further inland, in areas around Lake Victoria, and it has since spread through many East and Central African countries, causing high yield losses and jeopardizing the food security of subsistence farmers. This recent re-emergence has attracted intense scientific interest, with studies shedding light on CBSD viral epidemiology, sequence diversity, host interactions and potential sources of resistance within the cassava genome. This review reflects on 80 years of CBSD research history (1936-2016) with a timeline of key events. We provide insights into current CBSD knowledge, management efforts and future prospects for improved understanding needed to underpin effective control and mitigation of impacts on food security.

Generic Amplicon Deep Sequencing to Determine Ilarvirus Species Diversity in Australian Prunus

The distribution of Ilarvirus species populations amongst 61 Australian Prunus trees was determined by next generation sequencing (NGS) of amplicons generated using a genus-based generic RT-PCR targeting a conserved region of the Ilarvirus RNA2 component that encodes the RNA dependent RNA polymerase (RdRp) gene. Presence of Ilarvirus sequences in each positive sample was further validated by Sanger sequencing of cloned amplicons of regions of each of RNA1, RNA2 and/or RNA3 that were generated by species specific PCRs and by metagenomic NGS. Prunus necrotic ringspot virus (PNRSV) was the most frequently detected Ilarvirus, occurring in 48 of the 61 Ilarvirus-positive trees and Prune dwarf virus (PDV) and Apple mosaic virus (ApMV) were detected in three trees and one tree, respectively. American plum line pattern virus (APLPV) was detected in three trees and represents the first report of APLPV detection in Australia. Two novel and distinct groups of Ilarvirus-like RNA2 amplicon sequences were also identified in several trees by the generic amplicon NGS approach. The high read depth from the amplicon NGS of the generic PCR products allowed the detection of distinct RNA2 RdRp sequence variant populations of PNRSV, PDV, ApMV, APLPV and the two novel Ilarvirus-like sequences. Mixed infections of ilarviruses were also detected in seven Prunus trees. Sanger sequencing of specific RNA1, RNA2, and/or RNA3 genome segments of each virus and total nucleic acid metagenomics NGS confirmed the presence of PNRSV, PDV, ApMV and APLPV detected by RNA2 generic amplicon NGS. However, the two novel groups of Ilarvirus-like RNA2 amplicon sequences detected by the generic amplicon NGS could not be associated to the presence of sequence from RNA1 or RNA3 genome segments or full Ilarvirus genomes, and their origin is unclear. This work highlights the sensitivity of genus-specific amplicon NGS in detection of virus sequences and their distinct populations in multiple samples, and the need for a standardized approach to accurately determine what constitutes an active, viable virus infection after detection by molecular based methods.

Analysis of intra-host genetic diversity of Prunus necrotic ringspot virus (PNRSV) using amplicon next generation sequencing

PCR amplicon next generation sequencing (NGS) analysis offers a broadly applicable and targeted approach to detect populations of both high- or low-frequency virus variants in one or more plant samples. In this study, amplicon NGS was used to explore the diversity of the tripartite genome virus, Prunus necrotic ringspot virus (PNRSV) from 53 PNRSV-infected trees using amplicons from conserved gene regions of each of PNRSV RNA1, RNA2 and RNA3. Sequencing of the amplicons from 53 PNRSV-infected trees revealed differing levels of polymorphism across the three different components of the PNRSV genome with a total number of 5040, 2083 and 5486 sequence variants observed for RNA1, RNA2 and RNA3 respectively. The RNA2 had the lowest diversity of sequences compared to RNA1 and RNA3, reflecting the lack of flexibility tolerated by the replicase gene that is encoded by this RNA component. Distinct PNRSV phylo-groups, consisting of closely related clusters of sequence variants, were observed in each of PNRSV RNA1, RNA2 and RNA3. Most plant samples had a single phylo-group for each RNA component. Haplotype network analysis showed that smaller clusters of PNRSV sequence variants were genetically connected to the largest sequence variant cluster within a phylo-group of each RNA component. Some plant samples had sequence variants occurring in multiple PNRSV phylo-groups in at least one of each RNA and these phylo-groups formed distinct clades that represent PNRSV genetic strains. Variants within the same phylo-group of each Prunus plant sample had ≥97% similarity and phylo-groups within a Prunus plant sample and between samples had less ≤97% similarity. Based on the analysis of diversity, a definition of a PNRSV genetic strain was proposed. The proposed definition was applied to determine the number of PNRSV genetic strains in each of the plant samples and the complexity in defining genetic strains in multipartite genome viruses was explored.

Deep Sequencing Analysis of RNAs from Citrus Plants Grown in a Citrus Sudden Death-Affected Area Reveals Diverse Known and Putative Novel Viruses

Citrus sudden death (CSD) has caused the death of approximately four million orange trees in a very important citrus region in Brazil. Although its etiology is still not completely clear, symptoms and distribution of affected plants indicate a viral disease. In a search for viruses associated with CSD, we have performed a comparative high-throughput sequencing analysis of the transcriptome and small RNAs from CSD-symptomatic and -asymptomatic plants using the Illumina platform. The data revealed mixed infections that included Citrus tristeza virus (CTV) as the most predominant virus, followed by the Citrus sudden death-associated virus (CSDaV), Citrus endogenous pararetrovirus (CitPRV) and two putative novel viruses tentatively named Citrus jingmen-like virus (CJLV), and Citrus virga-like virus (CVLV). The deep sequencing analyses were sensitive enough to differentiate two genotypes of both viruses previously associated with CSD-affected plants: CTV and CSDaV. Our data also showed a putative association of the CSD-symptomatic plants with a specific CSDaV genotype and a likely association with CitPRV as well, whereas the two putative novel viruses showed to be more associated with CSD-asymptomatic plants. This is the first high-throughput sequencing-based study of the viral sequences present in CSD-affected citrus plants, and generated valuable information for further CSD studies.

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.

Characterisation of siRNAs derived from new isolates of bamboo mosaic virus and their associated satellites in infected ma bamboo (Dendrocalamus latiflorus)

We characterised the virus-derived small interfering RNAs (vsiRNA) of bamboo mosaic virus (Ba-vsiRNAs) and its associated satellite RNA (satRNA)-derived siRNAs (satsiRNAs) in a bamboo plant (Dendrocalamus latiflorus) by deep sequencing. Ba-vsiRNAs and satsiRNAs of 21-22 nt in length, with both (+) and (-) polarity, predominated. The 5'-terminal base of Ba-vsiRNA was biased towards A, whereas a bias towards C/U was observed in sense satsiRNAs, and towards A in antisense satsiRNAs. A large set of bamboo genes were identified as potential targets of Ba-vsiRNAs and satsiRNAs, revealing RNA silencing-based virus-host interactions in plants. Moreover, we isolated and characterised new isolates of bamboo mosaic virus (BaMV; 6,350 nt) and BaMV-associated satRNA (satBaMV; 834 nt), designated BaMV-MAZSL1 and satBaMV-MAZSL1, respectively.

Comparative analysis of virus-derived small RNAs within cassava (Manihot esculenta Crantz) infected with cassava brown streak viruses

•

The 21-nt virus-derived small RNAs were predominant, followed by the 22-nt class.

•

Susceptible cassava genotypes accumulated higher CBSV- than UCBSV-derived small RNA.

•

Tolerant cassava genotype accumulated high CBSV- and low UCBSV-derived small RNAs.

•

AGO2, DCL2 and DCL4 were differentially regulated in CBSV/UCBSV-infected plants.

•

CBSV and UCBSV interact differently in the same host genetic background

Infection of plant cells by viral pathogens triggers RNA silencing, an innate antiviral defense mechanism. In response to infection, small RNAs (sRNAs) are produced that associate with Argonaute (AGO)-containing silencing complexes which act to inactivate viral genomes by posttranscriptional gene silencing (PTGS). Deep sequencing was used to compare virus-derived small RNAs (vsRNAs) in cassava genotypes NASE 3, TME 204 and 60444 infected with the positive sense single-stranded RNA (+ssRNA) viruses cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV), the causal agents of cassava brown streak disease (CBSD). An abundance of 21–24 nt vsRNAs was detected and mapped, covering the entire CBSV and UCBSV genomes. The 21 nt vsRNAs were most predominant, followed by the 22 nt class with a slight bias toward sense compared to antisense polarity, and a bias for adenine and uracil bases present at the 5′-terminus. Distribution and frequency of vsRNAs differed between cassava genotypes and viral genomes. In susceptible genotypes TME 204 and 60444, CBSV-derived sRNAs were seen in greater abundance than UCBSV-derived sRNAs. NASE 3, known to be resistant to UCBSV, accumulated negligible UCBSV-derived sRNAs but high populations of CBSV-derived sRNAs. Transcript levels of cassava homologues of AGO2, DCL2 and DCL4, which are central to the gene-silencing complex, were found to be differentially regulated in CBSV- and UCBSV-infected plants across genotypes, suggesting these proteins play a role in antiviral defense. Irrespective of genotype or viral pathogen, maximum populations of vsRNAs mapped to the cytoplasmic inclusion, P1 and P3 protein-encoding regions. Our results indicate disparity between CBSV and UCBSV host-virus interaction mechanisms, and provide insight into the role of virus-induced gene silencing as a mechanism of resistance to CBSD.

Seed Transmission of Soybean vein necrosis virus: The First Tospovirus Implicated in Seed Transmission

Soybean vein necrosis virus (SVNV; genus Tospovirus; Family Bunyaviridae) is a negative-sense single-stranded RNA virus that has been detected across the United States and in Ontario, Canada. In 2013, a seed lot of a commercial soybean variety (Glycine max) with a high percentage of discolored, deformed and undersized seed was obtained. A random sample of this seed was planted in a growth room under standard conditions. Germination was greater than 90% and the resulting seedlings looked normal. Four composite samples of six plants each were tested by reverse transcription polymerase chain reaction (RT-PCR) using published primers complimentary to the S genomic segment of SVNV. Two composite leaflet samples retrieved from seedlings yielded amplicons with a size and sequence predictive of SVNV. Additional testing of twelve arbitrarily selected individual plants resulted in the identification of two SVNV positive plants. Experiments were repeated by growing seedlings from the same seed lot in an isolated room inside a thrips-proof cage to further eliminate any external source of infection. Also, increased care was taken to reduce any possible PCR contamination. Three positive plants out of forty-eight were found using these measures. Published and newly designed primers for the L and M RNAs of SVNV were also used to test the extracted RNA and strengthen the diagnosis of viral infection. In experiments, by three scientists, in two different labs all three genomic RNAs of SVNV were amplified in these plant materials. RNA-seq analysis was also conducted using RNA extracted from a composite seedling sample found to be SVNV-positive and a symptomatic sample collected from the field. This analysis revealed both sense and anti-sense reads from all three gene segments in both samples. We have shown that SVNV can be transmitted in seed to seedlings from an infected seed lot at a rate of 6%. To our knowledge this is the first report of seed-transmission of a Tospovirus.

Characterization of New Isolates of Apricot vein clearing-associated virus and of a New Prunus-Infecting Virus: Evidence for Recombination as a Driving Force in Betaflexiviridae Evolution

Double stranded RNAs from Prunus samples gathered from various surveys were analyzed by a deep-sequencing approach. Contig annotations revealed the presence of a potential new viral species in an Azerbaijani almond tree (Prunus amygdalus) and its genome sequence was completed. Its genomic organization is similar to that of the recently described Apricot vein clearing associated virus (AVCaV) for which two new isolates were also characterized, in a similar fashion, from two Japanese plums (Prunus salicina) from a French germplasm collection. The amino acid identity values between the four proteins encoded by the genome of the new virus have identity levels with those of AVCaV which fall clearly outside the species demarcation criteria. The new virus should therefore be considered as a new species for which the name of Caucasus prunus virus (CPrV) has been proposed. Phylogenetic relationships and nucleotide comparisons suggested that together with AVCaV, CPrV could define a new genus (proposed name: Prunevirus) in the family Betaflexiviridae. A molecular test targeting both members of the new genus was developed, allowing the detection of additional AVCaV isolates, and therefore extending the known geographical distribution and the host range of AVCaV. Moreover, the phylogenetic trees reconstructed with the amino acid sequences of replicase, movement and coat proteins of representative Betaflexiviridae members suggest that Citrus leaf blotch virus (CLBV, type member of the genus Citrivirus) may have evolved from a recombination event involving a Prunevirus, further highlighting the importance of recombination as a driving force in Betaflexiviridae evolution. The sequences reported in the present manuscript have been deposited in the GenBank database under accession numbers KM507061-KM504070.

Comparison of Next-Generation Sequencing Versus Biological Indexing for the Optimal Detection of Viral Pathogens in Grapevine

A bioassay is routinely used to determine the viral phytosanitary status of commercial grapevine propagation material in many countries around the world. That test is based on the symptoms developed in the field by specific indicator host plants that are graft-inoculated from the vines being tested. We compared the bioassay against next-generation sequencing (NGS) analysis of grapevine material. NGS is a laboratory procedure that catalogs the genomic sequences of the viruses and other pathogens extracted as DNA and RNA from infected vines. NGS analysis was found to be superior to the standard bioassay in detection of viruses of agronomic significance, including virus infections at low titers. NGS was also found to be superior to the bioassay in its comprehensiveness, the speed of its analysis, and for the discovery of novel, uncharacterized viruses.

Deep sequencing reveals the complete genome and evidence for transcriptional activity of the first virus-like sequences identified in Aristotelia chilensis (Maqui Berry)

Here, we report the genome sequence and evidence for transcriptional activity of a virus-like element in the native Chilean berry tree Aristotelia chilensis. We propose to name the endogenous sequence as Aristotelia chilensis Virus 1 (AcV1). High-throughput sequencing of the genome of this tree uncovered an endogenous viral element, with a size of 7122 bp, corresponding to the complete genome of AcV1. Its sequence contains three open reading frames (ORFs): ORFs 1 and 2 shares 66%–73% amino acid similarity with members of the Caulimoviridae virus family, especially the Petunia vein clearing virus (PVCV), Petuvirus genus. ORF1 encodes a movement protein (MP); ORF2 a Reverse Transcriptase (RT) and a Ribonuclease H (RNase H) domain; and ORF3 showed no amino acid sequence similarity with any other known virus proteins. Analogous to other known endogenous pararetrovirus sequences (EPRVs), AcV1 is integrated in the genome of Maqui Berry and showed low viral transcriptional activity, which was detected by deep sequencing technology (DNA and RNA-seq). Phylogenetic analysis of AcV1 and other pararetroviruses revealed a closer resemblance with Petuvirus. Overall, our data suggests that AcV1 could be a new member of Caulimoviridae family, genus Petuvirus, and the first evidence of this kind of virus in a fruit plant.

Comparative analysis of virus-specific small RNA profiles of three biologically distinct strains of Potato virus Y in infected potato (Solanum tuberosum) cv. Russet Burbank

Deep sequencing technology has enabled the analysis of small RNA profiles of virus-infected plants and could provide insights into virus-host interactions. Potato virus Y is an economically important viral pathogen of potato worldwide. In this study, we investigated the nature and relative levels of virus-derived small interfering RNAs (vsiRNAs) in potato cv. Russet Burbank infected with three biologically distinct and economically important strains of PVY, the ordinary strain (PVY-O), tobacco veinal-necrotic strain (PVY-N) and tuber necrotic strain (PVY-NTN). The analysis showed an overall abundance of vsiRNAs of 20-24nt in PVY-infected plants. Considerable differences were present in the distribution of vsiRNAs as well as total small RNAs. The 21nt class was the most prevalent in PVY-infected plants irrespective of the virus strain, whereas in healthy potato plants, the 24nt class was the most dominant. vsiRNAs were derived from every position in the PVY genome, though certain hotspots were identified for each of the PVY strains. Among the three strains used, the population of vsiRNAs of different size classes was relatively different with PVY-NTN accumulating the highest level of vsiRNAs, while PVY-N infected plants had the least population of vsiRNAs. Unique vsiRNAs mapping to PVY genome in PVY-infected plants amounted to 3.13, 1.93 and 1.70% for NTN, N and O, respectively. There was a bias in the generation of vsiRNAs from the plus strand of the genome in comparison to the negative strand. The highest number of total vsiRNAs was from the cytoplasmic inclusion protein gene (CI) in PVY-O and PVY-NTN strains, whereas from PVY-N, the NIb gene produced maximum total vsiRNAs. These findings indicate that the three PVY strains interact differently in the same host genetic background and provided insights into virus-host interactions in an important food crop.

Current impact and future directions of high throughput sequencing in plant virus diagnostics

The ability to provide a fast, inexpensive and reliable diagnostic for any given viral infection is a key parameter in efforts to fight and control these ubiquitous pathogens. The recent developments of high-throughput sequencing (also called Next Generation Sequencing - NGS) technologies and bioinformatics have drastically changed the research on viral pathogens. It is now raising a growing interest for virus diagnostics. This review provides a snapshot vision on the current use and impact of high throughput sequencing approaches in plant virus characterization. More specifically, this review highlights the potential of these new technologies and their interplay with current protocols in the future of molecular diagnostic of plant viruses. The current limitations that will need to be addressed for a wider adoption of high-throughput sequencing in plant virus diagnostics are thoroughly discussed.

Differential expression of tomato spotted wilt virus-derived viral small RNAs in infected commercial and experimental host plants

BACKGROUND

Viral small RNAs (vsiRNAs) in the infected host can be generated from viral double-stranded RNA replicative intermediates, self-complementary regions of the viral genome or from the action of host RNA-dependent RNA polymerases on viral templates. The vsiRNA abundance and profile as well as the endogenous small RNA population can vary between different hosts infected by the same virus influencing viral pathogenicity and host response. There are no reports on the analysis of vsiRNAs of Tomato spotted wilt virus (TSWV), a segmented negative stranded RNA virus in the family Bunyaviridae, with two of its gene segments showing ambisense gene arrangement. The virus causes significant economic losses to numerous field and horticultural crops worldwide.

PRINCIPAL FINDINGS

Tomato spotted wilt virus (TSWV)-specific vsiRNAs were characterized by deep sequencing in virus-infected experimental host Nicotiana benthamiana and a commercial, susceptible host tomato. The total small (s) RNA reads in TSWV-infected tomato sample showed relatively equal distribution of 21, 22 and 24 nt, whereas N. benthamiana sample was dominated by 24 nt total sRNAs. The number of vsiRNA reads detected in tomato was many a magnitude (~350:1) higher than those found in N. benthamiana, however the profile of vsiRNAs in terms of relative abundance 21, 22 and 24 nt class size was similar in both the hosts. Maximum vsiRNA reads were obtained for the M RNA segment of TSWV while the largest L RNA segment had the least number of vsiRNAs in both tomato and N. benthamiana. Only the silencing suppressor, NSs, of TSWV recorded higher antisense vsiRNA with respect to the coding frame among all the genes of TSWV.

SIGNIFICANCE

Details of the origin, distribution and abundance of TSWV vsiRNAs could be useful in designing efficient targets for exploiting RNA interference for virus resistance. It also has major implications toward our understanding of the differential processing of vsiRNAs in antiviral defense and viral pathogenicity.