Next Article in Journal
Detection of Ticks and Tick-Borne Pathogens of Urban Stray Dogs in South Africa
Next Article in Special Issue
Unlocking the Hidden Genetic Diversity of Varicosaviruses, the Neglected Plant Rhabdoviruses
Previous Article in Journal
Safety and Efficacy of Direct Antiviral Agents for Hepatitis C in Patients with Malignancies Other Than Liver Cancer: A Case Series
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
College of Biotechnology, Sardar Vallabhbhai Patel University of Agriculture and Technology, Meerut 250110, Uttar Pradesh, India
2
Division of Plant Pathology, ICAR-Indian Agricultural Research Institute, New Delhi 110012, India
3
Division of Genetics and Tree Improvement, Institute of Forest Biodiversity (ICFRE), Hyderabad 500100, Telangana, India
4
College of Agriculture, Sardar Vallabhbhai Patel University of Agriculture and Technology, Meerut 250110, Uttar Pradesh, India
5
College of Horticulture, Sardar Vallabhbhai Patel University of Agriculture and Technology, Meerut 250110, Uttar Pradesh, India
6
College of Veterinary and Animal Science, Sardar Vallabhbhai Patel University of Agriculture and Technology, Meerut 250110, Uttar Pradesh, India
*
Author to whom correspondence should be addressed.
Pathogens 2022, 11(8), 861; https://doi.org/10.3390/pathogens11080861
Submission received: 4 July 2022 / Revised: 20 July 2022 / Accepted: 22 July 2022 / Published: 30 July 2022
(This article belongs to the Special Issue Plant Virus Genome Diversity in Plant Hosts and Insect Vectors)

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 G1-2A6–7 motif at positions 2604–2859 nts and 2625–2888 nts, respectively. This motif resembled the highly conserved G1-2A6–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 256H-(X)8-D-(X)33-34-S298 in P1; 349K-I-T-C352 (aphid transmission), 494C-D-N-Q-L-D499 (symptomatology), and 607P-T-K609 (aphid transmission) in HC-Pro; 1455V-A-T-N-I-I-E-N-G-V-T-L1466 (potential helicase activity) and 1506G-R-V-G-R1510 in CI; 1237G-(X)2-G-X-G-K-S1244, 2466F-T-A-A-P2470, 2480C-V-D-D2483, 2584G-N-N-S-G-Q-P-S-T-V-V-D-N-S-L-M-V2600 (RNA-dependent polymerase activity), 26288GDD2630 (RNA-dependent polymerase activity), and 2670W-F-M-S2673 in NIb; 2817DAG2819 (aphid transmission) and 3023QMKAAA3028 in CP. Similarly, in TuMMoV, the predicted motifs in mature proteins are 206H-(X)8-D-(X)33-34-S250 in P1; 350K-I-T-C353 (aphid transmission), 495C-D-N-Q-L-D500 (symptomatology) in HC-Pro; 1456V-A-T-N-I-I-E-N-G-V-T-L1467 (potential helicase activity) and 1507G-R-V-G-R1511 in CI; 2466F-T-A-A-P2470, 2480C-V-D-D2483, 2584G-N-N-S-G-Q-P-S-T-V-V-D-N-S-L-M-V2600 (RNA-dependent polymerase activity), 2628GDD2630 (RNA-dependent polymerase activity), and 2670W-F-M-S2673 in NIb; 2817DAG2819 (aphid transmission) and 3033QMKAAA3038 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.

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.

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.

Funding

This research received no external funding.

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.

