Discovery of Viruses and Virus-Like Pathogens in Pistachio using High-Throughput Sequencing

Molecular characterization of a divergent genetic variant of wheat Eqlid mosaic virus from a Texas wheat field

Investigations conducted during the spring 2020 season to diagnose the associated viral agent of a severe mosaic disease of wheat in a Texas Panhandle field revealed the presence of wheat Eqlid mosaic virus (WEqMV; genus Tritimovirus, family Potyviridae) in the analyzed samples. The complete genome sequences of two WEqMV isolates were determined, and each was found to be 9,634 nucleotides (nt) in length (excluding the polyA tail) and to contain 5' and 3' untranslated regions of 135 nt and 169 nt, respectively, based on rapid amplification of cDNA ends (RACE) assays. Both sequences contained an open reading frame (ORF) of 9,330 nt encoding a polyprotein of 3,109 amino acids (aa). The ORF sequences of the two isolates were 100% identical to each other, but only 74.7% identical to that of the exemplar WEqMV-Iran isolate, with 85.7% aa sequence identity in the encoded polyprotein. The Texas WEqMV isolates also diverged significantly from WEqMV-Iran in the individual proteins at the nt and aa levels. This is the first report of WEqMV in the United States and the first report of this virus outside of Iran, indicating an expansion of its geographical range.

High-Throughput Sequencing of Grapevine in Mexico Reveals a High Incidence of Viruses including a New Member of the Genus Enamovirus

This is the first viral metagenomic analysis of grapevine conducted in Mexico. During the summer of 2021, 48 plants displaying virus-like symptoms were sampled in Queretaro, an important grapevine-producing area of Mexico, and analyzed for the presence of viruses via high-throughput sequencing (HTS). The results of HTS were verified by real-time RT-PCR following a standardized testing scheme (Protocol 2010). Fourteen different viruses were identified, including grapevine asteroid mosaic-associated virus (GAMaV), grapevine Cabernet Sauvignon reovirus (GCSV), grapevine fanleaf virus (GFLV), grapevine fleck virus (GFkV), grapevine Pinot gris virus (GPGV), grapevine red globe virus (GRGV), grapevine rupestris stem pitting-associated virus (GRSPaV), grapevine rupestris vein feathering virus (GRVFV), grapevine Syrah virus 1 (GSyV-1), grapevine virus B (GVB), and grapevine leafroll-associated viruses 1, 2, 3, 4 (GLRaV1, 2, 3, 4). Additionally, divergent variants of GLRaV4 and GFkV, and a novel Enamovirus-like virus were discovered. This is the first report of GAMaV, GCSV, GLRaV4, GPGV, GRGV, GRVFV, and GSyV-1 infecting grapevines in Mexico; the impact of these pathogens on production is unknown.

Identifying the Gut Virome of Diaphorina citri from Florida Groves

Asian citrus psyllid (Diaphorina citri) transmits the bacterial pathogen Candidatus Liberibacter asiaticus (CLas), the putative causative agent of citrus Huanglongbing disease (HLB). Insect-specific viruses can act against insects as their natural enemies, and recently, several D. citri-associated viruses were discovered. The insect gut plays an important role as not only a pool for diverse microbes but also as a physical barrier to prevent the spread of pathogens such as CLas. However, there is little evidence of the presence of D. citri-associated viruses in the gut and of the interaction between them and CLas. Here, we dissected psyllid guts collected from five growing regions in Florida, and the gut virome was analyzed by high throughput sequencing. Four insect viruses, including D. citri-associated C virus (DcACV), D. citri densovirus (DcDV), D. citri reovirus (DcRV), and D. citri flavi-like virus (DcFLV), were identified, and their presence in the gut, including an additional D. citri cimodo-like virus (DcCLV), were confirmed with PCR-based assays. Microscopic analysis showed that DcFLV infection leads to morphological abnormalities in the nuclear structure in the infected psyllid gut cells. The complex and diverse composition of microbiota in the psyllid gut suggests a possible interaction and dynamics between CLas and the D. citri-associated viruses. Our study identified various D. citri-associated viruses that localized in the psyllid gut and provided more information that helps to evaluate the potential vectors for manipulating CLas in the psyllid gut.

G, Rivera Y, Mavrodieva VA, Nakhla MK, Roy A. High-throughput sequencing application in the detection and discovery of viruses associated with the regulated citrus leprosis disease complex

