The sediment sample was collected from approximately 10 cm depth of the intertidal zone (37°31′ 36′′ N, 122°00′ 58′′ E) during September 2019. To obtain DNA for microbial community analysis, 5 mL of samples were centrifuged at 12,000 × g for 10 min, the supernatant was discarded, and the pellet was resuspended in 250 μL of sterile Milli-Q water. Then, DNA was extracted using the FastDNASpin Kit for Soil (MP Biomedicals, OH). The primer set used in the study was 338F (5′-ACTCCTACGGGAGGCAGCAG-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′) (Mori et al., 2013), which was designed to amplify the V3–V4 region of the 16S rRNA gene and was demonstrated in silico to be universal for nearly all bacterial taxa. Sequencing was performed on a MiSeq platform at Majorbio Bio-Pharm Technology Co., Ltd. (Shanghai, China).

The strains Oceanipulchritudo coccoides CK1056^T (from our lab), Coraliomargarita sinensis WN38^T (from our lab), Coraliomargarita akajimensis KCTC 12865^T, Pelagicoccus mobilis KCTC 13126^T, and Cerasicoccus arenae KCTC 12870^T that were obtained from the Korean Collection for Type Cultures (KCTC) and Puniceicoccus vermicola JCM 14086^T and Ruficoccus amylovorans JCM 31066^T that were obtained from the Japan Collection of Microorganisms (JCM) were tested for antimicrobial susceptibility. The standard disk diffusion technique was used against 20 different antimicrobial disks, according to the Clinical and Laboratory Standards Institute (CLSI) (CLSI, 2018). Suspensions of the isolates equivalent to 0.5 McFarland standard turbidity were prepared, and MA plates were incubated at the optimal temperature before aseptic placement of antibiotic disks onto the surface. The antimicrobial activity was determined by measuring the diameters of the zones of inhibition formed around each paper disk. The tested antibiotics were categorized as sensitive, intermediate sensitive, or resistant, according to CLSI.

According to the results of antimicrobial susceptibility testing of strains of Puniceicoccaceae, all the studied strains (the related strains) were resistant to ofloxacin and norfloxacin. Therefore, the strains of phylum Verrucomicrobia were cultured under the selective pressure of ofloxacin and norfloxacin. Isolation was performed using marine agar 2216 (MA; BD) supplemented with ofloxacin (100 μg l^–1) and norfloxacin (200 μg l^–1). The sample was diluted with sterilized artificial seawater. Serial dilutions were spread onto MA plates and incubated under aerobic conditions at 30°C for 7 days. After incubation, a white colony was selected from the plate, and the NFK12^T strain was obtained. The discrete colonies formed on the plate were selected and subjected to the identification procedure based on the analysis of the 16S rRNA gene sequence. Following purification with repeated streaking, the strain was routinely cultivated on MA at 30°C for 3 days. The strain was stored in sterile 15% (v/v) glycerol supplemented with 1% (w/v) NaCl suspensions at −80°C. The type strain, Pelagicoccus albus JCM 23202^T, which was obtained from the Japan Collection of Microorganisms (JCM), was used as a related strain for comparative purposes.

DNA was prepared using a bacterial DNA isolation kit (Takara). The 16S rRNA gene was amplified using primers 27F and 1492R (Lane, 1991). Sequencing was performed by the BGI Co. Ltd. (Qingdao, China), which resulted in a 1,406-bp almost complete 16S rRNA gene sequence. The cloned 16S rRNA gene sequence of the strain NFK12^T was preliminarily identified by searching for matches in the EzBioCloud database¹ and NCBI databases². Comparative sequence analyses of the 16S rRNA gene sequences were used to determine the phylogenetic position of the NFK12^T strain and related bacteria.

