Complete Genomic RNA Sequence of Tuberose Mild Mosaic Virus and Tuberose Mild Mottle Virus Acquired by High-Throughput Sequencing

Abstract

Tuberose (Polianthes tuberosa) is an ornamental flowering crop of the Amaryllidaceae family. Tuberose mild mosaic virus (TuMMV) and tuberose mild mottle virus (TuMMoV), members of the genus Potyvirus, are ubiquitously distributed in most tuberose growing countries worldwide with low biological incidence. Here, we report the first coding-complete genomic RNA of TuMMV and TuMMoV obtained through high-throughput sequencing (HTS) and further, the presence of both the viruses were confirmed using virus-specific primers in RT-PCR assays. Excluding the poly (A) tail, the coding-complete genomic RNA of TuMMV and TuMMoV was 9485 and 9462 nucleotides (nts) in length, respectively, and contained a single large open reading frame (ORF). Polyprotein encoded by both the viral genomes contained nine putative cleavage sites. BLASTn analysis of TuMMV and TuMMoV genomes showed 72.40–76.80% and 67.95–77% nucleotide sequence similarities, respectively, with the existing potyviral sequences. Phylogenetic analysis based on genome sequences showed that TuMMV and TuMMoV clustered in a distinct clade to other potyviruses. Further studies are required to understand the mechanism of symptom development, distribution, genetic variability, and their possible threat to tuberose production in India.

1. Introduction

Polianthes tuberosa, also known as tuberose, is an important perennial herb within Amaryllidaceae family. In the country’s tropical regions, it is frequently grown as an ornamental blooming plant for the development of long-lasting flower spikes. It is also known as “Rajanigandha” or “Nishigandha” in popular parlance. It is believed to have originated in Mexico. In India, it is popularly grown in Andhra Pradesh, Assam, Gujarat, Karnataka, Punjab, Rajasthan, Tamil Nadu, West Bengal, and parts of Uttar Pradesh. It is rich in fragrance constituents as well as secondary metabolites, including polyphenols, nonpolar volatiles, and benzenoid derivatives. Furthermore, the flower and bulb extracts exhibit anti-inflammatory, antispasmodic, diuretic, and emetic effects [1].

Potyviridae contains about 200 plant virus species that are now allocated to the twelve genera Arepavirus, Bevemovirus, Brambyvirus, Bymovirus, Celavirus, Ipomovirus, Macluravirus, Poacevirus, Roymovirus, Rymovirus, Tritimovirus, and Potyvirus [2]. Members of the family are distinguished by host range, vector, genomic features, and phylogeny, with species demarcation frequently based on sequence identity of the main ORF (or, if necessary, the CP-coding region) being <76% (nt) and <82% (aa) [3].

Potyviridae members can be monopartite or bipartite (flexuous filamentous particles 650–950 nm long and 11–20 nm wide), single-stranded positive-sense RNA genomes with sizes ranging from 8.2 kb to 11.3 kb, with an average size of 9.7 kb. Genomic RNA is translated into polyproteins that must be proteolytically processed to provide ten mature and one fusion protein required for replication and motility: P1 (modulator of replication and translation), helper component proteinase HC-Pro (aphid transmission and silencing suppression), P3 (movement virus and replication), P3N-PIPO (cell-to-cell movement), 6K1 (formation of replication vesicles), cytoplasmic inclusion, CI (replication and helicase involved in virus movement), 6K2 (formation of replication vesicles), genome-linked protein VPg (replication, translation, and virus movement), NIa-Pro (polyprotein processing), NIb (RNA-dependent RNA polymerase), and CP (aphid transmission, virion formation and virus movement) [3].

Some members have a narrow host range, whilst others can infect a wide range of hosts [4]. Several of the members cause disease symptoms in commercially significant crops, rendering this genus one of the largest thoroughly researched virus groups by plant virologists. Potyviruses are primarily spread by aphids in a non-persistent manner as well as by mechanical means. Vertical transmission of potyviruses through seeds has also been documented [5]. Considering the economic importance of tuberose, identifying and characterizing the disease-causing viruses of tuberose is a prerequisite for devising suitable management strategies.

The rapid development of high-throughput sequencing (HTS) technology in recent years has not only enabled the investigation of diseases with unknown etiology but also facilitated the identification and characterization of known/novel viruses [6,7,8,9,10]. Unlike traditional antibody and nucleic acid-based detection techniques, such as enzyme-linked immunosorbent assay (ELISA), hybridization and PCR assays, HTS technologies do not warrant any virus-specific primers/probes/antibodies for detection of target pathogens. Through bioinformatic analysis of HTS reads, viruses associated with a particular sample can be unveiled and their genomes be recovered [11]. Besides, HTS facilitates detection of multiple infection of a plant species, thereby, helping us to comprehensively capture the entire virome and viral variability [8,9,10,11,12]. Using HTS technology, this study reports the first coding-complete genomic RNA sequence of tuberose mild mosaic virus (TuMMV) and tuberose mild mottle virus (TuMMoV) from diseased samples in India.

2. Results and Discussion

