Quality Assessment and Validation of High-Throughput Sequencing for Grapevine Virus Diagnostics.

Development of High-Throughput Sequencing (HTS), also known as next generation sequencing, revolutionized diagnostic research of plant viruses. HTS outperforms bioassays and molecular diagnostic assays that are used to screen domestic and quarantine grapevine materials in data throughput, cost, scalability, and detection of novel and highly variant virus species. However, before HTS-based assays can be routinely used for plant virus diagnostics, performance specifications need to be developed and assessed. In this study, we selected 18 virus-infected grapevines as a test panel for measuring performance characteristics of an HTS-based diagnostic assay. Total nucleic acid (TNA) was extracted from petioles and dormant canes of individual samples and constructed libraries were run on Illumina NextSeq 500 instrument using a 75-bp single-end read platform. Sensitivity was 98% measured over 264 distinct virus and viroid infections with a false discovery rate (FDR) of approximately 1 in 5 positives. The results also showed that combining a spring petiole test with a fall cane test increased sensitivity to 100% for this TNA HTS assay. To evaluate extraction methodology, these results were compared to parallel dsRNA extractions. In addition, in a more detailed dilution study, the TNA HTS assay described here consistently performed well down to a dilution of 5%. In that range, sensitivity was 98% with a corresponding FDR of approximately 1 in 5. Repeatability and reproducibility were assessed at 99% and 93%, respectively. The protocol, criteria, and performance levels described here may help to standardize HTS for quality assurance and accreditation purposes in plant quarantine or certification programs.

Response of bitter and sweet Chenopodium quinoa varieties to cucumber mosaic virus: Transcriptome and small RNASeq perspective.

Saponins are secondary metabolites with antiviral properties. Low saponin (sweet) varieties of quinoa (Chenopodium quinoa) have been developed because seeds high in saponins taste bitter. The aim of this study was to elucidate the role of saponin in resistance of quinoa to Cucumber mosaic virus (CMV). Differential gene expression was studied in time-series study of CMV infection. High-throughput transcriptome sequence data were obtained from 36 samples (3 varieties × +/- CMV × 1 or 4 days after inoculation × 3 replicates). Translation, lipid, nitrogen, amino acid metabolism, and mono- and sesquiterpenoid biosynthesis genes were upregulated in CMV infections. In 'Red Head' (bitter), CMV-induced systemic symptoms were concurrent with downregulation of a key saponin biosynthesis gene, TSARL1, four days after inoculation. In local lesion responses (sweet and semi-sweet), TSARL1 levels remained up-regulated. Known microRNAs (miRNA) (81) from 11 miR families and 876 predicted novel miRNAs were identified. Differentially expressed miRNA and short interfering RNA clusters (24nt) induced by CMV infection are predicted to target genomic and intergenic regions enriched in repetitive elements. This is the first report of integrated RNASeq and sRNASeq data in quinoa-virus interactions and provides comprehensive understanding of involved genes, non-coding regions, and biological pathways in virus resistance.

First Report of Apricot vein clearing-associated virus Infecting flowering apricot (Prunus mume) in the United States.