To definite its exact taxonomic status, we did a further phylogenetic and taxonomic investigation. The 16S rRNA gene sequences of the novel strain NFK12^T and cultured type strains were aligned by using MAFFT version 7 (Katoh and Standley, 2013). The automated removal of spurious sequences and poorly aligned regions from the multiple sequence alignment was performed using trimAl (Capella-Gutierrez et al., 2009). The phylogenetic analysis was constructed by FastTree (Price et al., 2010) using GTR + CAT parameters and IQ-TREE (Jana et al., 2016) using the GTR + F + I + G4 model. Besides, these sequences were aligned by using MUSCLE (Edgar, 2004). A phylogenetic tree was reconstructed by using the neighbor-joining (Saitou and Nei, 1987) algorithm implemented in the software package MEGA (version 7.0) (Kumar et al., 2016). Phylogenetic trees were also generated with the maximum-likelihood (Felsenstein, 1981) and maximum-parsimony (Fitch, 1971) algorithms. Evolutionary distances were calculated using the Kimura 2-parameter model Bootstrap analysis was performed with 1,000 replications to provide confidence estimates for tree topologies.

In addition, 16S rRNA gene sequence of Puniceicoccaceae analysis was performed using the sequences available at the NCBI and EzBioCloud databases with completeness >90%. The 16S rRNA gene sequences of all of the strains of Puniceicoccaceae, including uncultured strains, were aligned using MAFFT version 7 (Katoh and Standley, 2013). The automated removal of spurious sequences and poorly aligned regions from the multiple sequence alignment was performed using trimAl (Capella-Gutierrez et al., 2009). The phylogenetic trees were constructed using IQ-TREE (Jana et al., 2016) by employing the GTR + F + I + G4 parameter and 1,000 bootstrap replications based on the 16S rRNA gene sequences from taxa with validly published names and the uncultured group containing almost 102 sequences. Bootstrap values were evaluated based on 1,000 replicates. Phylogenetic analysis was performed using MEGA version 7.0. The tree was rooted in Desulfocarbo indianensis SCBM^T (KF032904), which was then pruned from the tree. The tree was collapsed and formatted using iTOL v4 (Letunic and Bork, 2019).

The draft genome sequences were determined using the Shanghai OE Biotech Co., Ltd. (Shanghai, China). Sequencing was performed on an Illumina HiSeq Xten platform (Illumina Inc., San Diego, CA, United States). An Illumina shotgun library using the Illumina TruSeq Nano DNA LT Sample Prep Kit was reconstructed and sequenced in paired ends using the Illumina HiSeq platform. Raw sequencing data were generated using the Illumina base-calling software CASAVA v1.8.2³, according to the manufacturer’s protocol. The sequenced reads were assembled using SOAPdenovo software (Li et al., 2009).

To determine its exact taxonomic status, further phylogenetic and taxonomic investigations were performed. Multilocus sequence analysis was performed for all the genomes reported in this study. All genomes were available from the NCBI and EzBioCloud databases. The genomes of the members of Puniceicoccaceae were identified using the Rapid Annotations using Subsystems Technology (RAST) server (Aziz et al., 2008), with a focus on EPS synthesis. Besides, gene content was annotated using the NCBI Prokaryotic Genome Annotation Pipeline. Gene genes involved in metabolic pathways were analyzed in detail by using information from the KEGG database (Kanehisa et al., 2016). The genomes were aligned using MUSCLE v.3.8.31. An up-to-date bacterial core gene set (UBCG) and a phylogenomics pipeline were utilized for phylogenetic tree construction. The phylogenetic analysis was constructed by FastTree (Price et al., 2010) using GTR + CAT parameters and IQ-TREE (Jana et al., 2016) using the GTR + F + I + G4 model. Phylogenetic analysis was performed with 1,000 bootstrap replicates based on 36 genomes. The genome of D. indianensis SCBM^T served as the outgroup. The protein sequences (2,133 bp) of all the genomes reported in this study were subjected to phylogenetic analysis.

As measures of relatedness between the isolate and closely related taxa, the average amino acid identity (AAI) values, percentage of conserved proteins (POCP) values, the average nucleotide identity (ANI) (Yoon et al., 2017) values and the digital DNA-DNA hybridization (dDDH) values were determined. The average amino acid identity (AAI) between genomes were calculated by CompareM⁴ (Rodriguez-R and Konstantinidis, 2014). The POCP was calculated as described by Qin et al. (2014). The average nucleotide identity (ANI) values were calculated using JSpeciesWS⁵. The genome-to-genome distance calculator (GGDC 2.1⁶) was used to calculate all dDDH (Goris et al., 2007) values.

