Highlights
The genome of the Bacillus megaterium strain BHS1 was fully sequenced.
Genes regulating metabolic activities to degrade halogenated compounds were characterized
Alkaliphiles can be reengineered into effective bioremediation system to treat halogen contaminated sites
Strain BHS1 genomic data offers the NCBI database with new genomic information.
This permits comparative genomic analysis studies with other B. megaterium species
Introduction
The soda lake, Mavi Gölü (Blue Lake), in Turkey is a famous tourist destination known for its uniquely beautiful turquoise water that is retained in the lakebed from June to December. It is formed in the Göksu Creek, the only carbonated water that flows into the Black Sea. Given the lake’s unique and rugged landscape, the identification/characterization of the lake’s culturable microbial communities is essential for expanding the repertoire of alkaliphilic microbes in fortifying available databases for alkaliphilic microbes. The data are particularly pertinent for basic research that elucidates various applications of alkaliphiles, such as in food industries, bioremediation, and medicine (Batumalaie et al. 2018; Bagherbaigi et al. 2013; Neelam et al. 2019; Kevbrin, 2019). Despite being unpolluted, Mavi Gölü may contain organobromines, which are produced naturally in the lake. The organobromine naturally produced by an array of the lake’s biological and chemical processes. The processes might trigger certain bacteria to produce dehalogenases (Gribble, 2000). Hence, we investigated the properties of dehalogenases produced by bacteria in heavily polluted areas as well as natural ecosystems that produce environmentally detrimental halogenated compounds.
Many bacteria utilize halogenated compounds as an energy source with the aid of dehalogenases (Abel et al. 2012; Adamu et al. 2016; Edbeib et al. 2016; Edbeib et al. 2017; Akcay and Kaya, 2019; Heidarrezaei et al. 2020; Oyewusi et al. 2020; Muslem et al. 2020). These bacteria generate dehalogenases that catalyze the cleavage of carbon-halogen bonds, which generally yields less toxic intermediates (Kurihara and Esaki, 2008). In particular, Rhizobium sp. RC1 remains to be the only bacterium known to produce three different dehalogenases, namely DehD, DehE, and DehL, each of which has unique substrate specificities (Adamu et al. 2016; Hamid et al. 2011; Huyop and Cooper, 2011; Huyop and Sudi, 2012; Sudi et al. 2014).
Because BHS1 produces a dehalogenase, it can utilize 2,2-dichloropropionic acid (2,2DCP) as a sole carbon source (Wahhab et al. 2020). Therefore, we sequenced the genome of strain BHS1 to understand the regulatory mechanisms for overexpression of dehalogenase genes and how dehalogenases contribute to BHS1 survival in an alkali-laden ecosystem. Previous literature reported that B. megaterium produces proteins of unknown function (Korneli et al. 2013) that could potentially be used for bioremediation. Thus, to consider its appropriate applications and understand its adaptability to alkali-laden environments such as soda lakes, the acquisition of a full genomic sequence for BSH1 is of interest.
Materials and methods
Sampling
During June and July 2019, water samples were collected at ~ 30 cm below the lakewater surface of Mavi Gölü (Dereli), Giresun, Black Sea Region, Turkey (latitude 40°38′08′′N north, longitude 38°22′56′′E east).
Preparation of pure genomic DNA
B. megaterium strain BHS1 was isolated from water samples taken from Mavi Gölü and grown aerobically in minimal medium containing 20 mM 2,2-DCP (pH 9.0) at 30 °C with rotation. After 18 h, the sample was harvested via centrifugation (13,000×g, 20 min, 4 °C). The QIAamp DNA Minikit (Qiagen, Germany) was used to extract and purify genomic DNA, the quality of which was validated using a Qubit 2.0. fluorometer (Apical Scientific, Malaysia). The purity of DNA samples (UVA260nm/A280nm) was assessed using a NanoDrop spectrophotometer (ThermoFisher).
Analysis of the 16S rRNA gene
Amplification of the 16S rRNA gene. Genomic DNA from strain BHS1 was extracted from bacterial cultures grown on minimal media containing 20 mM 2,2-DCP using the Wizard® Genomic DNA Purification kit. PCR was used to amplify the target DNA fragments using the universal primers fP1 (5′-AGAGTTTGATCCTGGCTCAG-3′) and rP1 (5′-ACGGTCATACCTTGTTACGACTT-3′) (Fulton and Cooper, 2005). PCR was carried out by 30 cycles of denaturation at 94 °C for 1 min and annealing at 55 °C for 1 min, with final extension at 72 °C for 10 min. Amplicons were purified using the QIAquick PCR purification kit (Qiagen, Germany) and sequenced by the 1st Base Laboratories Sdn Bhd. (Malaysia).
Evolutionary distances were computed using the Maximum Composite Likelihood method (Tamura et al. 2013), with units being the number of base substitutions per site. This analysis involved 15 sequences for bacterial 16S rDNA.
Genome sequencing, de novo assembly, and scaffolding
The genome of strain BHS1 was sequenced using the HiSeq platform. The insert library with ≤ 500 bp was prepared using standard protocols as stated by New England Biolabs Inc. and sequenced using the NGS-Illumina HiSeq TM 2000 Platform (Apical Scientific, Malaysia). All sequences were uploaded to the UCLA CNSI Hoffman 2 computer cluster lab located at the Malaysia Genome Institute for gene assembly. The read maps were based on the Pipeline of Bioinformatics Analysis afforded by Apical Scientific. Data processing, involving the original optical density data obtained by high-throughput sequencing (Illumina platform), were transformed into raw sequenced reads by CASAVA (Hosseini, et al. 2010) and stored in FASTQ format for both the read sequences and the corresponding sequencing-quality information. Quality control was performed, and low-quality sequences were removed.
De novo genome assembly was accomplished using Velvet assembler v1.2.10 at the Malaysia Genome Institute (Zerbino and Birney, 2008). High-quality, short paired-end Illumina reads were assembled using VelvetOptimiser.pl v2.2.5. An optimal parameter value of k-mer lengths was identified, and the best assembly was chosen based on N50, number of contigs, and assembly size values. Scaffolding was performed using SSPACE-Standard v3.0 (Boetzer et al. 2010). SSPACE revealed distance information for the paired-end Illumina reads to assess the order, distance, and orientation of the assembled contigs, which were combined into scaffolds. The assembled genome was further refined using the multi-draft-based scaffolder MeDuSa v1.6 (Bosi et al. 2015). This tool allows for the use of multiple reference genomes during scaffolding.
The statistics for final assembly were successfully obtained as one continuous scaffold with 5.37 Mb of total genome size, GC content 38%, and 0.05% gaps (contains N bases). For data accessibility, the B. megaterium BHS1 (L1) genome was deposited into NCBI under BioProject PRJNA637885. The Fasta genome submission was under accession number CP058255 (https://submit.ncbi.nlm.nih.gov).
Genome data analysis
Following genome assembly, genetic prediction or annotation was performed to identify DNA regions containing protein-coding genes. To determine the function of genes in the B. megaterium BHS1 genome, suitable software and databases, such as BLAST (Basic Local Alignment Search Tool), InterProScan, KEGG, and Blast2Go, were used. BLAST is a program that compares protein or DNA sequences with sequences listed inside various databases, i.e., nr, ref-seq, or SWISS-PROT.
The BLAST search tools used in this study were Blast2GO high-performance cloud server, BLAST Search, and the stand-alone NCBI-BLAST+ with DIAMOND high-performance analysis with a local database. BLAST analysis aided in finding sequence similarity with proteins in the non-redundant database GenBank with the default parameter set BLOSUM62. Automated annotation using Blast2GO version 5.2.5. (2018) [http://www.blast2go.com] was conducted to characterize the coding regions in B. megaterium strain BHS1. OmicsBox version 1.3 (https://www.biobam.com/omicsbox) was used for mapping, annotation, and visualization as well as the quantitative and statistical analyses for the scaffold of strain BHS1. This program was utilized because of its advanced features compared with Blast2GO (Götz et al., 2008).
Depending on the Kyoto Encyclopedia of Genes and Genomes (KEGG) (Galperin et al. 2014; Zhang et al. 2019) with BLASTP, several databases were utilized for genome annotation and gene function prediction, including Gene Ontology, Clusters of Orthologous Groups (COG), and “shn.” Java version 1.8.0_152 (2019), which is the OS version for Java 10.0 (http://java.sun.com), was used for accessing most databases. Comparison against other reference genomes was performed using a progressive Mauve tool; Mauve is a system for constructing multiple genome alignments that provide a basis for comparative genomics and the study of genome-wide evolutionary dynamics (Darling et al. 2010).
Whole-cell analysis of maximum chloride released from 2,2DCP, D,L-2-chloropropionic acid (D,L-2CP), L-2-chloropropionic acid (L-2CP), and D-2CP
B. megaterium BHS1 was incubated with individual selected substrates (2,2-DCP, D,L-2CP, L-2CP, and D-2CP) under specific growth conditions as previously described (Wahhab et al. 2020). A pure grade of the substrates above were purchased from Sigma-Aldrich or Merck (USA). The chloride ion released was determined colorimetrically for each substrate (Bergmann and Sanik, 1957). The color was allowed to develop for 10 min at 20 °C and measured at A460 nm. After incubation at 20 °C for 10 min, chloride-ion release for each substrate was determined colorimetrically (A460 nm) (Bergmann and Sanik, 1957).
Results and discussion
Organism information
The general features of strain BHS1 have been described by Wahhab et al. (2020). Its full 16S rRNA gene sequence was deposited in GenBank under accession number MT883351. The 16S rRNA gene sequence of BHS1 (GenBank accession number MT883351) was analyzed to determine its evolutionary relationships with other bacterial species. The final dataset contained 1591 bases. Evolutionary analyses were conducted with MEGA X (Kumar et al. 2018). As shown in Figure 1, BHS1 is closely related to B. megaterium strains NCT-2 and WSH-002. Hence, it was confirmed that BHS1 belongs to the species B. megaterium (homotypic synonym of Priestia megaterium) (Gupta et al. 2020) (Fig. 1).
Fig. 1 [Images not available. See PDF.]
Phylogenetic tree highlighting the relative position of the Bacillus genus with another closely related genus/species and strains. The evolutionary history was inferred using the Neighbor-Joining method (Saitou and Nei, 1987), where the optimal tree with the sum of branch length is shown to be 0.287. The percentages of replicate trees in which the associated taxa clustered together in the bootstrap test (500 replicates) were sited next to the branches (Felsenstein, 1985). The tree is drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree
Comparison with other complete reference genome sequences
A draft genome for BHS1 was generated based on a comparison with complete reference genomes of B. megaterium strains available in the NCBI database. The reference strains WSH-002 (NC_01738) and NCT-2 (NZ_CP032527) agreed with the 16S rRNA gene analysis. Figure 2 shows the alignment of the whole genome, where each genome is laid out horizontally. Homologous genome segments in the BHS1 draft genome are shown as colored blocks connected by lines to similarly colored blocks of the reference genomes for the other two strains. Regions that are entirely white could not be aligned and probably contain sequence elements that were recently acquired and thus are unique to the particular strain.
Fig. 2 [Images not available. See PDF.]
Comparison of assembled strain BHS1 draft genome (Sample L1.L350) against two complete reference genomes. Line 1: Strain BHS1 draft genome (Sample L1.L350), Line 2: B. megaterium strain WSH-002, and Line 3: B. megaterium strain NCT-2
Whole-genome sequencing information for strain BHS1
Analysis of one scaffold revealed that B. megaterium strain BHS1 has a 5,376,285-bp circular chromosome with an average G+C content of 38%. Analysis of the nuclear genome of strain BHS1 predicted ~ 5509 protein-coding genes, 79 tRNA genes, 67 rRNA genes, and 6 non-coding (mRNA) genes. Figure 3 shows the circular map of the genome, which was generated using CGview based on open reading frames with COG information. Table 1 summarizes other general features of the genome, and Table 2 lists the functional categories in COG for the full genome sequence of strain BHS1. The main categories are cell wall/membrane/envelope biogenesis (M) ~ 28.3%, transcription (K) 8.0%, translation, ribosomal structure and biogenesis (J) 3.6%, and transport and metabolism of nucleotides (F) ~3%.
Fig. 3 [Images not available. See PDF.]
Circular chromosome of B. megaterium strain BHS1 (from inner to outer circles)
Table 1. Genome features of B. megaterium strain BHS1
Features | Chromosome (value) | % of total |
|---|---|---|
Genome size (bp) | 5,376,285 bp (~ 5.37 Mb) | 100% |
Protein-coding genes (CDS). | 5509 | 100% |
Number of bases for both strand genes (CDS features) without a pseudo count. | 4,484,403 | 81.4% |
G+ C content (%) | - | 38% |
tRNA genes (1353) | 79 | 1.4% |
Number of non-coding (mRNA genes) | 6 | 0.1% |
Number of rRNA genes | 67 | 1.2% |
DNA scaffolds | 1 | – |
Number of gaps | 39 | 0.7% |
Repeats regions | 22 | 0.399% ~ 0.4% |
Pseudo genes | 9 | 0.16% |
Coding percentage (including introns) | 83.4 | 1.5% |
Table 2. Number of genes associated with general COG functional categories of the full genome sequences of B. megaterium strain BHS1
Code | Functional annotation | Strain BHS1 (number of genes) | % * |
|---|---|---|---|
J | Translation, ribosomal structure, and biogenesis | 200 | 3.6% |
K | Transcription | 443 | 8.0% |
L | Replication, recombination, and repair | 139 | 2.5% |
D | Cell cycle control, cell division, chromosome partitioning | 59 | 1.0% |
O | Posttranslational modification, protein turnover, chaperones | 57 | 1.0% |
M | Cell wall/membrane/envelope biogenesis | 1559 | 28.29%~ 28.3% |
N | Cell motility | 17 | 0.3% |
P | Inorganic ion transport and metabolism | 48 | 0.87 % |
T | Signal transduction mechanisms | 117 | 2.1% |
C | Energy production and conversion | 8 | 0.1% |
G | Carbohydrate transport and metabolism | 16 | 0.29%~0.3 % |
E | Amino acid transport and metabolism | 32 | 0.58 %~0.6 % |
F | Nucleotide transport and metabolism | 164 | 2.97 %~3% |
H | Coenzyme transport and metabolism | 16 | 0.29%~0.3% |
I | Lipid transport and metabolism | 46 | 0.8% |
Q | Secondary metabolites biosynthesis, transport, and catabolism | 47 | 0.85% |
S | Function unknown | 1902 | 34.5% |
- | No. IPS (not in COGs). | 469 | 8.5% |
A | RNA processing and modifications | 116 | 2.1% |
B | Chromatin structure and dynamics | 0 | 0% |
V | Defense mechanism | 1 | 0.018% |
U | Intracellular trafficking | 0 | 0% |
Total | Total number of genes | 5509 |
*The total based on the total number of genes in the full genome
With respect to potential biotechnology applications on the industrial aspect, the genome for BHS1 harbors several important industrial enzymes such as α-2 and β-amylase 1. Moreover, although plants, algae, fungi, and yeast have been the traditional primary sources for urease, this enzyme is also produced by strain BHS1. Other biotechnologically relevant enzymes encoded in the BHS1 genome include thirteen proteases, nine ureases, seven aminopeptidases, five lipases, four aminomutases, three serine proteases, three enzymes for l-lysine synthesis, two glucose dehydrogenases, two pullulanases, and one mutarotase. Thus, it is clear that alkaliphilic Bacillus strains are quite important and exciting to explore. This includes the production of biotechnologically important enzymes such as the annotated extracellular hydrolases.
Dehalogenase genes and dehalogenase subunits for adaptation of strain BHS1 to an alkaline environment
The germane finding of the genomic annotation for strain BHS1 was the identification of a metabolic pathway that includes two dehalogenase families that were inferred based on the ability of the cells to optimally grow on minimal medium containing 2,2-DCP as the sole carbon source, as reported by Wahhab et al. (2020).
Thus, genes linked to these two dehalogenases with a Locs tag were found in the genome of strain BHS1: the first gene encodes a Cof-type haloacid dehalogenase IIB family hydrolase (Nemati et al. 2013), and the second encodes a haloacid dehalogenase type II (Table 3).
Table 3. The genes code for enzymes of dehalogenase degradation pathway in strain BHS1
Enzyme | Locs tag |
|---|---|
Cof-type HAD-IIB family hydrolase | Bacillus_3231. WP_013082899.1 |
Haloacid dehalogenase type II | Bacillus_4574. WP_045293204.1 |
The existence of more than a single dehalogenase gene in a bacterial genome has been reported only for Rhizobium sp. RC1, and this bacterium has the ability to grow solely on 2,2-DCP and D,L-2-chloropropionate as the only sources of both carbon and energy (Allison et al. 1983; Leigh et al. 1986). To date, no study has explored the possibility that dehalogenase genes could contribute to adaptation to an alkali-laden environment (Ismail et al. 2017).
The 2-haloacid dehalogenases act on 2-haloacids, liberating a halide ion(s) and the corresponding 2-hydroxy acid. The 2-haloacid dehalogenases have been categorized phylogenetically into two groups, namely I and II. Group II, or L-type enzymes, include L-2-haloacid dehalogenases (i.e., L-DEXs) that explicitly target L-2-haloacids. The group II enzymes are more common than group I dehalogenases, which catalyze on both D and L-haloacids (Wang et al. 2018).
Co-transport mechanisms
Many cation/proton antiporters have been identified in bacteria (Krulwich et al. 2009; Tsujii et al. 2020). Cation/proton antiporters contribute to maintaining the appropriate intracellular pH homeostasis particularly in alkaline environments that require import of hydrogen ions from a comparatively proton-poor external environment. Three families of cation/proton antiporters of the Gene Ontology Membrane Transporter Group have been identified in bacteria and perform significant roles in maintaining cation equilibrium in the cytoplasm in extreme environmental settings. For instance, monovalent cation/proton antiporters are integral membrane proteins that facilitate the efflux of cytoplasmic sodium, potassium, or lithium ions in exchange for external hydrogen ions (protons) (Morino et al. 2014). The genome of strain BHS1 which was isolated from the alkaline lake, encodes all seven cation/proton antiporter subunits (A-G) in a single gene. These seven subunits with their Locs tag are shown in Table 4. Na+/H+ antiporters have numerous functions, namely the creation of an essential electrochemical gradient of Na+ across the plasma membrane for the purpose Na+-driven flagellar rotation, the efflux of Na+ and Li+ (both of which are toxic if accumulated at high concentrations within cells), and regulation of intracellular pH under alkaline conditions (Wai Liew et al. 2007). Interestingly, strain BHS1 encodes several types of Na+/H+ antiporters (Table 4) that might be vital for maintaining intracellular pH during growth in alkaline environments, e.g., pH 10.5, similar to the previously reported Na+/H+ antiporter gene (g1-nhaC) of alkaliphilic Bacillus sp. G1 when expressed in Escherichia coli (Wai Liew et al. 2007). The full genome sequence of Bacillus may not only reveal individual bacilli functions but also play a role in identifying unique properties that are essential for the survival and adaptation of bacilli in alkaline ecosystems.
Table 4. The genetic codes for enzymes of the Na+/H+ adaptation pathway of strain BHS1 in an alkaline environment
Enzyme | Locs tag |
|---|---|
Na+/H+ antiporter subunit A | Bacillus_354. WP_050688978.1 |
Na+/H+ antiporter subunit B | Bacillus_353. WP_013085301.1 |
Na+/H+ antiporter subunit C | Bacillus_352. WP_013059592.1 |
Na+/H+ antiporter subunit D | Bacillus_351. WP_028412608.1 |
Na+/H+ antiporter subunit E | Bacillus_350. WP_013059594.1 |
Na+/H+ antiporter subunit F1 | Bacillus_349. WP_013059595.1 |
Na+/H+ antiporter subunit G | Bacillus_348. WP_013085304.1 |
Organization of dehalogenase genes and a putative regulatory gene
The regulation of dehalogenase gene expression is poorly understood from the perspective of whole-genome function (Huyop and Cooper, 2011). Hence, analysis of the full genome of BHS1 may shed light on this aspect via the identification of an operon consisting of a functional gene cluster, i.e., regulatory gene and structural genes, the expression of which depends on their location within the operon. A putative haloacid dehalogenase operon comprised of L-type dehalogenase (group II haloacid dehalogenase or L-specific dehalogenases) was identified in the BHS1 genome [4,389,524_4,390,183; 659 bp], with downstream genes consisting of a hypothetical protein [4,390,313_4,390,738; 425 bp], a capsular biosynthesis protein [439,171_439,333;1562 bp], and the polyglutamate capsule biosynthesis protein CAPE [4,393,457_4,393,660; 203 bp]. Conversely, the upstream region of the operon encodes the glycoside hydrolase family 10 [4,387,790_4,388,125; 335 bp] and a helix-turn-helix (HTH)-type transcriptional regulator [4,385,958_4,386,812; 854 bp]. Thus, the acquisition of the full genomic sequence for BHS1 is crucial for determining the most appropriate applications of these gene clusters for biotechnology purposes (Fig. 4). Toward that end, we are currently using cloning analysis to decipher the regulatory mechanisms for these genes, especially given the existence of a HTH-type transcriptional regulator gene.
Fig. 4 [Images not available. See PDF.]
Physical map of putative L-2-haloacid dehalogenase gene (encoded putative DehLBSH1 protein) and its putative gene HTH-type transcriptional regulator together with other genes location
Analysis of a putative L-type dehalogenase (group II haloacid dehalogenase) and HTH-type transcriptional regulator
The size of L-type dehalogenase was 657 bp equivalent to 219 amino acid residues. The L-type dehalogenase of strain BHS1 was later designated as DehLBHS1. The program ProtParam (Gasteiger et al. 2005) indicated that DehLBHS1 has a molecular mass of 25,636.47 Da and a theoretical pI of 6.66, i.e., close to neutral. The sequence is comprised of 28% negatively charged residues (Asp, Glu) and 27% positively charged residues (Arg, Lys). The grand average of hydropathicity index (GRAVY) (Kyte and Doolittle, 1981) of −0.42 suggests that, overall, DehLBHS1 is a hydrophilic protein. Amino-acid sequence alignment of DehLBHS1 and L-DEX YL (Hisano et al. 1996; Li et al. 1998) revealed three key residues, namely D (Asp), R (Arg), and S (Ser), that may play key roles in the catalytic mechanism of DehLBHS1. The functional residues of L-DEX YL were determined to be D10, R41, S118, and D180, which are comparatively equivalent to the BHS1 residues D10, R41, S127, and D180, respectively (Fig. 5). Although the amino-acid sequences of the L-haloacid dehalogenases from BHS1 and Pseudomonas sp. YL (L-DEX YL) are only 46% identical, the key residues are 100% identical and thus may have similar functions. Analysis of the HTH-type transcriptional regulator, designated as DehRBHS1 (a BHS1 strain, dehalogenase regulator protein) revealed very low sequence identity (< 17%) compared with the well-established regulator protein (DehR) of Rhizobium sp. RC1 (data not shown). This suggests that DehRBHS1 may have a similar regulatory function in the dehalogenase operon system that controls dehalogenases.
Fig. 5 [Images not available. See PDF.]
The pairwise sequence alignment results of DehLBHS1 and L-DEX YL from Pseudomonas sp. strain YL (Nardi-Dei et al., 1997). Mark in bold were three key amino acids (D,Asp; R,Arg; S,Ser), which play an important role in the catalytic function of L-DEX YL haloacid dehalogenase. The three key amino acids (in bold) also were located in DehLBHS1 in the same position except for S(Ser). * symbol indicates a stop codon
Functional characterization of dehalogenase-producing strain BHS1
Dehalogenase-producing B. megaterium strain BHS1 was incubated with 5 mM of each substrate (2,2-DCP, D,L-2CP, L-2CP, and D-2CP) under proper growth conditions, and the maximum amount of chloride ion released was monitored. In the presence of 2,2-DCP or D,L-2CP, 50% of the chloride was released (Fig. 6). In the presence of L-2CP alone, essentially all chloride ions were released (4.5 μmol Cl–/ml), whereas, with d-2CP, no chloride ion was detected, as expected, suggesting the stereospecificity of the enzyme. A control experiment, L-2CP without any cells, showed no dehalogenation reaction, affirming the substrate’s non-autodegradation (Fig. 6).
Fig. 6 [Images not available. See PDF.]
Resolution of racemic D,L-2-chloropropionic acid, L-2CP, D-2CP by B. megaterium strain BHS1 expressing stereospecific DehLBHS1
Conclusions
The study performed here is the first full genome sequence analysis of the functional haloacid-degrading gene from the alkaliphilic B. megaterium strain BHS1 from Mavi Gölü with the potential to degrade haloalkanoic acid in an alkaline environment that amenable to converting haloacids for carbon and energy sources. The genomic data acquired during this study may facilitate the discovery of new dehalogenases and its regulatory functions. The knowledge may also provide insights into the bacterium’s genetic and metabolic regulatory pathways, of which will be vital for potential application of BHS1 to the bioremediation of contaminated alkali-laden ecosystems.
Acknowledgements
The authors would like to acknowledge the Scientific Technological Research Council of Turkey (TUBITAK), Program 2221-Fellowship for Visiting Scientists and Scientists on Sabbatical Leave with Grant number 1059B11800242. FH and BHAW thanks the Fundamental Research Grant Scheme (FRGS)-Ministry of Education R.J130000.7854.5F189/FRGS/1/2019/STG05/UTM/01/1 and Q.J130000.2414.08G59 for partial financial assistance. Finally, AHMS, BHAW, and HAO also thank Universiti Malaysia Sabah (UMS) for a joint university research project between the Faculty of Science UTM and the Faculty of Food Science and Nutrition UMS.
Authors’ contributions
All authors contributed equally to the work. The first author is a Ph.D. candidate and responsible for the data collections. FH is the main supervisor and principal investigator for this research. All materials collected from YK, the principal researcher at Ondokuz Mayis University, Turkey, were initially analyzed and prepared by MFE. NHS expert in genomic data analysis and HAO contributed partially to data analysis. KIRN and RAW read the overview of the manuscript and data organization. AAAH helped to perform the phylogenetic analysis. AHMS and MSMAN approved special funding related to the research work and publication. All authors read and approved the final manuscript.
Funding
Program 2221-Fellowship for Visiting Scientists and Scientists on Sabbatical Leave with Grant number 1059B11800242; (FRGS)-Ministry of Education R.J130000.7854.5F189/FRGS/1/2019/STG05/UTM/01/1; RUGQ.J130000.2414.08G59; University Malaysia Sabah Special Grants.
Availability of data and materials
Strictly for research without profit-making interest.
Ethics approval and consent to participate
Ethics approval and consent to participate: This article does not contain any studies with human participants or animals performed by any authors. Informed consent was obtained from all individual participants included in the study.
Consent for publication
All authors agreed to participate and agreed on the submission.
Competing interests
The authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
References
Abel, E; Pakingking, RV, Jr; Gregoria, G; Wint, MT; Huyop, F. Characteristics of dehalogenase from bacteria isolated from the Gut of Pond-reared Rohu (Labeo rohita) Juveniles in Myanmar. Adv Biosci Biotechnol; 2012; 3,
Adamu, A; Wahab, RA; Huyop, F. L-2-haloacid dehalogenase (DehL) from Rhizobium sp. RC1. SpringerPlus; 2016; 5,
Akcay, K; Kaya, Y. Isolation, characterization and molecular identification of a halotolerant Bacillus megaterium CTBmeg1 able to grow on halogenated compounds. Biotechnol Biotechnol Equip; 2019; 33,
Allison, N; Skinner, AJ; Cooper, RA. The dehalogenases of a 2,2-dichloropropionate-degrading bacterium. J Gen Microbiol; 1983; 129, pp. 1283-1293.[COI: 1:CAS:528:DyaL3sXks1Ogs70%3D]
Bagherbaigi, S; Gicana, RG; Lamis, RJ; Nemati, M; Huyop, F. Characterisation of Arthrobacter sp. S1 that can degrade α and β-haloalkanoic acids isolated from contaminated soil. Ann Microbiol; 2013; 63, pp. 1363-1369.[COI: 1:CAS:528:DC%2BC3sXhvVahtb%2FF] [DOI: https://dx.doi.org/10.1007/s13213-012-0595-4] https://doi.org/10.1007/s13213-012-0595-4
Batumalaie, K; Edbeib, MF; Mahat, NA; Huyop, F; Wahab, RA. In silico and empirical approaches toward understanding the structural adaptation of the alkaline-stable lipase KV1 from Acinetobacter haemolyticus. J Biomol Struct Dynamics; 2018; 36,
Bergmann, J; Sanik, JJ. Determination of trace amounts of chlorine in naphtha. Anal Chem; 1957; 29,
Boetzer, M; Henkel, CV; Jansen, HJ; Butler, D; Pirovano, W. Scaffolding pre-assembled contigs using SSPACE. Bioinformatics; 2010; 27,
Bosi, E; Donati, B; Galardini, M; Brunetti, S; Sagot, M-F; Lió, P; Crescenzi, P; Fani, R; Fondi, M. MeDuSa: a multi-draft based scaffolder. Bioinformatics; 2015; 31,
Darling, AE; Mau, B; Perna, NT. Progressive Mauve: multiple genome alignment with gene gain, loss and rearrangement. PloS One; 2010; 5,
Edbeib, MF; Wahab, RA; Huyop, F. Characterization of an α-haloalkanoic acid–degrading Pseudomonas aeruginosa MX1 isolated from contaminated seawater. Bioremed J; 2016; 20,
Edbeib, MF; Wahab, RA; Kaya, Y; Huyop, F. In silico characterization of a novel dehalogenase (DehHX) from the halophile Pseudomonas halophila HX isolated from Tuz Gölü Lake, Turkey: insights into a hypersaline-adapted dehalogenase. Ann Microbiol; 2017; 67,
Felsenstein, J. Confidence limits on phylogenies: an approach using the bootstrap. Evolution; 1985; 39,
Fulton, CK; Cooper, RA. Catabolism of sulfamate by Mycobacterium sp. CF1. Environ Microbiol; 2005; 7, pp. 378-381.[COI: 1:CAS:528:DC%2BD2MXisV2qur0%3D] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15683398][DOI: https://dx.doi.org/10.1111/j.1462-2920.2005.00719.x] https://doi.org/10.1111/j.1462-2920.2004.00719.x
Galperin, MY; Makarova, KS; Wolf, YI; Koonin, EV. Expanded microbial genome coverage and improved protein family annotation in the COG database. Nucleic Acids Res; 2014; 43,
Gasteiger E, Hoogland C, Gattiker A, Duvaud S, Wilkins MR, Appel RD, Bairoch A (2005) Protein identification and analysis tools on the ExPASy Server. In: Walker JM (ed) The Proteomics Protocols Handbook. Humana Press, pp 571–607
Götz, S; Garcia-Gomez, JM; Terol, J; Williams, TD; Nagaraj, SH; Nueda, MJ; Robles, M; Talon, M; Dopazo, J; Conesa, A. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic acids Res; 2008; 36,
Gribble, GW. The natural production of organobromine compounds. Environ Sci Poll Res; 2000; 7,
Gupta, RS; Patel, S; Saini, N; Chen, S. Robust demarcation of 17 distinct bacillus species clades, proposed as novel bacillaceae genera, by phylogenomics and comparative genomic analyses: description of Robertmurraya Kyonggiensis sp. Nov. And Proposal for an Emended Genus Bacillus Limiting It Only to the Members of the Subtilis and Cereus Clades of Species. Int J Syst Evol Microbiol; 2020; 70,
Hamid, THTA; Hamid, AAA; Huyop, F. A review on non-stereospecific haloalkanoic acid dehalogenases. Afr J Biotechnol; 2011; 10,
Heidarrezaei, M., Shokravi, H., Huyop, F., Rahimian Koloor, S.S., Petrů, M. (2020) Isolation and characterization of a novel bacterium from the marine environment for trichloroacetic acid bioremediation. Appl Sci, 10(13):4593. https://doi.https://doi.org/10.3390/app10134593
Hisano, T; Hata, Y; Fujii, T; Liu, J-Q; Kurihara, T; Esaki, N; Soda, K. Crystal structure of L-2-haloacid dehalogenase from Pseudomonas sp. YL AN α/β hydrolase structure that is different from the α/β hydrolase fold. J Biol Chem; 1996; 271,
Hosseini, P; Tremblay, A; Matthews, BF; Alkharouf, NW. An efficient annotation and gene-expression derivation tool for Illumina Solexa datasets. BMC Res notes; 2010; 3, 183. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20598141][PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2908109][DOI: https://dx.doi.org/10.1186/1756-0500-3-183] [COI: 1:CAS:528:DC%2BC3cXos12mu7s%3D]https://doi.org/10.1186/1756-0500-3-183
Huyop, F; Cooper, RA. Regulation of dehalogenase E (Dehe) and expression of dehalogenase regulator gene ( Dehr ) from Rhizobium sp. RC1 in E. coli. Biotechnol Biotechnol Equip; 2011; 25,
Huyop, F; Sudi, IY. D-specific dehalogenases, a review. Biotechnol Biotechnol Equip; 2012; 26,
Ismail SNF, Wahab RA, Huyop F (2017) Microbial isolation and degradation of selected haloalkanoic aliphatic acids by locally isolated bacteria: A review. Malaysian J Microbiol 13(3):261–272
Kevbrin, VV. Mamo, G; Mattiasson, B. Isolation and Cultivation of Alkaliphiles. Alkaliphiles in Biotechnology. Advances in Biochemical Engineering/Biotechnology; 2019; Cham, Springer: pp. 195-220. https://doi.org/10.1007/10_2018_84
Korneli, C; David, F; Biedendieck, R; John, D; Wittmann, C. Getting the big beast to work—Systems biotechnology of Bacillus megaterium for novel high-value proteins. J Biotechnol; 2013; 163,
Krulwich, TA; Hicks, DB; Ito, M. Cation/proton antiporter complements of bacteria: why so large and diverse?. Mol Microbiol; 2009; 74,
Kumar, S; Stecher, G; Li, M; Knyaz, C; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol; 2018; 35,
Kurihara, T; Esaki, N. Bacterial hydrolytic dehalogenases and related enzymes: occurrences, reaction mechanisms, and applications. Chem Record.; 2008; 8,
Leigh, JA; Skinner, AJ; Cooper, RA. Isolation and partial characterization of dehalogenase-deficient mutants of a Rhizobium sp. FEMS Microbiol Letters; 1986; 36, pp. 136-166. [DOI: https://dx.doi.org/10.1111/j.1574-6968.1986.tb01688.x]
Li, Y-F; Hata, Y; Fujii, T; Hisano, T; Nishihara, M; Kurihara, T; Esaki, N. Crystal structures of reaction intermediates ofl-2-haloacid dehalogenase and implications for the reaction mechanism. J Biol Chem; 1998; 273,
Morino, M; Suzuki, T; Ito, M; Krulwich, TA. Purification and functional reconstitution of a seven-subunit Mrp-Type Na+/H+ Antiporte. J Bacteriol; 2014; 196,
Muslem, WH; Edbeib, MF; Aksoy, HM; Kaya, Y; Hamid, AAA; Hood, MHM; Wahab, RA; Huyop, F. Biodegradation of 3-chloropropionic acid (3-CP) by Bacillus cereus WH2 and its in silico enzyme-substrate docking analysis. J Biomol Struct Dynamics; 2020; 38,
Nardi-Dei, V; Kurihara, T; Chung, P; Esaki, N; Soda, K. Bacterial D,L-2-haloacid dehalogenase from Pseudomonas sp. strain 113: Gene cloning and structural comparison with D and L-2 haloacid dehalogenases. J Bacteriol; 1997; 179,
Neelam, DK; Agrawal, A; Tomer, AK; Bandyopadhayaya, S; Sharma, A; Jagannadham, MV; Mandal, CC; Dadheech, PK. A Piscibacillus sp. isolated from a soda lake exhibits anticancer activity against breast cancer MDA-MB-231 cells. Microorganisms; 2019; 7,
Nemati, M., Edbeib, M.F., Hamid, A.A.A., Gicana, R.G., Ibrahim, N., Huyop, F. (2013) Identification of putative Cof-like hydrolase associated with dehalogenase in Enterobacter cloacae MN1 isolated from the contaminated sea-side area of the Philippines. Malaysian J Microbiol, 9(3):253-259. https://doi: https://doi.org/10.21161/mjm.51113
Oyewusi, HA; Wahab, RA; Kaya, Y; Edbeib, MF; Huyop, F. Alternative bioremediation agents against haloacids, haloacetates and chlorpyrifos using novel halogen-degrading bacterial isolates from the hypersaline lake Tuz. Catalysts; 2020; 10,
Saitou, N; Nei, M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol; 1987; 4, pp. 406-425.[COI: 1:STN:280:DyaL1c7ovFSjsA%3D%3D]https://doi.org/10.1093/oxfordjournals.molbev.a040454
Sudi, IY; Hamid, AAA; Shamsir, MS; Jamaluddin, H; Wahab, RA; Huyop, F. Insights into the stereospecificity of the D-specific dehalogenase from Rhizobium sp. RC1 toward d- and l-2-chloropropionate. Biotechnol Biotechnol Equip; 2014; 28,
Tamura, K; Stecher, G; Peterson, D; Filipski, A; Kumar, S. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol Biol Evol; 2013; 30,
Tsujii, M; Tanudjaja, E; Uozumi, N. Diverse Physiological Functions of Cation Proton Antiporters across Bacteria and Plant Cells. Int J Mol Sci; 2020; 21,
Wahhab, BHA; Anuar, NFSK; Wahab, RA; Al Nimer, MS; Samsulrizal, NH; Hamid, AAA; Edbeib, MF; Kaya, Y; Huyop, F. Identification and characterization of a 2, 2-dichloropropionic acid (2, 2-DCP) degrading alkalotorelant bacterium strain BHS1 isolated from Blue Lake, Turkey. J Trop Life Sci; 2020; 10,
Wai Liew, C; Illias, RM; Muhammad Mahadi, N; Najimudin, N. Expression of the Na+/H+ antiporter gene (g1-nhaC) of alkaliphilic Bacillus sp. G1 in Escherichia coli. FEMS Microbiol Letters; 2007; 276,
Wang, Y; Feng, Y; Cao, X; Liu, Y; Xue, S. Insights into the molecular mechanism of dehalogenation catalyzed by D-2-haloacid dehalogenase from crystal structures. Sci Rep; 2018; 8,
Zerbino, DR; Birney, E. Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res; 2008; 18,
Zhang, H; Dong, S; Lou, T; Wang, S. Complete genome sequence unveiled cellulose degradation enzymes and secondary metabolic potentials in Streptomyces sp. CC0208. J Basic Microbiol; 2019; 59,
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© The Author(s) 2021. This work is published under http://creativecommons.org/licenses/by/4.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
Purpose
Bacillus megaterium strain BHS1, isolated from an alkaline water sample taken from Mavi Gölü (Blue Lake, Turkey), can grow on minimal medium containing 2,2-dichloropropionic acid. We characterized this bacterium at the genomic level.
Methods
The HiSeq platform was used to carry out genome sequencing, de novo assembly, and scaffolding with strain BHS1. Next, genome data were analyzed to demarcate DNA regions containing protein-coding genes and determine the function of certain BHS1 genes. Finally, results from a colorimetric chloride ion–release assay demonstrated that strain BHS1 produces dehalogenase.
Results
De novo assembly of the BHS1 genomic sequence revealed a genome size of ~ 5.37 Mb with an average G+C content of 38%. The predicted nuclear genome harbors 5509 protein-coding genes, 1353 tRNA genes, 67 rRNA genes, and 6 non-coding (mRNA) genes. Genomic mapping of strain BHS1 revealed its amenability to synthesize two families of dehalogenases (Cof-type haloacid dehalogenase IIB family hydrolase and haloacid dehalogenase type II), suggesting that these enzymes can participate in the catabolism of halogenated organic acids. The mapping identified seven Na+/H+ antiporter subunits that are vital for adaptation of the bacterium to an alkaline environment. Apart from a pairwise analysis to the well-established L-2-haloacid dehalogenases, whole-cell analysis strongly suggested that the haloacid dehalogenase type II might act stereospecifically on L-2-chloropropionic acid, D,L-2-chloropropionic acid, and 2,2-dichloropropionic acid. Whole-cell studies confirmed the utilization of these three substrates and the gene’s role in dehalogenation.
Conclusions
To our knowledge, this is the first report of the full genome sequence for strain BHS1, which enabled the characterization of selected genes having specific metabolic activities and their roles in the biodegradation of halogenated compounds.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details
1 Faculty of Medicine, Al-Mustansiriyah University, Department of Microbiology, Baghdad, Iraq (GRID:grid.411309.e); Universiti Teknologi Malaysia, Department of Biosciences, Faculty of Science, Johor Bahru, Malaysia (GRID:grid.410877.d) (ISNI:0000 0001 2296 1505)
2 International Islamic University Malaysia, Department of Plant Science, Kulliyah of Science, Kuantan, Malaysia (GRID:grid.440422.4) (ISNI:0000 0001 0807 5654)
3 Bani Walid University, Department of Medical Laboratories, Faculty of Medical Technology, Bani Walid, Libya (GRID:grid.440422.4)
4 Universiti Teknologi Malaysia, Department of Chemistry, Faculty of Science, Johor Bahru, Malaysia (GRID:grid.410877.d) (ISNI:0000 0001 2296 1505)
5 University Medical Hawler, Department of Pharmacology, Faculty of Pharmacy, Erbil, Iraq (GRID:grid.410877.d)
6 International Islamic University Malaysia, Department of Biotechnology, Kulliyah of Science, Kuantan, Malaysia (GRID:grid.440422.4) (ISNI:0000 0001 0807 5654)
7 Universiti Teknologi Malaysia, Department of Biosciences, Faculty of Science, Johor Bahru, Malaysia (GRID:grid.410877.d) (ISNI:0000 0001 2296 1505)
8 Kyrgyz-Turkish Manas University, Department of Biology, Faculty of Science, Bishkek, Kyrgyzstan (GRID:grid.444269.9); Ondokuz Mayis University, Department of Agricultural Biotechnology, Faculty of Agriculture, Samsun, Turkey (GRID:grid.411049.9) (ISNI:0000 0004 0574 2310)
9 University of Santo Tomas España, Faculty of Medicine and Surgery, Manila, Philippines (GRID:grid.412775.2) (ISNI:0000 0004 1937 1119)
10 Universiti Malaysia Sabah, Faculty of Food Science and Nutrition, Kota Kinabalu, Malaysia (GRID:grid.265727.3) (ISNI:0000 0001 0417 0814)