Apricot vein clearing-associated virus is the type species of genus Prunevirus, family Betaflexiviridae. The virus was first discovered from an Italian apricot tree (Prunus armeniaca) showing leaf vein clearing and mottling symptoms (Elbeaino et al. 2014). Since then, apricot vein clearing-associated virus (AVCaV) has been reported in symptomatic and asymptomatic plants from other countries (Marais et al. 2015; Kinoti et al. 2017; Kubaa et al. 2014). In 2018, a domestic selection of a flowering apricot (P. mume cv. Peggy Clarke) (PC01) with no discernible foliar virus-like symptoms was received for inclusion in the Foundation Plant Services (UC-Davis) collection. The plant originated from a private Prunus collection located in California. Total nucleic acids (TNA) were isolated from PC01 leaves using MagMax Plant RNA Isolation Kit (Thermo Fisher Scientific). The TNA were analyzed for a panel of 15 Prunus-infecting viruses by reverse-transcription quantitative PCR (RT-qPCR) (Diaz-Lara et al. 2020). In addition, to screen for sap-transmissible viruses, young leaves of PC01 were homogenized in inoculation buffer and were rubbed onto leaves of herbaceous indicator plants, Chenopodium amaranticolor, C. quinoa, Cucumis sativus, and Nicotiana clevelandii (Rowhani et al. 2005). The source PC01 tested negative for the 15 screened viruses. Interestingly, vein clearing symptoms were observed on leaves of C. quinoa and C. amaranticolor plants (Figure S1). These results suggested the presence of a mechanically transmissible virus in PC01. To determine the identity of mechanically transmissible viral agent, symptomatic C. quinoa and PC01 plant were advanced for high throughput sequencing analysis. Aliquots of TNA from PC01 and C. quinoa were rRNA-depleted and used for cDNA library preparation with TruSeq Stranded Total RNA kit (Illumina). The raw reads were trimmed, de novo assembled, and subsequently were annotated using tBLASTx algorithm (Al Rwahnih et al. 2018). A total of 47,261,138 and 8,812,296 single-end reads were obtained from cDNA libraries of PC01 and C. quinoa, respectively. The de novo assembly generated near-complete contigs resembling AVCaV genome ) from both PC01 and C. quinoa, which were 99.8% identical at the nucleotide level. The longest contig (8,342 nucleotides, 73.5x coverage depth) obtained from PC01 was further completed using SMARTer RACE 5'/3' kit (Takara Bio). The complete genome sequence of AVCaV-PC01 is 8,364 nucleotides long (GenBank: MK170158). The full-length virus genome was compared with GenBank database using BLASTn, which the best hit corresponded to KY132099 with 98% identity. Additionally, AVCaV infection was confirmed in both PC01 selection and the symptomatic C. quinoa by RT-PCR as previously described (Marais et al. 2015). Lastly, symptomatic leaves of C. quinoa were used in leaf dip method to visualize virus particles by transmission electron microscope. As a result, flexuous rod-shaped virions were observed from leaf dips of symptomatic C. quinoa plants (Figure S2). Therefore, our results represent the first report of AVCaV in California, USA. Furthermore, mechanical transmission of an AVCaV isolate infecting flowering apricot to herbaceous hosts was confirmed. Field surveys and biological studies are underway to determine the prevalence of AVCaV in commercial orchards and assess its effect on tree performance.

Cultivated and Wild Grapevines in Tennessee Possess Overlapping but Distinct Virus Populations.

Viruses and viroids prevalent in a population of 42 wild grapevines (i.e., free-living, uncultivated grapevines; Vitis spp.) were compared with those in a population of 85 cultivated grapevines collected in Tennessee, United States by RNA sequencing analysis of pools of ribosomal RNA-depleted total RNA. The sequences of 10 viruses (grapevine fleck virus, grapevine leafroll-associated virus 2, grapevine rupestris stem pitting-associated virus, grapevine Syrah virus 1, grapevine vein-clearing virus, grapevine virus B, grapevine virus E, tobacco ringspot virus, tomato ringspot virus, and a novel nano-like virus) and two viroids (hop stunt viroid and grapevine yellow speckle viroid 1) were detected in both grapevine populations. Sequences of four viruses (grapevine associated tymo-like virus, grapevine leafroll-associated virus 3, grapevine red blotch virus, and grapevine virus H) were identified only from cultivated grapevines. High, moderate, and low numbers of sequence reads were identified only from wild grapevines for a novel caulimovirus, an enamovirus, and alfalfa mosaic virus, respectively. The presence of most virus sequences and both viroids was verified independently in the original samples by reverse-transcription PCR followed by Sanger sequencing. Comparison of viral sequences shared by both populations showed that cultivated and wild grapevines harbored distinct sequence variants, which suggests that there was limited virus movement between the two populations. Collectively, this study represents the first unbiased survey of viruses and viroids in both cultivated and wild grapevines within a defined geographic region.

First report of Rose leaf rosette-associated virus infecting rose (Rosa spp.) in California, USA.

