Deep sequencing analysis of viral short RNAs from an infected Pinot Noir grapevine

Viromes of Plants Determined by High-Throughput Sequencing of Virus-Derived siRNAs

Plants growing in open airfields can be infected by several viruses even as a multiple infection. Virus infection in crops can lead to a serious damage to the harvest. In addition, virus presence in grapevine, fruit trees, and tuberous vegetables, propagated vegetatively affects the phytosanitary status of the propagation material (both the rootstock and the variety) having profound effect on the lifetime and health of the new plantations. The fast evolution of sequencing techniques provides a new opportunity for metagenomics-based viral diagnostics. Small interfering (si) RNAs produced by the RNA silencing-based host immune system during viral infection can be sequenced by high-throughput techniques and analyzed for the presence of viruses, revealing the presence of all known viral pathogens in the sample and therefore opening new avenues in virus diagnostics. This method is based on Illumina sequencing and bioinformatics analysis of virus-derived siRNAs in the host. Here we describe a protocol for this challenging technique step by step with notes, to ensure success for every user.

Strategies for Engineering Virus Resistance in Potato

Potato (Solanum tuberosum L.) is an important vegetable crop that plays a pivotal role in the world, especially given its potential to feed the world population and to act as the major staple food in many developing countries. Every year, significant crop loss is caused by viral diseases due to a lack of effective agrochemical treatments, since only transmission by insect vectors can be combated with the use of insecticides, and this has been an important factor hindering potato production. With the rapid development of molecular biology and plant genetic engineering technology, transgenic approaches and non-transgenic techniques (RNA interference and CRISPR-cas9) have been effectively employed to improve potato protection against devastating viruses. Moreover, the availability of viral sequences, potato genome sequences, and host immune mechanisms has remarkably facilitated potato genetic engineering. In this study, we summarize the progress of antiviral strategies applied in potato through engineering either virus-derived or plant-derived genes. These recent molecular insights into engineering approaches provide the necessary framework to develop viral resistance in potato in order to provide durable and broad-spectrum protection against important viral diseases of solanaceous crops.

Interference of Small RNAs in Fusarium graminearum through FgGMTV1 Infection

Small RNA (sRNA) plays a central role in RNA silencing in fungi. The genome of Fusarium graminearum gemytripvirus 1 (FgGMTV1) is comprised of three DNA segments: DNA-A, DNA-B, and DNA-C. DNA-A and DNA-B are associated with fungal growth and virulence reduction. To elucidate the role of RNA silencing during the interactions of fungi and viruses, the sRNA profiles of F. graminearum in association with FgGMTV1 were established, using an FgGMTV1-free library (S-S), a library for infection with the DNA-A and DNA-B segments (S-AB), and a library for infection with the DNA-A, DNA-B, and DNA-C segments (S-ABC). A large amount of virus-derived sRNA (vsiRNA) was detected in the S-AB and S-ABC libraries, accounting for 9.9% and 13.8% of the total sRNA, respectively, indicating that FgGMTV1 triggers host RNA silencing. The total numbers of sRNA reads differed among the three libraries, suggesting that FgGMTV1 infection interferes with host RNA silencing. In addition, the relative proportions of the different sRNA lengths were altered in the S-AB and S-ABC libraries. The genome distribution patterns of the mapping of vsiRNA to DNA-A and DNA-B in the S-AB and S-ABC libraries were also different. These results suggest the influence of DNA-C on host RNA silencing. Transcripts targeted by vsiRNAs were enriched in pathways that included flavin adenine dinucleotide binding, protein folding, and filamentous growth.

Use of deep sequencing to profile small RNAs derived from tomato spotted wilt orthotospovirus and hippeastrum chlorotic ringspot orthotospovirus in infected Capsicum annuum

Virus-derived small RNAs are one of the key factors of RNA silencing in plant defence against viruses. We obtained virus-derived small interfering RNA profiles from Tomato spotted wilt orthotospovirus and Hippeastrum chlorotic ringspot orthotospovirus infected Capsicum annuum XX19 and XY11 by deep sequencing one day after inoculation. The vsiRNAs data were mapped to the TSWV and HCRV genomes, and the results showed that the vsiRNAs measured 19-24 nucleotides in length. Most of the vsiRNAs were mapped to the S segment of the viral genome. For XX19 and XY11 infected with HCRV, the distribution range of vsiRNAs in S RNA was 52.06-55.20%, while for XX19 and XY11 infected with TSWV, the distribution range of vsiRNAs in S RNA was 87.76-89.07%. The first base at the 5' end of the siRNA from TSWV and HCRV was primarily biased towards A, U, or C. Compared with mock-inoculated XX19 and XY11, the expression level of CaRDR1 was upregulated in TSWV- and HCRV-inoculated XX19 and XY11. CaAGO2 and CaAGO5 were upregulated in XY11 against HCRV infection, and CaRDR2 was downregulated in TSWV-infected XY11 and XX19. The profile of HCRV and TSWV vsiRNA verified in this study could be useful for selecting key vsiRNA such as those in disease-resistant varieties by artificially synthesizing amiRNA.

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.

Small RNA profiling analysis of two recombinant strains of potato virus Y in infected tobacco plants

Plant viral infections lead to accumulation of virus-derived small interfering RNAs (vsiRNAs) as a result of host defense mechanisms. High-throughput sequencing technology enables vsiRNA profiling analyses from virus infected plants, which provide important insights into virus-host interactions. Potato virus Y (PVY) is a detrimental plant pathogen that can infect a variety of solanaceous crops, e.g., potato, tobacco, tomato, and pepper. We analyzed and characterized vsiRNAs derived from Nicotiana tabacum cv. Samsun infected with two recombinant PVY strains, N-Wi and NTN. We observed that the average percentage of vsiRNAs derived from plants infected with N-Wi was higher than from plants infected with NTN, indicating that N-Wi invokes a stronger host response than NTN in tobacco. The size distribution pattern and polarity of vsiRNAs were similar between both virus strains with the 21 and 22 nucleotide (nt) vsiRNA classes as most predominant and the sense/antisense vsiRNAs ratio nearly equal in the 20-24 nt class. However, the percentage of sense vsiRNAs was significantly higher in the 25-26 nt long vsiRNAs. Distinct vsiRNA hotspots, identifying highly abundant reads of different unique vsiRNA sequences, were observed in both viral genomes. Previous studies found an A or U bias at the 5' terminal nucleotide position of 21 nt vsiRNAs; in contrast, our analysis revealed a C and U nucleotide bias. This study provides insights that will help further elucidate differential processing of vsiRNAs in plant antiviral defense.

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.

sRIS: A Small RNA Illustration System for Plant Next-Generation Sequencing Data Analysis

