Phenotypic and Genomic Properties of Brachybacterium vulturis sp. nov. and Brachybacterium avium sp. nov.

Two strains, VM2412^T and VR2415^T, were isolated from the feces of an Andean condor (Vultur gryphus) living in Seoul Grand Park, Gyeonggi-do, South Korea. Cells of both strains were observed to be Gram-stain positive, non-motile, aerobic, catalase positive and oxidase negative. Growth was found to occur at 10-30^∘C, showing optimum growth at 30^∘C. The strains could tolerate up to 15% (w/v) NaCl concentration and grow at pH 6-9. The strains shared 99.3% 16S rRNA gene sequence similarity to each other but were identified as two distinct species based on 89.0-89.2% ANIb, 90.3% ANIm, 89.7% OrthoANI and 38.0% dDDH values calculated using whole genome sequences. Among species with validly published names, Brachybacterium ginsengisoli DCY80^T shared high 16S rRNA gene sequence similarities with strains VM2412^T (98.7%) and VR2415^T (98.4%) and close genetic relatedness with strains VM2412^T (83.3–83.5% ANIb, 87.0% ANIm, 84.3% OrthoANI and 27.8% dDDH) and VR2415^T (82.8–83.2% ANIb, 86.7% ANIm, 83.9% OrthoANI and 27.2% dDDH). The major fatty acid of the two strains was identified as anteiso-C_15:0 and the polar lipids consisted of phosphatidylglycerol, diphosphatidylglycerol, presumptively phosphatidylethanolamine and three unidentified glycolipids. Strain VR2415^T also produced an unidentified phospholipid. The cell walls of the two strains contained meso-diaminopimelic acid as diagnostic diamino acid and the whole cell sugars were ribose, glucose, and galactose. The strains contained MK-7 as their predominant menaquinone. The genomes of strains VM2412^T, VR2415^T, and B. ginsengisoli DCY80^T were sequenced in this study. The genomic G+C contents of strains VM2412^T and VR2415^T were determined to be 70.8 and 70.4 mol%, respectively. A genome-based phylogenetic tree constructed using an up-to-date bacterial core gene set (UBCG) showed that the strains formed a clade with members of the genus Brachybacterium, supporting their taxonomic classification into the genus Brachybacterium. Based on phenotypic and genotypic analyses in this study, strains VM2412^T and VR2415^T are considered to represent two novel species of the genus Brachybacterium and the names Brachybacterium vulturis sp. nov. and Brachybacterium avium sp. nov. are proposed for strains VM2412^T (=KCTC 39996^T = JCM 32142^T) and VR2415^T (=KCTC 39997^T = JCM 32143^T), respectively.

Introduction

The genus Brachybacterium was first described by Collins et al. (1988). At the time of writing, the genus contains 19 species with validly published names: B. faecium (Collins et al., 1988), B. nesterenkovii (Gvozdyak et al., 1992), B. conglomeratum, B. paraconglomeratum, B. rhamnosum (Takeuchi et al., 1995), B. alimentarium, B. tyrofermentans (Schubert et al., 1996), B. fresconis, B. sacelli (Heyrman et al., 2002), B. muris (Buczolits et al., 2003), B. phenoliresistens (Chou et al., 2007), B. zhongshanense (Zhang et al., 2007), B. saurashtrense (Gontia et al., 2011), B. squillarum (Park et al., 2011), B. ginsengisoli (Hoang et al., 2014), B. huguangmaarense (Liu et al., 2014), B. hainanense (Liu et al., 2015), B. aquaticum (Kaur et al., 2016) and B. horti (Singh et al., 2016). The genus Brachybacterium is affiliated to the family Dermabacteraceae with the three genera Dermabacter, Devriesea, and Helcobacillus. Currently, the genus Dermabacter is composed of three species: Dermabacter hominis (Jones and Collins, 1988), D. vaginalis (Chang et al., 2016), and D. jinjuensis (Park et al., 2016), while the genera Devriesea and Helcobacillus consist of a single species, Devriesea agamarum (Martel et al., 2008) and Helcobacillus massiliensis (Renvoise et al., 2009), respectively.

Strains of the genus Dermabacter have been classified as Centers of Disease Control and Prevention (CDC) group 3 and group 5 coryneform bacteria (Funke et al., 1994) and have been isolated from diverse clinical specimens such as cerebral abscess (Bavbek et al., 1998), blood cultures (Gómez-Garcés et al., 2001), peritoneal fluid (Radtke et al., 2001), bone deposits from an osteomyelitis patient (Van Bosterhaut et al., 2002) and fatal septicemia (Lee et al., 2011). Therefore, they are considered opportunistic pathogens. Strains of the genus Devriesea are involved in the pathogenesis of dermatitis, cheilitis and septicemia in lizards as facultative pathogenic bacteria (Hellebuyck et al., 2009) and strains of the genus Helcobacillus were isolated from human cutaneous discharge samples, showing the potential for toxicity (Renvoise et al., 2009). However, members within the genus Brachybacterium have not been reported to be pathogenic before now. They have been isolated from environmental samples such as coastal sand, sediment, plants, soil and seawater and from foods including milk products, cheese and salt-fermented seafood. The type strain of B. muris was isolated from murine liver (Buczolits et al., 2003).

Isolation of strains of Brachybacterium species from wide-ranging sources demonstrates that they have adapted to diverse environmental conditions using characteristics obtained in the process of diverging from other genera. Genotypic analysis based on whole genome sequences can indicate which genes are involved in the divergence and can contribute to more accurate taxonomic classification. In present study, two strains were isolated and found to belong to the genus Brachybacterium based on 16S rRNA gene sequences. The whole genome sequences were obtained for the newly isolated strains and from B. ginsengisoli DCY80^T as a closely related species. Although whole genome sequences of some type strains of Brachybacterium species are available in GenBank/EMBL/DDBJ databases, the genome sequences have not been used previously for taxonomic analysis or subjected to comparative genomic studies. Therefore, we used the available genome sequences of Brachybacterium species along with those of Dermabacter and Devriesea species to investigate the taxonomic status of the newly isolated strains within the genus Brachybacterium and to identify intra- and inter-genus differences based on genotypic characteristics.