Phenotypic characterizations, biochemical characterization, and chemotaxonomy were performed using cells grown on MA at 33°C for 2 days. Cell morphology and size were examined by light microscopy (E600, Nikon), transmission electron microscopy (JEM-1200, Jeol), and scanning electron microscopy (Nova NanoSEM450, FEI). Physiological and biochemical experimental details have been described in previous studies (Feng et al., 2020). Motility was determined by the hanging-drop method and tested by inoculating the bacteria on the 0.5% agar. The Gram reaction was determined as described by Tindall et al. (2007). The temperature range for growth was determined on MA at 4, 10, 15, 20, 25, 28, 30, 33, 37, and 40. The NaCl concentrations for growth were determined by incubating the bacteria on modified marine broth 2216 made with 0.5% peptone, 0.1% yeast extract, artificial seawater (consisting 0.32% MgSO₄, 0.22% MgCl₂, 0.12% CaCl₂, 0.07% KCl, 0.02% NaHCO₃, w/v) and in the presence of 0.0–10.0% (w/v) NaCl at intervals of 0.5%. The effects of pH were determined by adding the appropriate buffers (Sangon), including 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, and 9.5) to MB with a concentration of 20 mM and the pH of control groups were checked after autoclaved. OD₆₀₀ values of the cultures were measured after incubation for 2 days at 33°C. Unless otherwise mentioned, all experiments were performed with three replicates. Catalase activity was determined by the observation of bubble formation in a 3% H₂O₂ solution. Oxidase activity was examined using an oxidase reagent kit (bioMérieux) according to the manufacturer’s instructions. Tests for the hydrolysis of starch, casein, alginate, carboxymethylcellulose, Tweens 20, 40, 60, and 80 were determined as described by Cowan (Cowan et al., 1966). The API 20E (bioMérieux) kits were performed according to the manufacturer’s instructions, using the biomass of strain NFK12 grown on MA at 33°C for 2 days. Production of other enzymes was assessed using API ZYM (bioMérieux) kits. Carbon source oxidation was checked in BIOLOG GEN III microplates⁷. All the API tests were performed according to the manufacturer’s instructions. All of the API and the BIOLOG tests were performed with two replicates. All tests were carried out simultaneously with the related type strain.

Polar lipids were determined using 2D thin-layer chromatography (TLC) (Minnikin et al., 1984). Isoprenoid quinones of the strain NFK12^T were analyzed using reverse-phase HPLC (Kroppenstedt, 1982). The preparation and extraction of fatty acid methyl esters from biomass and their subsequent separation and identification by gas chromatography were performed as described elsewhere (Athalye et al., 2010). The fatty acids were extracted according to the standard protocol of MIDI (Sherlock Microbial Identification System, version 6.1). The fatty acids were methylated and analyzed by an Agilent 6890N gas chromatograph. Cellular fatty acids were identified using the TSBA40 database of the microbial identification system.

The strains were subjected to Congo red staining to detect the presence of EPS. The strains NFK12^T and P. albus JCM 23202^T were cultured at 33°C for 3 days on Congo red agar plates (supplemented with 0.8 g/L Congo red) with MA (Freeman et al., 1989) or in MB with Congo red solution (0.3 g/L). Based on the results of genomic and phylogenomic analyses of Puniceicoccaceae, other strains of the family were also suspected to produce EPS. The type strains, O. coccoides CK1056^T, C. akajimensis KCTC 12865^T, P. mobilis KCTC 13126^T, C. arenae KCTC 12870^T, P. vermicola JCM 14086^T, and R. amylovorans JCM 31066^T, were inoculated on the Congo red agar plates. The cells were harvested, washed three times, and then resuspended in sterile artificial seawater (OD₆₀₀ = 3.5). Congo red (150 μg/mL) and Trypan blue (20 μg/mL) were mixed evenly in a total volume of 125 μL with the suspensions of the cell strains. The samples were incubated at room temperature (RT) in the dark for 30 min and centrifuged at 12,000 × g for 10 min. Then, 200 μL of the supernatant was added to 96-well plates, and the optical density was measured at 490 and 585 nm for Congo red and Trypan blue, respectively.