Small RNA (sRNA), such as microRNA (miRNA) and short interfering RNA, are well-known to control gene expression based on degradation of target mRNA in plants. A considerable amount of research has applied next-generation sequencing (NGS) to reveal the regulatory pathways of plant sRNAs. Consequently, numerous bioinformatics tools have been developed for the purpose of analyzing sRNA NGS data. However, most methods focus on the study of sRNA expression profiles or novel miRNAs predictions. The analysis of sRNA target genes is usually not integrated into their pipelines. As a result, there is still no means available for identifying the interaction mechanisms between host and virus or the synergistic effects between two viruses. For the present study, a comprehensive system, called the Small RNA Illustration System (sRIS), has been developed. This system contains two main components. The first is for sRNA overview analysis and can be used not only to identify miRNA but also to investigate virus-derived small interfering RNA. The second component is for sRNA target prediction, and it employs both bioinformatics calculations and degradome sequencing data to enhance the accuracy of target prediction. In addition, this system has been designed so that figures and tables for the outputs of each analysis can be easily retrieved and accessed, making it easier for users to quickly identify and quantify their results. sRIS is available at http://sris.itps.ncku.edu.tw/.

Metagenomic Analysis of Riesling Grapevine Reveals a Complex Virome Including Two New and Divergent Variants of Grapevine leafroll-associated virus 3

The virome of a major white wine grape of cultivar Riesling showing decline and leafroll disease symptoms was analyzed through high-throughput sequencing (HTS) using total RNAs as templates and the Illumina HiSeq 2500 platform. Analysis of HTS data revealed the presence of five viruses and three viroids in the infected vine. These viruses are Grapevine leafroll-associated virus 1 (GLRaV-1) and GLRaV-3 (genus Ampelovirus, family Closteroviridae) and three viruses of the family Betaflexiviridae (namely, Grapevine virus A [GVA], Grapevine virus B, and Grapevine rupestris stem pitting-associated virus [GRSPaV]). We also show that multiple distinct strains of three viruses (GLRaV-3, GVA, and GRSPaV) were present in this diseased grapevine. The complete genomes of two novel and highly divergent isolates of GLRaV-3 were determined using the draft genomes derived from HTS data and two independent rapid amplification of cDNA ends (RACE) strategies to obtain sequences at both the 5' and the 3' termini of the viral genomes. Questionable genome regions of both isolates were also verified through cloning of reverse transcription polymerase chain reaction products and Sanger sequencing. These two isolates are vastly divergent from all other isolates of GLRaV-3 whose genome sequences are available in GenBank. Isolate ON8415A has up to 76% nucleotide sequence identities to other isolates representing existing variant groups. We also revealed high degrees of variation in both length and sequence in the terminal untranslated regions (UTRs) of GLRaV-3 variants. The 5'-UTR of most GLRaV-3 isolates whose complete genomes have been sequenced contain tandem repeats of 65 nucleotides, a highly unusual feature rarely observed in (+)single-stranded RNA viruses. Mechanisms for the biogenesis of these tandem repeats and their function in virus replication and pathogenesis require investigation. Findings of this research add to the genetic diversity, evolutionary biology, and diagnostics of GLRaV-3 that afflicts the global grape wine industry.

Identification and characterization of water chestnut Soymovirus-1 (WCSV-1), a novel Soymovirus in water chestnuts (Eleocharis dulcis)

BACKGROUND

A disease of unknown etiology in water chestnut plants (Eleocharis dulcis) was reported in China between 2012 and 2014. High throughput sequencing of small RNA (sRNA) combined with bioinformatics, and molecular identification based on PCR detection with virus-specific primers and DNA sequencing is a desirable approach to identify an unknown infectious agent. In this study, we employed this approach to identify viral sequences in water chestnut plants and to explore the molecular interaction of the identified viral pathogen and its natural plant host.

RESULTS

Based on high throughput sequencing of virus-derived small RNAs (vsRNA), we identified the sequence a new-to-science double-strand DNA virus isolated from water chestnut cv. 'Tuanfeng' samples, a widely grown cultivar in Hubei province, China, and analyzed its genomic organization. The complete genomic sequence is 7535 base-pairs in length, and shares 42-52% nucleotide sequence identity with viruses in the Caulimoviridae family. The virus contains nine predicated open reading frames (ORFs) encoding nine hypothetical proteins, with conserved domains characteristic of caulimoviruses. Phylogenetic analyses at the nucleotide and amino acid levels indicated that the virus belongs to the genus Soymovirus. The virus is tentatively named Water chestnut soymovirus-1 (WCSV-1). Phylogenetic analysis of the putative viral polymerase protein suggested that WCSV-1 is distinct to other well established species in the Soymovirus genus. This conclusion was supported by phylogenetic analyses of the amino acid sequences encoded by ORFs I, IV, VI, or VII. The sRNA bioinformatics showed that the majority of the vsRNAs are 22-nt in length with a preference for U at the 5'-terminal nucleotide. The vsRNAs are unevenly distributed over both strands of the entire WCSV-1 circular genome, and are clustered into small defined regions. In addition, we detected WCSV-1 in asymptomatic and symptomatic water chestnut samples collected from different regions of China by using PCR. RNA-seq assays further confirmed the presence of WCSV-1-derived viral RNA in infected plants.

CONCLUSIONS

This is the first discovery of a dsDNA virus in the genus Soymovirus infecting water chestnuts. Data presented also add new information towards a better understanding of the co-evolutionary mechanisms between the virus and its natural plant host.

High-throughput small RNA sequencing for evaluation of grapevine sanitation efficacy

This study describes the application of high-throughput sequencing of small RNA analysis of the efficacy of using Ribavirin to eliminate Grapevine leafroll-associated virus 1, Grapevine fleck virus and Grapevine rupestris stem pitting-associated virus from Vitis vinifera cv. Riesling. The original plant used for sanitation by Ribavirin treatment was one naturally infected with all the viruses mentioned above as confirmed by RT-PCR. A tissue cultures of the plant were established and plantlets obtained were sanitized using Ribavirin. Three years after sanitation, a small RNA sequencing method for virus detection, targeting 21, 22 and 24 nt-long viral small RNAs (vsRNAs), was used to analyze both the mother plant and the sanitized plants. The results showed that the mother plant was infected by the three mentioned viruses and additionally by two viroids - Hop stunt viroid and Grapevine yellow speckle viroid 1. After Ribavirin treatment, the plants contained only the two viroids, with the complete elimination of all the viruses previously present.

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.

A Short Indel-Lacking-Resistance Gene Triggers Silencing of the Photosynthetic Machinery Components Through TYLCSV-Associated Endogenous siRNAs in Tomato

