ABSTRACT

A comprehensive microbial surveillance was conducted at NASA’s Mars 2020 spacecraft assembly facility (SAF), where whole-genome sequencing (WGS) of 110 bacterial strains was performed. One isolate, designated 179-BFC-A-HST, exhibited less than 80% average nucleotide identity (ANI) to known species, suggesting a novel organism. This strain demonstrated high-level resistance [minimum inhibitory concentration (MIC) >256 mg/L] to third-generation cephalosporins, including ceftazidime, cefpodoxime, combination ceftazidime/avibactam, and the fourth-generation cephalosporin cefepime. The results of a comparative genomic analysis revealed that 179-BFC-A-HST is most closely related to Virgibacillus halophilus 5B73CT, sharing an ANI of 78.7% and a digital DNA-DNA hybridization (dDDH) value of 23.5%, while their 16S rRNA gene sequences shared 97.7% nucleotide identity. Based on these results and the recent recognition that the genus Virgibacillus is polyphyletic, strain 179-BFC-A-HST is proposed as a novel species of a novel genus, Tigheibacillus jepli gen. nov., sp. nov (type strain 179-BFC-A-HST = DSM 115946T = NRRL B-65666T), and its closest neighbor, V. halophilus, is proposed to be reassigned to this genus as Tigheibacillus halophilus comb. nov. (type strain 5B73CT = DSM 21623T = JCM 21758T = KCTC 13935T). It was also necessary to reclassify its second closest neighbor Virgibacillus soli, as a member of a novel genus Paracerasibacillus, reflecting its phylogenetic position relative to the genus Cerasibacillus, for which we propose Paracerasibacillus soli comb. nov. (type strain CC-YMP-6T = DSM 22952T = CCM 7714T). Within Amphibacillaceae (n = 64), P. soli exhibited 11 antibiotic resistance genes (ARG), while T. jepli encoded for 3, lacking any known β-lactamases, suggesting resistance from variant penicillin-binding proteins, disrupting cephalosporin efficacy. P. soli was highly resistant to azithromycin (MIC >64 mg/L) yet susceptible to cephalosporins and penicillins.

IMPORTANCE

The significance of this research extends to understanding microbial survival and adaptation in oligotrophic environments, such as those found in SAF. Whole-genome sequencing of several strains isolated from Mars 2020 mission assembly cleanroom facilities, including the discovery of the novel species Tigheibacillus jepli, highlights the resilience and antimicrobial resistance (AMR) in clinically relevant antibiotic classes of microbes in nutrient-scarce settings. The study also redefines the taxonomic classifications within the Amphibacillaceae family, aligning genetic identities with phylogenetic data. Investigating ARG and virulence factors (VF) across these strains illuminates the microbial capability for resistance under resource-limited conditions while emphasizing the role of human-associated VF in microbial survival, informing sterilization practices and microbial management in similar oligotrophic settings beyond spacecraft assembly cleanrooms such as pharmaceutical and medical industry cleanrooms.

INTRODUCTION

Phylogenetic investigations into the Bacillaceae family have revealed a polyphyletic structure, signaling the necessity for taxonomic re-evaluation to accurately depict evolutionary lineages (1). This has been partially addressed by the recently proposed family Amphibacillaceae (2), encompassing halophilic genera like Virgibacillus (3), renowned for its robustness in saline environments (47). This family, which includes 16 genera, contains many bacterial taxa that have halophilic or halotolerant traits, with members commonly isolated from marine habitats, salt lakes, saline soils, and seafood, suggesting significance in food safety (8, 9).
Aerobic spore-forming bacteria, with their diverse ecological adaptations, are instrumental in biotechnological applications and have implications for human health due to their secondary metabolite production. These metabolites range from antimicrobial agents, such as bacitracin produced by Bacillus licheniformis (10) and polymyxin E produced by Paenibacillus polymyxa (11). Amphibacillaceae members like Oceanobacillus (6) not only demonstrate antimicrobial compound production but also present a fertile ground for discovering novel metabolites (12).
While typically non-pathogenic, some Bacillaceae members, such as Bacillus cereus and Bacillus anthracis, possess VF enabling opportunistic infections (13, 14). Conversely, Amphibacillaceae pathogens are rare; however, any bacterium can become opportunistic, especially in immunocompromised hosts. For example, a recent publication identified an Oceanobacillus spp. infection in an immunocompetent patient, with the pathogenesis mechanisms still largely unidentified, highlighting the knowledge gaps in the virulence of halophilic species (15).
Environmental Amphibacillaceae strains also harbor ARG, indicating that resistance can evolve even without direct clinically associated antibiotic exposure. The discovery of the β-lactamase OIH-1 in Oceanobacillus iheyensis that confers resistance to several penicillins but no other β-lactams exemplifies the independent evolution of AMR in remote environments like the oceanic depths (16, 17). Virgibacillus species have shown resistance to aminoglycosides and other antibiotics, with cephalosporins maintaining their clinical relevance against Gram-positive spore formers (18). Despite the availability of first-line treatments for Gram-positive spore-former-induced infections, such as vancomycin, erythromycin, and ciprofloxacin (19), the role of cephalosporins in treating such infections remains clinically relevant. These findings underscore the underappreciated role of natural environments in shaping AMR and necessitate vigilance against emerging ARG threats.
As part of a microbial surveillance study (6 months) associated with the NASA Mars 2020 mission, spacecraft assembly environments were assessed for microbial contamination (20). Bacterial strains (n = 110) were isolated from cleanroom samples that were subjected to a heat-shock (80°C; 15 min) procedure and cultured aerobically as per the NASA standard spore assay (21), and their genomes were sequenced (22).
The primary aim of this study was to characterize a bacterial isolate retrieved from the Mars 2020 mission assembly facility, identified as belonging to Amphibacillaceae through traditional microbiological methods and genomic analyses. Conserved marker genes were analyzed, and the phylogenetic affiliations of the isolate were delineated. A secondary objective was to determine the relative abundance of the novel isolate on the floors of the SAF cleanrooms (n = 236 samples) where the Mars 2020 mission rover had been assembled. Detailed genetic profiles of the species described herein were created, including ARG and VF profiles, and their distinct phenotypic traits were predicted. Their AMR phenotypes were assayed, and MIC were defined. Lastly, a predictive functional analysis was conducted, revealing genes encoding for putative bioactive compounds. This approach provides insights into the survival mechanisms of spore-forming bacteria in the harsh and nutrient-deficient environments of SAF cleanrooms, which could pose challenges to future NASA life detection missions.

RESULTS

Genome assembly

The genome size of the novel strain 179-BFC-A-HST is 3,877,640 bp, and G+C content is 41.64%. Strain 5B73CT has a 4,199,006-bp genome with 41.59% G+C content, and the genome of CC-YMP-6T strain is 3,876,281 bp with 36.13% G+C value. Other genome statistics are summarized in Table 1 and File S1.
TABLE 1
TABLE 1 Genomic information of sequenced isolatesb
SpeciesStrainNCBI accession no.Isolation locationNo. of contigsGenome size (bp)N50 (bp)Average coverageG + C content (%)Filtered short reads used for assemblyFiltered long reads used for assemblyCoding sequences16S rRNA genesa
Tigheibacillus jepli179-BFCA-HSJAROCA000000000.2Clean room floor, JPL spacecraft assembly facility93,939,6462,866,0021,12141.7919,068,144176,0004,40116
Tigheibacillus halophilus5B73C(T)JAWDIP000000000.1Field soil, Kakegawa, Shizuoka, Japan44,199,0062,848,45430941.59N/A261,0005,9637
Paraceracibacillus soliCC-YMP-6JAWDIQ000000000.1Yang-Ming Mountain, Taiwan73,876,2811,950,47371436.138,485,365152,0005,1127
Virgibacillus pantothenticus179-F W3.5 NHSJAYDXX000000000Clean room floor, JPL spacecraft assembly facility894,850,066115,79131637.1610,376,860N/A4,3421
Virgibacillus pantothenticus179-I 4A3 NHSJAYDXW000000000Clean room floor, JPL spacecraft assembly facility834,852,173115,79123937.167,840,196N/A4,3431
Virgibacillus pantothenticus179-I 6B2 NHSJAYDXV000000000Clean room floor, JPL spacecraft assembly facility844,851,827118,78420437.166,708,850N/A4,3451
Virgibacillus pantothenticus179_I 6B3 NHSJAYDXU000000000Clean room floor, JPL spacecraft assembly facility824,850,152129,45724937.168,174,846N/A4,3451
Virgibacillus pantothenticus179-I 6B4 NHSJAYDXT000000000Clean room floor, JPL spacecraft assembly facility904,849,929116,03319037.166,257,450N/A4,3421
a
Owing to the draft status of these genomes, only 16S rRNA gene sequences >1,400 base pairs were counted. For V. pantothenticus genomes, the actual count of 16S rRNA genes likely surpasses the single copy indicated in this table due to genome incompleteness.
b
Assembly statistics for T. jepli 179-BFC-A-HST and the genomes of type strains of the reclassified T. halophilus 5B73CT and P. soli CC-YMP-6T described in this study. The genomic information of V. pantothenticus (n = 5), also isolated from JPL-SAF, sequenced and characterized for the purposes of this study is also included.

Genotypic characterization of ARG and VF

