1. Introduction
Seaweed is an abundant resource in the ocean. It contains functional components such as polysaccharides, polyunsaturated fatty acids, proteins, fiber, vitamins, amino acids, peptides, lipids, and minerals [1,2,3,4] that have various applications. Edible seaweed represents a bountiful marine drug source with health advantages, such as reducing the risk of chronic diseases [5,6,7]. The consumption of seaweed has become increasingly popular. In addition to direct consumption, a significant portion of harvested seaweed is utilized for processed food or in manufacturing of active substances or drugs with remarkable biological properties, such as polysaccharides carrageenan, agar, alginate, and some oligosaccharides [5]. There is increasing interest in using seaweed for global food security and other products, and thus the global demand for seaweed has grown exponentially over the last 25 years [8].
In the wake of the attention attracted by seaweed, the seaweed-associated microbes that form functional symbionts with seaweed also are of interest [1,9]. The surface of seaweed supports an abundant diverse microbiota. The interplay between seaweed and its microbiome can affect seaweed in diverse ways, such as in the growth, the defense, or destruction mechanisms, especially in terms of the degradation of seaweed and the production of biologically active metabolites [10,11,12]. Strains that can degrade seaweed efficiently are highly important in the development of marine algae drugs.
During our previous studies of seaweed microbial populations, a number of bacteria were cultured, including a novel strain 5k15T. In this study, this new taxon was classified with polyphasic taxonomy, including biochemical, chemotaxonomic, and phylogenetic analysis methods. The polyphasic evidence indicated that the strain 5K15T represents a novel phylogenetic species of the genus Oceaniferula. The genus Oceaniferula was first described by Jin et al. in 2021 [13] and contains only one recognized species, Oceaniferula marina, up to now. The type strain of O. marina was N1E253T, which was isolated from marine sediment. N1E253T possesses a draft genome size of 5.07 Mbp. The strain has an assimilatory sulfate reduction pathway that can assimilate sulfate into hydrogen sulfide. It also can dissimilate nitrate into nitrite. Most of the coding sequences in the genome were classified as unknown function, indicating that this genus has a lot of potentially unknown features. The genus Oceaniferula belongs to the family Verrucomicrobiaceae order Verrucomicrobiales class Verrucomicrobiae, phylum Verrucomicrobiota. Phylum Verrucomicrobiota is a core group of the PVC superphylum (Planctomycetes–Verrucomicrobiota–Chlamydiae, as well as Lentisphaerae, Kirimatiellaeota, and some uncultured candidate phyla) [14]. Although the phylum Verrucomicrobiota was reported to be ubiquitous in different habitats, it belongs to difficult-to-culture bacteria. They are important and unusual from many perspectives [15], especially in terms of their brilliant potential to degrade the sulfated polysaccharides, which are an important component of seaweed [16,17]. The sulfated polysaccharides in seaweeds such as fucoidans are probably recalcitrant and more difficult to degrade. No sufficient and convenient sources with enzyme activity are available yet for the degradation of fucalean algae. The abundance of Verrucomicrobiota was even coupled with the algae bloom onset [17]. In our previous study, we also found that Verrucomicrobiota was highly abundant in kelp. This indicated that Verrucomicrobiota is a potential source for degrading various sulfated polysaccharides.
The strain 5k15T in this study is a new member of the genus Oceaniferula phylum Verrucomicrobiota isolated from the seaweed habitat. The present study was designed to determine the taxonomic status, ecological characteristics, and potential application of the Oceaniferula strain, isolate 5K15T, recovered from kelp collected from Li Island, RongCheng, China.
2. Results and Discussion
2.1. Isolation and Maintenance of Strain 5K15T
Strain 5K15T, showing a yellow-colored, circular, and smooth colony, was isolated from kelp cultivated on Li Island, Rongcheng, China (37.2661° N, 122.6138° E). The purified culture was stored at −80 °C in sterile 15% (v/v) glycerol supplemented with 1% (w/v) NaCl. Oceaniferula marina N1E253T (= KCTC 72800T = MCCC 1H00405T) was used as a related type strain for the comparison of phenotypic and chemotaxonomic characteristics.
2.2. 16S rRNA Gene Sequence and Phylogenetic Analysis
Analysis of the 16S rRNA gene sequences (1462 bp) clearly showed that the most closely related species to strain 5K15T was O. marina N1E253T (96.4%), followed by Persicirhabdus sediminis KCTC 22039T (91.9%), Rubritalea sabuli KCTC 22127T (91.8%), and R. profundi KCTC 52460T (91.7%). The topological structure of the phylogenetic maximum likelihood tree (Figure 1) clearly illustrated that strain 5K15T was affiliated to the genus Oceaniferula, supported by a high bootstrap value of 100%. The tree topology was also supported by FastTree and IQTree (Figures S1 and S2). Therefore, it would appear that, on the basis of the sequence identity and phylogenetic data, strain 5K15T may represent a novel species of the genus Oceaniferula.
2.3. Morphological, Physiological, and Biochemical Characteristics
Cells of 5K15T were Gram-stain-negative, facultatively anaerobic, non-flagellated, non-motile, rod-shaped, and approximately 0.3–0.5 µm in width and 0.7–2.0 µm in length (Figure S3). Gliding motility was not observed. Growth occurred at 15–37 °C (optimum, 33 °C), pH 6.0–8.5 (optimum, 7.5), and in the presence of 1.0–5.5% (w/v) NaCl (optimum, 2.0–3.5%). No growth occurred at pH 5.5 or pH 9.0. It was positive for nitrate reduction, but negative for catalase, oxidase, DNase, and hydrolysis of starch; casein; cellulose; and Tweens 20, 40, 60, and 80. Strain 5K15T was susceptible to ceftriaxone (30 µg), cefotaxime (30 µg), ampicillin (10 µg), carbenicillin (100 µg), penicillin (30 µg), and lincomycin (2 µg), but resistant to rifampin (5 µg), chloramphenicol (30 µg), clarithromycin (15 µg), erythromycin (15 µg), ofloxacin (5 µg), gentamycin (10 µg), tetracycline (30 µg), kanamycin (30 µg), tobramycin (10 µg), vancomycin (30 µg), neomycin (30 µg), norfloxacin (10 µg), and streptomycin (10 µg).
The strain was found to be positive for the Voges–Proskauer reaction, but negative for citrate utilization; arginine dihydrolase; lysine decarboxylase; ornithine decarboxylase; urease; tryptophan deaminase; and gelatinase, H2S, and indole production. Activities of esterase (C4), leucine arylamidase, valine arylamidase, acid phosphatase, and naphthol-AS-BI-phosphphydrolase were present, but activities of alkaline phosphatase, esterase lipase (C8), lipase (C14), cystine arylamidase, trypsin, α-chymotrypsin, α-galactosidase, β-galactosidase, β-glucuronidase, α-glucosidase, β-glucosidase, N-acetyl-β-glucosaminidase, α-mannosidase, and β-fucosidase were absent. Acid was produced from
2.4. Chemotaxonomic Characterization
The sole respiratory quinone of strain 5K15T was MK-9, which was similar to O. marina N1E253T. A list of the fatty acid compositions of strains 5K15T and N1E253T are detailed in Table 2. The fatty acid profile of strain 5K15T was similar to that of O. marina N1E253T, although there were differences in the proportions of some components. The major fatty acids (>10% of the total) of strain 5K15T were iso-C14:0 (28.5%), anteiso-C15: 0 (11.8%), and C16:0 (10.0%). The major polar lipids of strain 5K15T were phosphatidylethanolamine (PE), diphosphatidylglycerol (DPG), and one unidentified phospholipid (PL). Minor amounts of phosphatidylglycerol (PG), one unidentified aminolipid (AL), and one unidentified glycolipid (GL) were also present. These polar lipids were also present in O. marina N1E253T but with different contents. In addition, three unidentified polar lipids (L1, L2, L3) were only present in strain O. marina N1E253T. Detailed polar lipid images of different strains are represented in Figure S4.
2.5. Genome insight of the Two Strains of Genus Oceaniferula
2.5.1. Genome Properties
The draft genome sequences of strain 5K15T were assembled into 51 contigs, and the total genome length was 3,953,077 bp, which was smaller than that of O. marina N1E253T (5,073,947 bp) (Table 3). The sequencing depth of coverage was 363×. The genome DNA G+C content of strain 5K15T was 55.4 mol%, which was higher than that of O. marina N1E253T (52.0 mol%).
The draft genome of 5K15T was predicted to contain 3354 genes, including 3293 protein-coding genes, 45 tRNAs, 3 rRNAs, and 4 ncRNAs by the NCBI Prokaryotic Genome Annotation Pipeline server online (PGAP). According to annotations from the KEGG database, 1426 genes (43.3%) could be assigned a putative function, including genetic information processing (162 genes), carbohydrate metabolism (152 genes), signaling and cellular processes (139 genes), amino acid metabolism (87 genes), metabolism of cofactors and vitamins (84 genes), environmental information processing (76 genes), energy metabolism (63 genes), lipid metabolism (62 genes), glycan biosynthesis and metabolism (41 genes), and drug resistance (47 genes). The comparison of 5K15T and O. marina N1E253T is shown in Figure 2.
The ANI value between strain 5K15T and O. marina N1E253T was 77.1%. This value is below the threshold value (95.0%) for species delineation [18]. The dDDH value of strain 5K15T in comparison to O. marina N1E253T was lower than 20.0%. This value is much lower than the threshold for species demarcation (70.0%) [19]. The POCP value between strain 5K15T and O. marina N1E253T was 50.6%, which is higher than the genus boundary cut-off of 50% [20]. The phylogenetic tree based on a concatenated alignment of 120 ubiquitous single-copy genes showed the evolutionary relationships of strain 5K15T and related strains (Figure 3). Combined with the 16S rRNA gene sequence identity, these data indicate that strain 5K15T represents a novel species of the genus Oceaniferula.
It is, therefore, proposed that isolate 5K15T be considered as the type strain of a novel Oceaniferula species that belongs to the phylum Verrucomicrobiota. The name proposed for this type species is Oceaniferula flavus sp. nov.
2.5.2. Metabolic Pathways
Analysis of metabolic pathways was performed using KEGG’s BlastKOALA service. Many complete metabolic modules were found in the genomes of strain 5K15T and O. marina N1E253T, such as the complete central carbohydrate metabolism, energy metabolism, lipid metabolism, nucleotide metabolism, amino acid metabolism, and cofactor and vitamin metabolism (Figure 4). For example, the genomes of strain 5K15T and O. marina N1E253T both contain a complete set of genes of the arginine and proline metabolism pathway (M00028, M00844, and M00015) (Figure 5A). Both strains contain genes argJ, argB, argC, argD, and argJ, which support transformation from glutamate to ornithine. The presence of genes argF, argG, and argH allows them to use ornithine to synthesize arginine. The strain 5K15T also has the ability to perform isoleucine biosynthesis through two pathways modules (M00535 and M00570), but only O. marina N1E253T can biosynthesize isoleucine through one module (M00570). In addition, both O. marina N1E253T and strain 5K15T have genes cysN, cysD, cysC, cysH, cysJ, and cysI, which can assimilate sulfate into hydrogen sulfide (M00176) for energy metabolism (Figure 5B). As for other energy metabolism pathways, these two strains also possess a carbon fixation ability through CAM (Crassulacean acid metabolism) pathways (M00168). Only O. marina N1E253T has a phosphate acetyltransferase-acetate kinase pathway (M00579), which is also a kind of carbon fixation pathway. For cofactor and vitamin metabolism, unlike O. marina N1E253T, strain 5K15T lacks genes for glucose-1-phosphate thymidylyltransferase (EC 2.7.7.24) and dTDP-4-dehydrorhamnose reductase (EC 1.1.1.133); furthermore, the dTDP-L-rhamnose biosynthesis pathway was found to be incomplete (M00793).
The functional annotation of the two strains were also conducted by searching against the COG database. Overall, except for the unknown function, the function of carbohydrate transport and metabolism, inorganic ion transport and metabolism, amino acid transport and metabolism, and cell wall/membrane/envelope biogenesis were the major functional categories (Figure S5). The ability of carbohydrate transport and metabolism makes the strains able to live in kelp, which is rich in carbohydrates. On the other hand, the decomposition and utilization of kelp carbohydrates can lead to the decay of fresh kelp, which is subversive for the cultivation of kelp. However, considering the deep processing of kelp, the decomposition and utilization of kelp by the bacteria can contribute to the production of kelp degradation products.
2.5.3. Prediction of Secondary Metabolites
Secondary metabolite analysis was performed by antiSMASH. The genome of strain 5K15T was found to have four putative biosynthetic gene clusters (BGCs), including two terpene biosynthetic gene clusters, one aryl polyene biosynthetic cluster, and one type III PKS (polyketide synthase) cluster (Table S1). Aryl polyenes are structurally similar to carotenoids but are biosynthesized by a PKS rather than a terpene synthase, and they have been predicted to be protective agents in response to photo-oxidative damage, just like carotenoids [21]. Type III PKSs have been found in bacteria and fungi, as well as plants. Type III PKSs in microorganisms are involved in the biosynthesis of some lipidic materials and diverse secondary metabolites, which have significant biological functions and/or important pharmaceutical activities [22]. The strain O. marina N1E253T has two terpene biosynthetic clusters and one phosphonate cluster (Table S1). The potential ability of multiple BGCs may help the bacteria grow better in a marine environment [23].
2.5.4. Estimation of Carbohydrate-Active Enzymes (CAZymes)
Considering the potential polysaccharide utilization ability of Verrucomicrobiota, the carbohydrate-active enzymes (CAZymes) of the two Oceaniferula genomes were estimated using the CAZy database. Strain 5K15T and O. marina N1E253T contained 106 and 203 carbohydrate-active enzymes, respectively (Figure 6). In the genome of strain O. marina N1E253T, more than 70% (142 enzymes) of the identified enzymes were assigned to the glycoside hydrolase (GH) family. Compared with O. marina N1E253T, glycosyltransferases (GT) are the predominant enzyme in the strain 5K15T. More than 41.5%, 29.2%, and 16.0% of the enzymes were assigned to the GT family (44 enzymes), GH family (31 enzymes), and carbohydrate esterases (17 enzymes) in the genome of strain 5K15T, respectively. The 31 glycoside hydrolases include fucosidases from the families GH29 and GH95, as well as fucoidanase from the families GH36, GH116, and GH117 with a yet unknown function and GH168 with activity of endo-1,3-fucanase. In addition, according to the NCBI Prokaryotic Genome Annotation Pipeline (PGAP) and the Database of Sulfatases (
2.6. Description of Oceaniferula flavus sp. nov. Oceaniferula flavus (L. masc. adj. flavus, Yellow, Referring to the Color of the Colonies)
Cells are Gram-stain-negative, facultatively anaerobic, rod-shaped and 0.3–0.5 µm wide, and 0.7–2.0 µm long. Neither cellular gliding movement nor mobile ability was observed. Colonies grown on 1/2-strength R2A agar with 75% artificial seawater were circular, convex, and yellow. Growth occurred at 15–37 ℃ (optimum, 33 ℃), pH 6.0–8.5 (optimum, 7.5), and 1.0–5.5% (w/v) NaCl (optimum, 2.0–3.5% (w/v)). It was positive for nitrate reduction, but negative for catalase, oxidase, and DNase, as well as hydrolysis of starch; casein; cellulose; and Tweens 20, 40, 60, and 80. It was positive for the Voges–Proskauer reaction, but negative for citrate utilization; arginine dihydrolase; lysine decarboxylase; ornithine decarboxylase; urease; tryptophan deaminase; and gelatinase, H2S, and indole production. Activities of esterase (C4), leucine arylamidase, valine arylamidase, acid phosphatase, and naphthol-AS-BI-phosphphydrolase were present, but activities of alkaline phosphatase, esterase lipase (C8), lipase (C14), cystine arylamidase, trypsin, α-chymotrypsin, α-galactosidase, β-galactosidase, β-glucuronidase, α-glucosidase, β-glucosidase, N-acetyl-β-glucosaminidase, α-mannosidase, and β-fucosidase were absent. Acid was produced from
The type strain was 5K15T (= KCTC 82748T = MCCC 1H00442T = SDUM 810003T), isolated from kelp of Li Island, Rongcheng, China (37.2661° N, 122.6138° E). Genomic analysis indicated that the strain has great potential to degrade sulfated polysaccharides and can be applied in the marine drug industry.
3. Materials and Methods
3.1. Sample Collection, Isolation, and Preservation
In the process of screening for seaweed microbial populations, strain 5K15T was isolated from kelp that was cultivated in Li Island, Rongcheng, China (37.2661° N, 122.6138° E). The fresh kelp sample was cut into small pieces and placed in sterile aged seawater followed by shaking at 300 rpm for 1 h. The 100 μL supernatant was spread on a self-made kelp powder medium (10 g kelp powder, 15 g agar, 1 L aged seawater). Plates were incubated in an aerobic environment at 30 °C for 7 days. A yellow-colored colony was isolated and further subcultivated on 1/2-strength R2A agar with 75% artificial seawater. The purified culture was stored at −80 ℃ in sterile 15% (v/v) glycerol supplemented with 1% (w/v) NaCl. Oceaniferula marina N1E253T (= KCTC 72800T = MCCC 1H00405T) was used as the related type strain for the comparison of phenotypic and chemotaxonomic characteristics.
3.2. Physiological and Biochemical Characteristics
The cell morphology and the presence of flagella and cell size were observed by using light microscopy (E600; Nikon) and scanning electron microscopy (model Nova NanoSEM450; FEI). The motility was tested with the hanging-drop method, and gliding motility was confirmed by inoculating the bacteria on the 0.5% agar, as described by Bowman [24]. Gram staining was conducted by using a Gram stain kit (bioMérieux). Growth at various temperatures (4, 10, 15, 20, 25, 28, 30, 33, 35, 37, 40, and 45 ℃) and at various NaCl concentrations (0–10%, at intervals of 0.5%) was tested on 1/2-strength R2A agar with 75% artificial seawater in triplicate. The pH range for growth was determined by adjusting the medium with 1 M HCl or NaOH, as well as the following buffers (Sangon, China): MES (pH 5.5 and 6.0), PIPES (pH 6.5 and 7.0), HEPES (pH 7.5 and 8.0), Tricine (pH 8.5), and CAPSO (pH 9.0, 9.5) with a concentration of 20 mM. Anaerobic growth was determined on 1/2-strength R2A agar with 75% artificial seawater with or without 0.1% (w/v) KNO3 incubated in an anaerobic jar (10% H2, 10% CO2, and 80% N2) for 14 days. Susceptibility to antibiotics was investigated using the disc-diffusion method as described previously [25]. The hydrolysis of starch, casein, alginate, carboxymethylcellulose, and Tweens (20, 40, 60, and 80) was investigated according to previously mentioned methods [26]. DNase activity was determined on DNase test agar (Hopebio, Qingdao, China). The catalase activity was determined by using the bubble formation in 3% H2O2 solution. The oxidase activity was tested with an oxidase reagent kit (bioMérieux, Lyon, France) according to the manufacturer’s instructions. The production of other enzymes was assessed using API ZYM (bioMérieux, Lyon, France) kits. Tests for acid production from different carbohydrates were performed using the API 50CHB system (bioMérieux, Lyon, France). Tests for the oxidation of carbon sources were conducted using the Biolog GEN III system. Tests for additional physiological and biochemical characteristics were conducted using the API 20E system (bioMérieux, Lyon, France). All the API and Biolog tests were carried out according to the manufacturer’s instructions, except that the concentration of the suspension medium was adjusted to 3.0% (w/v) NaCl. All tests were implemented simultaneously with strains 5K15T and O. marina N1E253T with two replicates.
3.3. Chemotaxonomy
The cell mass of strain 5K15T was obtained from cultures grown for 3 days on 1/2-strength R2A agar with 75% artificial seawater at 33 °C. The biomass of strain O. marina N1E253T was obtained from cultures grown on marine agar 2216 (MA) for 3 days at 30 °C. The fatty acids were saponified, methylated, and extracted according to the standard protocol of MIDI (Sherlock Microbial Identification System, version 6.1). Cellular fatty acids were analyzed using an Agilent 6890N gas chromatograph and identified using the TSBA40 database of the microbial identification system [27]. Respiratory quinones were extracted from 300 mg freeze-dried cells [28] and purified with a silica gel TLC plate (Merck Kieselgel 60 F254), then analyzed via HPLC (high-performance liquid chromatography) [29]. Polar lipids were determined using 2D thin-layer chromatography (TLC) [30].
3.4. DNA Sequencing and Phylogenetic Analyses
Universal bacterial primers 27F and 1492R [31] were used for the amplification of the 16S rRNA gene. The PCR product was purified and cloned into the pMD18-T vector (Takara, Dalian, China), and recombinant plasmids were reproduced in Escherichia coli DH5α cells (Trans-Gen Biotech, Beijing, China). Sequencing was performed by BGI Tech (Qingdao, China). To identify the taxonomic status of strain 5K15T, the cloned sequence was submitted to the GenBank databases, and sequence identities were calculated from alignments obtained by the BLAST algorithm (
3.5. Genomic DNA Sequencing, Assembly, Annotation, and Genome Comparison
The genomic DNA of strain 5K15T was extracted using the bacterial genomic DNA Kit (Takara, Dalian, China). The draft whole-genome sequencing of the strain 5K15T was carried out by Novogene (Beijing, China), using the Illumina Hiseq PE150 platform. All the good-quality sequenced reads were assembled using SOAPdenovo software into several scaffolds [34]. Genes were annotated by NCBI, KEGG, and COG databases. For analysis of gene clusters encoding secondary metabolites, the antiSMASH 5.0 database was used [35]. The DNA G+C content of the strain was calculated from the genome sequence. The average amino acid identity (AAI) value (
4. Conclusions
A new member of the phylum Verrucomicrobiota, strain 5K15T, was recovered from kelp taken from Li Island in China, being shown to be most closely related to the type strain Oceaniferula marina N1E253T. The characterization of strain 5K15T using polyphasic methods suggests that it belongs to a new Oceaniferula species, which we named Oceaniferula flavus sp. nov. For general use, we propose the name Oceaniferula flavus due to the yellow color of the colonies. The genome sequence of the strain (genome size 3.95 Mbp) showed the presence of four putative biosynthetic gene clusters (BGCs), including two terpene biosynthetic gene clusters, one arylpolyene biosynthetic cluster, and one type III PKS cluster. Strain 5K15T and O. marina N1E253T contained 106 and 203 carbohydrate-active enzymes, including 17 and 13 carbohydrate esterases (CEs), 79 and 137 sulfatase encoded genes, 24 and 52 sulfatase-like hydrolase/transferase-encoded genes, and 25 and 36 arylsulfatase encoded genes, respectively, which are essential for dismantling the heavily decorated structure. The sulfated polysaccharides are the main components of seaweed and are challenging to degrade for bacteria compared to other monomers. The genus Oceaniferula of phylum Verrucomicrobiota provides candidates for the degradation of highly sulfated polysaccharides in the production of active metabolites from seaweeds. On the basis of the analysis of this research, we will perform experiments to confirm the degrading ability of genus Oceaniferula in the future. The structure and function of the seaweed degradation products by this genus will be illuminated deeply for the development of active drugs.
M.-Q.Y. designed the experiment and offered experiment guidance. M.-Q.Y. and C.-B.J. performed the experimental operation and data collection, and finished the manuscript. X.-J.L. and Y.-Q.Y. performed the material preparation and experiment. X.-Y.T. helped perform comparative genome analysis. Z.-J.D. offered experiment guidance and critical revision of manuscripts. All authors contributed to the article and approved the submitted version. All authors have read and agreed to the published version of the manuscript.
Not applicable.
The GenBank accession numbers for the 16S rRNA gene sequence and draft genome sequence of strain 5K15T were MW260263 and JAFBGL000000000, respectively.
The scanning electron microscopy was supported by the Physical-Chemical Materials Analytical and Testing Center of Shandong University at Weihai.
The authors declare that they do not have any conflict of interest.
Footnotes
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Enlarge this image.Figure 1. Maximum likelihood phylogenetic tree based on 16S rRNA gene sequences showing the phylogenetic position of strain 5K15T among closely related taxa. Bootstrap values (expressed as percentages of 1000 replications) of >70% are shown at branch points. Escherichia coli JCM 1649T (AB242910) was used as an outgroup. Bar, 0.05 nucleotide substitutions per position.
Enlarge this image.Figure 2. Bar diagram of annotated KEGG functional categories of strains 5K15T and O. marina N1E253T.
Enlarge this image.Figure 3. The IQ-TREE based on 120 ubiquitous single-copy genes. Escherichia coli JCM 1649T was used as the outgroup. Bar, 0.20 substitutions per nucleotide position.
Enlarge this image.Figure 4. Heat maps of complete and incomplete metabolic pathways in the genomes of strains 5K15T and O. marina N1E253T.
Enlarge this image.Figure 5. Amino acid metabolism (A) and assimilatory sulfate reduction pathway (B) of the genus Oceaniferula.
Enlarge this image.Figure 6. Histogram of predicted carbohydrate-active enzymes in strains 5K15T and O. marina N1E253T. AAs, auxiliary activities; CBMs, carbohydrate-binding module; CEs, carbohydrate esterases; GHs, glycoside hydrolases; GTs, glycosyltransferases; PLs, polysaccharide lyases.
Differential characteristics of strain 5K15T and the related type strain.
| Characteristic | 5K15T | O. marina N1E253T |
|---|---|---|
| Cell size (μm) | 0.3–0.5 × 0.7–2.0 | 0.5–0.8 × 1.5–2.0 a |
| Temperature range (°C) | 15–37 | 15–37 a |
| Growth pH | 6.0–8.5 | 7.0–9.0 a |
| Tolerance to NaCl (%) | 1.0–5.5 | 1.5–5.0 a |
| Catalase | – | + a |
| Oxidase | – | + a |
| Hydrolysis of starch | – | + a |
| API 20E test: | ||
| Tryptophan deaminase | – | + |
| Voges–Proskauer reaction | + | + |
| Oxidation of (Biolog GEN III): | ||
| – | + | |
| + | – | |
| + | – | |
| + | – | |
| Enzyme activities (API ZYM): | ||
| Leucine arylamidase | + | + |
| N-acetyl-β-glucosaminidase | – | + |
| Acid production (API 50CHB): | ||
| – | + | |
| – | + | |
| + | – | |
| + | – | |
| DNA G+C content (mol%) | 55.4 | 52.0 |
All data are from this study unless otherwise indicated. +, positive; –, negative. Data from: a Jin et al. (2021) [
Fatty acid compositions (%) of strain 5K15T and the related strain.
| Fatty Acid | 5K15T | O. marina N1E253T |
|---|---|---|
| Branched fatty acids | ||
| Iso-C12:0 | 6.7 | - |
| Anteiso-C13: 0 | 1.4 | - |
| Iso-C14:0 | 28.5 | 33.3 |
| Anteiso-C15: 0 | 11.8 | 4.9 |
| Iso-C16:0 | 9.5 | 7.3 |
| Straight-chain fatty acids | ||
| C12:0 | Tr | 1.2 |
| C14:0 | Tr | 1.4 |
| C16:0 | 10.0 | 12.9 |
| C17:0 | 4.0 | 7.0 |
| C18:0 | Tr | 1.3 |
| Unsaturated fatty acids | ||
| C15: 1ω6c | 2.2 | 1.5 |
| C15: 1ω8c | - | 1.9 |
| C17: 1ω6c | Tr | 4.2 |
| C17: 1ω8c | 1.2 | 6.2 |
| Hydroxy fatty acids | ||
| C12: 0 3-OH | 2.0 | Tr |
| Iso-C12: 0 3-OH | 2.3 | - |
| Iso-C14: 0 3-OH | Tr | 3.8 |
| Summed features * | ||
| Summed feature 3 | 9.3 | 7.0 |
Tr, traces (<1.0%); -, not detected. Fatty acids amounting to (<1.0%) are not shown; dominant fatty acids (≥10.0%) are highlighted in bold. All data are from this study. * Summed features are groups of two or three fatty acids that are treated together for the purpose of evaluation in the MIDI system and include both peaks with discrete equivalent chain lengths (ECLs) as well as those where the ECLs are not reported separately. Summed feature 3 was listed as C16:1ω7c and/or C16:1ω6c.
Genome statistics of strain 5K15T and the related type strain.
| Characteristic | 5K15T | O. marina N1E253T |
|---|---|---|
| GenBank accession number | JAFBGL000000000 | JACBAZ000000000 |
| Genome size (bp) | 3,953,077 | 5,073,947 |
| N50 Value (bp) | 346,557 | 785,939 |
| Contigs (no.) | 51 | 88 |
| G+C content (mol/%) | 55.4 | 52 |
| Annotation Results by the NCBI Prokaryotic Genome Annotation Pipeline (PGAP) | ||
| Genes (no.) | 3354 | 3927 |
| Protein-coding genes | 3293 | 3855 |
| tRNAs (no.) | 45 | 46 |
| rRNAs (no.) | 3 | 7 |
| ncRNAs (no.) | 4 | 3 |
Supplementary Materials
The following supporting information can be downloaded at
References
1. Del Olmo, A.; Picon, A.; Nuñez, M. The microbiota of eight species of dehydrated edible seaweeds from North West Spain. Food Microbiol.; 2018; 70, pp. 224-231. [DOI: https://dx.doi.org/10.1016/j.fm.2017.10.009]
2. Xu, S.; Huang, X.; Cheong, K. Recent advances in marine algae polysaccharides: Isolation, structure, and activities. Mar. Drugs; 2017; 15, 388. [DOI: https://dx.doi.org/10.3390/md15120388]
3. Echave, J.; Otero, P.; Garcia-Oliveira, P.; Munekata, P.E.S.; Pateiro, M.; Lorenzo, J.M.; Simal-Gandara, J.; Prieto, M.A. Seaweed-derived proteins and peptides: Promising marine bioactives. Antioxidants; 2022; 11, 176. [DOI: https://dx.doi.org/10.3390/antiox11010176]
4. Rawiwan, P.; Peng, Y.; Paramayuda, I.G.P.B.; Quek, S.Y. Red seaweed: A promising alternative protein source for global food sustainability. Trends. Food Sci. Technol.; 2022; 123, pp. 37-56. [DOI: https://dx.doi.org/10.1016/j.tifs.2022.03.003]
5. Brown, E.M.; Allsopp, P.J.; Magee, P.J.; Gill, C.I.; Nitecki, S.; Strain, C.R.; McSorley, E.M. Seaweed and human health. Nutr. Rev.; 2014; 72, pp. 205-216. [DOI: https://dx.doi.org/10.1111/nure.12091] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24697280]
6. Meinita, M.D.N.; Harwanto, D.; Choi, J. Seaweed exhibits therapeutic properties against chronic diseases: An overview. Appl. Sci.; 2022; 12, 2638. [DOI: https://dx.doi.org/10.3390/app12052638]
7. Shah, M.D.; Venmathi Maran, B.A.; Shaleh, S.R.M.; Zuldin, W.H.; Gnanaraj, C.; Yong, Y.S. Therapeutic potential and nutraceutical profiling of north bornean seaweeds: A review. Mar. Drugs; 2022; 20, 101. [DOI: https://dx.doi.org/10.3390/md20020101]
8. Msuya, F.E.; Bolton, J.; Pascal, F.; Narrain, K.; Nyonje, B.; Cottier-Cook, E.J. Seaweed farming in Africa: Current status and future potential. J. Appl. Phycol.; 2022; 34, pp. 985-1005. [DOI: https://dx.doi.org/10.1007/s10811-021-02676-w]
9. Ren, C.G.; Liu, Z.Y.; Wang, X.L.; Qin, S. The seaweed holobiont: From microecology to biotechnological applications. Microb. Biotechnol.; 2022; 15, pp. 738-754. [DOI: https://dx.doi.org/10.1111/1751-7915.14014]
10. Ramanan, R.; Kang, Z.; Kim, B.; Cho, D.; Jin, L.; Oh, H.; Kim, H.S. Phycosphere bacterial diversity in green algae reveals an apparent similarity across habitats. Algal Res.; 2015; 8, pp. 140-144. [DOI: https://dx.doi.org/10.1016/j.algal.2015.02.003]
11. Grant, M.A.; Kazamia, E.; Cicuta, P.; Smith, A.G. Direct exchange of vitamin B12 is demonstrated by modelling the growth dynamics of algal–bacterial cocultures. ISME J.; 2014; 8, pp. 1418-1427. [DOI: https://dx.doi.org/10.1038/ismej.2014.9]
12. Bengtsson, M.M.; Sjotun, K.; Lanzen, A.; Ovreas, L. Bacterial diversity in relation to secondary production and succession on surfaces of the kelp Laminaria hyperborea. ISME J.; 2012; 6, pp. 2188-2198. [DOI: https://dx.doi.org/10.1038/ismej.2012.67]
13. Jin, C.B.; Feng, X.; Zou, Q.H.; Ye, M.Q.; Du, Z.J. Oceaniferula marina gen. nov., sp. nov., an anti-fluoroquinolone bacterium isolated from marine sediment. Antonie Leeuwenhoek; 2021; 114, pp. 1855-1865. [DOI: https://dx.doi.org/10.1007/s10482-021-01645-0] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34468960]
14. Devos, D.P.; Ward, N.L. Mind the PVCs. Environ. Microbiol.; 2014; 16, pp. 1217-1221. [DOI: https://dx.doi.org/10.1111/1462-2920.12349]
15. Fuerst, J.A. Phylum Verrucomicrobia. Reference Module in Life Sciences; Roitberg, B.D. Elsevier: Amsterdam, The Netherlands, 2019; pp. 551-563.
16. Sichert, A.; Corzett, C.H.; Schechter, M.S.; Unfried, F.; Markert, S.; Becher, D.; Fernandez-Guerra, A.; Liebeke, M.; Schweder, T.; Polz, M.F. et al. Verrucomicrobia use hundreds of enzymes to digest the algal polysaccharide fucoidan. Nat. Microbiol.; 2020; 5, pp. 1026-1046. [DOI: https://dx.doi.org/10.1038/s41564-020-0720-2]
17. Orellana, L.H.; Francis, T.B.; Ferraro, M.; Hehemann, J.; Fuchs, B.M.; Amann, R.I. Verrucomicrobiota are specialist consumers of sulfated methyl pentoses during diatom blooms. ISME J.; 2021; 16, pp. 630-641. [DOI: https://dx.doi.org/10.1038/s41396-021-01105-7] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34493810]
18. Richter, M.; Rosselló-Móra, R. Shifting the genomic gold standard for the prokaryotic species definition. Proc. Natl. Acad. Sci. USA; 2009; 106, pp. 19126-19131. [DOI: https://dx.doi.org/10.1073/pnas.0906412106] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19855009]
19. Goris, J.; Konstantinidis, K.T.; Klappenbach, J.A.; Coenye, T.; Vandamme, P.; Tiedje, J.M. DNA-DNA hybridization values and their relationship to whole-genome sequence similarities. Int. J. Syst. Evol. Microbiol.; 2007; 57, pp. 81-91. [DOI: https://dx.doi.org/10.1099/ijs.0.64483-0]
20. Qin, Q.L.; Xie, B.B.; Zhang, X.Y.; Chen, X.L.; Zhou, B.C.; Zhou, J.; Oren, A.; Zhang, Y.Z. A proposed genus boundary for the prokaryotes based on genomic insights. J. Bacteriol.; 2014; 196, pp. 2210-2215. [DOI: https://dx.doi.org/10.1128/JB.01688-14]
21. Jones, C.V.; Jarboe, B.G.; Majer, H.M.; Ma, A.T.; Beld, J. Escherichia coli Nissle 1917 secondary metabolism: Aryl polyene biosynthesis and phosphopantetheinyl transferase crosstalk. Appl. Microbiol. Biotechnol.; 2021; 105, pp. 7785-7799. [DOI: https://dx.doi.org/10.1007/s00253-021-11546-x]
22. Katsuyama, Y.; Ohnishi, Y. Chapter Sixteen—Type III Polyketide Synthases in Microorganisms. Methods in Enzymolog; Hopwood, D.A. Academic Press: Salt Lake City, UT, USA, 2012.
23. Lim, Y.; Kang, I.; Cho, J. Genome characteristics of Kordia antarctica IMCC3317T and comparative genome analysis of the genus Kordia. Sci. Rep.; 2020; 10, 14715. [DOI: https://dx.doi.org/10.1038/s41598-020-71328-9]
24. Bowman, J.P. Description of Cellulophaga algicola sp. nov., isolated from the surfaces of Antarctic algae, and reclassification of Cytophaga uliginosa (ZoBell and Upham 1944) Reichenbach 1989 as Cellulophaga uliginosa comb. nov. Int. J. Syst. Evol. Microbiol.; 2000; 50, pp. 1861-1868. [DOI: https://dx.doi.org/10.1099/00207713-50-5-1861]
25. Bauer, A.; Kirby, W.; Sherris, J.; Turck, M. Antibiotic susceptibility testing by a standardized single disk method. Am. J. Clin. Pathol.; 1966; 45, pp. 493-496. [DOI: https://dx.doi.org/10.1093/ajcp/45.4_ts.493]
26. Cowan, S.T.; Steel, K.J. Bacterial characters and characterization. Cowan and Steel’s Manual for the Identification of Medical Bacteria; 2nd ed. Cowan, S.T. Cambridge University Press: Cambridge, UK, 1974.
27. Sasser, M. Identification of Bacteria by Gas Chromatography of Cellular Fatty Acids, MIDI Technical Note 101; Microbial MIDI Inc.: Newark, DE, USA, 1990.
28. Minnikin, D.E.; O’Donnell, A.G.; Goodfellow, M.; Alderson, G.; Athalye, M.; Schaal, A.; Parlett, J.H. An integrated procedure for the extraction of bacterial isoprenoid quinones and polar lipids. J. Microbiol. Meth.; 1984; 2, pp. 233-241. [DOI: https://dx.doi.org/10.1016/0167-7012(84)90018-6]
29. Kroppenstedt, R.M. Separation of bacterial menaquinones by HPLC using reverse phase (RP18) and a silver loaded ion exchanger as stationary phases. J. Liq. Chromatogr.; 1982; 5, pp. 2359-2367. [DOI: https://dx.doi.org/10.1080/01483918208067640]
30. Bligh, E.G.; Dyer, W.J. A rapid method of total lipid extraction and purification. Can J. Biochem. Physiol.; 1959; 37, pp. 911-917. [DOI: https://dx.doi.org/10.1139/y59-099]
31. Liu, Q.Q.; Wang, Y.; Li, J.; Du, Z.J.; Chen, G.J. Saccharicrinis carchari sp. nov., isolated from a shark, and emended descriptions of the genus Saccharicrinis and Saccharicrinis fermentans. Int. J. Syst. Evol. Microbiol.; 2014; 64, pp. 2204-2209. [DOI: https://dx.doi.org/10.1099/ijs.0.061986-0] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24695056]
32. Luke, S. Two fast tree-creation algorithms for genetic programming. IEEE Trans. Evol. Comput.; 2000; 4, pp. 274-283. [DOI: https://dx.doi.org/10.1109/4235.873237]
33. Trifinopoulos, J.; Nguyen, L.; von Haeseler, A.; Minh, B.Q. W-IQ-TREE: A fast online phylogenetic tool for maximum likelihood analysis. Nucleic Acids Res.; 2016; 44, pp. W232-W235. [DOI: https://dx.doi.org/10.1093/nar/gkw256]
34. Li, R.Q.; Yu, C.; Li, Y.R.; Lam, T.W.; Yiu, S.M.; Kristiansen, K.; Wang, J. SOAP2: An improved ultrafast tool for short read alignment. Bioinformatics; 2009; 25, pp. 1966-1967. [DOI: https://dx.doi.org/10.1093/bioinformatics/btp336] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19497933]
35. Blin, K.; Shaw, S.; Steinke, K.; Villebro, R.; Ziemert, N.; Lee, S.Y.; Medema, M.H.; Weber, T. antiSMASH 5.0: Updates to the secondary metabolite genome mining pipeline. Nucleic Acids Res.; 2019; 47, pp. W81-W87. [DOI: https://dx.doi.org/10.1093/nar/gkz310]
36. Chaumeil, P.A.; Mussig, A.J.; Hugenholtz, P.; Parks, D.H. GTDB-Tk: A toolkit to classify genomes with the Genome Taxonomy Database. Bioinformatics; 2019; 36, pp. 1-3. [DOI: https://dx.doi.org/10.1093/bioinformatics/btz848] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31730192]
37. Nguyen, L.T.; Schmidt, H.A.; von Haeseler, A.; Minh, B.Q. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol.; 2015; 32, pp. 268-274. [DOI: https://dx.doi.org/10.1093/molbev/msu300] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25371430]
38. Valencia, R.; Gonzalez, V.; Undabarrena, A.; Zamora-Leiva, L.; Ugalde, J.A.; Camara, B. An integrative bioinformatic analysis for keratinase detection in marine-derived Streptomyces. Mar. Drugs; 2021; 19, 286. [DOI: https://dx.doi.org/10.3390/md19060286] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34063876]
39. Konstantinou, D.; Popin, R.V.; Fewer, D.P.; Sivonen, K.; Gkelis, S. Genome reduction and secondary metabolism of the marine sponge-associated Cyanobacterium Leptothoe. Mar. Drugs; 2021; 19, 298. [DOI: https://dx.doi.org/10.3390/md19060298] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34073758]
© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Details
1 Marine College, Shandong University, Weihai 264209, China; Weihai Research Institute of Industrial Technology, Shandong University, Weihai 264209, China
2 Marine College, Shandong University, Weihai 264209, China
3 Marine College, Shandong University, Weihai 264209, China; Weihai Research Institute of Industrial Technology, Shandong University, Weihai 264209, China; State Key Laboratory of Microbial Technology, Shandong University, Qingdao 266237, China