The NFK12^T strain was cultured at 33°C with shaking at 200 rpm for 3 days to induce EPS production. The crude EPS from the bacterial isolates was obtained as previously described (Liu et al., 2013). The culture medium was heated in boiling water for 10 min to inactivate the enzymes and then cooled down to 25°C. Then, the medium was centrifuged (20 min, 10,000 g, and 4°C) to remove cells and coagulated proteins, and the supernatant was collected. Next, 80% (w/v) trichloroacetic acid was added to the supernatant to a final concentration of 4% (w/v) with gentle stirring and kept at 4°C for 10 h. Protein precipitates were removed by centrifugation, and then EPS were precipitated from the supernatant with three volumes of cold ethanol followed by overnight incubation at 4°C. After centrifugation, the precipitate was resuspended in deionized water and dialyzed (molecular weight cut-off: 14 kDa) for 3 days. The retentate was centrifuged to remove the insoluble material. The supernatant was lyophilized, and white crude polysaccharides (LCP) were obtained (Ai et al., 2008). Total EPS (expressed as mg l^–1) was estimated by the phenol-sulfuric acid method using glucose as a standard (Dubois et al., 1956), and the amount of EPS produced by each isolate was calculated. The type species of the family Puniceicoccaceae were also cultured for EPS production, using the same method as that for the NFK12^T strain.

Scanning electron microscopy-energy dispersive X-ray (SEM-EDX) (model Nova NanoSEM450; FEI) was employed to analyze the film cross-sections of the crude EPS from the bacterial isolates. Films were placed on a metallic stub and were covered with gold under vacuum in an argon atmosphere.

The NFK12^T strain was cultured at 33°C for 2 days on MA with or without EPS produced by P. enzymogenes (300 μg/mL). Since the process of extracting EPS from the strain NFK12^T was complicated and chitosan as the only cationic marine polysaccharide is a random copolymer obtained from the alkaline deacetylation of chitin (Muxika et al., 2017), we consider whether it can be replaced by commercial chitosan for other related experiments. Therefore, the strain NFK12^T was cultured at 33°C for 2 days on MA with chitosan (Sangon Biotech, 300 μg/mL). Samples were then removed and serially diluted to 10^–6 in sterilized seawater and spotted on the modified MA agar plates for culturing at 33°C. Colony formation unit (CFU) counting was performed at 48 h. R. amylovorans JCM 31066^T and Escherichia coli DH5α were cultured on modified NB agar plates at 25°C and on modified LB agar plates at 37°C, respectively. For the strain NFK12^T, CFU counting was performed at 16 and 109 h using the same method. All experiments were performed in triplicate.

Antimicrobial susceptibility testing revealed that the strain NFK12^T was also resistant to ofloxacin (5 μg) and norfloxacin (10 μg). To investigate the MICs of ofloxacin and norfloxacin, the strain NFK12^T was inoculated into a 250-mL Erlenmeyer flask containing 100 mL of MB and incubated for 37 h at 33°C and 200 rpm until reaching the late logarithmic phase. The cells were harvested, washed three times, and then resuspended in sterile artificial seawater (OD₆₀₀ = 1.0). The sample was inoculated into a 96-well culture plate filled with MB medium. An equal volume of ofloxacin and norfloxacin at different concentrations (ofloxacin: 0–150 μg/mL at intervals of 10.0 μg/mL; norfloxacin: 0–300 μg/mL at intervals of 20.0 μg/mL) was mixed with the cell suspensions of the strain NFK12^T in 96-well culture plates. Samples were incubated at 33°C in the dark for 2 days. Then, 200 μL of bacterial cultures were added to 96-well plates, and growth was measured by recording the optical density (OD) at 600 nm. For R. amylovorans JCM 31066^T, the MIC was determined using the same method.

The effect of EPS on the antibiotic tolerance of the strain NFK12^T was studied by comparing biomass of EPS producing strains against non-EPS producing strains. Samples were then removed and serially diluted to 10^–7 in sterilized seawater and spotted on MA agar plates supplemented with 180 μg/mL norfloxacin or 50 μg/mL ofloxacin with or without 300 μg/mL EPS or 300 μg/mL chitosan at 33°C. Colony formation unit counting was performed at 48 h. In addition, the strain NFK12^T was diluted by 1.0% in MB containing 50 μg/mL with or without 300 μg/mL EPS or 300 μg/mL chitosan and incubated in a shaker (200 rpm) for 0, 6, 12, 24, 36, and 48 h at 30°C. Samples were then removed and serially diluted to 10^–7 in sterilized seawater and spotted on MA agar plates for culturing at 30°C. CFU counting was performed at 48 h. R. amylovorans JCM 31066^T, which had a very low yield of EPS production, were also cultured to verify the effect of EPS on antibiotic tolerance. The method was the same as that of the strain NFK12^T. All experiments were performed in triplicate.

