Xueqin Ran
Zhongmei ZhuHong LongQun TianLongjiang YouXingdiao WuQin Liu
Shihui Huang
Sheng Li
Xi Niu
Jiafu Wang*- Key Laboratory of Plant Resource Conservation and Germplasm Innovation in Mountainous Region (Ministry of Education), College of Animal Science/Institute of Agro-Bioengineering, Guizhou University, Guiyang, China
The mechanism of bacterial adaption to manganese-polluted environments was explored using 50 manganese-tolerant strains of bacteria isolated from soil of the largest manganese mine in China. Efficiency of manganese removal by the isolated strains was investigated using atomic absorption spectrophotometry. Bacillus safensis strain ST7 was the most effective manganese-oxidizing bacteria among the tested isolates, achieving up to 82% removal at a Mn(II) concentration of 2,200 mg/L. Bacteria-mediated manganese oxide precipitates and high motility were observed, and the growth of strain ST7 was inhibited while its biofilm formation was promoted by the presence of Mn(II). In addition, strain ST7 could grow in the presence of high concentrations of Al(III), Cr(VI), and Fe(III). Genome-wide analysis of the gene expression profile of strain ST7 using the RNA-seq method revealed that 2,580 genes were differently expressed under Mn(II) exposure, and there were more downregulated genes (n = 2,021) than upregulated genes (n = 559) induced by Mn stress. KAAS analysis indicated that these differently expressed genes were mainly enriched in material metabolisms, cellular processes, organism systems, and genetic and environmental information processing pathways. A total of twenty-six genes from the transcriptome of strain ST7 were involved in lignocellulosic degradation. Furthermore, after 15 genes were knocked out by homologous recombination technology, it was observed that the transporters, multicopper oxidase, and proteins involved in sporulation and flagellogenesis contributed to the removal of Mn(II) in strain ST7. In summary, B. safensis ST7 adapted to Mn exposure by changing its metabolism, upregulating cation transporters, inhibiting sporulation and flagellogenesis, and activating an alternative stress-related sigB pathway. This bacterial strain could potentially be used to restore soil polluted by multiple heavy metals and is a candidate to support the consolidated bioprocessing community.
Introduction
Manganese (Mn), one of the most abundant metals in the Earth’s crust, exerts a significant impact on the biogeochemical cycles of other metals, sulfur, and carbon and is considered the strongest naturally occurring oxidizing agent in the environment (Tebo et al., 2005). Furthermore, Mn is an essential trace element for all living organisms, acting as a cofactor for many enzymes, such as oxidoreductases, transferases, ligases, and hydrolases (Wang et al., 2019). Mn also participates in oxidative stress resistance, metabolism and decomposition of reactive oxygen species (ROS) (Bosma et al., 2021), and is required for the development and function of nerve and immune cells, control of blood sugar and vitamin levels in animals, and photosynthesis and respiration processes in plants (Gao et al., 2014). However, consumption of food or drinking water containing a high level of Mn has undesirable effects on human health (O’Neal and Zheng, 2015). Excessive Mn can accumulate in the human pancreas, bone, kidney, liver, adrenal, and pituitary glands, with toxic effects and a half-life of nine years in bones (O’Neal and Zheng, 2015). Mn toxicity is associated with dopaminergic dysfunction, and Parkinson’s and Alzheimer’s diseases (Bjørklund et al., 2019).
Excessive Mn in animals and plants predominantly results from a high level of Mn in the soil environment. There are seven valence states of Mn, and the dominant states in natural environments are Mn(II), Mn(III), and Mn(IV) oxides. Mn(II) is the main form toxic to life, and because it is soluble in water so it can accumulate in human cells via food chains. Mn in other valence states in natural environments can be transformed into Mn(II), which can easily permeate into ground water and contaminate the aquatic environment. Many countries, including Australia, Bangladesh, China, and the United States, have reported that the Mn concentration in ground water exceeds the permitted limit (Ali et al., 2017; McMahon et al., 2019; Adeyeye et al., 2020; Ghosh et al., 2020).
Mining and smelting are the two main source of contamination of ground water, rivers, streams, and lakes. In China, some of the largest Mn industries are built in Songtao county, northeast of Guizhou. The average amount of Mn(II) in the surface water of the Songtao Manganese Mine is 39.67-times higher than the safety standard for drinking water in China (0.1 mg/L), and almost half of the water sources in this area fail to meet the drinking water standard. Moreover, Mn pollution has impacted the growth and development of children living nearby. For remediation, chemical methods such as filtering with manganese sand are often used to reduce the high concentrations of Mn in wastewater from mine factories. However, not all Mn(II) is oxidized into Mn(III) or Mn(IV) by oxygen in chemical oxidation processes, and it is impossible for Mn metal to be completely eliminated via abiotic sedimentation and aeration (Li et al., 2019). Abiotic Mn(II) oxidation is thermodynamically favored at a higher pH, but many bacteria are capable of transforming soluble Mn(II) to its oxide precipitates of Mn(IV) oxides in neutral or acidic pH conditions (Bohu et al., 2015; Piazza et al., 2019). Manganese-oxidizing bacteria (MOB) harbor multicopper oxidase (MCO)-type enzymes, manganese oxidase, or heme peroxidase to catalyze Mn(II) oxidization (Anderson et al., 2009; Su et al., 2013; Das et al., 2020; Tang et al., 2021). The Mn(II) removal efficiency by enzymes of MOB is approximately five-times higher than that of oxidation reaction by O2 in natural environments (Katsoyiannis and Zouboulis, 2004; Su et al., 2013). Furthermore, some MOBs can survive in media with a Mn(II) concentration as high as 3,000 mg/L (Wan et al., 2017). In contrast, monkeys displayed abnormal cognition, motor functioning, and behavior when the Mn(II) concentration was only 29.4–73.7 μg/L in blood (Schneider et al., 2006). However, it remains to be elucidated how MOBs adapt to environments featuring a high concentration of toxic Mn. In the present study, MOBs were isolated from the soil of a manganese mine, and the response mechanism of these bacterial strains to manganese stress were explored.
Materials and Methods
Isolation and Screening of Manganese-Tolerant Bacterial Strains
Samples were collected from 0- to 20-cm depths of the subsoil located in the Songtao manganese mine (108°51′50″ E, 27°58′30″ N), Guizhou province, China. Soil in this area is heavily polluted with manganese as a result of anthropogenic mining activities; the total Mn in the sampled soils was in the range of 2,000–2,300 mg/kg. Soil samples were suspended in sterile water at a ratio of 1:9, and the mixtures were vigorously vortexed for 10 min. Approximately 80 μL supernatant was plated on PYCM solid medium [peptone 0.8 g/L, yeast extract 0.2 g/L, K2HPO4 0.1 g/L, MgSO4⋅7H2O 0.2 g/L, NaNO3 0.2 g/L, CaCl2 0.1 g/L, (NH4)2CO3 0.1 g/L, 1.5% agar, pH 7.0] (Tian et al., 2017) containing 200 mg/L MnCl2, and plates were incubated at 28°C for 7 days. Single colonies were selected and screened by plating on solid PYCM media containing MnCl2 concentrations from 700 mg/L up to 2,200 mg/L.
Identification of Isolated Strains
Total genomic DNA was extracted from each isolated strain using a TIANamp Bacteria DNA Kit [Tiangen Biotech (Beijing) Co., China] according to the manufacturer’s protocol. The 16S rRNA and gyrA genes were amplified by PCR using the bacterial genomic DNA as templates (Weisburg et al., 1991). Purified PCR products were inserted into the pMD19-T vector (Takara Co., Japan), transformed into Escherichia coli DH5α and sequenced by the Sanger method. Consensus phylogenetic trees were reconstructed by the UPGMA method using the MEGA v7 program (Kumar et al., 2016b) with 1,000 bootstrap repetitions.
Measurement of Manganese Removal Efficiency of Bacteria
Bacterial suspensions in exponential growth phase were inoculated in a 1:20 ratio into PYCM broth containing 250 mg/L MnCl2 and were cultured at 180 rpm and 28°C for 7 days. The culture supermatant from day 1 to day 7 was sampled and the Mn(II) in the samples was quantitated in triplicate using inductively coupled plasma-optical emission spectrometry (ICP-OES) (PQ9000, Analytik Jena, Germany). Mn removal efficiency was calculated via a previously reported equation (Shu et al., 2019).
Growth of Bacterial Strain Under Mn(II) Stress
The isolated bacterial strain was cultured in PYCM liquid medium containing MnCl2 at concentrations of 0, 250, 500, 1,000, 1,500, and 2,200 mg/L. The optical density at 600 nm (OD600) and pH values of the bacterial suspension were recorded at various time points from 0 to 64 h. The biofilm formation of the strain was detected according to a previous description (Auger et al., 2006).
Growth of Bacterial Strain in Media Containing Heavy Metals
The growth and tolerance of the selected bacterial strain to five heavy metals were measured according to the method reported by Oladipo et al. (2018). Briefly, stock solutions of heavy metal salts were prepared, including CdSO4, Al2(SO4), CuSO4, FeCl3, and K2Cr2O7. The bacterial strain was incubated at 180 rpm and 28°C for 72 h in PYCM liquid medium containing the different metal solutions at concentrations of 0, 100, 250, 500, and 1,000 mg/L. All tests were repeated four-times. The OD600 of the bacterial cultures was determined at 4, 8, 12, 16, 20, 24, 36, 48, and 72 h.
Motility of Bacterial Strains Under Mn(II) Stress
A puncture was created in a plate of PYCM solid media, and the isolated bacterial strain was inoculated to determine bacterial motility under Mn stress. The plate contained 0.3% agar to assess swimming motility, while 0.6% agar was used to test swarming motility. An equivalent number of bacteria in exponential growth phase were inoculated into PYCM soft-agar plates supplemented with or without 250 mg/L Mn(II) and were incubated for 8–32 h. The diameter of the spreading colony was measured, and swimming/swarming areas were calculated according to a previous report (Hariri et al., 2016). Motility detection in mutants with flagellar genes knocked out utilized the same method as that of wild strain ST7 without Mn supplied in the media. All assays were performed in quadruplicate.
Mn(II) Oxidation Activity Assays
The isolated strains were spread on solid-agar PYCM plates and cultured at 28°C for 7 days. Mn oxide products, including Mn(III) and Mn(IV), on plates or in liquid cultures were monitored by using the colorimetric dye solution of leucoberbelin blue (LBB; 0.04%, w/v) (Krumbein and Altmann, 1973). Manganese oxidation activity of strains was quantified by employing the LBB method and using KMnO4 solution for the standard curve as previously described (Johnson and Tebo, 2008).
Visualization of Mn(II) Oxides by Scanning Electron Microscopy
Bacteria cultured at 28°C for 7 days in PYCM media (refer to Section “Isolation and screening of manganese-tolerant bacterial strains”) with 2,200 mg/L MnCl2 were collected and washed in phosphate-buffered saline (PBS, pH 7.4). Bacterial pellets were prepared and observed using scanning electron microscopy (SEM) as previously reported (Bohu et al., 2015). Attachments on the bacterial surface were observed via a Hitachi S-3400N scanning electron microscope with 20,000 V acceleration voltage.
Analysis of Transcriptome Profile of Bacterial Strain Under Mn(II) Stress
The isolated bacterial strain was incubated in PYCM liquid medium in the absence or presence of 250 mg/L MnCl2. Three samples were collected from each culture: samples L01–L06 were from mid-exponential growth phase cultures while samples L07–L12 were from the onset of stationary phase cultures. Samples L01–L03, obtained in the absence of Mn(II), were collected after culture for 8 h [colony-forming unit (cfu)/mL = 3.183 × 108], while samples L04–L06 were obtained following culture in the presence of Mn(II) for 16 h (cfu/mL = 3.182 × 108). Samples L07–L09 were collected from media without Mn(II) after culture for 16 h (cfu/mL = 4.635 × 108), while samples L10–L12 were obtained following culture in the presence of Mn(II) for 24 h (cfu/mL = 4.573 × 108). Based on the protocol of Genedenovo Biotechnology Co., Ltd. (Guangzhou, China), 12 cDNA libraries were constructed and sequenced via an Illumina HiSeqTM 2500 platform.
The raw data were filtered using the fastp Toolkit1 to remove low-quality sequences, which were sequences with > 10% unidentified nucleosides (N), sequences with > 50% bases containing a phred quality score < 20, or sequences containing only the barcode adapter. The remaining sequences were then aligned with the reference genome of Bacillus safensis KCTC 12796BP (assembly no. GCF_001895885.1, containing 4,058 genes) using the STAR v2.7 program2 and allowing no mismatches, and the sequences mapping to rRNA were removed. Expression patterns of genes were calculated according to the value of counts per million (CPM) to standardize the gene expression levels. The DESeq2 and edgeR packages in the R v4.0.4 platform3 were used to identify differentially expressed genes (DEGs) across groups. Taking the FDR (false discovery rate) value of 0.05 as the threshold, up- or downregulated genes were recognized if the value of |log2 (fold change)| was ≥ 1. These data were further displayed in a volcano plot created with the SangerBox program4 taking the FDR value and log2FC as coordinate axes. Nucleotide and amino acid sequences of the DEGs were generated from the assembled transcripts by Trinity v2.85 (Haas et al., 2013). Amino acid sequences of the DEGs were inputted for Kyoto Encyclopedia of Genes and Genomes (KEGG) orthology via the KEGG Automatic Annotation Server (KAAS)5. Gene lists from 30 species of the genus Bacillus containing 148,323 complete genome sequences were used as references. The expressed genes were further examined to find those related with lignocellulose degradation based on KAAS hierarchical categories and previously archived literature.
Validation of Differentially Expressed Genes
The same aliquot of total RNA used for RNA-seq was also used to validate the DEGs by RT-qPCR. Gene-specific primers were designed using Primer 5.0 software (Supplementary Table S1) and using the 16S rRNA gene as the internal reference gene. The qPCR reaction was performed according to a previous study (Livak and Schmittgen, 2001). The corresponding RT-qPCR efficiency (E) was qualified in the range of 90.1–101.9% (Supplementary Figure S1). All assays were performed in triplicate.
Gene Knockout by Homologous Recombination Technology
Homologous recombination technology was used to knock out 15 genes by single-cross integration as previously reported (Liu et al., 2019). In brief, the partial fragment F1 without the stop codon of the target gene was amplified from genomic DNA of wild-type B. safensis ST7. The complete kanamycin gene (fragment F2) was amplified from plasmid pPIC9K (Invitrogen Life Technologies) by PCR using primers km-F and km-R (Supplementary Table S1). The combined fragment F3, comprising the incomplete target gene (fragment F1) fused upstream of the kanamycin gene (fragment F2), was generated by PCR (Liu et al., 2019). After digestion with BamHI, the F3 fragment formed a circular DNA and was electroporated into wild-type B. safensis ST7 cells.
The gene-knockout colonies were selected on solid PYCM media plates supplemented with 5 μg/mL kanamycin. The deletion of the target gene was verified by PCR tests and sequencing. The Mn(II) oxidation ability of each mutant (Δgene) was assessed by Mn(II) oxidation activity assays using the LBB and ICP-OES methods, as previously described (refer to Section “Mn(II) oxidation activity assays”), with the wild-type strain as a control. Mutants harboring knocked out flagellar genes were further used to detect motility, as described above (refer to Section “Motility of bacterial strains under Mn(II) stress”). Spore production was assayed in mutants containing deletions of sporulation-related genes according to the protocol of a spore staining kit (Shanghai Solarbio Bioscience & Technology Co., China, Cat no. G1132), using malachite green as the dye.
Results
Isolation and Identification of Manganese-Oxidizing Bacteria Strains
Soil was sampled from an abandoned manganese mine located in Songtao county, Guizhou, China, that had not performed smelting for 7 years, and 50 bacterial strains were isolated on media supplemented with 200 mg/L Mn(II) (Supplementary Table S2). The isolates were further tested in broth media containing different concentrations of Mn(II), and it was determined that seven strains were able to tolerate exposure to 1,800–2,200 mg/L Mn(II). These seven strains were characterized as Rhizobium sp., Arthrobacter oxydans, and B. safensis based on nucleotide sequence similarities of the 16S rRNA gene and gyrase subunit A gene (gyrA) (Supplementary Figures S2, S3).
Determination of Manganese Absorption and Oxidation Ability of the Isolated Strains
The isolated bacterial strains were investigated for their manganese removal efficiencies in PYCM media. After incubation for 7 days, the bacterial isolates presented manganese removal efficiencies ranging from 11 to 82% (Supplementary Table S2), with strain ST7 displaying the highest Mn removal rate. Higher efficiencies of manganese absorption were detected from day 3 until day 7 (Supplementary Figure S4). The Mn(II) oxides generated by isolate ST7 were further measured by LBB solution based on the characteristics of LBB dye being specifically oxidized by Mn(III) and Mn(IV) (Krumbein and Altmann, 1973).
A blue color appeared in an ST7 colony cultured on solid media supplemented with 250 mg/L MnCl2, and the color became much deeper on the plate containing 2,200 mg/L MnCl2 (Supplementary Figure S5). The oxidation capacity of isolate ST7 was high, with a manganese oxide biomass of 7.75 μM/d. Numerous irregular precipitates covering the surface of isolate ST7 were visible by SEM after incubation for 7 days in 2,200 mg/L MnCl2 liquid media (Supplementary Figure S5). Strain ST7 was thus characterized as a MOB with a high capacity for Mn removal. Furthermore, the swimming motility of strain ST7 was stimulated by Mn(II) in a 0.3% soft-agar plate supplemented with 250 mg/L Mn(II) (Supplementary Figure S6).
Growth of Strain ST7 Under Mn(II) Stress
Strain ST7 grew rapidly in PYCM medium without manganese (Supplementary Figure S7), reaching mid-exponential phase in 8 h and stationary phase in 16 h. When the strain was cultured in media supplemented with 2,200 mg/L MnCl2, the time to stationary phase increased to 44 h. The growth rate of strain ST7 in 250 mg/L MnCl2 was reduced to half of that of the control. Thus, 250 mg/L MnCl2 was considered as the half-effective dosage (ED50) for strain ST7. Furthermore, in PYCM media with a small amount of peptone and yeast extracts as an energy source, an increase in the pH value of cultures was observed during the exponential growth phase of strain ST7, which later returned to the initial pH of 7 at the stationary growth phase for all three cultures with or without Mn(II) stress (Supplementary Figure S8). In addition, the biofilm formation of strain ST7 was obviously promoted by the presence of Mn(II). After culturing in LB media for 72 h in a microtiter plate, the OD595 value of strain ST7 under Mn(II) stress was 3.340 ± 0.699, and was 0.715 ± 0.036 in the control wells lacking Mn(II) (Supplementary Figure S9).
Maximum Tolerance of Strain ST7 to Heavy Metals
The growth of strain ST7 was measured in minimal PYCM media containing heavy metals at various concentrations. Strain ST7 grew in media supplemented with Al(III), Cr(VI), or Fe(III), but not in media containing Cd(II) or Cu(II). The maximum metal tolerance of strain ST7 to Al(III), Cr(VI), and Fe(III) was 500 mg/L, 250 mg/L, and 250 mg/L, respectively (Supplementary Figure S10).
Expression Profile of Bacillus safensis ST7 During Mn(II) Exposure
To investigate the response mechanism of strain ST7 to manganese stress, 12 cDNA libraries were constructed from samples cultured in the presence or absence of 250 mg/L Mn(II). Whole-transcriptome gene expression profiles were analyzed using an Illumina HiSeq X-ten platform. All cDNA libraries generated 159.1 × 106 clean sequences after quality control and filtering. There were similar matching percentages for all samples, with an average of 80.77% of sequences mapping onto the reference genome of B. safensis, and the Q30 was greater than 89.34% (Table 1).
Table 1. Statistical statement of transcriptome sequencing of B. safensis strain ST7.
A total of 3,755 genes was annotated from the 12 libraries. The ratios of the expressed genes were 85.74–90.18% of the total reference genes (4,067). The normalized CPM data are listed in Supplementary Table S3. The expression levels for most genes were lower than 10,000 CPM (Figure 1). The expression profiles of strain ST7 differed in the absence and presence of Mn(II) (Figure 1). In total, 3,668 genes were expressed in the absence of Mn(II) at the mid-exponential growth phase. However, in the presence of Mn(II), the expression of 53 genes that functioned in the process of antimicrobial resistance was completely inhibited. Additionally, the expression of 64 genes enriched in the processes of transporting and DNA repair was induced by Mn(II) (Figures 2, 3). During the stationary growth phase, the numbers of expressed genes decreased to 3,487 and 3,547 in the absence and presence of Mn(II), respectively. There were 142 genes expressed in the presence of Mn(II) in mid-exponential growth phase that were not detected in the stationary growth phase. These included genes linked to the systems of membrane transport and secretion signaling processes.
Figure 1. Gene expression patterns of strain ST7 in the presence of Mn(II) at exponential (stage 1) and stationary (stage 2) growth phases.
Figure 2. Venn diagram showing the overlap of numbers of expressed genes in strain ST7 in the presence or absence of 250 mg/L Mn(II). Numbers of expressed genes in (A) exponential growth phase (EGP) and (B) stationary growth phase (SGP) are indicated.
Figure 3. Enriched pathways of genes affected by manganese exposure. (A) Enriched pathways of genes expressed during the exponential growth phase (EGP) of strain ST7 in the presence or absence of Mn(II). (B) Enriched pathways of genes expressed during stationary growth phase (SGP) of strain ST7 in the presence or absence of Mn(II).
Two software packages, edgeR and DESeq2, were used to analyze the DEGs. In total, 2,580 genes were differentially expressed in the presence and absence of Mn(II) (Supplementary Table S4). Based on KAAS analysis, the DEGs were enriched in six hierarchies including metabolisms, genetic information processing, environmental information processing, cellular processes, organism systems, and human disease-related processes (Figure 4). The enrichment categories of highly expressed DEGs with CPM above 1,000 are listed in Supplementary Table S5.
Figure 4. Enrichment hierarchies of differentially expressed genes (DEGs) in strain ST7 under manganese stress. The horizontal ordinate represents the number of enriched genes in each hierarchy. Numbers of enriched DEGs that were upregulated, downregulated, and without difference in exponential growth phase (EGP) are indicated as EGP-up (the upregulated DEGs at EGP), EGP-down (the downregulated DEGs at EGP), and EGP-no diff (genes with no significance at EGP), respectively. Numbers of DEGs that were upregulated and downregulated in stationary growth phase (SGP) are shown as SGP-up (the upregulated DEGs at SGP) and SGP-down (the downregulated DEGs at SGP), respectively.
At the mid-exponential growth phase, 1,205 DEGs were induced by Mn(II) and exhibited | logFC| ≥ 1, and of these genes, 599 were upregulated and 606 were downregulated (Figure 5A). The range of log2FC values varied from −3.82 to 5.99. The upregulated DEGs were enriched in pathways of energy metabolism, nucleotide metabolism, translation, and flagellar assembly (Supplementary Table S5), while the downregulated genes were clustered in quorum sensing, peptidoglycan biosynthesis, and propanoate metabolism. In the stationary growth phase (Supplementary Table S4), many more DEGs (2,133) with greater changes of log2FC values from −13.76 to 8.62 were identified compared with the DEGs from the mid-exponential growth phase. Of the DEGs in the stationary growth phase, 1,078 genes were upregulated and 1,055 genes were downregulated in the presence of Mn(II) (Figure 5B). Genes related to the metabolism of pyruvate and propanoate were upregulated, while those connected with ribosome function were largely inhibited by Mn(II) exposure (Supplementary Table S5).
Figure 5. Volcano plots of DEGs in strain ST7 under manganese stress. (A) DEGs in mid-exponential growth phase. (B) DEGs in the onset of stationary growth phase. Green and blue dots represent the downregulated genes, while red and orange dots represent the upregulated genes.
The top 20 genes with the highest expression levels were analyzed further (Tables 2, 3). Of these, ten genes did not show significant difference in the absence or presence of Mn(II) at the exponential growth phase, and these genes were involved in processes such as carbohydrate metabolism, lipid metabolism, cellular community in quorum sensing, and transcription. However, the expression of four genes related to carbohydrate and lipid metabolisms, and cellular community in quorum sensing was significantly increased by Mn(II) stress, and six genes involved in sporulation and protein metabolism were significantly downregulated during Mn(II) exposure. When the ST7 growth reached its stationary phase, most of the top 20 genes were upregulated, including 23S rDNA and genes enriched in the processes of transcription, amino acid metabolism, carbohydrate metabolism, lipid metabolism, exosomes, and transporters. Only three genes of ncRNAs and transfer-messenger RNA were downregulated during the stationary period.
Table 2. Top 20 highly expressed genes in the mid-exponential growth phase of strain ST7.
Table 3. Top 20 highly expressed genes in the stationary growth phase of strain ST7.
It was interesting that 26 genes from the transcriptome of strain ST7 were predicted to be related to lignocellulosic degradation (Supplementary Table S6). Of these, genes coding for three main types of cellulases were expressed at a high level, especially in media without Mn(II) stress, which included endoglucanases, β-glucosidases, and cellobiohydrolase. The other 15 genes participated in cellobiose transportation, and production of ethanol, acetate, and lactate from cellulose. It was ascertained that approximately six genes were hemicellulases such as xylanase, according to the hierarchies of KAAS enrichment analysis. Furthermore, two genes coding for manganese catalase and multicopper oxidase domain-containing protein (MCO) might participate in lignin degradation, which was the case for manganese peroxidase and bacterial laccase based on a homology blast with Bacillus pumilus (Kumar et al., 2016a; Ogonda et al., 2021).
Validation of RNA-Seq Data by RT-qPCR
Seventeen expressed genes from the RNA-seq data were selected for validation by RT-qPCR. The expression patterns obtained by RT-qPCR for these genes were consistent with the trends of transcriptomic abundance changes based on the RNA-seq data (Figure 6). To demonstrate the function of the DEGs in Mn(II) oxidation, 15 DEGs with high CPM values or high differential fold-changes were knocked out from the wild-type ST7 genome using homologous recombination technology by a single exchange (Figure 7 and Table 4). First, sporulation and Mn(II) oxidation abilities were detected in mutants of five sporulation-related genes. Sporulation and Mn(II) oxidation were completely abolished in two mutants, ΔRS10620 (glycosyltransferase family 2 protein) and ΔRS10635 (glycosyltransferase) (Figure 7 and Supplementary Figure S11), while the Mn(II) oxidation ability decreased by 55.78% in mutant ΔRS03010 (deletion of MCO-domain-containing protein). However, the sporulation and Mn(II) oxidation abilities were not affected by deletions of the RS05775 or RS11215 genes, which encode spore coat protein and stage IV sporulation protein A, respectively.
Figure 6. RT-qPCR validation of expressed genes in strain ST7 cultured in media supplemented with 250 mg/L Mn(II).
Figure 7. Manganese oxidation ability of mutants. Mutants with deleted genes and the wild-type strain were cultivated in PYCM liquid media supplemented with 250 mg/L Mn(II) for 7 days. (A) Manganese removal rates of mutants detected by ICP-OES method. (B) Manganese oxidation abilities of mutants detected by adding 150 μL of 0.04% LBB solution into the culture supernatants.
Table 4. Information of fifteen genes knocked out via homologous recombination technology.
Next, manganese oxidation activity and motility were detected in the mutants of six flagellar genes (Supplementary Figure S12). The swimming abilities of mutants ΔRS07925, ΔRS07865, and ΔRS07880 decreased by 87.13, 81.83, and 68.78% respectively, demonstrating that these flagellar genes were crucial for the motility of isolate ST7. Four mutants (ΔRS07865, ΔRS07895, ΔRS07920, and ΔRS07980) lost their Mn(II) oxidation abilities. These genes encode the flagellar basal body M-ring protein, flagellar hook-length control protein, flagellar motor switch protein, and the chemotaxis protein, respectively. No obvious changes were detected in the Mn(II) oxidation abilities of the other two mutants, ΔRS07880 and ΔRS07925.
Finally, Mn(II) oxidation abilities were detected in mutants of four transporter genes. Mn(II) oxidation was completely abolished in mutant ΔRS02135, while mutant ΔRS17055 exhibited decreased efficiency in Mn(II) oxidation from 82% to 17.11%. Disruption of the other two genes, RS06145 and RS16735, encoding phosphotransferase system (PTS) galactitol transporter subunit IIC and iron ATP-binding cassette (ABC) transporter permease, respectively, did not have any effect on Mn(II) oxidation.
Discussion
Isolation of Manganese-Oxidizing Bacteria From Soil
Toxic Mn(II) can be removed from different aqueous environments by MOB such as Bacillus sp. strain SG-1 (de Vrind et al., 1986), Leptothrix discophora strain SS-1 (Corstjens et al., 1997), Flavobacterium spp. (Hou et al., 2020), Exiguobacterium sp. CNU020 (Lee et al., 2009), and Pseudomonas putida strains MnB1 and GB-1 (Okazaki et al., 1997). In addition, several MOB—namely, Pseudomonas sp. nov., Rhizobium spp. (Moy et al., 2003), Bacillus cereus CP133 (Wang et al., 2016), Bacillus megaterium 1Y31 (Zhang et al., 2015), Acinetobacter spp. (Ren et al., 2019), and Citrobacter freundi (Tang et al., 2014)—were isolated from alfisol soil or mine soil. In the current study, 50 manganese-tolerant bacterial strains were isolated from the soil of the Songtao manganese mine (Supplementary Table S2), the largest mine in Asia. Most of the isolated strains were bacteria that are commonly distributed in soil, and the dominant genus was Bacillus. Some of the isolated bacterial taxa had previously been recognized for their abilities and characteristics related to manganese oxidation; the exceptions were Arthrobacter oxydans and B. safensis. The isolated bacteria showed diversity in Mn(II) removal efficiency even if the isolates were classified in the same genus.
Of the strains isolated in the current study, B. safensis strain ST7 was the most effective MOB. The maximum Mn(II) concentration tolerated by strain ST7 was 2,200 mg/L, and thus, up to 82% removal of Mn(II) was achieved, with production of 7.75 μM/d biomass of manganese oxides [Mn(III)/Mn(IV)]. The Mn(II) tolerance ability of strain ST7 was close to that of Arthrobacter sp. strain HW-16 (oxidation ratio of 66.28% in 3,000 mg/L Mn(II) media) (Wan et al., 2017), and much higher than that of B. cereus strain P1 (Fan et al., 2016) and Streptomyces spinoverrucosus strain NB-7 (Wang et al., 2013).
The current study demonstrated that strain ST7 can also tolerate multi-metal stress. The manganese mine region in Songtao county is polluted with numerous heavy metals, such as Cd(II), Al(III), Cr(VI), Cu(II), Fe(III). Therefore, the growth of B. safensis strain ST7 was explored in media containing these heavy metals. Isolate ST7 grew in the presence of three of these heavy metals, tolerating Al(III) at 500 mg/L, and Cr(VI) and Fe(III) at 250 mg/L. The maximum metal tolerance of strain ST7 was much higher than the threshold records in soil (Kabata-Pendias, 2010), as it existed in sandy soil with 1,000–2,223 μg/L iron at pH 2.5–4.5, and in soil with 5,700 μg/L Al at pH 4.4. The background content of Cr in soil ranged from 14 to 109 mg/kg, and the trigger action value (TAV) for Cr(VI) is 3–25 mg/kg. Moreover, the tolerance capacity of strain ST7 was higher than that of B. cereus and Bacillus kochii for Al(III) and Fe(III) from mining sites (Oladipo et al., 2018), and higher than that of Bacillus sp. S3 for Cr(VI) from soil (Zeng et al., 2020). Hexavalent chromium [Cr(VI)] is a strong oxidant, teratogen, and carcinogen that easily crosses biomembranes and is toxic to the proteins, genomic DNA, and lipids of bacteria (Pushkar et al., 2021). The solubility of aluminum is increased in acidic soil, and Al(III) toxicity severely reduced plant growth and productivity globally (Gu et al., 2021). Strain ST7 might therefore be considered as a valuable resource for bioremediation of multiple heavy metals.
Gene Expression Profiles of Strain ST7 Under Mn(II) Stress
In a previous study (Hoover et al., 2010), it was proposed that the different capabilities for Mn removal might be related to diverse genetic mechanisms and patterns of gene expression, especially for those proteins related to manganese oxidation such as MCO, which oxidizes Mn(II) into Mn oxides in Pseudomonas (Brouwers et al., 1999) and Bacillus pumilus (Su et al., 2013). However, the genes required to eliminate the toxicity of Mn(II) vary in the known MOBs. Moreover, for specific microorganisms that inhabit the soil in extreme environments, there is limited information available regarding their mechanisms of adaptation to high concentrations of Mn(II).
To investigate the adaptation mechanism of bacteria under Mn exposure, B. safensis strain ST7 was used for RNA-seq analysis of gene expression profiles in toxic Mn(II) media. A total of 3,755 genes from 12 cDNA libraries was detected, which included most genes in the reference transcriptome of B. safensis. There were 3,732 genes expressed in the exponential phase of growth, which decreased to 3,577 genes in the stationary phase. Mn(II) exposure specifically induced 64 genes and completely inhibited 53 genes in strain ST7. The genes induced by Mn(II) exposure in the exponential growth phase of strain ST7 were involved in the processes of transport and DNA repair, while the genes inhibited by Mn(II) exposure in this growth phase were involved in the processes of DNA replication and cell division. The patterns of specifically expressed genes might account for the slow growth of strain ST7 when encountering such a high concentration of Mn(II) (Supplementary Figure S7). Furthermore, in the stationary phase of growth, expression of genes related to sporulation was inhibited, while genes involved in membrane transport and secretion signaling were induced by Mn(II) exposure (Figure 3). This was similar to the profile of Escherichia coli MB266 obtained using iTRAQ-based proteomics analysis, whereby the proteins for ion and protein export and secretion systems were activated by Mn(II) (Wang et al., 2017).
It was recognized that 26 genes from the transcriptome of strain ST7were associated with lignocellulosic degradation (Supplementary Table S6). Lignocellulose in plant cell walls is composed of three components, namely cellulose, hemicellulose and lignin in the proportions of 35–50%, 25–35%, and 10–25%, respectively (Mazzoli, 2020). Cellulose is degraded into small molecules of chemicals mainly by cellobiohydrolases, endoglucanase, and β-glucosidases. Endoxylanase and β-xylosidase are two main enzymes involved in hemicellulose hydrolysis. The biodegradation of lignin requires oxidase and peroxidase, such as manganese peroxidase and laccase (Ma et al., 2020). Although no Bacillus members can ferment lignocellulosic material without prior pre-treatment of biomass saccharification, genes from ST7 transcriptome coding for cellulases, hemicellulases, and enzymes to decompose lignin were expressed at a high level according to the hierarchies of the KAAS enrichment analysis, and homology blasts with other Bacillus species (Kumar et al., 2016a; Ogonda et al., 2021). The other genes were related to the production of high-value chemicals such as ethanol, acetate, and lactate from cellulose. This indicated that strain ST7 has the potential to degrade lignocellulose while its capacity was inhibited by high concentrations of Mn(II).
Differentially Expressed Genes Under Manganese Stress
Previous studies in yeast (Andreeva et al., 2017) and Bacillus subtilis (Mhatre et al., 2016) suggested that some DEGs based on RNA-Seq data may be associated with adaptation to excessive Mn(II). In the current study, 2,580 DEGs were identified in strain ST7. Many DEGs were highly expressed in ST7 upon Mn(II) exposure. These highly expressed genes might be pivotal to the basic processes for bacterial survival, such as carbohydrate and lipid metabolism, cellular community in quorum sensing, and transcription and translation, which would be necessary in the response to Mn stress. For instance, gene RS09610 (encoding aconitate hydratase AcnA) was upregulated by Mn(II) with an increase greater than three-fold. AcnA modulates the efficiency of sporulation in B. subtilis (Craig et al., 1997), and a protein from marine spores can oxidize manganese (van Waasbergen et al., 1996). Thus, AcnA of ST7 might regulate the response to the toxic Mn(II) environment.
Gene RS03010 of strain ST7 encoded an MCO-domain-containing protein belonging to the MCO family. The enzyme MCO catalyzes the transformation of Mn(II) into Mn(III/IV) (Bohu et al., 2015). Mn oxidation in bacteria is an important requirement for MOB to survive in extremely toxic environments. Mn oxidation in MOB indirectly expends ROS in the creation of Mn(II) oxide and relieves ROS stress on bacterial cells. Furthermore, biological Mn(III/IV) oxides can catalyze the degradation of refractory compounds into low-molecular-weight substrates as nutrients for bacteria and protect bacteria against toxic heavy metals such as arsenic (Sunda and Kieber, 1994; Tebo et al., 2005). The protein encoded by RS03010 in strain ST7 showed 99% amino acid similarity (504/510) to the reference genome of B. safensis (Supplementary Figure S13).
Amino acid similarities of strain ST7 were much more diverse when compared with proteins from bacteria outside of the genus Bacillus that had been experimentally shown to harbor manganese oxidase (Corstjens et al., 1997; Brouwers et al., 1999; Hullo et al., 2001; Ridge et al., 2007; Dick et al., 2008; Geszvain et al., 2013; Su et al., 2013; Butterfield and Tebo, 2017). However, the four conserved domains for copper binding were identified in the MCO protein of strain ST7 (Supplementary Figures S13, S14). MCO, generally regarded as the indispensable enzyme in the process of Mn oxidation, resides on the outer surface of bacteria and accommodates four copper ions as cofactors (Su et al., 2013; Tang et al., 2021). The copper ions, which bind to conserved domains of MCO, are involved in electron transfer during the oxidation of Fe(II) or Mn(II) (van Waasbergen et al., 1996). Furthermore, in the experimental group with Mn(II) supplied at the exponential phase of growth, the expression level of the gene RS03010 was 1.72-times higher than that in the absence of Mn(II) (Supplementary Table S4). When RS03010 was knocked out, sporulation ceased, and the Mn oxidation capacity decreased to 26.22% in the mutant strain (Figure 7).
CotA, a spore-associated protein in the coat of prespores, is widely regarded as an MCO in B. subtilis (Hullo et al., 2001) and is stimulated under conditions of elevated temperature or pH (Baril et al., 2012). It was shown that proteins encoded by MCO genes from Bacillus sp. GZB, B. pumilus WH4, and Brevibacillus panacihumi MK-8 oxidize Mn(II) into Mn(IV) oxides using a heterologous expression system or site-directed mutagenesis method (Su et al., 2013; Das et al., 2020; Tang et al., 2021). These observations, together with the data from the current study, indicated that gene RS03010 of strain ST7 is important for spore production and might encode the MCO protein needed to oxidize soluble manganese compounds. Bacterial MCOs are recognized as significant contributors to the Mn(II) oxidation process. Bacteria may have developed many strategies to adapt to environments containing extremely toxic quantities of manganese (Bai et al., 2021), and Mn(II) oxidation is a crucial mechanism used by MOB to facilitate survival in high Mn surroundings. The presence of Mn oxides outside of the bacteria provide protection against ultraviolet (UV) ionizing radiation, viral attack, and heavy metal toxicity (Daly et al., 2004; Ghosal et al., 2005).
Clusters of transporter genes with CPM > 1,000 were also upregulated by Mn(II) in the current study, based on KAAS enrichment analysis (Supplementary Table S5). Manganese transport systems encoded by mntABCD and mntH operons have been found from various bacteria such as B. subtilis, Escherichia coli, Staphylococcus aureus, Salmonella enterica, and Corynebacterium diphtheria (Kliegman et al., 2006). The mntABCD and mntH operons are controlled by manganese transport regulator MntR. However, except for MntR, these manganese transport genes were not annotated in the ST7 transcripts, and this may due to the imperfect annotation of the reference genome of B. safensis (assembly no. GCA_001895885.1). Clean data from the libraries were therefore realigned using B. subtilis strain 168 as a reference (GCA_000009045.1). Four genes, mntB, mntD, mntH, and mntP, were annotated in strain ST7, and all of them were markedly decreased under Mn(II) stress. Thus, the mntABCD and mntH systems might not function when B. safensis strain ST7 encounters high concentrations of toxic Mn(II). This suggests that different transporter systems in B. safensis might help extrude excessive Mn(II).
Nine genes in strain ST7 were annotated in transporter systems. Two of them were classified as ABC transporters in the cell membrane, including ABC transporter ATP-binding protein and iron ABC transporter permease (encoded by RS16735). Four of the nine genes were substrate-binding proteins (SBPs) that assist ABC transporters in the uptake of metal ions; these were ABC transporter substrate-binding protein, iron (3+)-hydroxamate-binding protein FhuD, iron-uptake system-binding protein, and heme ABC transporter substrate-binding protein IsdE. Amino acid permease (encoded by RS17055) predominantly uptakes amino acids into the bacteria cell, and some amino acid permeases are ABC-type transporters (Quintero et al., 2001).
ATP-binding cassette transporters comprise two domains—the intracellular nucleotide-binding domain (NBD) and the transmembrane domain (TMD) (Lorca et al., 2007). ATP binds with the NBD and provides energy, while the TMD provides a path for cargo to pass through the cell membrane. Additionally, ABC transporters require a SBP to carry the substrate outside the membrane and reach the TMD (Qasem-Abdullah et al., 2017). The mechanism of PTS transporters is markedly different from ABC-type and natural resistance-associated macrophage protein (NRAMP) transporters. The phosphoenolpyruvate-dependent carbohydrate transport system (a PTS encoded by gene RS06145I in strain ST7) performs the translocation together with concomitant phosphorylation of sugars, galactitol, and hexitols (Volpon et al., 2006).
The copper-binding protein encoded by gene RS02135 in strain ST7 is the carrier that binds with Cu(I). Copper-binding protein, CopL, on the cell surface holds four Cu(I) ions on the outer surface of the cell and contributes to the Cu efflux effect of protein CopA or CopB (Rosario-Cruz et al., 2019). One copper-binding protein, prion protein (PrP), binds with both copper and manganese (Samorodnitsky and Nicholson, 2018). In the current study, the mutants with transporter genes RS02135 and RS17055 knocked out lost their Mn oxidation ability (Figure 7). Thus, these transporters might be manganese efflux pumps that remediate high concentrations of Mn(II), similar to the gene mntE in Streptococcus mutans (O’Brien et al., 2020) and gsp gene in Pseudomonas putida (de Vrind et al., 2003).
Several processes were significantly downregulated in strain ST7 under Mn exposure, including sporulation and flagellogenesis. In the sporulation pathway, the downregulated genes encoded enzymes and structural proteins required for spore assembly (Tables 2, 3). Species of the genus Bacillus produce dormant and metabolically inert spores to resist environmental stresses such as extreme pH and temperature conditions (Nicholson et al., 2000). The reported effects of Mn on bacterial sporulation are very diverse. Mn significantly enhanced sporulation in B. cereus (Ryu et al., 2005), but sporulation of B. subtilis strain B68 was completely suppressed under a Mn concentration of 7.2 mM/L (Song et al., 2007). This demonstrates that the regulation of sporulation is tightly controlled and diverse in bacteria.
Genes RS08160 and RS08165 in strain ST7 encode dipicolinic acid synthetase subunits A and B, which participate in the synthesis of dipicolinic acid (DPA) during the sporulation period in the mother cell, and DPA then accumulates in the dormant spores of Bacillus (Paidhungat et al., 2000). Genes RS03530, RS10620, and RS10635 in strain ST7 encode glycosyltransferase, which catalyzes the transference and modification of monosaccharides in the crust during the late stage of sporulation (Shuster et al., 2019). Gene RS15875 encodes KapD, which inhibits spore formation in Bacillus thuringiensis (Fan et al., 2017). The decrease in expression of these genes under Mn exposure could potentially block spore growth of strain ST7 in toxic Mn(II) media. Mutants ΔRS10620 and ΔRS10635 were completely defective in both sporulation and Mn(II) oxidation (Figure 7 and Supplementary Figure S11).
Although the relationship between glycosyltransferase and Mn has yet to be elucidated, it has been reported that spores from several species of the genus Bacillus contain proteins that oxidize Mn(II) into Mn(IV), including MnxG from the marine spores of Bacillus sp. SG-1 (Dick et al., 2008) and Bacillus sp. PL-12 (Butterfield and Tebo, 2017). In B. subtilis, a hallmark of the initiation of sporulation is the activation of Spo0A by phosphorylation through a phosphor-relay system (Rodriguez Ayala et al., 2020). When this system is activated by a starvation signal, a family of autophosphorylate kinases (KinA–KinE) initiates a cascade of phosphate transfers from sigma factor Spo0F to Spo0A. The production and accumulation of phosphorylated Spo0A (Spo0A-Pi) is the master regulator of sporulation. In addition, AbrB is a global gene regulator involved in the transition from exponential to stationary growth phase and antagonizes the sporulation process.
In the current study, AbrB was upregulated by Mn(II) in the stationary growth phase of strain ST7 (Supplementary Table S7), which inhibited the expression of SigH (Rodriguez Ayala et al., 2020). The SigH factor thus might not play a role in promoting the expression of both Spo0F and Spo0A. Furthermore, the high concentration of Mn(II) stimulated expression of kipI (logFC value 5.02), which then repressed the transcript of kinA. KinA is required for the phosphorylation and activation of Spo0F (Supplementary Table S7). In addition, the gene rapA was upregulated and is responsible for dephosphorylation of Spo0F-Pi (Rodriguez Ayala et al., 2020). These genes ultimately act to decrease the accumulation of Spo0A-Pi so that fewer spores are produced and sporulation is much slower in the harsh Mn(II) environment compared with the control environment lacking Mn(II).
Another stress-related pathway, the SigB pathway, was upregulated at the start of the stationary growth phase in strain ST7. It was reported that spore formation is the final survival strategy because it is a process that consumes large quantities of energy, and spore germination and outgrowth need to be tightly regulated to recover the planktonic growth in B. subtilis (Rodriguez Ayala et al., 2020). Furthermore, the SigB pathway is a link that interconnects the two dominant and mutually exclusive adaptive responses—sporulation and general stress response (GSR)—in the surviving regulatory network of B. subtilis and is needed for the bacteria to cope with very harsh conditions (Rodriguez Ayala et al., 2020).
Similar to B. subtilis, the SigB gene was upregulated 2.85-times in strain ST7 cultured in Mn(II)-supplemented media (Supplementary Table S7). The SigB pathway participates in regulation of the stress response under extreme conditions such as heat, salt, pH, and manganese in a B. subtilis model (Rodriguez Ayala et al., 2020). Congruent with this, a high concentration of Mn(II) in the culture media induced expression of the kinase RsbT in strain ST7, which may phosphorylate RsbR analogs and RsbS (Supplementary Table S7). The released RsbT protein can bind another phosphatase, RsbU, to dephosphorylate RsbV-Pi (Rodriguez Ayala et al., 2020). Dephosphorylated RsbV then binds to RsbW, releasing SigB, which in turn combines with RNA polymerase and activates its target genes. For example, the mRNA level of ctc (the general stress protein, similar to ribosomal protein L25) was upregulated approximately 7.39-fold in strain ST7 exposed to Mn(II) (Supplementary Table S7).
In the flagellum pathway, 32 genes were upregulated in strain ST7 at the exponential phase of growth following Mn exposure, but only 10 flagellum genes were expressed at the stationary phase of growth, and nine of these genes were downregulated (Table 5). The flagellar structure of Bacillus is complex and contains three basic structures—the hook, the basal body, and the long filament (Terahara et al., 2020). Of the 32 genes upregulated by Mn exposure in strain ST7 at the exponential phase of growth, genes encoding structural proteins were expressed at a high level. The proteins encoded by these genes were necessary for the assembly of flagella, including FlaD and flagellar capping protein located in the hook; MotA and MotB in the stator; FliD, FlgK, and FlgL in the filament; FlgB, FlgC, FliE, and FliF in the basal body; FliG in the rotor; and FliM and FliY, which are used to control the movement direction of flagellum.
Table 5. Flagellar genes expressed in exponential and stationary growth phases of strain ST7 exposed to manganese.
The functions of flagella change with the growth states of bacteria. B. subtilis exists in two states, including cell chains joined end-to-end and a single bacterium moving alone during the exponential growth phase (Chai et al., 2010). Flagellar genes of B. subtilis are expressed in the exponential growth phase, while the sessile biofilm-forming-related genes are induced in the stationary phase, during which the transition from motile cells to sessile growth is controlled by phosphorylation of modulator degU (Kobayashi, 2007). In strain ST7, swimming motility was stimulated but swarming motility was not affected by Mn(II) (Supplementary Figure S6). The gene RS17705 in strain ST7, corresponding to degU, was upregulated in the exponential phase of growth but was not detected in the stationary phase. Additionally, the biofilm formation of ST7 was promoted by Mn(II). This suggested that motility might be controlled by the degU gene of strain ST7. Furthermore, it has been reported that Mn(II) exposure increases the expression of flagellar genes (e.g., fliI, encoding the sigma factor for the flagellar operon) and the motility of gut bacteria in mice (Chi et al., 2017). Interestingly, disruption of Mn oxidation and a decrease in motility was observed in mutants ΔRS07865 (encoding FliF), ΔRS07895 (encoding flagellar hook-length control protein), ΔRS07920 (encoding FliM), and ΔRS07980 (encoding chemotaxis protein CheA) in the current study (Figure 7). This confirmed that expression of flagellar genes directly affected the capacity for Mn oxidation in strain ST7.
Conclusion
The capacities for Mn removal by the bacterial strains isolated from soil originating from a Mn mine area were very diverse. B. safensis strain ST7 was the most effective MOB among the isolated strains. The adaptation mechanism in B. safensis strain ST7 when confronted by an environment with a high concentration of Mn(II) is multifaceted, involving extracellular precipitation, enzymatic oxidation, and active efflux pumps. Increased expression of an MCO gene to oxidize excessive Mn(II), stimulation of transporter genes to pump excess Mn outside of the cell, and inhibition of expression of sporogenesis genes and flagellar genes to save energy and direct flux were detected. Many metabolic processes were also significantly affected by exposure to Mn. It was determined that an alternative stress-related pathway, the SigB pathway, was upregulated by Mn(II) at the stationary growth phase of strain ST7. The multi-metal tolerance traits of the isolated strain coupled with the high removal rate of Mn suggest that strain ST7 could be an effective bioremediation candidate for the restoration of ecological environments damaged by Mn and other heavy metal pollutants. The use of such a bacterial strain for bioremediation minimizes the use of chemical oxidants and is a potentially cost-effective technology. Additionally, the strain ST7 might degrade lignocellulosic biomass with other microbial collaborators to produce high-value chemicals with high efficiency (such as ethanol, acetate, and lactate) through the microbial fermentation pathway as an alternative to petroleum-based routes.
Data Availability Statement
The sequencing data for this study can be found in the SRA database of National Center for Biotechnology Information (NCBI) (Bioproject number PRJNA734365). The nucleotide sequence of 16S rDNA and gyrase subunit A (gryA) gene of strain ST7 are available in the NCBI under accession number of MT378374 and MT449449.
Author Contributions
XR, SL, and JW designed the work and the experiments. HL, QT, ZZ, LY, XW, QL, SH, and XN performed the experiments and data analyses. XR analyzed the data and wrote the manuscript. JW revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This project was funded by the Guizhou Province "Hundred" Innovative Talents Project (2016–4012), the National Natural Science Foundation of China (31672390 and 31960641), and the Guizhou Province Science and Technology Innovation Team Building Special (QKHPTRC[2019]5615). The funders monitored the ethics of the experiments and provided instruments and financial funds for all the research featured in this manuscript.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors are grateful to the owners for the permission to take soil samples from the Songtao manganese mine, Guizhou province, China.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.758889/full#supplementary-material
Footnotes
- ^ https://github.com/OpenGene/fastp
- ^ https://github.com/alexdobin/STAR
- ^ http://www.r-project.org/
- ^ http://sangerbox.com/Tool
- ^ https://www.genome.jp/tools/kaas/
References
Adeyeye, O., Xiao, C., Zhang, Z., and Liang, X. (2020). State, source and triggering mechanism of iron and manganese pollution in groundwater of Changchun, Northeastern China. Environ. Monit. Assess 192:619. doi: 10.1007/s10661-020-08571-0
Ali, A., Strezov, V., Davies, P., and Wright, I. (2017). Environmental impact of coal mining and coal seam gas production on surface water quality in the Sydney basin, Australia. Environ. Monit. Assess. 189:408. doi: 10.1007/s10661-017-6110-4
Anderson, C., Johnson, H., Caputo, N., Davis, R., Torpey, J., and Tebo, B. (2009). Mn(II) oxidation is catalyzed by heme peroxidases in "Aurantimonas manganoxydans" strain SI85-9A1 and Erythrobacter sp. strain SD-21. Appl. Environ. Microbiol. 75, 4130–4138. doi: 10.1128/AEM.02890-08
Andreeva, N., Kulakovskaya, E., Zvonarev, A., Penin, A., Eliseeva, I., Teterina, A., et al. (2017). Transcriptome profile of yeast reveals the essential role of PMA2 and uncharacterized gene YBR056W-A (MNC1) in adaptation to toxic manganese concentration. Metallomics 9, 175–182. doi: 10.1039/c6mt00210b
Auger, S., Krin, E., Aymerich, S., and Gohar, M. (2006). Autoinducer 2 affects biofilm formation by Bacillus cereus. Appl. Environ. Microbiol. 72, 937–941. doi: 10.1128/AEM.72.1.937-941.2006
Bai, Y., Su, J., Wen, Q., Huang, T., Chang, Q., and Ali, A. (2021). Characterization and mechanism of Mn(II)-based mixotrophic denitrifying bacterium (Cupriavidus sp. HY129) in remediation of nitrate (NO3–N) and manganese (Mn(II)) contaminated groundwater. J. Hazard Mater. 408:124414. doi: 10.1016/j.jhazmat.2020.124414
Baril, E., Coroller, L., Couvert, O., El Jabri, M., Leguerinel, I., Postollec, F., et al. (2012). Sporulation boundaries and spore formation kinetics of Bacillus spp. as a function of temperature, pH and a(w). Food Microbiol. 32, 79–86. doi: 10.1016/j.fm.2012.04.011
Bjørklund, G., Hofer, T., Nurchi, V., and Aaseth, J. (2019). Iron and other metals in the pathogenesis of Parkinson’s disease: toxic effects and possible detoxification. J. Inorg. Biochem. 199:110717. doi: 10.1016/j.jinorgbio.2019.110717
Bohu, T., Santelli, C., Akob, D., Neu, T., Ciobota, V., Rösch, P., et al. (2015). Characterization of pH dependent Mn(II) oxidation strategies and formation of a bixbyite-like phase by Mesorhizobium australicum T-G1. Front. Microbiol. 6:734. doi: 10.3389/fmicb.2015.00734
Bosma, E. F., Rau, M. H., van Gijtenbeek, L. A., and Siedler, S. (2021). Regulation and distinct physiological roles of manganese in bacteria. FEMS Microbiol. Rev. [Epub ahead of print]. doi: 10.1093/femsre/fuab028
Brouwers, G., de Vrind, J., Corstjens, P., Cornelis, P., Baysse, C., and de Vrind-de Jong, E. (1999). cumA, a gene encoding a multicopper oxidase, is involved in Mn2+ oxidation in Pseudomonas putida GB-1. Appl. Environ. Microbiol. 65, 1762–1768. doi: 10.1128/AEM.65.4.1762-1768.1999
Butterfield, C. N., and Tebo, B. M. (2017). Substrate specificity and copper loading of the manganese-oxidizing multicopper oxidase Mnx from Bacillus sp. PL-12. Metallomics 9, 183–191. doi: 10.1039/c6mt00239k
Chai, Y., Kolter, R., and Losick, R. (2010). Reversal of an epigenetic switch governing cell chaining in Bacillus subtilis by protein instability. Mol. Microbiol. 78, 218–229. doi: 10.1111/j.1365-2958.2010.07335.x
Chi, L., Gao, B., Bian, X., Tu, P., Ru, H., and Lu, K. (2017). Manganese-induced sex-specific gut microbiome perturbations in C57BL/6 mice. Toxicol. Appl. Pharmacol. 331, 142–153. doi: 10.1016/j.taap.2017.06.008
Corstjens, P., Vrind, J., Goosen, T., Vrind, J., de Vrind, J., Goosen, T., et al. (1997). Identification and molecular analysis of the Leptothrix discophora SS-1 mofA gene, a gene putatively encoding a manganese-oxidizing protein with copper domains. Geomicrobiol. J. 14, 91–108. doi: 10.1080/01490459709378037
Craig, J., Ford, M., Blaydon, D., and Sonenshein, A. (1997). A null mutation in the Bacillus subtilis aconitase gene causes a block in Spo0A-phosphate-dependent gene expression. J. Bacteriol. 179, 7351–7359. doi: 10.1128/jb.179.23.7351-7359.1997
Daly, M. J., Gaidamakova, E. K., Matrosova, V. Y., Vasilenko, A., Zhai, M., Venkateswaran, A., et al. (2004). Accumulation of Mn(II) in Deinococcus radiodurans facilitates gamma-radiation resistance. Science 306, 1025–1028. doi: 10.1126/science.1103185
Das, R., Liang, Z., Li, G., and An, T. (2020). A non-blue laccase of Bacillus sp. GZB displays manganese-oxidase activity: a study of laccase characterization, Mn(II) oxidation and prediction of Mn(II) oxidation mechanism. Chemosphere 252:126619. doi: 10.1016/j.chemosphere.2020.126619
de Vrind, J., de Groot, A., Brouwers, G. J., Tommassen, J., and de Vrind-de Jong, E. (2003). Identification of a novel Gsp-related pathway required for secretion of the manganese-oxidizing factor of Pseudomonas putida strain GB-1. Mol. Microbiol. 47, 993–1006. doi: 10.1046/j.1365-2958.2003.03339.x
de Vrind, J., de Vrind-de Jong, E., de Voogt, J., Westbroek, P., Boogerd, F., and Rosson, R. (1986). Manganese oxidation by spores and spore coats of a marine Bacillus species. Appl. Environ. Microbiol. 52, 1096–1100. doi: 10.1128/AEM.52.5.1096-1100.1986
Dick, G., Torpey, J., Beveridge, T., and Tebo, B. (2008). Direct identification of a bacterial manganese(II) oxidase, the multicopper oxidase MnxG, from spores of several different marine Bacillus species. Appl. Environ. Microbiol. 74, 1527–1534. doi: 10.1128/AEM.01240-07
Fan, Q., Zhang, S., Gong, Y., and He, J. (2017). Regulation of sporulation by two-component system YvcPQ in Bacillus thuringiensis. Wei Sheng Wu Xue Bao 57, 121–130.
Fan, X., Wang, S., and Li, C. (2016). The isolation and identification of an efficient Fe/Mn oxidizing bacterial strain P1, and the optimization of its oxidizing conditions. Biotechnol. Bull. 32, 175–183.
Gao, X., Sandberg, M., Quach, M., Bodin, B., Johansson, L., and Jansson, L. (2014). Effects of Mn-DPDP and manganese chloride on hemodynamics and glucose tolerance in anesthetized rats. Acta Radiol. 55, 328–334. doi: 10.1177/0284185113499327
Geszvain, K., McCarthy, J., and Tebo, B. (2013). Elimination of manganese(II,III) oxidation in Pseudomonas putida GB-1 by a double knockout of two putative multicopper oxidase genes. Appl. Environ. Microbiol. 79, 357–366. doi: 10.1128/AEM.01850-12
Ghosal, D., Omelchenko, M. V., Gaidamakova, E. K., Matrosova, V. Y., Vasilenko, A., Venkateswaran, A., et al. (2005). How radiation kills cells: survival of Deinococcus radiodurans and Shewanella oneidensis under oxidative stress. FEMS Microbiol. Rev. 29, 361–375. doi: 10.1016/j.femsre.2004.12.007
Ghosh, G. C., Khan, M., Chakraborty, T. K., Zaman, S., Kabir, A., and Tanaka, H. (2020). Human health risk assessment of elevated and variable iron and manganese intake with arsenic-safe groundwater in Jashore, Bangladesh. Sci. Rep. 10:5206. doi: 10.1038/s41598-020-62187-5
Gu, X., Jiang, Y., Wang, X., Jia, H., Li, J., Cui, Y., et al. (2021). Differences in aluminum tolerance and immobilization between two indigenous ectomycorrhizal fungi Lactarius deliciosus and Pisolithus tinctorius from Southwest China’s forest stands. Ecotoxicol. Environ. Saf. 213:112042. doi: 10.1016/j.ecoenv.2021.112042
Haas, B., Papanicolaou, A., Yassour, M., Grabherr, M., Blood, P., Bowden, J., et al. (2013). De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat. Protoc. 8, 1494–1512. doi: 10.1038/nprot.2013.084
Hariri, B. M., Payne, S. J., Chen, B., Mansfield, C., Doghramji, L. J., Adappa, N. D., et al. (2016). In vitro effects of anthocyanidins on sinonasal epithelial nitric oxide production and bacterial physiology. Am. J. Rhinol. Allergy. 30, 261–268. doi: 10.2500/ajra.2016.30.4331
Hoover, S., Xu, W., Xiao, W., and Burkholder, W. (2010). Changes in DnaA-dependent gene expression contribute to the transcriptional and developmental response of Bacillus subtilis to manganese limitation in Luria-Bertani medium. J. Bacteriol. 192, 3915–3924. doi: 10.1128/JB.00210-10
Hou, D., Zhang, P., Wei, D., Zhang, J., Yan, B., Cao, L., et al. (2020). Simultaneous removal of iron and manganese from acid mine drainage by acclimated bacteria. J. Hazard. Mater. 396:122631. doi: 10.1016/j.jhazmat.2020.122631
Hullo, M., Moszer, I., Danchin, A., and Martin-Verstraete, I. (2001). CotA of Bacillus subtilis is a copper-dependent laccase. J. Bacteriol. 183, 5426–5430. doi: 10.1128/jb.183.18.5426-5430.2001
Johnson, H., and Tebo, B. (2008). In vitro studies indicate a quinone is involved in bacterial Mn(II) oxidation. Arch. Microbiol. 189, 59–69. doi: 10.1007/s00203-007-0293-y
Kabata-Pendias, A. (2010). Trace Elements in Soils and Plants, The fourth Edn. Boca Raton, FL: CRC press.
Katsoyiannis, I., and Zouboulis, A. (2004). Biological treatment of Mn(II) and Fe(II) containing groundwater: kinetic considerations and product characterization. Water Res. 38, 1922–1932. doi: 10.1016/j.watres.2004.01.014
Kliegman, J., Griner, S., Helmann, J., Brennan, R., and Glasfeld, A. (2006). Structural basis for the metal-selective activation of the manganese transport regulator of Bacillus subtilis. Biochemistry 45, 3493–3505. doi: 10.1021/bi0524215
Kobayashi, K. (2007). Gradual activation of the response regulator DegU controls serial expression of genes for flagellum formation and biofilm formation in Bacillus subtilis. Mol. Microbiol. 66, 395–409. doi: 10.1111/j.1365-2958.2007.05923.x
Krumbein, W., and Altmann, H. (1973). A new method for the detection and enumeration of manganese oxidizing and reducing microorganisms. Helgoländer Wissenschaftliche Meeresuntersuchungen 25, 347–356.
Kumar, S., Stecher, G., and Tamura, K. (2016b). MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870–1874. doi: 10.1093/molbev/msw054
Kumar, S., Jain, K. K., Rani, S., Bhardwaj, K. N., Goel, M., and Kuhad, R. C. (2016a). In-vitro refolding and characterization of recombinant laccase (cotA) from Bacillus pumilus MK001 and its potential for phenolics degradation. Mol. Biotechnol. 58, 789–800. doi: 10.1007/s12033-016-9978-2
Lee, D., Oh, K., and Kahng, H. (2009). Molecular analysis of antioxidant genes in the extremohalophile marine bacterium Exiguobacterium sp. CNU020. Biotechnol. Lett. 31, 1245–1251. doi: 10.1007/s10529-009-0008-3
Li, N., Huang, T., Mao, X., Zhang, H., Li, K., Wen, G., et al. (2019). Controlling reduced iron and manganese in a drinking water reservoir by hypolimnetic aeration and artificial destratification. Sci. Total Environ. 685, 497–507. doi: 10.1016/j.scitotenv.2019.05.445
Liu, H., Liu, H., Yang, S., Wang, R., and Wang, T. (2019). Improved expression and optimization of trehalose synthase by regulation of Pglv in Bacillus subtilis. Sci. Rep. 9:6585. doi: 10.1038/s41598-019-43172-z
Livak, K., and Schmittgen, T. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods 25, 402–408. doi: 10.1006/meth.2001.1262
Lorca, G., Barabote, R., Zlotopolski, V., Tran, C., Winnen, B., Hvorup, R., et al. (2007). Transport capabilities of eleven gram-positive bacteria: comparative genomic analyses. Biochim. Biophys. Acta 1768, 1342–1366. doi: 10.1016/j.bbamem.2007.02.007
Ma, C., Lo, P. K., Xu, J., Li, M., Jiang, Z., Li, G., et al. (2020). Molecular mechanisms underlying lignocellulose degradation and antibiotic resistance genes removal revealed via metagenomics analysis during different agricultural wastes composting. Bioresour. Technol. 314:123731. doi: 10.1016/j.biortech.2020.123731
Mazzoli, R. (2020). Metabolic engineering strategies for consolidated production of lactic acid from lignocellulosic biomass. Biotechnol. Appl. Biochem. 67, 61–72. doi: 10.1002/bab.1869
McMahon, P. B., Belitz, K., Reddy, J. E., and Johnson, T. D. (2019). Elevated manganese concentrations in United States groundwater, role of land surface-soil-aquifer connections. Environ. Sci. Technol. 53, 29–38. doi: 10.1021/acs.est.8b04055
Mhatre, E., Troszok, A., Gallegos-Monterrosa, R., Lindstädt, S., Hölscher, T., Kuipers, O., et al. (2016). The impact of manganese on biofilm development of Bacillus subtilis. Microbiology 162, 1468–1478. doi: 10.1099/mic.0.000320
Moy, Y., Neilan, B., Foster, L., Madgwick, J., and Rogers, P. (2003). Screening, identification and kinetic characterization of a bacterium for Mn(II) uptake and oxidation. Biotechnol. Lett. 25, 1407–1413. doi: 10.1023/a:1025043326629
Nicholson, W., Munakata, N., Horneck, G., Melosh, H., and Setlow, P. (2000). Resistance of Bacillus endospores to extreme terrestrial and extraterrestrial environments. Microbiol. Mol. Biol. Rev. 64, 548–572. doi: 10.1128/mmbr.64.3.548-572.2000
O’Brien, J., Pastora, A., Stoner, A., and Spatafora, G. (2020). The S. mutans mntE gene encodes a manganese efflux transporter. Mol. Oral Microbiol. 35, 129–140. doi: 10.1111/omi.12286
Ogonda, L. A., Saumonneau, A., Dion, M., Muge, E. K., Wamalwa, B. M., Mulaa, F. J., et al. (2021). Characterization and engineering of two new GH9 and GH48 cellulases from a Bacillus pumilus isolated from Lake Bogoria. Biotechnol. Lett. 43, 691–700. doi: 10.1007/s10529-020-03056-z
Okazaki, M., Sugita, T., Shimizu, M., Ohode, Y., Iwamoto, K., de Vrind-de Jong, E., et al. (1997). Partial purification and characterization of manganese-oxidizing factors of Pseudomonas fluorescens GB-1. Appl. Environ. Microbiol. 63, 4793–4799. doi: 10.1128/AEM.63.12.4793-4799.1997
Oladipo, O. G., Ezeokoli, O. T., Maboeta, M. S., Bezuidenhout, J. J., Tiedt, L. R., Jordaan, A., et al. (2018). Tolerance and growth kinetics of bacteria isolated from gold and gemstone mining sites in response to heavy metal concentrations. J. Environ. Manage. 212, 357–366. doi: 10.1016/j.jenvman.2018.01.038
O’Neal, S., and Zheng, W. (2015). Manganese toxicity upon overexposure: a decade in review. Curr. Environ. Health Rep. 2, 315–328. doi: 10.1007/s40572-015-0056-x
Paidhungat, M., Setlow, B., Driks, A., and Setlow, P. (2000). Characterization of spores of Bacillus subtilis which lack dipicolinic acid. J. Bacteriol. 182, 5505–5512. doi: 10.1128/JB.182.19.5505-5512.2000
Piazza, A., Ciancio, C., Pacini, V., Sanguinetti, G., Ottado, J., and Gottig, N. (2019). Environmental bacteria involved in manganese(II) oxidation and removal from groundwater. Front. Microbiol. 10:119. doi: 10.3389/fmicb.2019.00119
Pushkar, B., Sevak, P., Parab, S., and Nilkanth, N. (2021). Chromium pollution and its bioremediation mechanisms in bacteria: a review. J. Environ. Manage. 287:112279. doi: 10.1016/j.jenvman.2021.112279
Qasem-Abdullah, H., Perach, M., Livnat-Levanon, N., and Lewinson, O. (2017). ATP binding and hydrolysis disrupt the high-affinity interaction between the heme ABC transporter HmuUV and its cognate substrate-binding protein. J. Biol. Chem. 292, 14617–14624. doi: 10.1074/jbc.M117.779975
Quintero, M., Montesinos, M., Herrero, A., and Flores, E. (2001). Identification of genes encoding amino acid permeases by inactivation of selected ORFs from the Synechocystis genomic sequence. Genome Res. 11, 2034–2040. doi: 10.1101/gr.196301
Ren, G., Yan, Y., Nie, Y., Lu, A., Wu, X., Li, Y., et al. (2019). Natural extracellular electron transfer between semiconducting minerals and electroactive bacterial communities occurred on the rock varnish. Front. Microbiol. 10:293. doi: 10.3389/fmicb.2019.00293
Ridge, J., Lin, M., Larsen, E., Fegan, M., McEwan, A., and Sly, L. (2007). A multicopper oxidase is essential for manganese oxidation and laccase-like activity in Pedomicrobium sp. ACM 3067. Environ. Microbiol. 9, 944–953. doi: 10.1111/j.1462-2920.2006.01216.x
Rodriguez Ayala, F., Bartolini, M., and Grau, R. (2020). The stress-responsive alternative sigma factor SigB of Bacillus subtilis and its relatives: an old friend with new functions. Front. Microbiol. 11:1761. doi: 10.3389/fmicb.2020.01761
Rosario-Cruz, Z., Eletsky, A., Daigham, N., Al-Tameemi, H., Swapna, G., Kahn, P., et al. (2019). The copBL operon protects Staphylococcus aureus from copper toxicity: CopL is an extracellular membrane-associated copper-binding protein. J. Biol. Chem. 294, 4027–4044. doi: 10.1074/jbc.RA118.004723
Ryu, J., Kim, H., and Beuchat, L. (2005). Spore formation by Bacillus cereus in broth as affected by temperature, nutrient availability, and manganese. J. Food Prot. 68, 1734–1738. doi: 10.4315/0362-028x-68.8.1734
Samorodnitsky, D., and Nicholson, E. (2018). Differential effects of divalent cations on elk prion protein fibril formation and stability. Prion 12, 63–71. doi: 10.1080/19336896.2017.1423187
Schneider, J., Decamp, E., Koser, A., Fritz, S., Gonczi, H., Syversen, T., et al. (2006). Effects of chronic manganese exposure on cognitive and motor functioning in non-human primates. Brain Res. 1118, 222–231. doi: 10.1016/j.brainres.2006.08.054
Shu, J., Wu, H., Chen, M., Peng, H., Li, B., Liu, R., et al. (2019). Fractional removal of manganese and ammonia nitrogen from electrolytic metal manganese residue leachate using carbonate and struvite precipitation. Water Res. 153, 229–238. doi: 10.1016/j.watres.2018.12.044
Shuster, B., Khemmani, M., Nakaya, Y., Holland, G., Iwamoto, K., Abe, K., et al. (2019). Expansion of the spore surface polysaccharide layer in Bacillus subtilis by deletion of genes encoding glycosyltransferases and glucose modification enzymes. J. Bacteriol. 201:e00321-19. doi: 10.1128/JB.00321-19
Song, K., Wang, X., and Zhang, R. (2007). Culture conditions affecting spore production of Bacillus subtilis B68. Zhongguo Sheng Wu Fang Zhi 23, 255–259.
Su, J., Bao, P., Bai, T., Deng, L., Wu, H., Liu, F., et al. (2013). CotA, a multicopper oxidase from Bacillus pumilus WH4, exhibits manganese-oxidase activity. PLoS One 8:e60573. doi: 10.1371/journal.pone.0060573
Sunda, W., and Kieber, D. (1994). Oxidation of humic substances by manganese oxides yields low molecular-weight substrates. Nature 367, 62–64. doi: 10.1038/367062a0
Tang, W., Liu, H., and Zeng, X. (2021). Structural and functional study on cysteine 495, coordinating ligand to T1Cu site in multicopper oxidase CopA. Chemosphere 281:130807. doi: 10.1016/j.chemosphere.2021.130807
Tang, W., Wu, L., Gong, J., Ye, G., and Zeng, X. (2014). Screening, identification, and removal dynamics of a novel iron-manganese removal strain. Biomed. Mater. Eng. 24, 2049–2056. doi: 10.3233/BME-141014
Tebo, B. M., Johnson, H. A., McCarthy, J. K., and Templeton, A. S. (2005). Geomicrobiology of manganese(II) oxidation. Trends Microbiol. 13, 421–428. doi: 10.1016/j.tim.2005.07.009
Terahara, N., Namba, K., and Minamino, T. (2020). Dynamic exchange of two types of stator units in Bacillus subtilis flagellar motor in response to environmental changes. Comput. Struct. Biotechnol. J. 18, 2897–2907. doi: 10.1016/j.csbj.2020.10.009
Tian, Q., Xu, Y., Yu, C., Wang, J., and Ran, X. (2017). Identification of manganese-resistant bacteria and their physiological response under manganese stress. J. Mountain Agric. Biol. 36, 23–28. doi: 10.15958/j.cnki.sdnyswxb.2017.03.004
van Waasbergen, L., Hildebrand, M., and Tebo, B. (1996). Identification and characterization of a gene cluster involved in manganese oxidation by spores of the marine Bacillus sp. strain SG-1. J. Bacteriol. 178, 3517–3530. doi: 10.1128/jb.178.12.3517-3530.1996
Volpon, L., Young, C., Matte, A., and Gehring, K. (2006). NMR structure of the enzyme GatB of the galactitol-specific phosphoenolpyruvate-dependent phosphotransferase system and its interaction with GatA. Protein Sci. 15, 2435–2441. doi: 10.1110/ps.062337406
Wan, W., Xue, Z., Zhang, Z., and He, D. (2017). Manganese oxidation characteristics and oxidation mechanism of a manganese-oxidizing bacterium Arthrobacter sp. HW-16. Huan Jing Ke Xue 38, 2036–2043. doi: 10.13227/j.hjkx.201611103
Wang, X., Qin, X., Hao, Z., Luo, H., Yao, B., and Su, X. (2019). Degradation of four major mycotoxins by eight manganese peroxidases in presence of a dicarboxylic acid. Toxins 11:566. doi: 10.3390/toxins11100566
Wang, Y., He, F., You, J., Chen, K., Xu, B., and Lin, Y. (2013). Screening and identification of a manganese resistant strain and its manganese biosorption conditions. Xibei Nong Ye Xue Bao 22, 193–199.
Wang, Y., Xu, P., and Dai, X. (2016). Isolation and characterization of a manganese-oxidizing bacterium from soils. Wei Sheng Wu Xue Bao 56, 1699–1708.
Wang, Z., Wang, J., Liu, J., Chen, H., Li, M., and Li, L. (2017). Mechanistic insights into manganese oxidation of a soil-borne Mn(II)-oxidizing Escherichia coli strain by global proteomic and genetic analyses. Sci. Rep. 7:1352. doi: 10.1038/s41598-017-01552-3
Weisburg, W., Barns, S., Pelletier, D., and Lane, D. (1991). 16S ribosomal DNA amplification for phylogenetic study. J. Bacteriol. 173, 697–703. doi: 10.1128/jb.173.2.697-703.1991
Zeng, W., Li, F., Wu, C., Yu, R., Wu, X., Shen, L., et al. (2020). Role of extracellular polymeric substance (EPS) in toxicity response of soil bacteria Bacillus sp. S3 to multiple heavy metals. Bioprocess Biosyst. Eng. 43, 153–167. doi: 10.1007/s00449-019-02213-7
Zhang, W., He, L., Wang, Q., and Sheng, X. (2015). Inoculation with endophytic Bacillus megaterium 1Y31 increases Mn accumulation and induces the growth and energy metabolism-related differentially-expressed proteome in Mn hyperaccumulator hybrid pennisetum. J. Hazard. Mater. 300, 513–521. doi: 10.1016/j.jhazmat.2015.07.049
Keywords: Bacillus safensis, manganese oxidation, transcriptome, soil, adaptation
Citation: Ran X, Zhu Z, Long H, Tian Q, You L, Wu X, Liu Q, Huang S, Li S, Niu X and Wang J (2021) Manganese Stress Adaptation Mechanisms of Bacillus safensis Strain ST7 From Mine Soil. Front. Microbiol. 12:758889. doi: 10.3389/fmicb.2021.758889
Received: 15 August 2021; Accepted: 21 October 2021;
Published: 25 November 2021.
Edited by:
Aicha Asma Houfani, University of British Columbia, CanadaReviewed by:
Facundo Rodriguez Ayala, Institute of Nanoscience and Nanotechnology (CNEA), Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), ArgentinaKaty Juarez, Universidad Nacional Autónoma de México, Mexico
Copyright © 2021 Ran, Zhu, Long, Tian, You, Wu, Liu, Huang, Li, Niu and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jiafu Wang, jfwang@gzu.edu.cn