Sampling and isolation of strains
The strain CQH2019T was isolated from pit mud obtained from a Nongxiang-flavor Baijiu company in Sichuan Province, China. Pit mud samples were collected from the bottom of cellar and immediately transferred to an anaerobic bag with an oxygen indicator (Hopebio HBYY007, Qingdao, China) and transported to the laboratory for bacterial separation.
For enrichment cultivation, 5 g homogenized pit mud was inoculated into 50 mL reinforced clostridial medium (RCM; ATCC Medium 1524) containing (per liter): 10.0 g beef extract, 3.0 g yeast extract, 3.0 g sodium acetate, 10.0 g tryptone, 3.0 g glucose, 1.0 g soluble starch, and 0.15 g CaCl₂·2H₂O. The medium pH was adjusted to 6.8 ± 0.1 using 1 mol/L NaOH. Enrichment cultures were maintained in COY vinyl anaerobic chambers (Coy Laboratory Products, Model 1000) with 85% N₂, 10% H₂, 5% CO₂ atmosphere at 37.0°C ± 0.5°C for 15 days. Subsequent isolation employed on solid RCM agar (add 1.5% w/v agar) with three successive streak-plate purification under anaerobic conditions. The selected pure culture (strain CQH2019T) with 25% (final v/v) glycerol was cryopreserved at -80°C.
Morphology
Use YCFA solid medium for the pure culture of strain CQH2019T and perform morphological observation of the strain. YCFA medium contained: tryptone, 10 g; yeast extract, 2.5 g; sodium acetate, 3.5 g; NaCl, 4.5g; CaCl2·6H2O, 0.09g; KH2PO4, 0.45g; K2HPO4, 0.45g; MgSO4·7H2O, 0.09g; L-Cysteine Hydrochloride, 0.8g; Sodium propionate, 0.865g; Sodium butyrate, 0.11g; Agar powder, 17.5g. The pH was adjusted to pH 7.0 with 10M NaOH. In an anaerobic chamber at 37°C for 48 h for observation. Cell morphology of the isolated strain was examined using a light microscope (Leica DM500, Leica Microsystems, Germany) and transmission electron microscope(Hitachi HT7800, Hitachi High-Tech, Japan).
Physiology and chemotaxonomy
To determine the salinity, pH, and temperature tolerance of strain CQH2019T, cultures were grown anaerobically in bottles containing 30 mL of Reinforced Clostridial Medium (RCM) broth. Salinity tolerance was assessed at NaCl concentrations ranging from 0 to 7% (w/v). The pH growth range was evaluated across pH 3.0–9.0 (adjusted at 0.5-unit intervals), and the temperature tolerance was tested from 4 to 50°C. All experiments were performed in triplicate biological replicates. Growth was monitored by measuring the optical density at 600 nm (OD₆₀₀) to determine the optimal conditions. To determine catalase activity, bacterial cells grown on RCM medium were exposed to a 3% hydrogen peroxide solution. Amylase activity was assessed by applying approximately 0.1 ml of Lugol's iodine solution around colonies grown on RCM agar plates containing 2% starch. Oxidase activity and nitrate reduction were identified using Oxyswab (bioMérieux, Marcy-l'Étoile, France) according to the manufacturer’s instructions.(Lan, Chen, Lin, Ye, Yan, Huang, Liu, Yang, 2016). Gram staining was conducted using a commercial kit (Hangzhou Tianhe Microbiology Reagent Co., Ltd.) following manufacturer instructions .
In RCM broth, after 7 days of anaerobic incubation at 37°C, the utilization of carbohydrates were determined according to the API 50CH (bioMérieux, Marcy-l'Étoile, France) instructions. and used HPLC (Shimadzu Corporation, Kyoto, Japan) to measure its fermentation end products(Duarte-Villanueva et al., 2025).The fatty acids were extracted following the instructions supplied with the MIDI / Hewlett Packard Microbial Identification System (MIDI Inc., Newark, DE, USA), and the subsequent separation and identification by Agilent 6850 gas chromatography (Agilent Technologies, Santa Clara, CA, USA) were performed (Athalye et al., 1985). The polar lipid composition was analyzed using two-dimensional thin-layer chromatography (TLC Merck, Germany)(Minnikin et al., 1979). The first-dimensional separation was performed with a chloroform-methanol-water solvent system (65:25:4, v/v), followed by second-dimensional development with chloroform-methanol-acetic acid-water (85:12:15:4, v/v). Total lipids were visualized by spraying with 10% (v/v) molybdatophosphoric acid solution. Isoprenoid quinones were extracted and analyzed according to the protocol established by Minnikin et al. (1984) (Minnikin et al., 1984). For quantification of short-chain fatty acids in fermentation products, including lactic acid, acetic acid, and butyric acid, identical analytical methodology was employed as described by Reyes et al. (2025) (Reyes et al., 2025).
Phylogeny analysis of 16S RNA gene
Genomic DNA was extracted from pure culture of strain CQH2019T using a commercial DNA extraction kit (Beijing SBS Genetech Co., Ltd). The nearly complete 16S rRNA gene sequence (1,458 bp) was amplified using universal primers 27F(5′-AGAGTTGATCCTGGCTCAG-3′) and 1492R (5′-GGTTACCTTGTTACGACTT-3′) (Ferdous, 2024). To ensure high accuracy of 16S rRNA gene, bidirectional Sanger sequencing was performed, and the final sequence was deposited in GenBank databases (accession number: OP925886). The 16S rRNA gene sequence was blasted using NCBI RefSeq database (https://www.ncbi.nlm.nih.gov/) and the EzBioCloud 16S database (https://www.ezbiocloud.net/), respectively. Related similar sequences were collected and phylogenetic trees were built by the neighbor-joining (NJ), minimum-evolution and maximum-likelihood (ML) algorithms using MEGA 7(Kimura, 1980; Kumar et al., 2016).
Genome Features
The whole-genome sequencing was conducted by Shanghai Meiji biomedical technology Co., Ltd (Shanghai, PR China) using Hiseq4000(Illumina)and Sequel II༈PacBio༉system according to the manufacturer’s protocols. The complete genome was assembled using SOAPdenovo2༈http://soap.genomics.org.cn/༉(Luo et al., 2012). Chromosomal completeness was confirmed by BUSCO (https://busco.ezlab.org/) with score of 99.2%. The genome sequences were annotated through the NCBI Prokaryotic Genome Annotation Pipeline (NCBI-PGAP) and deposited in the DDBJ, ENA, and GenBank databases under the INSDC framework (https://www.insdc.org/). The whole-genome sequencing project has been assigned with the BioProject accession number PRJNA908277. Transfer RNA (tRNA) genes were identified using tRNAscan-SE v1.21 and ARAGORN v1.2.41(Chan et al., 2019; Laslett et al., 2004), while ribosomal RNA (rRNA) genes were predicted with RNAmmer v1.2(Lagesen et al., 2007). The complete genome sequence was deposited in the NCBI GenBank database.
The digital DNA-DNA hybridization (dDDH) values were determined using Formula 2 of the Genome-to-Genome Calculator (http://ggdc.dsmz.de/ggdc.php) (Meier-Kolthoff et al., 2013). The ANI values between the genomes of strain CQH2019T and the close strains were calculated using the ANI calculator tool from EzBioCloud (Yoon et al., 2017) and amino acid identity (AAI) values were calculated by the AAI calculator (http://enve-omics.ce.gatech.edu/aai/) (Rodriguezr et al., 2016).
The taxonomic position of the strain was further confirmed by whole-genome phylogeny. A phylogenomic tree was constructed using the Up-to-date Bacterial Core Gene (UBCG) method, which leverages a set of 92 core genes (Na et al., 2018).