A novel strain exhibiting the defining features of the members of the genus Pelagicoccus was obtained based on its resistance to ofloxacin and norfloxacin. To culture marine Verrucomicrobia strains, sediment samples were collected from Xiaoshi Island (37°31 36′′ N, 122°00 58′′ E), Weihai, China, which is a national marine nature reserve. The samples were subjected to 16S high-throughput analysis, which revealed that the sequences of verrucomicrobia belonged to three families, Puniceicoccaceae, Pedosphaeraceae, and Rubritaleaceae, and two unclassified groups (Supplementary Appendix Table A). Among these, the species of Puniceicoccaceae were the most abundant and the genus Pelagicoccus was detected (Table 1). Furthermore, the antimicrobial susceptibility testing of the type strains of the family Puniceicoccaceae showed that they were resistant to ofloxacin (5 μg) and norfloxacin (10 μg) (Supplementary Appendix Table B).

Due to the results of the resistance to fluoroquinolones based on genomic analysis and antimicrobial susceptibility testing of the type strains of the family Puniceicoccaceae, we attempted to culture the strains of Puniceicoccaceae under the selective pressure of ofloxacin and norfloxacin. The strain NFK12^T was isolated from marine agar plates containing ofloxacin (100 μg l^–1) and norfloxacin (200 μg l^–1). Based on the pairwise comparison of 16S rRNA gene sequences of species available at GenBank databases, the strain NFK12^T exhibited the highest similarity to the most closely validly published species P. albus JCM 23202^T with 98.5% similarity. Based on the 16S rRNA gene sequence similarity threshold value (98.65%) suggested for species delineation (Kim et al., 2014), the strain NFK12^T represents a novel species within the Pelagicoccus clade. An evolutionary tree was constructed using the neighbor-joining method based on the 16S rRNA gene sequences. According to the results, the strain NFK12^T formed an independent linage within the genus Pelagicoccus (Figure 1 and Supplementary Figure 1). The distinct phylogenomic lineage of strain NFK12^T was also supported by the POCP values and average amino acid identity (AAI) values (Table 2). AAI was higher than the threshold proposed to include an organism in a given genus (60.0%). The POCP values were higher than 50% which the 50.0–60.0% POCP value proposed as a genus boundary. The ANI value of the strain NFK12^T with the related strains were lower than 74.6%. In addition, the dDDH value (identities/HSP length) of the strain NFK12^T with its related strains were lower than 20.3%. According to the proposed and generally accepted species boundary for dDDH values (i.e., <70% and ANI < 95%), the strain NFK12^T was identified as a new species distinguishable from the related strain.

After 2 days, cells of the strain NFK12^T exhibited a coccus shape with a diameter of 0.5–1.5 μm (Figures 2A–C). Through transmission electron microscopy (TEM), flagella were observed, but the number of flagella per cell of strain NFK12^T is unclear (Figure 2A). A faint layer of extracellular slime was seen around the cells (Figure 2A), and cells of the strain NFK12^T had randomly distributed finger-like extensions from which the fimbriae appendages projected, measuring more than 500 nm in diameter, as viewed under a scanning electron microscope (SEM) (Figures 2B,C). The electron micrographs indicate that the strain NFK12^T produces peripheral EPS, which has a three-dimensional web-like morphology (Figure 2C), as described in Lentisphaera araneosa (Cho et al., 2004).

The strains were subjected to Congo red staining to detect the presence of EPS. A positive result was indicated by red-stained colonies when the bacterial communities reached the late exponential stage of growth. The color of the strain NFK12^T colonies was checked after 3 days, where the red-stained colonies indicated EPS positivity. The positive colonies were identified by the appearance of a red halo around the colonies (Supplementary Figure 2A). In addition, the red-stained cells gathered at the bottom, while the samples cultured in MB were centrifuged at 3,700 rpm for 10 min (10 min, 3,700 g, RT) (Supplementary Figure 2B). The results indicated that the strain NFK12^T produced EPS, and we speculated that EPS production was one of the common features of Puniceicoccaceae. On the other hand, the results of Congo red and Trypan blue staining showed that the NFK12^T strain and its related strains were EPS-producing isolates (Supplementary Figure 3). The EPS yield of the strain NFK12^T under our experimental conditions was higher than that of its related strains. Among them, the yield of EPS produced by R. amylovorans JCM 31066^T was meager (Figures 3A,B).

