Draft Genome Sequences Resources of Mulberry Dwarf Phytoplasma Strain MDGZ-01 Associated with Mulberry Yellow Dwarf (MYD) Diseases.
Longhui Luo,Xingnan Zhang,Fang Meng,Yaqin Wang,Yinan Zhou,Jiping Liu
DOI: https://doi.org/10.1094/pdis-12-21-2671-a
IF: 4.5
2022-01-01
Plant Disease
Abstract:HomePlant DiseaseVol. 106, No. 8Draft Genome Sequences Resources of Mulberry Dwarf Phytoplasma Strain MDGZ-01 Associated with Mulberry Yellow Dwarf (MYD) Diseases PreviousNext RESOURCE ANNOUNCEMENT OPENOpen Access licenseDraft Genome Sequences Resources of Mulberry Dwarf Phytoplasma Strain MDGZ-01 Associated with Mulberry Yellow Dwarf (MYD) DiseasesLonghui Luo, Xingnan Zhang, Fang Meng, Yaqin Wang, Yinan Zhou, and Jiping LiuLonghui Luohttps://orcid.org/0000-0002-1050-1031College of Animal Science, Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding, South China Agriculture University, Wushan Road, Guangzhou, Guangdong, 510642, China, Xingnan ZhangCollege of Animal Science, Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding, South China Agriculture University, Wushan Road, Guangzhou, Guangdong, 510642, China, Fang MengCollege of Animal Science, Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding, South China Agriculture University, Wushan Road, Guangzhou, Guangdong, 510642, China, Yaqin WangCollege of Animal Science, Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding, South China Agriculture University, Wushan Road, Guangzhou, Guangdong, 510642, China, Yinan ZhouCollege of Animal Science, Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding, South China Agriculture University, Wushan Road, Guangzhou, Guangdong, 510642, China, and Jiping Liu†Corresponding author: J. Liu; E-mail Address: liujiping@scau.edu.cnCollege of Animal Science, Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding, South China Agriculture University, Wushan Road, Guangzhou, Guangdong, 510642, ChinaAffiliationsAuthors and Affiliations Longhui Luo Xingnan Zhang Fang Meng Yaqin Wang Yinan Zhou Jiping Liu † College of Animal Science, Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding, South China Agriculture University, Wushan Road, Guangzhou, Guangdong, 510642, China Published Online:21 Jun 2022https://doi.org/10.1094/PDIS-12-21-2671-AAboutSectionsPDF ToolsAdd to favoritesDownload CitationsTrack Citations ShareShare onFacebookTwitterLinked InRedditEmailWechat Genome Announcement‘Candidatus Phytoplasma asteris’ is the causal agent of mulberry yellow dwarf (MYD) disease, which is present in many sericulture regions and cause great economic losses to the sericulture industry. Due to the limitation of genomic resources, little is known about the molecular pathogenesis of ‘Ca. P. asteris’. Mulberry dwarf phytoplasma strain MDGZ-01, the pathogen of MYD disease, was identified from the leaves of infected mulberry in Guangdong, in the southern part of China. Here, we used Illumina NovaSeq sequencing technology to sequence and analyze the whole genome of MDGZ-01. The complete genome sequence of MDGZ-01 was characterized by a 622,358-base circular chromosome with a 29.10% GC content. This genome provides a resource to explore the pathogenetic mechanisms of MYD disease. To our knowledge, this is the first genome announcement of a mulberry-associated Phytoplasma.Mulberry (Morus alba L.) is a woody plant belonging to the family Moraceae; it is widespread in Asia, Africa, and Europe and has diverse beneficial characteristics (Dhanyalakshmi and Nataraja 2018). In recent years, MYD disease was observed increasingly in most mulberry yards in China, which led to huge economic losses in sericulture regions. The earliest records of phytoplasma diseases can be traced back to the Song Dynasty, about 1,000 years ago (Wang and Mamamorosch 1988). Doi et al. (1967) reported that phytoplasma was first observed in phloem cells of infected mulberry by electron microscopy. In 2004, the International Research Program on Comparative Mycoplasmology (IRPCM) issued a guide for ‘Candidatus Phytoplasma’, and pointed out that phytoplasma identification should be based on 16 ribosomal RNA (16 rRNA) sequences and biological characteristics (IRPCM 2004). It is usually required that the 16 rRNA gene of plant phytoplasma (≥1,200 nucleotides), homology ≥97.5%, is classified as the same species. Currently, 33 taxa and more than 120 subgroups have been assigned the temporary taxonomic status of ‘Candidatus’ (Harrison et al. 2018). In 2020, a sample of mulberry leaves exhibiting MYD was collected from Guangzhou, Guangdong, China. The presence of phytoplasmas was confirmed by detection and sequencing of the 16S rRNA gene amplicons from the mulberry samples. The strain was designated MDGZ-01 and the voucher specimens were deposited in the herbarium of the Asia-Pacific Sericulture Training Center (Guangzhou, China) with sample number MDGZ. Here, we reported the first complete genome sequence of Mulberry dwarf phytoplasma strain MDGZ-01, providing a resource for comparative genomics of phytoplasmas.Genomic DNA of strain MDGZ-01 was extracted from symptomatic mulberry leaves using an EZ-10 Spin Column Plant Genomic DNA Purification Kit (Sangon Biotech Co., Ltd., Shanghai, China). DNA quantity was assessed using a NanoDrop 1000 spectrophotometer (Thermo Fisher Scientific) and a Qubit 2.0 fluorometer (Life Technologies). DNA integrity was assessed using the the quantitative PCR method and Agilent 2100 BioAnalyzer system (Agilent). Bacterial DNA was enriched using the NEBNext microbiome DNA enrichment kit (New England Biolabs, Ipswich, MA, U.S.A.) that was further amplified using the REPLI-g Mini Kit (Qiagen). For paired-end sequencing, whole-genome DNA was used to prepare a DNA library with a 350-bp insert. Whole-genome sequencing was conducted with manufactured 2 × 150-bp alignment reads on an Illumina NovaSeq System by Science Corporation of Gene Technology Co., Ltd. (Guangzhou, China). The completeness of genome was evaluated using benchmarking universal single-copy ortholog (BUSCO) (V5.2.2) software, with mollicutes_odb10 database and default parameters (Simão et al. 2015).Genome sequences were assembled using SPAdes, version 3.5.0 (Bankevich et al. 2012). The genome was automatically annotated using the NCBI Prokaryotic Genome Annotation Pipeline (Haft et al. 2018) and the RAST tool kit (Overbeek et al. 2014). Functional analysis was done by using clusters of orthologous groups (COGs) of proteins (Galperin et al. 2015), gene ontology (GO) (Ashburner et al. 2000), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways (Kanehisa et al. 2016) annotation. The chromosomal genome of Mulberry dwarf phytoplasma strain MDGZ-01 is 622,358 bp in length with 29.10% GC content, and is characterized by 76 contigs and 67 scaffolds, with a contig N50 of 32,555 bp and a scaffold N50 of 32,555 bp. The resulting genome contains 580 genes in total, with 524 protein-coding genes; 15 pseudogenes; 2, 2, and 2 rRNA-coding sequences (5S, 16S, and 23S rRNAs, respectively); 32 transfer RNAs; and 3 noncoding RNAs. By using functional protein alignments, 388 proteins were assigned to the COG database, 110 proteins were assigned to GO terms, and 78 proteins were mapped to KEGG pathways. The characteristics of genome sequencing and BUSCO evaluation result are summarized in Table 1.Table 1. Genome of Mulberry dwarf phytoplasma strain MDGZ-01aFeaturesNumbers (%)Genome size (bp)6,22,358Genome coverage1,139.0×GC content (%)29.10Gene (total)580Gene (coding)524Transfer RNAs32Ribosomal RNAs6Noncoding RNAs3Pseudogenes15Hypothetical proteins246Proteins with COG assignments388Proteins with KEGG assignments78Proteins with GO assignments110Genome assembly Contig76 N50 (bp)24,801 N90 (bp)2,598 Scaffold67 N50 (bp)32,555 N90 (bp)2,700BUSCOs Complete BUSCOs144 (95.4) Complete and single-copy BUSCOs143 (94.7) Complete and duplicated BUSCOs1 (0.7) Fragmented BUSCOs1 (0.7) Missing BUSCOs6 (4) Total BUSCO groups searched151aAbbreviations: COG = clusters of orthologous groups, KEGG = Kyoto Encyclopedia of Genes and Genomes, GO = gene ontology, and BUSCO = benchmarking universal single-copy ortholog.Table 1. Genome of Mulberry dwarf phytoplasma strain MDGZ-01aView as image HTML To validate the phylogenetic position of MDGZ-01, the closest reference and representative genomes, registered by NCBI, were identified by Mash/MinHash (Ondov et al. 2016) and selected to construct a phylogenetic tree. The single-copy orthologous genes of MDGZ-01 and another 12 phytoplasmas were selected by OrthoFinder (v2.0) software (Emms and Kelly 2019), then connected into a single super gene. We used jModelTest (v2.1.10) (Darriba et al. 2012) to test the nucleic acid model of the super gene DNA sequence, and selected the minimum value of Akaike Information Criterion as the best model for tree framing. The super gene was exported to MEGA 11 software (v 11.0.10) (Kumar et al. 2018) to build maximum-likelihood (ML) phylogenetic trees. ML tree searches were performed under the GTRGAMMA model with 1,000 pseudoreplicates. In the ML tree (Fig. 1), phytoplasma can be divided into at least four groups, including ‘Candidatus Phytoplasma’ 16Sr I group, ‘Candidatus Phytoplasma’ 16Sr XII (Stolbur group), ‘Candidatus Phytoplasma’ 16Sr X (host: apple), and ‘Candidatus Phytoplasma’ 16Sr V (host: jujube) group. Mulberry dwarf phytoplasma strain MDGZ-01 was be classified into the ‘Candidatus Phytoplasma’ 16Sr I group. The phylogenetic tree showed that MDGZ-1 was closest to ‘Ca. P. asteris’ OY-M (AP006628) and had the highest similarity (99% genome-wide similarity and 90% coverage). High bootstrap and posterior probability values indicate that the presented relationship is stable. In addition, MDGZ-01 and ‘Ca. P. asteris’ CR-DV (NZ_CP035949), ‘Ca. P. asteris’ MBS-M3 (NZ_CP015149), ‘Ca. P. asteris’ ROL-LD1 (NZ_MIEP01000004), ‘Ca. P. asteris’ WBD (NZ_AVAO 01000003), and ‘Ca. P. asteris’ AYWB (NC_007716) were clustered in the same branch and belonged to the ‘Candidatus Phytoplasma’ 16Sr I group. Therefore, it is necessary to consider whether there is the possibility of cross infection of phytophores in mulberry intercropped or interplanted cash crops. In a word, our data provide a reference for the study of phytoplasma diseases.Fig. 1. Phylogenetic relationships among 13 phytoplasma genomes, including ‘Candidatus Phytoplasma asteris’. A maximum-likelihood (GTRGAMMA model) tree, based on 71 single-copy orthologous genes, from Mulberry dwarf phytoplasma isolates was obtained from Morus spp. and reference isolates from GenBank. Bootstrap support values above 50% (left) are shown at the nodes.Download as PowerPointData AvailabilityThe genome sequence of Mulberry dwarf phytoplasma strain MDGZ-01 has been deposited in NCBI GenBank in BioProject number PRJNA773635 under accession number CP085837.1.The author(s) declare no conflict of interest.Literature CitedAshburner, M., Ball, C. A., Blake, J. A., Botstein, D., Butler, H., Cherry, J. M., Davis, A. P., Dolinski, K., Dwight, S. S., Eppig, J. T., Harris, M. A., Hill, D. P., Issel-Tarver, L., Kasarskis, A., Lewis, S., Matese, J. C., Richardson, J. E., Ringwald, M., Rubin, G. M., and Sherlock, G. 2000. Gene Ontology: Tool for the unification of biology. Nat. Genet. 25:25-29. https://doi.org/10.1038/75556 Crossref, ISI, Google ScholarBankevich, A., Nurk, S., Antipov, D., Gurevich, A. A., Dvorkin, M., Kulikov, A. S., Lesin, V. M., Nikolenko, S. I., Pham, S., Prjibelski, A. D., Pyshkin, A. V., Sirotkin, A. V., Vyahhi, N., Tesler, G., Alekseyev, M. A., and Pevzner, P. A. 2012. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19:455-477. https://doi.org/10.1089/cmb.2012.0021 Crossref, ISI, Google ScholarDarriba, D., Taboada, G. L., Doallo, R., and Posada, D. 2012. jModelTest 2: More models, new heuristics and parallel computing. Nat. Methods 9:772. https://doi.org/10.1038/nmeth.2109 Crossref, ISI, Google ScholarDhanyalakshmi, K. H., and Nataraja, K. N. 2018. Mulberry (Morus spp.) has the features to treat as a potential perennial model system. Plant Signal. Behav. 13:e1491267. https://doi.org/10.1080/15592324.2018.1491267 ISI, Google ScholarDoi, Y., Teranaka, M., and Yora, K. 1967. Mycoplasma- or PLT group-like microorganisms found in the phloem elements of plants infected with mulberry dwarf, potato witches’ broom, aster yellows, or paulownia witches’ broom. Jpn. J. Phytopathol. 33:259-266. https://doi.org/10.3186/jjphytopath.33.259 Crossref, Google ScholarEmms, D. M., and Kelly, S. 2019. OrthoFinder: Phylogenetic orthology inference for comparative genomics. Genome Biol. 20:238. https://doi.org/10.1186/s13059-019-1832-y Crossref, ISI, Google ScholarGalperin, M. Y., Makarova, K. S., Wolf, Y. I., and Koonin, E. V. 2015. Expanded microbial genome coverage and improved protein family annotation in the COG database. Nucleic Acids Res. 43:D261-D269. https://doi.org/10.1093/nar/gku1223 Crossref, ISI, Google ScholarHaft, D. H., DiCuccio, M., Badretdin, A., Brover, V., Chetvernin, V., O’Neill, K., Li, W., Chitsaz, F., Derbyshire, M. K., Gonzales, N. R., Gwadz, M., Lu, F., Marchler, G. H., Song, J. S., Thanki, N., Yamashita, R. A., Zheng, C., Thibaud-Nissen, F., Geer, L. Y., Marchler-Bauer, A., and Pruitt, K. D. 2018. RefSeq: An update on prokaryotic genome annotation and curation. Nucleic Acids Res. 46:D851-D860. https://doi.org/10.1093/nar/gkx1068 Crossref, ISI, Google ScholarHarrison, N. A., Gundersen-Rindal, D., and Davis, R. E. 2018. Candidatus phytoplasma. In: Bergey’s Manual of Systematics of Archaea and Bacteria. https://doi.org/10.1002/9781118960608.gbm01259.pub2 Google ScholarIRPCM. 2004. ‘Candidatus Phytoplasma’, a taxon for the wall-less, non-helical prokaryotes that colonize plant phloem and insects. Int. J. Syst. Evol. Microbiol. 54:1243-1255. https://doi.org/10.1099/ijs.0.02854-0 Crossref, ISI, Google ScholarKanehisa, M., Sato, Y., Kawashima, M., Furumichi, M., and Tanabe, M. 2016. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 44:D457-D462. https://doi.org/10.1093/nar/gkv1070 Crossref, ISI, Google ScholarKumar, S., Stecher, G., Li, M., Knyaz, C., and Tamura, K. 2018. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 35:1547-1549. https://doi.org/10.1093/molbev/msy096 Crossref, ISI, Google ScholarOndov, B. D., Treangen, T. J., Melsted, P., Mallonee, A. B., Bergman, N. H., Koren, S., and Phillippy, A. M. 2016. Mash: Fast genome and metagenome distance estimation using MinHash. Genome Biol. 17:132. https://doi.org/10.1186/s13059-016-0997-x Crossref, ISI, Google ScholarOverbeek, R., Olson, R., Pusch, G. D., Olsen, G. J., Davis, J. J., Disz, T., Edwards, R. A., Gerdes, S., Parrello, B., Shukla, M., Vonstein, V., Wattam, A. R., Xia, F., and Stevens, R. 2014. The SEED and the rapid annotation of microbial genomes using subsystems technology (RAST). Nucleic Acids Res. 42:D206-D214. https://doi.org/10.1093/nar/gkt1226 Crossref, ISI, Google ScholarSimão, F. A., Waterhouse, R. M., Ioannidis, P., Kriventseva, E. V., and Zdobnov, E. M. 2015. BUSCO: Assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31:3210-3212. https://doi.org/10.1093/bioinformatics/btv351 Crossref, ISI, Google ScholarWang, M. Q., and Maramorosch, K. 1988. Earliest historical record of a tree mycoplasma disease: Beneficial effect of mycoplasma-like organisms on peonies. Pages 349-356 in: Mycoplasma Diseases of Crops. K. Maramorosch and S. P. Raychaudhuri, eds. Springer, New York, NY, U.S.A. https://doi.org/10.1007/978-1-4612-3808-9_19 Crossref, Google ScholarFunding: This work was supported by the Agriculture Research System of China MOF and MARA (CARS-18-ZJ0304).The author(s) declare no conflict of interest.DetailsFiguresLiterature CitedRelated Vol. 106, No. 8 August 2022SubscribeISSN:0191-2917e-ISSN:1943-7692 Download Metrics Article History Issue Date: 29 Jul 2022Published: 21 Jun 2022Accepted: 6 Apr 2022 Pages: 2239-2242 Information© 2022 The American Phytopathological SocietyFundingAgriculture Research System of ChinaGrant/Award Number: CARS-18-ZJ0304Keywords‘Candidatus Phytoplasma asteris’genomeIllumina NovaSeqMulberry dwarf phytoplasmaMulberry yellow dwarf diseaseThe author(s) declare no conflict of interest.PDF download