Plant viruses modify gene expression in infected tissues by altering the micro (mi)RNA-mediated regulation of genes. Among conserved miRNA targets there are transcripts coding for transcription factors, RNA silencing core, and disease-resistance proteins. Paralogs in these gene families are widely present in plant genomes and are known to respond differently to miRNA-mediated regulation during plant virus infections. Using genome-wide approaches applied to Solanum lycopersicum infected by a nuclear-replicating virus, we highlighted miRNA-mediated cleavage events that could not be revealed in virus-free systems. Among them we confirmed miR6024 targeting and cleavage of RX-coiled-coil (RX-CC), nucleotide binding site (NBS), leucine-rich (LRR) mRNA. Cleavage of paralogs was associated with short indels close to the target sites, indicating a general functional significance of indels in fine-tuning gene expression in plant-virus interaction. miR6024-mediated cleavage, uniquely in virus-infected tissues, triggers the production of several 21-22 nt secondary siRNAs. These secondary siRNAs, rather than being involved in the cascade regulation of other NBS-LRR paralogs, explained cleavages of several mRNAs annotated as defence-related proteins and components of the photosynthetic machinery. Outputs of these data explain part of the phenotype plasticity in plants, including the appearance of yellowing symptoms in the viral pathosystem.

Accumulation of 24 nucleotide transgene-derived siRNAs is associated with crinivirus immunity in transgenic plants

RNA silencing is a conserved antiviral defence mechanism that has been used to develop robust resistance against plant virus infections. Previous efforts have been made to develop RNA silencing-mediated resistance to criniviruses, yet none have given immunity. In this study, transgenic Nicotiana benthamiana plants harbouring a hairpin construct of the Lettuce infectious yellows virus (LIYV) RNA-dependent RNA polymerase (RdRp) sequence exhibited immunity to systemic LIYV infection. Deep sequencing analysis was performed to characterize virus-derived small interfering RNAs (vsiRNAs) generated on systemic LIYV infection in non-transgenic N. benthamiana plants as well as transgene-derived siRNAs (t-siRNAs) derived from the immune-transgenic plants before and after LIYV inoculation. Interestingly, a similar sequence distribution pattern was obtained with t-siRNAs and vsiRNAs mapped to the transgene region in both immune and susceptible plants, except for a significant increase in t-siRNAs of 24 nucleotides in length, which was consistent with small RNA northern blot results that showed the abundance of t-siRNAs of 21, 22 and 24 nucleotides in length. The accumulated 24-nucleotide sequences have not yet been reported in transgenic plants partially resistant to criniviruses, and thus may indicate their correlation with crinivirus immunity. To further test this hypothesis, we developed transgenic melon (Cucumis melo) plants immune to systemic infection of another crinivirus, Cucurbit yellow stunting disorder virus (CYSDV). As predicted, the accumulation of 24-nucleotide t-siRNAs was detected in transgenic melon plants by northern blot. Together with our findings and previous studies on crinivirus resistance, we propose that the accumulation of 24-nucleotide t-siRNAs is associated with crinivirus immunity in transgenic plants.

First natural crossover recombination between two distinct species of the family Closteroviridae leads to the emergence of a new disease

Lettuce chlorosis virus-SP (LCV-SP) (family Closteroviridae, genus Crinivirus), is a new strain of LCV which is able to infect green bean plants but not lettuce. In the present study, high-throughput and Sanger sequencing of RNA was used to obtain the LCV-SP full-length sequence. The LCV-SP genome comprises 8825 nt and 8672 nt long RNA1 and RNA2 respectively. RNA1 of LCV-SP contains four ORFs, the proteins encoded by the ORF1a and ORF1b are closely related to LCV RNA1 from California (FJ380118) whereas the 3´ end encodes proteins which share high amino acid sequence identity with RNA1 of Bean yellow disorder virus (BnYDV; EU191904). The genomic sequence of RNA2 consists of 8 ORFs, instead of 10 ORFs contained in LCV-California isolate. The distribution of vsiRNA (virus-derived small interfering RNA) along the LCV-SP genome suggested the presence of subgenomic RNAs corresponding with HSP70, P6.4 and P60. Results of the analysis using RDP4 and Simplot programs are the proof of the evidence that LCV-SP is the first recombinant of the family Closteroviridae by crossover recombination of intact ORFs, being the LCV RNA1 (FJ380118) and BnYDV RNA1 (EU191904) the origin of the new LCV strain. Genetic diversity values of virus isolates in the recombinant region obtained after sampling LCV-SP infected green bean between 2011 and 2017 might suggest that the recombinant virus event occurred in the area before this period. The presence of LCV-SP shows the role of recombination as a driving force of evolution within the genus Crinivirus, a globally distributed, emergent genus.

Small RNA NGS Revealed the Presence of Cherry Virus A and Little Cherry Virus 1 on Apricots in Hungary

Fruit trees, such as apricot trees, are constantly exposed to the attack of viruses. As they are propagated in a vegetative way, this risk is present not only in the field, where they remain for decades, but also during their propagation. Metagenomic diagnostic methods, based on next generation sequencing (NGS), offer unique possibilities to reveal all the present pathogens in the investigated sample. Using NGS of small RNAs, a special field of these techniques, we tested leaf samples of different varieties of apricot originating from an isolator house or open field stock nursery. As a result, we identified Cherry virus A (CVA) and little cherry virus 1 (LChV-1) for the first time in Hungary. The NGS results were validated by RT-PCR and also by Northern blot in the case of CVA. Cloned and Sanger sequenced viral-specific PCR products enabled us to investigate their phylogenetic relationships. However, since these pathogens have not been described in our country before, their role in symptom development and modification during co-infection with other viruses requires further investigation.

Use of siRNAs for Diagnosis of Viruses Associated to Woody Plants in Nurseries and Stock Collections

Woody perennial plants like grapevine and fruit trees can be infected by several viruses even as multiple infections. Since they are propagated vegetatively, the phytosanitary status of the propagation material (both the rootstock and the variety) can have a profound effect on the lifetime and health of the new plantations. The fast evolution of sequencing techniques provides a new opportunity for metagenomics-based viral diagnostics. Viral derived small RNAs produced by the host immune system during viral infection can be sequenced by next-generation techniques and analyzed for the presence of viruses, revealing the presence of all known viral pathogens in the sample. This method is based on Illumina sequencing of short RNAs and bioinformatics analysis of virus-derived small RNAs in the host. Here we describe a protocol for this challenging technique step by step with notes, in order to ensure success for every user.

Application of next generation sequencing toward sensitive detection of enteric viruses isolated from celery samples as an example of produce