References

  1. Verma, C.; Kashyap, S.; Kerketta, A. Major Diseases of Tuberose and Their Management. Pop. Kheti 2020, 8, 74–76. [Google Scholar]
  2. Pasin, F.; Daròs, J.A.; Tzanetakis, I.E. Proteome expansion in the Potyviridae evolutionary radiation. FEMS Microbiol. Rev. 2022, 46, 1–22. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar]
  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. [Google Scholar] [CrossRef]
  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, 259–264. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Hadidi, A.; Flores, R.; Candresse, T.; Barba, M. Next-generation sequencing and genome editing in plant virology. Front. Microbiol. 2016, 7, 1325. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  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. [Google Scholar] [CrossRef]
  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, 96–105. [Google Scholar] [CrossRef]
  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, 716–725. [Google Scholar] [CrossRef] [PubMed]
  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, 427–432. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef]
  13. Kainthura, P.; Srivastava, R. Induction of genetic variability and isolation of mutants in tuberose (Polianthes tuberosa L.). Trop. Agric. Res. 2015, 26, 721–732. [Google Scholar] [CrossRef] [Green Version]
  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, 356–363. [Google Scholar] [CrossRef]
  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, 741–749. [Google Scholar] [CrossRef]
  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, 719–724. [Google Scholar] [CrossRef] [PubMed]
  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, 870–872. [Google Scholar]
  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, 93–95. [Google Scholar]
  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. [Google Scholar] [CrossRef]
  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, 1107–1116. [Google Scholar] [CrossRef] [PubMed]
  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, 1397–1401. [Google Scholar] [CrossRef] [PubMed]
  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, 291–293. [Google Scholar] [CrossRef] [PubMed]
  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, 5897–5902. [Google Scholar] [PubMed] [Green Version]
  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, 995–1004. [Google Scholar] [CrossRef]
  25. Nigam, D.; LaTourrette, K.; Souza, P.F.; Garcia-Ruiz, H. Genome-wide variation in potyviruses. Front. Plant Sci. 2019, 10, 1439. [Google Scholar] [CrossRef] [Green Version]
  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, 3099–3102. [Google Scholar] [CrossRef] [PubMed]
  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, 181–194. [Google Scholar] [CrossRef] [PubMed]
  28. Revers, F.; García, J.A. Molecular biology of potyviruses. Adv. Virus Res. 2015, 92, 101–199. [Google Scholar] [PubMed]
  29. Gibbs, A.J.; Mackenzie, A.M.; Wei, K.J.; Gibbs, M.J. The potyviruses of Australia. Arch. Virol. 2008, 153, 1411–1420. [Google Scholar] [CrossRef] [PubMed]
  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, 1909–1914. [Google Scholar] [CrossRef]
  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, 277–289. [Google Scholar] [CrossRef] [PubMed]
  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, 830–833. [Google Scholar] [CrossRef] [PubMed]
  33. Hall, T.A. BioEdit Sequence Alignment Editor 7.0. 1; Isis Pharmaceuticals: Carlsbad, CA, USA, 2004. [Google Scholar]
  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. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Transmission electron micrograph of symptomatic leaves tuberose.