Citrus leprosis (CiL) is one of the destructive emerging viral diseases of citrus in the Americas. Leprosis syndrome is associated with two taxonomically distinct groups of Brevipalpus-transmitted viruses (BTVs), that consist of positive-sense Cilevirus, Higrevirus, and negative-sense Dichorhavirus. The localized CiL symptoms observed in multiple citrus species and other alternate hosts indicates that these viruses might have originated from the mites and eventually adopted citrus as a secondary host. Genetic diversity in the genomes of viruses associated with the CiL disease complex have complicated current detection and diagnostic measures that prompted the application of High-Throughput Sequencing (HTS) protocols for improved detection and diagnosis. Two cileviruses are known to infect citrus, and among them only citrus leprosis virus C2 (CiLV-C2) hibiscus strain (CiLV-C2H) has been reported in hibiscus and passion fruit in the US. Based on our current CiL disease complex hypothesis, there is a high probability that CiL disease is associated with more viruses/strains that have not yet been identified but exist in nature. To protect the citrus industry, a Ribo-Zero HTS protocol was utilized for detection of cileviruses infecting three different hosts: Citrus spp., Swinglea glutinosa, and Hibiscus rosa-sinensis. Real-time RT-PCR assays were used to identify plants infected with CiLV-C2 or CiLV-C2H or both in mixed infection in all the above-mentioned plant genera. These results were further confirmed by bioinformatic analysis using HTS generated data. In this study, we utilized HTS assay in confirmatory diagnostics to screen BTVs infecting Dieffenbachia sp. (family: Araceae), Passiflora edulis (Passifloraceae), and Smilax auriculata (Smilacaceae). Through the implementation of HTS and downstream data analysis, we detected not only the known cileviruses in the studied hosts but also discovered a new strain of CiLV-C2 in hibiscus from Colombia. Phylogenetically, the new hibiscus strain is more closely related to CiLV-C2 than the known hibiscus strain, CiLV-C2H. We propose this strain to be named as CiLV-C2 hibiscus strain 2 (CiLV-C2H2). The findings from the study are critical for citrus growers, industry, regulators, and researchers. The possible movement of CiLV-C2H2 from hibiscus to citrus by the Brevipalpus spp. warrants further investigation.

Molecular characterization of divergent isolates of Citrus bent leaf viroid (CBLVd) from citrus cultivars of Punjab, Pakistan

Citrus viroid infection is emerging as a serious threat because of its efficient systemic movement within the host plant and its quick spread due to contaminated pruning tools. A survey was conducted to investigate the primary distribution and molecular characterization of Citrus bent leaf viroid (CBLVd) and its variants in different citrus cultivars. A total of 154 symptomatic citrus samples were collected and detected by RT‒PCR with newly designed specific primers with the incidence of 36.33%. During biological indexing study on Etrog citron, expressions of reduced leaf size, yellowing with a light green pattern, and bending were observed. Amplified products were sequenced and analyzed using a nucleotide BLAST search, which showed 98% homology with other CBLVd isolates. The results of the phylogenetic tree analysis showed the presence of two main groups (A and B), with the predominant variants of CBLVd, i.e., CVd-I-LSS (Citrus viroid Low Sequence Similarity) sequences, clustering in subgroup A1 along with newly detected CVd-I-LSS from Palestinian sweet lime (Citrus limettioides), which has been identified as a new host of CVd-I-LSS in Pakistan. Further analysis of the sequences in subgroup A1 showed that the variant of CVd-I-LSS infecting citrus cultivars had a close relationship with isolates reported from China, Japan, and Iran, which may have resulted from the exchange of planting material. This study also unveiled the variability in nucleotide sequences of CBLVd, which made it unable to be detected by old primers. The results of this study indicate that the widespread presence of divergent variants of CBLVd is a major concern for the citrus industry in Pakistan and other countries where virulent isolates of CBLVd are prevalent. These findings suggest the need for future research on effective management and quarantine measures to stop the spread of CBLVd.

Genetic Diversity of Viral Populations Associated with Ananas Germplasm and Improvement of Virus Diagnostic Protocols

Pineapple (Ananas comosus L. [Merr.]) accessions from the U.S. Tropical Plant Genetic Resources and Disease Research (TPGRDR) in Hilo, Hawaii were subjected to RNA-sequencing to study the occurrence of viral populations associated with this vegetatively propagated crop. Analysis of high-throughput sequencing data obtained from 24 germplasm accessions and public domain transcriptome shotgun assembly (TSA) data identified two novel sadwaviruses, putatively named "pineapple secovirus C" (PSV-C) and "pineapple secovirus D" (PSV-D). They shared low amino acid sequence identity (from 34.8 to 41.3%) compared with their homologs in the Pro-pol region of the previously reported PSV-A and PSV-B. The complete genome (7485 bp) corresponding to a previously reported partial sequence of the badnavirus, pineapple bacilliform ER virus (PBERV), was retrieved from one of the datasets. Overall, we discovered a total of 69 viral sequences representing ten members within the Ampelovirus, Sadwavirus, and Badnavirus genera. Genetic diversity and recombination events were found in members of the pineapple mealybug wilt-associated virus (PMWaV) complex as well as PSVs. PMWaV-1, -3, and -6 presented recombination events across the quintuple gene block, while no recombination events were found for PMWaV-2. High recombination frequency of the RNA1 and RNA2 molecules from PSV-A and PSV-B were congruent with the diversity found by phylogenetic analyses. Here, we also report the development and improvement of RT-PCR diagnostic protocols for the specific identification and detection of viruses infecting pineapple based on the diverse viral populations characterized in this study. Given the high occurrence of recombination events, diversity, and discovery of viruses found in Ananas germplasm, the reported and validated RT-PCR assays represent an important advance for surveillance of viral infections of pineapple.

Survey for Virus Diversity in Common Bean (Phaseolus vulgaris) Fields and the Detection of a Novel Strain of Cowpea polerovirus 1 in Zambia

