Discovery and characterization of two highly divergent variants of a novel potyvirus species infecting Madagascar periwinkle (Catharanthus roseus L.).

During summer 2022, a cluster of Madagascar periwinkle plants with white and mauve flowers were observed with foliar mild yellow mosaic symptoms on a private property in Harlingen, Cameron County, Texas. The symptoms were reproduced on mechanically inoculated periwinkle and Nicotiana benthamiana plants. Virions of 776 to 849 nm in length and 11.7 to 14.8 nm in width were observed in transmission electron microscopy of leaf dip preparations made from symptomatic periwinkle leaves. Highthroughput sequencing (HTS) analysis of total RNA extracts from symptomatic leaves revealed the occurrence of two highly divergent variants of a novel Potyvirus species as the only virus-like sequences present in the sample. The complete genomes of both variants were independently amplified via RT-PCR, cloned, and Sanger sequenced. The 5' and 3' of the genomes were acquired using RACE methodology. The assembled virus genomes were 9,936 and 9,944 nucleotides (nt) long and they shared 99.9-100% identities with the respective HTS-derived genomes. Each genome encoded hypothetical polyprotein of 3,171 amino acids (aa) (362.6 kDa) and 3,173 aa (362.7 kDa), respectively, and they shared 77.3%/84.4% nt/aa polyproteins identities, indicating that they represent highly divergent variants of the same Potyvirus species. Both genomes also shared below species threshold polyprotein identity levels with the most closely phylogenetically related known potyviruses thus indicating that they belong to a novel species. The name periwinkle mild yellow mosaic virus (PwMYMV) is given to the potyvirus with complete genomes of 9,936 nt for variant 1 (PwMYMV-1) and 9,944 nt for variant 2 (PwMYMV-2). We propose that PwMYMV be assigned into the genus Potyvirus (family Potyviridae).

First report of Nerine yellow stripe virus infecting crinium lily (Crinium sp.) in Texas.

Crinum sp. (family Amaryllidaceae) is an ornamental flower bulb that is commonly called crinum lily, cape lily, cemetery plant, spider lily, and swamp lily. In April 2023, two plants of Crinum sp. var. Maiden's Blush with yellow stripe symptoms (Fig. S1) were submitted to the Texas Plant Virus Diagnostic Laboratory, Weslaco, TX for virus diagnosis. Due to the resemblance of the observed symptoms to those described for potyviruses infecting ornamental flower bulbs (Pearson et al. 2009), total RNA extracts were made from each sample using the SpectrumTM Plant Total RNA Kit (Sigma-Aldrich, USA), according to the manufacturer's protocol. Complementary DNA (cDNA) was synthesized from 2 µg total RNA per sample with Oligo(dT) primers using the PrimeScript™ 1st strand cDNA Synthesis Kit (Takara Bio, USA) as recommended by the manufacturer. A 2µL aliquot of each cDNA template was initially subjected to PCR using the generic primer pair CIFor/CIRev (Ha et al., 2008) that targets a fragment of the cylindrical inclusion (CI) body of potyviruses. The expected ~700 bp DNA band was amplified from both samples using the Taq DNA polymerase, dNTPack kit (Sigma-Aldrich). The amplicons were cloned and sequenced (three recombinant clones per sample) as described by Hernandez et al. (2021) and the BLASTX analyses of the consensus sequence (GenBank acc. no. OR137018) returned significant hits only to nerine yellow stripe virus (NeYSV; Potyvirus, Potyviridae) at 100% query coverage. To further confirm the results, another pair of universal primers (Jordan et al. 2011) was used to amplify the expected ∼1,600 bp product specific to the partial nuclear inclusion body (NIb), coat protein (CP) cistron, and 3' untranslated region of potyviruses from the same samples. The amplicons were similarly cloned, and a consensus sequence obtained (OR137019). In pairwise comparisons, the partial CI sequence of NeYSV from Texas (NeYSV-TX; OR137018) shared 83% nucleotide (nt)/93% amino acids (aa) identities with the corresponding sequences of NeYSV isolate 63 (MT396083) from the United Kingdom. The partial (649 nt) NIb sequences of NeYSV-TX (OR137019) and the complete CP (OR137019) of NeYSV-TX shared 77-94%/88-94% and 83-99%/89-98% nt/aa identities with the corresponding sequences of global NeYSV isolates that were retrieved from GenBank. Phylogenetic analysis revealed a closer relationship between NeYSV-TX and the isolates Stenomesson (EU042758) and DC (MG012805) from the Netherlands and USA, respectively based on the partial NIb and CP cistrons (Fig. S2), suggesting that NeYSV-TX may have been introduced from foreign and/or domestic sources. NeYSV has been documented previously from the United Kingdom, the Netherlands, Australia, New Zealand, and India; its first report from the United States was a decade ago from Amaryllis belladonna in California (Guaragna et al. 2013). To the best of our knowledge, this is the first report of NeYSV in Texas, thus expanding the geographical range of the virus in the USA. Anecdotal information from the sample submitter implicated infected crinum lily bulbs as the likely source of NeYSV introduction into the property, with subsequent vegetative propagation of plants resulting in 100% incidence of symptomatic lilies (n>100) over time. Thus, the results underscore the importance of ensuring that only virus-free vegetative plant materials are distributed and propagated by florists to curtail virus spread.

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 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.