Materials and Methods

Bacterial Strain Isolation

Strains VM2412^T and VR2415^T were isolated from fecal sample of an Andean condor living in Seoul Grand Park, Gyeonggi-do, South Korea (37°25^′39.7^′′N 127°01^′01.2^′′E). Strains VM2412^T and VR2415^T were isolated using a marine agar (MA) plate and a Reasoner’s 2A (R2A) medium, respectively, which were inoculated with 10^-4 diluted fecal sample and incubated at 20°C. Tryptic soy broth (TSB) and tryptic soy agar (TSA) were used for culture at 30°C.

16S rRNA Gene Sequence Analysis

Genome-derived 16S rRNA gene sequences of strains VM2412^T and VR2415^T were 1,513 bp. The sequences were compared with those of other type strains using EzBioCloud server (Yoon et al., 2017). The 16S rRNA gene sequences of the two strains and type strains of closely related taxa were aligned using BioEdit (Hall, 2011). Phylogenetic trees were constructed by neighbor-joining (NJ) (Saitou and Nei, 1987), maximum likelihood (ML) (Felsenstein, 1981) and maximum parsimony (MP) (Kluge and Farris, 1969) algorithms in the MEGA6 program (Tamura et al., 2013).

Phenotypic Characterization

Based on the 16S rRNA gene sequences, the strains were observed to share high sequence similarities (>98%) with B. ginsengisoli DCY80^T and B. faecium DSM 4810^T. Therefore, these type strains were obtained from the Korean Collection for Type Cultures (B. ginsengisoli KCTC 29226^T) and the Japan Collection of Microorganisms (B. faecium JCM 11609^T) and used for comparison of phenotypic characteristics with the isolated strains.

Growth temperature range was determined in TSB by incubating at 4, 10, 15, 20, 25, 30, 37, 45, 55, and 70°C. The effect of salt concentration on bacterial growth was determined in TSB with 0, 1, 2, 3, 4, 5, 6, 8, 10, 12, 15, 20, 22, and 25% (w/v) NaCl. Growth in TSB at pH 4, 5, 6, 7, 8, 9, 10, and 11 was investigated for identifying pH range and optimal pH for bacterial growth. Bacterial growth was evaluated by measuring the absorbance at 600 nm (OD₆₀₀) after incubation for 24 h, 48 h, and 1 week in the different temperature, salinity and pH conditions.

The two isolates were cultured on TSA [4% (w/v) NaCl, pH 7] and incubated at 30°C for optimal growth. After 48 h, cell morphology was identified using a light microscope (ECLIPSE 50i; Nikon) and colony morphology was determined macroscopically. Gram-staining was conducted with a Gram-staining kit (bioMérieux) according to the manufacturer’s instructions. Cell motility was identified by bacterial migration from an inoculation line in semi-solid agar with 0.4% (w/v) agar (Tittsler and Sandholzer, 1936). Bacterial survival under anaerobic conditions was determined by checking bacterial growth after incubation for 1 week in an anaerobic chamber with 90% N₂, 5% CO₂ and 5% H₂.

Catalase activity was identified by bubble formation after adding drops of 3% (v/v) hydrogen peroxide solution and oxidase activity was confirmed by color transition to deep purple or blue after adding drops of 1% (w/v) tetramethyl-p-phenylenediamine solution (bioMérieux). Bacterial enzyme activities were examined using API ZYM and API 20 NE test strips (bioMérieux). API 20 NE test strips were also used to check the ability of bacteria to assimilate carbohydrates. Bacterial ability to ferment various carbohydrates was tested with API 50 CH test strips (bioMérieux) using API 50 CHB/E media.

For analyses of fatty acids, cell wall components, polar lipids and menaquinone, culture was performed on TSA [4% (w/v) NaCl, pH 7] at 30°C for 48 h. Bacterial fatty acid methyl esters were extracted by saponification and methylation of bacterial cells according to the protocol of the Sherlock Microbial Identification Systems (MIDI, 1999). The composition of fatty acid methyl esters was determined by gas chromatography using a calibration standard (Microbial ID). For analysis of cell wall amino acids, peptidoglycan was purified by the method of Schumann (2011). The purified peptidoglycan was hydrolyzed using 4M HCl and the hydrolysate was spotted onto a cellulose thin layer chromatography (TLC) plate with alanine, lysine, glutamic acid, glycine, aspartic acid, and meso-diaminopimelic acid as reference amino acids. One-dimensional TLC was performed using methanol-water-6M HCl-pyridine (80:26:4:10, v/v) as solvent. The plate was air-dried and sprayed with ninhydrin reagent. Spots were visualized by heating at 100°C for 5 min. For analysis of cell wall sugars, bacterial cells were hydrolyzed according to the protocol of Cummins and Harris (1956). The whole cell hydrolysate was spotted onto a cellulose TLC plate with ribose, arabinose, glucose, galactose, rhamnose, xylose and mannose for spot comparison. The solvent used for one-dimensional TLC development was n-butanol-water-pyridine-toluene (10:6:6:1, v/v). After drying, the plate was sprayed with aniline phthalate reagent and heated at 120°C for 5 min for spot detection. Polar lipids were extracted from bacterial cells and dissolved in chloroform-methanol (2:1, v/v) solution (Minnikin et al., 1984). The composition of polar lipids was determined by two-dimensional TLC using chloroform-methanol-water (65:24:4, v/v) for one dimension and chloroform-acetic acid-methanol-water (80:15:12:4, v/v) for second dimension. Total lipids were detected by spraying with 10% (v/v) ethanolic molybdophosphoric acid, then heating at 180°C for 15 min. Ninhydrin and α-naphthol reagents were used for detection of aminolipids and glycolipids, respectively. Spots were visualized by heating at 110°C for 15 min. Bacterial menaquinone was extracted from bacterial cells using chloroform-methanol (2:1, v/v) solution (Collins and Jones, 1981). The composition of bacterial menaquinones was determined by high-performance liquid chromatography (HPLC) using methanol-isopropanol (2:1, v/v) for elution.