The production of common bean (Phaseolus vulgaris L.) is adversely affected by virus-like diseases globally, but little is known about the occurrence, distribution, and diversity of common bean-infecting viruses in Zambia. Consequently, field surveys were conducted during the 2018 season in 128 fields across six provinces of Zambia and 640 common bean leaf tissue samples were collected with (n = 585) or without (n = 55) symptoms. The prevalence of symptomatic fields was 100%, but incidence of symptomatic plants ranged from 32 to 67.5%. Metagenomic analyses of nine composite samples and a single plant sample of interest revealed the occurrence of isolates of Bean common mosaic necrosis virus, Bean common mosaic virus, Cowpea aphid-borne mosaic virus, Peanut mottle virus, Southern bean mosaic virus (SBMV), Cucumber mosaic virus, Phaseolus vulgaris alphaendornavirus 1 (PvEV-1), PvEV-2, Ethiopian tobacco bushy top virus (ETBTV), and a novel strain of Cowpea polerovirus 1 (CPPV1-Pv) of 5,902 nt in length. While CPPV1-Pv was consistently detected in mixed infection with ETBTV and its satellite RNA molecule, based on results of mechanical transmission assays it does not appear to be involved in disease etiology, suggesting that its role may be limited to being a helper virus for the umbravirus. Screening of the survey samples by real-time PCR for the viruses detected by high-throughput sequencing revealed the prevalence of single (65.2% or 417/640) over mixed (1.9% or 12/640) infections in the samples. SBMV was the most frequently detected virus, occurring in ∼29.4% (188/640) of the samples and at a prevalence rate of 58.6% (75/128) across fields. The results showed that diverse virus species are present in Zambian common bean fields and the information will be useful for the management of common bean viral diseases.

Preserving plant samples from remote locations for detection of RNA and DNA viruses

Viral diseases in plants have a significant impact on agricultural productivity. Effective detection is needed to facilitate accurate diagnosis and characterization of virus infections essential for crop protection and disease management. For sensitive polymerase chain reaction (PCR)-based methods, it is important to preserve the integrity of nucleic acids in plant tissue samples. This is especially critical when samples are collected from isolated areas, regions distant from a laboratory, or in developing countries that lack appropriate facilities or equipment for diagnostic analyses. RNAlater^® provides effective, reliable sample storage by stabilizing both RNA and DNA in plant tissue samples. Our work indicated that total RNA or DNA extracted from virus-infected leaf samples preserved in RNAlater^® was suitable for reverse transcription polymerase chain reaction (RT-PCR), PCR, Sanger sequencing, high-throughput sequencing (HTS), and enzyme-linked immunosorbent assay (ELISA)-based diagnostic analyses. We demonstrated the effectiveness of this technology using leaf tissue samples from plants with virus symptoms grown in farmers’ fields in Bangladesh. The results revealed that RNAlater^® technology was effective for detection and characterization of viruses from samples collected from remote areas and stored for extended periods. Adoption of this technology by developing countries with limited laboratory facilities could greatly increase their capacity to detect and diagnose viral infections in crop plants using modern analytical techniques.

Analysis of Citrus Tristeza Virus Incidences within Asian Citrus Psyllid (Diaphorina citri) Populations in Florida via High-Throughput Sequencing

The destructive citrus disease, Huanglongbing (HLB) or citrus greening, continues to devastate Florida’s citrus industry. A hemipteran insect, the Asian citrus psyllid (ACP), disperses Candidatus Liberibacter asiaticus, one of the putative bacterial pathogens of HLB. This study builds upon ongoing research utilizing high-throughput sequencing to analyze the virome of ACP populations collected from citrus groves throughout Florida. Following the widespread detection of sequences aligning to the genome of citrus tristeza virus (CTV) across consecutive years in the Florida ACP virome, we continued to detect a pervasive amount of CTV in Florida ACPs during subsequent years. Simultaneously, we also detected mixed infections of CTV strains in pooled ACPs from different Florida regions. Predating the HLB epidemic, CTV has been present in Florida for many years and our results confirm its widespread and diverse persistence in Florida citrus groves through a unique lens, the ACP. CTV presence in the ACP likely results from feeding on CTV-infected citrus trees in Florida citrus groves, which may help to understand an overlapping presence of CTV and HLB, both endemic citrus pathosystems in the state, and their role in future integrated pest management strategies.

Diaphorina citri flavi-like virus localization, transmission, and association with Candidatus Liberibacter asiaticus in its psyllid host

Huanglongbing is caused by Candidatus Liberibacter asiaticus (CLas) and transmitted by Diaphorina citri. D. citri harbors various insect-specific viruses, including the Diaphorina citri flavi-like virus (DcFLV). The distribution and biological role of DcFLV in its host and the relationship with CLas are unknown. DcFLV was found in various organs of D. citri, including the midgut and salivary glands, where it co-localized with CLas. CLas-infected nymphs had the highest DcFLV titers compared to the infected adults and CLas-free adults and nymphs. DcFLV was vertically transmitted to offspring from female D. citri and was temporarily detected in Citrus macrophylla and grapefruit leaves from greenhouse and field. The incidences of DcFLV and CLas were positively correlated in field-collected D. citri samples, suggesting that DcFLV might be associated with CLas in the vector. These results provide new insights on the interactions between DcFLV, the D. citri, and CLas.

Olea europaea geminivirus is present in a germplasm repository and in California and Texas olive (Olea europaea L.) groves

