search for




 

Draft genome sequence of Weizmannia coagulans BA375 isolated from saltern soil in Korea
Korean J. Microbiol. 2021;57(2):119-121
Published online June 30, 2021
© 2021 The Microbiological Society of Korea.

Ki-Hong Yoon*

Food Science & Biotechnology Major, Woosong University, Daejeon 34606, Republic of Korea
Correspondence to: E-mail: ykh@wsu.ac.kr; Tel.: +82-42-630-9742; Fax: +82-42-630-9389
Received April 21, 2021; Revised June 1, 2021; Accepted June 1, 2021.
Abstract
The strain Weizmannia coagulans BA375, which was isolated from soil of a saltern, produced 1,4-β-xylanase and acetylxylan esterase. The draft genome sequence of the strain W. coagulans BA375 was determined using Illumina MiSeq platform. The genome comprises 82 contigs with genome size of 3,060,890 bp long and GC content of 46.9%. The draft genome contains 2,785 protein-coding genes, 130 pseudogenes and 117 RNA genes including 33 ribosomal RNA genes, 80 transfer RNA (tRNA) genes and 4 non-coding RNA (ncRNA) genes. The strain W. coagulans BA375 has polyprenyl glycosylphosphotransferase gene which is involved in biosynthesis of exopolysaccharides.
Keywords : Weizmannia coagulans, draft genome sequence, exopolysaccharide, Illumina MiSeq, xylanase
Body

Weizmannia coagulans is the type species of the genus Weizmannia which was recently transferred from the Coagulans clade consisting of three Bacillus species, B. coagulans, B. acidiproducens, and B. ginsengihumi in the genus Bacillus (Gupta et al., 2020). Weizmannia coagulans is facultatively aerobic, Gram-stain-positive, rod shaped, endospore-forming bacteria, and has been used as probiotics being more stable than Lactobacillus strains in heat treatment. Besides, it was known to produce lactic acid and food enzymes such as α-amylase, α-galactosidase, β-galactosidase and lipase (Konuray and Erginkaya, 2018). Several strains of W. coagulans have been isolated from canned food, tomato juice, gelatin, milk, pharmaceuticals and silage as well as soil (Gupta et al., 2020). Xylanase is widely used in food and pulp bleaching and several studies have been performed to produce xylanase by W. coagulans (Heck et al., 2005; Chauhan et al., 2006) and exopolysaccharide have found diverse applications in various food and pharmaceutical industries. In this study, we present the draft genome sequence and annotation of the strain W. coagulans BA375.

The strain was isolated by diluting a soil sample in sterile 0.85% (w/v) NaCl solution and plating on tryptic soy agar (TSA; BD) supplemented with 3% (w/v) NaCl and the adjustment of pH to 7.0 using 1 M NaOH. The strain was subsequently purified three times by plating on TSA supplemented with 3% (w/v) NaCl at 37°C for 2 days and maintained on the same medium.

Genomic DNA of strain BA375 was extracted using Genomic DNA Prep Kit (Biofact) and used to sequence the 16S rRNA gene, according to the previous study (Lee and Yoon, 2008). The comparison of the 16S rRNA gene sequence showed the highest similarity to Bacillus coagulans ATCC 7050T (98.9%). The whole-genome sequencing and analysis of the strain BA375 was performed using the paired-end sequencing method on the MiSeq platform (Illumina) at ChunLab. The estimation of genome completeness and quality was performed by CheckM v1.0.12 (Parks et al., 2015). The resulting short sequencing reads were assembled in contigs using the SPAdes 3.13.0 program (Bankevich et al., 2012). The genome sequence was annotated with the Rapid Annotation of microbial genomes using Subsystems Technology (RAST) (Aziz et al., 2008; Overbeek et al., 2014), and the protein-coding sequences were determined using the NCBI Prokaryotic Genome Annotation Pipeline (PGAP) (Tatusova et al., 2016). The genome sequences of the members of the species W. coagulans which are available in the NCBI database were taken for comparative analysis. The average nucleotide identity (ANI) values between genomes of the strain BA375 and other W. coagulans strains were calculated using the ANI Calculator (www.ezbiocloud.net/ tools/ani) (Yoon et al., 2017a, 2017b). To determine genomic relatedness, digital DDH (dDDH) values were computed using Genome-to-Genome Distance Calculator (GGDC) version 2.1 (Meier-Kolthoff et al., 2013).

The draft genome of W. coagulans BA375 composed of 82 contigs with a chromosome length of 3,060,890 bp (N50 value 90,601 bp) with 572× sequencing coverage, and 46.9% GC content. The result of CheckM estimation indicated that genome completeness at 100.0% with 0.0% contamination and strain heterogeneity. The NCBI PGAP annotation revealed that the genome contains 2,785 protein-coding genes, 130 pseudogenes and 117 RNA genes including 33 rRNA genes (10 5S ribosomal RNA, 11 16S ribosomal RNA, 12 23S ribosomal RNA), 80 transfer RNA (tRNA) genes and 4 non-coding RNA (ncRNA) genes (Table 1). Two conserved signature indels (CSIs), acetate kinase and O-methyltransferase, that are shared exclusively by the members of the genus Weizmannia (Gupta et al., 2020) were found in the genome of strain BA375. The average nucleotide identity (ANI) and dDDH value between genomes of strain BA375 and the closest type strain, W. coagulans ATCC 7050T was 97.9% and 80.6%, respectively. The PGAP analysis showed that the genome of the strain BA375 contains genes encoding 1,4-β-xylanase (J1C53_RS14485) and acetylxylan esterase (J1C53_RS14455). Moreover, we also found genes involved in the production of exopolysaccharides, polyprenyl glycosylphosphotransferase (J1C53_RS10835). Among genomes of 43 W. coagulans strains, only strain 36D1 had endo-1,4-β-xylanase (BCOA_RS21565), while the genes of acetylxylan esterase and polyprenyl glycosylphosphotransferase were not found in any other genomes except the strain BA375. These findings signify that W. coagulans BA375 has potential for applications in food and feed industries.