Whole Genome Sequence Analysis

Genomic DNA was extracted using a MG Genomic DNA Purification kit (MGmed), checked by electrophoresis in 1% agarose gel and quantified with a NanoDrop spectrophotometer (Thermo Scientific) and a Qubit fluorometer (Thermo Scientific). To generate a 20 kb library, genomic DNA was sheared with g-TUBE (Covaris) and purified using AMPure PB magnetic beads (Beckman Coulter). A total library was prepared using a DNA Template Prep Kit 1.0 (Pacific Biosciences). DNA/Polymerase Binding Kit P6 (Pacific Biosciences) was used for binding DNA template libraries to DNA polymerase P6. Whole genome sequencing was performed using the PacBio RS II sequencing platform. DNA Sequencing Kit 4.0 (Pacific Biosciences) and PacBio RS II single-molecule real-time (SMRT) Cells 8 Pac (Pacific Biosciences) were used for SMRT sequencing on the PacBio system. Genomes were assembled by RS HGAP Assembly v3.0 and polished with Quiver. Complete genome sequences were annotated using the Rapid Annotation using Subsystem Technology (RAST).

Genomic relatedness was evaluated by different algorithms for genome-to-genome comparison. Average nucleotide identity (ANI) based on BLAST (ANIb) and MUMmer (ANIm) values were calculated in the JSpeciesWS online server (Richter et al., 2016). Digital DNA–DNA hybridization (dDDH) values were computed with Genome-to-Genome Distance Calculator (GGDC) version 2.1 (Meier-Kolthoff et al., 2013). OrthoANI values were determined using the Orthologous Average Nucleotide Identity Tool (OAT) (Lee et al., 2016). Genome-based phylogenetic tree was constructed using an up-to-date bacterial core gene set (UBCG) consisting of 92 genes (Na et al., 2018).

Results and Discussion

16S rRNA Phylogeny

Genome-derived 16S rRNA gene sequences of strains VM2412^T and VR2415^T were found to share high sequence similarities with those of B. ginsengisoli DCY80^T (98.7 and 98.4%, respectively) and B. faecium DSM 4810^T (98.3% for both strains). The two strains formed a monophyletic clade with B. ginsengisoli DCY80^T in the phylogenetic tree based on 16S rRNA gene sequences (Figure 1). The clade containing type strains of Brachybacterium species was distinct from those composed of the members of other genera within the family Dermabacteraceae. Furthermore, the four genera in the family Dermabacteraceae formed a clade distinct from representative members of the family Micrococcaceae. These results suggest that the 16S rRNA tree accurately reflects the phylogenetic relationships.

FIGURE 1

FIGURE 1. Neighbor-joining tree based on the 16S rRNA gene sequences. Timonella senegalensis JC301^T was used as an out-group. Nodes overlapping in neighbor-joining, maximum likelihood and maximum parsimony trees are indicated by filled diamonds. Bootstrap values (NJ/ML/MP) >70% based on 1000 replications are presented at nodes. Bar, 0.005 substitutions per nucleotide.

Phenotypic Characteristics

The two strains were found to grow at 10-30°C, with optimum growth at 30°C. They could tolerate salinity up to 15% (w/v) and grow at pH 6-9. Their cells were observed to be coccoid, Gram-stain positive, non-motile and aerobic. Colonies were observed to be 0.5-1 mm in diameter, ivory-colored, circular and low-convex with entire margins after incubation for 48 h on TSA.

Strains VM2412^T and VR2415^T were found to be catalase positive and oxidase negative. Further biochemical characteristics of strains VM2412^T, VR2415^T, B. ginsengisoli KCTC 29226^T and B. faecium JCM 11609^T were investigated using API kits.

Strain VM2412^T had the following enzyme activities: tryptophanase, urease, protease, esterase (C4), esterase lipase (C8), leucine arylamidase, valine arylamidase, crystine arylamidase, acid phosphatase, naphthol-AS-BI-phosphohydrolase, β-galactosidase, β-glucuronidase, α-glucosidase, β-glucosidase, N-acetyl-β-D-glucosaminidase and α-mannosidase. The strain could ferment glycerol, D-arabinose, D-ribose, D-xylose, D-galactose, D-glucose, D-fructose, D-mannose, dulcitol, D-mannitol, D-sorbitol, N-acetylglucosamine, amygdalin, arbutin, esculin, salicin, D-cellobiose, D-maltose, D-lactose, D-melibiose, D-saccharose, D-melezitose, D-raffinose, amidon, xylitol, gentiobiose, D-turanose, D-lyxose, D-tagatose, L-fucose, potassium gluconate and potassium 5-ketogluconate.

Strain VR2415^T had tryptophanase, protease, esterase (C4), esterase lipase (C8), leucine arylamidase, valine arylamidase, crystine arylamidase, acid phosphatase, naphthol-AS-BI-phosphohydrolase, β-galactosidase, α-glucosidase, β-glucosidase, N-acetyl-β-D-glucosaminidase and α-mannosidase enzyme activities. The strain could ferment glycerol, erythritol, D-arabinose, L-arabinose, D-ribose, D-xylose, D-adonitol, methyl β-D-xylopyranoside, D-galactose, D-glucose, D-fructose, D-mannose, D-mannitol, D-sorbitol, N-acetylglucosamine, amygdalin, arbutin, esculin, salicin, D-cellobiose, D-maltose, D-lactose, D-melibiose, D-saccharose, D-trehalose, D-melezitose, D-raffinose, amidon, glycogen, xylitol, gentiobiose, D-turanose, D-lyxose, D-tagatose, L-fucose, D-arabitol, L-arabitol, potassium gluconate, and potassium 5-ketogluconate.

