Rapid Identification of Jasmine Virus H Infecting Ixora coccinea by Nanopore Metatranscriptomics

Article information

Plant Pathol J. 2023;39(3):303-308
1Department of Applied Biology, Institute of Environmentally Friendly Agriculture, Chonnam National University, Gwangju 61185, Korea
2Horticultural and Herbal Crop Environment Division, National Institute of Horticultural and Herbal Science, Rural Development Administration, Wanju 55365, Korea
*Corresponding author. Phone) +82-62-530-2075, FAX) +82-62-530-2069, E-mail) jraed2@jnu.ac.kr
Handling Editor : Ju-Yeon Yoon
Received 2023 March 6; Revised 2023 April 10; Accepted 2023 April 23.

Abstract

The global climate change and international trade have facilitated the movement of plants across borders, increasing the risk of introducing novel plant viruses in new territories. Ixora coccinea exhibited virus-like foliar symptoms, including mosaic and mild mottle. An Oxford Nanopore Technologies-based compact and portable MinION platform was used to identify the causal viral pathogen. The complete genome sequence of jasmine virus H (JaVH; 3867 nt, JaVH-CNU) was determined and found to share 88.4–90.3% nucleotide identity with that of Jasminum sambac JaVH isolate in China. Phylogenetic analysis based on the complete amino acid sequences of RNA-dependent RNA polymerase and coat protein revealed that JaVH-CNU was grouped separately with other JaVH isolates. This is the first report of a natural JaVH infection of I. coccinea. The application of rapid nanopore sequencing for plant virus identification was demonstrated and is expected to provide accurate and rapid diagnosis for virus surveillance.

Ixora coccinea, commonly known as jungle geranium, is a flowering shrub native to the countries of Eastern South Asia, such as India, Bangladesh, and Sri Lanka, and it is now cultivated in both tropical and subtropical countries. In addition, it is used as a medicinal plant for skin disease, dysentery, ulcers, and gonorrhea (Shreelakshmi et al., 2021). With the increasing demand, the import and cultivation of I. coccinea has been growing in South Korea, and this is mostly propagated by cuttings (Lin et al., 2020). Although several studies reported natural infections of fungal pathogens, including Mycoleptodiscus indicus and Colletotrichum aeschynomenes, in I. coccinea (Banerjee et al., 2018; Li et al., 2021), no viral pathogens have been reported to date.

Recently, Illumina high-throughput sequencing-based diagnostic platforms have been widely utilized to identify various unknown viral pathogens in crops (Ibaba and Gubba, 2020; Rivarez et al., 2021). In addition to its advantages, this technology also has some limitations in viral diagnosis, including costly equipment, long sequencing times, low portability, and complex data analysis (Fitzpatrick et al., 2021). Nanopore sequencing introduced by Oxford Nanopore Technology (ONT), known as third-generation sequencing, has several advantages compared to Illumina-based sequencing, including reduced cost, high portability, single-molecule (DNA or RNA) sequencing, long sequencing reads, faster sequencing, and easy data analysis that can be performed on a personal computer, which facilitate its use in plant viral diagnostics (Lee et al., 2022; Sun et al., 2022). The biggest drawbacks of nanopre sequencing have been the lower throughput of sequence data and an approximate 10% high error rate. However, recent research has achieved validation of the accuracy of nanopore sequencing long reads through improvements in technology (Lu et al., 2016).

In September 2022, two out of the 20 I. coccinea plants showing foliar disease symptoms with mosaic and mild mottle were collected from the streets of Jeju Island, South Korea, where they had been displayed (Fig. 1). Total RNA was extracted from symptomatic leaf samples (~200 mg) using a Clear-S Total RNA extraction kit (Invirustech Co., Gwangju, Korea) according to the manufacturer’s instructions. DNA removal and RNA purification for library construction were performed as described by Lee et al. (2022). The purified RNA was used for cDNA synthesis using a Maxima H Minus Double-Stranded cDNA Synthesis Kit (Thermo Fisher Scientific, Waltham, MA, USA). For library construction, dA tailing and adapter ligation were performed using a NEBNext Ultra II End Repair/dA-Tailing Module (New England Biolabs, Ipswich, MA, USA) and a Blunt/TA Ligase Master Mix (New England Biolabs) with an Adapter Mix (SQK-DCS109, Oxford Nanopore Technologies, Oxford, UK), respectively, generating the nanopore library with a final concentration of 3.3 ng/μl. The procedures of sample loading into Flongle flow cell (R9.4.1), sequencing of the nanopore library, and subsequent analyses of the reads were performed as described by Lee et al. (2022). For metatranscriptomics analysis, the FASTQ files were converted to FASTA files, and the adapter sequences were removed using the tool Geneious Prime (version 21.1.1, Biomatters, Auckland, New Zealand). The FASTQ data files were deposited in the National Center for Biotechnology Information (NCBI) Sequence Read Archive repository under SRR23561814. To identify genome sequences of plant viruses, ONT EPI2ME WIMP workflow (version 3.4.2) and BLASTn search on the NCBI Viral Genome data (version 21.11.4, Viral RefSeq Database, VirDB) were utilized as described by Lee et al. (2022).