In the Amphibacillaceae family (n = 64 species), 52 distinct ARG across 11 genera (Fig. 1A and B) were identified, with over 80% coverage and 70% identity to the NCBI’s ARG database. These ARG confer resistance to multiple antibiotic classes. Strain 5B73CT harbors 11 ARG, resistant to aminoglycosides, diaminopyrimidines, macrolide-lincosamide-streptogramin B (MLSB), tetracyclines, and glycopeptides. Fifteen species possess a single ARG, with strain 179-BFC-A-HST exhibiting three ARG related to vancomycin resistance.
Fig 1
Fig 1 ARG and phenotypic resistance to antibiotics: the mean resistome (± stdev) of the newly described T. jepli 179-BFC-A-HST and representative genomes of type species of the family Amphibacillaceae. (A) Number of ARG per genus analyzed (n = 11). (B) Breakdown of ARG per genome. Color gradient reflects percentage identity with respective gene in the NCBI database. Only hits with >70% identity and coverage are depicted. White color indicates absence of the gene. Tigheibacillus spp. is highlighted in red text, and Paracerasibacillus soli is highlighted in blue text. (C) MIC (mg/L) values for T. jepli 179-BFC-A-HST and P. soli CC-YMP-3T for 11 antibiotics interpreted with EUCAST clinical breakpoints under the “non-species” category: imipenem (IMI), ampicillin (AMP), ceftazidime-avibactam (CZA), ampicillin-sulbactam (SAM), cefotaxime (CT), tigecycline (TGC), ceftazidime (TZ), ciprofloxacin (CI), azithromycin (AZ), amikacin (AK), and aztreonam (AT).
Following a bacteremia case in South Korea caused by an Oceanobacillus species, a VF gene analysis across related genera revealed 17 human-associated VF genes (Fig. 2), aligning with the virulence factor database (VFDB) (23). Pseudogracilibacillus contained the most VF (n = 8), while Tigheibacillus had the fewest (n = 3). EF-Tu- and ATP-dependent Clp protease genes appeared in all 64 tested genomes. In contrast, hasC, bspD/E/F, and lspA were unique to single genomes: Cerasibacillus terrae (hasC), Oceanobacillus zhagkaii (bspE/F), and Ornithinibacillus scapharcae (lspA), respectively. These results underscore the varied resistance and virulence potentials within Amphibacillaceae, indicative of an environmental influence on emerging AMR and virulence.
Fig 2
Fig 2 Virulence factors. (A) The virulome of the newly described T. jepli 179-BFC-A-HST and genomes from the type species of associated genera (n = 64). The number of VF per genome is colored by associated genera. (B) Breakdown of VF per genome. Color gradient reflects percentage identity with respective gene in the VFDB. Only hits >70% identity are depicted. White color indicates absence of the gene. (C) An overview of mean number (± stdev) of identified VF per genus. Tigheibacillus spp. is highlighted in red text, and Paracerasibacillus soli is highlighted in blue text.

Phenotypic characterization of antimicrobial resistance properties