Four strains used for the biochemical analysis were found to have tryptophanase, esterase (C4), esterase lipase (C8), leucine arylamidase, valine arylamidase, acid phosphatase, naphthol-AS-BI-phosphohydrolase, β-galactosidase, α-glucosidase, β-glucosidase, N-acetyl-β-D-glucosaminidase and α-mannosidase enzyme activities. All these strains could assimilate mannose, mannitol and potassium gluconate and ferment glycerol, D-arabinose, D-ribose, D-xylose, D-galactose, D-glucose, D-fructose, D-mannose, D-mannitol, D-sorbitol, N-acetylglucosamine, amygdalin, arbutin, esculin, salicin, D-cellobiose, D-maltose, D-lactose, D-melibiose, D-saccharose, D-melezitose, D-raffinose, amidon, xylitol, gentiobiose, D-turanose, D-lyxose, L-fucose, potassium gluconate and potassium 5-ketogluconate. Biochemical characteristics that differ between strains are shown in Table 1.

TABLE 1

TABLE 1. Biochemical characteristics of Brachybacterium vulturis VM2412^T (1), B. avium VR2415^T (2), B. ginsengisoli KCTC 29226^T (3), and B. faecium JCM 11609^T (4).

Strain VM2412^T was found to have anteiso-C_15:0 (48.9%), anteiso-C_17:0 (22.3%) and iso-C_16:0 (13.5%) as major fatty acids (>10%). By contrast, strain VR2415^T was found to have anteiso-C_15:0 (57.7%) and iso-C_16:0 (13.7%) as predominant fatty acids (>10%). The fatty acid compositions of strains VM2412^T, VR2415^T and closely related taxa are given in Table 2. The major fatty acid of all strains was identified as anteiso-C_15:0, accounting for 42.9-57.7% of the total composition.

TABLE 2

TABLE 2. Fatty acid profiles of B. vulturis VM2412^T (1), B. avium VR2415^T (2), B. ginsengisoli KCTC 29226^T (3), and B. faecium JCM 11609^T (4).

The peptidoglycan of the two strains was found to contain alanine, glutamic acid, glycine and meso-diaminopimelic acid. Strain VR2415^T also contained aspartic acid within the peptidoglycan. Ribose, glucose, and galactose were found to be present as whole cell sugars. Polar lipids included phosphatidylglycerol (PG), diphosphatidylglycerol (DPG), phosphatidylethanolamine (PE, presumptively) and three unidentified glycolipids (Supplementary Figures 1A,B). In addition, strain VR2415^T produced an unidentified phospholipid. Both strains were found to have MK-7 as their predominant menaquinone.

Along with the closely related taxa B. ginsengisoli DCY80^T and B. faecium NCIB 9860^T, phenotypic characteristics of the type strains of B. squillarum and B. nesterenkovii were investigated because they were used for genotypic analyses using their available genome sequences. All type strains of Brachybacterium, Devriesea and Dermabacter spp. are non-motile and contain meso-diaminopimelic acid within their cell wall. Differential characteristics between the type strains of the three genera are presented in Table 3. PE was tentatively detected only in strains VM2412^T and VR2415^T.

TABLE 3

TABLE 3. Differential characteristics of type strains of Brachybacterium, Devriesea, and Dermabacter species: 1, B. vulturis VM2412^T; 2, B. avium VR2415^T; 3, B. ginsengisoli DCY80^T (Hoang et al., 2014); 4, B. faecium NCIB 9860^T (Collins et al., 1988); 5, B. squillarum M-6-3^T (Park et al., 2011); 6, B. nesterenkovii 35^T (Gvozdyak et al., 1992); 7, Devriesea agamarum IMP2^T (Martel et al., 2008); 8, Dermabacter vaginalis AD1-86^T (Chang et al., 2016); 9, D. hominis S69^T (Jones and Collins, 1988).

Genome Properties and Genetic Relatedness

Whole genome sequences were obtained from strains VM2412^T, VR2415^T, and B. ginsengisoli DCY80^T in this study. The genome properties of nine strains used in this study are presented in Table 4. The genome of strain VM2412^T is comprised of one circular chromosome of 3,796,663 bp long with 70.8 mol% G+C content. The genome of strain VR2415^T, with three assembled contigs, is 3,568,502 bp long with 70.4 mol% G+C content. The G+C contents (%) of the two strains are the lowest among all strains of current Brachybacterium species. Members of the genus Brachybacterium have higher G+C contents and larger genomes than the members of the genera Devriesea and Dermabacter. Strains VM2412^T, VR2415^T, B. ginsengisoli DCY80^T, and B. faecium DSM 4810^T each have a total of 9 ribosomal RNA-encoding genes (three 5S rRNAs, three 16S rRNAs and three 23S rRNAs). The genomes of all type strains of the Brachybacterium spp. analyzed have 50 tRNA-encoding genes except strain VM2412^T.

TABLE 4

TABLE 4. Genomic traits of type strains of Brachybacterium, Devriesea and Dermabacter species used in this study.

ANIb, ANIm, dDDH, and OrthoANI values were calculated to identify the genomic similarities of strains VM2412^T and VR2415^T to the seven strains of Brachybacterium, Devriesea, and Dermabacter species with available genome sequences. ANIb, ANIm, and dDDH values are presented in Table 5 and OrthoANI values are shown in Figure 2. Strains VM2412^T and VR2415^T share 99.3% sequence similarity with each other based on 16S rRNA gene sequences whereas ANIb, ANIm, OrthoANI, and dDDH values between the two strains are 89.0-89.2, 90.3, 89.7, and 38.0%, respectively. The two strains show high 16S rRNA gene sequence similarities with B. ginsengisoli DCY80^T (98.7% for VM2412^T and 98.4% for VR2415^T) and also have the closest genetic relatedness with B. ginsengisoli DCY80^T among the analyzed type strains of Brachybacterium species: 83.3-83.5% ANIb, 87.0% ANIm, 84.3% OrthoANI and 27.8% dDDH for VM2412^T; 82.8–83.2% ANIb, 86.7% ANIm, 83.9% OrthoANI, and 27.2% dDDH for VR2415^T. These findings support the classification of the strains as two distinct species according to cut-off values for species demarcation (95-96% for ANI and 70% for dDDH).

TABLE 5

TABLE 5. ANIb, ANIm, and dDDH values between pairs of type strains of Brachybacterium, Devriesea, and Dermabacter species.

FIGURE 2

