First Report of Beet Western Yellows Virus Infecting Vasconcellea x heilbornii (Caricaceae) in Ecuador.

Abstract

Vasconcellea x heilbornii, known as babaco, is a hybrid native to Ecuador grown in small orchards in sub-tropical highland regions. Over the last decade, several viruses have been identified in babaco using high-throughput sequencing (HTS) (Cornejo-Franco et al. 2020, (Reyes-Proaño et al. 2023). In 2021, total RNA from a babaco plant showing distinctive leaf yellowing was extracted using the PureLink RNA Mini Kit (Thermo Fischer Scientific, USA) and subjected to HTS on an Illumina NovaSeq6000 system as 150 paired-end reads (Macrogen Inc., South Korea). Library construction was done using the TruSeq Stranded Total RNA Sample kit with Plant Ribo-Zero, as described (Villamor et al. 2022). Reads were processed using BBDuk and de novo assembled using SPAdes 3.15. both implemented in Geneious 2022. Contig analysis was done by BLASTx using the NCBI viral sequence database (as of November 2022). HTS generated 54 million reads, of which 12% assembled into contigs corresponding to genomes of previously reported babaco viruses including babaco virus Q (BabVQ), babaco nucleorhabdovirus 1 (BabRV1) and babaco ilarvirus 1 (BabIV1). Interestingly, 144 reads (0.0003%) assembled into seven contigs ranging from 100 to 480 nucleotides (nt) in length. These contigs showed homology, with 97% amino acid (aa) identity (100% query coverage), to regions of the RNA-dependent-RNA-polymerase (RdRp) of beet western yellows virus (BWYV, Acc. No. NC_004756), a member of the Polerovirus genus. To confirm the occurrence of BWYV in babaco, double-stranded RNA (dsRNA) was extracted from 15 g of leaf tissue from the original sample as described (Dodds et al. 1984) and used as template for reverse-transcription (RT)-PCR using overlapping primers designed to span all short contigs. RT-PCR amplified fragments were cloned into a pGEM®T-easy vector (Promega, USA) and sequenced by the Sanger method (Macrogen Inc., South Korea). The sequences were assembled into a single 2.7 kbp BWYV genome fragment comprising the complete protein 1 (P1) and partial RdRp gene (GenBank Acc. No. PP480670). Sequence alignments between the partially sequenced genome of the babaco isolate and its corresponding fragment from the closest BWYV isolate (NC_004756) revealed 94% and 97% identities at the nt and aa levels, respectively. To assess the prevalence of BWYV in babaco, 30 leaf samples showing yellowing symptoms from Pichincha (n=15) and Azuay (n=15) provinces were tested by RT-PCR using total RNA. Total RNA extraction and reverse transcription were done using the methodology described by Halgren et al. (2007). For RT-PCR, the primer set BWYV_Bab_F: 5'-CAGTGTCCTCCAAGTGCAACAT-3' / BWYV_Bab_R: 5'GGTTCCTTCCCAGTTTGGTGGT-3', which amplifies a 235 nt-long P1 region, was used. Three RT-PCR products from each positive sample were purified using the GeneJET PCR clean-up kit (Thermo Scientific, USA) and sequenced. BWYV was confirmed in 9 out of 15 samples (60%) from Pichincha, and in 10 out of 15 samples (64%) from Azuay. Samples were also tested for additional babaco viruses as described (Reyes-Proaño et al. 2023). All BWYV-infected plants turned out positive for papaya ringspot virus (PRSV), babaco mosaic virus (BabMV), BabVQ, and BabIV1. Hence, the impact of BWYV infection on babaco plants in single and mixed infections warrants further investigation. To the best of our knowledge, this is the first report of BWYV in a crop in Ecuador, and the first time it has been found in a Caricaceae species.