Next generation sequencing (NGS) holds promise as a single application for both detection and sequence identification of foodborne viruses; however, technical challenges remain due to anticipated low quantities of virus in contaminated food. In this study, with a focus on data analysis using several bioinformatics tools, we applied NGS toward amplification-independent detection and identification of norovirus at low copy (<10³ copies) or within multiple strains from produce. Celery samples were inoculated with human norovirus (stool suspension) either as a single norovirus strain, a mixture of strains (GII.4 and GII.6), or a mixture of different species (hepatitis A virus and norovirus). Viral RNA isolation and recovery was confirmed by RT-qPCR, and optimized for library generation and sequencing without amplification using the Illumina MiSeq platform. Extracts containing either a single virus or a two-virus mixture were analyzed using two different analytic approaches to achieve virus detection and identification. First an overall assessment of viral genome coverage for samples varying in copy numbers (1.1×10³ to 1.7×10⁷) and genomic content (single or multiple strains in various ratios) was completed by reference-guided mapping. Not unexpectedly, this targeted approach to identification was successful in correctly mapping reads, thus identifying each virus contained in the inoculums even at low copy (estimated at 12 copies). For the second (metagenomic) approach, samples were treated as "unknowns" for data analyses using (i) a sequence-based alignment with a local database, (ii) an "in-house" k-mer tool, (iii) a commercially available metagenomics bioinformatic analysis platform cosmosID, and (iv) an open-source program Kraken. Of the four metagenomics tools applied in this study, only the local database alignment and in-house k-mer tool were successful in detecting norovirus (as well as HAV) at low copy (down to <10³ copies) and within a mixture of virus strains or species. The results of this investigation provide support for continued investigation into the development and integration of these analytical tools for identification and detection of foodborne viruses.

Deep sequencing analysis reveals a TMV mutant with a poly(A) tract reduces host defense responses in Nicotiana benthamiana

Tobacco mosaic virus (TMV) possesses an upstream pseudoknotted domain (UPD), which is important for replication. After substituting the UPD with an internal poly(A) tract (43 nt), a mutant TMV-43A was constructed. TMV-43A replicated slower than TMV and induced a non-lethal mosaic symptom in Nicotiana benthamiana. In this study, deep sequencing was performed to detect the differences of small RNA profiles between TMV- and TMV-43A-infected N. benthamiana. The results showed that TMV-43A produced lesser amount of virus-derived interfering RNAs (vsiRNAs) than that of TMV. However, the distributions of vsiRNAs generation hotspots between TMV and TMV-43A were similar. Expression of genes related to small RNA biogenesis in TMV-43A-infected N. benthamiana was significantly lower than that of TMV, which leads to generation of lesser vsiRNAs. The expressions of host defense response genes were up-regulated after TMV infection, as compared to TMV-43A-infected plants. Host defense response to TMV-43A infection was lower than that to TMV. The absence of UPD might contribute to the reduced host response to TMV-43A. Our study provides valuable information in the role of the UPD in eliciting host response genes after TMV infection in N. benthamiana. (187 words).

Nicotiana benthamiana plants asymptomatically infected by Pelargonium line pattern virus show unusually high accumulation of viral small RNAs that is neither associated with DCL induction nor RDR6 activity

Pelargonium line pattern virus (PLPV, Tombusviridae) normally establishes systemic, low-titered and asymptomatic infections in its hosts. This type of interaction may be largely determined by events related to RNA silencing, a major antiviral mechanism in plants. This mechanism is triggered by double or quasi double-stranded (ds) viral RNAs which are cut by DCL ribonucleases into virus small RNAs (vsRNAs). Such vsRNAs are at the core of the silencing process as they guide sequence-specific RNA degradation Host RNA dependent-RNA polymerases (RDRs), and particularly RDR6, strengthen antiviral silencing by promoting biosynthesis of secondary vsRNAs. To approach PLPV-host relationship, here we have characterized the vsRNAs that accumulate in PLPV-infected Nicotiana benthamiana. Such accumulation was found unprecedented high despite DCLs were not induced in infected tissue and neither vsRNA generation nor PLPV infection was apparently affected by RDR6 impairment. From the obtained data, triggers and host factors likely involved in anti-PLPV silencing are proposed.

Differential expression of miRNAs and associated gene targets in grapevine leafroll-associated virus 3-infected plants

MicroRNAs (miRNAs) are a class of endogenous small non-coding RNAs (sRNA) that play an essential role in the regulation of target mRNAs expressed during plant development and in response to stress. MicroRNA expression profiling has helped to identify miRNAs that regulate a range of processes, including the plant's defence response to pathogens. In this study, differential miRNA expression in own-rooted Vitis vinifera cv. Cabernet Sauvignon plants infected with grapevine leafroll-associated virus 3 was investigated with microarrays and next-generation sequencing (NGS) of sRNA and mRNA. These high-throughput approaches identified several differentially expressed miRNAs. Four miRNAs, identified by both approaches, were validated by stemloop RT-PCRs. Three of the predicted targets of the differentially expressed miRNAs were also differentially expressed in the transcriptome data of infected plants, and were validated by RT-qPCR. Identification of these miRNAs and their targets can lead to a better understanding of host-pathogen interactions involved in grapevine leafroll disease and the identification of possible targets for virus resistance.

Comprehensive Virus Detection Using Next Generation Sequencing in Grapevine Vascular Tissues of Plants Obtained from the Wine Regions of Bohemia and Moravia (Czech Republic)

Comprehensive next generation sequencing virus detection was used to detect the whole spectrum of viruses and viroids in selected grapevines from the Czech Republic. The novel NGS approach was based on sequencing libraries of small RNA isolated from grapevine vascular tissues. Eight previously partially-characterized grapevines of diverse varieties were selected and subjected to analysis: Chardonnay, Laurot, Guzal Kara, and rootstock Kober 125AA from the Moravia wine-producing region; plus Müller-Thurgau and Pinot Noir from the Bohemia wine-producing region, both in the Czech Republic. Using next generation sequencing of small RNA, the presence of 8 viruses and 2 viroids were detected in a set of eight grapevines; therefore, confirming the high effectiveness of the technique in plant virology and producing results supporting previous data on multiple infected grapevines in Czech vineyards. Among the pathogens detected, the Grapevine rupestris vein feathering virus and Grapevine yellow speckle viroid 1 were recorded in the Czech Republic for the first time.

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.

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.

Deep sequencing of banana bract mosaic virus from flowering ginger (Alpinia purpurata) and development of an immunocapture RT-LAMP detection assay

Banana bract mosaic virus (BBrMV) has never been reported in banana plants in Hawaii. In 2010, however, it was detected in a new host, flowering ginger (Alpinia purpurata). In this study, we characterize the A. purpurata isolate and study its spread in flowering ginger in Hawaii. A laboratory study demonstrated that BBrMV could be transmitted from flowering ginger to its natural host, banana, therefore raising a serious concern about the potential risk to the rapidly growing banana industry of Hawaii. To quickly monitor this virus in the field, we developed a robust immunocapture reverse transcription loop-mediated isothermal amplification (IC-RT-LAMP) assay. Deep sequencing of the BBrMV isolate from A. purpurata revealed a single-stranded RNA virus with a genome of 9,713 nt potentially encoding a polyprotein of 3,124 aa, and another predicted protein, PIPO, in the +2 reading-frame shift. Most of the functional motifs in the Hawaiian isolate were conserved among the genomes of isolates from one found in the Philippines and India. However, the A. purpurata isolate had an amino acid deletion in the Pl protein that was most similar to the Philippine isolate. Phylogenetic analysis of an eastern Pacific subpopulation that included A. purpurata was closest in genetic distance to a Southeast Asian subpopulation, suggesting frequent gene flow and supporting the hypothesis that the A. purpurata isolate arrived in Hawaii from Southeast Asia.