FIGURE 2. OrthoANI values between type strains of Brachybacterium, Devriesea, and Dermabacter species.

The UBCG phylogenetic tree verified the phylogenetic relationships derived from the 16S rRNA phylogeny and genetic relatedness, showing a monophyletic clade comprising strains VM2412^T and VR2415^T which is closely related to B. ginsengisoli DCY80^T (Figure 3).

FIGURE 3

FIGURE 3. Genome-based phylogenetic tree constructed using an up-to-date bacterial core gene (UBCG) set. Arthrobacter globiformis NBRC 12137^T was used as an out-group. Numbers on nodes indicate the gene support index.

Comparative Genomic Analysis

Genome annotation was performed using the RAST annotation system. The following are comparisons of the genomic features derived from the genomes of strains VM2412^T and VR2415^T. Common features shared by all nine strains used in the comparative genomic study are mentioned at the end of each section. Differences among characteristics derived from genome sequences of the members of the three genera are given in Table 6.

TABLE 6

TABLE 6. Different genome-derived features of type strains of Brachybacterium, Devriesea, and Dermabacter species: 1, B. vulturis VM2412^T; 2, B. avium VR2415^T; 3, B. ginsengisoli DCY80^T; 4, B. faecium DSM 4810^T; 5, B. squillarum M-6-3^T; 6, B. nesterenkovii CIP104813^T; 7, Devriesea agamarum IMP2^T; 8, Dermabacter vaginalis AD1-86^T; 9, D. hominis NBRC 106157^T.

Defense Against Foreign DNA

Gene editing using CRISPR/Cas and restriction-modification systems is used for protection from foreign DNA. Strain VM2412^T has pre-crRNA processing-related genes cse1-4, 5e and cas1, 3 (Gesner et al., 2011). Genes encoding the type III restriction-modification system are present in both strains, but the type I restriction-modification system is present only in strain VM2412^T.

Motility

Genes belonging to the motility subsystem were not detected in the genomes of the two strains, consistent with the result of bacterial motility test using semi-solid agar. Similarly, the other strains analyzed do not encode genes that involved in the motility subsystem. Indeed, non-motility is one of the common features of the members of the family Dermabacteraceae (Stackebrandt, 2014).

Peptidoglycan Biosynthesis

Both strains contain gene murE, which catalyzes the addition of meso-diaminopimelic acid into peptidoglycan. This confirms the result of cell wall analysis using one-dimensional TLC. All the other strains analyzed also encode gene murE within their genomes. Members of the family Dermabacteraceae are known to have A4γ-type peptidoglycan, which contains meso-diaminopimelic acid (Stackebrandt, 2014).

Glycerophospholipid Metabolism

The enzyme for DPG synthesis is encoded in the genomes of the two strains. In both strains, enzymes for the biosynthesis of PG including phosphatidate cytidylyltransferase and CDP-diacylglycerol–glycerol-3-phosphate 3-phosphatidyltransferase are encoded but phosphatidylglycerophosphate is absent. Genes involved in the synthesis of phosphatidylcholine (PC), phosphatidylinositol (PI), phosphatidylserine (PS), and phosphatidylethanolamine (PE) are not present in the two strains while PE was detected by two-dimensional TLC. This latter identification needs to be confirmed; it is an atypical polar lipid in the majority of Brachybacterium species although PE has been detected in B. aquaticum (Kaur et al., 2016). All nine strains analyzed have genes for the synthesis of DPG and PG. The genome of D. agamarum IMP2^T includes the cls gene encoding cardiolipin (DPG) synthetase and the pcs gene encoding PC synthetase but DPG and PC were not detected in a previous characterization (Martel et al., 2008). This indicates potential discrepancies between genome-derived features and those obtained from some chemotaxonomic methods that will await detailed re-evaluation.

Sulfur Metabolism

In relation to metabolic pathways of sulfur compounds, strain VM2412^T possesses aryl sulfatase. However, none of the nine strains analyzed have genes for sulfur oxidation.

Anaerobic Respiration

Ferredoxin reductase, one of the anaerobic respiratory reductases, is present in the two strains. Furthermore, strain VM2412^T has arsenate reductase. In spite of the presence of anaerobic respiratory reductases, the strains are not able to grow under anaerobic conditions. Both strains possess respiratory nitrate reductase, which converts nitrate into nitrite. The ability to reduce nitrate was also detected using an API test kit. However, denitrification-related enzymes including nitrite reductase, nitric oxide reductase and nitrous oxide reductase are not present in any of the nine genomes analyzed. Nitrogenase involved in nitrogen fixation is also not encoded in the genomes of these nine strains.

Stress Response

Both strains contain enzymes involved in the biosynthesis of betaine from choline, including choline dehydrogenase (betA) and betaine aldehyde dehydrogenase (betB). Genes for ectoine biosynthesis (ectB, ectA, and ectC) are also present in the two strains. Osmolytes such as betaine and ectoine can protect cells against osmotic stress. Oxidative stress is suppressed by catalase, ferroxidase and manganese-dependent superoxide dismutase in the two strains. A synthetase for glutathione, known to be an important antioxidant, is encoded in the genome of strain VM2412^T. Catalase and manganese-dependent superoxide dismutase-encoding genes are also present in the other strains analyzed.

Resistance to Antibiotics

Both strains encode β-lactamase that may provide resistance to β-lactam antibiotics. Strain VR2415^T has a gene encoding VanW, a vancomycin B-type resistance protein. All the other strains analyzed here also contain genes encoding β-lactamase.

Conclusion

This study indicates that strains VM2412^T and VR2415^T represent two novel species of the genus Brachybacterium based on the results of phenotypic and genotypic analyses. Phylogenetic trees constructed using 16S rRNA gene sequences and core gene sets from whole genome sequences reveal the phylogenetic relationships between the isolates and other type strains of the genus Brachybacterium. In addition, this is the first study to reveal the taxogenomic relationships between members of the genus Brachybacterium. The following are descriptions of the newly isolated strains. The Digital Protologue database¹ TaxoNumbers for strains VM2412^T and VR2415^T are TA00602 and TA00603.