First report of Ethiopian tobacco bushy top virus and its associated satellite RNA infecting common bean (Phaseolus vulgaris L.) in Zambia.

During surveys for common bean viruses in Central Province of Zambia in April 2018, symptoms of bushy top, deep green curled branches and patchy leaf chlorosis were observed on five plants in a 2-ha farmer's field. Total RNA was isolated from symptomatic leaf samples using the CTAB method (Chang et al. 1993). The RNA from one sample (CP414-1) was used to construct a cDNA library with the Illumina TruSeq RNA Library Prep Kit (Illumina, San Diego, CA), followed by high-throughput sequencing (HTS) on the Illumina MiSeq platform that generated ~3.1M single-end raw reads of ~300 nucleotides (nt) each. A total of 355,885 reads showed hits to Ethiopian tobacco bushy top virus (ETBTV; Umbravirus), ETBTV satellite RNA (satRNA) and peanut mottle virus (PeMoV, Potyvirus) based on BLASTn analysis. The full-length genomes of ETBTV (4239-nt; MT225089), its satRNA (521-nt; MT225092) and PeMoV (9,643-nt) were assembled from the HTS reads using Geneious R11.1.2 (Biomatters, Auckland, New Zealand). The obtained complete genome sequences of ETBTV (MT225089) and ETBTV satRNA (MT225092) shared 88% and 95% nt identities, respectively with the corresponding viral (KJ918748) and satRNA (KJ918747) sequences of isolate 18-2 (Abraham et al. 2014). The near complete PeMoV genome was 89% identical to isolate Liaoning (MH270528). The HTS results were validated by two-step RT-PCR analyses of the five field-collected samples using newly designed primer pairs (data not shown). All five samples gave the expected 988-bp ETBTV-specific and 521-bp satRNA-specific DNA bands while three samples produced the expected 2100-bp PeMoV-specific fragment. The virus specificities of the agent specific PCR fragments were ascertained by Sanger sequencing (ETBTV: MT225090-91; ETBTV satRNA: MT225093-94; PeMoV: MT900843-44) and they shared 98-100% identities with their corresponding HTS-derived sequences. To further probe for the presence of an ETBTV helper virus, the samples were screened by RT-PCR with the degenerate primer pair Lu1-mod-F/C2R3 that was modified from Robertson et al. (1991). The expected 245-bp DNA bands was obtained from all five samples, indicating the presence of a possible luteovirus or polerovirus target in these samples. The BLASTn analyses of the two Sanger sequenced gel-eluted products (MT900845-46) showed that they shared 100% identity with each other and 96% nt identity with cowpea polerovirus 1 (CPPV1, KX599163). Leaf tissue extracts from a common bean plant that was confirmed by RT-PCR to be positive for all four agents were rub-inoculated onto Nicotiana occidentalis and common bean (Sutter Pink) plants (n=5 each) at the three fully expanded leaf stage, with a buffer inoculation as control. Systemic foliar symptoms consisting of leaf deformation, stunting and leaf bushy top were observed on all ten plants, 10 days post-inoculation whereas the control plants remained symptomless. All the test plants were screened by RT-PCR as described above. The results showed that all five N. occidentalis plants were positive for ETBTV+ETBTVsatRNA, the five common bean plants tested positive for ETBTV+satRNA+PeMoV, and all 10 plants of both species were negative for CPPV1. To the best of our knowledge, this is the first report of ETBTV, ETBTV satRNA and CPPV1 infecting common bean in Zambia, and the first molecular based confirmation of PeMoV occurrence in the country. Ongoing studies are focused on determining the extent of the disease spread and assessment of its economic impact.