Generation of a high resolution map of sRNAs from Fusarium graminearum and analysis of responses to viral infection

Previously, we characterized F. graminearum hypovirus 1 (FgHV1) and F. graminearum hypovirus 2 (FgHV2), which are the only two hypoviruses in F. graminearum that are closely related to Cryphonectria hypovirus 1 (CHV1) and Cryphonectria hypovirus 2 (CHV2) in the Hypoviridae family. In this study, we preliminarily elucidated the RNA silencing mechanism of the F. graminearum/hypovirus system from a small RNA (sRNA) perspective by using HiSeq deep sequencing. The length distributions of F. graminearum sRNA were altered by hypoviral infection. Potential microRNA-like (milRNA) candidates were differentially expressed between the hypovirus-free and hypovirus-infected library types. Extensive virus-derived small interfering RNAs (vsiRNAs) were also principally defined. The 1,831,081 and 3,254,758 total reads generated from the FgHV1 and FgHV2 genomes in F. graminearum yielded the first high-resolution sRNA maps of fungal viruses. In addition, extensive bioinformatics searches identified a large number of transcripts that are potentially targeted by vsiRNAs, several of which were effectively down-regulated. In particular, the RNA silencing-related genes FgDicer1 and FgRdRp5 were predicted targets of FgHV1- and FgHV2-derived siRNAs, possibly revealing a novel anti-RNA silencing strategy employed by mycoviruses.

Mechanisms, applications, and perspectives of antiviral RNA silencing in plants

Viral diseases of plants cause important economic losses due to reduction in crop quality and quantity to the point of threatening food security in some countries. Given the reduced availability of natural sources, genetic resistance to viruses has been successfully engineered for some plant-virus combinations. A sound understanding of the basic mechanisms governing plant-virus interactions, including antiviral RNA silencing, is the foundation to design better management strategies and biotechnological approaches to engineer and implement antiviral resistance in plants. In this review, we present current molecular models to explain antiviral RNA silencing and its application in basic plant research, biotechnology and genetic engineering.

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.

In vivo study of Dicer-2-mediated immune response of the small interfering RNA pathway upon systemic infections of virulent and avirulent viruses in Bombus terrestris

Recent studies suggest a potent role of the small interfering RNA (siRNA) pathway in the control of bee viruses and its usefulness to tackle these viral diseases. However, the involvement of the siRNA pathway in the defense against different bee viruses is still poorly understood. Therefore, in this report, we comprehensively analyzed the response of the siRNA pathway in bumblebees of Bombus terrestris to systemic infections of the virulent Israeli acute paralysis virus (IAPV) and the avirulent slow bee paralysis virus (SBPV). Our results showed that IAPV and SBPV infections induced the expression of Dicer-2. IAPV infections also triggered the production of predominantly 22 nt-long virus-derived siRNAs (vsiRNAs). Intriguingly, these 22 nt-long vsiRNAs showed a high proportion of antigenomic IAPV sequences. Conversely, these predominantly 22 nt-long vsiRNAs of SBPV were not detected in SBPV infected bees. Furthermore, an "RNAi-of-RNAi" experiment on Dicer-2 did not result in altered genome copy numbers of IAPV (n = 17-18) and also not of SBPV (n = 11-12). Based on these results, we can speculate about the importance of the siRNA pathway in bumblebees for the antiviral response. During infection of IAPV, this pathway is probably recruited but it might be insufficient to control viral infection in our experimental setup. The host can control SBPV infection, but aside from the induction of Dicer-2 expression, no further evidence of the antiviral activity of the siRNA pathway was observed. This report may also enhance the current understanding of the siRNA pathway in the innate immunity of non-model insects upon different viral infections.

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.

A highly effective and versatile technology for the isolation of RNAs from grapevines and other woody perennials for use in virus diagnostics

BACKGROUND

Isolation of pure RNA from woody perennials, especially fruit crops such as grapevine rich in complex secondary metabolites, has remained very challenging. Lack of effective RNA isolation technology has resulted in difficulties in viral diagnosis and discovery as well as studies on many biological processes of these highly important woody plants. It is imperative to develop and refine methodologies with which large amounts of pure nucleic acids can be readily isolated from woody perennials.

METHODS

We compared five commonly used RNA isolation kits in isolating total RNA from twelve species of woody perennials. We made modifications to select RNA isolation systems to simplify and improve their efficiency in RNA isolation. The yield and quality of isolated RNAs were assessed via gel electrophoresis and spectrophotometric measurement. We also performed RT-PCR and RT-qPCR to detect several major viruses from grapevines.

RESULTS

Two of the kits were shown to be the best in both the yield and quality of the isolated RNA from all twelve woody species. Using disposable extraction bags for tissue homogenization not only improved the yield without affecting quality, but also made the RNA isolation technology simpler, less costly, and suitable for adoption by many potential users with facility limitations. This system was successfully applied to a wide range of woody plants, including fruit crops, ornamentals and timber trees. Inclusion of polyvinylpyrrolidone in the extraction buffer drastically improved the performance of the system in isolating total RNA from old grapevine leaves collected later in the season. This modification made our system highly effective in isolating quality RNA from grapevine leaves throughout the entire growing season. We further demonstrated that the resulting nucleic acid preparations are suitable for detection of several major grapevine viruses with RNA or DNA genomes using PCR, RT-PCR and qPCR as well as for assays on plant microRNAs.

CONCLUSIONS

This improved RNA isolation system would have wide applications in viral diagnostics and discovery, studies on gene expression and regulation, transcriptomics, and small RNA biology in grapevines. We believe this system will also be useful in diverse applications pertaining to research on many other woody perennials and recalcitrant plant species.

Small RNA profiles of wild-type and silencing suppressor-deficient tomato spotted wilt virus infected Nicotiana benthamiana

Tospoviruses are plant-infecting viruses belonging to the family Bunyaviridae. We used a collection of wild-type, phylogenetically distinct tomato spotted wilt virus isolates and related silencing-suppressor defective mutants to study the effects on the small RNA (sRNA) accumulation during infection of Nicotiana benthamiana. Our data showed that absence of a functional silencing suppressor determined a marked increase of the total amount of viral sRNAs (vsRNAs), and specifically of the 21 nt class. We observed a common under-representation of vsRNAs mapping to the intergenic region of S and M genomic segments, and preferential mapping of the reads against the viral sense open reading frames, with the exception of the NSs gene. The NSs-mutant strains showed enrichment of NSm-derived vsRNA compared to the expected amount based on gene size. Analysis of 5' terminal nucleotide preference evidenced a significant enrichment in U for the 21 nt- and in A for 24 nt-long endogenous sRNAs in all the samples. Hotspot analysis revealed a common abundant accumulation of reads at the 5' end of the L segment, mostly in the antiviral sense, for the NSs-defective isolates, suggesting that absence of the silencing suppressor can influence preferential targeting of the viral genome.

