Multiplex detection, distribution, and genetic diversity of Hop stunt viroid and Citrus exocortis viroid infecting citrus in Taiwan
© Lin et al.; licensee BioMed Central. 2015
Received: 10 October 2014
Accepted: 20 January 2015
Published: 3 February 2015
Two citrus viroids, Citrus exocortis viroid (CEVd) and Hop stunt viroid (HSVd), have been reported and become potential threats to the citrus industry in Taiwan. The distributions and infection rates of two viroids have not been investigated since the two diseases were presented decades ago. The genetic diversities and evolutionary relationships of two viroids also remain unclear in the mix citrus planted region.
Multiplex RT-PCR was used to detect the two viroids for the first time in seven main cultivars of citrus. Multiplex real-time RT-PCR quantified the distributions of two viroids in four citrus tissues. Sequence alignment and phylogenetic analysis were performed using the ClustalW and MEGA6 (neighbor-joining with p-distance model), respectively.
HSVd was found more prevalent than CEVd (32.2% vs. 30.4%). Both CEVd and HSVd were commonly found simultaneously in the different citrus cultivars (up to 55%). Results of the multiplex quantitative analysis suggested that uneven distributions of both viroids with twig bark as the most appropriate material for studies involving viroid sampling such as quarantine inspection.
Sequence alignment against Taiwanese isolates, along with analysis of secondary structure, revealed the existence of 10 and 5 major mutation sites in CEVd and HSVd, respectively. The mutation sites in CEVd were located at both ends of terminal and variability domains, whereas those in HSVd were situated in left terminal and pathogenicity domains. A phylogenetic analysis incorporating worldwide viroid isolates indicated three and two clusters for the Taiwanese isolates of CEVd and HSVd, respectively.
Moderately high infection and co-infection rates of two viroids in certain citrus cultivars suggest that different citrus cultivars may play important roles in viroid infection and evolution. These data also demonstrate that two multiplex molecular detection methods developed in the present study provide powerful tools to understand the genetic diversities among viroid isolates and quantify viroids in citrus host. Our field survey can help clarify citrus-viroid relationships as well as develop proper prevention strategies.
KeywordsCitrus exocortis viroid (CEVd) Hop stunt viroid (HSVd) Multiplex RT-PCR Multiplex real-time RT-PCR Co-infection Uneven distribution
Viroids are small, circular, single-stranded noncoding RNAs that only infect plants. With tiny genome sizes (246–401 nt) and simple structures, viroids are the smallest known agents that infect hosts and cause disease. Because they lack gene-encoded proteins to provide specific functions, viroids depend on host-encoded factors and enzymes for replication [1-3]. Viroid replication occurs in specific subcellular compartments and trafficking throughout the plant, leading to complete systemic infection [4-6]. Viroids are classified into two families: Pospiviroidae and Avsunviroidae. Pospiviroidae species have a central conserved region (CCR) and do not contain hammerhead ribozymes, whereas the Avsunviroidae lack a CCR but can self-cleave through hammerhead ribozymes .
Citrus species are natural hosts of at least seven viroids in the family Pospiviroidae: Citrus exocortis viroid (CEVd, genus Pospiviroid), Citrus bent leaf viroid (CBLVd, CVd-I-b, genus Apscaviroid), Hop stunt viroid (HSVd, CVd-II, genus Hostuviroid), Citrus dwarfing viroid (CDVd, CVd-III, genus Apscaviroid), Citrus bark cracking viroid (CBCVd, CVd-IV, genus Cocadviroid), Citrus viroid V (CVd-V, genus Apscaviroid), and Citrus viroid VI (CVd-VI, genus Apscaviroid; initially named Citrus viroid original sample, CVd-OS) [8-11]. CEVd induces initial bark shelling and produces subsequent sloughing symptoms on trifoliate orange (Poncirus trifoliata [L.] Raf.), Troyer citrange, and Rangpur lime (Citrus × limonia Osb.), all widely used as rootstocks in commercial orchards [12,13]. HSVd variants with corresponding disease being known as cachexia induces discoloration, gumming, browning of phloem tissue, wood pitting, bark cracking, and stunting symptoms in mandarin (C. reticulata Blanco), clementine (C. clementina Hort. ex Tan.), satsuma (C. unshiu [Macf.] Marc.), alemow (C. macrophylla Webster), Rangpur lime, kumquat (Fortunella spp.), and mandarin hybrids such as tangelo (C. paradisi Macf. × C. tangerina Hosrt. ex Tan.) .
In Taiwan, citrus diseases caused by CEVd (365–475 nt) and HSVd (294–303 nt) were reported decades ago , but no further studies have been carried out during recent years. Following up on previous findings, in this study we simultaneously detected CEVd and HSVd isolates for the first time from seven citrus cultivars. We found that HSVd is more prevalent than CEVd and that co-infection with the two viroids is common in different citrus cultivars. In addition, we determined that twig bark is the best material for sampling viroids for assays using multiplex real-time RT-PCR . Furthermore, taking into account the diversity and complexity of citrus cultivars, we aligned sequences of 18 Taiwanese isolates of CEVd and HSVd with worldwide isolates and plotted their divergences against their type species to understand the diversity and evolution of the two main viroids in Taiwan, a complex mixed-citrus planted region.
Identification of two citrus viroids by three detection methods
Uneven distribution of viroids in different citrus tissues
CEVd and HSVd infection rates of citrus cultivars in Taiwan
Multiplex RT-PCR detection of CEVd and HSVd in samples collected from different regions in Taiwan
Locations (Number of samples)
Cultivars of rootstocks
Percentages of CEVd/HSVd infection rates (%/%)
Percentages of CEVd-infected/total samples (%)
Percentages of HSVd-infected/total samples (%)
Percentages of CEVd + HSVd co-infected/total samples (%)
Genetic diversity and secondary structure of CEVd and HSVd isolates
Nine samples each of CEVd and HSVd from various cultivars were cloned and sequenced using CEVd dF28/dR27 and HSVd F1/R1  primers. Eighteen sequences of CEVd and HSVd were deposited in GenBank under the following accession numbers: KC290927, KC290928, KJ956798-KJ956800, KJ956802-KJ956803, and KJ810543-KJ810544 for CEVd and KC290929, KJ956804-KJ956808, KJ810548, and KJ810552-KJ810553 for HSVd.
Current status of citrus diseases caused by infection with the two main viroids in Taiwan
Since the first reports of citrus diseases caused by CEVd or HSVd several decades ago, no studies have been performed on viroid diseases in Taiwan. The initial aim of our study was to fill in knowledge gaps and better understand distributions of these viroids in the citrus industry. Our field observations suggest that incorrect diagnoses are frequently made, as most viroid-infected citrus plants had latent infections without obvious external symptoms, except in sensitive rootstock 5 years old or younger. Based on our preliminary test, we tested four RNA extract methods as TRIzol buffer isolation, total nucleic acid extraction, TENS buffer extraction and purification by commercial kit (MACHEREY-NAGEL, GmbH & Co. KG, Düren, Germany). The TRIzol buffer isolation method is the best extraction method to purify viroids RNA (data not shown). When mass quantities of quarantined samples need to be assayed, multiplex RT-PCR and real-time RT-PCR are more efficient and specific for viroid detection than are field observations. In addition, the location of pathogens in infected hosts is the main issue when quarantine is required. Certain pathogens showed an uneven distribution in infected hosts; for example, Candidatus Liberibacter, which causes Huanglongbing in citrus, is limited to the phloem , and two beet viruses are distributed in different tissues of sugarbeet seedlings . Another previous study also indicated that the concentrations of several viroids as CEVd, HSVd or ASSVd may be varying and uneven when infect different hosts such as citrus, pear and hop . We estimated viroid distribution by real-time RT-PCR. This analysis revealed that both CEVd and HSVd were present in the lower parts of citrus plants, including roots and rootstocks where exocortis symptoms were found. Although roots and rootstock bark are the most abundant viroid-containing tissues, collection of these materials is time-consuming and difficult and the wounds may be contaminated with other pathogens. Taking into account sampling efficiency and convenience, we recommend using twig bark for viroid detection. In contrast, neither single-step RT-PCR nor qPCR could be successfully used to detect the two viroids present in high concentrations and copy numbers in leaf tissue. We found that viroids were more abundant in young tomato and citrus leaves than in mature ones (data not shown). Previous research has shown that members of different viroid genera, such as Potato spindle tuber viroid and Peach latent mosaic viroid, have opposite reactions when infecting shoot apical organs [20,21]. The cited studies demonstrate that mechanisms such as RNA silencing in leaves may operate to prevent viroid replication or transmission; nevertheless, additional evidence is required to support this hypothesis. Previous studies using multiplex RT-PCR for the field detection of citrus-infecting viroids have exploited a mixture of several primer sets in different concentrations [22,23]. Our method, which was focused on two main viroids, was developed to use a less complicated set of conditions. Our preliminary investigation yielded several discoveries. First, different cultivars and rootstocks may play an important role in viroid-host interactions and can increase citrus resistance or sensitivity, as they exhibited a range of disease scales upon viroid infection. For example, citrus plants grafted onto Rangpur lime rootstock displayed severe exocortis and stunting in Yunlin County, whereas pummelo in Hualien County showed no symptoms. Second, because mechanical inoculation is crucial to viroid transmission, viroid disease severity varied significantly according to crop management and field sanitation practices. In our experience, late-season field management was associated with the greatest severity of viroid disease. Third, most infected citrus trees showed latent infections and no symptoms unless the farmer used sensitive cultivars as rootstock. Another interesting finding was that simultaneous infection of individual citrus hosts with various viruses and viroids is common and severe in Taiwan. This situation is due to mixed crop plantings, where unknown interactions between each virus and viroid complicate detection of single diseases. Finally, slight differences in viroid sequences uncovered among citrus cultivars indicate that viroids use nucleotide and conformational alterations to adapt to different hosts. We found evidence of rapid mutation in viroids, with 2–5% divergence uncovered even within one host. Preliminary findings from our field survey were that nearly 35% of plants were infected with CEVd and HSVd viroids. The potential threat posed by viroid infection in citrus has been greatly underestimated. Regular field surveys are therefore necessary.
Genetic variation and diversity in Taiwanese CEVd and HSVd
Although CEVd and HSVd have existed in Taiwan for decades, no molecular analyses have been performed. In this study, we randomly selected nine isolates each of CEVd and HSVd from different citrus cultivars in various areas. One isolate each of the two viroids from Yunlin County (KC290927 and KC290929) were collected as Taiwanese type isolates for further analysis. According to the alignment analysis (data not shown), CEVd isolates from different cultivars had 93–99.4% identity, whereas HSVd isolates from various cultivars were 91.2-100% identical. Among analyzed CEVd isolates, the greatest sequence divergence was observed between Murcott orange and kumquat cultivars. Surprisingly, even when CEVd isolates were collected from an area in which both Murcott and blood oranges were growing, a high level of divergence was found between isolates. In contrast, HSVd isolates were most diverged between Tankan mandarin and Murcott orange. Interestingly, identities among HSVd isolates were lower on average than those among CEVd isolates, suggesting that the two viroids may experience different mutation rates. Previous studies have shown that viroids have higher mutation rates than other microbes . Even on a single citrus host, viroid populations are preserved as a genetic pool. High mutation rates may thus be triggered by different hosts or geographic factors, leading to more difficult, unreliable, and unstable detection conditions than for other pathogens. While the results of our sequence analysis might support this view, further investigations must be carried out, such as on viroid evolution in single hosts and viroid mutation among different citrus cultivars.
With respect to plotting of predicted secondary structures, we aligned all generated or downloaded CEVd and HSVd sequences against their Taiwanese type isolates. We identified 35 mutation sites, of which 10 were present more than three times among these isolates. The 10 major mutation sites were visualized on the plotted secondary structures, where they were distributed in both end terminals and variability domains as described in a previous study . Interestingly, most of the CEVd sequence variants exhibited five nucleotide substitutions or indels relative to the type sequence (KC290927): U19C, U27A, A32C, −74G, and G130A. In contrast, 20 mutation sites or substitutions were present among the HSVd sequences, but only five sites were repeated more than three times among them. Five sites were located on left terminal and pathogenicity domains of the secondary structure. Almost all of the HSVd sequence variants featured three unique nucleotide substitutions or indels compared with the type sequence (KC290929), namely G26AA, G53-, and U246G. The presence of these mutations, especially those in the pathogenicity domain, is interesting, as their occurrence may be responsible for the fact that HSVd Taiwanese isolates cause scarcely any symptoms in certain hosts. Additional experiments are required to confirm which mutated position is able to affect symptoms of HSVd infection. Finally, the insertion -288GGATCC was observed in four of the aligned HSVd isolates; this additional sequence was derived from the BamHI restriction enzyme cleavage site and is an experimental artifact of the cloning procedure.
In the phylogenetic analysis, nine CEVd isolates and eight HSVd isolates were analyzed along with selected worldwide isolates or variants. The CEVd Taiwanese isolates were distributed among three clusters, whereas the HSVd Taiwanese isolates fell into two separate clusters. Of the three CEVd clusters, one was most similar to Canadian (DQ094296) and Australian (S67437) strains, one was most closely related to a USA strain (GU295988) and the third was similar to a Uruguayan strain (AF148717). Interestingly, two isolates collected from kumquat in Yilan County were placed in separate clusters, suggesting that CEVd isolates may be quite different even when present in the same host. In the phylogenetic tree of HSVd sequences, one cluster was identical to a Chinese strain (EU382208) and the other to a Colombian strain (EU519454). A previous study found that AGVd isolates from India were more closely related to Chinese strains, suggesting that geographic proximity is an important factor shaping genetic similarity  and another studies also addressed the species-dependency of the population structure of six viroids while infecting grapevines in China and Japan . In our study, however, no relationship between geographical region and genetic variation was apparent. To clarify the evolutionary origin of these viroids, additional sequences of worldwide isolates should be aligned and analyzed.
We used three methods to detect two major viroids infecting citrus plants in Taiwan. Multiplex real-time RT-PCR revealed that the two viroids had similarly uneven distributions in different citrus tissues and suggested that twig bark may be the best material for sampling. Our multiplex RT-PCR-based field survey also gave preliminary estimates of CEVd and HSVd infection of nearly 30-35% in all citrus orchards in Taiwan. Sequence alignment against Taiwanese type isolates and secondary structure analysis revealed that 10 mutation sites were distributed on both ends of CEVd terminal and variability domains while five mutation sites existed in HSVd left terminal and pathogenicity domains. Finally, phylogenetic analyses incorporating worldwide viroid isolates indicated that Taiwanese isolates of CEVd are separated into three clusters and HSVd isolates form two clusters. Citrus diseases caused by viroids are a rising threat in the citrus industry in Taiwan. The results of this study should provide a solid basis for exploration of the two citrus viroids and for understanding the relationship between viroids and citrus plants. The fact that nearly 35% of citrus plants in Taiwan are infected by a combination of CEVd and HSVd suggests that viroid diseases will become an increasingly important and urgent global issue. Further studies must be performed in the near future and should focus on the ecology and epidemiology of viroids.
Materials and methods
Plant materials and viroid sources
The three repeat indicator plants Etrog citron Arizona 861-S and G. aurantiaca for CEVd were planted for stem cutting and grafting, respectively; the indicator plant C. sativus for HSVd was grown by germination from seeds (KNOWN-YOU Seed Co., Taiwan). For mechanical inoculation, infected citrus leaves from the field were ground in a 0.5 M K-P buffer (0.4 M K2HPO4, 0.08 M KH2PO4, 0.01 M EDTA-2Na, 1 mM NA2SO4, pH 7.5) as inoculum sap with carborundum was spotted onto three mature leaves of each plant by use of cotton swabs. For graft inoculation, plants with at least two infected buds, shoots and petioles were grafted with one healthy 861-S plant. Inoculated and healthy plants were grown in a greenhouse with 16 h light (28°C) and 8 h dark (24°C) cycles. Plants were regularly watered with commercial plant nutrients and checked for symptoms.
TRIzol RNA isolation method
This method was described by  with some modification. Citrus tissue (100 mg) was homogenized using pestle and liquid nitrogen in 1 mL TRIzol reagent (0.8 M guanidine thiocyanate; 0.4 M ammonium thiocyanate; 0.1 M sodium acetate, at pH 5.0; 5% glycerol; and 38% phenol in a saturated buffer). After centrifugation (13,500 × g) for 10 min, the supernatant was added a 1/5 volume of chloroform, shaken vigorously for 15 s and placed at room temperature for 3 min. The supernatant was collected and added a 1/2 volume of isopropanol and 1/2 volume of 0.8 M sodium citrate/1.2 M NaCl and then placed at room temperature for 10 min. The RNA pellet was collected by high-speed centrifugation (17,000 × g) at 4°C for 15 min and then washed with 70% ethanol and resuspended in 50 μL DEPC-water.
Sequences of primer pairs for sequencing and multiplex detection methods
Details of the primers used for the molecular detections and quantifications of CEVd and HSVd
Primer pairs for multiplex RT-PCR
Product size (nt)
Full-length cDNA for sequencing and qPCR
Primers/Probes for multiplex real-time RT-PCR
Product size (nt)
Amplification and sequencing of viroid cDNA
Complementary DNAs were synthesized by M-MLV reverse transcriptase enzyme (Invitrogen, Carlsbad, CA, USA) using primers as described previously according to the manufacturer’s description. The fragment of viroid cDNA was then amplified by polymerase chain reaction (PCR) using Taq polymerase (Gibco BRL) in each primer pairs as described previously. The error rate of Taq polymerase is 2.2 × 10−5. Amplification products were analyzed on 1.6% agarose gel electrophoresis. The DNA fragments were cloned into pCR2.1 TOPO (Invitrogen, Carlsbad, CA, USA) and 5–10 clones of each fragment were sequenced with the use of the ABI PRISM automated DNA sequencer (Applied Biosystems). The nucleotide sequences were analyzed and blasted against each other by the Clustal W method using MEGA 6 software .
Multiplex detection methods
Multiplex RT-PCR was performed in a single tube with 25 μL of the reaction mixture containing 45 mM Tris–HCl (pH 8.3), 80 mM KCl, 4 mM MgCl2, 0.2 mM dNTP, 5 mM 1, 4-dithiothreitol (DTT), 5 pmol of CEVd and HSVd primer pairs as mentioned above, 50 units of reverse transcriptase (Superscript II RNase H Reverse Transcriptase), 25 units of Taq DNA polymerase (Gibco BRL), and a 300 ng template of the nucleic acid. The thermal cycle conditions were one cycle at 50°C for 35 min; one cycle at 94°C for 2 min; 40 cycles at 94°C for 30 s, 58°C for 30 s, 72°C for 60 s; and then 72°C extension for 8 min. Reactions were carried out in a DNA Thermal Cycler 2720 (Applied Biosystems).
Multiplex real-time RT-PCR
Primers and probes
All primers and probes of the two viroids are listed in Table 2. Conserved regions of probes were identified and sent to Applied Biosystems Co. (Foster City, USA) for primer and probe designations. The 5′ terminal reporter dyes were FAM (6-carboxyfluorescein) for CEVd probe; VIC for HSVd probe, and the 3′ quencher dyes were NFQ (non-fluorescent quencher) and MGB (minor groove binder).
Standard curve and TaqMan assay setup
The standard curves for absolute quantification of CEVd and HSVd were established as described by the manufacturer (Applied Biosystems). The final concentration ranges of standard curves of CEVd and HSVd were measured from 1010 to 10 copies by serial dilution. The full-length of CEVd and HSVd were amplified using the CEVd dF28/dR27 and HSVd F1/R1 (Table 2) primer pairs and cloned into the pcr2.1-TOPO vector (Invitrogen, Life technologies, Carlsbad, CA, USA), respectively. The measurement of plasmid concentration was recorded as the absorbance at 260 nm and the copy number was calculated using the plasmid’s molecular weight and Avogadro’s constant. Synthesis of first strand cDNA for real-time PCR was followed by methods described in the above mentioned. The TaqMan PCR reactions (total volume 20 μL) were set up in 48-well reaction plates using PCR Master Mix reagent kits (Applied Biosystems). The real-time PCR conditions followed Applied Biosystems’ instruction with slight modification: 9 μL (100 ng) of cDNA, 10 μL 2X TaqMan® Gene Expression Master Mix, 1 μL 20X assay mix for each CEVd and HSVd (900 nM forward primer, 900 nM reverse primer, 250 nM TaqMan® MGB probes), and ROX passive reference dye (concentration undisclosed by manufacturer). The assays were carried out on an ABI StepOne Real-Time PCR System (Applied Biosystems) at cycling conditions (50°C /2 min, 95°C /10 min and 40 cycles of 95°C /15 min, 60°C /1 min).
Secondary structure prediction and phylogenetic analysis
All the sequences of CEVd or HSVd were aligned with each type isoaltes (GenBank Acc. No. KC290927 and KC290929). Possible secondary structures of two viroids were predicted using the software tool as Mfold (http://mfold.rna.albany.edu/?q=mfold) and sequence variants in Taiwan were plotted against each secondary structure. All sequences were aligned with the others of randomly selected CEVd or HSVd sequences which submitted in the GenBank database using the ClustalW method and phylogenetic analysis were performed using neighbor-joining with p-distance model of MEGA6 software .
We appreciate the grants received from the Council of Agriculture (nos. 100AS-1.1.1-BQ-B1, 100AS-9.1.1-BQ-B2) and the National Science Council (no. 98-2313-B-002-040-MY3), Executive Yuan, Taiwan. We would also like to thank Dr. Hsin-Hung Yeh and Dr. Chin-Cheng Yang for generously revising the manuscript and providing many valuable opinions; Dr. Cheng-En Chen for kindly sharing experiences and skills and the farmers from Taiwan for providing citrus cultivars and samples to complete the field survey.
- Diener TO, Maramorosch K, Murphy FA, Shatkin AJ. The viroid: biological oddity or evolutionary fossil? Adv Virus Res. 2001;57:137–84.View ArticlePubMedGoogle Scholar
- Flores R, Hernandez C, de Alba AE M, Daròs JA, Di Serio F. Viroids and viroid-host interactions. Annu Rev Phytopathol. 2005;43:117–39.View ArticlePubMedGoogle Scholar
- Ding B, Itaya A, Zhong X. Viroid trafficking: a small RNA makes a big move. Curr Opin Plant Biol. 2005;8:606–12.View ArticlePubMedGoogle Scholar
- Ding B, Itaya A. Control of directional macromolecular trafficking across specific cellular boundaries: a key to integrative plant biology. J Integr Plant Biol. 2007;49:1227–34.View ArticleGoogle Scholar
- Ding B, Itaya A. Viroid: a useful model for studying the basic principles of infection and RNA biology. MPMI. 2007;20:7–20.View ArticlePubMedGoogle Scholar
- Ding B. The biology of viroid-host interactions. Annu Rev Phytopathol. 2009;47:105–31.View ArticlePubMedGoogle Scholar
- Flores R, Delgado S, Gas ME, Carbonell A, Molina D, Gago S, et al. Viroids: the minimal non-coding RNAs with autonomous replication. FEBS Lett. 2004;567:42–8.View ArticlePubMedGoogle Scholar
- Duran-Vila N, Roistacher CN, Rivera-Bustamante R, Semancik JS. A definition of citrus viroid groups and their relationship to the exocortis disease. J Gen Virol. 1988;69:3069–80.View ArticleGoogle Scholar
- Foissac X, Duran-Vila N. Characterisation of two citrus apscaviroids isolated in Spain. Arch Virol. 2000;145:1975–83.View ArticlePubMedGoogle Scholar
- Ito T, Ieki H, Ozaki K, Ito T. Characterization of a new citrus viroid species tentatively termed citrus viroid OS. Arch Virol. 2001;146:975–82.View ArticlePubMedGoogle Scholar
- Serra P, Barbosa CJ, Daròs JA, Flores R, Duran-Vila N. Citrus viroid V: molecular characterization and synergistic interactions with other members of the genus Apscaviroid. Virol. 2008;370:102–12.View ArticleGoogle Scholar
- Duran-Vila N, Flores R, Semancik JS. Characterization of viroidlike RNAs associated with the citrus exocortis syndrome. Virol. 1986;150:75–84.View ArticleGoogle Scholar
- Fawcett HS, Klotz LJ. Exocortis of trifoliate orange. Citrus Leaves. 1948;28:8.Google Scholar
- Hsu YH, Chen W, Owens RA. Nucleotide sequence of a hop stunt viroid variant isolated from citrus growing in Taiwan. Virus Genes. 1995;9:193–5.View ArticlePubMedGoogle Scholar
- Papayiannis LC, Papayiannis LC. Diagnostic real-time RT-PCR for the simultaneous detection of Citrus exocortis viroid and Hop stunt viroid. J Virol Methods. 2014;196:93–9.View ArticlePubMedGoogle Scholar
- Bernad L, Duran-Vila N. A novel RT-PCR approach for detection and characterization of citrus viroids. Mol Cellul Probes. 2006;20:105–13.View ArticleGoogle Scholar
- Jagoueix S, Bove JM, Garnier M. The phloem-limited bacterium of greening disease of citrus is a member of the a subdivision of the Proteobacteria. Intl J System Bacterial. 1994;44:379–86.View ArticleGoogle Scholar
- Kaufmann A, Koenig R, Lesemann DE. Tissue print-immunoblotting reveals an uneven distribution of beet necrotic yellow vein and beet soil-borne viruses in sugarbeets. Arch Virol. 1992;126:329–35.View ArticlePubMedGoogle Scholar
- Li SF, Onodera S, Sano T, Yoshida K, Wang GP, Shikata E. Gene diagnosis of viroids: comparisons of return-PAGE and hybridization using DIG-labeled DNA and RNA probes for practical diagnosis of Hop stunt, Citrus exocortis and Apple scar skin viroids in their natural host plants. Ann Phytopathol Soc Jpn. 1995;61:381–90.View ArticleGoogle Scholar
- Zhu Y, Green L, Woo YM, Owens RA, Ding B. Cellular basis of Potato spindle tuber viroid systemic movement. Virol. 2001;279:69–77.View ArticleGoogle Scholar
- Rodio ME, Delgado S, De Stradis A, Gòmez MD, Flores R, Di Serio F. A viroid RNA with a specific structural motif inhibits chloroplast development. Plant Cell. 2007;19:3610–26.View ArticlePubMed CentralPubMedGoogle Scholar
- Ito T, Namba N, Ito T. Distribution of citrus viroids and Apple stem grooving virus on citrus trees in Japan using multiplex reverse transcription polymerase chain reaction. J Gen Plant Pathol. 2003;69:205–7.View ArticleGoogle Scholar
- Wang XF, Zhou CY, Tang KZ, Zhou Y, Li ZG. A rapid one-step multiplex RT-PCR assay for the simultaneous detection of five citrus viroids in China. Eur J Plant Pathol. 2009;124:175–80.View ArticleGoogle Scholar
- Gago S, Elena SF, Flores R, Sanjuan R. Extremely high mutation rate of a hammerhead viroid. Science. 2009;323:1308.View ArticlePubMedGoogle Scholar
- Keese P, Symons RH. Domains in viroids: evidence of intermolecular RNA rearrangements and their contribution to viroid evolution. Proc Natl Acad Sci U S A. 1985;82:4582–6.View ArticlePubMed CentralPubMedGoogle Scholar
- Adkar-Purushothama CR, Kanchepalli PR, Yanjarappa SM, Zhang Z, Sano T. Detection, distribution, and genetic diversity of Australian grapevine viroid in grapevines in India. Virus Genes. 2014;49:304–11.View ArticlePubMedGoogle Scholar
- Jiang D, Sano T, Tsuji M, Araki H, Sagawa K, Purushothama CRA, et al. Comprehensive diversity analysis of viroids infecting grapevine in China and Japan. Virus Res. 2012;169:237–45.View ArticlePubMedGoogle Scholar
- Chomczynski P, Sacchi N. Single step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. Anal Biochem. 1987;162:156–9.View ArticlePubMedGoogle Scholar
- Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.View ArticlePubMed CentralPubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.