Phylogenetic characteristics
The 16S rRNA gene sequence (1511 bp; GenBank accession number MZ695036) showed that the closely related species of strain D16T were genus or members of the genera Eudoraea adriatica DSM 19308T (93.0 %), Poritiphilus flavus R33T (92.5 %), Muriicola (91.9-92.4 %), Zeaxanthinibacter (91.7-92.3 %), Robiginitalea (91.9-92.3 %), Cellulophaga (91.2-91.7), Maribacter (90.9-92.0 %) and Arenibacter (90.9-91.5 %) in the family Flavobacteriaceae. All the similarities between strain D16T and closely related genera were lower than 94 %, which was considered a reasonable lower cut-off window for describing a new genus (Yarza et al. 2014; Varghese et al. 2015; Konstantinidis et al. 2017; Johnson et al. 2019). The neighbor-joining phylogenetic tree (Supplementary Figure 1) showed that strain D16T belonged to the clade represented by the closely related genera listed above. Still, strain D16T formed a single lineage in this clade with a high bootstrap value (97 %) (Fig. 1). In the maximum-likelihood tree (Supplementary Figure 2) based on single-copy orthologous clusters (concatenated protein sequences), like in neighbor-joining, maximum-likelihood and maximum-parsimony trees (Supplementary Figure 3) based on 16S rRNA gene sequences, strain D16T formed a separate branch, at the root of the branch formed by member of the genus Flavobacterium, within the cluster of the family Flavobacteriaceae, supporting its affiliation with a novel genus. A phylogenomic tree built based on genome sequences strain D16T and relative genera showed the same phylogenetic position of strain D16T (Supplementary Figure 4).
Genome composition and Genomes comparison
The genome size of strain D16T was 3,225,638 bp. According to the NCBI PGAP, strain D16T was predicted 2,981 genes in the genome, including 2,916 protein-coding genes. The amino acid sequence of strain D16T and close relatives were analyzed by KEGG and 1374 genes (46.3%) could be assigned a putative function (Supplementary Figure 5). Strain D16T contains 184 strain-specific genes. The numbers of gene functional categories comparison of strain D16T and Eudoraea adriatica DSM 19308T, Robiginitalea biformata HTCC2501T, Zeaxanthinibacter enoshimensis TD-ZE3T, Cellulophaga tyrosinoxydans EM41T and Muriicola jejuensis EM44T were shown in Supplementary Figure 6.
The genomic DNA G+C content of strain D16T was 42.8 mol %, which was in the range (37.8–55.3 mol%) of DNA G+C contents in the family Flavobacteriaceae (Table 1). The ANI, dDDH, AAI values of strain D16T compared to similar genera ranged from 67.3 to 70.2 %, 16.2 to 19.1 % and 68.3 to 70.1 %, respectively (Table 2), below the threshold for recommended species division of ANI 95–96%, dDDH 70% and AAI 95% (Lachance et al. 2020; Kim et al. 2021; Meier-Kolthoff et al. 2022). The heat map results of DNA-DNA hybridization (dDDH) and average nucleotide identity (ANI) values of genomes between D16T and related strains are presented in Supplementary Figure 7. The AAI-profiler comparison of amino acids of D16T results suggested that the closest sequenced proteome was from Eudoraea adriatica DSM 19308T, with 65% matched fraction and 68.32% average amino acid identity values (Supplementary Figure 8). These data indicated that strain D16T represents a novel species. Furthermore, POCP values between strain D16T and similar genera were less than 50.0%, consistent with the description of a new genus (Qin et al. 2014). Further supported that strain D16T was a novel genus of the family Flavobacteriaceae.
Table 1. Genome statistics of strain D16T and the related genera in family Flavobacteriaceae.
Strains: 1, D16T; 2, Eudoraea adriatica DSM 19308T(Alain et al. 2008); 3, Zeaxanthinibacter enoshimensis TD-ZE3T(Asker et al. 2007); 4, Robiginitalea biformata HTCC2501T(JC and SJ 2004); 5, Muriicola jejuensis EM44T(Kahng et al. 2010); 6, Cellulophaga tyrosinoxydans EM41T(Kahng et al. 2009); 7, Poritiphilus flavus R33T(Wang et al. 2020)
characteristics
|
1
|
2
|
3
|
4
|
5
|
6
|
7
|
Size of genome (mb)
|
3.1
|
3.91
|
3.34
|
3.53
|
3.31
|
3.56
|
4.99
|
N50 value(bp)
|
281,706
|
562,897
|
709,771
|
3,530,383
|
1,788,650
|
377,537
|
550,540
|
G+C content (mol%)
|
37.8
|
38.3
|
46.4
|
55.3
|
48.3
|
33.2
|
44.5
|
Contigs(no.)
|
20
|
22
|
13
|
1
|
26
|
18
|
19
|
Annotation by the NCBI Prokaryotic Genome Annotation Pipeline (PGAP)
|
|
Genes(no.)
|
2981
|
3476
|
2953
|
3070
|
2921
|
3166
|
4,191
|
Protein-coding genes (no.)
|
2916
|
3401
|
2897
|
3131
|
2994
|
3106
|
4118
|
Gene specific to genome
|
183
|
79
|
48
|
64
|
34
|
120
|
56
|
tRNA(no.)
|
40
|
36
|
38
|
41
|
39
|
33
|
3
|
rRNA(no.)
|
3
|
9
|
3
|
6
|
3
|
6
|
38
|
ncRNAs (no.)
|
4
|
4
|
4
|
4
|
4
|
4
|
4
|
DDBJ/ENA/GenBank accession number
|
JAJTJC000000000
|
KB907546
|
NZ_SNYI00000000
|
NC_013222.1
|
NZ_JAABOP000000000.1
|
NZ_FWXO01000000
|
NZ_WXYO01000000
|
Table 2 AAI, ANI, POCP, DDH values among strain D16T and relative genera of the family Flavobacteriaceae
Strain
|
Accession no.
|
AAI/%
|
ANI/%
|
POCP/%
|
DDH/%
|
Reference
|
D16T
|
JAJTJC000000000
|
100.0
|
100.0
|
100.0
|
100.0
|
This study
|
Robiginitalea biformata HTCC2501T
|
KB907546
|
71.4
|
69.6
|
71.4
|
16.5
|
(JC and SJ, 2004)
|
Muriicola jejuensis EM44T
|
NZ_SNYI00000000
|
70.1
|
68.9
|
67.9
|
17.0
|
(Kahng et al., 2010)
|
Eudoraea adriatica DSM 19308T
|
NC_013222.1
|
68.3
|
68.4
|
66.3
|
19.2
|
(Alain et al., 2008)
|
Zeaxanthinibacter enoshimensis TD-ZE3T
|
NZ_JAABOP000000000.1
|
70.4
|
69.6
|
75.1
|
16.8
|
(Asker et al., 2007)
|
Cellulophaga tyrosinoxydans EM41T
|
NZ_FWXO01000000
|
68.7
|
68.4
|
62.4
|
16.2
|
(Kahng et al., 2009)
|
Saonia flava DSM 29762T
|
NZ_JAATJJ000000000
|
70.1
|
69.2
|
62.6
|
16.8
|
(Fagervold et al., 2017)
|
Maribacter polysiphoniae KCTC 22021T
|
NZ_JACWLN000000000
|
70.3
|
69.3
|
22.7
|
16.8
|
(Nedashkovskaya et al., 2007)
|
Robiginitalea sediminis O458T
|
NZ_NGNR00000000
|
68.2
|
68.4
|
71.3
|
16.4
|
(Zhang et al., 2018b)
|
Robiginitalea myxolifaciens DSM 21019T
|
NZ_FOYQ00000000
|
67.5
|
68.2
|
69.2
|
18.0
|
(Manh et al., 2008)
|
Maribacter aquivivus DSM 16478T
|
NZ_FQZX00000000
|
69.2
|
68.0
|
53.2
|
16.6
|
(Nedashkovskaya et al., 2004)
|
Maribacter algicola PoM-212T
|
NZ_QUSX00000000
|
68.7
|
68.5
|
62.3
|
16.9
|
(Khan et al., 2020)
|
Arenibacter troitsensis DSM 19835T
|
NZ_FXAO00000000
|
69.3
|
68.9
|
55.4
|
17.3
|
(Nedashkovskaya et al., 2003)
|
Maribacter luteus RZ05T
|
NZ_WKJH01000030
|
69.9
|
69.2
|
59.6
|
17.7
|
(Liu et al., 2020)
|
Maribacter arenosus CAU 1321T
|
NZ_JABTCG010000000
|
70.0
|
69.0
|
65.0
|
18.0
|
(Thongphrom et al., 2016)
|
Maribacter aurantiacus KCTC 52409T
|
NZ_VBUK00000000
|
68.9
|
68.5
|
60.9
|
17.1
|
(Khan et al., 2020)
|
Maribacter orientalis DSM 16471T
|
NZ_FNZN00000000
|
68.1
|
68.0
|
61.2
|
17.2
|
(Nedashkovskaya et al., 2004)
|
Cellulophaga algicola DSM 14237T
|
NC_014934
|
68.1
|
68.1
|
51.0
|
18.4
|
(Bowman, 2000)
|
Maribacter flavus KCTC 42508T
|
NZ_VUOE00000000
|
69.0
|
68.7
|
61.3
|
17.3
|
(Tang et al., 2015)
|
Maribacter vaceletii DSM 25230T
|
NZ_RBIQ00000000
|
68.9
|
68.9
|
57.1
|
16.3
|
(Jackson et al., 2015)
|
Arenibacter aquaticus GUOT
|
NZ_RQPJ00000000
|
69.5
|
68.9
|
54.4
|
18.1
|
(Guo et al., 2020)
|
Arenibacter nanhaiticus CGMCC 1.8863T
|
NZ_SNZW00000000
|
68.0
|
68.2
|
60.4
|
17.5
|
(Sun et al., 2010)
|
Maribacter caenipelagi CECT 8455T
|
NZ_FNBD00000000
|
68.4
|
68.3
|
41.8
|
16.5
|
(Jung et al., 2014)
|
Cellulophaga baltica DSM 24729T
|
PRJNA599484
|
71.2
|
69.8
|
61.2
|
17.5
|
(Johansen et al., 1999)
|
Poritiphilus flavus R33T
|
GCA_009901585.1
|
71.1
|
68.3
|
61.2
|
17.5
|
(Wang et al., 2020)
|
Genome Function prediction
The RAST annotation indicated that gene functions of strain D16T covered all aspects of cellular metabolism, with the most diverse gene types and number associated with biological process (917 metabolic pathways). Statistics of metabolic function of strain D16T are shown in Supplementary Figure 9. About 110 carbon metabolism related pathways was annotated in genome of strain D16T, which also includes relavite genes and enzymes of CO2 fixation, such as Ss,RBCS,RBCl,ClCP,CA. There are13 metabolism pathways associated with nitrogen metabolism was annotated in genome of D16T, which including denitrification pathway, denitrification-related genes, such as narH and narI, which was consistent with the results of anaerobic experiments. Furthermore, 33 metabolism pathways associated with cellular respiration were found in the strain D16T. [NiFe], One of the special key enzymes, also is closely related with carbon fixation. the strain D16T may be involved in carbon and nitrogen cycling in marine sediments, predict from the above analysis. CO2 fixation、denitrification and [NiFe] was also found in other relative genera of family Flavobacteriaceae [Supplementary table 2].
Metabolic pathway comparison
Some complete metabolic pathways have been found in strain D16T and close relatives (Supplementary Figure 10). For Fatty acid biosynthesis, beta-Oxidation, acyl-CoA synthesis, biosynthesis of Inosine monophosphate, adenine ribonucleotide, guanine ribonucleotide and pyrimidine deoxyribonuleotide. These pathways were found in all these strains. Conversely, reductive pentose phosphate cycle, uridine monophosphate biosynthesis, ethylene biosynthesis and urea cycle present in all these strains were found to be incomplete.
Some metabolic pathways exhibited striking differences among these Flavobacteriaceae strains (Figure 2). For instance, Strain D16T, Maribacter polysiphoniae KCTC 22021T and Arenibacter troitsensis DSM 19835T lack complete Pentose phosphate pathway; only Eudoraea adriatica DSM 19308T and Arenibacter troitsensis DSM 19835T has a complete D-Galacturonate degradation; only Eudoraea adriatica DSM 19308T and Muriicola jejuensis EM44T has complete Glycogen biosynthesis, for Nucleotide sugar biosynthesis, only Maribacter polysiphoniae KCTC 22021T has a complete pathway; Zeaxanthinibacter enoshimensis TD-ZE3T, Robiginitalea biformata HTCC2501T, Muriicola jejuensis EM44T and Poritiphilus flavus R33T has a complete Propanoyl-CoA metabolism; only Eudoraea adriatica DSM 19308T and Cellulophaga tyrosinoxydans EM41T has complete Formaldehyde assimilation; For Amino acid metabolism, Valine biosynthesis, Leucine biosynthesis and Lysine biosynthesis are found in all these strains; Maribacter polysiphoniae KCTC 22021T lacks a complete Lysine biosynthesis; Strain D16T, Maribacter aquivivus DSM 16478T and Arenibacter troitsensis DSM 19835T has complete Methionine degradation; D16T lacks complete Methionine degradation, Zeaxanthinibacter enoshimensis TD-ZE3T, Robiginitalea biformata HTCC2501T, Cellulophaga tyrosinoxydans EM41T and Arenibacter troitsensis DSM 19835T has complete Proline biosynthesis, only Maribacter polysiphoniae KCTC 22021T lacks Histidine biosynthesis/degradation. Complete metabolism pathways which were only found in the strain D16 were TCA cycle, second carbon oxidation, Methionine degradation, Isoleucine biosynthesis, NAD biosynthesis and Heme biosynthesis.
Secondary Metabolite Biosynthesis Clusters Analysis
Clusters related to secondary metabolite biosynthesis were identified using the antiSMASH program and found differences in gene cluster abundance and diversity between strain D16T and close relative genera (Supplementary Figure 11). Common gene clusters encoding of strain D16T, Eudoraea adriatica DSM 19308T, Robiginitalea biformata HTCC2501T, Muriicola jejuensis EM44T and Zeaxanthinibacter enoshimensis TD-ZE3T were T3PKS: Chal_sti_synt_C, T3PKS: Chal_sti_synt_N, terpene: Lycopene_cycl, terpene: phytoene_synt, 2OG-FeII_Oxy, adh_short_C2, Alpha-amylase, PP-binding, RimK, t2fas, short-chain dehydrogenase/reductase SDR, aminotransferase class-III, oxidoreductase, Polyprenyl synthetase, BC transporter ATP-binding protein, TetR family transcriptional regulator, ND family efflux transporter MFP subunit and dehydrogenase. Gene clusters found in strain D16T were arylpolyene: APE_KS2, 8-amino-7-oxononanoate synthase, GATase_7, Aminotran_3, Peptidase_M16_C, PF04055 and PF06968. Based ResFinder builds upon assembly and BLAS (Clausen et al. 2016), the tigecycline-resistance tet(X) gene was found in the genome of the strain D16T. Strain D16T was susceptible to tetracycline, Which was consistent with that Flavobacteriaceae is a potential ancestral source of tigecycline resistance gene tet (X) (Zhang et al. 2020).
Carbohydrate-Active Enzymes (CAZymes)
Studies have shown that polysaccharides utilization was found in many Bacteroidota bacteria of marine(Francis et al. 2019). Genomes of D16T and closely related genera were analyzed by the CAZy database. Strain D16T, Eudoraea adriatica DSM 19308T, Zeaxanthinibacter enoshimensis TD-ZE3T, Robiginitalea biformata HTCC 2501T, Muriicola jejuensis EM44T, Maribacter vaceletii DSM 25230T and Cellulophaga tyrosinoxydans EM41T contain 160, 129, 109, 112, 87, 231 and 237 carbohydrate-active enzymes, respectively (Figure 4). carbohydrate-active enzymes including 46 GHs were found in strain D16T, which was less than Eudoraea adriatica DSM 19308T, Zeaxanthinibacter enoshimensis TD-ZE3T and Robiginitalea biformata HTCC 2501T, Muriicola jejuensis EM44T and Cellulophaga tyrosinoxydans EM41T. The experimental results show strain D16T, Eudoraea adriatica DSM 19308T, Zeaxanthinibacter enoshimensis TD-ZE3T and Robiginitalea biformata HTCC 2501T, Muriicola jejuensis EM44T and Cellulophaga tyrosinoxydans EM41T can utilize various carbon sources (Table1 ), which was consistent with the genome annotation results.
Table 3. Differential characteristics of strain D16T (Aoguangibacterium sediminis gen. nov., sp. nov.) and members of related genera of the family Flavobacteriaceae.
Strains: 1, D16T(data from this study); 2, Eudoraea adriatica DSM 19308T(Alain et al. 2008); 3, Muriicola jejuensis EM44T(Kahng et al. 2010); 4, Zeaxanthinibacter enoshimensis TD-ZE3T(Asker et al. 2007); 5, Robiginitalea biformata HTCC2501T(JC and SJ 2004); 6, Maribacter arcticus KOPRI 20941T (Cho et al. 2008); 7, Cellulophaga tyrosinoxydans EM41T (Kahng et al. 2009); 8, Arenibacter troitsensis DSM 19835T(Nedashkovskaya et al. 2003); 9, Poritiphilus flavus R33T(Wang et al. 2020).
+, Positive; 2, negative; ND, not determined; W, weak
Characteristic
|
1
|
2
|
3
|
4
|
5
|
6
|
7
|
8
|
9
|
Pigmentation*
|
R/O
|
-
|
Y/O
|
Y
|
Y/O
|
-
|
-
|
-
|
+
|
Flexirubin
|
+
|
-
|
-
|
-
|
-
|
N
|
+
|
-
|
-
|
Gliding motility
|
-
|
-
|
-
|
+
|
-
|
+
|
+
|
-
|
-
|
Reaction to oxygen!
|
A
|
A
|
A
|
A
|
A
|
A
|
A
|
A
|
A
|
Growth at 42 ℃
|
+
|
-
|
-
|
-
|
+
|
-
|
-
|
-
|
-
|
Nitrate reduction
|
-
|
-
|
+
|
-
|
-
|
+
|
+
|
-
|
+
|
Oxidase/catalase production
|
-/+
|
+/+
|
+/+
|
+/+
|
+/+
|
+/+
|
+/+
|
+/+
|
+/+
|
β-galactosidase
|
-
|
+
|
+
|
-
|
-
|
-
|
+
|
+
|
+
|
Glucose acidification
|
-
|
-
|
+
|
+
|
-
|
-
|
+
|
-
|
+
|
Degradation of:
|
|
|
|
|
|
|
|
|
|
Starch
|
+
|
-
|
-
|
+
|
+
|
-
|
+
|
+
|
+
|
Aesculin
|
-
|
+
|
+
|
+
|
+
|
-
|
N
|
|
N
|
Gelatin
|
-
|
-
|
-
|
+
|
-
|
-
|
+
|
+
|
+
|
DNA
|
-
|
N
|
N
|
-
|
-
|
N
|
-
|
N
|
N
|
Cellulose
|
+
|
-
|
-
|
-
|
-
|
N
|
-
|
-
|
N
|
Urea
|
+
|
-
|
-
|
-
|
-
|
-
|
-
|
-
|
+
|
Arginine
|
-
|
-
|
N
|
N
|
-
|
N
|
|
-
|
+
|
16S rRNA gene sequence similarity (%)#
|
100
|
93
|
91.9
|
92.3
|
92.3
|
91.9
|
91.5
|
90.9
|
92.5
|
*Y, Yellow; O, orange; DO, dark orange.
!O, Obligate aerobe; A, facultative anaerobe.
#Similarities calculated in reference to the 16S rRNA gene sequence of strain D16T
Table 4. Phenotypic and genotypic characteristics of strain D16T (Aoguangibacterium sediminis gen. nov., sp. nov.)
Characteristic
|
Strain D16T
|
Temperature range for growth (°C) (optimum)
|
10-42 (30-33)
|
NaCl range for growth (%, w/v) (optimum)
|
0-5 (1)
|
pH range for growth (optimum)
|
5.5-9.5 (7.0)
|
Biochemical properties (API 20E)
|
|
o-nitrophenyl-β-D-galactopyranoside, Urease
|
+
|
arginine dihydrolase, lysine decarboxylase, ornithine decarboxylase
|
-
|
citrate utilization, H2S production, tryptophan deaminase, indole production
|
|
Voges–Proskauer reaction, gelatinase, glucose, mannitol, nositol, sorbitol, rhamnol
|
|
sucrose, melibiose, amygdalin, arabinose
|
|
Enzymic activities (API ZYM)
|
|
alkaline phosphatase, esterase (C4), esterase lipase (C8), leucine arylamidase, valine arylamidase
|
+
|
trypsin, α-chymotrypsin, acid phosphatase, naphthol-AS-BI-phosphohydrolase,
|
|
β-glucosidase, N-acetyl-β-glucosaminidase
|
|
lipase (C14), cystine arylamidase, α-galactosidase, β-galactosidase, β-glucuronidase,
|
-
|
α-glucosidase, α-mannosidase, α-fucosidase
|
|
Substrate assimilation
|
|
1% Sodium Lactate, Fusidic Acid, D-Serine, Troleandomycin, Rifamycin SV ,Minocycline, Lincomycin, Glucuronamide, Nalidixic Acid, Potassium Tellurite, Acetoacetic Acid, Sodium Bromate
|
+
|
Vancomycin, L-Malic Acid, Bromo-Succinic Acid, Lithium Chloride
|
Acetic Acid, Formic Acid, Aztreonam, Sodium Butyrate
|
|
Tween 20, Tween 40, Tween 60, Tween 80
|
-
|
Degration of macromolecules
|
|
Amylase, Casein
|
+
|
Cellulase, Agar, Algin, Casein, Dnase, starch, urea
|
-
|
Susceptibility to antibiotics
|
|
Tetracycline, ampicillin, rifampin, carbenicillin, vancomycin.
|
+
|
Polar lipids
|
DPG、PE、AGL、GL、AL、PGS
|
Quinones
|
MK-6
|
DNA G+C content (mol%)
|
42.8
|
Phenotypic and Biochemical Characteristics
Strain D16T formed aerobic, non-motile, red, mesophilic and slightly halotolerant colonies with diameters of 1-3 mm after growing for 96 hours on modified MA at 33 °C. The cell was Gram-stain-negative, motile by gliding and spherical (0.2-0.4 μm wide and 1.3-2.0 μm long) (Figure 2). Strain D16T could adapt to growth at higher temperatures (Table 3). The strain D16T was able to hydrolysis cellulose, urea and starch but was negative for degradation of arginine, Tweens (80, 60, 40, and 20), DNA, gelatin, aesculin agar, alginate and casein. Strain D16T produced Poly-β-hydroxybutyrate, similar to the related strains. Strain D16T was susceptible to tetracycline, ampicillin, rifampin, carbenicillin, vancomycin. The results are summarized in Table 4. Comparative analysis of physiological and biochemical strain D16T and related strains were summarized in Supplementary Figure 12.
Chemotaxonomic Analysis
Major components of fatty acids (>10 %) of strain D16T were iso-C15:0 ( 32.0 %) and C16:0 (10.5 %) , iso-C15 : 0 was also detected as the major fatty acids in type species of all type genus in Flavobacteriaceae. Listed in Table 5 and iso-C15 : 1 was also the major component in type strains of Eudoraea adriatica DSM 19308T, Muriicola jejuensis EM44T, Robiginitalea biformata HTCC 2501T, Maribacter arcticus KOPRI 20941T and Arenibacter antarcticus R18H21T. The five most closely phylogenetic relatives. However, C16:0, one of the major fatty acids in strain D16T, was detected only in trace amounts in type strains of Muriicola jejuensis EM44T, Robiginitalea biformata HTCC 2501T and Arenibacter antarcticus R18H21T (Table 5). Polar lipids of strain D16T included Diphosphatidylglycerol (DPG)、Phosphatidylethanolamine (PE)、Aminoglycolipid (AGL)、glycolipid (GL)、Aminolipid (AL) and one unidentified lipid PGS (Supplementary Figure 13). The respiratory quinone of strain D16T was menaquinone-6 (MK-6), typical features of the family Flavobacteriaceae (J.-F. Bernardet et al., 2002).
Environmental Distribution
The 16S rRNA sequence of the strain D16T was added to the SILVA SSU database ((Quast et al. 2013); SSU Ref NR release 138.1) and aligned by the ARB version 7.0 software package (Ludwig et al. 2004). The phylogenetic position of strain D16T and relative genera were determined by comparative analyses of these sequences in the ARB version 7.0 software package. Compare 16S rRNA sequence the strain D16T, those sequences of strains whose sequence similarity between D16T was higher 95% was selected from SILVA SSU 138.1. The environmental distribution of these strains and relative genera was analyzed by the isolation source form environmental sequences affiliated with relative genera in the European Nucleotide Archive (ENA)8. Supplementary Table S1 lists all the 16S rRNA gene sequences used.
The difference between the proposed Aoguangibacterium genus and other related genera was also reflected in different environmental isolation source. The majority of representatives of the genus Muriicola, Eudoeaea, Maribacter, Cellulophaga and Zeaxanthinibacter inhabit marine environments, soil or fresh water and are expected to have a metabolism adapted to aerobic or microaerophilic conditions. as shown in Supplementary Figure 14, the proposed Aoguangibacterium genus of the family Flavobacteriaceae mainly was isolated from hypoxic or microhypoxic environments, like marine sediments or hypersaline sediments. Only a few representatives were found in aerobic environments such as sea water or freshwater. Strain D16T, isolated from marine sediments, is similar to isolation source of the proposed Aoguangibacterium genus of the family Flavobacteriaceae, which was consistent with results of anaerobic experiments and denitrification pathway annotated in genome of D16T. Specific analysis of the isolation source of members of the proposed Aoguangibacterium genus and related genera were presented in Supplementary Table S1.
Discussion
The majority of representatives in the family Flavobacteriaceae were distributed in the marine environment and are expected to adapt to micro-aerobic and anaerobic environments. One of the family Flavobacteriaceae was isolated from human blood (Leyer et al. 2020). Candidatus Endobryopsis, one of the family Flavobacteriaceae, can live in symbiosis with the alga Bryopsis sp (Zan et al. 2019). Marine sediments, especially offshore sediments, which was a habitat where posses extremely active material transformation and energy flow. The offshore area, where the water depth is less than 50M no more than 2% of the total ocean area, transported to the ocean almost 48% of the global inflows of organic carbon fluxes into the oceans (Arndt et al. 2013).
In the genome of strain D16T, 1373 genes were predicted, including 1340 protein-coding genes,which include genes associated with CO2 fixation and denitrification, which can be inferred that the strain D16T was also autotrophic. Some differences in key metabolic pathways also supported the separate evolution of the strain D16T and other related genera of the family Flavobacteriaceae. (succinate dehydrogenase/fumarate reductase (sdhC/ frdC) catalyzes the generation of Fumarate and Succinate in the Citrate cycle (TCA cycle, Krebs cycle) and functions as a key switch in the regulation of carbon flux distribution (Akram 2014). Both the substrates and products of sdhC/frdC are involved in the tricarboxylic acid cycle, anaplerosis and energy anabolism. sdhC/frdC was present in the genomes of the proposed genus D16T. However, it was missing in the genome of other relative genera.
Furthermore, it is involved in the second carbon oxidation of the Citrate cycle, which was also only present in D16T. Only strain D16T has a complete tricarboxylic acid (TCA) cycle, which is related to the metabolism of the three major substances and is also the hub of energy metabolism (Steffens et al. 2021). In addition, for the heme biosynthesis pathway, uroporphyrinogen-III synthase (hemD, UROS) was present in D16T, which was missing in other relative genera. The cyclic tetrapyrrole heme is used as a prosthetic group in a wide variety of different proteins in almost all organisms. It is often essential for vital biochemical processes such as aerobic and anaerobic respiration and photosynthesis (Layer 2021), which is consistent with the that the D16T has CAM (Crassulacean acid) metabolism) metabolic pathways. This indicates it can perform photosynthesis (Yurkov and Beatty 1998; Béjà et al. 2002; Karl 2002; Oz et al. 2005; Fuchs et al. 2007; Suzuki and Béjà 2007).
The evolutionary divergence among members the strain D16T and other relative genera was also reflected by different growth performances under aerobic conditions. The strain D16T cannot grow under microaerobic/aerobic conditions. However, other relative genera can also grow under microaerobic/aerobic conditions, relative to the environmental distribution of members of the proposed genus Aoguangibacterium and other relative genera. (JC and SJ 2004; Asker et al. 2007; Alain et al. 2008; Hu et al. 2015; Wang et al. 2020).
As a new genus of the family Flavobacteriaceae, strain D16T has many features that differentiate it from other bacterial groups of Flavobacteriaceae. Such as cellulose degradation and flexirubin biosynthesis, absent in other related genera, was found in the strain D16T. Compared to the high levels of iso-C15:0 and 3-OH i17:0 of most members of Flavobacteriaceae, the strain D16T only has a low proportion of iso-C15:0 and 3-OH i17:0. Most members can degrade lipid (Tweens20 40, 60, and 80). However, the strain D16T can not degrade Tweens 20, 40, 60, and 80. Many family Flavobacteriaceae were positive oxidase production, but oxidase was absent in the strain D16T.
Taxonomic Conclusion
In this study, a novel strain D16T with carbon fixation and denitrification was isolated from marine sediment from China. Phylogenetic trees based on genomes and 16S rRNA genes illustrated that the strain D16T formed a separate branch with related genera (Eudoraea, Poritiphilus, Zeaxanthinibacter, Robiginitalea and Muriicola) and was distinctly different from other genera of family Flavobacteriaceae (Figures 1–2). Compare strain D16T with closest relative genera; the ANI, POCP, DDH and AAI values didn’t exceed the recommended critical point for genus or species division (Table 5), revealing that the strain D16T represents a novel genus. When comparing the strains D16T and close relatives, Their physiological and biochemical characteristics were markedly different. Therefore, strain D16T indicated a new genus in the family Flavobacteriaceae.
Description of Aoguangibacterium sediminis gen. nov., sp. nov.
Description of Aoguangibacterium gen. nov.
Aoguangibacterium (Ao.guang.i.bac.te’ri.um. N.L. masc. n. bacterium a rod; N.L. neut. n. Aoguangibacterium, a rod named after Ao Guang, Dragon King of the Eastern Sea in Chinese folklore).
Gram-stain-negative, non-motile, Aerobic, non-spore-forming, red, mesophilic, denitrifying, carbon-fixing and slightly halotolerant. Pigments are Red or Orange. Oxidase-negative and catalase-positive. Neutrophilic. Polar lipids include (DPG、AGL、 GL 、 AL 、PGS). The predominant quinone is MK-6. Major cellular fatty acids were (iso-C15:0, 32.0 %) and (C16:0, 10.5 %). The G+C content of the DNA is 42.8 mol %. Phylogenetically, the genus Aoguangibacterium belongs to the family Flavobacteriaceae, phylum Bacteroidota, showing a distant relatedness to the marine genera Eudoraea, Zeaxanthinibacter, Robiginitale, Cellulophaga, Muriicola and Poritiphilus. The type species is D16T.
Aoguangibacterium sediminis sp. nov.
Aoguangibacterium sediminis (se.di’mi.nis. L. gen. neut. n. sediminis, of sediment).
Displays the following characteristics in addition to those given in the genus description. Cells are 0.2-0.4 μm in width and 1.3-2.0 μm in length; On MA, colonies (1.0-3.0 mm in diameter) are round, red, smooth, non-motile and shiny. Growth occurs at 10-42 °C (optimum,30-33 °C), at pH 5.5-10.0 (optimum,7.0) and at 0–5% (optimum, 1.0%) NaCl; Starch, cellulose and urea are hydrolyzed, but aesculin, gelatin, DNA and arginine are not. O-nitrophenyl-β-D-galactopyranoside and Urease are produced. arginine dihydrolase, lysine decarboxylase, ornithine decarboxylase, citrate utilization, H2S production, tryptophan deaminase, indole production, Voges–Proskauer reaction, gelatinase, glucose, mannitol, inositol, sorbitol, rhamnol, sucrose, melibiose, amygdalin and arabinose are not produced. Alkaline phosphatase, esterase (C4), esterase lipase (C8), leucine arylamidase, valine arylamidase, trypsin, α-chymotrypsin, acid phosphatase, naphthol-AS-BI-phosphohydrolase, β-glucosidase, N-acetyl-β-glucosaminidase activities are positive. lipase (C14), cystine arylamidase, α-galactosidase, β-galactosidase, β-glucuronidase, α-glucosidase, α-mannosidase, α-fucosidase activities are negative. Acid is produced from 1 % Sodium Lactate, Fusidic Acid, D-Serine, Troleandomycin, Rifamycin SV, Minocycline, Lincomycin, Glucuronamide, Vancomycin, L-Malic Acid, Bromo-Succinic Acid, Nalidixic Acid, Lithium Chloride, Potassium Tellurite, Acetoacetic Acid, Acetic Acid, Formic Acid, Aztreonam, Sodium Butyrate, Sodium Bromate. Predominant cellular fatty acids (representing 10 % of the total fatty acids) are (iso-C15:0, 32.0 %) and (C16:0, 10.5 %). The DNA G+C content of the type strain is 42.8 mol %.
The type strain, D16T (=MCCC 1H00463T= KCTC 82746T), was isolated from coastal sediment from Xiaoshi Island, Weihai, China. The DNA G + C content of the type strain is 42.8%. Accession numbers are MZ695036 (16S rRNA gene) and JAJTJC000000000 (whole genome).