General genomic features of W. coagulans BA375

Features W. coagulans BA375
Genome length (bp) 3,060,890
GC-content (%) 46.94
Contigs 82
Total genes 2,915
Protein-coding genes 2,785
Pseudo-genes 130
rRNA genes (5S, 16S, 23S) 33 (10,11,12)a
tRNA genes 80
ncRNA genes 4

a rRNA genes including 10 (5S) complete rRNAs, and 11, 12 (16S, 23S) partial rRNAs.



Nucleotide sequence and strain accession numbers

The 16S rRNA gene sequence and draft genome sequence of Weizmannia coagulans BA375 has been deposited at GenBank under the accession LC140744 and JAFMQA010000000. The strain BA375 has been deposited at Korean Agricultural Culture Collection (KACC 18681) and at NITE Biological Resource Center in Japan (NBRC 112273).

적 요

염전 토양에서 분리한 Weizmannia coagulans BA375는 1,4-β-xylanase 및 acetylxylan esterase를 생산한다. BA375 균주의 전장유전체는 3,060,890 bp 크기로 82개의 contig로 구성되어 있으며 GC 함량은 46.9% 이다. 유전체는 2,785개의 단백질 코딩 유전자로 되어 있으며, 130개의 유사 유전자 및 33개의 리보솜 RNA 유전자, 80개의 운반 RNA (tRNA) 유전자 및 4개의 비 코딩 RNA 유전자 (ncRNA) 등 117개의 RNA 유전자를 포함하고 있다. BA375 균주는 균체외 다당체의 생합성에 관여하는 polyprenyl glycosylphosphotransferase를 가지고 있는 것으로 확인되었다.

Acknowledgments

This research is based on the support of 2020 Woosong University Academic Research Funding.

Conflict of Interest

The author declares that I have no conflicts of interests.

References
  1. Aziz RK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, Formsma K, Gerdes S, Glass EM, and Kubal MKubal M, et al. 2008. The RAST server: rapid annotations using subsystems technology. BMC Genomics 9, 75.
    Pubmed KoreaMed CrossRef
  2. Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, Lesin VM, Nikolenko SI, Pham S, and Prjibelski ADPrjibelski AD, et al. 2012. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19, 455-477.
    Pubmed KoreaMed CrossRef
  3. Chauhan S, Choudhury B, Singh SN, and Ghosh P. 2006. Application of xylanase enzyme of Bacillus coagulans as a prebleaching agent on non-woody pulps. Process Biochem. 41, 226-231.
    CrossRef
  4. Gupta RS, Patel S, Saini N, and Chen S. 2020. Robust demarcation of 17 distinct Bacillus species clades, proposed as novel Bacillaceae genera, by phylogenomics and comparative genomic analyses: description of Robertmurraya kyonggiensis sp. nov. and proposal for an emended genus Bacillus limiting it only to the members of the Subtilis and Cereus clades of species. Int. J. Syst. Evol. Microbiol. 70, 5753-5798.
    Pubmed CrossRef
  5. Heck JX, Flôres SH, Hertz PF, and Ayub MAZ. 2005. Optimization of cellulase-free xylanase activity produced by Bacillus coagulans BL69 in solid-state cultivation. Process Biochem. 40, 107-112.
    CrossRef
  6. Konuray G and Erginkaya Z. 2018. Potential use of Bacillus coagulans in the food industry. Foods 7, 92.
    Pubmed KoreaMed CrossRef
  7. Lee JC and Yoon KH. 2008. Paenibacillus woosongensis sp. nov., a xylanolytic bacterium isolated from forest soil. Int. J. Syst. Evol. Microbiol. 58, 612-616.
    Pubmed CrossRef
  8. Meier-Kolthoff JP, Auch AF, Klenk HP, and Göker M. 2013. Genome sequence-based species delimitation with confidence intervals and improved distance functions. BMC Bioinformatics 14, 60.
    Pubmed KoreaMed CrossRef
  9. Overbeek R, Olson R, Pusch GD, Olsen GJ, Davis JJ, Disz T, Edwards RA, Gerdes S, Parrello B, and Shukla MShukla M, et al. 2014. The SEED and the rapid annotation of microbial genomes using subsystems technology (RAST). Nucleic Acids Res. 42, D206-D214.
    Pubmed KoreaMed CrossRef
  10. Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, and Tyson GW. 2015. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 25, 1043-1055.
    Pubmed KoreaMed CrossRef
  11. Tatusova T, DiCuccio M, Badretdin A, Chetvernin V, Nawrocki EP, Zaslavsky L, Lomsadze A, Pruitt KD, Borodovsky M, and Ostell J. 2016. NCBI prokaryotic genome annotation pipeline. Nucleic Acids Res. 44, 6614-6624.
    Pubmed KoreaMed CrossRef
  12. Yoon SH, Ha SM, Kwon S, Lim J, Kim Y, Seo H, and Chun J. 2017a. Introducing EzBioCloud: A taxonomically united database of 16S rRNA gene sequences and whole-genome assemblies. Int. J. Syst. Evol. Microbiol. 67, 1613-1617.
    Pubmed KoreaMed CrossRef
  13. Yoon SH, Ha SM, Lim J, Kwon S, and Chun J. 2017b. A large-scale evaluation of algorithms to calculate average nucleotide identity. Antonie van Leeuwenhoek 110, 1281-1286.
    Pubmed CrossRef


June 2021, 57 (2)