Complete Genome Resource of a Commonly Used Laboratory Substrain of Xanthomonas oryzae pv. oryzae PXO99(A)
Yiming Li,Xue Dan,Tianjiao Li,Qiaozhi Wu,Zihao Wang,Lidan Xu,Xiuguo Ma,Fengzhi Bie,Shimin Lei,Sheng Huang,Wei Jiang,Yong-Qiang He
DOI: https://doi.org/10.1094/PDIS-10-21-2176-A
IF: 4.5
2022-01-01
Plant Disease
Abstract:HomePlant DiseaseVol. 106, No. 3Complete Genome Resource of a Commonly Used Laboratory Substrain of Xanthomonas oryzae pv. oryzae PXO99A PreviousNext RESOURCE ANNOUNCEMENT OPENOpen Access licenseComplete Genome Resource of a Commonly Used Laboratory Substrain of Xanthomonas oryzae pv. oryzae PXO99AYiming Li, Xue Dan, Tianjiao Li, Qiaozhi Wu, Zihao Wang, Lidan Xu, Xiuguo Ma, Fengzhi Bie, Shimin Lei, Sheng Huang, Wei Jiang, and Yong-Qiang HeYiming LiNational Demonstration Center for Experimental Plant Science Education, College of Agriculture, Guangxi University, Nanning 530004, Guangxi, ChinaState Key Laboratory for Conservation and Utilization of Subtropical Agro-bioresources and College of Life Science and Technology, 100 Daxue Road, Nanning, Guangxi 530004, ChinaSearch for more papers by this author, Xue DanNational Demonstration Center for Experimental Plant Science Education, College of Agriculture, Guangxi University, Nanning 530004, Guangxi, ChinaSearch for more papers by this author, Tianjiao LiNational Demonstration Center for Experimental Plant Science Education, College of Agriculture, Guangxi University, Nanning 530004, Guangxi, ChinaState Key Laboratory for Conservation and Utilization of Subtropical Agro-bioresources and College of Life Science and Technology, 100 Daxue Road, Nanning, Guangxi 530004, ChinaSearch for more papers by this author, Qiaozhi WuState Key Laboratory for Conservation and Utilization of Subtropical Agro-bioresources and College of Life Science and Technology, 100 Daxue Road, Nanning, Guangxi 530004, ChinaSearch for more papers by this author, Zihao WangState Key Laboratory for Conservation and Utilization of Subtropical Agro-bioresources and College of Life Science and Technology, 100 Daxue Road, Nanning, Guangxi 530004, ChinaSearch for more papers by this author, Lidan XuState Key Laboratory for Conservation and Utilization of Subtropical Agro-bioresources and College of Life Science and Technology, 100 Daxue Road, Nanning, Guangxi 530004, ChinaSearch for more papers by this author, Xiuguo MaNational Demonstration Center for Experimental Plant Science Education, College of Agriculture, Guangxi University, Nanning 530004, Guangxi, ChinaSearch for more papers by this author, Fengzhi BieState Key Laboratory for Conservation and Utilization of Subtropical Agro-bioresources and College of Life Science and Technology, 100 Daxue Road, Nanning, Guangxi 530004, ChinaSearch for more papers by this author, Shimin LeiNational Demonstration Center for Experimental Plant Science Education, College of Agriculture, Guangxi University, Nanning 530004, Guangxi, ChinaSearch for more papers by this author, Sheng HuangState Key Laboratory for Conservation and Utilization of Subtropical Agro-bioresources and College of Life Science and Technology, 100 Daxue Road, Nanning, Guangxi 530004, ChinaSearch for more papers by this author, Wei Jiang†Corresponding authors: Y.-Q. He; E-mail Address: yqhe@gxu.edu.cn, and W. Jiang; E-mail Address: weijiang@gxu.edu.cnState Key Laboratory for Conservation and Utilization of Subtropical Agro-bioresources and College of Life Science and Technology, 100 Daxue Road, Nanning, Guangxi 530004, ChinaSearch for more papers by this author, and Yong-Qiang He†Corresponding authors: Y.-Q. He; E-mail Address: yqhe@gxu.edu.cn, and W. Jiang; E-mail Address: weijiang@gxu.edu.cnhttps://orcid.org/0000-0001-6803-0431National Demonstration Center for Experimental Plant Science Education, College of Agriculture, Guangxi University, Nanning 530004, Guangxi, ChinaState Key Laboratory for Conservation and Utilization of Subtropical Agro-bioresources and College of Life Science and Technology, 100 Daxue Road, Nanning, Guangxi 530004, ChinaSearch for more papers by this author AffiliationsAuthors and Affiliations Yiming Li1 2 Xue Dan1 Tianjiao Li1 2 Qiaozhi Wu2 Zihao Wang2 Lidan Xu2 Xiuguo Ma1 Fengzhi Bie2 Shimin Lei1 Sheng Huang2 Wei Jiang2 † Yong-Qiang He1 2 † 1National Demonstration Center for Experimental Plant Science Education, College of Agriculture, Guangxi University, Nanning 530004, Guangxi, China 2State Key Laboratory for Conservation and Utilization of Subtropical Agro-bioresources and College of Life Science and Technology, 100 Daxue Road, Nanning, Guangxi 530004, China Published Online:8 Mar 2022https://doi.org/10.1094/PDIS-10-21-2176-AAboutSectionsView articlePDFSupplemental ToolsAdd to favoritesDownload CitationsTrack Citations ShareShare onFacebookTwitterLinked InRedditEmailWechat View articleGenome AnnouncementThe phytopathogenic bacterium Xanthomonas oryzae pv. oryzae is the causal agent of rice bacterial leaf blight (BB), one of the most destructive rice diseases in major rice-growing regions worldwide (Niño-Liu et al. 2006). The interactions between X. oryzae pv. oryzae and the rice plant follow the gene-for-gene hypothesis to some extent, making the X. oryzae pv. oryzae–rice pathosystem an excellent model platform in plant sciences (Zhang and Wang 2013). Among the identified X. oryzae pv. oryzae strains, PXO99A, also known as Philippine race 6, is widely acknowledged as the most representative type strain of X. oryzae pv. oryzae and important for differentiating X. oryzae pv. oryzae strains in laboratories and in rice disease resistance breeding. In 2008, Salzburg and his colleagues published the first version of the complete genome of PXO99A, which contains a unique, near-perfect tandem repeat of 212 kb close to the replication terminus (Salzburg et al. 2008). In 2015, Booher and his colleagues confirmed that there are two genomic structures in the progenies from the original PXO99A, with or without the 212-kb duplication, by using long-range sequencing and specific PCR verification (Booher et al. 2015). The genome of the PXO99A duplication version, designated the PXO99A-s accession, has widely been used as a benchmarking reference genome in relevant research around the world (Booher et al. 2015; Salzburg et al. 2008). In 2010, we obtained a substrain of PXO99A, named X. oryzae pv. oryzae PXO99A-GX (shortened to GX99), from the College of Plant Protection at Nanjing Agriculture University, the first institute which had imported the strain PXO99A used in China. In our mutational gene manipulations in PXO99A-GX, we found that the phenotypic changes in mutants of genes in the putative 212-kb duplication region had no significant differences from that of the genes in the other genomic regions (data not shown), which suggested that PXO99A-GX is not a PXO99A-s accession (Booher et al. 2015). In order to clarify the relationship between this substrain and the original strain PXO99A, we sequenced the complete genome of PXO99A-GX using PacBio and Illumina sequencing.X. oryzae pv. oryzae PXO99A-GX was grown in PSG medium (peptone at 10 g liter−1, sucrose at 10 g liter−1, and glutamic acid at 1 g liter−1) for 36 h, about to the early stationary phase. Genomic DNA was isolated using HiPure Bacterial DNA Kits (Magen, Guangzhou, China). The quality of the isolated genomic DNA was assessed using Qubit and Nanodrop accordingly and was fragmented with G-tubes. The fragment size (>10 kb) was selected by blue pippin system to prepare SMRTbell DNA template libraries. DNA sequencing was performed on the Pacific Biosciences Sequel (PacBio, Menlo Park, CA, U.S.A.) by Genedenovo Biotechnology Co., Ltd. (Guangzhou, China). In all, 145,826 reads and 1,162,120,421 bp were obtained and the N50 length was 9,258 bp. Sequencing depth of PXO99A-GX was around 230×. Continuous long reads obtained from SMRT sequencing were assembled de novo using Hifiasm software (version 0.16.0) (Cheng et al. 2021). The reads from the Illumina platform were used to correct the genome sequences to improve the quality of the assembly by Pilon software (version 1.23) (Walker et al. 2014). A circular chromosome of 5,028,051 bp (Fig. 1) was attained with 63.65% G+C content. The genome completeness and contamination were evaluated by CheckM v1.1.3 (Parks et al. 2015). The genome annotation of PXO99A-GX was carried out by the NCBI prokaryotic genome annotation pipeline (Tatusova et al. 2016) and GeneMarkS-2 software (Lomsadze et al. 2018). In total, 4,493 protein coding genes, 6 ribosomal RNAs, 53 transfer RNAs, and 147 soluble RNAs were predicted by the pipeline (Table 1).Fig. 1. Comparison of genomic structures and transcription activator-like (tal) gene clusters between Xanthomonas oryzae pv. oryzae strains PXO99A and PXO99A-GX. Block diagram showing the syntenic relationships between two genomes. The white sector is the 210-kb duplication present in PXO99A but absent in PXO99A-GX. The tal gene clusters were mapped to the genomes according to their relative positions on each chromosome.Download as PowerPointTable 1. Comparison of the genomic features of Xanthomonas oryzae pv. oryzae strains PXO99A and PXO99A-GXGenome featuresPXO99APXO99A-GXPXO86IXO1088AccessionyNC_010717.2CP073820.1NZ_CP007166.1NZ_CP040687.1Genome size (bp)5,238,5555,028,0515,016,6235,093,052G+C content (%)63.6063.6563.7263.67Coding density (%)86.7086.4586.7986.49Genes4,8774,6994,6724,758 Coding sequences (CDSs)4,6664,4934,4704,554 Transfer RNA54535353 Ribosomal RNA6666 Noncoding RNA151147143145 TALE genesz19171819 Average CDS length (bp)957953957950 Pseudogenes759737720745yGenome data of the other X. oryzae pv. oryzae strains in this table were obtained from NCBI (https://www.ncbi.nlm.nih.gov/genome/).zTALE = transcription activator-like effector.Table 1. Comparison of the genomic features of Xanthomonas oryzae pv. oryzae strains PXO99A and PXO99A-GXView as image HTML Comparative genomic analysis between X. oryzae pv. oryzae strains PXO99A-s and PXO99A-GX showed that both chromosome structures shared good collinearity, seemingly identical if the PXO99A-s duplication was removed from the genome of PXO99A (Booher et al. 2015). In total, there are 43 sites of genomic variations, including 8 single-nucleotide polymorphisms and 35 insertion-deletions (Supplementary Table S1), which resulted in 210,504 bp and 178 fewer genes in PXO99A-GX than in PXO99 A (Table 1). Of the genomic variation sites, 16 were found in the intergenic regions and 25 in the protein coding genes. Most of the annotated genes in PXO99A-GX had identical allelic counterparts in PXO99A, and even the pseudogenes in PXO99A-GX mostly kept the same mutated sites as those in PXO99A. Transcription activator-like effectors (TALEs) are the most important virulence factors in X. oryzae pv. oryzae, which act as transcription factors in the plant cells by binding to certain DNA targets (Moscou and Bogdanove 2009). The tal genes for TALEs consist of highly repetitive DNA sequences, which are presumed to be highly dynamic, involving recombination and deletions (Erkes et al. 2017). To our surprise, PXO99A-GX encodes 17 tal genes which are identical to the orthologs in PXO99A, except for the 2 tal genes in the duplication region, with no one alteration in the repeat variable diresidue sequence in the reciprocal TALEs (Grau et al. 2016) (Table 2). Because the sequenced sample of PXO99A-GX has been kept at −80°C in an ultralow temperature refrigerator, it is tempting to attribute the genome stability to the infrequent subculturing, which posed few stresses or selection pressures on PXO99A-GX. Although the genome sizes of the two genomes are significantly distinct, PXO99A has only three strain-specific genes, including PXO_RS12260, which encodes an atypical TALE (Supplementary Tables S1 and S2); and PXO99A-GX has one unique gene (i.e., KFK26_14970), which encodes a dipeptidylpeptidase IV N-terminal domain-containing protein (Supplementary Table S1). Assays on the impacts of the genomic variations on the biological traits in PXO99A-GX are in process.Table 2. Repeat variable diresidue (RVD) sequences in transcription activator-like effectors (TALEs) of Xanthomonas oryzae pv. oryzae strains PXO99A-GX and PXO99AClassaPXO99A-GXPXO99AFunctionsRVD sequencesTalBMtal1tal1pthXo7NI NG NI NI N* NN HD HD N* NI NI NI NG HD HG NN NS NN HD HD NG N*TalBXtal2btal2bpthXo1NN HD NI HG HD NG N* HD HD NI NG NG NI HD NG NN NG NI NI NI NI N* NS N*TalAGtal2atal2a−NI NG NN NG NK NG NI NN NI NN NI HD N* NS N*TalAWtal3atal3aiTalNS HD NG NG NG HD HD NG HD NN NG HD NN HD NG HD NI N*TalAItal3btal3biTalNS HD NG NG NG NG NG HD HD HD NN HD NG HD HD HD HD N*TalAFtal4tal4−NI NN NN NI NI NI HD NS HG NN NN NN NI NI HG HDTalARtal5btal5b−NI H* NI NN NN NN NN NN HD NI HD HG HD NI N* NS NI NI HG HD NS NS NGTalBAtal5atal5a−NI NS HD HG NS NN HD H* NG NN NN HD HD NG HD NGTalAHtal6atal6apthXo6NI N* NI NS NN NG NN NS N* NS NN NS N* NI HG HD NI HD HD NGTalANtal6btal6b−NI HG NI HG NI NI NI HD NN HD NS NG SS HD NI NI NN NI NN NI NGTalAB−tal7a−NI HG NI NI NI NN HD NS NN NS NN HD NN NI HD NN NS NG−tal7b−NI HG NI NN NGTalABtal7atal8a−NI HG NI NI NI NN HD NS NN NS NN HD NN NI HD NN NS NGTalAAtal7btal8b−NI HG NS HG HG HD NS NG HD NN NG HG NG HD HG HD HD NI NN NGTalAPtal8atal9apthXo8HD HD HD NG N* NN HD HD N* NI NI NN HD HI ND HD NI HD NG NGTalAQtal8btal9bavrXa23HD HD NN NN NG NG HD NS HG HD NG N* HD HD HD N* NN NI NN HD HI ND HD HG NN HG N*TalAOtal8ctal9cavrXa27NI NN N* NG NS NN NN NN NI NN NI N* HD HD NI NG NGTalAEtal8dtal9d−NI NN NI HG HG NV HG HD HG HD HD HD NGTalADtal8etal9e−NN HD NS NG HD NN N* NI HD NS HD NN HD NN HD NN NN NN NN NN NN NN HD NGaTALE classes are based on AnnoTALE (Grau et al. 2016); − indicates no data available.Table 2. Repeat variable diresidue (RVD) sequences in transcription activator-like effectors (TALEs) of Xanthomonas oryzae pv. oryzae strains PXO99A-GX and PXO99AView as image HTML In this sequencing project, we have sequenced the complete genome of a substrain of PXO99A—namely, PXO99A-GX—which has been used as the lab strain in many laboratories in China. We speculated that PXO99A-GX is similar to PXO99A-B or PXO99A-L-’13 (Booher et al. 2015), representing a type of genome state of PXO99A derivatives. In addition to the genome of PXO99A-s, the genome of PXO99A-GX might serve as an alternative high-quality reference genome resource to facilitate studies of the X. oryzae pv. oryzae–rice pathosystem.This whole-genome PacBio sequencing project has been deposited in NCBI under the accession number CP073820, with the raw sequencing data at BioProject PRJNA725623 and Sequence Read Archive SRR15390483 and SRR15390484.The author(s) declare no conflict of interest.Literature CitedBooher, N. J., Carpenter, S. C. D., Sebra, R. P., Wang, L., Salzberg, S. L., Leach, J. E., and Bogdanove, A. J. 2015. Single molecule real-time sequencing of Xanthomonas oryzae genomes reveals a dynamic structure and complex TAL (transcription activator-like) effector gene relationships. Microb. Genomics 1:e000032. https://doi.org/10.1099/mgen.0.000032 Crossref, ISI, Google ScholarCheng, H., Concepcion, G. T., Feng, X., Zhang, H., and Li, H. 2021. Haplotype-resolved de novo assembly using phased assembly graphs with hifiasm. Nat. Methods 18:170-175. https://doi.org/10.1038/s41592-020-01056-5 Crossref, ISI, Google ScholarErkes, A., Reschke, M., Boch, J., and Grau, J. 2017. Evolution of transcription activator-like effectors in Xanthomonas oryzae. Genome Biol. Evol. 9:1599-1615. https://doi.org/10.1093/gbe/evx108 Crossref, ISI, Google ScholarGrau, J., Reschke, M., Erkes, A., Streubel, J., Morgan, R. D., Wilson, G. G., Koebnik, R., and Boch, J. 2016. AnnoTALE: Bioinformatics tools for identification, annotation, and nomenclature of TALEs from Xanthomonas genomic sequences. Sci. Rep. 6:21077. https://doi.org/10.1038/srep21077 Crossref, ISI, Google ScholarLomsadze, A., Gemayel, K., Tang, S., and Borodovsky, M. 2018. Modeling leaderless transcription and atypical genes results in more accurate gene prediction in prokaryotes. Genome Res. 28:1079-1089. https://doi.org/10.1101/gr.230615.117 Crossref, ISI, Google ScholarMoscou, M. J., and Bogdanove, A. J. 2009. A simple cipher governs DNA recognition by TAL effectors. Science 326:1501. https://doi.org/10.1126/science.1178817 Crossref, ISI, Google ScholarNiño-Liu, D. O., Ronald, P. C., and Bogdanove, A. J. 2006. Xanthomonas oryzae pathovars: Model pathogens of a model crop. Mol. Plant Pathol. 7:303-324. https://doi.org/10.1111/j.1364-3703.2006.00344.x Crossref, ISI, Google ScholarParks, D. H., Imelfort, M., Skennerton, C. T., Hugenholtz, P., and Tyson, G. W. 2015. CheckM: Assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 25:1043-1055. https://doi.org/10.1101/gr.186072.114 Crossref, ISI, Google ScholarSalzberg, S. L., Sommer, D. D., Schatz, M. C., Phillippy, A. M., Rabinowicz, P. D., Tsuge, S., Furutani, A., Ochiai, H., Delcher, A. L., Kelley, D., Madupu, R., Puiu, D., Radune, D., Shumway, M., Trapnell, C., Aparna, G., Jha, G., Pandey, A., Patil, P. B., Ishihara, H., Meyer, D. F., Szurek, B., Verdier, V., Koebnik, R., Dow, J. M., Ryan, R. P., Hirata, H., Tsuyumu, S., Lee, S. W., Ronald, P. C., Sonti, R. V., Van Sluys, M.-A., Leach, J. E., White, F. F., and Bogdanove, A. J. 2008. Genome sequence and rapid evolution of the rice pathogen Xanthomonas oryzae pv. oryzae PXO99A. BMC Genomics 9:204. https://doi.org/10.1186/1471-2164-9-204 Crossref, ISI, Google ScholarTatusova, T., DiCuccio, M., Badretdin, A., Chetvernin, V., Nawrocki, E. P., Zaslavsky, L., Lomsadze, A., Pruitt, K. D., Borodovsky, M., and Ostell, J. 2016. NCBI prokaryotic genome annotation pipeline. Nucleic Acids Res. 44:6614-6624. https://doi.org/10.1093/nar/gkw569 Crossref, ISI, Google ScholarWalker, B. J., Abeel, T., Shea, T., Priest, M., Abouelliel, A., Sakthikumar, S., Cuomo, C. A., Zeng, Q., Wortman, J., Young, S. K., and Earl, A. M. 2014. Pilon: An integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS One 9:e112963. https://doi.org/10.1371/journal.pone.0112963 Crossref, ISI, Google ScholarZhang, H., and Wang, S. 2013. Rice versus Xanthomonas oryzae pv. oryzae: A unique pathosystem. Curr. Opin. Plant Biol. 16:188-195. https://doi.org/10.1016/j.pbi.2013.02.008 Crossref, ISI, Google ScholarY. Li and X. Dan contributed equally to this work.Funding: The work is supported by the National Key R&D Program of China (2018YFD0200302), the Natural Science Foundation of China (32060600), and the Natural Science Foundation of Guangxi (2019GXNSFAA245052).The author(s) declare no conflict of interest.DetailsFiguresLiterature CitedRelated Vol. 106, No. 3 March 2022SubscribeISSN:0191-2917e-ISSN:1943-7692 Download Metrics Downloaded 433 times Article History Issue Date: 30 Mar 2022Published: 8 Mar 2022Accepted: 7 Dec 2021 Pages: 1045-1048 Information© 2022 The American Phytopathological SocietyFundingNational Key R&D Program of ChinaGrant/Award Number: 2018YFD0200302Natural Science Foundation of ChinaGrant/Award Number: 32060600Natural Science Foundation of GuangxiGrant/Award Number: 2019GXNSFAA245052KeywordsgenomePacBio sequencingPXO99(^A)substrainXanthomonas oryzae pv. oryzaeThe author(s) declare no conflict of interest.PDF download