Rose leaf rosette-associated virus (RLRaV) is a member of genus Closterovirus, family Closteroviridae. The virus was first discovered in China in 2015 from a mixed infected wild rose (Rosa multiflora Thunb.) showing small leaf rosettes on branches, dieback and severe decline symptoms (He et al. 2015). In 2013, a rose plant (cv. Roses Are Red) was introduced to Foundation Plant Services (FPS, UC-Davis) rose collection. The plant was originated from a private rose breeder collection located in California. In 2019, total nucleic acids (TNA) were isolated from leaf tissues of one asymptomatic plant (Roses Are Red plant) using MagMax Plant RNA Isolation Kit (Thermo Fisher Scientific, USA). Extracted TNA were screened by reverse-transcription quantitative PCR (RT-qPCR) for six common viruses infecting roses, including prunus necrotic ringspot virus (PNRSV), apple mosaic virus (ApMV), rose spring dwarf associated virus (RSDaV), rose yellow vein virus (RYVV), rose rosette virus (RRV), and blackberry chlorotic ringspot virus (BCRV); however, the results were negative. Therefore, the sample was subjected to high throughput sequencing (HTS). Briefly, TNA was depleted of rRNA and advanced for cDNA library preparation using TruSeq Stranded Total RNA kit (Illumina, USA). HTS was performed on Illumina NextSeq 500 platform. The raw reads were trimmed, de novo assembled, and subsequently were annotated using tBLASTx algorithm (Al Rwahnih et al. 2018). HTS generated 23.6 million 75 nucleotide (nt) single-end raw data reads. De novo assembly generated a contig (16,528 nts) resembling RLRaV reference sequence (KJ748003) with 74% identity at the nucleotide level. Putative coat protein and heat shock protein 70-like protein were identified based on >90% identity with RLRaV genes. To confirm HTS results, RT-PCR was performed using two primer sets, 1) Clo-F4916 (5'-GGTGTTCCAACGCTATCGTG-3') and Clo-R5215 (5'- TGTCCTCAAACCGCCTACAT-3') targeting nucleotide sequences of putative polyprotein 1a, and 2) Clo-F10006 (5'-GATTCCGCGGACGAATTAAT-3') and Clo-R10311 (5'-GGTAACCGAAAGGTAAAGTATTC-3') targeting nucleotide sequences of putative protein p25. The RLRaV amplicons with expected size of 300 nt were confirmed using bidirectional Sanger sequencing. The near-complete sequence of the new RLRaV isolate was deposited in GenBank under accession number MW056181. In addition, HTS analysis showed that RLRaV was in mixed infection with two mycoviruses (rose cryptic virus with 8,267 mapped reads and rose partitivirus with 7,283 mapped readss). To our knowledge, this is the first report of RLRaV affecting roses in California. Further research is needed to determine the prevalence of RLRaV in California as well as evaluation of RLRaV effect on rose performance.

Transcriptome profiles of Quercus rubra responding to increased O3 stress.

BACKGROUND: Climate plays an essential role in forest health, and climate change may increase forest productivity losses due to abiotic and biotic stress. Increased temperature leads to the increased formation of ozone (O3). Ozone is formed by the interaction of sunlight, molecular oxygen and by the reactions of chemicals commonly found in industrial and automobile emissions such as nitrogen oxides and volatile organic compounds. Although it is well known that productivity of Northern red oak (Quercus rubra) (NRO), an ecologically and economically important species in the forests of eastern North America, is reduced by exposure to O3, limited information is available on its responses to exogenous stimuli at the level of gene expression. RESULTS: RNA sequencing yielded more than 323 million high-quality raw sequence reads. De novo assembly generated 52,662 unigenes, of which more than 42,000 sequences could be annotated through homology-based searches. A total of 4140 differential expressed genes (DEGs) were detected in response to O3 stress, as compared to their respective controls. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of the O3-response DEGs revealed perturbation of several biological pathways including energy, lipid, amino acid, carbohydrate and terpenoid metabolism as well as plant-pathogen interaction. CONCLUSION: This study provides the first reference transcriptome for NRO and initial insights into the genomic responses of NRO to O3. Gene expression profiling reveals altered primary and secondary metabolism of NRO seedlings, including known defense responses such as terpenoid biosynthesis.