Cultivation-based analysis revealed that the novel strain 179-BFC-A-HST exhibited resistance to multiple antibiotics, including third-generation cephalosporins, ceftazidime, and cefpodoxime; fourth-generation cephalosporin cefepime; third-generation cephalosporin-β-lactamase inhibitor combination ceftazidime/avibactam and monobactam aztreonam (all with MIC >256 mg/L) but remained susceptible to penicillins [ampicillin and ampicillin/sulbactam (MIC <0.047 mg/L) (Fig. 1C). The strain CC-YMP-6T exhibited high-level resistance to macrolide azithromycin (MIC >64 mg/L) and monobactam aztreonam (MIC >256 mg/L) while being susceptible to all tested penicillins, cephalosporins, and cephamycins.

Phylogeny and genomic relatedness

The comparative genomic analysis indicates that the strain 179-BFC-A-HST shares considerable genomic similarity with strain 5B73CT, with an ANI of 78.7% (File S1) and a dDDH value of 23.5%. A heat map of ANI, including representative genomes of type species of the family Amphibacillaceae (n = 64), is shown in File S2. The gyrB sequence similarity between these strains is 79.13%, and the 16S rRNA gene sequence of strain 179-BFC-A-HST shares a 97.7% similarity with that of 5B73CT. The next closest genetic relatives based on 16S rRNA gene identity are Virgibacillus marismortui (95.77%) and Virgibacillus salaries (95.71%). Phylogenetic trees, including representative genomes of all validly published species of the family Amphibacillaceae, were reconstructed using 16S rRNA and gyrB gene sequences. These trees reveal a distinct clade, uniquely identifying the new strain 179-BFC-A-HST (File S3). Further average amino acid identity (AAI) comparisons indicate that strains 179-BFC-A-HST and 5B73CT share a 69.82% AAI, substantiating their classification within the same genus. In contrast, the AAI values between 179-BFC-A-HST and CC-YMP-6T, and between 5B73CT and CC-YMP-6T, are 58.53% and 55.88%, respectively, which support their categorization in different genera from CC-YMP-6T.
Recent analysis based on the Genome Taxonomy Database (GTDB) has reclassified many genera within Bacillaceae, including Virgibacillus into the new family Amphibacillaceae (2). Genomic analysis reveals Virgibacillus is polyphyletic, interspersed with genera like Oceanobacillus, Lentibacillus, Paucisalibacillus, and Cerasibacillus (Fig. 3). Consequently, strain 179-BFC-A-HST is proposed as a new genus and species, Tigheibacillus jepli gen. nov., sp. nov. V. halophilus 5B73CT is reclassified as Tigheibacillus halophilus comb. nov. and V. soli CC-YMP-6T as Paracerasibacillus soli comb. nov., positioned between Cerasibacillus and Tigheibacillus. The phylogenetic analysis indicates the genus Virgibacillus within Amphibacillaceae differentiates into clades E, F, and G, suggesting the need for further taxonomic studies to determine if these should be distinct genera.
Fig 3
Fig 3 Maximum-likehood phylogenomic tree based on 120 conserved marker genes showing the placement of 179-BFC-A-HST with the closely related members of Amphibacillaceae family. The tree was rooted on Thalassobacillus devorans CECT 7046 and Pontibacillus chungwhensis BH030062 as outgroup (not shown). Ultrafast bootstrap values >90% are indicated below the branches.

Biochemical and chemotaxonomic characteristics

Cell growth for T. jepli type strain 179-BFC-A-HST was measured at temperatures ranging from 4°C to 40°C, with an optimum at 35°C, and in 0.5%‒15% NaCl concentrations, optimally at 4%‒5%. Biochemical test results using the BioLog GNIII MicroPlate system and the Vitek 2 GP card are detailed in File S4. Fatty acid analysis revealed C 15:0 iso (37.8%) as the predominant component, followed by C 15:0 anteiso (24.2%), C 17:0 anteiso (13.3%), C 16:0 (9.0%), and C 16:0 iso (7.1%) (File S5). This contrasts with Virgibacillus pantothenticus CN8028T, which predominantly produces C 15:0 anteiso (47.4%). The polar lipid profile of strain 179-BFC-A-HST includes diphosphatidylglycerol (DPG), phosphatidylglycerol (PG), phosphatidylcholine (PC), phosphatidylmonomethylethanolamine (PME), phosphatidylethanolamine (PE), phosphatidylserine (PS), two unidentified aminolipids (AL1‒2), an unidentified phospholipid (PL1), and an unidentified aminophospholipid (APL1), as shown in File S6. Menaquinone MK-7 is the primary respiratory quinone, and A1γ with meso 2,6-diaminopimelic acid is the cell wall peptidoglycan’s diagnostic diamino acid.

Differential characteristics of Tigheibacillus genus

Table 2 highlights distinctive features of Tigheibacillus, differentiating it from other Amphibacillaceae genera strains. Our study compared phenotypes of type strains from five genera: Tigheibacillus, Virgibacillus, Oceanobacillus, Ornithinibacillus, and Lentibacillus, comprising 96 recognized species. T. jepli 179-BFC-A-HST grows between 4°C and 40°C (optimal at 35°C) and tolerates up to 15% NaCl, indicating moderate heat and high salt tolerance. This contrasts with Lentibacillus, which prefers cooler temperatures (30°C) and tolerates up to 23% salt. T. jepli 179-BFC-A-HST can produce acid from D-mannose and D-trehalose, a trait shared with Virgibacillus but not Ornithinibacillus or Lentibacillus. Its iso-C15:0/anteiso-C15:0 fatty acid ratio is 1.6, differing from Virgibacillus and Lentibacillus, implying variations in cell membrane structure. Tigheibacillus has a unique A1γ peptidoglycan type with meso 2,6-DAP, distinct from other genera’s peptidoglycan types. While sharing phosphatidylglycerol as a common lipid with other genera, its specific MK-7 quinone and DNA G+C content of 41.6% further set it apart biochemically.
TABLE 2
TABLE 2 Differential characteristics of T. jepli 179-BFC-A-HST and closely related type strains of the type species of the family Amphibacillaceaej
Characteristic1234e5h
TigheibacillusVirgibacillus (n = 34)Oceanobacillus (n = 29)Ornithinibacillus (n = 10)Lentibacillus (n = 19)
T. jepli 179-BFC-A HSTV. pantothenticusO. iheyensisO. bavariensisL. salicampi
OriginSAF/JPL,
USA
Soil,
UKa
Deep-sea sediment,
Japanb
Pasteurized milk, GermanySalt field,
Korea
Growth temperature range (°C)4‒40 (24)37c15‒42b15‒45 (25)15–40 (26)
NaCl tolerance (%)15 (4‒5)4c0‒21b10 (0.5‒4)2‒23
pH range for growth6–9 (7.5)7c6.5‒10b7‒106‒8
Acid production from:     
 D-mannose++c+b
 D-trehalose++cb+
 iso-C15:0/anteiso-C15:0 ratio1.60.3d<1e1.960.8
 Peptidoglycan typeA1γ with meso 2,6-DAPmeso DAP directfmeso DAP directeA4β (L-Orn←D-Asp)meso DAP
Polar lipids     
 IdentifiedDPG, PG, PME, PE, PC, PSDPGd, PGd, PEdDPGg, PGg, PEg, PCgDPG, PGDPG, PL
 UnidentifiedPL1, AL1‒2, APL1APLd, PLd, GLdPLg, Lg, ALgPL, L, APLUPL
Phosphatidylglycerol contentMajorMajorMajor to moderateMinorUnknown
Major quinoneMK-7MK-7dMK-7bMK-7MK-7
gDNA G+C content (%)41.6i37.1d35.8b36.444
a
Proom and Knight (27).
b
Lu et al. (6).
c
Heyndrickx et al. (26).
d
Wainø et al. (28).
e
Mayr et al. (9).
f
Claus and Berkeley (29).
g
Long et al. (data were not from type species) (30).
h
Yoon et al. (31).
i
Deduced from genome sequence data; SAF, Spacecraft Assembly Facility; JPL, Jet Propulsion Laboratory; DAP, Diaminopimelic Acid; +, Positive; −, Negative; ND, Not Determined.
j
Taxa: 1: T. jepli 179-BFC-A-HST; 2, Virgibacillus pantothenticus CN8028T (26); 3, Oceanobacillus iheyensis HTE831T (6); 4, Ornithinibacillus bavariensis WSBC 24001T (9); 5, Lentibacillus salicampi SF-20T (31).

Description of Tigheibacillus gen. nov.

Tigheibacillus (Tig.he.i.ba.cil’lus. L. masc. n. bacillus, a small rod; N.L. masc. n. Tigheibacillus a rod named to honor Scott Tighe, an American microbiologist, for his contribution in all areas of microbiology and extremophiles).
Cells are Gram-staining positive, strictly aerobic, non-spore forming, chemoheterotrophic, and mesophilic. They exhibit positive oxidase reaction and are catalase negative. The cells are flagellated, thick, slender rods with rounded ends, and endospores (Fig. 4). Their major fatty acids include 15:0 iso, 15:0 anteiso, and 17:0 anteiso. The predominant polar lipids are DPG, PG, PME, PE, and PS, with MK-7 being the major isoprenoid quinone. A1γ with meso 2,6-diaminopimelic acid characterizes the cell wall peptidoglycans. Based on the phylogenetic analyses the genus is placed within the family Amphibacillaceae, with Tigheibacillus jepli as the type species.
Fig 4
Fig 4 Scanning electron microscopy (SEM) images ([A] 10.0 μm and [B] 20 μm zoom) showcasing cellular structure and overall morphology, including presence of flagellum structure and spatial arrangement of T. jepli 179-BFC-A-HST.

Description of Tigheibacillus jepli sp. nov.

Tigheibacillus jepli (jep’li. N.L. gen. n. jepli, arbitrary name derived from the abbreviation JPL, meaning of or pertaining to the NASA’s Jet Propulsion Laboratory, where the type strain was isolated).
The following properties are observed in addition to those given for the genus description. Cells are 2.4–3.8 µm in length and 0.3–0.4 mm in width. On trypticase soy agar (TSA), after 1–2 days of incubation at 30°C, colonies are circular with regular margins, convex, 0.5–1.0 mm in diameter, and motile. Cell growth occurs at 4°C‒40°C (35°C optimum), 0.5%‒15% NaCl (4%‒5% optimum), and at pH 6.0–9.0 (optimum, pH 7.5). In addition to the major polar lipids listed in the genus description, PC, AL1‒2, PL1, and APL1 are produced in moderate-to-minor amounts.
The type strain of T. jepli 179-BFC-A-HST (DSM 115946T = NRRL B-65666T), isolated from spacecraft assembly cleanrooms.

Description of Tigheibacillus halophilus comb. nov.

Tigheibacillus halophilus (ha.lo.phi’lus. Gr. masc. n. hals, salt; Gr. masc. adj. philos, loving; N.L. masc. adj. halophilus, salt-loving).
Basonym: Virgibacillus halophilus An et al. (32).
The description is identical to that given by An et al. (32). The type strain, 5B73CT = DSM 21623T = JCM 21758T = KCTC 13935T, was isolated from field soil in Kakegawa, Shizuoka, Japan.

Description of Paracerasibacillus gen. nov.

Paracerasibacillus (Pa.ra.ce.ra.si.ba.cil’lus. Gr. prep. para, beside, near; N.L. masc. n. Cerasibacillus, a bacterial genus name; N.L. masc. n. Paracerasibacillus beside the genus Cerasibacillus, referring to the close but distinct relationship to the genus Cerasibacillus).
The reclassification of Virgibacillus soli to Paracerasibacillus soli is based on the phylogenetic analyses of 120 conserved protein genes. The description of the genus is identical to that given by Kämpfer et al. (33). The key properties of this species, now representing the genus Paracerasibacillus, include Gram-positive bacilli, predominant isoprenoid quinone of menaquinone MK-7, a polar lipid profile featuring diphosphatidylglycerol, phosphatidylglycerol, and phosphatidylethanolamine, and a polyamine pattern dominated by spermidine. Additionally, physiological and biochemical tests differentiate this strain from other Virgibacillus species. These distinguishing features support the reclassification of Virgibacillus soli as the type species of the new genus Paracerasibacillus, named Paracerasibacillus soli.

Description of Paracerasibacillus soli comb. nov.

Paracerasibacillus soli (so’li. L. gen. n. soli, of soil).
Basonym: Virgibacillus soli Kämpfer et al. (33).
The description is identical to that given by Kämpfer et al. (33). The type strain, CC-YMP-6T = DSM 22952T = CCM 7714T, was isolated from soil samples collected from Yang-Ming Mountain, Taiwan.

In silico metabolic predictions

In silico metabolic predictions for T. jepli, T. halophilus, P. soli, and 61 other genomes corresponding to type species representing genera within the Amphibacillaceae were conducted using the distilled and refined annotation of metabolism (DRAM) software (Fig. 5). The genome of T. jepli 179-BFC-A-HST encoded for parts 2 and 3 of the acetate to methane pathway, nitrate to nitrite and nitrite to nitrate metabolic pathways, arsenic and mercury reduction pathways, and as part 2 and part 1 of the butyrate and acetate pathways, respectively.
Fig 5
Fig 5 Phenotypic predictions based on genomic data using the DRAM tool. The chart showcases predicted metabolic capabilities and associated pathways, providing insight into the potential physiological roles and adaptations of T. jepli 179-BFC-A-HST and associated representative genomes of type species within the Amphibacillaceae family. (A) Electron transport chain complexes green color-coded gradient based on their completeness, ranging from 0 (absent) to 1 (fully complete). (B) Metabolic predictions are visualized with specific focus on carbohydrate-active enzyme, nitrogen metabolism, sulfur metabolism, various other reductases, photosynthesis, methanogenesis, methanotrophy, short-chain fatty acid conversions, and alcohol conversions. Metabolic features are color coded to indicate whether a function is present (purple) or absent (white).
The presence and characteristics of various biosynthetic gene clusters (BGCs) within the genome of T. jepli 179-BFC-A-HST were predicted using antiSMASH v.7.0.0 with a “strict” detection system. A loosely related (8% identity) BGC corresponding to a type III polyketide synthase cluster and specifically to legonindolizidine A6 was identified (Table S3). Furthermore, the production of a yet-to-be-characterized terpene cluster-associated metabolite was also predicted.

Metagenomics-based mapping to reads generated from JPL-SAF

The genome of the newly characterized T. jepli was tracked in metagenomic data sets gathered from the JPL-SAF during a period spanning 6 months in 2016, coinciding with the preparation phase for the Mars 2020 rover’s components. Analysis of 236 paired-end shotgun metagenomic data sets, with a threshold of greater than 1% genome coverage, revealed genomic evidence of T. jepli in 15 samples as depicted in Fig. 6. These samples were collected over six different dates ranging from 15 March to 28 June 2016. The findings suggest that T. jepli is present in low relative abundance and demonstrates limited spatial-temporal distribution within the SAF environment.
Fig 6
Fig 6 Mapping of the novel species to SAF metagenomic reads: box plots showing the breadth of coverage of the consensus genome were constructed from mapped reads aligned to T. jepli 179-BFC-A-HST (reported percentage breadth of coverage cutoff >1%). Reads mapping to the genome of T. jepli 179-BFC-A-HST were identified for 15 samples retrieved from 6 different collection dates.

DISCUSSION

The genomic landscape within the Bacillaceae family has undergone substantial revaluation, particularly with the reclassification of numerous genera following the rank normalization process by GTDB (2). Our comparative genomic analysis emphasizes this taxonomic fluidity, revealing that strain 179-BFC-A-HST is not closely aligned with the current genus Virgibacillus. Although ANI values below the 95%–96% threshold typically suggest a novel species (34), a universally accepted threshold for genus-level distinction is still elusive. To confirm our genomic findings, we conducted an in-depth phylogenetic analysis, using marker genes (16S rRNA gene and gyrB) and genome-based phylogeny involving 120 conserved protein markers. The average AAI value between P. soli and both T. jepli and T. halophilus genomes is substantially below the 63.43% suggested as a threshold for genus delineation (35). In addition, the dDDH analysis further supports the classification of isolate 179-BFC-A HST as a novel genus, with values below the 70%-mark indicative of a distinct species or potentially novel genera (36). While genomic data play a pivotal role in the delineation of species (ANI) and genera (AAI), the core attributes of an organism are frequently manifested through its phenotypic traits and the makeup of its cell wall. Consequently, we adopted an integrative methodology, combining genomic insights with phenotypic characteristics and chemotaxonomic data to achieve a comprehensive understanding. Based on both phenotypic and genotypic evidence and adhering to the rules of the International Code of Nomenclature of Prokaryotes (ICNP) (37), we propose the establishment of a new species within a new genus Tigheibacillus for the strain 179-BFC-A-HST.
Our phylogenetic analysis reinforces the polyphyletic nature of Virgibacillus (26), as members of this genus are separated by other distinct genera such as Oceanobacillus and Lentibacillus, suggesting a need for re-evaluation of its genus boundaries (6, 38). In contemporary research, genome-based phylogeny is widely recognized as the primary reference point for genus delineation despite the lack of uniform standards in this area. However, the absence of universally applicable thresholds presents significant challenges in making comparative analyses across taxa within the same hierarchical rank (39). This has been addressed by the introduction of the relative evolutionary divergence approach in GTDB (40), resulting in rank-normalized taxonomy for bacterial and archaeal taxa. The divergence of strains 179-BFC-A-HST and 5B73CT, leading to their reassignment to the newly proposed genus Tigheibacillus, reflects their genetic distinction within the Amphibacillaceae family. Moreover, based on the phylogenetic placement and relative evolutionary divergence, we suggest to place V. soli into its own genus Paracerasibacillus, to reflect its distinction from Virgibacillus (41). Furthermore, our study calls for more investigation into the clades E, F, and G within Virgibacillus, identified as potentially distinct genera.
Species within the Amphibacillaceae, including Virgibacillus and related genera, are generally not associated with human diseases. Yet, O. oncorhynchi subsp. incaldanensis has been reported to infect immunocompetent individuals (15) showing growth preference for alkaline environments like the male reproductive system (42). Furthermore, there are documented instances of Virgibacillus senegalensis and Lentibacillus sp. CBA3610 being isolated from the human gut (4, 24), indicating possible transient colonization with unknown systemic effects. Our VF analysis of the Amphibacillaceae family revealed 17 human-associated VF gene homologs. Notably, the hasC gene homolog in C. terrae’s genome is associated with biofilm formation and commonly found in pathogens like group B streptococci (43). The type IV secretion system effector BspE homolog in O. zhagkaii and the signal peptidase II encoding lspA homolog in O. scapharcae are associated with the human pathogens Brucella species and Rickettsia typhi, respectively (44, 45), highlighting a diverse range of VF. This range of human-associated VF in the genomes of environmental microorganisms signifies the importance of further understanding the ecological role of VF genes as well as the need for an in vitro evaluation of the pathogenic potential of such microorganisms.
T. jepli 179-BFC-A-HST exhibits an unusual antibiotic resistance phenotypic profile, showing high-level resistance to third- and fourth-generation cephalosporins (MIC >256 mg/L), even when these are combined with the β-lactamase inhibitor avibactam while remaining susceptible to ampicillin. This selective resistance pattern is atypical and not fully explained by the known type resistome as no known β-lactamases were identified. Albeit atypical, similar resistance patterns have been observed in Bacillus species within the Bacillaceae family. Adamski et al. (46) reported susceptibility to ampicillin and resistance to cefotaxime among a subset of isolates of B. cereus, B. pumilus, and B. licheniformis isolated from raw milk.
β-Lactam antibiotics function by binding to penicillin-binding proteins (PBPs), inhibiting cell wall synthesis (47). In Gram-positive bacteria, PBPs are more accessible due to the lack of an outer membrane (48), and resistance often involves PBPs with low β-lactam affinity (25, 49). The resistance phenotype in T. jepli 179-BFC-A-HST may be attributed to modifications in PBPs, diminishing their affinity for cephalosporins relative to penicillin antibiotics. The β-lactam-inducible PBPs of T. jepli 179-BFC-A-HST show substancial divergence (<63.7% amino acid identity) from those in closely related species, suggesting a distinct resistance profile. The isolation of cephalosporin-resistant bacteria like T. jepli 179-BFC-A-HST in cleanrooms signals a need for the pharmaceutical and medical industries to strengthen microbial reduction protocols and invest in new decontamination technologies. Such measures are vital to safeguard product integrity, comply with stricter regulatory standards, and mitigate health risks, ensuring the reliability of the global supply chain and maintaining industry reputation against the backdrop of rising antibiotic resistance (50).
Conversely, P. soli CC-YMP-6T presented high-level resistance to azithromycin (MIC > 64 mg/L), a clinically significant macrolide antibiotic for treating infections caused by Gram-positive bacteria (51). This resistance is likely due to the presence of plasmid-mediated erythromycin methyltransferase genes (ermG/Y/C), which confer resistance to MLSB antibiotics (52). Isolated from mountainous soil in Yang-Ming Mountain, Taiwan, P. soli CC-YMP-6T’s possession of erm might offer uncharacterized adaptive advantages or indicate anthropogenic environmental contamination. The detection of erm genes in soil samples, particularly in urban areas with high anthropogenic activity, varies widely (53), underscoring the need to understand the environmental and evolutionary factors influencing AMR in natural settings.
T. jepli 179-BFC-A-HST genomic analysis predicts a metabolic profile involved in methanogenesis, nitrogen cycling, and reductive reactions, notably in the latter stages of acetate-to-methane conversion. It also exhibits pathways for nitrogen utilization, reducing heavy metals like arsenate and mercury, while it lacks pathways for carbohydrate-active enzymes (CAZy) and short-chain fatty acid (SCFA)/alcohol transformations. The predicted heavy metal metabolic capacity of T. jepli 179-BFC-A-HST aligns with the members of the genus Bacillus known for their bioremediation potential in reducing environmental concentrations of heavy metals (54). The metabolic traits of T. jepli indicate its putative capacity to utilize energy and nutrients from limited sources, neutralize heavy metals, and participate in biogeochemical cycles, suggesting a resilient phenotype suitable for nutrient-limited environments such as controlled clean rooms. The resilience of T. jepli in oligotrophic environments, mirroring traits of related Amphibacillaceae genera like Oceanobacillus and Virgibacillus, highlights their efficient nutrient utilization and stress resistance (6, 33). Such traits are vital for survival in nutrient-poor habitats, suggesting a possible role in nutrient cycling and ecosystem stability (55). Studying these related genera enhances our understanding about the ecological role of T. jepli.
Metagenomic analysis in JPL-SAF revealed the infrequent occurrence of T. jepli 179-BFC-A-HST during the assembly phase of the Mars 2020 rover subsystems. Over 6 months in 2016, genomic signatures of T. jepli were identified in only 15 of 236 samples. This intermittent presence implies that T. jepli may have a niche-specific ecology or be a transient organism in the SAF environment. Generally, clean room microbiomes are dominated by gram-positive, spore-forming species, with Bacillus species comprising 10%–13% of isolates (56). Previous culture-based analyses in the JPL-SAF have identified 16 bacterial genera, frequently detecting Bacillus subtilis and V. pantothenticus (20). Bacillaceae and Amphibacillaceae species are notable for their persistence in such environments. Understanding the low prevalence of certain Amphibacillaceae species, including T. jepli, can inform contamination control strategies for space missions (57) and contribute to planetary protection efforts.
The presence of a high 16S rRNA gene copy number (16S GCN) in T. jepli is an intriguing aspect of the findings that merits further investigation. The 16S GCN in bacteria usually ranges from 1 to 15 and is often correlated with ecological strategies: oligotrophs, which thrive in nutrient-poor environments, usually have a lower 16S GCN, while copiotrophs, which flourish in nutrient-rich settings, have a higher 16S GCN, reflecting their adaptive strategies to their respective environment (58). The high full length 16S GCN in T. jepli (n = 16), found in the oligotrophic environment of the SAF, suggests it might confer an advantage for growth when conditions become favorable (59). The observed high 16S GCN copy number may also confer ecological versatility, allowing T. jepli to rapidly adjust to sporadic nutrient availability, which could be a result of evolutionary adaptations (60). While a high 16S GCN is advantageous for initial growth, it is suggested that the ability to endure extended periods of starvation involves more complex survival mechanisms that go beyond just the number of 16S GCN (61). Furthermore, the presence of 27 dormancy, sporulation, and spore germination-associated genes in the genome of T. jepli indicates their role in enabling survival in nutrient-poor environments, employing spore formation as an adaptive mechanism for survival.
The robustness of T. jepli in oligotrophic environments, coupled with its potential for heavy metal detoxification, mirrors similar capabilities observed in Bacillus species (54), indicating its potential suitability for bioremediation in nutrient-scarce, contaminated areas. Additionally, atypical resistance profile of T. jepli, alongside its occurrence in human-associated settings, necessitates further investigation into its pathogenicity. Studies focusing on its interactions with human cells and serum resistance are essential for evaluating potential health risks in medical and industrial contexts.

Conclusion

In conclusion, microbial surveillance conducted as part of the Mars 2020 mission revealed a novel bacterial strain, 179-BFC-A-HST, classified as T. jepli which represents the type strain of a novel species of a novel genus named herein. The high-level resistance of T. jepli to third- and fourth-generation cephalosporins, contrasted with its susceptibility to penicillins, along with its unique genetic profile, underscores the importance of characterizing resistance patterns, even in the absence of clinically associated environmental pressures. Furthermore, several human pathogenicity-associated VF were identified in the genome of T. jepli. The metabolic versatility of T. jepli, particularly its involvement in pathways such as methanogenesis, nitrogen cycling, reductive reactions, and heavy metal detoxification, highlights the resilience and adaptability of microorganisms in nutrient-limited environments. These findings are not only crucial for space exploration but also have significant implications for maintaining stringent microbial control standards in other sensitive industries, such as semiconductor manufacturing, pharmaceutical production, and medical device fabrication, where microbial contamination can have profound impacts.

MATERIALS AND METHODS

Sample collection, isolation, and preliminary identification of the strain 179-BFC-A-HST using 16S rRNA gene phylogeny were presented elsewhere (20). Additionally, type strains of CC-YMP-6T (formerly Virgibacillus soli DSM 22952T) and 5B73CT (formerly Virgibacillus halophilus DSM 21623T) were purchased from the Leibniz Institute DSMZ—the German Collection of Microorganisms and Cell Cultures GmbH. Their genomes were sequenced, and extensive AMR phenotypic analyses were performed.

DNA extraction: second- and third-generation sequencing

DNA from overnight-grown (30°C, TSA) cultures was extracted using a GeneJET Genomic DNA Purification Kit (Thermo Scientific, USA). Samples with >10 ng/µL concentration, A260/A280 ratios of 1.8–2.0, and A260/A230 ratios of 1.8–2.2 were selected for sequencing, and libraries prepared using the SQK-RBK114.24 Kit (Oxford Nanopore Technologies, UK) were sequenced on a MinION Mk1C system. Base calling was completedusing MinKnow/Guppy, adapter trimming with Porechop (v.0.2.4) (62), and read filtering by Filtlong v.0.2.1 with parameters --min_length 1000, --keep_percent 90. Hybrid genome assemblies for CC-YMP-6T (ASM3400633v1) and 179-BFC-A-HST (ASM3073200v2) used Unicycler (v.0.5.0) (63), while strain 5B73CT (ASM3400637v1) had a long-read assembly via Flye v.2.9.1 (64). Assembly quality was assessed using Quast (v.5.2.0) (65). For V. pantothenticus isolates 179-FW3.5N-HS, 179-I4A3N-HS, 179-I6B2N-HS, 179-I6B3N-HS, and 179-6B4N-HS, DNA was extracted, sequenced, and assembled as previously described (22).

Phenotypic and genotypic assessment of antimicrobial resistance properties and genotypic prediction of virulence potential

ARG and VF-encoding genes were classified using Abricate v. 1.0.0 against National Center for Biotechnology Information (NCBI) (66) and VFDB (23) databases, considering only hits with >80% nucleotide coverage and >70% identity. Data analysis and visualization were carried out using R programming. Antibiotic susceptibility testing was conducted on 11 antibiotics (imipenem, ampicillin, ceftazidime-avibactam, ampicillin-sulbactam, cefotaxime, tigecycline, ceftazidime, ciprofloxacin, azithromycin, amikacin, and aztreonam) using MIC Test Strips (Liofilchem, Italy) on Mueller-Hinton II agar (Fannin Ltd., Ireland). Susceptibility categories were assigned based on the EUCAST clinical breakpoints (v. 13.0, 2023) using PK-PD interpretive criteria (67).

Phylogeny and genome annotations

A two-step phylogenomic analysis was conducted on isolate 179-BFC-A-HST using 120 conserved single-copy marker proteins (40) as follows: initially using GTDB-Tk v2.3.0 (68) by running (i) the ani_rep command to compare genome sequence from strain 179-BFC-A-HST against all GTDB representative genomes (release 08-RS214) and (ii) classify_wf to find out the maximum-likelihood placement in the GTDB reference tree. Based on the output from the classify workflow (classify_wf), we selected a total of 61 genome sequences representing closely related type strains of all genera from the family Amphibacillaceae according to GTDB taxonomy (08-RS214), including genomes sequenced as part of this study from P. soli CC-YMP-6T and T. halophilus 5B73CT strains. Second, we inferred genome tree from the multiple sequence alignment of 61 selected genomes including the genome of strain 179-BFC-A-HST and two genomes representing an outgroup using FastTree v.2.1.11 (69) with the WAG model executed via GTDB-Tk v2.3.0. The inferred tree was further used as a starting tree in the IQ-TREE v.2.1.2 (70) executed via the following command: iqtree -s < alignment > m LG + C10+F + G -ft < starting_tree> -bb 1000. The tree was visualized, rooted, and pruned for illustration purposes in ARB v.6.0.6 and beautified in iTOL (71).
We assessed the genetic relatedness of 64 genomes in this study by creating an ANI heatmap with ANIclustermap (v1.2.0), which utilizes fastANI and seaborn algorithms. Genomic annotations were performed using PROKKA v1.14.5 with the --compliant extension (72). AAI values were determined using protein sequences as input based on the AAI calculator software developed by Kostas Lab. We extracted the gyrB and 16S rRNA gene sequences of T. jepli 179-BFC-A-HST and compared their percentage identities with their three closest sequences using BLASTn. Additionally, the dDDH value was estimated using formula 2 of the Genome-to-Genome Distance Calculator v.3.0 and the BLAST+ alignment tool (73). For the 61 species of Amphibacillaceae family, we extracted their 16S rRNA gene sequences from NCBI and aligned them with P. soli CC-YMP-6T, T. jepli 179-BFC-A-HST, and T. halophilus 5B73CT using MUSCLE v3.8.1551. A phylogenetic tree was constructed using IQTREE v.2.2.0.3 (74) with the GTR + G model and 1,000 ultrafast bootstrap replicates, while trees were visualized in iTOL v.6.7 (75).

Light microscopy and scanning electron microscopy

Spore formation was checked using malachite green and safranin followed by light microscopy using a DP25 camera and cellSens software (76). SEM analyses were carried out as previously published (77), and images were collected with an FEI Quanta 200F electron microscope at the California Institute of Technology’s Kavli Nanoscience Institute.

Biochemical and phenotypic characteristics

We used the Vitek 2 GP ID (bioMérieux) card for biochemical testing as per the manufacturer’s guidelines. Colonies from TSA were aseptically transferred into saline (0.45%–0.50% NaCl, pH 4.5–7.0) to create a 0.5-McFarland suspension (~1.5 × 108 cfu/mL), verified with VITEK 2 DensiCHEK Plus. This suspension, along with the Vitek 2 GP ID card, was incubated in a cassette at 37°C. Cassette loading, data entry, and retrieval of results (within 10 hours of inoculation) followed the VITEK instrument manual.
For phenotypic fingerprinting, a GNIII MicroPlate (BioLog) was used, preparing inoculum from TSA colonies in solution A (Cat # 72401; BioLog) to achieve a McFarland of 0.5. A volume of 100 µL of this bacterial suspension was loaded into each well of a GNIII MicroPlate and incubated at 37°C for 24 hours. OmniLog readings (A590–A750) were taken with a FLUO star Omega MicroPlate reader (BMG Labtech, Germany) within 24 hours of incubation.
Temperature tolerance and optimal growth temperatures were tested by incubating cells in tryptic soy broth at various temperatures (4°C, 10°C, 20°C, 27°C, 35°C, 40°C, 50°C, 70°C, 100°C). Salt tolerance and optimal concentration were assessed by growing cells in R2A broth with varying NaCl concentrations (0%–18% [wt/vol], in increments of 0.5% up to 6% and larger increments thereafter). For pH tolerance and optimal pH, cells were cultured across a pH range of 4.0–10.0 (in 1.0 pH unit intervals) using specific buffer systems: 0.1 M citric acid/0.1 M trisodium citrate for pH 4.0–5.0; 0.2 M Na2HPO4/0.2 M NaH2PO4 for pH 6.0–8.0; and 0.1 M NaHCO3/0.1 M Na2CO3 for pH 9.0–10.0. The pH was verified and adjusted post-sterilization of the media (78).

Chemotaxonomic analysis

Cells grown on TSA at 30°C for 48 hours to mid-exponential phase were used for fatty acid methyl esters (FAMEs) analysis. FAMEs were extracted from cell biomass via sequential saponification, methylation, and extraction (79) and then separated using a gas chromatograph (Agilent 7890A) with a flame ionization detector. Identification of FAMEs utilized the Microbial Identification System (MIDI) (80) and the Aerobe database (RTSBA6) in Sherlock version 6.0, adhering to the standard protocol (81).
Cells grown on TSA for 3 days at 30°C were used for extracting polar lipids, quinones, and peptidoglycans. Polar lipids and quinones were analyzed using 2D thin-layer chromatography (TLC) (82). TLC plates were sprayed with various reagents for lipid identification: 10% ethanolic phosphomolybdic acid for total lipids, 0.2% ninhydrin in butanol for aminolipids, Dittmer and Lester’s Zinzadze reagent for phospholipids, and alpha-naphthol for glycolipids. Peptidoglycans characterization followed established protocols by Staneck and Roberts (83).

Genotypic prediction of metabolic traits

The DRAM software v.1.4.6 (84) was used to analyze the genomic assemblies of all three strains tested during this study and 61 associated representative genomes of type species within the Amphibacillaceae family, for electron transport chain complexes and key metabolic functions, assigning a completeness scale from 0 (absent) to 1 (fully complete). This analysis included CAZy profiles, nitrogen and sulfur metabolism, reductases, photosynthesis pathways, methanogenesis, methanotrophy, and SCFA and alcohol conversions. Features were color coded for presence or absence in the genomic data. For strain 179-BFC-A-HST, antiSMASH v.7.0.1 was additionally used with a “strict” detection method (85) to predict BGCs associated with secondary metabolites, aligning them with the MIBiG v.3.1 database for functional annotations (85).

Mapping of the novel species to SAF metagenomic reads

Shotgun metagenomics sequence data sets (236 samples) from the NASA’s-JPL SAF were tested to evaluate the prevalence of strain 179-BFC-A-HST (20). Quality control and adapter trimming were conducted using Fastp v.0.23.4 with default settings. To map and quantify the metagenome reads relative to the novel species, MetaCompass v.2.0 was used with a maximum read length of 150, employing a reference-guided assembly approach. This included using Bowtie2 for read alignment and MEGAHIT for consensus genome construction. Samples without sequence data mapping the reference genome were excluded from further analysis.

Footnote

This article is a direct contribution from Nikos C. Kyrpides, a Fellow of the American Academy of Microbiology, who arranged for and secured reviews by Fumito Maruyama, Office of Academic Research and Industry-Government Collaboration, Academy of Hiroshima University, and Alexandre Rosado, KAUST, King Abdullah University of Science and Technology.

SUPPLEMENTAL MATERIAL

Supplemental material - mbio.00181-24-s0001.docx
Supplemental figures and tables.
ASM does not own the copyrights to Supplemental Material that may be linked to, or accessed through, an article. The authors have granted ASM a non-exclusive, world-wide license to publish the Supplemental Material files. Please contact the corresponding author directly for reuse.

REFERENCES

1.
Patel S, Gupta RS. 2020. A phylogenomic and comparative genomic framework for resolving the polyphyly of the genus Bacillus: proposal for six new genera of Bacillus species, Peribacillus gen. nov., Cytobacillus gen. nov., Mesobacillus gen. nov., Neobacillus gen. nov., Metabacillus gen. nov. and Alkalihalobacillus gen. nov. Int J Syst Evol Microbiol 70:406–438.
2.
Chuvochina M, Mussig AJ, Chaumeil P-A, Skarshewski A, Rinke C, Parks DH, Hugenholtz P. 2023. Proposal of names for 329 higher rank taxa defined in the genome taxonomy database under two prokaryotic codes. FEMS Microbiol Lett 370:fnad071.
3.
Heyrman J, Logan NA, Busse H-J, Balcaen A, Lebbe L, Rodriguez-Diaz M, Swings J, De Vos P. 2003. Virgibacillus carmonensis sp. nov., Virgibacillus necropolis sp. nov. and Virgibacillus picturae sp. nov., three novel species isolated from deteriorated mural paintings, transfer of the species of the genus Salibacillus to Virgibacillus, as Virgibacillus marismortui comb. nov. and Virgibacillus salexigens comb. nov., and emended description of the genus Virgibacillus. Int J Syst Evol Microbiol53:501–511.
4.
Ahn J, Hayes RB. 2021. Environmental influences on the human microbiome and implications for noncommunicable disease. Annu Rev Public Health 42:277–292.
5.
Qu J-H, Fu Y-H, Yue Y-F, Li H-F. 2019. Description of Ornithinibacillus gellani sp. nov., a halophilic bacterium isolated from lake sediment, and emended description of the genus Ornithinibacillus. Int J Syst Evol Microbiol 69:2632–2637.
6.
Lu J, Nogi Y, Takami H. 2001. Oceanobacillus iheyensis gen. nov., sp. nov., a deep-sea extremely halotolerant and alkaliphilic species isolated from a depth of 1050 m on the Iheya Ridge. FEMS Microbiol Lett 205:291–297.
7.
Konate S, Camara A, Lo CI, Tidjani Alou M, Hamidou Togo A, Niare S, Armstrong N, Djimdé A, Thera MA, Fenollar F, Raoult D, Million M. 2021. Virgibacillus doumboii sp. nov., a halophilic bacterium isolated from the stool of a healthy child in Mali. New Microbes New Infect 42:100890.
8.
Akuzawa S, Nagaoka J, Kanekatsu M, Kanesaki Y, Suzuki T. 2016. Draft genome sequence of Oceanobacillus picturae Heshi-B3, isolated from fermented rice bran in a traditional Japanese seafood dish. Genome Announc 4:e01621-15.
9.
Mayr R, Busse H-J, Worliczek HL, Ehling-Schulz M, Scherer S. 2006. Ornithinibacillus gen. nov., with the species Ornithinibacillus bavariensis sp. nov. and Ornithinibacillus californiensis sp. nov. Int J Syst Evol Microbiol 56:1383–1389.
10.
Haavik HI. 1975. Bacitracin production by the neotype; bacillus licheniformis ATCC 14580. Acta Pathol Microbiol Scand Suppl 83B:534–540.
11.
Yu Z, Sun Z, Yin J, Qiu J. 2018. Enhanced production of polymyxin E in Paenibacillus polymyxa by replacement of glucose by starch. Biomed Res Int 2018:1–10.
12.
Asharaf S, Chakraborty K. 2022. Pharmacological potential of seaweed-associated heterotrophic Firmicutes. Lett Appl Microbiol 75:1042–1054.
13.
Owusu-Kwarteng J, Parkouda C, Adewumi GA, Ouoba LII, Jespersen L. 2022. Technologically relevant Bacillus species and microbial safety of West African traditional alkaline fermented seed condiments. Crit Rev Food Sci Nutr 62:871–888.
14.
Enosi Tuipulotu D, Mathur A, Ngo C, Man SM. 2021. Bacillus cereus: epidemiology, virulence factors, and host–pathogen interactions. Trends Microbiol 29:458–471.
15.
Kim M, Choi JH, Shin H, Lee WK, Byun J-H. 2022. First case of bacteremia and epididymo-orchitis caused by Oceanobacillus oncorhynchi subspecies incaldanensis in an immunocompetent patient. Ann Lab Med 42:290–292.
16.
Toth M, Smith C, Frase H, Mobashery S, Vakulenko S. 2010. An antibiotic-resistance enzyme from a deep-sea bacterium. J Am Chem Soc 132:816–823.
17.
Matobole RM, van Zyl LJ, Parker-Nance S, Davies-Coleman MT, Trindade M. 2017. Antibacterial activities of bacteria isolated from the marine sponges Isodictya compressa and Higginsia bidentifera collected from Algoa Bay, South Africa. Mar Drugs 15:47.
18.
Muriuki SW, Neondo JO, Budambula NLM. 2020. Detection and profiling of antibiotic resistance among culturable bacterial isolates in vended food and soil samples. Int J Microbiol 2020:6572693.
19.
Aygun FD, Aygun F, Cam H. 2016. Successful treatment of Bacillus cereus bacteremia in a patient with propionic acidemia. Case Rep Pediatr 2016:6380929.
20.
Hendrickson R, Urbaniak C, Minich JJ, Aronson HS, Martino C, Stepanauskas R, Knight R, Venkateswaran K. 2021. Clean room microbiome complexity impacts planetary protection bioburden. Microbiome 9:238.
21.
NASA. 2010. Handbook for the microbiological examination of space hardware. NASA-HDBK-6022, Washington, D.C., USA.
22.
Wood JM, Peres JB, Bateh J, Singh NK, Aronson HS, Hendrickson R, Mason CE, Venkateswaran K. 2021. Draft genome sequences of heat shock-tolerant microbes isolated from a spacecraft assembly facility. Microbiol Resour Announc 10:e0065321.
23.
Liu B, Zheng D, Zhou S, Chen L, Yang J. 2022. VFDB 2022: a general classification scheme for bacterial virulence factors. Nucleic Acids Res 50:D912–D917.
24.
Seck E, Rathored J, Khelaifia S, Croce O, Robert C, Couderc C, Di Pinto F, Sokhna C, Raoult D, Lagier J-C. 2015. Virgibacillus senegalensis sp. nov., a new moderately halophilic bacterium isolated from human gut. New Microbes New Infect 8:116–126.
25.
Malouin F, Bryan LE. 1986. Modification of penicillin-binding proteins as mechanisms of beta-lactam resistance. Antimicrob Agents Chemother 30:1–5.
26.
Heyndrickx M, Lebbe L, Kersters K, De Vos P, Forsyth G, Logan NA. 1998. Virgibacillus: a new genus to accommodate Bacillus pantothenticus (Proom and Knight 1950). emended description of Virgibacillus pantothenticus. Int J Syst Bacteriol48:99–106.
27.
Proom H, Knight B. 1950. Bacillus pantothenticus (n.sp.). J Gen Microbiol 4:539–541.
28.
Wainø M, Tindall BJ, Schumann P, Ingvorsen K. 1999. Gracilibacillus gen. nov., with description of Gracilibacillus halotolerans gen. nov., sp. nov.; transfer of Bacillus dipsosauri to Gracilibacillus dipsosauri comb. nov., and Bacillus salexigens to the genus Salibacillus gen. nov., as Salibacillus salexigens comb. nov. Int J Syst Bacteriol 49 Pt 2:821–831.
29.
Claus D, Berkeley RC. 1986. Genus Bacillus Cohn 1872. In Bergey’s manual of systematic bacteriology. Williams & Wilkins, Baltimore.
30.
Long X, Ye R, Zhang S, Liu B, Zhang Y, Zeng Z, Tian Y. 2015. Oceanobacillus damuensis sp. nov. and Oceanobacillus rekensis sp. nov., isolated from saline alkali soil samples. Antonie Van Leeuwenhoek 108:731–739.
31.
Yoon J-H, Kang KH, Park Y-H. 2002. Lentibacillus salicampi gen. nov., sp. nov., a moderately halophilic bacterium isolated from a salt field in Korea. Int J Syst Evol Microbiol 52:2043–2048.
32.
An S-Y, Asahara M, Goto K, Kasai H, Yokota A. 2007. Virgibacillus halophilus sp. nov., spore-forming bacteria isolated from soil in Japan. Int J Syst Evol Microbiol 57:1607–1611.
33.
Kämpfer P, Arun AB, Busse H-J, Langer S, Young C-C, Chen W-M, Syed AA, Rekha PD. 2011. Virgibacillus soli sp. nov., isolated from mountain soil. Int J Syst Evol Microbiol 61:275–280.
34.
Goris J, Konstantinidis KT, Klappenbach JA, Coenye T, Vandamme P, Tiedje JM. 2007. DNA–DNA hybridization values and their relationship to whole-genome sequence similarities. Int J Syst Evol Microbiol 57:81–91.
35.
Park M-J, Kim YJ, Park M, Yu J, Namirimu T, Roh Y-R, Kwon KK. 2022. Establishment of genome based criteria for classification of the family Desulfovibrionaceae and proposal of two novel genera, Alkalidesulfovibrio gen. nov. and Salidesulfovibrio gen. nov. Front Microbiol 13.
36.
Stackebrandt E, Frederiksen W, Garrity GM, Grimont PAD, Kämpfer P, Maiden MCJ, Nesme X, Rosselló-Mora R, Swings J, Trüper HG, Vauterin L, Ward AC, Whitman WB. 2002. Report of the ad hoc committee for the re-evaluation of the species definition in bacteriology. Int J Syst Evol Microbiol 52:1043–1047.
37.
International code of nomenclature of prokaryotes. 2019. Int J Syst Evol Microbiol 69:S1–S111.
38.
Yoon J-H, Kang KH, Park Y-H. 2002. Lentibacillus salicampi gen. nov., sp. nov., a moderately halophilic bacterium isolated from a salt field in Korea. Int J Syst Evol Microbiol 52:2043–2048.
39.
Avise JC, Johns GC. 1999. Proposal for a standardized temporal scheme of biological classification for extant species. Proc Natl Acad Sci U S A 96:7358–7363.
40.
Parks DH, Chuvochina M, Waite DW, Rinke C, Skarshewski A, Chaumeil P-A, Hugenholtz P. 2018. A standardized bacterial taxonomy based on genome phylogeny substantially revises the tree of life. Nat Biotechnol 36:996–1004.
41.
Nakamura K, Haruta S, Ueno S, Ishii M, Yokota A, Igarashi Y. 2004. Cerasibacillus quisquiliarum gen. nov., sp. nov., isolated from a semi-continuous decomposing system of kitchen refuse. Int J Syst Evol Microbiol 54:1063–1069.
42.
Mishra AK, Kumar A, Swain DK, Yadav S, Nigam R. 2018. Insights into pH regulatory mechanisms in mediating spermatozoa functions. Vet World 11:852–858.
43.
Schindler Y, Rahav G, Nissan I, Treygerman O, Prajgrod G, Attia BZ, Raz R, Valenci GZ, Tekes-Manova D, Maor Y. 2023. Group B streptococcus virulence factors associated with different clinical syndromes: asymptomatic carriage in pregnant women and early-onset disease in the newborn. Front Microbiol 14:1093288.
44.
Rahman MS, Ceraul SM, Dreher-Lesnick SM, Beier MS, Azad AF. 2007. The lspA gene, encoding the type II signal peptidase of Rickettsia typhi: transcriptional and functional analysis. J Bacteriol 189:336–341.
45.
Byndloss MX, Tsolis RM. 2016. Brucella spp. virulence factors and immunity. Annu Rev Anim Biosci 4:111–127.
46.
Adamski P, Byczkowska-Rostkowska Z, Gajewska J, Zakrzewski AJ, Kłębukowska L. 2023. Prevalence and antibiotic resistance of Bacillus sp. isolated from raw milk. Microorganisms 11:1065.
47.
Rice LB. 2012. Mechanisms of resistance and clinical relevance of resistance to β-lactams, glycopeptides, and fluoroquinolones. Mayo Clin Proc 87:198–208.
48.
Sauvage E, Terrak M. 2016. Glycosyltransferases and transpeptidases/penicillin-binding proteins: valuable targets for new antibacterials. Antibiotics (Basel) 5:12.
49.
Fisher JF, Mobashery S. 2016. β-lactam resistance mechanisms: gram-positive bacteria and Mycobacterium tuberculosis. Cold Spring Harb Perspect Med 6:a025221.
50.
Mahnert A, Moissl-Eichinger C, Zojer M, Bogumil D, Mizrahi I, Rattei T, Martinez JL, Berg G. 2019. Man-made microbial resistances in built environments. Nat Commun 10:968.
51.
Dinos GP. 2017. The macrolide antibiotic renaissance. Br J Pharmacol 174:2967–2983.
52.
Beirne C, McCann E, McDowell A, Miliotis G. 2022. Genetic determinants of antimicrobial resistance in three multi-drug resistant strains of Cutibacterium acnes isolated from patients with acne: a predictive in silico study. Access Microbiol 4.
53.
Yi X, Wang M, Zhou Z. 2019. The potential impact of naturally produced antibiotics, environmental factors, and anthropogenic pressure on the occurrence of erm genes in urban soils. Environmental Pollution 245:282–289.
54.
Wróbel M, Śliwakowski W, Kowalczyk P, Kramkowski K, Dobrzyński J. 2023. Bioremediation of heavy metals by the genus Bacillus. Int J Environ Res Public Health20:4964.
55.
Haruta S, Kanno N. 2015. Survivability of microbes in natural environments and their ecological impacts. Microbes Environ 30:123–125.
56.
Sandle T. 2011. A review of cleanroom microflora: types, trends, and patterns. PDA J Pharm Sci Technol 65:392–403.
57.
McDonagh F, Cormican M, Morris D, Burke L, Singh NK, Venkateswaran K, Miliotis G. 2023. Medical astro-microbiology: current role and future challenges. J Indian Inst Sci 103:771–796.
58.
Roller BRK, Stoddard SF, Schmidt TM. 2016. Exploiting rRNA operon copy number to investigate bacterial reproductive strategies. Nat Microbiol 1:16160.
59.
Klappenbach JA, Dunbar JM, Schmidt TM. 2000. rRNA operon copy number reflects ecological strategies of bacteria. Appl Environ Microbiol 66:1328–1333.
60.
Puigbò P, Wolf YI, Koonin EV. 2010. The tree and net components of prokaryote evolution. Genome Biol Evol 2:745–756.
61.
Chen Y, Neilson JW, Kushwaha P, Maier RM, Barberán A. 2021. Life-history strategies of soil microbial communities in an arid ecosystem. ISME J 15:649–657.
62.
Wick RR, Judd LM, Gorrie CL, Holt KE. 2017. Completing bacterial genome assemblies with multiplex MinION sequencing. Microb Genom 3:e000132.
63.
Wick RR, Judd LM, Gorrie CL, Holt KE. 2017. Unicycler: resolving bacterial genome assemblies from short and long sequencing reads. PLoS Comput Biol 13:e1005595.
64.
Kolmogorov M, Yuan J, Lin Y, Pevzner PA. 2019. Assembly of long, error-prone reads using repeat graphs. Nat Biotechnol 37:540–546.
65.
Gurevich A, Saveliev V, Vyahhi N, Tesler G. 2013. QUAST: quality assessment tool for genome assemblies. Bioinformatics 29:1072–1075.
66.
Feldgarden M, Brover V, Haft DH, Prasad AB, Slotta DJ, Tolstoy I, Tyson GH, Zhao S, Hsu C-H, McDermott PF, Tadesse DA, Morales C, Simmons M, Tillman G, Wasilenko J, Folster JP, Klimke W. 2019. Validating the AMRFinder tool and resistance gene database by using antimicrobial resistance genotype-phenotype correlations in a collection of isolates. Antimicrob Agents Chemother 63:e00483-19.
67.
EUCAST. 2020. The European committee on antimicrobial susceptibility testing. breakpoint tables for interpretation of MICs and zone diameters. Version10.0. European Committee on Antimicrobial Susceptibility Testing
68.
Chaumeil P-A, Mussig AJ, Hugenholtz P, Parks DH. 2022. GTDB-Tk v2: memory friendly classification with the genome taxonomy database. Bioinformatics 38:5315–5316.
69.
Price MN, Dehal PS, Arkin AP. 2010. FastTree 2 – approximately maximum-likelihood trees for large alignments. PLoS One 5:e9490.
70.
Nguyen L-T, Schmidt HA, von Haeseler A, Minh BQ. 2015. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol 32:268–274.
71.
Letunic I, Bork P. 2019. Interactive tree of life (iTOL) V4: recent updates and new developments. Nucleic Acids Res 47:W256–W259.
72.
Seemann T. 2014. Prokka: rapid prokaryotic genome annotation. Bioinformatics 30:2068–2069.
73.
Auch AF, Klenk H-P, Göker M. 2010. Standard operating procedure for calculating genome-to-genome distances based on high-scoring segment pairs. Stand Genomic Sci 2:142–148.
74.
Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, von Haeseler A, Lanfear R. 2020. IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol Biol Evol 37:1530–1534.
75.
Letunic I, Bork P. 2021. Interactive tree of life (iTOL) V5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res 49:W293–W296.
76.
Schaeffer AB, Fulton MD. 1933. A simplified method of staining endospores. Science 77:194.
77.
Simpson AC, Sengupta P, Zhang F, Hameed A, Parker CW, Singh NK, Miliotis G, Rekha PD, Raman K, Mason CE, Venkateswaran K. 2023. Phylogenomics, phenotypic, and functional traits of five novel (earth-derived) bacterial species isolated from the international space station and their prevalence in metagenomes. Sci Rep 13:19207.
78.
Xu P, Li W-J, Tang S-K, Zhang Y-Q, Chen G-Z, Chen H-H, Xu L-H, Jiang C-L. 2005. Naxibacter alkalitolerans gen. nov., sp. nov., a novel member of the family ‘Oxalobacteraceae’ isolated from China. Int J Syst Evol Microbiol 55:1149–1153.
79.
Miller LT. 1982. Single derivatization method for routine analysis of bacterial whole-cell fatty acid methyl esters, including hydroxy acids. J Clin Microbiol 16:584–586.
80.
Sasser M. 1990. Bacterial identification by gas chromatographic analysis of fatty acid methyl esters (GC-FAME). MIDI-Inc., technical note #101
81.
Paisley R. 1996. MIS whole cell fatty acid analysis by gas chromatography training manual, Newark
82.
Minnikin DE, Minnikin SM, Parlett JH, Goodfellow M, Magnusson M. 1984. Mycolic acid patterns of some species of Mycobacterium. Arch Microbiol 139:225–231.
83.
Staneck JL, Roberts GD. 1974. Simplified approach to identification of aerobic actinomycetes by thin-layer chromatography. Appl Microbiol 28:226–231.
84.
Shaffer M, Borton MA, McGivern BB, Zayed AA, La Rosa SL, Solden LM, Liu P, Narrowe AB, Rodríguez-Ramos J, Bolduc B, Gazitúa MC, Daly RA, Smith GJ, Vik DR, Pope PB, Sullivan MB, Roux S, Wrighton KC. 2020. DRAM for distilling microbial metabolism to automate the curation of microbiome function. Nucleic Acids Res 48:8883–8900.
85.
Terlouw BR, Blin K, Navarro-Muñoz JC, Avalon NE, Chevrette MG, Egbert S, Lee S, Meijer D, Recchia MJJ, Reitz ZL, et al. 2023. MIBiG 3.0: a community-driven effort to annotate experimentally validated biosynthetic gene clusters. Nucleic Acids Res 51:D603–D610.

Information & Contributors

Information

Published In

cover image mBio
mBio
Volume 15Number 410 April 2024
eLocator: e00181-24
Editor: Igor B. Jouline, The Ohio State University, Columbus, Ohio, USA
PubMed: 38477597

History

Received: 24 January 2024
Accepted: 29 January 2024
Published online: 13 March 2024

Keywords

  1. cephalosporin
  2. phylogenetic analysis
  3. Virgibacillus taxonomy

Contributors

Authors

Antimicrobial Resistance and Microbial Ecology Group, School of Medicine, University of Galway, Galway, Ireland
Centre for One Health, Ryan Institute, University of Galway, Galway, Ireland
Author Contributions: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing – original draft, and Writing – review and editing.
Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai, Tamil Nadu, India
Center for Integrative Biology and Systems mEdicine (IBSE), Indian Institute of Technology Madras, Chennai, Tamil Nadu, India
Robert Bosch Centre for Data Science and Artificial Intelligence (RBCDSAI), Indian Institute of Technology Madras, Chennai, Tamil Nadu, India
Author Contributions: Data curation, Formal analysis, Investigation, Software, Validation, Visualization, Writing – original draft, and Writing – review and editing.
Division of Microbiology and Biotechnology, Yenepoya Research Centre, Yenepoya (Deemed to be University), Mangalore, Karnataka, India
Author Contributions: Formal analysis, Resources, Validation, Visualization, Writing – original draft, and Writing – review and editing.
The University of Queensland, School of Chemistry and Molecular Biosciences, Australian Centre for Ecogenomics, Brisbane, Australia
Author Contributions: Data curation, Formal analysis, Methodology, Software, Visualization, Writing – original draft, and Writing – review and editing.
Antimicrobial Resistance and Microbial Ecology Group, School of Medicine, University of Galway, Galway, Ireland
Author Contributions: Data curation, Formal analysis, Investigation, Methodology, Software, Visualization, Writing – original draft, and Writing – review and editing.
Anna C. Simpson
Biotechnology and Planetary Protection Group, Jet Propulsion Laboratory, California Institute of Technology, Pasadena, California, USA
Author Contributions: Formal analysis, Investigation, Methodology, Supervision, Visualization, Writing – original draft, and Writing – review and editing.
Ceth W. Parker
Biotechnology and Planetary Protection Group, Jet Propulsion Laboratory, California Institute of Technology, Pasadena, California, USA
Author Contributions: Validation, Visualization, Writing – original draft, and Writing – review and editing.
Nitin K. Singh
Biotechnology and Planetary Protection Group, Jet Propulsion Laboratory, California Institute of Technology, Pasadena, California, USA
Author Contributions: Conceptualization, Data curation, Formal analysis, Methodology, Supervision, Writing – original draft, and Writing – review and editing.
Division of Microbiology and Biotechnology, Yenepoya Research Centre, Yenepoya (Deemed to be University), Mangalore, Karnataka, India
Author Contributions: Methodology, Project administration, Resources, Visualization, Writing – original draft, and Writing – review and editing.
Dearbháile Morris
Antimicrobial Resistance and Microbial Ecology Group, School of Medicine, University of Galway, Galway, Ireland
Centre for One Health, Ryan Institute, University of Galway, Galway, Ireland
Author Contributions: Investigation, Supervision, Writing – original draft, and Writing – review and editing.
Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai, Tamil Nadu, India
Center for Integrative Biology and Systems mEdicine (IBSE), Indian Institute of Technology Madras, Chennai, Tamil Nadu, India
Robert Bosch Centre for Data Science and Artificial Intelligence (RBCDSAI), Indian Institute of Technology Madras, Chennai, Tamil Nadu, India
Author Contributions: Investigation, Validation, Visualization, Writing – original draft, and Writing – review and editing.
US Department of Energy Joint Genome Institute, Berkeley, California, USA
Environmental Genomics and Systems Biology Division, Lawrence Berkeley National Laboratory, Berkeley, California, USA
Author Contributions: Investigation, Methodology, Writing – original draft, and Writing – review and editing.
The University of Queensland, School of Chemistry and Molecular Biosciences, Australian Centre for Ecogenomics, Brisbane, Australia
Author Contributions: Conceptualization, Investigation, Methodology, Supervision, Writing – original draft, and Writing – review and editing.
Biotechnology and Planetary Protection Group, Jet Propulsion Laboratory, California Institute of Technology, Pasadena, California, USA
Author Contributions: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, and Writing – review and editing.

Editor

Igor B. Jouline
Editor
The Ohio State University, Columbus, Ohio, USA

Notes

The authors declare no conflict of interest.

Metrics & Citations

Metrics

Note:

  • For recently published articles, the TOTAL download count will appear as zero until a new month starts.
  • There is a 3- to 4-day delay in article usage, so article usage will not appear immediately after publication.
  • Citation counts come from the Crossref Cited by service.

Citations

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. For an editable text file, please select Medlars format which will download as a .txt file. Simply select your manager software from the list below and click Download.

View Options

Figures

Tables

Media

Share

Share

Share the article link

Share with email

Email a colleague

Share on social media

American Society for Microbiology ("ASM") is committed to maintaining your confidence and trust with respect to the information we collect from you on websites owned and operated by ASM ("ASM Web Sites") and other sources. This Privacy Policy sets forth the information we collect about you, how we use this information and the choices you have about how we use such information.
FIND OUT MORE about the privacy policy