Detection of a New Luteovirus in Imported Nectarine Trees: A Case Study to Propose Adoption of Metagenomics in Post-Entry Quarantine

In spring 2013, 5-year-old nectarine (Prunus persica) trees, grafted on peach rootstock Nemaguard, were found stunted in a propagation block in California. These trees had been propagated from budwood of three nectarine cultivars imported from France and cleared through the post-entry quarantine procedure. Examination of the canopy failed to reveal any obvious symptoms. However, examination of the trunks, after stripping the bark, revealed extensive pitting on the woody cylinder. To investigate the etiological agent, double-stranded RNA was extracted from bark scrapings from the scion and rootstock portions, and a cDNA library was prepared and sequenced using the Illumina platform. BLAST analysis of the contigs generated by the de novo assembly of sequence reads indicated the presence of a novel luteovirus. Complete sequence of the viral genome was determined by sequencing of three overlapping cDNA clones generated by reverse transcription-polymerase chain reaction (RT-PCR) and by rapid amplification of the 5'- and 3'-termini. The virus genome was comprised of 4,991 nucleotides with a gene organization similar to members of the genus Luteovirus (family Luteoviridae). The presence of the virus, tentatively named Nectarine stem pitting-associated virus, was confirmed in symptomatic trees by RT-PCR. Discovery of a new virus in nectarine trees after post-entry quarantine indicates the importance of including (i) metagenomic analysis by next-generation sequencing approach as an essential tool to assess the plant health status, and (ii) examination of the woody cylinders as part of the indexing process.

Detection and molecular characterisation of Grapevine Syrah virus-1 isolates from Central Europe

Grapevine Syrah virus-1 (GSyV-1) was identified by small-RNA deep sequencing in Slovak grapevine co-infected by several other viruses. The RT-PCR assays developed in this work substantially improved the virus detection and allowed the identification of GSyV-1 in tested grapevine samples from Slovakia and the Czech Republic at an unexpectedly high rate (ca. 30 %). Subsequently, complete genome sequences of 3 GSyV-1 isolates (2 Slovak and 1 Czech) were determined by Sanger sequencing, showing a typical marafivirus genome organization. Analyses of complete genome sequences showed a higher intra-group diversity among these 3 central European GSyV-1 isolates (differences reaching 7.1 % at the nucleotide level) in comparison to 3 previously characterized North American isolates (only 1.2 % intra-group divergence). A substantially higher divergence among central European isolates and their clustering into two major phylogenetic groups was further confirmed by the partial genome analysis of additional 26 isolates. The CP-centered study did not support the geography-based clustering among central European and American isolates. Nevertheless, the sequence data of the highly variable 5'-proximal portion of the genome obtained for few additional isolates from Slovakia and Czech Republic showed the presence of both, "European-" and "north American-like", GSyV-1 isolates in the analyzed grapevine samples.

Genetic Variability of Grapevine Pinot gris virus and Its Association with Grapevine Leaf Mottling and Deformation

The role of Grapevine Pinot gris virus (GPGV) in the etiology of grapevine leaf mottling and deformation was investigated by biological and molecular assays. A survey on different cultivars from the Trentino Region in Italy showed a widespread distribution of GPGV, which was associated with symptomatic (79%) but also with symptomless (21%) vines. Symptomatic and GPGV-infected 'Pinot gris' vines induced symptoms on grafted vines of healthy Pinot gris or 'Traminer', whereas GPGV-infected but symptomless vines did not. High-throughput sequencing of small RNA (sRNA) populations of two infected Pinot gris accessions confirmed the existence of nearly overlapping viromes in vines with or without symptoms but phylogenetic analyses of the genomes of seven GPGV isolates from Italy and the Czech and Slovak Republics clearly differentiated those infecting symptomatic vines. The involvement of Grapevine rupestris vein feathering virus (GRVFV) in the disease, which was only infecting the symptomatic vine, was ruled out by reverse-transcription polymerase chain reaction studies. Maximum likelihood and Bayesian phylogenetic analysis of two GPGV genomic regions, encompassing part of the movement protein (MP) and coat protein gene sequences and the RNA-dependent RNA polymerase domain of the replicase gene, showed that isolates from symptomatic vines form a lineage distinct from that of symptomless vines. Moreover, the presence or lack of the MP stop codon identified in viral isolates from symptomatic or symptomless vines, respectively, is likely responsible for an MP six amino acids longer in symptomless isolates.

A novel grapevine badnavirus is associated with the Roditis leaf discoloration disease

Roditis leaf discoloration (RLD), a graft-transmissible disease of grapevine, was first reported in Greece in the 1980s. Even though various native grapevine viruses were identified in the affected vines, the etiology of the disease remained unknown. In the present study, we used an NGS platform for sequencing siRNAs from a twenty-year old Roditis vine showing typical RLD symptoms. Analysis of the NGS data revealed the presence of various known grapevine viruses and viroids as well as a hitherto uncharacterized DNA virus. The circular genome of the new virus was fully reassembled. It is 6988 nts long and includes 4 open reading frames (ORFs). ORF1, ORF2 and ORF4 code for proteins with unknown functions while ORF3 encodes a polyprotein with motifs related to the replication, encapsidation and movement of the virus. Phylogenetic analysis classified the novel virus within the genus Badnavirus, with closest relationship to Fig badnavirus 1. Further studies showed that the new badnavirus is closely related with the RLD disease and the provisional name grapevine Roditis leaf discoloration-associated virus (GRLDaV) is proposed. Our findings extend the number of DNA viruses identified in grapevine, further drawing attention to the potential importance of this virus group on grapevine pathology.

Community-onset Klebsiella pneumoniae pneumonia in Taiwan: clinical features of the disease and associated microbiological characteristics of isolates from pneumonia and nasopharynx

As virus diseases cannot be controlled by traditional plant protection methods, the risk of their spread have to be minimized on vegetatively propagated plants, such as grapevine. Metagenomic approaches used for virus diagnostics offer a unique opportunity to reveal the presence of all viral pathogens in the investigated plant, which is why their application can reduce the risk of using infected material for a new plantation. Here we used a special branch, deep sequencing of virus-derived small RNAs, of this high-throughput method for virus diagnostics, and determined viromes of vineyards in Hungary. With NGS of virus-derived small RNAs we could detect not only the viruses tested routinely, but also new ones, which had never been described in Hungary before. Virus presence did not correlate with the age of the plantation, moreover phylogenetic analysis of the identified virus isolates suggests that infections are mostly caused by the use of infected propagating material. Our results, validated by other molecular methods, raised further questions to be answered before this method can be introduced as a routine, reliable test for grapevine virus diagnostics.

