- Research article
- Open Access
Comparative genome analysis of jujube witches’-broom Phytoplasma, an obligate pathogen that causes jujube witches’-broom disease
BMC Genomics volume 19, Article number: 689 (2018)
JWB phytoplasma is a kind of insect-transmitted and uncultivable bacterial plant pathogen causeing a destructive Jujube disease. To date, no genome information about JWB phytoplasma has been published, which hindered its characterization at genomic level. To understand its pathogenicity and ecology, the genome of a JWB phytoplasma isolate jwb-nky was sequenced and compared with other phytoplasmas enabled us to explore the mechanisms of genomic rearrangement.
The complete genome sequence of JWB phytoplasma (jwb-nky) was determined, which consisting of one circular chromosome of 750,803 bp with a GC content of 23.3%. 694 protein-encoding genes, 2 operons for rRNA genes and 31 tRNA genes as well as 4 potential mobile units (PMUs) containing clusters of DNA repeats were identified. Based on PHIbaes analysis, a large number of genes were genome-specific and approximately 13% of JWB phytoplasma genes were predicted to be associated with virulence. Although transporters for maltose, dipeptides/oligopeptides, spermidine/putrescine, cobalt, Mn/Zn and methionine were identified, KEGG pathway analysis revealed the reduced metabolic capabilities of JWB phytoplasma. Comparative genome analyses between JWB phytoplasma and other phytoplasmas shows the occurrence of large-scale gene rearrangements. The low synteny with other phytoplasmas indicated that the expansion of multiple gene families/duplication probably occurred separately after differentiation.
In this study, the complete genome sequence of a JWB phytoplasma isolate jwb-nky that causing JWB disease was reported for the first time and a number of species-specific genes were identified in the genome. The study enhanced our understandings about genomic basis and the pathogenicity mechanism of this pathogen, which will aid in the development of improved strategies for efficient management of JWB diseases.
Pytoplasma was a group of plant pathogenic bacteria that belong to the class of Mollicutes. Phylogenetic analysis, virtual restriction fragment length polymorphism (RFLP) and online software iPhyClassifier analysis of 16S rRNA gene sequences have become the widely accepted phytoplasma classification schemes . Pytoplasmas can affect a wide range of plant hosts, including agriculturally and economically important plants, such as fruit tree, landscape plant and flowers. Phytoplasma infected plants showing various symptoms ranging from mild yellowing to death. As wall-less obligate parasites, phytoplasma propagate intracellularly in both insects and plant phloem and transmit from plant to plant by insect vectors . Although axenic cultivation of phytoplasmas was reported , most of the phytoplasma strains are not successfully cultivated in vitro owe to the strict anaerobic conditions and complex media.
Whole-genome projects can provide insight into the organism biology and help understand host-pathogen interactions. However, the difficult cultivated in vitro of phytoplasmas poses a challenge to the collection of sufficiently high-quality phytoplasma genomic DNA and hinders attempts to study structural features and genome organization. Thus far, only six complete phytoplasma genomes have been fully sequenced [4,5,6,7,8,9], in which, three from 16Sr I group phytoplasmas (yellows mild strain OY-M, aster yellows witches’-broom isolate AY-WB and maize bushy stunt phytoplasma MBSP), two from 16SrXII phytoplasmas (Australian isolate PAa and strawberry lethal yellows isolate SLY of ‘Candidatus Phytoplasma australiense’) and one from 16Sr X phytoplasma (apple proliferation phytoplasma AT of ‘Ca. P. mali’). In addition, fifteen partial draft genomes, including members of 16SrI, 16SrII, 16SrIII, 16SrIX, 16SrXI and 16SrXII groups, have been reported [10,11,12,13,14,15,16,17,18,19,20,21]. Based on these sequences, more and more genomic information of phytoplasma were obtained. The genome length of phytoplasmas is between 530 kb and 1350 kb, and the size of the closely related strains varies significantly . The genomes have a relatively low GC content (23.0–29.5 Mol%), two rRNA operons, a small amount of tRNAs and a limited set of genes encoding metabolic enzymes .
Chinese Jujube (Ziziphus jujuba Mill.), a native and economically important fruit tree of China, is one of the most important members in the Rhamnaceae family. The complex jujube genome and more and more function genes were identified [23,24,25]. Jujube witches’-broom (JWB) disease caused by JWB phytoplasma widespread in China and causes serious problems to the industry. Until the phytoplasmal cells were demonstrated through transmission electron microscopy , JWB disease was described as a graft-transmissible viral disease . JWB phytoplasma now classified as ‘Ca. P. ziziphi’, belongs to elm yellows group (16SrV) subgroup B based on analyses of 16S rDNA sequences [28, 29]. The infected plants exhibit a variety of symptoms, such as little leaves, phyllody, witches’ broom, leaves etiolating and eventually death within a few years of infection . Previous studies mainly focused on the investigations of JWB, whereas very little knowledge is known about its pathogenicity mechanisms. Some progresses have been made in understanding the regulation of plant-phytoplasma interactions, and several candidate genes and proteins that might be involved in the interaction of Mexican lime trees with the phytoplasma were identified [31, 32]. Recently, microRNAs that may participate in the JWB phytoplasma response of jujube have been identified . Furthermore, combination of iTRAQ proteomics and RNA-seq transcriptomics reveals multiple levels of regulation in phytoplasma-infected jujube . However, the molecular mechanism of JWB disease is not fully understood. In this study, the complete genome sequence of JWB phytoplasma was reported, as well as and the relationship with other five phytoplasmas and the potential key genes were determined by comparative gene analysis. This work is of great importance for a comprehensive understanding of the interactions between JWB phytoplasmas and their hosts.
Samples of Z. jujuba (Jujube) with witches’-broom were collected from Shandong, China. The phytoplasma strain was then maintained and propagated in an insect-proof glasshouse by tissue culture. The transmitted JWB phytoplasma strain was confirmed by PCR using non-specific primers 16S rDNA-F (5’-TAAAAAGGCATCTTTTTGTT-3′) and 16S rDNA-R (5’-AATCCGGACTAAGACTGT-3′) that amplify the 16S rDNA sequence .
Genome DNA preparation and whole genome sequencing
JWB phytoplasma chromosomal DNA was prepared as described previously . JWB phytoplasma pellets were purified by differential centrifugation. Phytoplasma DNA were separated by pulsed-field gel electrophoresis (PFGE) after digestion by proteinase K at 50 °C for 72 h. Yeast chromosomes (New England Biolabs, Frankfurt, Germany) were used as a molecular size marker. The DNA was then electroeluted from the excised PFGE agars slice and concentrated by ethanol precipitation. To improve the accuracy of the genome sequences, Illumina HiSeq 4000 platform and PacBio RS II platform (10 kb inserts library; Pacific Biosciences, Menlo Park, CA, USA) were used for the JWB-nky sequencing at the Beijing Genomics Institute (BGI, Shenzhen, China). Four SMRT cells Zero-Mode Waveguide arrays of sequencing were used by the PacBio platform to generate the subreads set. PacBio subreads (length < 1 kb) were removed.
Genome assembly and annotation
PacBio sequencing reads were de novo assembled with hierarchical genome-assembly process (HGAP) . Sequencing data were combined, and further polished with Pion  to obtain the complete sequences. Removed genomic gene of plant based on the whole jujube genome sequence . Prodigal v2.60  was used to predict open reading frames (ORFs) for the JWB phytoplasma genome sequence. UGA was used as a stop codon, which was consistent with the ORF prediction for other phytoplasmas. tRNAscan-SE was applied to identify tRNAs  and RNAmmer  was used to predict the locations of rRNA genes. JWB phytoplasma genes were named according to their homologous genes and OrthoMCL  was used to identify homologous genes between the JWB phytoplasma genome and other complete phytoplasma genomes. Similarity searches were carried out using BLASTP  against the UniProt database. Those genes with no homology to the other complete phytoplasma genomes were searched by BLASTP against the NCBI nr database. The rest of the genes, which had no homology with the phytoplasma complete genomes or the nr database, were presumed to be putative protein-coding genes only when they were longer than 100 amino acids or had a confidence score of more than 10 from Prodigal. COG (Clusters of Orthologous Groups) , KEGG (Kyoto Encyclopedia of Genes and Genomes, http://www.genome.jp/kegg/), NR (Non-Redundant Protein Database) and GO (Gene Ontology) databases were used for general function annotation. The program signalP V4.0  was used to predict SP cleavage in all proteins. For SP-containing proteins, TMHMM V2.0  was used to predict the transmembrane regions after eliminating the SP. Virulence factors and resistance gene were identified based on the core dataset in VFDB (Virulence Factors of Pathogenic Bacteria) database, PHI (Pathogen Host Interactions) database and CAZy (Carbohydrate-Active enZYmes Database).
Comparative genome analysis
For comparative genomic analysis, five complete phytoplasma genomes (Ca.P. australiense [‘Ca. P. australiense’, GenBank accession number: GCF_000069925.1], Candidatus Phytoplasma mali [‘Ca. P. mali’, GenBank accession number: GCF_000026205.1], Candidatus Phytoplasma pruni [‘Ca. P. pruni’, GenBank accession number: GCA_001277135.1], Onion yellows phytoplasma [OY-M, GenBank accession number: GCA_000009845.1], Peanut witches’-broom phytoplasma [PnWB NTU2011, GenBank accession number: GCA_000364425.1]) were chosen for comparison. OrthoMCL was used for gene clustering analysis. All genes conserved among the six phytoplasma genomes were aligned using MUSCLE  with the default settings. The resulting multiple sequence alignment was used to infer the species phylogeny using the maximum likelihood program PhyML . The visualization of syntenic blocks was done in Web-based Genome Synteny Viewer . BLASTP was used for the alignment of protein sequences.
Nucleotide sequence accession number
The JWB phytoplasma chromosome sequence was deposited in DDBJ/EMBL/GenBank nucleotide sequence databases under accession number CP025121.
Results and discussion
General features of the JWB phytoplasma genome
The JWB phytoplasma genome consists of 750,803 bp with GC content of 23.31%. Based on the annotation, 694 protein-coding genes, two operons for rRNA genes and 31 tRNA genes were found in the genome (Table 1, Fig. 1). Genes in jwb-nky were annotated with different databases (See Additional file 1). 513 of the 694 annotated protein coding genes in jwb-nky can be annotated to Non-Redundant Protein Database based on homologous alignment with E-value< 10− 5. The genes homologous to Ca. Phytoplasma are the most abundant, accounting for 38.52% of total genes. 377 genes were categorized into cellular components, molecular functions, and biological processes according to GO (gene ontology) analysis (Fig. 2, Table 2). The most predominant cellular components were located in the cell (111) and cell part (111). The main molecular functions were binding activity (264) and catalytic activity (196). The most common biological processes were metabolic process (240), cellular process (236) and single-organism process (96). Other major functional categories were membrane (44), macromelocular complex (65), orgenalle (56) and structural molecular activity (52). Pathway analysis categorized 328 genes into 4 levels (Fig. 3, Table 2). It is worth noting that most genes in metabolism were associated with global and overview maps (61), carbohydrate metabolism (14) and energy metabolism (11); in genetic information processing were replication and repair (39) and translation (74); and in environmental information processing were membrane transport (20). 369 genes had specific functional ditributions according to the COG categories (Fig. 4, Table 2) and these 369 genes were assigned to 19 functional categories. The most abundant functional category was COG category J (translation, ribosomal structure and biogenesis) (120), accounting for 17.29% of the JWB genes, mainly due to 2 ribosomal proteins and 31 tRNA synthetases found in the genome. Followed by the COG category L (replication, recombination and repair), which contained 56 genes. A total of 537 protein-coding genes (77.37% of protein coding genes) were annotated to NR/GO/KEGG/COG database (See Additional file 1). As described in “Ca. P. asteris” strains OY-M and AY-WB and Ca. P. australiense , UGA was also used as a stop codon for ORF prediction.
Comparative genomic analysis of phytoplasmas
Clustering of gene family and core gene analysis
JWB phytoplasma has more multiple-copy orthologs with Ca. P. australiense and OY-M, and the genome size of these phytoplasmas is larger than that of Ca. P. mali, Ca. P. pruni and PnWB (Fig. 5), suggesting that expansion of multiple gene families/duplication occurred in these genomes and may play an important role in pathogenesis.
70 core genes were identified in all of these genomes: 12 from JWB, 14 from ‘Ca. P. australiense’, 9 from Ca. P. mali, 12 from Ca. P. pruni, 9 from OY-M and 14 from PnWB (See Additional file 2). In addition to these homologous genes, JWB phytoplasma also contained 307 specific genes, with no assigned function.
A total of 150 single copy gene proteins, common to all 6 phytoplasma genomes were identified (see Additional file 3). The phylogenetic tree based on the concatenated alignment of these genes showed that JWB phytoplasma were grouped closely with Ca. P. pruni and remotely with Ca. P. australiense (Fig. 6), indicating that the expansion of multiple gene families/duplication occurred separately after divergence.
Genomic synteny analysis between the JWB phytoplasma genome and other complete phytoplasma genomes were performed. Synteny between these genomes is quite poor. The alignment among the genomes of JWB phytoplasma, Ca. P. australiense, Ca. P. mali, OY-M and PnWB showed some of synteny regions (Fig. 7b, c, d, e). However, the genome alignments between JWB and Ca. P. pruni showed more synteny regions than other phytoplasmas (Fig. 7a). The result from synteny analysis is consistant with the phylogenetic relationships in which JWB and Ca. P. pruni were in the same group (Fig. 6).
The metabolic capabilities of many pathways are greately reduced in phytoplasmas, as is in the JWB phytoplasma. Investigation of the genome sequences showed that the pathways for oxidative phosphorylation, amino acid and fatty acid biosynthesis, pentose phosphate pathway and tricarboxylic acid cycle are not completed in JWB phytoplasma. It lacks genes involved in most amino acid sythesis, and only one gene involved in the reversible conversions of L-serine to glycine (jwb-nkyGM000266) was identified. And as is in other phytoplasmas, the pathways for ATP synthesis and fatty acid metabolism in JWB phytoplasma are incomplete. Only one gene involved in oxidative phosphorylation (inorganic pyrophosphatase, jwb-nkyGM000439) and seven genes related to glycerolipid and glycerophospholipid metabolism (gspA, jwb-nkyGM000052; plsY, jwb-nkyGM000048, pssA, jwb-nkyGM000102; psd, jwb-nkyGM000103; plsX, jwb-nkyGM000108; acpS, jwb-nkyGM000692 and cdsA, jwb-nkyGM000069) were identified. It also lacks a vitamin B6 metabolism gene (pdxK) and folate biosynthesis genes (folKPC), which were found in the OY-M phytoplasma genome.
None of the ATP-synthase subunits were identified in the JWB phytoplasma genome, while genes involved in glycolysis were detected (pgi, jwb-nkyGM000319; pfkA, jwb-nkyGM000444; fba, jwb-nkyGM000447; gapA, jwb-nkyGM000446; pgk, jwb-nkyGM000445; gpmI, jwb-nkyGM000443; eno, jwb-nkyGM000442; pykF, jwb-nkyGM000551; tpiA, jwb-nkyGM000540). suggesting that JWB phytoplasma depends on glycolysis for its energy generation. The phosphoenol pyruvate dependent sugar phosphotransferase (PTS) systems responsible for sugar importation and phosphorylation were not identified, whereas three genes encoding ABC-type maltose transporters (malF, jwb-nkyGM000324; malG, jwb-nkyGM000325; malK, jwb-nkyGM000323) were found in JWB phytoplasma. The maltose ABC transporter could also recognize sucrose and trehalose in some bacteria , which are major sugars in plant phloem and insect hemolymph, respectively. As in WBD phytoplasma , enzymes responsible for phosphorylation of glucose to produce glucose 6-phosphate were not identified in the JWB phytoplasma genome. Although a sucrose phosphorylase (jwb-nkyGM000581) catalyzing the conversion of sucrose to D-fructose and α-D-glucose-1-phosphate was found, no phosphoglucomutase gene was identified. A putative malate/citrate symporter (citS, jwb-nkyGM000689) and a malate dehydrogenase (sfcA, jwb-nkyGM000280) were found in other metabolic pathways that enabled the JWB phytoplasma to use malate, and this carbon source for ATP generation in phytoplasmas . No enzymes for the conversion of pyruvate to acetyl-CoA and acetyl-P/ADP to acetate/ATP were found. Although no phosphate acetyltransferase was found to catalyze acetyl-CoA to produce acetul-p, a putative phosphate propanoyltransferase (pduL, jwb-nkyGM000388) was identifiedthat may catalyze this reaction in phytoplasma.
Transporter and secretion systems
Phytoplasmas have strongly reduced metabolic capabilities, so they must absorb metabolites from their plant and insect hosts. The JWB phytoplasma genome contains 13 ORFs that encode transport and 264 ORFs encode binding proteins. Two ABC transporter genes (evbG, jwb-nkyGM000115 and mdlA, jwb-nkyGM000116) involved in multidrug resistance were identified. Others including ATP-binding cassette (ABC) transporter systems for dipeptides/oligopeptides (oppA, jwb-nkyGM000118; dppB, jwb-nkyGM000257 and dppD, jwb-nkyGM000255), spermidine/putrescine (potA, jwb-nkyGM000481; potB, jwb-nkyGM000480; potC, jwb-nkyGM000479; potD, jwb-nkyGM000478), cobalt (cbiO, jwb-nkyGM000001), Mn/Zn (znuB1, jwb-nkyGM000084; znuB2, jwb-nkyGM000085; znuC, jwb-nkyGM000086) and D-methionine (metN, jwb-nkyGM000558, metQ, jwb-nkyGM000559, metI jwb-nkyGM000560) were identified. However, amino acid transporter systems except for D-methionine such as GlnQ and ArtIMPQ that are present in the AY strains  are lacking in JWB phytoplasma, which is similar to Ca. P. mali . The methionine transport system may be also involved in other amino acids transportation.The components of protein translocation system were identified, including secA (jwb-nkyGM000372), secY (jwb-nkyGM000009), yidC (jwb-nkyGM000141), ffh (jwb-nkyGM000385), ftsY (jwb-nkyGM000384), dnaJ (jwb-nkyGM000334), dnaK (jwb-nkyGM000335), grpE (jwb-nkyGM000336), groEL (jwb-nkyGM000114) and signal peptidase SPaseI (jwb-nkyGM000063), which suggested that a functional sec-dependent protein translocation system exists in JWB phytoplasma.
Plant-pathogenic bacteria can secrete a subset of carbohydrate-activated enzymes (CAZymes) to degrade plant cell walls. The JWB phytoplasma genome encodes three CAZymes: jwb-nkyGM000581 (GH13), jwb-nkyGM000537 (GH4) and jwb-nkyGM000495 (GT78) (See Additional file 4). jwb-nkyGM000537 and jwb-nkyGM000495 are specific to JWB phytoplasma, while jwb-nkyGM000581 has a homologue in ‘Ca. P. australiense’ with 70.37% identify. These enzymes may be involved in phytoplasma virulence.
Virulence factors such as hemolysins and adhesion-related proteins are thought to be involved in bacterial pathogenicity. JWB phytoplasma genome contains two genes encoding putative hemolysins (jwb-nkyGM000059, jwb-nkyGM000572), this is similar to the results of hemolysine related proteins identified in AY-WB phytoplasma and Ca. p. australiense. The virulence factor TENGU , effectors SAP11  and SAP54  have been shown to alter plant morphology. 28 candidate secreted JWB proteins (SAP) were identified according Bai et al. (data no show) . 87 gene products in JWB phytoplasma have homologues in Virulence Factor Database (VFDB, http://www.mgc.ac.cn/VFs/main.htm) (See Additional file 5) and 92 gene products have homologues in Pathogen Host Interactions Database (PHI-base, http://www.phi-base.org/) using Blastp (E value< 10− 5) in which 13 gene products were in the class PHI:1566|GzOB006 and 11 genes were in the class PHI:2042|ABC3 (See Additional file 6). Among these factors, the SOD protein encoded by sodA (jwb-nkyGM000089), which was described as a virulence factor for Mycoplasma pneumoniae  was found. it has also been found in OY-W and WBD phytoplasma that SodA may play an important role in the phytoplasma colonization of plants and insects through degradation of reactive oxygen species (ROS) produced by the host during infection . In addition, molecular chaperones DnaK and DnaJ have been reported to be necessary for pathogen survival in the host cell . The multi-copy hflB (or ftsH) genes that encode membrane associated ATP-dependent Zn proteases also present in this group. These proteins are conserved among bacteria and are involved in protein secretion  and membrane protein assembly  as well as adaptations to nutritional conditions and osmotic stress [56, 58]. Other virulence-related factors in this group include AAA+ ATPases  and nucleases .
Besides virulence facotors, a large-conductance mechanosensitive channel MscL protein encoded by jwb-nkyGM000669 was identified. MscL involves in osmotic shock defense  and may help JWB phytoplasma adapt to changes in osmotic pressure between plant and insect cells. The presence of MscL, SodA, as well as DnaJ/DnaK suggests that JWB phytoplasma adapts to complex environments by different mechanisms.
Genome sequencing has revealed that many phytoplasma genomes contain large repetitive regions called potential mobile units (PMUs) that appear to have promoted intrachromosomal recombination [5, 7]. PMUs contain several genes involved in gene duplication and transposition, and some contain putative ‘virulence genes’. Genome alignment suggests that PMUs are involved in the instability and recombination of phytoplasma genome. Four PMUs were identified in JWB phytoplasma genome (Fig. 8, Table 3) and they showed the highest similar structure with AY-WB PMU1. The phylogenetic tree based on PMUs gene support different relationships between the JWB PMUs and homologs in the PnWB PMU and ‘Ca. P. asteris’ OY-M and AY-WB genomes, suggesting that potential mobile units are transferred among these phytoplasmas. JWB PMUs are more likely to have been transferred from the 16SrI group, due to PnWB PMU is originated from one of the PMUs present in a close relative of 16SrI group OY-M . PMUs translocation within the JWB genome resulted in the coexistence of different but homologous PMUs (PMU1–4), as in ‘Ca. P. australiense’ . Due to the close association of effectors and PMUs , PMU-mediated horizontal gene transfers among different phytoplasmas may have an impact on the evolution of pathogenicity .
Jujube witches’-broom (JWB) disease, caused by JWB phytoplasma, is one of the destructive diseases affecting worldwide jujube orchards. Some progresses have been made in understanding the regulation of Jujube-phytoplasma interactions. However, genomic basis and mechanism of its pathogenicity on Jujube remains to be researched. In this study, the genome of a JWB phytoplasma isolate jwb-nky was sequenced. Comparative genomic analysis revealed that, although jwb-nky is closely related to ‘Ca. P. pruni’ according to phylogenetic analysis, the gene syntety between the two phytoplasmas is low. And metabolic analysis suggested that some essential pathways of JWB phytoplasma are completely missing and others are greatly reduced.
Previous study based on 16S rDNA sequences analysis also indicates that the JWB phytoplasmas represents a novel taxon, ‘Ca. Phytoplasma ziziphi’ . There is few report on genome sequence of other JWB phytoplasma isolates and the information about them are not sufficient to access the envolutionary relationships extensively. Therefore, in order to have a better understanding about the genomic basis, evolution and pathogenicity mechanism of this phytoplasmas group, more genomes of JWB phytoplasma isolates need to be sequenced in the future.
Basic local alignment search tool-Protein
Clusters of Orthologous Groups
Hierarchical genome-assembly process
Kyoto Encyclopedia of Genes and Genomes
National Center for Biotechnology Information
Non-Redundant Protein Database
Ppen reading frames
Pulsed-field gel electrophoresis
Pathogen Host Interactions
Potential mobile unit
Virulence Factors of Pathogenic Bacteria
Zhao Y, Wei W, Lee IM, Shao J, Suo X, Davis RE. Construction of an interactive online phytoplasma classification tool, iPhy-classifier, and its application in analysis of the peach X-disease phytoplasma group (16SrIII). Int J Syst Evol Microbiol. 2009;59:2582–93.
Weintraub PG, Beanland L. Insect vectors of phytoplasmas. Annu Rev Entomol. 2006;51:91–111.
Contaldo N, Satta E, Zambon Y, Paltrinieri S, Bertaccini A. Development and evaluation of different complex media for phytoplasma isolation and growth. J Microbiol Meth. 2016;127:105–10.
Oshima K, Kakizawa S, Nishigawa H, Jung HY, Wei W, Suzuki S, et al. Reductive evolution suggested from the complete genome sequence of a plant-pathogenic phytoplasma. Nat Genet. 2004;36(1):7–29.
Bai X, Zhang J, Ewing A, Miller SA, Jancso Radek A, Shevchenko DV, et al. Living with genome instability: the adaptation of phytoplasmas to diverse environments of their insect and plant hosts. J Bacteriol. 2006;188(10):3682–96.
Tran-Nguyen LT, Kube M, Schneider B, Reinhardt R, Gibb KS. Comparative genome analysis of “Candidatus Phytoplasma australiense” (subgroup tuf-Australia I; rp-a) and “Ca. Phytoplasma asteris” strains OY-M and AY-WB. J Bacteriol. 2008;190(11):3979–91.
Andersen MT, Liefting LW, Havukkala I, Beever RE. Comparison of the complete genome sequence of two closely related isolates of’Candidatus Phytoplasma australiense’ reveals genome plasticity. BMC Genomics. 2013;14:529.
Kube M, Schneider B, Kuhl H, Dandekar T, Heitmann K, Migdoll AM, Reinhardt R, et al. The linear chromosome of the plant-pathogenic mycoplasma’Candidatus Phytoplasma mali’. BMC Genomics. 2008;9:306.
Orlovskis Z, Canale MC, Haryono M, Lopes JRS, Kuo C-H, Hogenhout SA. A few sequence polymorphisms among isolates of maize bushy stunt phytoplasma associate with organ proliferation symptoms of infected maize plants. Ann Bot. 2017;119:869–84.
Saccardo F, Martini M, Palmano S, Ermacora P, Scortichini M, Loi N, Firrao G. Genome drafts of four phytoplasma strains of the ribosomal group 16SrIII. Microbiology. 2012;158(Pt 11):2805–14.
Chung WC, Chen LL, Lo WS, Lin CP, Kuo CH. Comparative analysis of the peanut witches’-broom phytoplasma genome reveals horizontal transfer of potential mobile units and effectors. PLoS One. 2013;8(4):e62770.
Chen W, Li Y, Wang Q, Wang N, Wu Y. Comparative genome analysis of wheat blue dwarf phytoplasma, an obligate pathogen that causes wheat blue dwarf disease in China. PLoS One. 2014;9(5):e96436.
Kakizawa S, Makino A, Ishii Y, Tamaki H, Kamagata Y. Draft genome sequence of “Candidatus Phytoplasma asteris” strain OY-V, an unculturable plant-pathogenic bacterium. Genome Announc. 2014;2(5):e00944–14.
Mitrović J, Siewert C, Duduk B, Hecht J, Mölling K, Broecker F, et al. Generation and analysis of draft sequences of “stolbur” phytoplasma from multiple displacement amplification templates. J Mol Microbiol Biotechnol. 2014;24:1–11.
Chang SH, Cho ST, Chen CL, Yang JY, Kuo CH. Draft genome sequence of a 16SrII-A subgroup phytoplasma associated with purple coneflower (Echinacea purpurea) witches’-broom disease in Taiwan. Genome Announc. 2015;3:e01398–15.
Lee IM, Shao J, Bottner-Parker KD, Gundersen-Rindal DE, Zhao Y, Davis RE. Draft genome sequence of’Candidatus Phytoplasma pruni’ strain CX, a plant-pathogenic bacterium. Genome Announc. 2015;3:e01117–5.
Quaglino F, Kube M, Jawhari M, Abou-Jawdah Y, Siewert C, Choueiri E, et al. Candidatus Phytoplasma phoenicium’ associated with almond witches’-broom disease: from draft genome to genetic diversity among strain populations. BMC Microbiol. 2015;15:148.
Fischer A, Santana-Cruz I, Wambua L, Olds C, Midega C, Dickinson M, et al. Draft genome sequenceof’Candidatus Phytoplasma oryzae’ strain Mbita1, the causative agent of Napier grass stunt disease in Kenya. Genome Announc. 2016;4:e00297–16.
Zamorano A, Fiore N. Draft genome sequence of 16SrIII-J phytoplasma, a plant pathogenic bacterium with a broad spectrum of hosts. Genome Announc. 2016;4:e00602–16.
Zhu Y, He Y, Zheng Z, Chen J, Wang Z, Zhou G. Draft genome sequence of rice orange leaf phytoplasma from Guangdong. China Genome Announc. 2017;5:e00430–17.
Sparks ME, Bottner-Parker KD, Gundersen-Rindal DE, Lee IM. Draft genome sequence of the New Jersey aster yellows strain of’Candidatus Phytoplasma asteris’. PLoS One. 2018;13(2):e0192379.
Marcone C, Neimark H, Ragozzino A, Lauer U, Seemüller E. Chromosome sizes of phytoplasmas composing major phylogenetic groups and subgroups. Phytopathology. 1999;89(9):805–10.
Liu MJ, Zhao J, Cai QL, Liu GC, Wang JR, Zhao ZH, et al. The complex jujube genome provides insights into fruit tree biology. Nat Commun. 2014;5:5315.
Huang J, Zhang C, Zhao X, Fei Z, Wan K, Zhang Z, et al. The jujube genome provides insights into genome evolution and the domestication of sweetness/acidity taste in fruit trees. PLoS Genet. 2016;12(12):e1006433.
Song S, Zhou H, Sheng S, Cao M, Li Y, Pang X. Genome-wide organization and expression profiling of the SBP-box gene family in Chinese jujube (Ziziphus jujuba mill.). Int J Mol Sci. 2017;18(8):1734.
Yi CK, La YJ. Mycoplasma-like bodies found in the phloem elements of jujube infected with witches’- broom disease. J Korean For Soc. 1973;20:111–4.
Kim CJ. Witches’-broom of jujube tree (Zizyphus jujuba Mill. Var. inermis Rehd.). Transmission by grafting. Korean J Microbiol. 1965;3:1–6.
Jung HY, Sawayanagi T, Kakizawa S, Nishigawa H, Wei W, Oshima K, et al. ‘Candidatus Phytoplasma ziziphi’, a novel phytoplasma taxon associated with jujube witches’-broom disease. Int J Syst Evol Micro. 2003;53(4):1037–41.
Lee IM, Martini M, Marcone C, Zhu SF. Classification of phytoplasma strains in the elm yellows group (16SrV) and proposition of ‘Candidatus Phytoplasma ulmi’ for the phytoplasma associated with elm yellows. Int J Syst Evol Micro. 2004;54(2):337–47.
Liu X, Liu M, Ning Q, Liu G. Reverse-cleft in vitro micrografting of Ziziphus jujuba mill. Infected with jujube witches’-broom (JWB). Plant Cell Tiss Org. 2012;108(2):339–44.
Monavarfeshani A, Mirzaei M, Sarhadi E, Amirkhani A, Khayam Nekouei M, Haynes PA, et al. Shotgun proteomic analysis of the Mexican lime tree infected with “Candidatus Phytoplasma aurantifolia”. J Proteome Res. 2013;12(2):785–95.
Zamharir MG, Mardi M, Alavi SM, Hasanzadeh N, Nekouei MK, Zamanizadeh HR, et al. Identification of genes differentially expressed during interaction of Mexican lime tree infected with “Candidatus Phytoplasma aurantifolia”. BMC Microbiol. 2011;11:1.
Shao F, Zhang Q, Liu H, Lu S, Qiu D. Genome-wide identification and analysis of microRNAs involved in witches’-broom phytoplasma response in ziziphus jujuba. PLoS One. 2016;11(11):e0166099.
Ye X, Wang H, Chen P, Fu B, Zhang M, Li J, et al. Combination of iTRAQ proteomics and RNA-seq transcriptomics reveals multiple levels of regulation in phytoplasma-infected Ziziphus jujuba mill. Hortic Res. 2017;4:17080.
Chin CS, Alexander DH, Marks P, Klammer AA, Drake J, Heiner C, Clum A, et al. Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data. Nat Methods. 2013;10(6):563–9.
Walker BJ, Abeel T, Shea T, Priest M, Abouelliel A, Sakthikumar S, et al. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS One. 2014;9(11):e112963.
Hyatt D, Chen GL, Locascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics. 2010;11:119.
Lowe TM, Eddy SR. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997;25(5):955–64.
Lagesen K, Hallin P, Rodland EA, Staerfeldt HH, Rognes T, Ussery DW. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007;35(9):3100–8.
Li L, Stoeckert CJ Jr, Roos DS. OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 2003;13(9):2178–89.
Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25(17):3389–402.
Tatusov RL, Fedorova ND, Jackson JD, Jacobs AR, Kiryutin B, Koonin EV, et al. The COG database: an updated version includes eukaryotes. BMC Bioinformatics. 2003;4:41.
Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8(10):785–6.
Krogh A, Larsson B, von Heijne G, Sonnhammer EL. Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J Mol Biol. 2001;305(3):567–80.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7.
Guindon S, Gascuel O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003;52(5):696–704.
Revanna KV, Chiu CC, Bierschank E, Dong Q. GSV: a web-based genome synteny viewer for customized data. BMC Bioinformatics. 2011;12:316.
Silva Z, Sampaio MM, Henne A, Böhm A, Gutzat R, Boos W, et al. The high-affinity maltose/trehalose ABC transporter in the extremely thermophilic bacterium Thermus thermophilus HB27 also recognizes sucrose and palatinose. J Bacteriol. 2005;187(4):1210–8.
Hoshi A, Oshima K, Kakizawa S, Ishii Y, Ozeki J, Hashimoto M, et al. A unique virulence factor for proliferation and dwarfism in plants identified from a phytopathogenic bacterium. Proc Natl Acad Sci. 2009;106(15):6416–21.
Sugio A, Kingdom HN, MacLean AM, Grieve VM, Hogenhout SA. Phytoplasma protein effector SAP11 enhances insect vector reproduction by manipulating plant development and defense hormone biosynthesis. Proc Natl Acad Sci. 2011;108(48):1254–63.
MacLean AM, Sugio A, Makarova OV, Findlay KC, Grieve VM, Tóth R, et al. Phytoplasma effector SAP54 induces indeterminate leaf-like flower development in Arabidopsis plants. Plant Physiol. 2011;157(2):831–41.
Bai X, Correa VR, Toruño TY, Ammar e-D, Kamoun S, Hogenhout SA. AY-WB phytoplasma secretes a protein that targets plant cell nuclei. Mol Plant-Microbe Interact. 2009;22(1):18–30.
Cohen G, Somerson NL. Glucose-dependent secretion and destruction of hydrogen peroxide by mycoplasma pneumoniae. J Bacteriol. 1969;98(2):547–51.
Miura C, Sugawara K, Neriya Y, Minato N, Keima T, Himeno M, et al. Functional characterization and gene expression profiling of superoxide dismutase from plant pathogenic phytoplasma. Gene. 2012;510(2):107–12.
Takaya A, Tomoyasu T, Matsui H, Yamamoto T. The DnaK/DnaJ chaperone machinery of Salmonella enterica serovar typhimurium is essential for invasion of epithelial cells and survival within macrophages, leading to systemic infection. Infect Immun. 2004;72(3):1364–73.
Deuerling E, Mogk A, Richter C, Purucker M, Schumann W. The ftsH gene of Bacillus subtilis is involved in major cellular processes such as sporulation, stress adaptation and secretion. Mol Microbiol. 1997;23(5):921–33.
Akiyama Y, Shirai Y, Ito K. Involvement of FtsH in protein assembly into and through the membrane II. Dominant mutations affecting FtsH functions. J Biol Chem. 1994;269(7):5225–9.
Lysenko E, Ogura T, Cutting SM. Characterization of the ftsH gene of Bacillus subtilis. Microbiol. 1997;143(3):971–8.
Seemüller E, Sule S, Kube M, Jelkmann W, Schneider B. The AAA+ ATPases and HflB/FtsH proteases of’Candidatus Phytoplasma mali’: phylogenetic diversity, membrane topology, and relationship to strain virulence. Mol Plant-Microbe Interact. 2013;26(3):367–76.
Bendjennat M, Blanchard A, Loutfi M, Montagnier L, Bahraoui E. Role of Mycoplasma penetrans endonuclease P40 as a potential pathogenic determinant. Infect Immun. 1999;67(9):4456–62.
Kloda A, Martinac B. Common evolutionary origins of mechanosensitive ion channels in archaea, Bacteria and cell-walled Eukarya. Archaea. 2002;1(1):35–44.
Ku C, Lo WS, Kuo CH. Horizontal transfer of potential mobile units in phytoplasma. Mob Genet Elments. 2013;3:e26145.
Sugio A, Hogenhout SA. The genome biology of phytoplasma: modulators of plants and insects. Curr Opin Microbiol. 2012;15(3):247–54.
Funds for National Natural Science Foundation of China (31401718 and 31701899); Excellent Young Scientists of Shandong Province (BS2013NY012); Young Talents Training Program of Shandong Academy of Agricultural Sciences and Youth Fund of Shandong Academy of Agricultural Sciences (2016YQN28). The funding body was not involved in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Availability of data and materials
Ethics approval and consent to participate
The diseased jujube trees were collected from an orchard in Shandong province with no permissions since it is not an endangered fruit trees in China.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Gene annotation of jwb-nky with different databases. (XLS 468 kb)
The identity of core genes of the six phytoplasmas (XLS 6 kb)
Single copy genes in six phytoplasmas. (XLS 41 kb)
Annotated CAZymes in jwb-nky. (XLS 1 kb)
87 putative effectors in JWB phytoplasma genome assigned to Virulence Factor Database. (XLS 32 kb)
92 genes in JWB phytoplasma genome assigned to Pathogen Host Interactions Database. (XLS 46 kb)
About this article
Cite this article
Wang, J., Song, L., Jiao, Q. et al. Comparative genome analysis of jujube witches’-broom Phytoplasma, an obligate pathogen that causes jujube witches’-broom disease. BMC Genomics 19, 689 (2018). https://doi.org/10.1186/s12864-018-5075-1
- Genome sequence
- Potential mobile unit (PMU)