Tuberose is an economically important, vegetatively propagated flower crop grown worldwide [13]. In 2021 and 2022, symptoms like chlorotic stripes in leaves that initiated from the base and yellowing of the whole leaf, characteristic of potyvirus infection were observed in tuberose plants grown at the Horticultural Research Centre, Sardar Vallabhbhai Patel University of Agriculture and Technology, Meerut. The disease incidence of 79% was recorded. To detect the presence of viral particles, twenty-four symptomatic leaf samples were collected and samples were observed under TEM, which revealed the presence of flexuous particles of length 778.89 nm (Figure 1). Morphology and size of the virus particle was similar to those of potyviruses [14]. An increasing number of viral isolates have been unraveled utilizing HTS around the world, allowing for major advancements in diagnostic tools [6,7,8,9,10,11,12,15]. Furthermore, HTS technology may provide detailed insights into the etiological and epidemiological correlations of viruses associated with diseased plant samples [8,9,10,11,12,15,16]. Thus, HTS was employed to identify the viral agent causing the observed disease in the present study.

RNA of pooled symptomatic leaf samples was sequenced using the Illumina HiSeq 2000 platform to identify viruses that may be associated with the symptoms. In the two libraries, the Illumina sequencing data generated was approximately 90 million 75 bp paired-end reads and after trimming, 90,380,608 reads (average length 75 bp) were acquired (Table S1). A total of 108,668 contigs were generated from samples of tuberose. BLASTx search against the nr database revealed the complete-genome sequence of tuberose mild mosaic virus (TuMMV) and tuberose mild mottle virus (TuMMoV). Copy numbers of TuMMV and TuMMoV were 11,011 and 108,535, respectively. In public domain, only five partial sequences of TuMMV isolates deposited from Taiwan, India, and USA [17,18], and eight partial sequences of TuMMoV isolates deposited from China, India, Mexico, Netherlands, and Taiwan [19,20] are available up to date.

In present study, the obtained TuMMV and TuMMoV genomes contained 9485 and 9462 nt, respectively excluding the poly(A) tail. TuMMV and TuMMoV genomic RNA contained a large ORF, encoding for a polyprotein of 3070 aa and 3081 aa, which began with an ATG codon at 138 and 164 nt and ended with UAG and UAA termination codon at 9350 and 9490 nt, respectively. The characteristic domains, conserved cleavage sites, and motifs specific to the genus Potyvirus [21,22,23] were predicted in the polyprotein sequence encoded by the genomes.

Furthermore, a PIPO ORF resulting from viral RNA polymerase slippage in the TuMMV and TuMMoV P3 cistrons was predicted downstream of the G_1-2A_6–7 motif at positions 2604–2859 nts and 2625–2888 nts, respectively. This motif resembled the highly conserved G_1-2A_6–7 motif found in other members of the family Potyviridae [24]. Like other potyviruses [25], TuMMV polyprotein is proteolytically processed into 10 mature proteins: P1, HC-Pro, P3, 6K1, CI, 6K2, NIa-VPg, NIa-Pro, Nib, and CP, at aa positions 298, 756, 1099, 1151, 1787, 1841, 2033, 2276, and 2798, respectively (Table S2). Similarly, TuMMoV polyprotein is proteolytically cleaved to yield 10 mature proteins: P1, HC-Pro, P3, 6K1, CI, 6K2, NIa-VPg, NIa-Pro, Nib, and CP, at aa positions 299, 757, 1100, 1152, 1788, 1841, 2033, 2276, and 2797, respectively (Table S3).

Furthermore, potyvirus-specific conserved motifs were discovered in proteins encoded by the genomes of both the viruses. The predicted motifs in TuMMV cleavage products are ²⁵⁶H-(X)₈-D-(X)_33-34-S²⁹⁸ in P1; ³⁴⁹K-I-T-C³⁵² (aphid transmission), ⁴⁹⁴C-D-N-Q-L-D⁴⁹⁹ (symptomatology), and ⁶⁰⁷P-T-K⁶⁰⁹ (aphid transmission) in HC-Pro; ¹⁴⁵⁵V-A-T-N-I-I-E-N-G-V-T-L¹⁴⁶⁶ (potential helicase activity) and ¹⁵⁰⁶G-R-V-G-R¹⁵¹⁰ in CI; ¹²³⁷G-(X)₂-G-X-G-K-S¹²⁴⁴, ²⁴⁶⁶F-T-A-A-P²⁴⁷⁰, ²⁴⁸⁰C-V-D-D²⁴⁸³, ²⁵⁸⁴G-N-N-S-G-Q-P-S-T-V-V-D-N-S-L-M-V²⁶⁰⁰ (RNA-dependent polymerase activity), ²⁶²⁸⁸GDD²⁶³⁰ (RNA-dependent polymerase activity), and ²⁶⁷⁰W-F-M-S²⁶⁷³ in NIb; ²⁸¹⁷DAG²⁸¹⁹ (aphid transmission) and ³⁰²³QMKAAA³⁰²⁸ in CP. Similarly, in TuMMoV, the predicted motifs in mature proteins are ²⁰⁶H-(X)₈-D-(X)_33-34-S²⁵⁰ in P1; ³⁵⁰K-I-T-C³⁵³ (aphid transmission), ⁴⁹⁵C-D-N-Q-L-D⁵⁰⁰ (symptomatology) in HC-Pro; ¹⁴⁵⁶V-A-T-N-I-I-E-N-G-V-T-L¹⁴⁶⁷ (potential helicase activity) and ¹⁵⁰⁷G-R-V-G-R¹⁵¹¹ in CI; ²⁴⁶⁶F-T-A-A-P²⁴⁷⁰, ²⁴⁸⁰C-V-D-D²⁴⁸³, ²⁵⁸⁴G-N-N-S-G-Q-P-S-T-V-V-D-N-S-L-M-V²⁶⁰⁰ (RNA-dependent polymerase activity), ²⁶²⁸GDD²⁶³⁰ (RNA-dependent polymerase activity), and ²⁶⁷⁰W-F-M-S²⁶⁷³ in NIb; ²⁸¹⁷DAG²⁸¹⁹ (aphid transmission) and ³⁰³³QMKAAA³⁰³⁸ in CP [26,27]. Few of these motifs are suggested to play critical roles in the transmission of potyviruses by aphids [28]. Despite the presence of all of these conserved motifs, aphid transmission necessitates species-specific interactions [29,30]. Further researches are needed to determine the aphid species transmitting TuMMV and TuMMoV.

In Blast analysis, coding-complete genomic RNA sequence of TuMMV shared 72.40–76.80% sequence identities with complete genomes of other potyviruses like pokeweed mosaic virus (KU133475), plum pox virus (OK562680), turnip mosaic virus (AB093615), leek yellow stripe virus (AJ307057), Japanese yam mosaic virus (KJ789139), Moroccan watermelon mosaic virus (MN688647), sunflower chlorotic mottle virus (GU181200), and sweet potato feathery mottle virus (MH763680) while that of TuMMoV shared 67.95–77.00% sequence identities with other viruses like Pennisetum mosaic virus (JX070146), turnip mosaic virus (AB093615), leek yellow stripe virus (AJ307057), jasmine virus T (KX398054), and achyranthes bidentata mosaic virus (MT648692). Pairwise sequence identities (%) of the coding-complete genomic RNA sequence of TuMMV (ON116187) shared nucleotide (nt) identities 48.10–59.0% and amino acid (aa) identities 41.30–61.80%, and TuMMoV (and ON219793) shared nt identities 47.40–59.0% and aa identities 42.0–62.30% with some other closely related potyvirus isolates reported globally (Table 1). According to ICTV species demarcation, these results indicate that they belong to the same species [3,4]. Since 76% nt sequence identity in the complete genome has been defined as the species demarcation criteria for potyviruses [31], the Indian isolates of TuMMV and TuMMoV can be considered members of a distinct potyvirus species.

To analyze the evolutionary relationship between TuMMV, TuMMoV, and 44 other potyviruses, a phylogenetic tree based on their complete genomes was constructed. The results showed that these potyviruses can be divided into two distinct groups, subgroup I and subgroup II. Further, subgroup I divided into three distinct groups, subgroup Ia, subgroup Ib, and subgroup Ic. TuMMV and TuMMoV genomes that shared 65.20% nt identity clustered together in a single subgroup Ia, and both the viruses were showed genetic relatedness with papaya ring spot virus and Moroccan watermelon mosaic virus (Figure 2).

RT-PCR assay using primers mentioned in Table 1 yielded an expected amplicon of ~700 bp (Figure S1) from 31 samples representing four cultivars, including samples used for HTS. The sequences obtained were deposited to NCBI GenBank with accession numbers ON116188-ON116195 and ON241024–ON241026. BLASTn analysis of the partial CP gene sequence of TuMMV showed 87.33–97.31% nt sequence identity with other TuMMV isolates available in GenBank while partial CP TuMMoV sequences of this study shared 87.33–97.31% nt identities with other TuMMoV isolates available in GenBank. Pairwise sequence identity comparison of the partial CP sequence of TuMMoV with similar sequences shared 95.80–98.40% nt and 96.0–99.56% aa identities with the global isolates (Figure 3B). The TuMMV isolate shared highest sequence identities (97.33% nt and 96.7% aa) with a Taiwanian (EF137178) isolate while TuMMoV shared highest sequence identities (98.40% nt and 99.56%) with a Chinese isolate (AJ581528) whereas the suggested threshold for species demarcation in the coat protein is 76% [3,4].

To better understand the genetic variability of the TuMMV isolates, a phylogenetic tree was constructed based on the partial CP coding sequences of TuMMV and TuMMoV from different geographical locations. In the phylogenetic tree, the TuMMV isolate of this study fell in sister clade to the TuMMV Taiwanian isolates while the TuMMoV isolate obtained in the current study formed a sister clade to the TuMMoV isolates of China (Figure 3A).

Phylogenetic analysis based on the genomic RNA sequence (Figure 2) and partial coat protein region (Figure 3) showed consistent clustering of isolates suggesting a close relationship between the Indian isolate and other TuMMV and TuMMoV isolates, supported by high posterior probability values.

3. Materials and Methods

3.1. Sample Collection and Electron Microscopy

Symptoms of viral disease were observed in the experimental field of tuberose cv. Mexican-single, Hyderabad-single, Rajni-phule, and Sikkim-selection during January 2021 at the Horticulture Research Center, Sardar Vallabhbhai Patel University of Agriculture and Technology, Meerut, India. Leaf samples exhibiting typical symptoms of viral diseases such as mosaic, puckering, and stunting were collected and stored at −80 °C till further analysis along with asymptomatic samples. Sap obtained from symptomatic leaf samples were examined under transmission electron microscope (TEM) after staining with 2% aqueous uranyl acetate (UA) by following the leaf dip method [32].

3.2. RNA Extraction, Library Preparation and HTS

Leaf samples of Mexican-single and Hyderabad-single cultivar were subjected to total RNA extraction using GeneJET RNA purification kit (Thermo Scientific, MA, USA) in accordance with the user instructions manual. The quality and quantity of total RNA were assessed using Qubit™ fluorometer (Thermo Fisher Scientific, MA, USA) and TapeStation (Agilent, CA, USA). From the isolated pool of total RNA, ribosomal RNA was depleted using the Ribo-Zero rRNA Removal Kit (Illumina, San Diego, CA, USA). Sequencing libraries were prepared from the ribodepleted RNA pool using truseq-stranded-total RNA library preparation kit (Illumina, CA, USA) by following the manufacturer’s instructions. Bioanalyzer (Agilent, CA, USA) and TapeStation were used to verify the quality of the prepared libraries and the libraries were sequenced using Illumina HiSeq 2000 platform at NxGenBio Life Sciences (New Delhi, India) to produce 2 × 75 bp paired-end reads.

3.3. Data Processing and Virus Analysis