Figure 1. Transmission electron micrograph of symptomatic leaves tuberose.
Pathogens 11 00861 g001
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 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.
Pathogens 11 00861 g002
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.
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.
Pathogens 11 00861 g003
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.
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 Identity
nts/aa
ON116187 Tuberose mild mosaicON219793 Tuberose mild mottle virusNC001785 Papaya ringspot virusNC009995 Moroccan watermelon mosaic virusNC001616 Potato virus YNC001517 Pepper mottle viruNC000947 Japanese yam mosaic virusNC011541 Narcissus yellow stripe virusNC002509 Turnip mosaic virusNC003605 Lettuce mosaic virusNC001445 Plum pox virusNC001841 Sweet potato feathery mottle virusNC004039 Potato virus ANC001768 Tobacco vein mottling virusNC004011 Leek yellow stripe virusNC004752 Yam mosaic virusNC001555 Tobacco etch virusNC040802 Lily yellow mosaic virusNC003536 Clover yellow vein virusNC003492 Bean yellow mosaic virus
ON116187 Tuberose mild mosaic virusID71.061.846.043.043.345.445.145.441.745.742.843.643.741.344.543.442.642.043.2
ON219793 Tuberose mild mottle virus65.2ID62.346.243.143.344.444.945.042.045.142.343.743.941.743.643.343.043.043.5
NC001785 Papaya ringspot virus58.157.9ID57.941.041.042.842.342.543.044.042.341.941.344.042.249.042.347.041.3
NC009995 Moroccan watermelon mosaic virus49.549.061.0ID43.843.745.045.345.843.847.644.343.844.142.945.844.344.143.944.5
NC001616 Potato virus Y54.055.048.149.1ID61.344.845.345.945.648.043.544.745.242.047.245.743.242.843.8
NC001517 Pepper mottle virus59.059.047.848.561.7ID45.245.746.345.047.643.644.644.342.146.345.843.343.443.7
NC000947 Japanese yam mosaic virus50.048.849.749.949.549.9ID58.858.646.950.146.947.046.843.448.146.345.245.246.1
NC011541 Narcissus yellow stripe virus49.248.348.749.349.849.261.9ID60.846.350.847.647.546.843.747.346.445.645.945.9
NC002509 Turnip mosaic virus49.248.649.051.049.949.962.063.0ID47.450.949.348.546.844.147.846.845.345.245.1
NC003605 Lettuce mosaic virus47.947.449.049.249.348.454.153.554.1ID47.745.344.143.942.446.143.742.743.042.9
NC001445 Plum pox virus59.050.051.051.251.159.052.951.752.551.4ID50.349.147.046.250.347.347.846.447.4
NC001841 Sweet potato feathery mottle virus48.948.954.056.049.149.751.451.451.154.054.2ID45.643.642.446.044.243.642.944.0
NC004039 Potato virus A48.449.348.549.449.148.558.051.351.149.652.255.8ID54.643.546.948.545.845.546.4
NC001768 Tobacco vein mottling virus49.649.447.848.353.049.249.149.148.748.949.849.652.1ID43.946.749.345.444.844.9
NC004011 Leek yellow stripe virus49.649.248.748.549.949.749.448.649.348.551.149.649.354.6ID44.043.247.645.045.7
NC004752 Yam mosaic virus48.847.949.049.549.149.057.049.454.049.352.051.551.348.948.5ID45.444.945.745.1
NC001555 Tobacco etch virus48.447.747.648.649.348.749.849.349.748.651.258.051.649.848.954.8ID44.845.245.6
NC040802 Lily yellow mosaic virus48.748.149.049.548.849.352.049.449.949.352.651.058.048.353.052.051.2ID49.349.4
NC003536 Clover yellow vein virus48.848.248.749.448.649.250.049.556.048.651.453.058.048.549.751.952.255.8ID68.9
NC003492 Bean yellow mosaic virus48.348.149.049.149.349.449.849.649.648.751.453.059.048.653.052.151.754.965.9ID
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Prajapati, M.R.; Manav, A.; Singhal, P.; Sidharthan, V.K.; Sirohi, U.; Kumar, M.; Bharti, M.K.; Singh, J.; Kumar, P.; Kumar, R.; et al. Complete Genomic RNA Sequence of Tuberose Mild Mosaic Virus and Tuberose Mild Mottle Virus Acquired by High-Throughput Sequencing. Pathogens 2022, 11, 861. https://doi.org/10.3390/pathogens11080861

AMA Style

Prajapati MR, Manav A, Singhal P, Sidharthan VK, Sirohi U, Kumar M, Bharti MK, Singh J, Kumar P, Kumar R, et al. Complete Genomic RNA Sequence of Tuberose Mild Mosaic Virus and Tuberose Mild Mottle Virus Acquired by High-Throughput Sequencing. Pathogens. 2022; 11(8):861. https://doi.org/10.3390/pathogens11080861

Chicago/Turabian Style

Prajapati, Malyaj R., Aakansha Manav, Pankhuri Singhal, Venkidusamy K. Sidharthan, Ujjwal Sirohi, Mukesh Kumar, Mahesh Kumar Bharti, Jitender Singh, Pankaj Kumar, Ravindra Kumar, and et al. 2022. "Complete Genomic RNA Sequence of Tuberose Mild Mosaic Virus and Tuberose Mild Mottle Virus Acquired by High-Throughput Sequencing" Pathogens 11, no. 8: 861. https://doi.org/10.3390/pathogens11080861

APA Style

Prajapati, M. R., Manav, A., Singhal, P., Sidharthan, V. K., Sirohi, U., Kumar, M., Bharti, M. K., Singh, J., Kumar, P., Kumar, R., Prakash, S., & Baranwal, V. K. (2022). Complete Genomic RNA Sequence of Tuberose Mild Mosaic Virus and Tuberose Mild Mottle Virus Acquired by High-Throughput Sequencing. Pathogens, 11(8), 861. https://doi.org/10.3390/pathogens11080861

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop