Thevetia ahouai L. is an evergreen shrub with shiny, dark green, ovate leaves and bright, lobed red fruits with a milky sap. The shrub belongs to the Apocynaceae family, and is native from Brazil; although it can be found all the way to Mexico [3]. T. ahouai has been used in traditional medicine as treatment for hemorrhoids, toothache, and rheumatism, and has also shown to have anti-promastigote activity against Leishmania [6]. In recent years, the ornamental use of T. ahouai has expanded to private gardens and public places, such as urban parks or sidewalks, due to its colorful appearance.
In June of 2020, virus-like symptoms including leaf white spot and fruit discoloration were observed in T. ahouai plants at two different locations of Guayaquil, a coastal city of Ecuador. Symptomatic leaves from two selected shrubs, one located in an urban park (sample 1) and the second located in Prosperina, a tropical dry forest in the western side of the city (sample 2), were collected for virus identification.
Due to the lack of reports on viruses infecting Thevetia sp., collected leaves were submitted to a partial virus purification protocol as described before [4]. Partially purified extracts were mounted on a carbon-coated formvar (1%) grid, and negatively stained with 2% PTA (phosphotungstic acid at pH 7.0). Grids were examined using a JEOL JEM-1400Plus transmission electron microscope hosted at the University of Minnesota Imaging Center. Flexuous filamentous virus-like particles of approximately 700 nm in length were identified in partially purified extracts from both samples (Online Resource 1).
High-throughput sequencing (HTS) was applied for virus identification using double-stranded RNA (dsRNA) for sample 1, or total RNA for sample 2, as initial template. The dsRNA was extracted from 15 g of fresh symptomatic leaf tissue following the protocol described by Morris and Dodds [7]. Total RNA was extracted from 100 mg of symptomatic leaf tissue using the RNeasy Plant Mini Kit (Qiagen, Germany) and subjected to plant ribosomal RNA (rRNA) depletion using Illumina’s Ribo-Zero kit prior to the generation of a complementary DNA (cDNA) library using the TruSeq library prep kit. Sequencing was done on NovaSeq6000 Illumina platform as 150bp paired-end reads (Macrogen, South Korea).
A total of 22.3 and 21.2 million sequence reads were obtained for samples 1 and 2, respectively. Sequence data sets were analyzed using HTS-processing tools available from Geneious Prime® 2022.0.1. Raw sequences were trimmed for adapter removal and quality using the BBDuk plugin and de novo assembled using SPAdes. Several thousands of contigs were assembled from each sequence set. Blastx search identified a 9,528 nt long contig from sample 1, and a 9,542 nt long contig from sample 2, both showing sequence homology to several members of the potyvirus genus. A closer examination of each contig revealed that 15,844 reads (0.07%) were assembled into the potyvirus contig from sample 1; whereas 454,852 reads (2.15 %) were assembled into the potyvirus contig from sample 2. Blast searches of the remaining contigs revealed their host origin, which is expected as the RNA ribosomal depletion step does not subtract the whole host RNA. As for the difference in coverage observed between the viral contigs obtained from the two template types, dsRNA and total RNA, it is known that potyviruses yield low amounts of dsRNA compared to other plant viruses [12]. Given the low coverage of contig 1, a series of overlapping primers were designed to amplify and re-sequence contig 1 by cloning each RT-PCR fragment using a pGEM-T-easy Vector System (Promega, USA) followed by Sanger sequencing.
The original template was used to obtain the 5’ and 3’ terminal sequences using the 5´/3´ RACE Kit, 2nd Generation (Roche, Germany), according to manufacturer instructions. For the dsRNA, an additional denaturation step (96 C for 10 min) was used prior the reverse-transcription reaction. RT-PCR amplicons for each terminus (n =5) were cloned as described above and sequenced in both directions.
The complete genomic sequence, excluding the poly (A) tail, consisted of 9,912 nt (acc. numb. OM263475) and 9,904 nt (acc. numb. OM263476), for potyvirus sequence 1 and 2, respectively. The identity percentages between the two sequences (80% nt level and 86.7% amino acid level) suggest that both sequences belong to isolates of the same potyvirus (hereafter isolate 1 and isolate 2) [15].
Blastn searches revealed that both isolates share up to ~ 73% nt identity (86% coverage) with their closest relative, a potyvirus sequence from weeds in a papaya orchard of Chiapas, Mexico (access. numb. MN203192). When the complete genome of this potyvirus was compared to both isolates of the thevetia potyvirus, the identity was 69%. According to the species demarcation criteria for potyviruses [15], at this genomic identity level, it can be inferred that the thevetia virus described in this study represents a new member of the potyvirus genus.
Symptomatic leaves infected with each virus isolate were used for mechanical inoculations onto T. ahouai virus-free seedlings as described [9]. Leaf white spots were observed at an average of 15 days post-inoculation, with no differences between the symptoms induced by each isolate. The presence of the virus in the inoculated plants was confirmed by RT-PCR and Sanger sequencing using the primers: Det_F: 5’-TCAGGAACGGTCTCGGTTCC-3’ and Det_R: 5’-CCATCATCACCCAAACTCCAT-3’, which amplify a 292 bp fragment of the virus coat protein (CP) gene. Inoculated plants were maintained under controlled conditions and symptoms were monitored for one year. Original symptoms including leaf white spots and fruit discoloration were reproduced in the inoculated plants, while no other virus-like sequence was found in the HTS data sets. Taken together, these findings suggest that the new potyvirus is the causal agent of the described symptomatology (Fig 1). Hence, the name thevetia white spot virus (ThWSV) is proposed and will be used hereafter.
In order to test seed-transmission, as reported for other potyviruses [11], five symptomatic fruits from a single virus-inoculated plant, were collected and all the twenty seeds (each fruit has 4 seeds) were potted in sterile germination medium. The third true leaf of each seedling was tested for the virus by RT-PCR as described above. None of the plants tested positive for the virus and no symptoms were observed during the study.
The genome organization of ThWSV is identical to those of recognized potyviruses, containing a long open reading frame (ORF) at nt positions 123-9,671 for isolate 1, and 118-9,663 for isolate 2. The hypothetical polyprotein precursor from isolate 1 has 3,183 amino acids with a predicted molecular mass of 363.5 kDa; while the polyprotein from isolate 2 has 3,182 aa (363.6 kDa), having a single lysine deletion (derived from codon AAA) at the N-terminus of the coat protein (CP) with respect to isolate 1.
Analyses of the hypothetical polyprotein of ThWSV identified the nine conserved proteolytic cleavage sites previously described for potyviruses [1], resulting in ten mature putative proteins (Fig 2A). In addition, typical potyvirus conserved motifs [14] were identified in the proteins of ThWSV, with slight differences in a few motifs between the two isolates (Table 1).
The putative small ORF termed PIPO was also identified in both isolates (Fig 2A), overlapping with the P3 coding region through the presence of the highly conserved motif G1–2A6–7 at the beginning of the PIPO ORF (isolate 1: 3324GAAAAAT3330; isolate 2: 3315GAAAAAT3325). In both isolates, the PIPO ORF is found in a non-frame fashion, suggesting its expression through a -1 ribosomal frameshifting from the P3 coding region, which would result in a fused protein (P3N-PIPO) as previously described [13].
Phylogenetic inferences were done using the complete polyprotein amino acid sequences of 36 representative potyvirus species and a member of the Ipomovirus genus (Potyviridae), which was used as outgroup. Multiple sequence alignment was done using MUSCLE [8] and the best fitted protein model (LG+G+F+I) was obtained. The phylogenetic tree was generated using the maximum-likelihood method with 500 bootstrap replicates in MEGAX [10]. The topology of the tree was consistent with previous Blast results, showing a most recent common ancestor for ThWSV, a potyvirus sequenced from an unknown weed in Mexico and asclepias virus A, a perennial herb in the same family (Apocynaceae) as T. ahouai. Other closely related potyviruses include pokeweed mosaic virus, tobacco vein mottling, potato virus A and potato virus B (Fig 2B). Amino acid identities among most closely related species ranged from 51% to 72%.
To the best of our knowledge, this is the first report of a virus infecting T. ahouai. ThWSV induced a range of symptoms including white spots on the leaves, darkening and black ringspots on the stems, and fruit discoloration (Fig 1). Further studies should be conducted to investigate the host range and natural vector of this new virus, especially due to its increased use as ornamental, which might pose a threat to cultivated plants. Based on the presence of aphid-transmission related motifs, such as KITC (KIAC in isolate 2) and PTK, in the helper component protein (HC-Pro), and DAG in the CP (Table 1), it is reasonable to speculate that natural transmission of ThWSV is mediated by aphids. However, epidemiology studies should focus on identifying those aphid species that are more common in T. ahouai and closely related species such as Catharanthus roseus (L.), which is ubiquitously found as ornamental in tropical regions. Sequence comparisons at nucleotide and amino acid level indicate that ThWSV is most closely related to a potyvirus found in an unknown weed growing within papaya orchards in Mexico [2] and asclepias virus A [5], isolated from Asclepias syriaca, a perennial herb in the same family (Apocynaceae) as T. ahouai.