Fig. 1

Ixora coccinea leaf with mosaic and mild mottle symptoms.

The nanopore sequencing generated a total of 434,106 raw reads with an average length of 317.1 nucleotides. After adapter removal and quality filtering steps (>Q7), a total of 349,069 reads with an average length of 44.6 nucleotides were obtained (Table 1). The identification of plant virus-associated sequences was performed using the EPI2ME WIMP workflow. A total of 796 (0.22%) virus-like reads were found and were BLASTed on the NCBI VirDB. All virus-like reads were mapped to jasmine virus H (JaVH) (Table 2).

Summary of Oxford Nanopore sequencing data

List of identified JaVH reads in Ixora coccinea sample

De novo assembly of MinION reads was performed with genome sequences of JaVH isolate (JaVH-CNU) (Fig. 2A). A total of 796 reads (0.22% of the total reads) were mapped with 95.5% nucleotide coverage and 95.1% nucleotide identity. To obtain the complete genome sequence of JaVH, the consensus sequences were finally mapped against a reference viral genome (NC055545; at 60× coverage depth and at least 20% of support fraction for base-call ambiguity). In addition, Sanger sequencing after reverse transcription polymerase chain reaction (RT-PCR) was performed for filling in any sequence gaps using virus-specific primers (Fwd: 5′-TCC AAG GCG AAT GCT CTC TCT GT-3′/Rev: 5′-GGA TTG TAG TGG AGC GGT GAA AAA CCC-3′), which amplify the fragments between the aligned reads, designed based on the consensus sequence. The RT-PCR reaction mixture (final volume 20 μl) contained 10 μl SuPrimerScript RT Premix (Genet Bio, Daejeon, Korea), 2 μl RNA template, 2 μl of each of the 10 μM forward and reverse primers, and diethylpyrocarbonate-treated water. Amplification reaction was conducted under the following conditions: 50°C for 30 min; 95°C for 5 min; 35 cycles of 95°C for 30 s, 62°C for 30 s, and 72°C for 30 s; and 72°C for 5 min. The 5′ and 3′ end sequences of JaVH were determined by rapid amplification of cDNA end (RACE) technique using a SMARTer RACE 5′/3′ Kit (Clontech, Mountain View, CA, USA). The complete genome consensus sequence of the JaVH infecting I. coccinea was deposited in GenBank (LC757512).

Fig. 2

Genome assembly of viral genome sequences determined by nanopore sequencing. (A) Genomic mapping of JaVH-CNU identified in Ixora coccinea. The blue graph shows the distribution of reads, and each read is represented with a black bar. (B) Genome structure of JaVH-CNU isolate from I. coccinea. (C) Confirmation of presence of JaVH in I. coccinea by reverse transcription polymerase chain reaction. M, 100-bp ladder; lane 1, JaVH-infected I. coccinea sample; lane 2, healthy sample. JaVH, jasmine virus H.

