- Short genome report
- Open Access
High-quality-draft genomic sequence of Paenibacillus ferrarius CY1T with the potential to bioremediate Cd, Cr and Se contamination
Standards in Genomic Sciencesvolume 12, Article number: 60 (2017)
Paenibacillus ferrarius CY1T (= KCTC 33419T = CCTCC AB2013369T) is a Gram-positive, aerobic, endospore-forming, motile and rod-shaped bacterium isolated from iron mineral soil. This bacterium reduces sulfate (SO4 2−) to S2−, which reacts with Cd(II) to generate precipitated CdS. It also reduces the toxic chromate [Cr(VI)] and selenite [Se(VI)] to the less bioavailable chromite [Cr(III)] and selenium (Se0), respectively. Thus, strain CY1T has the potential to bioremediate Cd, Cr and Se contamination, which is the main reason for the interest in sequencing its genome. Here we describe the features of strain CY1T, together with the draft genome sequence and its annotation. The 9,184,169 bp long genome exhibits a G + C content of 45.6%, 7909 protein-coding genes and 81 RNA genes. Nine putative Se(IV)-reducing genes, five putative Cr(VI) reductase and nine putative sulfate-reducing genes were identified in the genome.
The genus Paenibacillus was established in 1993 with Paenibacillus polymyxa as the type species [1, 2]. The common characteristics of the Paenibacillus members are aerobic, Gram-positive, rod-shaped and endospore-forming . Some Paenibacillus strains have the ability for plant growth promotion, biocontrol, manufacturing process and bioremediation, which making them very important in agricultural, industrial and medical applications . A variety of industrial wastes including crude oil, diesel fuel, textile dyes, aliphatic and aromatic organic pollutants could be degraded by Paenibacillus strains [5,6,7,8,9,10,11]. However, the bioremediation of heavy metal(loids) contamination by Paenibacillus strains are rarely reported.
Paenibacillus ferrarius CY1T is a multi-metal(loids) resistant bacterium isolated from iron mineral soil in Hunan Province, China . During cultivation, it could efficiently reduce sulfate (SO4 2−) to S2−, which could precipitate with cadmium [Cd(II)] to generate CdS . In addition, it also reduces the more toxic chromate [Cr(VI)] and selenite [Se(VI)] to the much less toxic chromite [Cr(III)] and selenium (Se0), respectively. Based on these interesting features, we propose that strain CY1T represents a promising candidate for bioremediation of Cd, Cr and Se contamination. To gain insight into the molecular mechanisms involved in sulfate/chromate/selenite reduction and metal(loids) resistance, and to enhance its biotechnological applications, we analyze the high quality draft genome of this bacterium.
Classification and features
P. ferrarius CY1T is a Gram-positive, endospore-forming, motile and aerobic bacterium. The rod-shaped cells are 0.5–0.8 mm in width and 4.2–5.7 mm in length with peritrichous flagella (Fig. 1). Colonies are yellowish to creamy-white, smooth and circular on NA agar plate . Growth occurs at temperature and pH range of 4–37 °C and pH 5.0–8.0, respectively . Optimal growth occurs at 28 °C and pH 6.0–7.0 (Table 1). Strain CY1T grows on NA/R2A/LB and TSA media, but cannot grow on MacConkey agar . The phylogenetic relationship of P. ferrarius CY1T with other members within the genus Paenibacillus is shown in a 16S rRNA based neighbor-joining tree, and strain CY1T is closely related to Paenibacillus marchantiophytorum R55 T (KP056549) (Fig. 2).
Physiological and biochemical analyses were performed using the API 20NE test (bioMérieux, France), ID 32GN text (bioMérieux, France) and traditional classification methods. Strain CY1T is positive for oxidase and catalase activities, hydrolysis of Tween 80 and aesculin and production of NH3 and H2S, but is negative for nitrate reduction, citrate utilization, egg yolk reaction, production of indole, and hydrolysis of starch, gelatin, casein, urea, L-tyrosine, arginine, Tween 20, DNA and CM-cellulose . The carbon sources, which can be used by strain CY1T, are shown in Table 1.
The resistance levels of P. ferrarius CY1T for multi-metal(loids) were tested with the minimal inhibition concentration on NA agar plates using Na3AsO3, K2Sb2(C4H2O6)2, Na2SeO3, K2CrO4, CdCl2, PbCl2, CuCl2 and MnCl2. The results showed that the MICs for As(III), Sb(III), Se(IV), Cr(VI), Cd(II), Pb(II), Cu(II) and Mn(II) are 2, 1, 8, 4, 0.08, 1, 0.5 and 100 mmol/L, respectively. In addition, the abilities of strain CY1T for Cd(II) removal, and Cr(VI) and Se(IV) reduction were tested. Strain CY1T was incubated in LB medium for Cd(II) removal and in NA medium for Cr(VI) and Se(IV) reduction, since NA medium can absorb some of the Cd(II). When OD600 reach 0.6-0.7, CdCl2 (50 μmol/L), K2CrO4 (200 μmol/L) and Na2SeO3 (200 μmol/L) were each added to the culture. At designated times, culture samples were taken for measuring the residual concentrations of Cd(II), Cr(VI) and Se(IV). The concentration of Cd(II) was measured by the atomic absorption spectrometry . The concentration of Cr(VI) was measured by the UV spectrophotometer (DU800, Beckman, CA, USA) with the colorimetric diphenylcarbazide method , and the concentration of Se(IV) was tested by HPLC-HG-AFS (Beijing Titan Instruments Co., Ltd., China) . The results showed that strain CY1T could remove nearly 50 μmol/L Cd(II) in 72 h (Fig. 3a) and reduce 200 μmol/L Cr(VI) and Se(IV) in 5 h and 6 h, respectively (Fig. 3b, c). The removed Cd(II) is presented as pellets that is most probably by the reaction of Cd(II) with H2S to produce precipitated CdS.
Genome sequencing information
Genome project history
Strain CY1T was selected for genome sequencing on the basis of its ability for Cd(II) removal, Cr(VI) and Se(IV) reduction, these characters made strain CY1T with great value for genetic study and for bioremediation of Cd, Cr and Se contamination. The draft genome sequence is deposited at DDBJ/EMBL/GenBank under the accession number MBTG00000000. The final genome consists of 73 scaffolds with 289.77 × coverage. A summary of the project information is shown in Table 2.
Growth conditions and genomic DNA preparation
Overnight cultures of strain CY1T was inoculated into 50 mL of NA medium at 28 °C with 120 rpm shaking. After incubation for 36 h, the bacterial cells were harvested through centrifugation (13,400×g for 5 min at 4 °C). Genomic DNA was extracted using the QiAamp kit (Qiagen, Germany). The quality and quantity of the DNA were determined by a spectrophotometer (NanoDrop 2000, Thermo). Then, 10 μg of DNA was sent to Bio-broad Technology Co., Ltd., Wuhan, China for sequencing.
Genome sequencing and assembly
Genome sequencing and assembly were performed by Bio-broad Technology Co., Ltd., Wuhan, China, and all original sequence data can be found at the NCBI Sequence Read Archive. An Illumina standard shotgun library was constructed and sequenced using an Illumina Hiseq2000 platform with pair-end sequencing strategy (300 bp insert size) . The following quality control steps were performed for removing low quality reads: 1) removed the adapter sequences of reads; 2) trimmed the ambiguous bases (N) in 5′ end and the reads with a quality score lower than 20; and 3) filtered the reads which contain N more than 10% or have the length less than 50 bp (without adapters and N in 5′ end). The assembly of CY1T genome is based on 20,189,278 quality reads totaling 3,000,798,615 bp, which provides a coverage of 289.77×. Subsequently, the reads were assembled into 75 contigs (> 200 bp) using SOAPdenovo v2.04 , and the gaps between the contigs were closed by GapCloser v1.12 .
The draft genome of strain CY1T was annotated through the RAST server version 2.0 and the NCBI Prokaryotic Genome Annotation Pipeline. Genes were identified using the gene caller GeneMarkS+ with the similarity-based gene detection approach . Pseudogenes were also predicted using the NCBI PGAP. Internal gene clustering was performed by OrthoMCL using Match cutoff of 50% and E-value Exponent cutoff of 1-e5 [21, 22]. The COGs functional categories were assigned by WebMGA server  with E-value cutoff of 1-e10. The translations of the predicted CDSs were used to search against the Pfam protein family database  and the KEGG database . The transmembrane helices and signal peptides were predicted by TMHMM v. 2.0  and SignalP 4.1 , respectively.
The whole genome of strain CY1T reveals a genome size of 9,184,169 bp and a G + C content of 45.6% (Table 3). The genome contains 8260 coding sequences, 19 rRNA, 58 tRNA, and 4 ncRNA. Among 7909 protein-coding genes, 4231 were assigned as putative function, while the other 3678 were designated as hypothetical proteins. In addition, 6632 genes were categorized into COGs functional groups. Information about the genome statistics is shown in Table 3 and the classification of genes into COGs functional categories is summarized in Table 4.
Insights from the genome sequence
P. ferrarius CY1T is a multi-metal(loids) resistant bacterium with the capability of SO4 2−, Cr(VI) and Se(IV) reduction, suggesting that it has developed a number of evolutionary strategies to adapt to heavy metal (or metalloids) contaminated environments. To identify pathways and enzymes involved in SO4 2−, Cr(VI) and Se(IV) reduction, high quality draft genome sequence of strain CY1T was generated. The map of the P. ferrarius CY1T genome is shown in Fig. 4.
KEGG analysis showed that strain CY1T contains a complete SO4 2− reduction pathway, which is consistent with the phenotype of H2S production. The genes responsible for SO4 2− reduction include sulfate ABC transporter CysPWA, sulfate adenylyltransferase CysD, adenylylsulfate kinase CysC, adenylylsulfate reductase CysH and sulfite reductase CysJI (Table 5). The S2− generated from SO4 2− reduction could react with Cd(II) to form the participated CdS , which may contribute to the Cd(II) removal. For Cr(VI) reduction, five NADPH-dependent FMN reductase which have the same conserved domain as the Cr(VI) reductases ChrR (from Pseudomonas putida ) and YieF (from Escherichia coli ) , were identified in the genome of strain CY1T (Table 5). It has been reported that thioredoxin reductase ThxR and NADH:flavin oxidoreductase could reduce Se(IV) in Pseudomonas seleniipraecipitans and Rhizobium selenitireducens , respectively [29,30,31]. According to the NCBI and RAST annotation, seven thioredoxin reductases and two NADH-dependent flavin oxidoreductases were found in the genome of strain CY1T (Table 5), and some of these proteins may responsible for Se(IV) reduction in strain CY1T.
Strain CY1T could tolerant multi-metal(loids), such as As(III), Sb(III), Cr(VI), Cd(II), Pb(II), Cu(II) and Mn(II). Expectably, various metal resistant genes were identified in its genome (Table 6). Several transporters were found to responsible for the efflux of these metal(loids). In addition, the transcriptional regulator ArsR and arsenite reductase ArsC were also found to be involved in the As(III)/Sb(III) resistance (Table 6) [32,33,34]. Recently, it has been reported that an oxidoreductase AnoA, which belongs to the short-chain dehydrogenase/reductase family, and catalase KatA, which is responsible for H2O2 degradation, are all involved in bacterial Sb(III) oxidation/resistance in Agrobacterium tumefaciens GW4 [35,36,37,38]. One AnoA homologue oxidoreductase gene and five catalase genes were identified in the genome of strain CY1T (Table 6), which may associate with Sb(III) oxidation/resistance.
The genome of P. ferrarius CY1T harbors various genes responsible for sulfate transport and reduction, chromate and selenite reduction and resistance of multi-metal(loids), which is consistent with its phenotypes. To date, the utilization of Paenibacillus species in immobilization of heavy-metals (or metalloids) is still limited and the genes and enzymes involves in Cr(VI) and Se(IV) reduction were poorly understood in Paenibacillus members. The genomic sequence of strain CY1T enriches the genome information of Paenibacillus strains. More importantly, the genome information provides basis for understanding molecular mechanisms of microbial redox transformations of metal(loids).
Ash C, Priest FG, Collins MD. Molecular identification of rRNA group 3 bacilli (Ash, Farrow, Wallbanks and Collins) using a PCR probe test. Proposal for the creation of a new genus Paenibacillus. Antonie Van Leeuwenhoek. 1993;64:253–60.
Ash C, Priest FG, Collins MD. Paenibacillus gen. nov. and Paenibacillus polymyxa comb. nov. In validation of the publication of new names and new combinations previously effectively published outside the IJSB, List no. 51. Int J Syst Bacteriol. 1994;44:852.
Zhou Y, Gao S, Wei DQ, Yang LL, Huang X, He J, et al. Paenibacillus thermophilus sp. nov., a novel bacterium isolated from a sediment of hot spring in Fujian province, China. Antonie Van Leeuwenhoek. 2012;102:601–9.
Grady EN, MacDonald J, Liu L, Richman A, Yuan ZC. Current knowledge and perspectives of Paenibacillus: a review. Microb Cell Factories. 2016;15:203.
Li O, Lu C, Liu A, Zhu L, Wang PM, Qian CD, et al. Optimization and characterization of polysaccharide-based bioflocculant produced by Paenibacillus elgii B69 and its application in wastewater treatment. Bioresour Technol. 2013;134:87–93.
Abbasian F, Lockington R, Mallavarapu M, Naidu RA. comprehensive review of aliphatic hydrocarbon biodegradation by bacteria. Appl Biochem Biotechnol. 2015;176:670–99.
Haritash A, Kaushik C. Biodegradation aspects of polycyclic aromatic hydrocarbons (PAHs): a review. J Hazard Mater. 2009;169:1–15.
Spadaro JT, Isabelle L, Renganathan V. Hydroxyl radical mediated degradation of azo dyes: evidence for benzene generation. Environ Sci Technol. 1994;28:1389–93.
Moosvi S, Kher X, Madamwar D. Isolation, characterization and decolorization of textile dyes by a mixed bacterial consortium JW-2. Dyes Pigments. 2007;74:723–9.
Choi K, Park C, Kim S, Lyoo W, Lee SH, Lee J. Polyvinyl alcohol degradation by Microbacterium barkeri KCCM 10507 and Paenibacillus amylolyticus KCCM 10508 in dyeing wastewater. J Microbiol Biotechnol. 2004;14:1009–13.
Zheng B, Zhang F, Dong H, Chai L, Shu F, Yi S, et al. Draft genome sequence of paenibacillus sp. Strain A2. Stand Genomic Sci. 2016;11:9.
Cao Y, Chen F, Li Y, Wei S, Wang G. Paenibacillus ferrarius sp. Nov., isolated from iron mineral soil. Int J Syst Evol Microbiol. 2015;65:165–70.
Pagnanelli F, Cruz Viggi C, Toro L. Isolation and quantification of cadmium removal mechanisms in batch reactors inoculated by sulphate reducing bacteria: biosorption versus bioprecipitation. Bioresour Technol. 2010;101:2981–7.
Liao S, Zhou J, Wang H, Chen X, Wang H, Wang G. Arsenite oxidation using biogenic manganese oxides produced by a deep-sea manganese-oxidizing bacterium, Marinobacter sp. MnI7-9. Geomicrobiol J. 2013;30(2):150–9.
Monteiro MI, Fraga IC, Yallouz AV, de Oliveira NM, Ribeiro SH. Determination of total chromium traces in tannery effluents by electrothermal atomic absorption spectrometry, flame atomic absorption spectrometry and UV-visible spectrophotometric methods. Talanta. 2002;58(4):629–33.
Zheng S, Su J, Wang L, Yao R, Wang D, Deng Y, et al. Selenite reduction by the obligate aerobic bacterium Comamonas testosteroni S44 isolated from a metal-contaminated soil. BMC Microbiol. 2014;14:204.
Bennett S. Solexa Ltd. Pharmacogenomics. 2004;5:433–8.
Li R, Li Y, Kristiansen K, Wang JSOAP. short oligonucleotide alignment program. Bioinformatics. 2008;24:713–4.
Li R, Zhu H, Ruan J, Qian W, Fang X, Shi Z, et al. De novo assembly of human genomes with massively parallel short read sequencing. Genome Res. 2010;20(2):265–72.
Besemer J, Lomsadze A, Borodovsky M. GeneMarkS: a self-training method for prediction of gene starts in microbial genomes. Implications for finding sequence motifs in regulatory regions. Nucleic Acids Res. 2001;29:2607–18.
Li L, Stoeckert CJ Jr, Roos DS. OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 2003;13:2178–89.
Fischer S, Brunk BP, Chen F, Gao X, Harb OS, Iodice JB, et al. Using OrthoMCL to assign proteins to OrthoMCL-DB groups or to cluster proteomes into new ortholog groups. Curr Protoc Bioinformatics. 2011;6:12–9.
Wu S, Zhu Z, Fu L, Niu B, Li W. WebMGA: a customizable web server for fast metagenomic sequence analysis. BMC Genomics. 2011;12:444.
Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, et al. Pfam: the protein families database. Nucleic Acids Res. 2014;42:222–30.
Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M. The KEGG resource for deciphering the genome. Nucleic Acids Res. 2004;32:277–80.
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:567–80.
Petersen TN, Brunak S, Heijne GV, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8:785–6.
Ackerley DF, Gonzalez CF, Park CH, Blake R, Keyhan M, Matin A. Chromatereducing properties of soluble flavoproteins from Pseudomonas putida and Escherichia coli. Appl Environ Microbiol. 2004;70:873–82.
Bjornstedt M, Kumar S, Bjorkhem L, Spyrou G, Holmgren A. Selenium and the thioredoxin and glutaredoxin systems. Biomed Environ Sci. 1997;10:271–9.
Tamura T, Sato K, Komori K, Imai T, Kuwahara M, Okugichi T, et al. Selenite reduction by the thioredoxin sysem: kinetics and identification of protein-bound selenide. Biosci Biotechnol Biochem. 2011;75:118–7.
Hunter WJA. Rhizobium selenitireducens protein showing selenite reductase activity. Curr Microbiol. 2014;68:311–6.
Xu C, Zhou T, Kuroda M, Rosen BP. Metalloid resistance mechanisms in prokaryotes. J Biochem. 1998;23:16–23.
Martin P, DeMel S, Shi J, Gladysheva T, Gatti DL, Rosen BP, et al. Insights into the structure, solvation, and mechanism of ArsC arsenate reductase, a novel arsenic detoxification enzyme. Structure. 2001;9:1071–81.
Suzuki K, Wakao N, Kimura T, Sakka K, Ohmiya K. Expression and regulation of the arsenic resistance operon of Acidiphilium multivorum AIU 301 plasmid pKW301 in Escherichia coli. Appl Environ Microbiol. 1998;64:411–8.
Li J, Wang Q, Li M, Yang B, Shi M, Guo W, et al. Proteomics and genetics for identification of a bacterial antimonite oxidase in Agrobacterium tumefaciens. Environ Sci Technol. 2015;49(10):5980–9.
Li J, Wang Q, Oremland RS, Kulp TR, Rensing C, Wang G. Microbial antimony biogeochemistry - enzymes, regulation and related metabolic pathways. Appl Environ Microbiol. 2016;82(18):5482–95.
Li J, Yang B, Shi M, Yuan K, Guo W, Wang Q, et al. Abiotic and biotic factors responsible for antimonite oxidation in Agrobacterium tumefaciens GW4. Sci Rep. 2017;7:43225.
Li J, Yang B, Shi M, Yuan K, Guo W, Li M, et al. Effects upon metabolic pathways and energy production by Sb(III) and As(III)/Sb(III)-oxidase gene aioA in Agrobacterium tumefaciens GW4. PLoS One. 2017;12(2):e0172823.
Woese CR, Kandler O, Wheelis ML. Towards a natural system of organisms: proposal for the domains Archaea, Bacteria, and Eucarya. Proc Natl Acad Sci U S A. 1990;87:4576–9.
Gibbons N, Murray R. Proposals concerning the higher taxa of bacteria. Int J Syst Bacteriol. 1978;28:1–6.
Garrity GM, Holt JG. The road map to the manual. In: Bergey’s Manual® of Systematic Bacteriology. New York: Springer; 2001. p. 119–66.
Murray R. The higher taxa, or, a place for everything. Bergey’s Man Syst Bacteriol. 1984;1:31–4.
Ludwig WSK, Whitman WB. Class I. Bacilli class nov. In: De Vos P, Garrity G, Jones D, Krieg NR, Ludwig W, Rainey FA, Schleifer KH, Whitman WB, editors. Bergey’s manual of systematic bacteriology, vol. 3. 2nd ed. New York: Springer; 2009. p. 19–20.
Skerman VBD, Mcgowan V, Sneath PHA. Approved lists of bacterial names. Int J Syst Bacteriol. 1980;30:255–420.
Hauduroy P, Ehringer G. Dictionnaire des bactéries pathogènes. Paris: Masson; 1953.
Euzeby J. List of new names and new combinations previously effectively, but not validly, published. Int J Syst Evol Microbiol. 2006;56:925–7.
Judicial Commission of the International Committee for Systematics of P. The type species of the genus Paenibacillus Ash et al. 1994 is Paenibacillus polymyxa. Opinion 77. Int J Syst Evol Microbiol. 2005;55:513.
Validation List no. 51. Validation of the publication of new names and new combinations previ-ously effectively published outside the IJSB. Int J Syst Bacteriol. 1994;44:852.
Shida O, Takagi H, Kadowaki K, Nakamura LK, Komagata K. Transfer of Bacillus alginolyticus, Bacillus chondroitinus, Bacillus curdlanolyticus, Bacillus glucanolyticus, Bacillus kobensis, and Bacillus thiaminolyticus to the genus Paenibacillus and emended description of the genus Paenibacillus. Int J Syst Bacteriol. 1997;47:289–98.
Behrendt U, Schumann P, Stieglmeier M, Pukall R, Augustin J, Sproer C, et al. Characterization of heterotrophic nitrifying bacteria with respiratory ammonification and denitrification activity–description of Paenibacillus uliginis sp. nov., an inhabitant of fen peat soil and Paenibacillus purispatii sp. nov., isolated from a spacecraft assembly clean room. Syst Appl Microbiol. 2010;33:328–36.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25:25–9.
We thank Mr. Xian Xia and Dr. Jing Huang for technical assistance. This study was supported by National key research and development program of China (2016YFD0800702).
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.