Description of Brachybacterium vulturis sp. nov.

Brachybacterium vulturis (vul.tu’ris. L. gen. n. vulturis from a vulture, the isolation source of the type strain).

Cells are coccoid, Gram-stain positive, non-motile, aerobic, catalase positive and oxidase negative. Colonies are 0.5-1 mm in diameter, ivory-colored, circular and low-convex with entire margins after 48 h incubation on TSA. Grows at 10-30°C, 0-15% (w/v) NaCl and pH 6-9. Optimal growth occurs at 30°C, 4% (w/v) NaCl and pH 7-8. Can reduce nitrates to nitrites and assimilate mannose, mannitol, N-acetylglucosamine and potassium gluconate. The major fatty acids are anteiso-C_15:0, anteiso-C_17:0 and iso-C_16:0. The cell wall contains meso-diaminopimelic acid, alanine, glutamic acid and glycine and whole cell sugars are ribose, glucose and galactose. The polar lipids are composed of phosphatidylglycerol, diphosphatidylglycerol, presumptively phosphatidylethanolamine and three unidentified glycolipids. The predominant menaquinone is MK-7. The genomic G+C content of the type strain is 70.8 mol%.

The type strain is VM2412^T (=KCTC 39996^T = JCM 32142^T), isolated from a fecal sample of an Andean condor (Vultur gryphus) in Seoul Grand Park, Gyeonggi-do, South Korea.

Description of Brachybacterium avium sp. nov.

Brachybacterium avium (a’vi.um. L. gen. pl. n. avium of birds).

Cells are coccoid, Gram-stain positive, non-motile, aerobic, catalase positive and oxidase negative. Colonies are ivory-colored, circular, low-convex and 0.5-1 mm in diameter after 48 h incubation on TSA. Grows at 10–30°C (optimal growth at 30°C), 0-15% (w/v) NaCl [optimal growth at 3-4% (w/v) NaCl] and pH 6-9 (optimal growth at pH 6). Can reduce nitrates to nitrites and assimilate glucose, arabinose, mannose, mannitol, N-acetylglucosamine, maltose and potassium gluconate. The major fatty acids are anteiso-C_15:0 and iso-C_16:0. Peptidoglycan contains alanine, glutamic acid, glycine, aspartic acid and meso-diaminopimelic acid and ribose, glucose, and galactose are present as whole cell sugars. The major polar lipids are phosphatidylglycerol, diphosphatidylglycerol, presumptively phosphatidylethanolamine, three unidentified glycolipids and an unidentified phospholipid. The predominant menaquinone is MK-7. The genomic G+C content of the type strain is 70.4 mol%.

The type strain is VR2415^T (=KCTC 39997^T = JCM 32143^T), isolated from a fecal sample of an Andean condor (Vultur gryphus) in Seoul Grand Park, Gyeonggi-do, South Korea.

Author Contributions

All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.

Funding

This work was supported by grants from the Mid-Career Researcher Program (NRF-2016R1E1A1A02921587) through the National Research Foundation of Korea (NRF), from the National Institute of Biological Resources (NIBR) funded by the Ministry of Environment (NIBR201701107) and from the Collaborative Genome Program for Fostering New Post-Genome industry (NRF-2015M3C9A2054299) through the National Research Foundation of Korea funded by the Ministry of Science ICT & Future Planning.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank staffs of Seoul Grand Park Zoo, especially J. R. Roh and S. D. Park for providing research permits to conduct our work and H. E. Shim for offering the avian feces.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.01809/full#supplementary-material

Footnotes

1. ^ http://imedea.uib-csic.es/dprotologue/

References

Bavbek, M., Caner, H., Arslan, H., Demirhan, B., Tunçbilek, S., and Altinörs, N. (1998). Cerebral Dermabacter hominis abscess. Infection 26, 181–183. doi: 10.1007/bf02771848

CrossRef Full Text | Google Scholar

Buczolits, S., Schumann, P., Weidler, G., Radax, C., and Busse, H. J. (2003). Brachybacterium muris sp. nov., isolated from the liver of a laboratory mouse strain. Int. J. Syst. Evol. Microbiol. 53, 1955–1960. doi: 10.1099/ijs.0.02728-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Chang, D. H., Rhee, M. S., and Kim, B. C. (2016). Dermabacter vaginalis sp. nov., isolated from human vaginal fluid. Int. J. Syst. Evol. Microbiol. 66, 1881–1886. doi: 10.1099/ijsem.0.000960

PubMed Abstract | CrossRef Full Text | Google Scholar

Chou, J. H., Lin, K. Y., Lin, M. C., Sheu, S. Y., Wei, Y. H., Arun, A. B., et al. (2007). Brachybacterium phenoliresistens sp. nov., isolated from oil-contaminated coastal sand. Int. J. Syst. Evol. Microbiol. 57, 2674–2679. doi: 10.1099/ijs.0.65019-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Collins, M. D., Brown, J., and Jones, D. (1988). Brachybacterium faecium gen. nov., sp. nov., a coryneform bacterium from poultry deep litter. Int. J. Syst. Evol. Microbiol. 38, 45–48. doi: 10.1099/00207713-38-1-45

CrossRef Full Text | Google Scholar

Collins, M. D., and Jones, D. (1981). Distribution of isoprenoid quinone structural types in bacteria and their taxonomic implication. Microbiol. Rev. 45, 316–354.

PubMed Abstract | Google Scholar

Cummins, C. S., and Harris, H. (1956). The chemical composition of the cell wall in some Gram-positive bacteria and its possible value as a taxonomic character. Microbiology 14, 583–600. doi: 10.1099/00221287-14-3-583

PubMed Abstract | CrossRef Full Text | Google Scholar

Felsenstein, J. (1981). Evolutionary trees from DNA sequences: a maximum likelihood approach. J. Mol. Evol. 17, 368–376. doi: 10.1007/BF01734359

CrossRef Full Text | Google Scholar

Funke, G., Stubbs, S., Pfyffer, G. E., Marchiani, M., and Collins, M. D. (1994). Characteristics of CDC group 3 and group 5 coryneform bacteria isolated from clinical specimens and assignment to the genus Dermabacter. J. Clin. Microbiol. 32, 1223–1228.