Raw reads obtained were imported into CLC Genomic workbench (20.0.4) and quality trimmed using the trimming tool to remove ambiguous and adaptor sequences. Trimmed reads were assembled using the de novo assembly tool in CLC Genomic workbench (20.0.4) to obtain longer contigs. Obtained contigs were subjected to BlastX analysis against non-redundant (NR) database in the workbench. The alignment against reference viral genomes in the database were made using OmicsBox 2.1 (https://www.biobam.com/omicsbox, accessed on 23 February 2022). The outcomes showed in a report have a number of assembled reads and total used reads. Open reading frames (ORFs) encoded by the putative viral genomes were determined by NCBI ORF Finder (https://www.ncbi.nlm.nih.gov/orffinder/, accessed on 2 March 2022). Conserved domains in the virus proteins were predicted using NCBI Conserved Domain-Search tool (https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi, accessed on 2 March 2022), and the genome organization of the identified viruses was constructed using the Bioedit 7.2 program [33].

3.4. Sequence Similarity and Phylogenetic Analysis

Genome sequences of identified viruses along with the available genome sequences of known potyviruses were retrieved and subjected to multiple sequence alignment using the CLUSTALW tool in MEGA-X [34]. Phylogenetic trees were constructed using neighbour-joining method and 1000 bootstrap replicates in MEGA-X. The nucleotide and amino acid sequence compositions and percent sequence identities were calculated using Bioedit 7.2 program [33] and visualized by the Sequence Demarcation Tool version 1.2 [35].

3.5. RT-PCR and Sanger Sequencing

Total RNA isolation was performed from four tuberose cvs. Mexican-single, Hyderabad-single, Rajni-phule, and Sikkim-selection as described in Section 3.2. cDNAs were synthesized from isolated RNA using random hexamers (1 µL) and RevertAid Reverse Transcriptase kit (Thermo Scientific, Waltham, MA, USA) according to the manufacturer’s protocol. PCR amplification of identified viruses was performed using virus-specific primers designed from the coat protein (CP) region of the obtained viral genomes (Table S4). The final 25 µL reaction volume contained 1x reaction buffer, 10 µM of each forward and reverse primers and 12.5 µL of DreamTaq PCR MasterMix (Thermo Scientific, Waltham, MA, USA). The reaction conditions were as follows: initial denaturation at 94 °C for 10 min, 35 cycles of 94 °C for 30 s, 52/55 °C for 45 s, 72 °C for 1 min followed by final extension at 72 °C for 10 min. PCR amplicons were purified and subjected to bidirectional Sanger sequencing at Centyle Biotech Pvt., Ltd., New Delhi, India.

4. Conclusions

Virus diseases limit tuberose production worldwide. With increased global trade, the possibility of introduction of tuberose viruses into newer areas grows by the day. As a result, rapid and reliable tuberose virus detection techniques are required to halt this inadvertent introduction and ensure virus-free tuberose production. We used HTS in this study and determined the first coding-complete genome sequences of TuMMV and TuMMoV. The genome sequences would be useful in providing a clearer picture of the virus evolution and in inferring biological properties of the viruses, which ultimately would help in devising management strategies. However, biological studies are required to understand the mode of transmission and effects of single/mixed viral infection on tuberose plants.

Author Contributions

Conceptualization J.S. and V.K.B.; software M.R.P. and V.K.S.; validation M.R.P., A.M. and P.S.; formal analysis M.R.P. and V.K.S.; investigation M.R.P. and J.S.; writing—original draft preparation M.R.P. and J.S.; writing—review and editing M.R.P., J.S., U.S., M.K., M.K.B., S.P., P.K. and R.K.; supervision J.S. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The HTS data used in this study have been deposited under the NCBI SRA under Bioproject: PRJNA810951 (SRA accessions: SRR18162840 and SRR18162840). The complete genomic RNA sequence of TuMMV and TuMMoV have been deposited in NCBI under accession number ON116187 and ON219793, respectively. Partial CP gene sequences have been deposited under accession numbers ON116188-ON116195 and ON241024-ON241026.

Acknowledgments

The authors acknowledge the Vice Chancellor, Sardar Vallabhbhai Patel University of Agriculture and Technology, Meerut-250110, Uttar Pradesh, Council of Science & Technology, Uttar Pradesh and Bioinformatics facility, Department of Biotechnology, India, for providing financial support and the facilities to carry out this research work.

Conflicts of Interest

All authors declare that they have no conflict of interest.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figures and Table

Figure 1. Transmission electron micrograph of symptomatic leaves tuberose.

Figure 2. Complete genome-based phylogenetic tree using the neighbor joining algorithm with p-distance matrix and 1000 bootstrap replicates. The scale bar represents a genetic distance of 0.1. Viral genome sequences of the current study is highlighted in red. Ryegrass mosaic virus was used as an outgroup.

Figure 3. Phylogenetic analysis of TuMMV and TuMMoV based on partial coat protein sequences. (A) Phylogenetic tree of TuMMV and TuMMoV isolates and other viruses worldwide were constructed using partial CP sequences using neighbour-joining algorithm and p-distance matrix with 1000 bootstrap replicates. (B) Pairwise identity plots of partial CP sequences of TuMMV and TuMMoV aligned by Clustal W and visualized using Sequence Demarcation Tool software. Colors on the scale denote the pairwise identities on colored heat map. Partial CP sequences of the current study is highlighted in red.

Table 1

Comparisons of nucleotide sequence (nts) and amino acid (aa) identity of pairwise combinations of complete genomic RNA sequences of TuMMV and TuMMoV with other potyviruses.

Sequence Identitynts/aa	ON116187 Tuberose mild mosaic	ON219793 Tuberose mild mottle virus	NC001785 Papaya ringspot virus	NC009995 Moroccan watermelon mosaic virus	NC001616 Potato virus Y	NC001517 Pepper mottle viru	NC000947 Japanese yam mosaic virus	NC011541 Narcissus yellow stripe virus	NC002509 Turnip mosaic virus	NC003605 Lettuce mosaic virus	NC001445 Plum pox virus	NC001841 Sweet potato feathery mottle virus	NC004039 Potato virus A	NC001768 Tobacco vein mottling virus	NC004011 Leek yellow stripe virus	NC004752 Yam mosaic virus	NC001555 Tobacco etch virus	NC040802 Lily yellow mosaic virus	NC003536 Clover yellow vein virus	NC003492 Bean yellow mosaic virus
ON116187 Tuberose mild mosaic virus	ID	71.0	61.8	46.0	43.0	43.3	45.4	45.1	45.4	41.7	45.7	42.8	43.6	43.7	41.3	44.5	43.4	42.6	42.0	43.2
ON219793 Tuberose mild mottle virus	65.2	ID	62.3	46.2	43.1	43.3	44.4	44.9	45.0	42.0	45.1	42.3	43.7	43.9	41.7	43.6	43.3	43.0	43.0	43.5
NC001785 Papaya ringspot virus	58.1	57.9	ID	57.9	41.0	41.0	42.8	42.3	42.5	43.0	44.0	42.3	41.9	41.3	44.0	42.2	49.0	42.3	47.0	41.3
NC009995 Moroccan watermelon mosaic virus	49.5	49.0	61.0	ID	43.8	43.7	45.0	45.3	45.8	43.8	47.6	44.3	43.8	44.1	42.9	45.8	44.3	44.1	43.9	44.5
NC001616 Potato virus Y	54.0	55.0	48.1	49.1	ID	61.3	44.8	45.3	45.9	45.6	48.0	43.5	44.7	45.2	42.0	47.2	45.7	43.2	42.8	43.8
NC001517 Pepper mottle virus	59.0	59.0	47.8	48.5	61.7	ID	45.2	45.7	46.3	45.0	47.6	43.6	44.6	44.3	42.1	46.3	45.8	43.3	43.4	43.7
NC000947 Japanese yam mosaic virus	50.0	48.8	49.7	49.9	49.5	49.9	ID	58.8	58.6	46.9	50.1	46.9	47.0	46.8	43.4	48.1	46.3	45.2	45.2	46.1
NC011541 Narcissus yellow stripe virus	49.2	48.3	48.7	49.3	49.8	49.2	61.9	ID	60.8	46.3	50.8	47.6	47.5	46.8	43.7	47.3	46.4	45.6	45.9	45.9
NC002509 Turnip mosaic virus	49.2	48.6	49.0	51.0	49.9	49.9	62.0	63.0	ID	47.4	50.9	49.3	48.5	46.8	44.1	47.8	46.8	45.3	45.2	45.1
NC003605 Lettuce mosaic virus	47.9	47.4	49.0	49.2	49.3	48.4	54.1	53.5	54.1	ID	47.7	45.3	44.1	43.9	42.4	46.1	43.7	42.7	43.0	42.9
NC001445 Plum pox virus	59.0	50.0	51.0	51.2	51.1	59.0	52.9	51.7	52.5	51.4	ID	50.3	49.1	47.0	46.2	50.3	47.3	47.8	46.4	47.4
NC001841 Sweet potato feathery mottle virus	48.9	48.9	54.0	56.0	49.1	49.7	51.4	51.4	51.1	54.0	54.2	ID	45.6	43.6	42.4	46.0	44.2	43.6	42.9	44.0
NC004039 Potato virus A	48.4	49.3	48.5	49.4	49.1	48.5	58.0	51.3	51.1	49.6	52.2	55.8	ID	54.6	43.5	46.9	48.5	45.8	45.5	46.4
NC001768 Tobacco vein mottling virus	49.6	49.4	47.8	48.3	53.0	49.2	49.1	49.1	48.7	48.9	49.8	49.6	52.1	ID	43.9	46.7	49.3	45.4	44.8	44.9
NC004011 Leek yellow stripe virus	49.6	49.2	48.7	48.5	49.9	49.7	49.4	48.6	49.3	48.5	51.1	49.6	49.3	54.6	ID	44.0	43.2	47.6	45.0	45.7
NC004752 Yam mosaic virus	48.8	47.9	49.0	49.5	49.1	49.0	57.0	49.4	54.0	49.3	52.0	51.5	51.3	48.9	48.5	ID	45.4	44.9	45.7	45.1
NC001555 Tobacco etch virus	48.4	47.7	47.6	48.6	49.3	48.7	49.8	49.3	49.7	48.6	51.2	58.0	51.6	49.8	48.9	54.8	ID	44.8	45.2	45.6
NC040802 Lily yellow mosaic virus	48.7	48.1	49.0	49.5	48.8	49.3	52.0	49.4	49.9	49.3	52.6	51.0	58.0	48.3	53.0	52.0	51.2	ID	49.3	49.4
NC003536 Clover yellow vein virus	48.8	48.2	48.7	49.4	48.6	49.2	50.0	49.5	56.0	48.6	51.4	53.0	58.0	48.5	49.7	51.9	52.2	55.8	ID	68.9
NC003492 Bean yellow mosaic virus	48.3	48.1	49.0	49.1	49.3	49.4	49.8	49.6	49.6	48.7	51.4	53.0	59.0	48.6	53.0	52.1	51.7	54.9	65.9	ID

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/pathogens11080861/s1. Table S1: HTS statistics for each sample before and after quality control showing the number of reads and the range of read lengths in base pairs (bp), Table S2: Putative genes along with nine protease cleavage sites in polyprotein of TuMMV, Table S3: Putative genes along with nine protease cleavage sites in polyprotein of TuMMoV, Table S4: List of virus-specific primers designed from the coat protein (CP) region of the obtained viral genomes, Figure S1: PCR detection of TuMMV and TuMMoV. Expected amplicon size was ~700 bp. Lane M, 1 kb DNA ladder (Cat. No. SM0311; Thermo Scientific); Lanes 2, 3, 4, 5, 6, 7, 8, 9, 10 and 11, PCR amplified products of TuMMV; Lane 12, 13, 14, 15, 16, 17, 18, 19, 20, 21 and 22, PCR amplified products of TuMMoV.

References

1. Verma, C.; Kashyap, S.; Kerketta, A. Major Diseases of Tuberose and Their Management. Pop. Kheti; 2020; 8, pp. 74-76.

2. Pasin, F.; Daròs, J.A.; Tzanetakis, I.E. Proteome expansion in the Potyviridae evolutionary radiation. FEMS Microbiol. Rev.; 2022; 46, pp. 1-22. [DOI: https://dx.doi.org/10.1093/femsre/fuac011] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35195244]

3. Inoue-Nagata, A.K.; Jordan, R.; Kreuze, J.; Li, F.; López-Moya, J.J.; Mäkinen, K.; Ohshima, K.; Wylie, S.J. ICTV Virus Taxonomy Profile: Potyviridae. J. Gen. Virol.; 2022; 103, 001738.

4. Wylie, S.J.; Adams, M.; Chalam, C.; Kreuze, J.; López-Moya, J.J.; Ohshima, K.; Praveen, S.; Rabenstein, F.; Stenger, D.; Wang, A. et al. ICTV virus taxonomy profile: Potyviridae. J. Gen. Virol.; 2017; 98, 352. [DOI: https://dx.doi.org/10.1099/jgv.0.000740]

5. Simmons, H.E.; Dunham, J.P.; Zinn, K.E.; Munkvold, G.P.; Holmes, E.C.; Stephenson, A.G. Zucchini yellow mosaic virus (ZYMV, Potyvirus): Vertical transmission, seed infection and cryptic infections. Virus Res.; 2013; 176, pp. 259-264. [DOI: https://dx.doi.org/10.1016/j.virusres.2013.06.016] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23845301]

6. Hadidi, A.; Flores, R.; Candresse, T.; Barba, M. Next-generation sequencing and genome editing in plant virology. Front. Microbiol.; 2016; 7, 1325. [DOI: https://dx.doi.org/10.3389/fmicb.2016.01325] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27617007]

7. Pecman, A.; Kutnjak, D.; Gutiérrez-Aguirre, I.; Adams, I.; Fox, A.; Boonham, N.; Ravnikar, M. Next generation sequencing for detection and discovery of plant viruses and viroids: Comparison of two approaches. Front. Microbiol.; 2017; 8, 1998. [DOI: https://dx.doi.org/10.3389/fmicb.2017.01998] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29081770]

8. Sidharthan, V.K.; Kalaivanan, N.S.; Baranwal, V.K. Discovery of putative novel viruses in the transcriptomes of endangered plant species native to India and China. Gene; 2021; 786, 145626. [DOI: https://dx.doi.org/10.1016/j.gene.2021.145626]

9. Prajapati, M.R.; Manav, A.; Singh, J.; Singh, M.K.; Ranjan, K.; Kumar, A.; Kumar, P.; Kumar, R.; Baranwal, V.K. Identification of Garlic virus A infecting Allium sativum L. through next generation sequencing technology. J. Hortic. Sci. Biotechnol.; 2022; 97, pp. 96-105. [DOI: https://dx.doi.org/10.1080/14620316.2021.1963848]

10. Villamor, D.E.V.; Ho, T.; Al Rwahnih, M.; Martin, R.R.; Tzanetakis, I.E. High throughput sequencing for plant virus detection and discovery. Phytopathology; 2019; 109, pp. 716-725. [DOI: https://dx.doi.org/10.1094/PHYTO-07-18-0257-RVW] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30801236]

11. Singh, J.; Truong, T.N.; An, D.; Prajapati, M.R.; Manav, A.; Quoc, N.B.; Ranjan, K.; Kumar, A.; Kumar, P.; Kumar, R. et al. Complete genome sequence and genetic organization of a Garlic virus D infecting garlic (Allium sativum) from northern India. Acta Virol.; 2020; 64, pp. 427-432. [DOI: https://dx.doi.org/10.4149/av_2020_405] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33112640]

12. Prajapati, M.R.; Manav, A.; Singh, J.; Kumar, P.; Kumar, A.; Kumar, R.; Prakash, S.; Baranwal, V.K. Identification and Characterization of a Garlic Virus E Genome in Garlic (Allium sativum L.) Using High-Throughput Sequencing from India. Plants; 2022; 11, 224. [DOI: https://dx.doi.org/10.3390/plants11020224]

13. Kainthura, P.; Srivastava, R. Induction of genetic variability and isolation of mutants in tuberose (Polianthes tuberosa L.). Trop. Agric. Res.; 2015; 26, pp. 721-732. [DOI: https://dx.doi.org/10.4038/tar.v26i4.8135]

14. Harrison, B.D.; Finch, J.T.; Gibbs, A.J.; Hollings, M.; Shepherd, R.J.; Valenta, V.; Wetter, C. Sixteen groups of plant viruses. Virology; 1971; 45, pp. 356-363. [DOI: https://dx.doi.org/10.1016/0042-6822(71)90336-9]

15. Adams, I.P.; Miano, D.W.; Kinyua, Z.M.; Wangai, A.; Kimani, E.; Phiri, N.; Reeder, R.; Harju, V.; Glover, R.; Hany, U. et al. Use of next-generation sequencing for the identification and characterization of M aize chlorotic mottle virus and S ugarcane mosaic virus causing maize lethal necrosis in Kenya. Plant Pathol.; 2013; 62, pp. 741-749. [DOI: https://dx.doi.org/10.1111/j.1365-3059.2012.02690.x]

16. Perotto, M.C.; Pozzi, E.A.; Celli, M.G.; Luciani, C.E.; Mitidieri, M.S.; Conci, V.C. Identification and characterization of a new potyvirus infecting cucurbits. Arch. Virol.; 2018; 163, pp. 719-724. [DOI: https://dx.doi.org/10.1007/s00705-017-3660-2] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29196817]

17. Kulshrestha, S.; Mehra, A.; Hallan, V.; Raikhy, G.; Ram, R.; Zaidi, A.A. Molecular evidence for occurrence of Tuberose mild mottle virus infecting tuberose (Polianthes tuberosa) in India. Curr. Sci.; 2005; 89, pp. 870-872.

18. Raj, S.K.; Snehi, S.K.; Kumar, S.; Ram, T.; Goel, A.K. First report of Tuberose mild mosaic potyvirus from tuberose (Polianthes tuberosa L.) in India. Australas. Plant Dis. Notes; 2009; 4, pp. 93-95.

19. Dey, K.K.; Melzer, M.J.; Li, C.; Sun, X.; Adkins, S. First report of tuberose mild mottle virus infecting tuberose (Polianthes tuberosa) in the United States. Plant Disease; 2018; 102, 461. [DOI: https://dx.doi.org/10.1094/PDIS-07-17-1024-PDN]

20. Lin, L.; Zheng, H.Y.; Chen, J.; Chen, J.P.; Zhang, Q.Y.; Zhao, M.F.; Antoniw, J.F.; Adams, M.J. A new potyvirus from tuberose (Polianthes tuberosa) in China. Arch. Virol.; 2004; 149, pp. 1107-1116. [DOI: https://dx.doi.org/10.1007/s00705-003-0289-0] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15168198]

21. Lucinda, N.; Da Rocha, W.B.; Inoue-Nagata, A.K.; Nagata, T. Complete genome sequence of pepper yellow mosaic virus, a potyvirus, occurring in Brazil. Arch. Virol.; 2012; 157, pp. 1397-1401. [DOI: https://dx.doi.org/10.1007/s00705-012-1313-z] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22527869]

22. Knierim, D.; Menzel, W.; Winter, S. Analysis of the complete genome sequence of euphorbia ringspot virus, an atypical member of the genus Potyvirus. Arch. Virol.; 2017; 162, pp. 291-293. [DOI: https://dx.doi.org/10.1007/s00705-016-3087-1] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27686070]

23. Chung, B.Y.W.; Miller, W.A.; Atkins, J.F. & Firth, A.E. An overlapping essential gene in the Potyviridae. Proc. Natl. Acad. Sci. USA; 2008; 105, pp. 5897-5902. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18408156]

24. Olspert, A.; Chung, B.Y.W.; Atkins, J.F.; Carr, J.P.; Firth, A.E. Transcriptional slippage in the positive-sense RNA virus family Potyviridae. EMBO Rep.; 2015; 16, pp. 995-1004. [DOI: https://dx.doi.org/10.15252/embr.201540509]

25. Nigam, D.; LaTourrette, K.; Souza, P.F.; Garcia-Ruiz, H. Genome-wide variation in potyviruses. Front. Plant Sci.; 2019; 10, 1439. [DOI: https://dx.doi.org/10.3389/fpls.2019.01439]

26. Lan, P.; He, P.; Zhang, Y.; Zhang, S.; Zhang, Z.; Chen, X.; Tan, S.; Luo, H.; Cao, M.; Li, F. Molecular characterization of a novel potyvirus infecting noni. Arch. Virol.; 2019; 164, pp. 3099-3102. [DOI: https://dx.doi.org/10.1007/s00705-019-04400-z] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31520219]

27. Worrall, E.A.; Hayward, A.C.; Fletcher, S.J.; Mitter, N. Molecular characterization and analysis of conserved potyviral motifs in bean common mosaic virus (BCMV) for RNAi-mediated protection. Arch. Virol.; 2019; 164, pp. 181-194. [DOI: https://dx.doi.org/10.1007/s00705-018-4065-6] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30302583]

28. Revers, F.; García, J.A. Molecular biology of potyviruses. Adv. Virus Res.; 2015; 92, pp. 101-199. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25701887]

29. Gibbs, A.J.; Mackenzie, A.M.; Wei, K.J.; Gibbs, M.J. The potyviruses of Australia. Arch. Virol.; 2008; 153, pp. 1411-1420. [DOI: https://dx.doi.org/10.1007/s00705-008-0134-6] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18566735]

30. Pelletier, Y.; Nie, X.; Giguère, M.A.; Nanayakkara, U.; Maw, E.; Foottit, R. A new approach for the identification of aphid vectors (Hemiptera: Aphididae) of potato virus Y. J. Econ. Entomol.; 2012; 105, pp. 1909-1914. [DOI: https://dx.doi.org/10.1603/EC12085]

31. Romay, G.; Lecoq, H.; Desbiez, C. Zucchini tigré mosaic virus is a distinct potyvirus in the papaya ringspot virus cluster: Molecular and biological insights. Arch. Virol.; 2014; 159, pp. 277-289. [DOI: https://dx.doi.org/10.1007/s00705-013-1798-0] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23979176]

32. Manav, A.; Prajapati, M.R.; Singh, J.; Kumar, A.; Kumar, P.; Pant, R.P.; Baranwal, V.K. First report of natural infection by two potyviruses on amaryllis (Hippeastrum hybridum) plants from India. VirusDisease; 2021; 32, pp. 830-833. [DOI: https://dx.doi.org/10.1007/s13337-021-00741-z] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34901330]

33. Hall, T.A. BioEdit Sequence Alignment Editor 7.0. 1; Isis Pharmaceuticals: Carlsbad, CA, USA, 2004.

34. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol.; 2018; 35, 1547. [DOI: https://dx.doi.org/10.1093/molbev/msy096] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29722887]

35. Muhire, B.M.; Varsani, A.; Martin, D.P. SDT: A virus classification tool based on pairwise sequence alignment and identity calculation. PLoS ONE; 2014; 9, e108277. [DOI: https://dx.doi.org/10.1371/journal.pone.0108277] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25259891]