Olea europaea geminivirus (OEGV) from olive accessions in Italy was characterized recently. OEGV was also detected during routine high-throughput sequencing screening of olive (cv. Leccino) material, and its complete bipartite genome segments were sequenced and shown to be 100% identical to those of the isolate from Italy. Using two pairs of newly designed primers targeting the AV1 and BV1 genes, OEGV was detected in randomly sampled olive trees from the U.S. Department of Agriculture National Clonal Germplasm Repository (USDA-NCGR) (21.4% or 6/28), commercial and residential settings in California (47.6% or 10/21), and an orchard in Texas (60% or 30/50). The cuttings for the USDA-NCGR-positive trees originated from the former Serbia and Montenegro, Spain, Italy, and Greece. Comparative analysis of the directly sequenced gene fragments from randomly selected samples showed that OEGV isolates from the different sources were 100% identical to each other. The results indicate that OEGV spread was likely facilitated by inadvertent movement of contaminated olive germplasm.

Partial Genome Sequence of a Novel Reo-Like Virus Detected in Asian Citrus Psyllid (Diaphorina citri) Populations from Florida Citrus Groves

This report describes the partial (nearly complete) genome sequence of a novel reo-like virus tentatively named Diaphorina citri Cimodo-like virus. This putative virus has 10 double-stranded RNA segments and was detected in Asian citrus psyllid (Diaphorina citri) populations collected from Florida commercial citrus groves.

The Sisal Virome: Uncovering the Viral Diversity of Agave Varieties Reveals New and Organ-Specific Viruses

Sisal is a common name for different plant varieties in the genus Agave (especially Agave sisalana) used for high-quality natural leaf fiber extraction. Despite the economic value of these plants, we still lack information about the diversity of viruses (virome) in non-tequilana species from the genus Agave. In this work, by associating RNA and DNA deep sequencing we were able to identify 25 putative viral species infecting A. sisalana, A. fourcroydes, and Agave hybrid 11648, including one strain of Cowpea Mild Mottle Virus (CPMMV) and 24 elements likely representing new viruses. Phylogenetic analysis indicated they belong to at least six viral families: Alphaflexiviridae, Betaflexiviridae, Botourmiaviridae, Closteroviridae, Partitiviridae, Virgaviridae, and three distinct unclassified groups. We observed higher viral taxa richness in roots when compared to leaves and stems. Furthermore, leaves and stems are very similar diversity-wise, with a lower number of taxa and dominance of a single viral species. Finally, approximately 50% of the identified viruses were found in all Agave organs investigated, which suggests that they likely produce a systemic infection. This is the first metatranscriptomics study focused on viral identification in species from the genus Agave. Despite having analyzed symptomless individuals, we identified several viruses supposedly infecting Agave species, including organ-specific and systemic species. Surprisingly, some of these putative viruses are probably infecting microorganisms composing the plant microbiota. Altogether, our results reinforce the importance of unbiased strategies for the identification and monitoring of viruses in plant species, including those with asymptomatic phenotypes.

Sequencing a Strawberry Germplasm Collection Reveals New Viral Genetic Diversity and the Basis for New RT-qPCR Assays

Viruses are considered of major importance in strawberry (Fragaria × ananassa Duchesne) production given their negative impact on plant vigor and growth. Strawberry accessions from the National Clonal Germplasm Repository were screened for viruses using high throughput sequencing (HTS). Analyses of sequence information from 45 plants identified multiple variants of 14 known viruses, comprising strawberry mottle virus (SMoV), beet pseudo yellows virus (BPYV), strawberry pallidosis-associated virus (SPaV), tomato ringspot virus (ToRSV), strawberry mild yellow edge virus (SMYEV), strawberry vein banding virus (SVBV), strawberry crinkle virus (SCV), strawberry polerovirus 1 (SPV-1), apple mosaic virus (ApMV), strawberry chlorotic fleck virus (SCFaV), strawberry crinivirus 4 (SCrV-4), strawberry crinivirus 3 (SCrV-3), Fragaria chiloensis latent virus (FClLV) and Fragaria chiloensis cryptic virus (FCCV). Genetic diversity of sequenced virus isolates was investigated via sequence homology analysis, and partial-genome sequences were deposited into GenBank. To confirm the HTS results and expand the detection of strawberry viruses, new reverse transcription quantitative PCR (RT-qPCR) assays were designed for the above-listed viruses. Further in silico and in vitro validation of the new diagnostic assays indicated high efficiency and reliability. Thus, the occurrence of different viruses, including divergent variants, among the strawberries was verified. This is the first viral metagenomic survey in strawberry, additionally, this study describes the design and validation of multiple RT-qPCR assays for strawberry viruses, which represent important detection tools for clean plant programs.

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.

First report of squash vein yellowing virus naturally infecting butternut squash (Cucurbita moschata) in Texas

