Phylogenetic and phylogenomic analysis
The query results against the GenBank database revealed that the 16S rRNA sequence of strain NTOU-MSR1T exhibits the highest similarity of 98.82% to O. marisflavi DSM 106032ᵀ (MG456888). Following closely is the similarity of 98.54% to O. doudoroffii DSM 7028ᵀ (AB021371.1). The phylogenetic tree based on 16S rRNA gene sequences showed that strain NTOU-MSR1T belongs to the Oceanimoans clade, with O. marisflavi DSM 106032ᵀ and O. doudoroffii DSM 7028ᵀ as its closest relative (Fig. 1). The phylogenomic tree was consistent with the phylogenetic trees based on 16S rRNA, clustering strain NTOU-MSR1T with O. marisflavi DSM 106032ᵀ and O. doudoroffii DSM 7028ᵀ (Fig. 2). The clustering results of the phylogenomic tree, constructed using whole genome sequences, closely resemble the findings of the phylogenetic analysis conducted using 16S rRNA gene sequences. The pairwise comparisons of digital DDH values between Oceanimonas pelagia NTOU-MSR1 and all Oceanimonas species are 32.80% (O. doudoroffii DSM 7028ᵀ), 25.90% (O. baumannii DSM 15594ᵀ), 26.60% (O. smirnovii ATCC BAA-899ᵀ), 35.40% (O. marisflavi DSM 106032ᵀ), 34.90% (O. sp. GK1), 32.50% (O. sp. MB9), and 33.10% (O. sp. NS1), individually, and these values are much lower than the recommended threshold of species (> 70%) of digital DDH [19]. Oceanimonas pelagia NTOU-MSR1 exhibits the highest ANI values of 87.88% with O. marisflavi DSM 106032ᵀ, then 87.80% (O. sp. GK1), 86.78 (O. sp. NS1), 86.66% (O. doudoroffii DSM 7028ᵀ), 86.50 (O. sp. MB9), 82.87% (O. smirnovii ATCC BAA-899ᵀ), 82.37% (O. baumannii DSM 15594ᵀ), respectively, which is below the 95% threshold established considering a species boundary. Based on the overall genome relatedness index (OGRI), the results indicate that this strain NTOU-MSR1T can be distinguished from its closest relatives, establishing it as a novel species within the Oceanimonas genus.
Previous research has demonstrated that the plastisphere-associated Oceanimonas sp. isolate can degrade low-density polyethylene (Joshi et al., 2022). In 2022, a phylogenetic analysis of three Oceanimonas sp. isolates (NIOT.PDB3- NIOT.PDB5, accession number: MZ708724-MZ708726) from the plastisphere was conducted. These isolates showed the closest genetic similarity to O. marisflavi DSM 106032ᵀ (strain 102-Na3ᵀ) with similarities of 98.84 to 98.98% (Joshi et al., 2022). However, with the addition of strain NTOU-MSR1T into the database, it became apparent that this novel species exhibited even higher similarity percentages, ranging from 99.23 to 99.87%. Consequently, it has been determined that all these isolates belong to this new species.
Genome analysis
The general features of strain NTOU-MSR1T are summarized in Table S1, and the genome is visualized in Supplementary Figure S1. The complete genome of the strain NTOU-MSR1 is comprised of a single circular chromosome with a length of 3,592,406-bp long. The GC content is 61.0%. No plasmid, 1 type I-F CRISPR–Cas system and 5 prophage regions are detected. Thirteen genetic islands, which include 295 genes, are identified. The genome is estimated to be 100% complete with contamination of 0.09% based on the CheckM analysis. Totally 3,321 protein-coding sequences (CDS) are predicted, and 94 tRNA, 25 rRNA (5S, 16S, 23S), and 4 ncRNA are detected (Table S1).
The genome of strain NTOU-MSR1 T contains several genes associated with the catechol branch of β- ketoadipate pathway involved in the degradation of PAHs (locus tag PU634_15375 to _15380, PU634_15400 to _15410, PU634_04665). These genes include essential enzymes such as catechol 1,2-dioxygenase (catA), muconate/chloromuconate family cycloisomerase, muconolactone delta-isomerase (catC), alpha/beta hydrolase, 3-oxoadipate CoA-transferase, 3-oxoadipyl-CoA thiolase (pcaF), and CoA-transferase. Additionally, the presence of genes encoding benzoate 1,2-dioxygenase (PU634_15415 to _15425) in the genome of strain NTOU-MSR1T indicates its capacity for aerobic benzoate degradation via the β-ketoadipate pathway. These findings highlight the metabolic capabilities of strain NTOU-MSR1T and suggest its potential for environmental remediation and applications in addressing environmental contamination caused by PAHs and benzoate-containing substances.
This genome also contains a number of genes that are related to resistance to heavy metals, encoding copper resistance system multicopper oxidase (PU634_03685, PU634_06750), copper resistance protein B (PU634_03690, PU634_06755), cupredoxin family protein (PU634_03695, PU634_06760), heavy metal sensor histidine kinase (PU634_06770), and heavy metal translocating P-type ATPase (PU634_02520) involved in copper homeostasis mechanisms (Giachino & Waldron, 2020). Other heavy metals resistant genes encoding chromate efflux transporter (chrA, PU634_02540), two nickel/cobalt efflux protein RcnA (PU634_03740 and PU634_16635), Ni(II)/Co(II)-binding transcriptional repressor RcnR (rcnR; PU634_03735, PU634_16640), magnesium/cobalt transporter CorA (corA, PU634_09050), and CNNM family magnesium/cobalt transport protein CorC (corC, PU634_11385) are also being detected. Besides, several genes related to arsenic resistance or arsenate detoxification are detected (PU634_07750 to PU634_07760, PU634_07770), which encode organoarsenical efflux MFS transporter ArsJ (arsJ), ArsJ-associated glyceraldehyde-3-phosphate dehydrogenase, ACR3 family arsenite efflux transporter (arsB), and metalloregulator ArsR/SmtB family transcription factor. The presence of these genes indicates the potential of strain NTOU-MSR1 to survive in environments contaminated by heavy metals and toxic compounds.
Furthermore, the genes related to the biosynthesis pathway of PHB are detected (PU634_06235, PU634_06240, PU634_06250), including the genes encodings acetyl-CoA C-acetyltransferase, acetoacetyl-CoA reductase (phbB) and class I poly(R)-hydroxyalkanoic acid synthase (phaC). The genome of strain NTOU-MSR1T also contains genes encoding a TIGR01841 family phasin protein (phaP, PU634_16465) and a hypothetical phasin family protein (PU634_06245), as depicted in Figure S2. Phasins, also known as granule-associated proteins, are surface proteins that stabilize the PHA storage granules in bacteria (Mezzina & Pettinari, 2016). A gene encoding polyhydroxyalkanoate synthesis repressor PhaR (phaR, PU634_16475), which represses and autoregulates the expression of genes encoding the phasin protein (PhaP) and PhaR itself, was also found. Interestingly, a peptidylprolyl isomerase-encoding gene (PU634_16470) is identified through NCBI PGAP; however, its relation to the PHB cluster, phasin genes, and regulatory genes remains unclear. The presence of these genes suggests that strain NTOU-MSR1 may be capable of synthesizing PHA.
The prediction of BGCs within the genome sequence is analyzed using antiSMASH 7.0, revealing the presence of 3 BGCs, which are responsible for the biosynthesis of ectoine, ribosomally synthesized and post-translationally modified peptides (RiPP)-like product, and non-ribosomal peptide synthetases (NRPS)-independent siderophore. The ectoine biosynthetic gene cluster (PU634_06015 to _06025) is responsible for the production of ectoine, which is a low-molecular-mass compatible solute. Ectoine acts as an osmoprotectant to counteract the increase in environmental salinity (Bursy et al., 2007). The genes detected in this cluster include diaminobutyrate acetyltransferase (ectA), diaminobutyrate-2-oxoglutarate transaminase (ectB), ectoine synthase (ectC), and a specialized aspartate kinase (PU634_06030) involved in ectoine biosynthesis. The RiPP-like cluster is detected based on the presence of the gene encoding 30S ribosomal protein S12 methylthiotransferase accessory factor YcaO (ycaO, PU634_02785). Additionally, the presence of a non-ribosomal peptide synthetase (NRPS)-independent, IucA/IucC-like (NI-) siderophore biosynthetic gene cluster was detected by antiSMASH 7.0, which shows high similarity (> 85%) to the vibrioferrin biosynthetic gene cluster from Azotobacter vinelandii CA (BGC0002527). It encodes two IucA/IucC family proteins (PU634_12740, PU634_12750), one MFS transporter (PU634_12745), one type III PLP-dependent enzyme (PU634_12755), and an additional hypothetical protein (PU634_12735). The presence of this cluster suggests the potential for siderophore production in order to acquire iron from the nutrient-poor marine environment. It is worth noting that this specific siderophore has not been previously reported in the Oceanimonas genus. However, further investigation is required to confirm the presence of this siderophore and identify its specific product.
Morphological, physiological, biochemical characteristics and chemotaxonomy
Strain NTOU-MSR1T was Gram-stain-negative, aerobic, motile, and rod-shaped. The size of the cell was approximately 0.5–0.6 µm in width and 1.8-2.0 µm in length, as illustrated in Figure S3. A single flagellum was observed on one end of the bacterium, indicating a monotrichous arrangement. Colonies were convex, smooth, round with a light pink color, and approximately 8.0 mm wide after three days of incubation at 30°C. Growth between 10–45°C was observed on TSB agar (Optimum at 30°C). The pH range for growth was 7.0–10.0 (optimal growth at pH 7.0–8.0). The strain grew in TSB broth containing 0%-12% (w/v). Catalase and oxidase activity were positive. The reduction of nitrates to nitrites was positive. Assimilations of malic acid and trisodium citrate were observed. The strain NTOU-MSR1T exhibited the ability to utilize D-aspartic acid, L-alanine, L-arginine, L-aspartic acid, L-glutamic acid, L-histidine, L-serine, methyl pyruvate, L-lactic acid, citric acid, L-malic acid, bromo-succinic acid, ɑ-hydroxy-butyric acid, β-hydroxy-D, L-butyric acid, ɑ -keto-butyric acid, propionic acid, acetic acid, and formic acid as single carbon. The differential morphological and physio-biochemical characteristics of strain NTOU-MSR1T and the experimental control strains were summarized in Table 1. Table S2 presented detailed information regarding the results of the BIOLOG GEN III Microplates analysis conducted on strain NTOU-MSR1T and other strains. The fatty acid identification result showed that the cell membrane of the strain NTOU-MSR1T primarily consisted of C16:1ω7c at 34.53% ± 1.22%, C18:1ω7c at 27.78% ± 1.38%, and C16:0 at 22.38% ± 0.77% (Table S3). The comparison of the cellular fatty acid profiles between strain NTOU-MSR1T and related species indicates that they share a similar composition of major fatty acids (Lee et al., 2018). Through in-silico analysis, key genes encoding the enzymes involved in the biosynthetic pathway were detected within the genome of the strain NTOU-MSR1T (Table S4). Previous reports have indicated that ubiquinol-8 (Q-8) was the predominant quinone in O. marisflavi DSM 106032T (strain 102-Na3T) (Lee et al., 2018). Based on the result of in silico analysis, it suggests the isoprenoid quinone type of strain NTOU-MSR1T might be ubiquinone.
Table 1
Difference of morphological, physio-biochemical characteristics of strain NTOU-MSR1T and closely related Oceanimonas species. Strains: 1, stain NTOU-MSR1T; 2, Oceanimonas doudoroffii DSM 7028T, and; 3, Oceanimonas baumannii DSM 15594T; 4, Oceanimonas marisflavi DSM 106032T. +, Positive reaction; -, Negative reaction; W, Weakly positive reaction.
Characteristics | 1 | 2 | 3 | 4 |
Colony Color | Pink | White | White | Pink |
Temperature Range for Growth (℃) | 10–45 | 10–45 | 10–45 | 10–45 |
pH Range for Growth | 7–10 | 8–10 | 6–10 | 7–10 |
Oxidase | + | + | + | + |
Catalase | + | + | + | + |
Enzyme Activities (API 20NE): | | | | |
Reduction of Nitrate | + | + | - | + |
Gelatin Hydrolysis | - | - | - | + |
Phenylacetate | - | + | - | - |
Carbon-source utilization (BIOLOG GEN III Microplate): | | | | |
D-Fructose | - | + | - | + |
D-Serine | - | + | + | - |
Glycyl-L-Proline | - | + | + | + |
D-Glucuronic Acid | - | w | - | + |
D-Lactic Acid Methyl Ester | - | + | + | + |
Citric Acid | + | + | + | + |
⍺-Keto-Glutaric Acid | - | + | + | - |
D-Malic Acid | - | + | + | + |
L-Malic Acid | + | + | + | + |
𝛾-Amino-Butryric Acid | - | + | + | + |
β-Hydroxy-D, L-Butyric Acid | + | + | + | + |
⍺-Keto-Butyric Acid | + | - | + | w |
Acetoacetic Acid | - | + | - | + |
Propionic Acid | + | - | + | - |
Acetic Acid | + | + | + | + |
Formic Acid | + | + | + | + |
DNA G + C content (mol%) | 61.0 | 60.0 | 55.5 | 59.5 |
Emulsifier’s activity and plastic degradation potential
The results of the emulsification test experiment were shown in Fig. 3A. A positive result was observed of the strain NTOU-MSR1T, with an \({\text{E}}_{24}\) value of 19.91 ± 2.23 after one day of culturing of the supernatant of the bacterial culture. This finding indicated the potential production of the biosurfactants that was secreted into the medium. However, a decline in biosurfactant production was observed during a 4-day incubation period (Fig. 3), suggesting possible instability of the biosurfactants or a decrease in the specific metabolic pathway responsible for biosurfactant generation. Further investigation is required to elucidate the chemical structure and precise characteristics of the potential biosurfactants. Our research provides preliminary evidence for the presence of biosurfactants in the novel marine species NTOU-MSR1T, highlighting its potential as a candidate for further studies related to biosurfactant research.
For detecting LDPE degradation was carried out by the clear zone method. We found that strain NTOU-MSR1T formed detectable clear zones, whereas Staphylococcus epidermidis F-1, isolated from marine animals, did not form any (Fig. 4A). O. baumannii DSM 15594T also exhibited a slight clear zone around the colony (Fig. 4B). The clear zone experiments indicated that strain NTOU-MSR1T and O. baumannii DSM 15594T possess plastic degradation potential. However, Staphylococcus epidermidis F-1 showed negative results. Previous research has demonstrated that the clear zone method can be utilized to identify bacterial consortia capable of degrading low-density polyethylene (Joshi et al., 2022). Therefore, the candidate strain NTOU-MSR1T merits further research in the context of plastic degradation.