Comparative analysis among the small RNA populations of source, sink and conductive tissues in two different plant-virus pathosystems

Background

In plants, RNA silencing plays a fundamental role as defence mechanism against viruses. During last years deep-sequencing technology has allowed to analyze the sRNA profile of a large variety of virus-infected tissues. Nevertheless, the majority of these studies have been restricted to a unique tissue and no comparative analysis between phloem and source/sink tissues has been conducted. In the present work, we compared the sRNA populations of source, sink and conductive (phloem) tissues in two different plant virus pathosystems. We chose two cucurbit species infected with two viruses very different in genome organization and replication strategy; Melon necrotic spot virus (MNSV) and Prunus necrotic ringspot virus (PNRSV).

Results

Our findings showed, in both systems, an increase of the 21-nt total sRNAs together with a decrease of those with a size of 24-nt in all the infected tissues, except for the phloem where the ratio of 21/24-nt sRNA species remained constant. Comparing the vsRNAs, both PNRSV- and MNSV-infected plants share the same vsRNA size distribution in all the analyzed tissues. Similar accumulation levels of sense and antisense vsRNAs were observed in both systems except for roots that showed a prevalence of (+) vsRNAs in both pathosystems. Additionally, the presence of overrepresented discrete sites along the viral genome, hot spots, were identified and validated by stem-loop RT-PCR. Despite that in PNRSV-infected plants the presence of vsRNAs was scarce both viruses modulated the host sRNA profile.

Conclusions

We compare for the first time the sRNA profile of four different tissues, including source, sink and conductive (phloem) tissues, in two plant-virus pathosystems. Our results indicate that antiviral silencing machinery in melon and cucumber acts mainly through DCL4. Upon infection, the total sRNA pattern in phloem remains unchanged in contrast to the rest of the analyzed tissues indicating a certain tissue-tropism to this polulation. Independently of the accumulation level of the vsRNAs both viruses were able to modulate the host sRNA pattern.

Electronic supplementary material

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

Community-onset Klebsiella pneumoniae pneumonia in Taiwan: clinical features of the disease and associated microbiological characteristics of isolates from pneumonia and nasopharynx

Klebsiella pneumoniae is an important cause of community-onset pneumonia in Asian countries and South Africa. We investigated the clinical characteristics of K. pneumoniae causing community-onset pneumonia, and the associated microbiological features between K. pneumoniae isolates from pneumonia and those from the nasopharynx in Taiwan. This study was conducted at the Taipei Veterans General Hospital during July, 2012 to February, 2014. The clinical characteristics in patients with community-onset K. pneumoniae pneumonia were analyzed. K. pneumoniae isolates from the nasopharynx of adults attending otorhinolaryngology outpatient clinics were collected to compare their microbiological features with those from pneumonia. Capsular genotypes, antimicrobial susceptibility, and multilocus sequence type (MLST) were determined among these strains. Ninety-one patients with community-onset K. pneumoniae pneumonia were enrolled. We found a high mortality (29.7%) among these patients. Capsular types K1, K2, K5, K20, K54, and K57 accounted for ∼70% of the K. pneumoniae isolates causing pneumonia, and ∼70% of all the K. pneumoniae strains isolated from the nasopharynx of patients in outpatient clinics. The MLST profiles further demonstrated the genetic relatedness between most pneumonia isolates and those from the nasopharynx. In conclusion, our results show that community-onset pneumonia caused by K. pneumoniae was associated with high mortality and could have a reservoir in the nasopharynx. To tackle this high-mortality disease, the distribution of capsular types in the nasopharynx might have implications for future vaccine development.

Bioinformatics approaches for viral metagenomics in plants using short RNAs: model case of study and application to a Cicer arietinum population

Over the past years deep sequencing experiments have opened novel doors to reconstruct viral populations in a high-throughput and cost-effective manner. Currently a substantial number of studies have been performed which employ next generation sequencing techniques to either analyze known viruses by means of a reference-guided approach or to discover novel viruses using a de novo-based strategy. Taking advantage of the well-known Cymbidium ringspot virus we have carried out a comparison of different bioinformatics tools to reconstruct the viral genome based on 21–27 nt short (s)RNA sequencing with the aim to identify the most efficient pipeline. The same approach was applied to a population of plants constituting an ancient variety of Cicer arietinum with red seeds. Among the discovered viruses, we describe the presence of a Tobamovirus referring to the Tomato mottle mosaic virus (NC_022230), which was not yet observed on C. arietinum nor revealed in Europe and a viroid referring to Hop stunt viroid (NC_001351.1) never reported in chickpea. Notably, a reference sequence guided approach appeared the most efficient in such kind of investigation. Instead, the de novo assembly reached a non-appreciable coverage although the most prominent viral species could still be identified. Advantages and limitations of viral metagenomics analysis using sRNAs are discussed.

Characterization by Deep Sequencing of Prunus virus T, a Novel Tepovirus Infecting Prunus Species

Double-stranded RNAs purified from a cherry tree collected in Italy and a plum tree collected in Azerbaijan were submitted to deep sequencing. Contigs showing weak but significant identity with various members of the family Betaflexiviridae were reconstructed. Sequence comparisons led to the conclusion that the viral isolates identified in the analyzed Prunus plants belong to the same viral species. Their genome organization is similar to that of some members of the family Betaflexiviridae, with three overlapping open reading frames (RNA polymerase, movement protein, and capsid protein). Phylogenetic analyses of the deduced encoded proteins showed a clustering with the sole member of the genus Tepovirus, Potato virus T (PVT). Given these results, the name Prunus virus T (PrVT) is proposed for the new virus. It should be considered as a new member of the genus Tepovirus, even if the level of nucleotide identity with PVT is borderline with the genus demarcation criteria for the family Betaflexiviridae. A reverse-transcription polymerase chain reaction detection assay was developed and allowed the identification of two other PrVT isolates and an estimate of 1% prevalence in the large Prunus collection screened. Due to the mixed infection status of all hosts identified to date, it was not possible to correlate the presence of PrVT with specific symptoms.

Drawing siRNAs of viral origin out from plant siRNAs libraries

Viruses are obligate intracellular entities that infect all forms of life. In plants, invading viral nucleic acids trigger RNA silencing machinery and it results in the accumulation of viral short interfering RNAs (v-siRNAs). The study of v-siRNAs population in biological samples has become a major part of many research projects aiming to identify viruses infecting them, including unknown viruses, even at extremely low titer. Currently, siRNA populations are investigated by high-throughput sequencing approaches, which generate very large data sets. The major difficulty in these studies is to properly analyze such huge amount of data. In this regard, easy-to-use bioinformatics tools to groom and decipher siRNA libraries and to draw out v-siRNAs are needed. Here we describe a workflow, which permit users with little experience in bioinformatics to draw out v-siRNAs from raw data sequences obtained by Illumina technology. Such pipeline has been released in the context of Galaxy, an open source Web-based platform for bioinformatics analyses.