The genome structure of JaVH-CNU was predicted using the ORFfinder software (http://www.ncbi.nlm.nih.gov/orffinder). The RNA genome of JaVH-CNU consisted of 3,867 nucleotides (nt), with a total of four open reading frames (ORFs) (Fig. 2B). The 5′ and 3′ untranslated regions were found to be 19 nt and 226 nt in length, respectively. The first ORF (ORF1) comprised 729 nt and has potential to read-through a downstream ORF, resulting in production of a 2,301 nt sequence encoding the RNA-dependent RNA polymerase (RdRP). The ORF2 (189 nt) and ORF3 (252 nt) encoded movement protein. The ORF4 (1020 nt) encoded a coat protein (CP) (Fig. 2B).

To confirm the presence of JaVH in I. coccinea, RT-PCR was performed using JaVH-specific diagnostic primer sets (Dey et al., 2018), and the RT-PCR product (657 bp) was sequenced by Sanger sequencing (Fig. 2C). The obtained sequences confirmed the presence of JaVH (data not shown).

The complete genome of JaVH-CNU exhibited sequence identities ranging from 88.4% (JaVH-Hunan, MH231179) to 90.3% (JaVH-Fujian, NC055545) with other JaVH isolates. In addition, to compare the amino acid sequences, the complete amino acid sequences of RdRP and CP of JaVH-CNU were aligned with those of other JaVH isolates available on GenBank using BioEdit software version 7.2.5 (Ibis Biosciences, Carlsbad, CA, USA). Pairwise distances were calculated using the PASC algorithm (NCBI, Bethesda, MD, USA) available on the GenBank database. JaVH-CNU shared amino acid sequence identities of 94.9–95.3% in RdRP (Fig. 3A) and 94.4–96.2% in CP (Fig. 3B).

Fig. 3

Pairwise comparisons of the complete amino acid sequences of RdRP (A) and CP genes (B) of JaVH-CNU. The accession numbers, geographic locations, and host plants of each JaVH isolate are listed, ordered by the taxon name. Each number represents pairwise identities between the corresponding isolates. JaVH, jasmine virus H; RdRP, RNA-dependent RNA polymerase; CP, coat protein.

To assess the phylogenetic relationship among JaVH-CNU and other JaVH isolates, the complete amino acid sequences of the RdRP and CP of JaVH-CNU, along with those of other JaVH isolates, were used in phylogenetic analysis. The phylogenetic tree was constructed using the MEGA 10.0 tool (Kumar et al., 2018). The alignments were used to infer neighbor-joining trees in MEGA 10.0 with Tamura-Nei model and 1,000 bootstrap replicates as described by Lee et al. (2022). Phylogenetic analysis of RdRP and CP revealed that the JaVH-CNU differed from the existing isolates in GenBank. The RdRP and CP of JaVH-CNU were clustered separately with different JaVH isolates, probably having different biological characteristics (Fig. 4).

Fig. 4

Phylogenetic analyses using the neighbor-joining method based on the complete amino acid sequences of RdRP (A) and CP (B). The reference sequences are indicated with their accession numbers and names. The genome of the JaVH-CNU is marked with black circles. The scale bar indicates genetic distance. JaVH, jasmine virus H; RdRP, RNA-dependent RNA polymerase; CP, coat protein.

The spread of plan viruses through global trade is a significant concern. The international movement of infected plant materials or their products can introduce viruses to new geographic locations and host plant species, potentially leading to the emergence of new disease. Moreover, the movement of insect vectors, such as aphids or whiteflies, can also facilitate the spread of plant viruses through global (Amari et al., 2021). To successfully respond to an unexpected outbreak of plant viruses, it has become increasingly important to establish rapid and reliable diagnostic methods for effective disease management.

With the advantages of nanopore sequencing, this study sequenced and identified the full genome of a JaVH isolate from I. coccinea, and it took only 48 h from RNA extraction to virus identification in the laboratory, providing rapid diagnosis for disease survey and management. The JaVH, which has been recently characterized from Jasminum sambac (Zhuo et al., 2017) and J. multiflorum (Dey et al., 2018), belongs to a member of the genus Pelarspovirus in the family Tombusviridae. Despite the JaVH isolates reported in other countries and hosts, the genome of JaVH seems to be highly variable based on phylogenetic analyses of complete amino acid sequences of RdRP and CP (Dey et al., 2018). It is worth noting that this is first report of a natural infection of JaVH in I. coccinea worldwide. Thus, it is urgently required to survey the incidence of JaVH in I. coccinea in Korea, and more importantly, to evaluate the disease risk in different types of crops and ecologies.

Notes

Conflicts of Interest

No potential conflict of interest relevant to this article was reported.

Acknowledgments

This work was carried out with the support of Cooperative Research Program for Agriculture Science and Technology Development (Project No. PJ014947032023) Rural Development Administration, Republic of Korea.

References

Amari K., Huang C., Heinlein M.. 2021;Potential impact of global warming on virus propagation in infected plants and agricultural productivity. Front. Plant Sci 12:649768.
Banerjee A., Mandal R., Nath P. S. 2018;First report of leaf blight and twig dieback caused Mycoleptodiscus indicus on Ixora coccinea. New Dis. Rep 38:5.
Dey K. K., Leite M., Hu J. S., Jordan R., Melzer M. J. 2018;Detection of Jasmine virus H and characterization of a second pelarspovirs infecting star jasmine (Jasminum multiflorum) and angelwing jasmine (J. nitidum) plants displaying virus-like symptoms. Arch. Virol 163:3051–3058.
Fitzpatrick A. H., Rupnik A., O’Shea H., Crispie F., Keaveney S., Cotter P.. 2021;High throughput sequencing for the detection and characterization of RNA viruses. Front. Microbiol 12:621719.
Ibaba J. D, Gubba A.. 2020;High-throughput sequencing application in the diagnosis and discovery of plant-infecting viruses in Africa, a decade later. Plants 9:1376.
Kumar S., Stecher G., Li M., Knyaz C., Tamura K.. 2018;Mega X: molecular evloutionalry genetics analysis across computing platforms. Mol. Biol. Evol 35:1547–1549.
Lee Y.-J., Cho I.-S, Jeong R.-D. 2022;Nanopore metagenomics sequencing for rapid diagnosis and characterization of lily viruses. Plant Pathol. J 38:503–512.
Li M., Zhang S., Gao Z., Wang Y., Zhang W., Gong D., Zhao C., Hu M.. 2021;First report of Colletotrichum aeschynomenes causing leaf anthracnose on Ixora coccinea in China. J. Plant Pathol 103:1337–1338.
Lin Y.-R., Lee S., Lu C.-H, Chu C.-C. 2020;Genetic phenotypic characterization of Xanthomonas axonopodis pv. maculifoliigardeniae causing bacteral leaf spot of Ixora in Taiwan. J. Phytopathol 168:478–489.
Lu H., Giordano F., Ning Z.. 2016;Oxford nanopore MinION sequencing and genome assembly. Genomics Proteomics Bioinformatics 14:265–279.
Rivarez M. P. S., Vučurovič A., Mehle N., Ranvikar M., Kutnjak D.. 2021;Global advances in tomato virome research: current status and the impact of high-throughput sequencing. Front. Microbiol 12:671925.
Shreelakshmi S. V., Chaitrashree N., Kumar S. S., Shetty N. P, Giridhar P.. 2021;Fruits of Ixora coccinea are a rich source of phytoconstituents, bioactives, exhibit antioxidant activity and cytotoxicity against human prostate carcinoma cells and development of RTS beverage. J. Food Process. Preserv 45:e15656.
Sun K., Liu Y., Zhou X., Yin C., Zhang P., Yang Q., Mao L., Shentu X., Yu X.. 2022;Nanopore sequencing technology and its application in plant virus diagnostics. Front. Microbiol 13:939666.
Zhuo T., Zhu L.-J., Lu C.-C., Jiang C.-Y., Chen Z.-Y., Zhang G., Wang Z.-H., Jovel J., Han Y.-H. 2017;Complete nucleotide sequence of jasmine virus H, a new member of the family Tombusviridae. Arch. Virol 163:731–735.

Article information Continued

Fig. 1

Ixora coccinea leaf with mosaic and mild mottle symptoms.

Fig. 2

Genome assembly of viral genome sequences determined by nanopore sequencing. (A) Genomic mapping of JaVH-CNU identified in Ixora coccinea. The blue graph shows the distribution of reads, and each read is represented with a black bar. (B) Genome structure of JaVH-CNU isolate from I. coccinea. (C) Confirmation of presence of JaVH in I. coccinea by reverse transcription polymerase chain reaction. M, 100-bp ladder; lane 1, JaVH-infected I. coccinea sample; lane 2, healthy sample. JaVH, jasmine virus H.

Fig. 3

Pairwise comparisons of the complete amino acid sequences of RdRP (A) and CP genes (B) of JaVH-CNU. The accession numbers, geographic locations, and host plants of each JaVH isolate are listed, ordered by the taxon name. Each number represents pairwise identities between the corresponding isolates. JaVH, jasmine virus H; RdRP, RNA-dependent RNA polymerase; CP, coat protein.

Fig. 4

Phylogenetic analyses using the neighbor-joining method based on the complete amino acid sequences of RdRP (A) and CP (B). The reference sequences are indicated with their accession numbers and names. The genome of the JaVH-CNU is marked with black circles. The scale bar indicates genetic distance. JaVH, jasmine virus H; RdRP, RNA-dependent RNA polymerase; CP, coat protein.

Table 1

Summary of Oxford Nanopore sequencing data

Sample ID Raw reads Trimmed reads SRA accession no.


Total reads Minimum read length (nt) Maximum read length (nt) Mean read length (nt) Total reads Minimum read length (nt) Maximum read length (nt) Mean read length (nt)
Iksora cococinea 434,106 89 6,990 317.1 349,069 20 903 44.6 SRR23561814

SRA, Sequence Read Archive.

Table 2

List of identified JaVH reads in Ixora coccinea sample

Virus name Read no. Maximum read length (nt) Mean read length (nt) Identity (%) Coverage (%) Accession no.
Jasmine virus H (JaVH) 796 51 11.8 95.1 95.5 LC757512