The yield of EPS production was characterized by gravimetric analysis of the lyophilized crude EPS. The strain NFK12^T could produce EPS at a yield of 65.5 mg/L by the phenol-sulfuric acid method in MA at an optimal temperature of 33°C. Regarding SEM analysis, film cross-sections of the strain NFK12^T are shown in Supplementary Figure 4. A compact structure was observed in all films, but the surface was not flat. Sulfated polysaccharides are prevalent in marine environments (William, 2017), and the EPS of NFK12^T are polyanionic for the presence of phosphate and sulfate based on the EDX analysis. The yields of EPS synthesized by all strains were the same. The detailed differences in the EPS profiles between the strain NFK12^T and its related strain are shown in Table 3. Among them, EPS production by R. amylovorans JCM 31066^T was not detected by the phenol-sulfuric acid method.

The strain NFK12^T, which produced the highest EPS, was selected for further analysis. For E. coli DH5α, only the colonies on the modified LB agar plates with chitosan were visible after 8 h, and more colonies on all plates gradually emerged over subsequent hours. However, the number of colonies was not significantly different on the plates with or without EPS/chitosan (Figure 4A). For R. amylovorans JCM 31066^T and the strain NFK12^T, the number of CFUs on the modified agar plates with chitosan or EPS was higher than that of the control. In particular, EPS and chitosan significantly promoted the growth of the NFK12^T strain. Therefore, the uncultured strains of the genus Pelagicoccus and other rare bacteria could be supplemented by chitosan or EPS for optimum growth.

In agreement with previous findings (Venkatesan et al., 2015), the growth of the strain NFK12^T with EPS in the presence of ofloxacin or norfloxacin was higher than that of the control. Because of antibiotic treatment, the strain NFK12^T could be observed at sublethal concentrations of 180 μg/mL norfloxacin or 50 μg/mL ofloxacin. After treatment with norfloxacin, the number of CFUs on the plates with chitosan or EPS increased, but it was less than that on the plates without norfloxacin. Since the EPS yield of R. amylovorans JCM 31066^T was meager, we speculated that the effect of chitosan or EPS on the strain might be different after antibiotic treatment. Therefore, we performed the same experiment after antibiotic treatment. The strain R. amylovorans JCM 31066^T could be observed at sublethal concentrations of 260 μg/mL norfloxacin or 110 μg/mL ofloxacin. The result of R. amylovorans JCM 31066^T was the same as that of the strain NFK12^T (Figure 4B). It was concluded that EPS could promote bacterial resistance to norfloxacin. However, we observed that cells did not survive on MA with 50 μg/mL ofloxacin when the CFUs were counted at 48 h. Moreover, some cells resumed growth (T_regrowth = 132 h) after ofloxacin did not work. Therefore, the strain NFK12^T was cultured in MB. After a bacterial population was treated with 50 μg/mL ofloxacin for 0, 6, 12, 24, 36, and 48 h, the surviving cells were cultured on fresh MA plates. With chitosan or EPS in MB, the colonies of the strain NFK12^T were visible after 24 h on MA, while the colonies of the strain NFK12^T were visible after 48 h in the case of control. More colonies gradually emerged over subsequent hours (Figure 4C). After the addition of chitosan or EPS, the number of CFUs of the NFK12^T strain was higher than that of the control. After incubation in MB with ofloxacin for a different time, the number of CFUs of the NFK12^T strain was significantly less than that of the control without ofloxacin treatment. This indicated that some cells were killed by ofloxacin. When chitosan or EPS was added, the number of CFUs increased. However, after incubation with ofloxacin for 36 and 48 h, no effect was observed with the addition of EPS. The number of CFUs was almost the same with or without EPS at 36 and 48 h. In other words, the EPS could tentatively improve the ofloxacin tolerance of the strain NFK12^T, but ofloxacin still worked.

