The family Flavobacteriaceae is one of the major branches of the phylum Bacteroidota Krieg et al. 2021 (Oren and Garrity, 2021). Over the past two decades, numerous novel genera and species have been reported as new member of the family Flavobacteriaceae. However, phylogenetic positioning based on 16S rRNA gene sequence of strains within the family remained somewhat ambiguous, especially in the branch including the genera Algibacter, Hyunsoonleella, Lacinutrix, etc. The genus Seonamhaeicola, adjacent to the genus Hyunsoonleella, also belong to the family Flavobacteriaceae and was first proposed by Park et al. (2014) and its description was subsequently amended by Zhou et al. (2016). Most members of the genus Seonamhaeicola are characterized as Gram-stain-negative, rod-shaped, facultatively anaerobic, heterotrophic, non-motile, non-spore-forming, yellowish- to orange- coloured bacteria, with menaquinone-6 (MK-6) as the major respiratory quinone, and iso-C15:0 as the major fatty acid. Members of the genus have been isolated from diverse habitats such as coral (Yoon et al., 2019), marine algae (Zhou et al., 2016; Fang et al., 2017), and sediments (Park et al., 2014; Cao et al., 2020; Zhang et al., 2020). A novel marine bacterial strain affiliated to the genus Seonamhaeicola has been isolated from a debris of Eisenia bicyclis and its taxonomic properties are reported here.
During the exploration of marine bacteria, strain MEBiC01930T (=DOKDO 020T) was isolated from a debris of marine algae (Eisenia bicyclis) sampled at Dokdo Island, Republic of Korea (37° 14' 34'' N, 131° 52' 00'' E) in Oct. 2004. For strains isolation, a bit of algal debris (approx. 0.1 g) was grinded and serially diluted with sterilized seawater, plated on marine agar 2216 medium (MA; BD) and cultivated at 25℃ for 3 days, then, distinctive colonies were selected for pure culture. After primary isolation and purification, strain MEBiC01930T was cultivated at 25°C on the same medium for biochemical and physiological characterization. This strain was also stored at −80°C in marine broth 2216 (MB; BD) with 20% (v/v) glycerol. For phenotypic comparisons, Seonamhaeicola aphaedonensis KCTC 32578T, S. acroporae KCTC 62713T, S. algicola KCTC 42396T, S. maritimus KCTC 72528T, S. sediminis KCTC 72085T, and S. marinus NBRC 112333T were purchased from KCTC (Korean Collection for Type Cultures) or NBRC (Biological Resource Center, NITE) and grown on MA at 25℃. Unless stated otherwise, strain MEBiC01930T was cultivated at 25℃ for 2 days.
Extraction of genomic DNA was carried out using a commercial DNA extraction kit (MoBio), and the 16S rRNA gene was amplified by using a bacterial primer set 27F and 1492R (Giovannoni, 1991), the detailed procedure and conditions were described in Lee et al. (2013). The 16S rRNA gene sequencing was performed using an ABI 3730xl automatic sequencer. The obtained 16S rRNA gene sequence was compared by BLAST pair-wise alignment with sequences in the EzBioCloud database (Kim et al., 2012). Phylogenetic analysis based on the almost-complete 16S rRNA gene sequence (1,407 bp) of strain MEBiC01930T was conducted with members of the genus Seonamhaeicola and closely related genera by using MEGA11 (Tamura et al., 2021). The phylogenetic tree with 1,000 replicated bootstrap analysis and the complete deletion option was reconstructed by using the neighbor-joining method [NJ; (Saitou and Nei, 1987)] with the Kimura 2 parameter model (Kimura, 1980) as well as the maximum-likelihood [ML; (Felsenstein,1981)] and maximum-parsimony [MP; (Fitch, 1971)] methods. The sequences of Flavobacterium aquatile DSM 1132T (AM230485) and Kordia algicida OT-1T (AY195836) were used as outgroups (Fig. 1).
The genome sequencing of strain MEBiC01930T was performed using a PacBio RS II system. Sequence assembly was performed by using HGAP version 3.0. Obtained genome sequence was annotated using the NCBI Prokaryotic Genome Annotation Pipeline (Tatusova et al., 2016). A phylogenomic analysis was conducted with the obtained genome sequence of strain MEBiC01930T with the closely related strains in the family Flavobacteriaceae (Table 1) using PhyloPhlAn (Segata et al., 2013). Average nucleotide identity (ANI) and Average amino acids identity (AAI) values were computed using the OrthoANIu algorithm (Yoon et al., 2017) and EzAAI (Kim et al., 2021), respectively. The genus delineation cut-off of AAI values was calculated according to Park et al. (2022).
Unless otherwise stated, the physiological and morphological characterization of the strain was conducted according to methods described by Yang et al. (2006). Transmission electron micrographs were taken using a LIBRA120 (Carl-Zeiss) electron microscope with fixed cells that were negatively stained with 2% phosphotungstic acid at pH 7.0. Additionally, scanning electron micrographs (SEM) were taken using a VEGAII LSU (Tescan) after dehydrated with series of ethanol. The range and optimal growth temperatures were determined in MB at 12 different temperatures (10, 16, 20, 23, 26, 29, 32, 36, 39, 42, 46 and 50°C) in a temperature gradient incubator (TVS126MA; Adaventec) for up to 3 days. The tolerance range for NaCl was tested in modified ZoBell 2216 broth (ZoBell, 1941) prepared with distilled water and supplemented with the NaCl (Sigma; 0, 0.5, 1, 2, 3, 3.5, 4, 6, 10, 15, and 20%, w/v) for 1 week. The tolerance range for pH was determined (pH 4, 5, 6, 6.5, 7, 7.5, 8, 9, and 10) in MB with the pH adjusted by 1 N HCl or 1 N NaOH solutions and using 10 mM MES (pH 4–6), 10 mM HEPES (pH 6–8) or 10 mM AMPSO (pH 8–10) as biological buffers for 1 week. The tolerance ranges for NaCl and pH were also determined in a temperature gradient incubator and the OD600 was monitor at 10 min intervals. The bacterial suspension used to inoculate into API 20E, 20NE, API ZYM kit (bioMérieux) and a Microlog GN2 system (Biolog) was prepared in a 2% sea salt (Sigma) solution. API and Microlog panels were recorded after incubation at 25°C for 2 days. To confirm anaerobic growth, cells were inoculated with 0.3% sodium nitrate and 0.5% yeast extract in the serum vial capped with aluminum seal, purged with deoxygenated nitrogen gas until color of indicator (final 0.2 mg/L of resazurin) was disappeared, and cultivated at 25°C for 5 days and observing growth.
The cellular fatty acid profiles of strain MEBiC01930T and reference strains were determined using the MIDI/Hewlett Packard Microbial Identification System (MIS) (Sasser, 1990) with Sherlock version 6.2 and the TSBA6 database according to the manufacturer’s instructions at KCCM (Korea Culture Center of Microorganisms). Briefly, strains were streaked on marine agar to produce four segments, and the third segment was sampled when growth was observed at the fourth segment. Then, the samples were saponified, methylated, extracted, purified, and analyzed by gas chromatograph. Polar lipids were extracted using a chloroform/methanol system and separated by two-dimensional TLC using silica gel 60 F254 aluminum-backed thin-layer plates (Merck) (Minnikin et al., 1984). The detailed procedure is described in Yang et al. (2013). After 2-dimensional development, each component was visualized using the following reagents; all lipids -10 % (w/v) molybdatophosphoric acid, free amino groups -0.2% ninhydrin solution (Consden and Gordon, 1948), phosphorus - Zinzadze reagent (Dittmer and Lester, 1964), sugar groups - α- naphthol reagent (Jacin and Mishkin, 1965), and phosphatidylcholine (PC) - Dragendorff’s reagent. The major respiratory quinone in strain MEBiC01930T was determined by HPLC analysis according to the method described by Collins (1985).
The GenBank/DDBJ/EMBL accession number for the 16S rRNA gene sequence of strain MEBiC01930T is DQ191181 and that of the genome sequence is CP104270.
After purification of the strain using standard agar plating method, a yellow-coloured, circular, convex, and opaque colony that was butyrous and entire edges was selected and designated as strain MEBiC01930T (=DOKDO020T). Strain MEBiC01930T was found to be closely related to Seonamhaeicola maritimus 1505T and S. algicola Gy8T with 97.84% and 97.71% 16S rRNA gene sequence similarity, respectively, and all other type species showing lower than 97% similarity. The 16S rRNA gene sequence similarity values were lower than the cut-off value (98.65%) of the novel species suggested by Kim et al. (2014). The NJ tree revealed that strain MEBiC01930T formed a coherent clade with members of the genus Seonamhaeicola but phylogenetic position pf species within the genus was not stable (Fig. 1). This result implied that strain MEBiC01930T could be a separate species in the genus Seonamhaeicola but genome based phylogenetic analysis was required for determining exact position.
The genome of strain MEBiC01930T was assembled into 1 contig with a size of 3.65 Mbp (3652651 bp). Calculated coverage was 4437×, and considering the self-mapping coverage and 16S rRNA gene sequences, genome was not contaminated by materials from other organisms. Briefly, the genome contains 3,102 protein-coding genes (CDSs) in addition to the 2 set of rRNAs, 49 tRNAs, 4 ncRNAs and 9 pseudogenes. The DNA G + C content is estimated to be 32.4 mol% (Table 1). Details on the genome contents can be found in Oh et al. (2024) and they reported that the strain possesses a set of alginate lyases such as guluronate-specific alginate lyases, mannuronate-specific alginate lyases, and an oligo-alginate lyase, but was not agarose, carragenase and cellulae. The 2 set of 16S rRNA genes encoded in the genome were identical but shows 1 bp difference against that by Sanger method. One or two copies of rRNA gene cluster are common in the closely related strains (Table 1). The copy number of the rRNA operon has been thought to be related to growth rate but no direct relationship was shown in nature. However, it can be related to the availability of limited nutrients (Valdivia-Anistro et al., 2016).
The phylogenomic tree inferred by PhyloPhlan showed that strain MEBiC01930T was well-clustered with the members of the genus Seonamhaeicola and this result demonstrated that strain MEBiC01930T is a member of the genus Seonamhaeicola (Fig. 2). The ANI used only for species delineation but AAI values widely used not only species but genus delineation (Park et al., 2022). The ANI values calculated were lower than 86.61% and demonstrated that strain MEBiC01930T is a novel species (Table 2). In case of the AAI values, higher than 76.44% within the genus Seonamhaeicola. Yang et al. (2024) reported that AAI value for genus delineation at adjacent groups of the family Flavobacteriaceae ranged between 76.0–78.7% and the above value (76.44%) is in this range.
Cell images obtain by TEM and SEM were rod shaped (0.3–0.5 × 1.6–2.1 mm) without flagella (Fig. 3), and this result is well matched with previously reported species of the genus Seonamhaeicola (Zhang et al., 2020). Staining the cell wall of strain was resulted as Gram-negative. Growth observed at temperature between 20.0 and 32.2°C (optimum 26–29°C), pH between 6.5 and 8.5 (optimum pH 7.5) and in the presence of 0–4% (w/v, optimum 2%) NaCl. No growth was observed at 18°C and 35°C, with 5% NaCl, and at pH 6.0 and 9.0 (Table 3). Anaerobic growth was observed by reducing nitrate to nitrite. Enzyme activities and utilization of carbon and energy sources, as measured using commercial kits, are summarized in Table 3 and in the description of the species and it can be briefly summarized as follows; strain MEBiC01930T was able to utilize a variety of substrates such as carbohydrates, short chain organic acids, and amino acids. Compared with the previously reported species, basic physiological properties were similar. However, cell size of strain MEBiC01930T was smaller than others and growth ranges for temperature, NaCl, and pH revealed to be lower than previously reported species (Table 3).
The dominant fatty acids of strain MEBiC01930T were determined to be iso-C15:1 G (17.2%), iso-C15:0 (14.8%), iso-C15:0 3-OH (9.8%), and iso-C17:0 3-OH (12.2%) and summed feature 3 (C16:1 ω6c and/or C16:1 ω7c; 4.0%) (Table 4) and were not quite different from other species in the genus Seonamhaeicola (Table 4). The predominant polar lipids of strain MEBiC01930T were phosphatidylethanolamine (PE), two unidentified aminolipids (ALs), and additional three unidentified lipids (Ls) (Fig. 4). PE as the major polar lipid is consistent with previously reported results for the genus Seonamhaeicola. The major respiratory quinone of strain MEBiC01930T was determined to be MK-6 by HPLC analysis according to Collins (1985). The DNA G + C content deduced from genome sequence was 32.4 mol%.
Strain MEBiC01930T could be classified into the genus Seonamhaeicola on the basis of the phylogenetic and phylogenomic trees (Figs. 1 and 2), fatty acids profile (Table 4), types of respiratory quinone, DNA G + C ratio, reduction of nitrate, enzyme activities determined by API ZYM kit (Table 3) etc. However, the isolate could be distinguished from other type strains by the existence of some short-chain fatty acids (Table 4), range and optimum of pH and NaCl for growth (Table 3), and the assimilation of some carbon sources such as gluconate and malate (Table 3). Additionally, ANI and AAI values of strain MEBiC01930T against members of the genus Seonamhaeicola were lower than the species delineation range (Kim et al., 2014) (Table 2). On the basis of this polyphasic taxonomical evidence, it was concluded that strain MEBiC01930T is a novel species in the genus Seonamhaeicola and proposed it as a name Seonamhaeicola citrea sp. nov.
Seonamhaeicola citrea (ci. tre′a. N.L. fem. adj. citrea lemon-yellow coloured)
Gram-stain-negative, oxidase and catalase positive, facultatively anaerobic, non-motile, and rod-shaped cells are 0.3–0.5 × 1.6–2.1 μm in size. The ranges and optimal conditions of temperature, pH and NaCl for growth are 20.0–32°C (optimum, 26–29°C), pH 6.5–8.5 (optimum, pH 7.5) and 0–4% (w/v; optimum, 2%), respectively. Nitrate is reduced to nitrite. Urease, β-galactosidase and β-glucosidase activities are present. Arginine dihydrolase, lysine decarboxylase and ornithine decarboxylase activities are not detected in API20E and 20NE kits. When assayed with the API ZYM system, alkaline phosphatase, esterase (C4), esterase lipase (C8), leucine arylamidase, acid phosphatase, naphtol-AS-BI phosphohydrolase, β-glucosidase, α-glucosidase and N-acetyl-β-glucosaminidase activities are present. Lipase (C14), valine arylamidase, trypsin, α-chymotrypsin, cystine arylamidase, α-galactosidase, β-galactosidase, β-glucuronidase, α-mannosidase and α-fucosidase activities are absent. In the API 50 CHB/E test, acid is produced from esculin, sucrose and 5-keto-gluconate, but not from others. Major fatty acids (>5.0%) are C12:0, iso-C15:1 G, iso-C15:0, iso-C15:0 3-OH and iso-C17:0 3-OH when grown at 25°C. The major respiratory quinone is MK-6. The polar lipids present are phosphatidylethanolamine, two unidentified aminolipids and three unidentified lipids. The DNA G + C content is 32.4 mol%.
The type strain MEBiC01930T (=KCCM 42330T =JCM 14619T) was isolated from a debris of marine alga Eisenia bicyclis collected from Dokdo Island, Republic of Korea. The GenBank/DDBJ/EMBL accession number for the 16S rRNA gene is DQ191181 and that of the genome is CP104270.
그람 음성의 통성 혐기성의 이동성이 없는 간균 (1.6 × 2.1 μm – 0.3 × 0.5 μm) MEBiC01930T (= DOKDO020T)은 2004년 10월 우리나라 독도(131°52'00''E, 37°14'34''N)에서 채취된 대황 잔해에서 분리되었다. 16S rRNA 유전자 서열 분석 결과 MEBiC01930T 균주는 Seonamhaeicola 속 균주들과 높은 유사성(최고 97.84%)을 보였다. MEBiC01930T 균주는 섭씨 20–32도 (최적 26–29도), pH 6.5–8.5 (최적 pH 7.5) 및 0–4% (최적 2%) NaCl 농도 조건에서 성장이 관찰되었다. 주요 세포 지방산은 iso-C15:1 G (17.2%), iso-C15:0 (14.8%), iso-C17:0 3-OH (12.2%), iso-C15:0 3-OH (9.8%), C12:0 (5.2%)였다. 균주의 유전체는 3.65 Mbp 크기에 3,170개의 유전자를 지니는 것으로 확인되었으며 DNA G + C 함량은 33.6 mol%이다. 유전체 정보 분석 결과 균주는 alginate lyase 종류를 여러 개 가지고 있음을 확인하였다. 주요 호흡기 퀴논은 MK-6이다. 이와 같은 다상분류 결과에 기초하여 균주 MEBiC01930T는 Seonamhaeicola 속의 새로운 종으로 분류되어야 함을 제안하며 이를 Seonamhaeicola citrea sp. nov.로 제안하였다. 표준균주는 MEBiC01930T (= KCCM 42330T = JCM 14619T)이다.
This work was supported by the KIOST in-house program (PEA0212) and the management of Marine Fishery Bioresources Center (2024) funded by the National Marine Biodiversity Institute of Korea (MABIK).
The authors declare that they have no conflict of interests.