Grapevine leafroll disease and associated viruses: a unique pathosystem

Grapevine leafroll is the most complex and intriguing viral disease of grapevine (Vitis spp.). Several monopartite closteroviruses (family Closteroviridae) from grapevines have been molecularly characterized, yet their role in disease etiology is not completely resolved. Hence, these viruses are currently designated under the umbrella term of Grapevine leafroll-associated viruses (GLRaVs). This review examines our current understanding of the genetically divergent GLRaVs and highlights the emerging picture of several unique aspects of the leafroll disease pathosystem. A systems biology approach using contemporary technologies in molecular biology, -omics, and cell biology aids in exploring the comparative molecular biology of GLRaVs and deciphering the complex network of host-virus-vector interactions to bridge the gap between genomics and phenomics of leafroll disease. In addition, grapevine-infecting closteroviruses have a great potential as designer viruses to pursue functional genomics and for the rational design of novel disease intervention strategies in this agriculturally important perennial fruit crop.

Identification of viruses and viroids by next-generation sequencing and homology-dependent and homology-independent algorithms

A fast, accurate, and full indexing of viruses and viroids in a sample for the inspection and quarantine services and disease management is desirable but was unrealistic until recently. This article reviews the rapid and exciting recent progress in the use of next-generation sequencing (NGS) technologies for the identification of viruses and viroids in plants. A total of four viroids/viroid-like RNAs and 49 new plant RNA and DNA viruses from 18 known or unassigned virus families have been identified from plants since 2009. A comparison of enrichment strategies reveals that full indexing of RNA and DNA viruses as well as viroids in a plant sample at single-nucleotide resolution is made possible by one NGS run of total small RNAs, followed by data mining with homology-dependent and homology-independent computational algorithms. Major challenges in the application of NGS technologies to pathogen discovery are discussed.

Diagnosis and discovery of fungal viruses using deep sequencing of small RNAs

Analysis of virus-derived small RNAs with high-throughput sequencing has been successful for detecting novel viruses in plants and invertebrates. However, the applicability of this method has not been demonstrated in fungi, although fungi were among the first organisms reported to utilize RNA silencing. Here, we used virus-infected isolates of the fungal species complex Heterobasidion annosum sensu lato as a model system to test whether mycovirus genome segments can be detected with small RNA deep sequencing. Species of the genus Heterobasidion are some of the most devastating forest pathogens in boreal forests. These fungi cause wood decay and are commonly infected with species of the family Partitiviridae and the unassigned virus species Heterobasidion RNA virus 6. Small RNA deep sequencing allowed the simultaneous detection of all eight double-stranded RNA virus strains known to be present in the tested samples and one putative mitovirus species (family Narnaviridae) with a single-stranded RNA genome, designated here as Heterobasidion mitovirus 1. Prior to this study, no members of the family Narnaviridae had been described as infecting species of Heterobasidion. Quantification of viral double- and single-stranded RNA with quantitative PCR indicated that co-infecting viral species and viruses with segmented genomes can be detected with small RNA deep sequencing despite vast differences in the amount of RNA. This is the first study demonstrating the usefulness of this method for detecting fungal viruses. Moreover, the results suggest that viral genomes are processed into small RNAs by different species of Heterobasidion.

Viruses in maize and Johnsongrass in southern Ohio

The two major U.S. maize viruses, Maize dwarf mosaic virus (MDMV) and Maize chlorotic dwarf virus (MCDV), emerged in southern Ohio and surrounding regions in the 1960s and caused significant losses. Planting resistant varieties and changing cultural practices has dramatically reduced virus impact in subsequent decades. Current information on the distribution, diversity, and impact of known and potential U.S. maize disease-causing viruses is lacking. To assess the current reservoir of viruses present at the sites of past disease emergence, we used a combination of serological testing and next-generation RNA sequencing approaches. Here, we report enzyme-linked immunosorbent assay and RNA-Seq data from samples collected over 2 years to assess the presence of viruses in cultivated maize and an important weedy reservoir, Johnsongrass (Sorghum halepense). Results revealed a persistent reservoir of MDMV and two strains of MCDV in Ohio Johnsongrass. We identified sequences of several other grass-infecting viruses and confirmed the presence of Wheat mosaic virus in Ohio maize. Together, these results provide important data for managing virus disease in field corn and sweet corn maize crops, and identifying potential future virus threats.

Construction of a synthetic infectious cDNA clone of Grapevine Algerian latent virus (GALV-Nf) and its biological activity in Nicotiana benthamiana and grapevine plants

BACKGROUND

Grapevine Algerian latent virus (GALV) is a tombusvirus first isolated in 1989 from an Algerian grapevine (Vitis spp.) plant and more recently from water samples and commercial nipplefruit and statice plants. No further reports of natural GALV infections in grapevine have been published in the last two decades, and artificial inoculations of grapevine plants have not been reported. We developed and tested a synthetic GALV construct for the inoculation of Nicotiana benthamiana plants and different grapevine genotypes to investigate the ability of this virus to infect and spread systemically in different hosts.

METHODS

We carried out a phylogenetic analysis of all known GALV sequences and an epidemiological survey of grapevine samples to detect the virus. A GALV-Nf clone under the control of the T7 promoter was chemically synthesized based on the full-length sequence of the nipplefruit isolate GALV-Nf, the only available sequence at the time the project was conceived, and the infectious transcripts were tested in N. benthamiana plants. A GALV-Nf-based binary vector was then developed for the agroinoculation of N. benthamiana and grapevine plants. Infections were confirmed by serological and molecular analysis and the resulting ultrastructural changes were investigated in both species.

RESULTS

Sequence analysis showed that the GALV coat protein is highly conserved among diverse isolates. The first epidemiological survey of cDNAs collected from 152 grapevine plants with virus-like symptoms did not reveal the presence of GALV in any of the samples. The agroinoculation of N. benthamiana and grapevine plants with the GALV-Nf binary vector promoted efficient infections, as revealed by serological and molecular analysis. The GALV-Nf infection of grapevine plants was characterized in more detail by inoculating different cultivars, revealing distinct patterns of symptom development. Ultrastructural changes induced by GALV-Nf in N. benthamiana were similar to those induced by tombusviruses in other hosts, but the cytopathological alterations in grapevine plants were less severe.

CONCLUSIONS

This is the first report describing the development of a synthetic GALV-Nf cDNA clone, its artificial transmission to grapevine plants and the resulting symptoms and cytopathological alterations.