The order Puniceicoccales, containing Puniceicoccaceae, and the order Opitutales, containing the family Opitutaceae, were proposed for the classification of species belonging to subdivision 4 (Cho et al., 2004). The genus Puniceicoccus (Choo et al., 2007) is the type genus of the family Puniceicoccaceae. Besides Puniceicoccus, the genera Cerasicoccus (Yoon et al., 2007a), Ruficoccus (Lin et al., 2017), Pelagicoccus (Yoon et al., 2007b), Coraliomargarita (Yoon et al., 2007c), and Oceanipulchritudo (Feng et al., 2020) have been described as members of the present family. In 2007, the genus Pelagicoccus was proposed by Yoon et al. (2007b) as a new member of the family Puniceicoccaceae. The strains of the genus Pelagicoccus have been isolated from the sea environment. As of December 2020, this genus contains four species: P. mobilis, P. albus, Pelagicoccus litoralis (Yoon et al., 2007b), and Pelagicoccus croceus (Yoon et al., 2007b). A few representatives of Puniceicoccaceae have been cultivated, most of which are yet to be cultured.

With the golden era of the sequencing technology revolution, the abundance, diversity, and ecology of microorganisms gained another dimension. A total of 102 16S rRNA gene sequences of Puniceicoccaceae were available from the NCBI and EzBioCloud databases with completeness >90%, and 80.4% of them belonged to uncultured strains. Cluster analysis of the 16S rRNA gene sequences, including the complete 16S rRNA gene sequence of the strain NFK12^T, was performed, and the phylogenetic trees were constructed using the IQ-TREE server (Price et al., 2010). Of the gene sequences, 67.6% belonged to strains isolated from the marine environment, which could be divided into six groups (groups I–VI) (Supplementary Figure 5). This result was consistent with the previous study that subdivisions 1 and 4 generally dominated marine bacterial communities (Freitas et al., 2012). Most of the cultured species were grouped into groups V and VI. The strains of each group were from different sampling sites and environmental conditions. For group III, there were no cultured type species. Members of Puniceicoccaceae appear especially significant, but isolates do not represent the known phylogenetic breadth from culture-independent studies. Therefore, there seems to be a vast wealth of verrucomicrobial species remaining to be cultured and described.

The draft genome size of the strain NFK12^T was determined to be 6,276,818 bp, arranged into 200 contigs. The G + C content was calculated to be 56.8 mol%, which was in the range of that of the genus Pelagicoccus (51.0–57.0 mol%). An average of 235× coverage depth was accomplished. According to the KEGG and RAST analyses, the membrane fusion protein of RND family multidrug efflux pump, multiple antibiotic resistance protein (marC), and multidrug resistance protein (MATE family) were found in the draft genome of strain the NFK12^T.

In this study, we found that the majority of the new strains brought into the pure culture are affiliated with subdivision 1. In particular, for Puniceicoccaceae, thus far, only two genomes from validly described species (C. sinensis and C. akajimensis) are available. Therefore, we obtained the type species of Puniceicoccaceae and sequenced the draft genome sequences. Nineteen genomes of the uncultured members of Puniceicoccaceae were available at the NCBI and EzBioCloud databases, and the 15 genomes were identified as close neighbors after calculation of AAI with the published genome sequences of the members of Puniceicoccaceae (AAI > 55.0%) (Rodriguez-R and Konstantinidis, 2014). The 15 genomes were utilized for phylogenetic tree construction and were consistently clustered together within the same clade. When the whole genome sequences of these strains were analyzed, nearly all the branches (including the deeper branches) in the protein sequence tree were well supported by bootstrap replicates (Figure 5). In addition, the genomic analysis revealed EPS biosynthesis in those strains, which was almost consistent with the phenotype. Exopolysaccharide biosynthesis protein was predicted in 27 genomes of the class Opitutae after annotation using the RAST server and the NCBI Prokaryotic Genome Annotation Pipeline (Supplementary Appendix Table D). Interestingly, the results indicated that P. albus JCM 23202^T consisted of EPS-producing isolates, but EPS biosynthesis protein was not predicted in the genome. However, it was less clear for us to arrive at that conclusion based on the draft genome. According to the phylogenetic trees of EPS biosynthesis protein sequences of the class Opitutae, the strains of Puniceicoccaceae and Opitutaceae clustered, respectively (Supplementary Figure 6).