PubMed Abstract | Google Scholar

Gesner, E. M., Schellenberg, M. J., Garside, E. L., George, M. M., and MacMillan, A. M. (2011). Recognition and maturation of effector RNAs in a CRISPR interference pathway. Nat. Struct. Mol. Biol. 18, 688–692. doi: 10.1038/nsmb.2042

PubMed Abstract | CrossRef Full Text | Google Scholar

Gómez-Garcés, J. L., Oteo, J., Garcíìa, G., Aracil, B., Alós, J. I., and Funke, G. (2001). Bacteremia by Dermabacter hominis, a rare pathogen. J. Clin. Microbiol. 39, 2356–2357. doi: 10.1128/JCM.39.6.2356-2357.2001

PubMed Abstract | CrossRef Full Text | Google Scholar

Gontia, I., Kavita, K., Schmid, M., Hartmann, A., and Jha, B. (2011). Brachybacterium saurashtrense sp. nov., a halotolerant root-associated bacterium with plant growth-promoting potential. Int. J. Syst. Evol. Microbiol. 61, 2799–2804. doi: 10.1099/ijs.0.023176-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Gvozdyak, O. R., Nogina, T. M., and Schumann, P. (1992). Taxonomic study of the genus Brachybacterium: Brachybacterium nesterenkovii sp. nov. Int. J. Syst. Evol. Microbiol. 42, 74–78. doi: 10.1099/00207713-42-1-74

PubMed Abstract | CrossRef Full Text | Google Scholar

Hall, T. (2011). BioEdit: an important software for molecular biology. GERF Bull. Biosci. 2, 60–61.

Google Scholar

Hellebuyck, T., Martel, A., Chiers, K., Haesebrouck, F., and Pasmans, F. (2009). Devriesea agamarum causes dermatitis in bearded dragons (Pogona vitticeps). Vet. Microbiol. 134, 267–271. doi: 10.1016/j.vetmic.2008.08.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Heyrman, J., Balcaen, A., De Vos, P., Schumann, P., and Swings, J. (2002). Brachybacterium fresconis sp. nov. and Brachybacterium sacelli sp. nov., isolated from deteriorated parts of a medieval wall painting of the chapel of Castle Herberstein (Austria). Int. J. Syst. Evol. Microbiol. 52, 1641–1646. doi: 10.1099/00207713-52-5-1641

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoang, V. A., Kim, Y. J., Nguyen, N. L., and Yang, D. C. (2014). Brachybacterium ginsengisoli sp. nov., isolated from soil of a ginseng field. Int. J. Syst. Evol. Microbiol. 64, 3063–3068. doi: 10.1099/ijs.0.058388-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Jones, D., and Collins, M. D. (1988). Taxonomic studies on some human cutaneous coryneform bacteria: description of Dermabacter hominis gen. nov., sp. nov. FEMS Microbiol. Lett. 51, 51–55. doi: 10.1111/j.1574-6968.1988.tb02967.x

CrossRef Full Text | Google Scholar

Kaur, G., Kumar, N., Mual, P., Kumar, A., Kumar, R. M., and Mayilraj, S. (2016). Brachybacterium aquaticum sp. nov., a novel actinobacterium isolated from seawater. Int. J. Syst. Evol. Microbiol. 66, 4705–4710. doi: 10.1099/ijsem.0.001414

PubMed Abstract | CrossRef Full Text | Google Scholar

Kluge, A. G., and Farris, J. S. (1969). Quantitative phyletics and the evolution of anurans. Syst. Biol. 18, 1–32. doi: 10.1093/sysbio/18.1.1

CrossRef Full Text | Google Scholar

Lee, H. J., Cho, C. H., Kwon, M. J., Nam, M. H., Lee, K. N., and Lee, C. K. (2011). A patient with fatal septicemia caused by a rare pathogen Dermabacter Hominis. Infect. Chemother. 43, 86–88. doi: 10.3947/ic.2011.43.1.86

CrossRef Full Text | Google Scholar

Lee, I., Kim, Y. O., Park, S. C., and Chun, J. (2016). OrthoANI: an improved algorithm and software for calculating average nucleotide identity. Int. J. Syst. Evol. Microbiol. 66, 1100–1103. doi: 10.1099/ijsem.0.000760

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Y., Xie, Q. Y., Shi, W., Li, L., An, J. Y., Zhao, Y. M., et al. (2014). Brachybacterium huguangmaarense sp. nov., isolated from Lake sediment. Int. J. Syst. Evol. Microbiol. 64, 1673–1678. doi: 10.1099/ijs.0.052464-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Y., Zhai, L., Yao, S., Cao, Y., Cao, Y., Zhang, X., et al. (2015). Brachybacterium hainanense sp. nov., isolated from noni (Morinda citrifolia L.) branch. Int. J. Syst. Evol. Microbiol. 65, 4196–4201. doi: 10.1099/ijsem.0.000559

PubMed Abstract | CrossRef Full Text | Google Scholar

Martel, A., Pasmans, F., Hellebuyck, T., Haesebrouck, F., and Vandamme, P. (2008). Devriesea agamarum gen. nov., sp. nov., a novel actinobacterium associated with dermatitis and septicaemia in agamid lizards. Int. J. Syst. Evol. Microbiol. 58, 2206–2209. doi: 10.1099/ijs.0.65478-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Meier-Kolthoff, J. P., Auch, A. F., Klenk, H. P., and Göker, M. (2013). Genome sequence-based species delimitation with confidence intervals and improved distance functions. BMC Bioinformatics 14:60. doi: 10.1186/1471-2105-14-60

PubMed Abstract | CrossRef Full Text | Google Scholar

MIDI (1999). Sherlock Microbial Identification System Operating Manual, Version 3.0. Newark, NJ: MIDI, Inc.

Google Scholar

Minnikin, D. E., O’Donnell, A. G., Goodfellow, M., Alderson, G., Athalye, M., Schaal, A., et al. (1984). An integrated procedure for the extraction of bacterial isoprenoid quinones and polar lipids. J. Microbiol. Methods. 2, 233–241. doi: 10.1016/0167-7012(84)90018-6