Virus diseases are major constraints to the production of cucurbits in the Texas Lower Rio Grande Valley. In September 2020, a ~8.1 ha butternut squash (Cucurbita moschata) field in Hidalgo County, Texas, was observed with virus-like symptoms of vein yellowing, leaf curl, mosaic, and foliar chlorosis. The proportion of plants with virus-like symptoms in this field was estimated at 30% and seven samples (symptomatic = 5; non-symptomatic = 2) were collected randomly for virus diagnosis. Initially, equimolar mixtures of total nucleic acid extracts (Dellaporta et. al. 1983) from two symptomatic samples from this field and extracts from 12 additional symptomatic samples from six other fields across south and central Texas was used to generate one composite sample for diagnosis by high throughput sequencing (HTS). The TruSeq Stranded Total RNA with Ribo-Zero Plant Kit (Illumina) was used to construct cDNA library from the composite sample, which was then sequenced on the Illumina NextSeq 500 platform. More than 26 million single-end HTS reads (75 nt each) were obtained and their bioinformatic analyses (Al Rwahnih et al. 2018) revealed several virus-like contigs belonging to different species (data not shown). Among them, 6 contigs that ranged in length from 429 to 3,834 nt shared 96 to 100% identities with isolates of squash vein yellowing virus (SqVYV), genus Ipomovirus, family Potyviridae. To confirm the HTS results, total nucleic acid extracts from the cucurbit samples from all seven fields (n = 46) were used for cDNA synthesis with random hexamers and the PrimeScript 1st strand cDNA Synthesis Kit (Takara Bio). A 1-μL aliquot of cDNA was used in 12.5-μL PCR reaction volumes with PrimeSTAR GXL DNA Polymerase (Takara Bio) and two pairs of SqVYV-specific primers designed based on the HTS derived contigs. The primer pairs SqYVV-v4762: 5'-CTGGATTCTGCTGGAAGATCA & SqYVV-c5512: 5'-CCACCATTAAGGCCATCAAAC and SqYVV-v8478: 5'-TTTCTGGGCAAACAAACATGG & SqYVV-c9715: 5'-TTCAGCGACGTCAAGTGAG targeted ~0.75 kb and ~1.2 kb fragments of the cylindrical inclusion (CI) and the complete coat protein (CP) gene sequences of SqVYV, respectively. The expected DNA band sizes were obtained only from the five symptomatic butternut squash samples from the Hidalgo Co. field. Two amplicons per primer pair from two samples were cloned into pJET1.2/Blunt vector (Life Technologies) and bidirectionally Sanger sequenced, generating 753 nt partial CI specific sequences (MW584341-342) and 1,238 nt that encompassed the complete CP (MW584343-344) of SqVYV. In pairwise comparisons, the partial CI sequences shared 100% nt/aa identity with each other and 98-99% nt/aa identity with corresponding sequences of SqVYV isolate IL (KT721735). The CP cistron of TX isolates shared 100% nt/aa identity with each other and 90-98% nt (97-100% aa) identities with corresponding sequences of several SqVYV isolates in GenBank, with isolates IL (KT721735) and Florida (EU259611) being at the high and low spectrum of nt/aa identity values, respectively. This is the first report of SqVYV in Texas, naturally occurring in butternut squash. SqVYV was first discovered in Florida (Adkins et al. 2007) and subsequently reported from few other states in the U.S. (Adkins et al. 2013; Egel and Adkins 2007; Batuman et al. 2015), Puerto Rico (Acevedo et al. 2013), and locations around the world. The finding shows an expansion of the geographical range of SqVYV and adds to the repertoire of cucurbit-infecting viruses in Texas. Further studies are needed to determine the prevalence of SqVYV in Texas cucurbit fields and an assessment of their genetic diversity.

First report of Cucurbit chlorotic yellows virus infecting cantaloupe (Cucumis melo L.) in Texas

Texas is a major producer of cucurbits such as cantaloupe (Cucumis melo L.), but outbreaks of virus-like diseases often adversely affect yields. Little is known about the identity of the causal or associated viruses. During studies conducted in fall 2020 to explore the virome of cucurbit fields in Texas, a commercial cantaloupe field (~4.1 ha) in Cameron County was observed with virus-like symptoms of interveinal chlorotic mottle and foliar chlorosis and disease incidence was estimated at 100%. Virus-like symptoms including mosaic and leaf curl were also observed in six additional fields across five south and central Texas counties of Atascosa, Hidalgo, Fort Bend, Frio, and Wharton. Forty-six plants, which included 32 symptomatic and 14 non-symptomatic, were sampled from these fields for virus diagnosis and each sample was subjected to total nucleic acid extraction according to Dellaporta et. al. (1983). Initially, equal amounts of nucleic acids from 14 symptomatic plants (two/field) were pooled into one composite sample for preliminary diagnosis by high throughput sequencing (HTS). The cDNA library obtained from the composite sample with a TruSeq Stranded Total RNA with Ribo-Zero Plant Kit (Illumina) was sequenced on the Illumina NextSeq 500 platform, generating ~26.3 M single-end HTS reads (75 nucleotides [nt] each). Analyses of the reads according to Al Rwahnih et al. (2018) revealed several virus-like contigs; among them 23 contigs (206 to 741 nt) shared 98 to 100% nt identities to isolates of cucurbit chlorotic yellows virus (CCYV), genus Crinivirus, family Closteroviridae. Three pairs of CCYV-specific primers were designed from the HTS contigs with primers CCYV-v1330: 5'-AGTCCCTTACCCTGAGATGAA/CCYV-c2369: 5'-CGGAGCATTCGACAACTGAATA targeting ~1 kb fragment of the ORF1a (RNA1), primers CCYV-v4881: 5'-ATAAGGCGGCGACCTAATC/CCYV-c5736: 5'-GATCACTTGACCATCTCCTTCT targeting a ~0.9 kb fragment encompassing the coat protein (CP) cistron of CCYV (RNA2), and primers CCYV-v6362: 5'-CACCTCTTCCAGAACCAGTTAAA/CCYV-c7423: 5'-TGGGAACAACTTATTTCTCCTAGC targeting ~1 kb spanning partial minor coat protein (CPm) and p26 sequences (RNA2). Total nucleic acid extracts of each of the 46 samples from the seven fields were tested by two-step reverse transcription polymerase chain reaction using all three CCYV-specific primer pairs and they yielded amplicons of expected sizes from all five symptomatic cantaloupe samples from the Cameron County field and one additional symptomatic butternut squash sample from a field in Hidalgo County. The DNA bands from three randomly chosen cantaloupe samples were cloned and sequenced as previously described (Oke et al. 2020). In pairwise comparisons, the obtained 1,040 nt ORF1a (MW584332-334), 753 nt complete CP (MW584335-337), and 1,062 nt CPm/p26 (MW584338-340) gene specific sequences from the three samples shared 100% nt identity with each other, and 99-100% nt identities with corresponding RNA1 (AB523788) and RNA2 (AB523788) sequences of the exemplar isolate of CCYV. This is the first report of CCYV in Texas, thus expanding the current geographical range of the virus in the U.S. that includes California (Wintermantel et al. 2019) and Georgia (Kavalappara et al. 2021). The abundance of whiteflies of the Bemisia tabaci species complex in south Texas and other major U.S. cucurbit production areas presents additional challenges to virus disease management.

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.