Based on the results of genomic analysis, we suspect that almost all bacteria of this class produce EPS. The EPS-producing strains were screened out quickly at the first stage by adding Congo red to the medium. Therefore, besides supplementing with ofloxacin (100 μg l^–1) and norfloxacin (200 μg l^–1), Congo red solution (0.3 g/L) was added to the culture medium. By optimizing the isolation methodology, a large number of Puniceicoccaceae species could be obtained. In this study, we obtained two strains of the genus Coraliomargarita (R2-22 MW255631 and E2-5 MW255632). In future studies, we will use this method to isolate more species of this family and try to verify this method.

Based on 16S rRNA gene sequences analysis, phenotypic, genomic and chemotaxonomic, we suggest that strain NFK12^T could be affiliated to genus Pelagicoccus currently. Consequently, we concluded that strain NFK12^T could be placed in a novel species of the genus Pelagicoccus, for which the name P. enzymogenes sp. nov. is proposed. Other phenotypic characteristics from classical experiments are listed in Table 4 and the species description, while the comparative analyses with members of phylogenetically related strain are given in the species description. The detailed differences in the polar lipid profiles between strain NFK12^T and the related strain are shown in Supplementary Figure 7. The fatty acid compositions of strain NFK12^T and the related strains were similar, but there were differences with strain NFK12^T in the proportions of some fatty acids, particularly summed feature 6 (C_19:1 ω11c and/or C_19:1 ω9c) (Supplementary Table 1).

Pelagicoccus enzymogenes (en.zy.mo’ge.nes. N.L. n. enzyma, enzyme; from Gr. mac. n. zumê, leaven; N.L. suff. -genes, producing; from Gr. v. gennaô, to produce; N.L.mac. adj. enzymogenes, enzyme-producing).

Cells are Gram-stain-negative, aerobic, cocci-shaped. Cells are motile by means of flagella, the mean cell size is 1.0 μm in diameter. Colonies on MA are round, white and glossy with entire margins, about 2.0 mm in diameter after incubation for 3 days at 33°C. Growth occurs at 10–37° (optimum, 33°), at pH 6.0–8.5 (optimum, 7.0–7.5) and in the presence of 1.0–8.5% (w/v) NaCl (optimum, 1.0–2.0%). Growth is not observed under anaerobic conditions on MA with or without 0.1% (w/v) KNO₃. Nitrate is not reduced to nitrite. Catalase-positive and oxidase-positive. The strain was positive for hydrolyzes of Tweens 20, 40, and 60, alginate, starch, casein and cellulose but negative for DNase. Positive for ornithine decarboxylase, urease, gelatinase and citrate utilization. The following carbon sources are oxidized: dextrin, D-mannose, α-D-lactose, D-fucose, β-methyl-D-glucoside, L-rhamnose, glucuronamide, acetic acid and L-aspartic acid. Acid is produced from L-arabinose, D-ribose, methyl β-D-xylopyranoside, D-galactose, D-xylose, glycogen, D-raffinose, D-melibiose, D-lactose, D-maltose, D-cellobiose, D-tagatose and 5-ketogluconate. Positive for esterase (C4), gelatinase, tryptophan deaminase, α-galactosidase, β-galactosidase, β-glucuronidase, alkaline phosphatase, N-acetyl-β-glucosaminidase, acid phosphatase, naphthol-AS-BI-phosphohydrolase. The predominant menaquinone is MK-7 and the major cellular polar lipids are aminophosphoglycolipid, phosphatidylethanolamine and one unidentified lipid. Main cellular fatty acids are anteiso-C_{15: 0}, C_16:0 and summed feature 3 (C_16:1 ω7c and/or C_16:1 ω6c).

The type strain is NFK12^T (=KCTC 72940^T = MCCC 1H00424^T), isolated from sediment of Xiaoshi Island, Weihai, China. The genomic DNA G + C content of the type strain is 56.8 mol%. The GenBank accession number for the 16S rRNA gene sequence of P. enzymogenes NFK12^T is MW023105 and the draft genome has been deposited in GenBank under the accession number JACYFG000000000.