3.1. Isolation and Morphological Characterization of Bacillus altitudinis DG4
The bacterial isolate obtained from the Damanganga River water samples displayed distinct colony morphologies on the nutrient agar plates. The colonies were medium to large (2–5 mm in diameter), circular with entire (smooth) edges, and exhibited raised or umbonate elevation. The surface appeared smooth and glistening with a moist or mucoid appearance. The colonies demonstrated a buttery or creamy texture with a soft and opaque consistency. The color ranged from cream to off-white or pale yellow with an opaque appearance and entire or smooth margins (Table 1).
Table 1
Morphological characteristics of Bacillus altitudinis DG4 colonies on nutrient agar.
| Morphology | Characteristic |
| Size | Medium to large in size |
| Diameter | 2–5 mm |
| Shape | Circular in shape with entire (smooth) edges |
| Elevation | Raised or umbonate elevation |
| Surface | Smooth and glistening, with a moist or mucoid appearance |
| Texture | Buttery or creamy |
| Consistency | Soft and opaque |
| Color | Cream, off-white or pale yellow |
| Opacity | Opaque |
| Margin | Entire or smooth |
Biochemical tests were performed on Bacillus altitudinis DG4, which tested positive for catalase, methyl red, and citrate utilization, but negative for oxidase, indole production, and Voges-Proskauer tests. These results indicate that B. altitudinis DG4 possesses catalase activity, can utilize citrate as a sole carbon source, and follows the mixed acid fermentation pathway (Methyl Red positive), but does not produce oxidase, indole, or acetoin. Gram staining of the pure isolate confirmed its classification as a Gram-positive bacterium.
These observations are consistent with the known attributes of Bacillus species, which frequently demonstrate noticeable morphological traits that facilitate identification [25].
Existing knowledge regarding Bacillus altitudinis indicates its significance in various biotechnological and environmental applications. Studies have established that members of this genus possess diverse metabolic capabilities that are useful for industrial microbiology and bioremediation [26].
In comparison to other Bacillus species, B. altitudinis DG4 exhibited unique biochemical characteristics, notably catalase positivity and oxidase negativity. This finding aligns with the general classifications of the genus, where catalase activity is a common trait among Bacillus strains. Conversely, the absence of indole production in B. altitudinis DG4 presents a notable differentiation from some other Bacillus species, which are typically indole-positive [27]. These biochemical traits are important for further understanding metabolic pathways and ecological adaptations.
B. altitudinis DG4 possesses catalase activity and can utilize citrate as the sole carbon source, reinforcing its classification within the Bacillus genus. Such abilities suggest its potential biotechnological applications, particularly for degrading environmental pollutants or serving as a biocontrol agent in agriculture contexts [28].
3.2. Naphthalene Degradation by Bacillus altitudinis DG4
Among the five isolates, strain DG5 demonstrated in Table 2 the highest growth rate with an OD₆₀₀ value of 0.6070, indicating robust proliferation in the presence of naphthalene. However, the highest naphthalene degradation efficiency was observed for strain DG4, which removed 90.0% of the initial naphthalene concentration, despite exhibiting a moderate growth rate (OD₆₀₀ = 0.4241). This suggests that strain DG4 possesses highly efficient naphthalene-degrading enzymatic machinery.
Table 2
Growth rates and naphthalene degradation efficiencies of bacterial isolates after 7 days of incubation
| No. | Organism | The growth rate (OD at λ600) | Napthelene degradation efficiency (%) |
| 1 | DG1 | 0.4092 | 80% |
| 2 | DG2 | 0.3476 | 79% |
| 3 | DG3 | 0.4343 | 82% |
| 4 | DG4 | 0.4241 | 90% |
| 5 | DG5 | 0.607 | 85% |
Strain DG3 showed the second-highest degradation efficiency (82.0%) with a growth rate (OD₆₀₀ = 0.4343) comparable to DG4. Strain DG5, despite having the highest growth rate, achieved a naphthalene degradation efficiency of 85.0%, ranking third among the isolates. Strains DG1 and DG2 exhibited relatively lower degradation efficiencies of 80.0% and 79.0%, respectively, with DG2 showing the lowest growth rate (OD₆₀₀ = 0.3476) among all isolates.
Supporting this, the data presented in this study confirm that strain DG4 operates efficiently, even with a moderate growth profile, focusing its metabolic resources on the degradation of naphthalene [29]. These findings pave the way for further investigations focused on the practical application of this strain in the bioremediation of naphthalene-contaminated environments, although there remains a limitation due to the absence of in situ field trials to corroborate laboratory results under natural environmental conditions [25].
3.3. Genome Assembly and Annotation of Bacillus altitudinis DG4
Whole-genome sequencing of B. altitudinis DG4 yielded a high-quality genome assembly. The assembled genome consisted of 19 contigs with a total length of 3,831,796 base pairs (bp) and an average G + C content of 41.03%. The genome assembly displayed a Contig L50 of 2 and a Contig N50 of 987,862 bp, indicating robust assembly quality (Table 3). Plasmids were not identified within the genome. The complete genome sequence of Bacillus altitudinis strain DG4 is available in the European Nucleotide Archive (ENA) under project accession number PRJEB76590.0.
Table 3
Assembly statistics of Bacillus altitudinis DG4 genome
| Assembly Parameter | Value |
| Contigs | 19 |
| GC Content | 41.03% |
| Plasmids | 0 |
| Contig L50 | 2 |
| Genome Length | 3,831,796 bp |
| Contig N50 | 987,862 |
| Chromosomes | 0 |
Genome annotation using the RAST toolkit (RASTtk) assigned a unique genome identifier (293387.487) to B. altitudinis DG4. The taxonomic classification placed the organism in the superkingdom Bacteria, phylum Bacillota, class Bacilli, order Bacillales, family Bacillaceae, genus Bacillus, and genus Bacillus altitudinis.
Genome annotation revealed 4,120 protein-coding sequences (CDS), 66 transfer RNA (tRNA) genes, and 3 ribosomal RNA (rRNA) genes (Table 4). Among the protein-coding sequences, 1,017 were annotated as hypothetical proteins and 3,103 were assigned functional annotations. Of these, 942 proteins received Enzyme Commission (EC) number assignments, 787 had Gene Ontology (GO) assignments, and 695 were mapped to KEGG pathways. Additionally, 3,819 proteins were classified into PATRIC genus-specific protein families (PLFams), and 3,924 proteins were assigned to PATRIC cross-genus protein families (PGFams).
Table 4
Annotated genome features of Bacillus altitudinis DG4
| Feature Type | Count |
| CDS | 4,120 |
| tRNA | 66 |
| Repeat Regions | 5 |
| rRNA | 3 |
| Partial CDS | 0 |
| Miscellaneous RNA | 0 |
| Hypothetical proteins | 1,017 |
| Proteins with functional assignments | 3,103 |
| Proteins with EC number assignments | 942 |
| Proteins with GO assignments | 787 |
| Proteins with Pathway assignments | 695 |
| Proteins with PATRIC genus-specific family (PLfam) assignments | 3,819 |
| Proteins with PATRIC cross-genus family (PGfam) assignments | 3,924 |
3.4. Subsystem Analysis and Specialty Genes
The complete genomic features of B. altitudinis DG4 was elucidated and is preseanted as a circular genome map in Fig. 2, whereby showing the arrangement of coding sequences, RNA genes, and other genomic features. Subsystem analysis of the B. altitudinis DG4 genome revealed a distribution of functional categories across the genome (Fig. 3). The analysis identified various functional subsystems including those related to metabolism, cellular processes, and stress responses.
Genome analysis also identified several specialty genes, including those homologous to known transporters, virulence factors, drug targets, and antibiotic resistance genes (Table 5). Specifically, 63 genes showed homology to transporters in the Transporter Classification Database (TCDB), 41 genes were associated with antibiotic resistance in the PATRIC database, and 21 genes were identified as potential drug targets in DrugBank.
Table 5
Specialty genes identified in the Bacillus altitudinis DG4 genome.
| Specialty Category | Source | Genes |
| Antibiotic Resistance | CARD | 2 |
| Antibiotic Resistance | NDARO | 2 |
| Antibiotic Resistance | PATRIC | 41 |
| Drug Target | DrugBank | 21 |
| Drug Target | TTD | 1 |
| Transporter | TCDB | 63 |
| Virulence Factor | PATRIC_VF | 2 |
| Virulence Factor | Victors | 2 |
Previous studies have demonstrated the significance of Bacillus altitudinis in adapting to high-altitude environments, hinting at the genetic underpinnings that facilitate resilience. The role of specific genetic factors, including metabolic pathways involved in stress tolerance, underscores the ecological versatility [27].
Compared with other Bacillus species, B. altitudinis DG4 possesses a distinctive set of transporters and antibiotic resistance genes. Specifically, our analysis identified a significant number of genes associated with these traits, indicating functional niches that may not be present in closely related species. However, comprehensive comparative genomic analyses of specific reference genes are necessary to establish these distinctions [30].
Our detailed analysis revealed that 63 genes showed homology to transporters documented in the Transporter Classification Database (TCDB), aligning with metabolic and transport functionalities essential for nutrient acquisition. Additionally, 41 genes associated with antibiotic resistance were identified in the PATRIC database, whereas 21 genes designated as potential drug targets were identified in DrugBank. These findings provide substantial evidence for the enhanced functional capacity of B. altitudinis DG4 [31], [32].
This study provides novel insights into the genetic landscape of B. altitudinis DG4, particularly concerning its specialized genes and functional subsystems that contribute to its environmental adaptability and potential uses in biotechnology [33].
3.5. Antimicrobial Resistance Gene Analysis
Genome annotation identified several genes associated with antimicrobial resistance (Table 6). These included genes encoding antibiotic inactivation enzymes (CatA6 family, FosB), antibiotic target protection proteins (BcrC), antibiotic target replacement proteins (fabL), and efflux pumps conferring antibiotic resistance (BceA, BceB, EbrA, and EbrB).
Table 6
Antimicrobial resistance genes identified in Bacillus altitudinis DG4.
| AMR Mechanism | Genes |
| Antibiotic inactivation enzyme | CatA6 family, FosB |
| Antibiotic target in susceptible species | Alr, Ddl, dxr, EF-G, EF-Tu, folA, Dfr, folP, gyrA, gyrB, inhA, fabI, Iso-tRNA, kasA, MurA, rho, rpoB, rpoC, S10p, S12p |
| Antibiotic target protection protein | BcrC |
| Antibiotic target replacement protein | fabL |
| Efflux pump conferring antibiotic resistance | BceA, BceB, EbrA, EbrB |
| Gene conferring resistance via absence | gidB |
| Protein altering cell wall charge conferring antibiotic resistance | GdpD, MprF, PgsA |
| Regulator modulating expression of antibiotic resistance genes | BceR, BceS, LiaF, LiaR, LiaS |
3.6. Comparative Genome Analysis
This analysis further confirmed the taxonomic classification of the isolate and demonstrated its close relationship with other members of the Bacillus pumilus group, while showing clear separation from other Bacillus species such as B. anthracis and B. psychrosaccharolyticus. As shown in Fig. 4, this analysis positioned B. altitudinis DG4 within the Bacillus altitudinis clade, highlighted in red in the phylogenetic tree. Phylogenetic analysis based on the PATRIC global protein families (PGFams) positioned B. altitudinis DG4 within the Bacillus altitudinis clade (Fig. 5). This analysis further confirmed the taxonomic classification of these isolates.
3.7. Pan-Genome Analysis
The pan-genome analysis of B. altitudinis DG4, conducted using the Integrated Prokaryotic Genome and Pan-genome Analysis (IPGA) service with 122 B. altitudinis genomes, revealed a comprehensive view of the gene content and variability within this strain (Fig. 6). The analysis identified 2,403 core genes that were present across all genomes, indicating their essential roles in basic cellular functions and survival. Additionally, 1,588 accessory genes were identified, which were present only in a subset of the genomes, suggesting their involvement in adaptive functions. Notably, 40 unique genes were found to be specific to B. altitudinis DG4, potentially indicating recent acquisitions through horizontal gene transfer or unique adaptations.
Analysis of accessory genes in B. altitudinis DG4 revealed genes involved in diverse metabolic pathways, such as accA_1 (propanoate metabolism), accC_2 (biotin carboxylase), and acsA_1 (acetyl-coenzyme A synthetase). Notably, genes associated with degradation pathways, including acyP (1,4-dichlorobenzene degradation) and dhaT (naphthalene and anthracene degradation), were identified, suggesting their potential xenobiotic degradation capabilities. Additionally, genes for various transport systems (btuD_1 and btuD_6) and transcriptional regulators (gabR_1 and rhaS_4) were present, indicating complex gene regulation mechanisms and nutrient acquisition systems (Fig. 6).
Functional Analysis successfully identified key accessory genes, including accA_1, accC_2, and acsA_1, which are implicated in a range of metabolic pathways. Understanding these genes is vital because accessory genes play a significant role in bacterial adaptability and survival in various environments [34]. This insight underscores the potential application of B. altitudinis DG4 in bioremediation, where environmental pollutants pose challenges to health and ecosystems [35].
Previous studies have established that accessory genes contribute to the metabolic diversity of bacteria, enabling them to exploit a broader spectrum of substrates. In comparison to other Bacillus species, B. altitudinis DG4 has shown a unique set of accessory genes connected to xenobiotic degradation, specifically acyP and dhaT. These genes may provide a strain with enhanced capacities for utilizing and degrading environmental pollutants, contradicting earlier findings that suggested limited xenobiotic degradation abilities within Bacillus species [36], [37]. The data presented in this study indicate promising metabolic versatility in B. altitudinis DG4, opening avenues for its utilization in bioremediation strategies.
3.8. Analysis of Residual Naphthalene with GC-MS by B. altitudinis DG4 (continued)
Gas Chromatography-Mass Spectrometry (GC-MS) analysis of naphthalene degradation by B. altitudinis DG4 confirmed the metabolism of naphthalene and formation of intermediate metabolites. The major metabolites were 1-naphthalenol (18.21 min), 1,8-naphthalic anhydride (19.48 min), 1-acenaphthanone (20.31 min), and benzoic acid (7.70 min) (Fig. 7). The detection of these metabolites suggests that B. altitudinis DG4 employs a specific naphthalene degradation pathway, likely involving initial methylation followed by carboxylation.
The degradation of naphthalene by B. altitudinis DG4 requires a detailed understanding of the metabolic byproducts generated during the process, especially using techniques such as Gas Chromatography-Mass Spectrometry (GC-MS). The findings from our GC-MS analysis identified several key metabolites, including 1-naphthalenol and benzoic acid, indicating the significant metabolic activity of naphthalene by B. altitudinis DG4. This aligns with existing literature suggesting naphthalene's metabolic pathways in various bacterial strains, particularly that certain metabolites can be produced during the degradation of naphthalene by microbes [38], [39].
The naphthalene degradation pathway, as illustrated in Fig. 8, was investigated through the identification of key metabolites using chromatographic and spectroscopic analysis. The primary intermediates detected at different retention times included 1-naphthalenol, 1,8-naphthalic anhydride, 1-acenaphthanone, and benzoic acid.
These intermediates suggested a sequential degradation pathway involving hydroxylation, oxidation, and ring cleavage. The presence of benzoic acid as a degradation product indicates the possibility of its further breakdown into simpler aliphatic acids, which can enter the tricarboxylic acid (TCA) cycle for complete mineralization.
The experimental findings confirmed the successful degradation of naphthalene through microbial and oxidative processes, with the formation of key intermediates leading to the mineralization of the parent compound. The identification of 1,8-naphthalic anhydride and 1-acenapthol suggests an alternative degradation pathway involving anhydride formation before ring cleavage.
These results align with those of previous studies on naphthalene biodegradation by Bordetella avium and Paraburkholderia aromaticivorans, which identified similar intermediates and confirmed the role of dioxygenase enzymes in initiating naphthalene catabolism [38]. Zeng et al. (2024) also supported this degradation pathway, demonstrating the role of radical oxidation in producing naphthalenic intermediates [21].
Elucidating this degradation pathway is important for bioremediation strategies that target environments contaminated by PAHs. Insights gained from the metabolic pathways of B. altitudinis DG4 could play a crucial role in developing effective bioremediation techniques tailored for naphthalene-contaminated sites, thereby enhancing the efficacy of environmental clean-up efforts [40]. Although our research contributes to the understanding of microbial degradation, caution is warranted when extending these findings to natural environments, as in vivo conditions can differ from laboratory settings [39].