First report of watermelon crinkle leaf-associated virus 1 (WCLaV-1) and WCLaV-2 infecting watermelon (Citrullus lanatus) in the United States

Watermelon (Citrullus lanatus) and other cucurbits are cultivated globally, and Texas ranks among its top 5 producers in the U.S. In July 2020, plants with virus-like disease symptoms consisting of mild leaf crinkling and yellow mosaic patterns were observed in a 174-ha watermelon field in Burleson Co., TX; disease incidence was visually estimated at 5%. Total nucleic acids were extracted from leaf tissues of 5 randomly sampled plants (Dellaporta 1983) and their equimolar amounts were made into a composite sample that was used for cDNA library construction with TruSeq Stranded Total RNA with Ribo-Zero Plant Kit (Illumina). The cDNA library was sequenced on the Illumina NextSeq 500 platform, generating ~37M single-end reads (each 75 nt), which were analyzed as per Al Rwahnih et al. (2018). Of these, 58,200 and 27,500 reads mapped to the genomes of watermelon crinkle leaf-associated virus 1 (WCLaV-1) and WCLaV-2 (Xin et al. 2017), respectively, along with 4 other virus-specific reads (data not shown). The near complete RNA1-RNA3 segments of WCLaV-1 (354-652X) and WCLaV-2 (144-258X) were generated from the mapped reads and they shared ≥96% nt identities with published RNA segments of both viruses. The results were verified by RT-PCR using newly designed primers WCLaV-1vRP: 5'-GGTGAGTTAGTGTGTCTGAAGG/WCLaV-1cRP: 5'-GAGGTTGCCTGAGGTGATAAG to target 881 bp of the RNA1-encoded RNA-dependent RNA polymerase (RdRP), WCLaV-1vMP: 5'-GAAGGTTTGCTCCCTTGAAATG/WCLaV-1cMP: 5'-GACTGTGGCTGAAGAGTCTATG target 538 bp of the RNA2-encoded movement protein (MP), and WCLaV-1vNP: 5'-CGAATAGACTCTGGAGGGTAGA/WCLaV-1cMP: 5'-GAAAGCAAGAAAGCTGGCTAAA target 786 bp of the RNA3-encoded nucleoprotein (NP). Similarly, the WWCLaV-2-specific primers WCLaV-2vRP: 5'-GTCTCACATTCCTGCACTAACT/WCLaV-2cRP: 5'-ATCGGTCCTGGGTTATTTGTATC target 968 bp of the RdRP, WCLaV-2vMP: 5'-GACTTCAGAACCTCAACATCCA/WCLaV-2cMP: 5'-CAAGGGAGAGTGCTGACAAA target 562 bp of the MP, and WCLaV-2vNP: 5'-ATTCCCAGTGAGAGCAACAA/WCLaV-2cMP: 5'-GAGGTGGAGGTAGGAAAGAAAG target 449 bp of the NP. Fresh cDNA synthesized from the 5 samples with PrimeScript First Strand cDNA synthesis kit (Takara Bio) were tested by PCR with all 6 primer pairs using the PrimeSTAR GXL DNA Polymerase kit (Takara Bio). Three of the 5 samples were positive for both viruses and one sample was positive for each virus. The obtained products from 4 samples were cloned individually into pJET1.2/Blunt vector (Thermo Scientific, USA), followed by bidirectional Sanger-sequencing of the plasmids with the GenElute Five-Minute Plasmid Miniprep kit (Sigma-Aldrich). In pairwise comparisons, the partial RNA1-RNA3 sequences of WCLaV-1 (GenBank accession nos. MW559074-82) shared 100% nt/aa identities with each other and with corresponding sequences of WCLaV-1 isolate KF-1 from China (KY781184-86). The partial RNA1-RNA3 sequences of WCLaV-2 (MW559083-91) shared 97-100% nt/96-100% aa identities with each other and with corresponding sequences of WCLaV-2 isolate KF-15 from China (KY781187-89). This is the first report of WCLaV-1 and WCLaV-2 in Texas and the first record of both viruses in the U.S. and elsewhere outside of China. Both negative-sense, single-stranded RNA viruses represent a novel taxon in the family Phenuiviridae (order Bunyavirales) (Xin et al. 2017). While aspects of the biology of both viruses are yet to be elucidated, our results expand their geographical range. The detection primers developed here will be useful for screening cucurbits germplasm to avert their spread.