CrossRef Full Text | Google Scholar

Na, S. I., Kim, Y. O., Yoon, S. H., Ha, S. M., Baek, I., and Chun, J. (2018). UBCG: up-to-date bacterial core gene set and pipeline for phylogenomic tree reconstruction. J. Microbiol. 56, 280–285. doi: 10.1007/s12275-018-8014-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, S. K., Kim, M. S., Jung, M. J., Nam, Y. D., Park, E. J., Roh, S. W., et al. (2011). Brachybacterium squillarum sp. nov., isolated from salt-fermented seafood. Int. J. Syst. Evol. Microbiol. 61, 1118–1122. doi: 10.1099/ijs.0.022517-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, Y. K., Lee, K. M., Lee, W. K., Cho, M. J., Lee, H. S., Cho, Y. G., et al. (2016). Dermabacter jinjuensis sp. nov., a novel species of the genus Dermabacter isolated from a clinical specimen. Int. J. Syst. Evol. Microbiol. 66, 2573–2577. doi: 10.1099/ijsem.0.001092

PubMed Abstract | CrossRef Full Text | Google Scholar

Radtke, A., Bergh, K., Øien, C. M., and Bevanger, L. S. (2001). Peritoneal dialysis-associated peritonitis caused by Dermabacter hominis. J. Clin. Microbiol. 39, 3420–3421. doi: 10.1128/JCM.39.9.3420-3421.2001

PubMed Abstract | CrossRef Full Text | Google Scholar

Renvoise, A., Aldrovandi, N., Raoult, D., and Roux, V. (2009). Helcobacillus massiliensis gen. nov., sp. nov., a novel representative of the family Dermabacteraceae isolated from a patient with a cutaneous discharge. Int. J. Syst. Evol. Microbiol. 59, 2346–2351. doi: 10.1099/ijs.0.003319-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Richter, M., Rosselló-Móra, R., Oliver Glöckner, F., and Peplies, J. (2016). JSpeciesWS: a web server for prokaryotic species circumscription based on pairwise genome comparison. Bioinformatics 32, 929–931. doi: 10.1093/bioinformatics/btv681

PubMed Abstract | CrossRef Full Text | Google Scholar

Saitou, N., and Nei, M. (1987). The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 4, 406–425. doi: 10.1093/oxfordjournals.molbev.a040454

PubMed Abstract | CrossRef Full Text | Google Scholar

Schubert, K., Ludwig, W., Springer, N., Kroppenstedt, R. M., Accolas, J. P., and Fiedler, F. (1996). Two coryneform bacteria isolated from the surface of french gruyère and beaufort cheeses are new species of the genus Brachybacterium: Brachybacterium alimentarium sp. nov. and Brachybacterium tyrofermentans sp. nov. Int. J. Syst. Evol. Microbiol. 46, 81–87. doi: 10.1099/00207713-46-1-81

PubMed Abstract | CrossRef Full Text | Google Scholar

Schumann, P. (2011). “5 - Peptidoglycan Structure,” in Methods in Microbiology, eds F. Rainey and A. Oren (Cambridge, MA: Academic Press), 101–129. doi: 10.1016/B978-0-12-387730-7.00005-X

CrossRef Full Text | Google Scholar

Singh, H., Du, J., Yang, J. E., Yin, C. S., Kook, M., and Yi, T. H. (2016). Brachybacterium horti sp. nov., isolated from garden soil. Int. J. Syst. Evol. Microbiol. 66, 189–195. doi: 10.1099/ijsem.0.000696

PubMed Abstract | CrossRef Full Text | Google Scholar

Stackebrandt, E. (2014). “The Family Dermabacteraceae,” in The Prokaryotes, eds E. Rosenberg, E. F. DeLong, S. Lory, E. Stackebrandt, and F. Thompson (Heidelberg: Springer), 289–299. doi: 10.1007/978-3-642-30138-4_170

CrossRef Full Text | Google Scholar

Takeuchi, M., Fang, C. X., and Yokota, A. (1995). Taxonomic Study of the Genus Brachybacterium: proposal of Brachybacterium conglomeratum sp. nov., nom. rev., Brachybacterium paraconglomeratum sp. nov., and Brachybacterium rhamnosum sp. nov. Int. J. Syst. Evol. Microbiol. 45, 160–168. doi: 10.1099/00207713-45-1-160

CrossRef Full Text | Google Scholar

Tamura, K., Stecher, G., Peterson, D., Filipski, A., and Kumar, S. (2013). MEGA6: molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 30, 2725–2729. doi: 10.1093/molbev/mst197

PubMed Abstract | CrossRef Full Text | Google Scholar

Tittsler, R. P., and Sandholzer, L. A. (1936). The use of semi-solid agar for the detection of bacterial motility. J. Bacteriol. 31, 575–580.

PubMed Abstract | Google Scholar

Van Bosterhaut, B., Boucquey, P., Janssens, M., Wauters, G., and Delmée, M. (2002). Chronic osteomyelitis due to Actinomyces neuii subspecies neuii and Dermabacter hominis. Eur. J. Clin. Microbiol. Infect. Dis. 21, 486–487. doi: 10.1007/s10096-002-0747-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoon, S. H., Ha, S. M., Kwon, S., Lim, J., Kim, Y., Seo, H., et al. (2017). Introducing EzBioCloud: a taxonomically united database of 16S rRNA gene sequences and whole-genome assemblies. Int. J. Syst. Evol. Microbiol. 67, 1613–1617. doi: 10.1099/ijsem.0.001755

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, G., Zeng, G., Cai, X., Deng, S., Luo, H., and Sun, G. (2007). Brachybacterium zhongshanense sp. nov., a cellulose-decomposing bacterium from sediment along the Qijiang River, Zhongshan City, China. Int. J. Syst. Evol. Microbiol. 57, 2519–2524. doi: 10.1099/ijs.0.64968-0

PubMed Abstract | CrossRef Full Text | Google Scholar