Reassessing species demarcation criteria in viroid taxonomy by pairwise identity matrices

With a small, circular and non-protein coding RNA genome, viroids are the smallest infectious agents. They invade plants, which in turn may develop symptoms. Since their discovery about 50 years ago, more than thirty viroids have been reported and classified using as species demarcation less than 90 per cent sequence identity on the overall genome and evidence of biological divergence with respect to the closest related viroids. In the last few years, new viroids have been identified that infect latently their (frequently) woody hosts and have a narrow experimental hosts range, complicating and slowing down studies on their biology. As a consequence, several viroids are still waiting for classification. Moreover, the number of new viroids is expected to increase in the next years due to the use of high-throughput sequencing technologies with diagnostics purposes. Therefore, establishment of reliable species demarcation criteria mainly based on molecular features of viroids is needed. Here, viroid classification is reassessed and a scheme based on pairwise sequence identity matrices is developed. After identifying a threshold pairwise identity score (PWIS) for each viroid genus, to be used as a species demarcation criterion, we show that most of those yet unclassified viroids can be assigned to a known or to a new species, thus limiting the need for additional biological evidence to only a few more complex situations. The advantages of this PWIS-based method are that the proposed identity thresholds for species demarcations are not arbitrarily established and evidence for biological divergence is not mandatory. Importantly, the current classification is not essentially modified. A protocol for a tentative fast classification of new viroids according to the proposed approach is also provided.

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.

Next-Generation Sequencing and the CRISPR-Cas Nexus: A Molecular Plant Virology Perspective

In recent years, next-generation sequencing (NGS) and contemporary Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-CRISPR-associated (Cas) technologies have revolutionized the life sciences and the field of plant virology. Both these technologies offer an unparalleled platform for sequencing and deciphering viral metagenomes promptly. Over the past two decades, NGS technologies have improved enormously and have impacted plant virology. NGS has enabled the detection of plant viruses that were previously undetectable by conventional approaches, such as quarantine and archeological plant samples, and has helped to track the evolutionary footprints of viral pathogens. The CRISPR-Cas-based genome editing (GE) and detection techniques have enabled the development of effective approaches to virus resistance. Different versions of CRISPR-Cas have been employed to successfully confer resistance against diverse plant viruses by directly targeting the virus genome or indirectly editing certain host susceptibility factors. Applications of CRISPR-Cas systems include targeted insertion and/or deletion, site-directed mutagenesis, induction/expression/repression of the gene(s), epigenome re-modeling, and SNPs detection. The CRISPR-Cas toolbox has been equipped with precision GE tools to engineer the target genome with and without double-stranded (ds) breaks or donor templates. This technique has also enabled the generation of transgene-free genetically engineered plants, DNA repair, base substitution, prime editing, detection of small molecules, and biosensing in plant virology. This review discusses the utilities, advantages, applications, bottlenecks of NGS, and CRISPR-Cas in plant virology.

Comprehensive Real-Time RT-PCR Assays for the Detection of Fifteen Viruses Infecting Prunus spp

Viruses can cause economic losses in fruit trees, including Prunus spp., by reducing yield and marketable fruit. Given the genetic diversity of viruses, reliable diagnostic methods relying on PCR are critical in determining viral infection in fruit trees. This study evaluated the broad-range detection capacity of currently available real-time RT-PCR assays for Prunus-infecting viruses and developed new assays when current tests were inadequate or absent. Available assays for 15 different viruses were exhaustively evaluated in silico to determine their capacity to detect virus isolates deposited in GenBank. During this evaluation, several isolates deposited since the assay was designed exhibited nucleotide mismatches in relation to the existing assay’s primer sequences. In cases where updating an existing assay was impractical, we performed a redesign with the dual goals of assay compactness and comprehensive inclusion of genetic diversity. The efficiency of each developed assay was determined by a standard curve. To validate the assay designs, we tested them against a comprehensive set of 87 positive and negative Prunus samples independently analyzed by high throughput sequencing. As a result, all the real-time RT-PCR assays described herein successfully detected the different viruses and their corresponding isolates. To further validate the new and updated assays a Prunus germplasm collection was surveyed. The sensitive and reliable detection methods described here will be used for the large-scale pathogen testing required to maintain the highest quality nursery stock.

The Detection and Surveillance of Asian Citrus Psyllid (Diaphorina citri)-Associated Viruses in Florida Citrus Groves

The plant pathogenic bacterium Candidatus Liberibacter asiaticus (CLas), the causal agent of the citrus disease Huanglongbing (HLB), and its insect vector, the Asian citrus psyllid (ACP; Diaphorina citri), have been devastating the Florida citrus industry. To restore the competitive production presence of Florida in the worldwide citrus market, effective and sustainable control of HLB and the ACP needs to be identified. As alternatives for resistance-inducing insecticides, viruses are currently being considered for biological control of the ACP. To identify possible biological control candidates, we conducted one of the most comprehensive surveys of natural ACP populations in major citrus production regions spanning 21 counties in Florida. By optimizing PCRs and RT-PCRs, we were able to successfully detect and monitor the prevalence of five previously identified ACP-associated RNA and DNA viruses throughout Florida citrus groves, which include: Diaphorina citri-associated C virus (DcACV), Diaphorina citri flavi-like virus (DcFLV), Diaphorina citri densovirus (DcDNV), Diaphorina citri reovirus (DcRV), and Diaphorina citri picorna-like virus (DcPLV). Adult and nymph ACP populations from 21 of Florida's major citrus-producing counties were collected each month during approximately 18 consecutive months. RNA extracts used for these viral screens were also regionally combined and subjected to High Throughput Sequencing (HTS) to reveal a more comprehensive picture of known and unknown viruses in Florida ACP populations. We discovered that DcACV was the most prevalent ACP-associated virus throughout nymph and adult ACP populations in Florida, detected in more than 60% of all samples tested, followed by DcPLV and DcFLV. HTS allowed us to identify a novel ACP-associated reo-like virus and a picorna-like virus. The putative reo-like virus, tentatively named Diaphorina citri cimodo-like virus, was later surveyed and detected back in seasonal adult and nymph ACP samples collected in Florida during this study. HTS generated data also revealed that the most abundant virus in Florida ACP populations was Citrus tristeza virus (CTV), which is not an ACP-associated virus, suggesting persistent presence of CTV infection in citrus throughout Florida groves. Collectively, information obtained from our study may be able to help guide the direction of biotechnological pest control efforts involving a number of viruses that were detected for the first time in Florida ACP populations, including two newly identified ACP-associated viruses.

Papaya Ringspot Virus Isolates From Papaya in Bangladesh: Detection, Characterization, and Distribution

Papaya ringspot virus (PRSV) is the major constraint to papaya (Carica papaya) production in Bangladesh. Disease symptoms occurred in 90 to 100% of the plants surveyed. Full-length genomes of PRSV strains from severely infected papaya plants were determined using the Illumina NextSeq 500 platform, followed by Sanger DNA sequencing of viral genomes obtained by reverse-transcription PCR(RT-PCR). The genome sequences of two distinct PRSV strains, PRSV BD-1 (10,300 bp) and PRSV BD-2 (10,325 bp) were 74 and 83% identical to each other, respectively, at the nucleotide and amino acid levels. PRSV BD-1 and PRSV BD-2 were 74 to 75% and 79 to 88% identical, respectively, to other full-length PRSV sequences at the nucleotide level. Based on phylogenetic analysis, PRSV BD-2 was most closely related to PRSV-Meghalaya (MF356497) from papaya in India. PRSV BD-1 formed a branch distinct from the other PRSV sequences based on nucleotide and amino acid sequence comparisons. Comparisons of the genome sequences of these two strains with other sequenced PRSV genomes indicated two putative recombination events in PRSV BD-2. One recombinant event contained a 2,766-nucleotide fragment highly identical to PRSV-Meghalaya (MF356497). The other recombinant event contained a 5,105-nucleotide fragment highly identical to PRSV-China (KY933061). The occurrence rates of PRSV BD-1 and PRSV BD-2 in the sampled areas of Bangladesh were approximately 19 and 69%, respectively. Plants infected with both strains (11%) exhibited more severe symptoms than plants infected with either strain alone. The full-length genome sequences of these new PRSV strains and their distribution provide important information regarding the dynamics of papaya ringspot virus infections in papaya in Bangladesh.

Two Novel Negative-Sense RNA Viruses Infecting Grapevine Are Members of a Newly Proposed Genus within the Family Phenuiviridae

Two novel negative-stranded (ns)RNA viruses were identified by high throughput sequencing in grapevine. The genomes of both viruses, named grapevine Muscat rose virus (GMRV) and grapevine Garan dmak virus (GGDV), comprise three segments with each containing a unique gene. Based on sequence identity and presence of typical domains/motifs, the proteins encoded by the two viruses were predicted to be: RNA-dependent RNA polymerase (RdRp), nucleocapsid protein (NP), and putative movement protein (MP). These proteins showed the highest identities with orthologs in the recently discovered apple rubbery wood viruses 1 and 2, members of a tentative genus (Rubodvirus) within the family Phenuiviridae. The three segments of GMRV and GGDV share almost identical sequences at their 5' and 3' termini, which are also complementary to each other and may form a panhandle structure. Phylogenetics based on RdRp, NP and MP placed GMRV and GGDV in the same cluster with rubodviruses. Grapevine collections were screened for the presence of both novel viruses via RT-PCR, identifying infected plants. GMRV and GGDV were successfully graft-transmitted, thus, they are the first nsRNA viruses identified and transmitted in grapevine. Lastly, different evolutionary scenarios of nsRNA viruses are discussed.

Grapevine virus M, a novel vitivirus discovered in the American hybrid bunch grape cultivar Blanc du Bois in Texas

A novel ssRNA (+) virus with molecular properties typical of members of the genus Vitivirus (family Betaflexiviridae; subfamily Trivirinae) was discovered by high-throughput sequencing in samples of the American hybrid bunch grape cultivar Blanc du Bois in Texas. The results were independently confirmed by Sanger sequencing of the virus isolate, whose genome length is 7,387 nt, excluding the polyA tail. The genome sequence contains five ORFs that are homologous and phylogenetically related to ORFs of grapevine-infecting vitiviruses. The name "grapevine virus M" is proposed for this new virus, whose sequence divergence exceeds the current ICTV species demarcation threshold for the genus Vitivirus.