Superficieibacter electus gen. nov., sp. nov., an Extended-Spectrum β-Lactamase Possessing Member of the Enterobacteriaceae Family, Isolated From Intensive Care Unit Surfaces

Two Gram-negative bacilli strains, designated BP-1(T) and BP-2, were recovered from two different Intensive Care Unit surfaces during a longitudinal survey in Pakistan. Both strains were unidentified using the bioMerieux VITEK MS IVD v2.3.3 and Bruker BioTyper MALDI-TOF mass spectrometry platforms. To more precisely determine the taxonomic identity of BP-1(T) and BP-2, we employed a biochemical and phylogenomic approach. The 16S rRNA gene sequence of strain BP-1(T) had the highest identity to Citrobacter farmeri CDC 2991-81(T) (98.63%) Citrobacter amalonaticus CECT 863(T) (98.56%), Citrobacter sedlakii NBRC 105722(T) (97.74%) and Citrobacter rodentium NBRC 105723(T) (97.74%). The biochemical utilization scheme of BP-1(T) using the Analytic Profile Index for Enterobacteriaceae (API20E) indicated its enzymatic functions are unique within the Enterobacteriaceae but most closely resemble Kluyvera spp., Enterobacter cloacae and Citrobacter koseri/farmeri. Phylogenomic analysis of the shared genes between BP-1(T), BP-2 and type strains from Kluyvera, Citrobacter, Escherichia, Salmonella, Kosakonia, Siccibacter and Shigella indicate that BP-1(T) and BP-2 isolates form a distinct branch from these genera. Average Nucleotide Identity analysis indicates that BP-1(T) and BP-2 are the same species. The biochemical and phylogenomic analysis indicate strains BP-1(T) and BP-2 represent a novel species from a new genus within the Enterobacteriaceae family, for which the name Superficieibacter electus gen. nov., sp. nov., is proposed. The type strain is BP-1(T) (= ATCC BAA-2937, = NBRC 113412).

Introduction

The hospital built environment is a key source of nosocomial pathogens, which can cause high mortality infections in vulnerable patient populations (Oberauner et al., 2013). The danger from these infections is exacerbated by the higher levels of antibiotic resistance harbored by nosocomial pathogens compared to their community-associated relatives (Horcajada et al., 2013). As many pathogens can persist on surfaces for extended periods of time, microbiological and molecular surveillance of high-contact surfaces in health-care settings is an important feature of effective infection control and prevention strategies (Russotto et al., 2015). The isolates reported herein were recovered as part of such a surveillance program of intensive care unit (ICU) surfaces at a tertiary care hospital in Pakistan in 2016.

The Enterobacteriaceae is a diverse family of Gammaproteobacteria that includes many common human pathogens (Paradis et al., 2005; Potter et al., 2016). Enterobacteriaceae taxonomy is complicated, as demonstrated by the use of multilocus sequence analysis to reclassify several Enterobacter species into the Lelliottia, Pluralibacter and Kosakonia genera (Brady et al., 2013). Next-generation sequencing has significantly improved our understanding of Enterobacteriaceae taxonomy by identifying similarities between species of different genera (e.g., Salmonella subterranean and Escherichia hermannii into Atlantibacter subterranea and Atlantibacter hermannii, respectively) or resolving closely related species within the same genus (e.g., Klebsiella pneumoniae and Klebsiella variicola) (Holt et al., 2015; Hata et al., 2016). Many Enterobacteriaceae strains are medically relevant not only as etiological agents of disease but also as reservoirs for transferable antibiotic resistance genes (Sirot et al., 1988). Thus, documenting the appearance of novel Enterobacteriaceae species in hospital settings is an important component in surveilling and preventing emerging infectious diseases (Iredell et al., 2016).

Materials and Methods

Bacterial Isolation

We isolated strain BP-1(T) (bedside rail) and BP-2 (nurse call button) from an intensive care room at a tertiary care hospital in Pakistan with the Eswab collection device and transport system (Becton Dickenson & Company, Franklin Lakes, NJ) in June 2016. Surfaces were swabbed in triplicate during an ongoing yearlong longitudinal sampling study. Swabs were sent to the US at room temperature and arrived in the US within ~3 days of sampling. Samples from Pakistan were received by the Burnham lab under approval of the Washington University IBC Protocol number 6,572. Swabs were cultured within 24 h of arrival in the US laboratory site. The isolates were recovered on MacConkey agar with 5 μg/ml cefotaxime (Hardy Diagnostics, Santa Maria, CA) and re-streaked on blood agar using the four-quadrant method to recover distinct single colonies. The isolates were unidentified by the MALDI-TOF VITEK MS IVD v2.3.3 (bioMerieux, Durham, NC) and Bruker BioTyper (Bruker Daltonics, Billerica, MA) mass spectrometry systems. For subsequent genomic analysis, frozen stocks of each isolate were made by creating a dense suspension of the isolate in Tryptic Soy Broth (Sigma Aldrich, St. Louis, MO) with 15% Glycerol (Hardy Diagnostics, Santa Maria, CA).

Biochemical Tests and Growth Assays

Glycerol stocks of strain BP-1(T) and BP-2 were streaked separately onto 5% sheep′s blood agar (Hardy Diagnostics, Santa Maria, CA), HardyCHROM™ ESBL Agar (Hardy Diagnostics, Santa Maria, CA), Hektoen Enteric Agar (Remel, Lenexa, KS) and MacConkey agar (Hardy Diagnostics, Santa Maria, CA) plates using four quadrant streaking. Plates were incubated for 35°C in air and then imaged. Additional streaks of blood agar plates were incubated at either 4, 25, 35, or 42°C. Blood agar plates were also stored in an anaerobic bag (BD, Franklin Lakes, NJ) at 35°C. Phenotypic assays to assess indole production, catalase, urease and motility were performed. The isolates were also evaluated on lysine iron agar (LIA) (Remel, Lenexa, KS), triple sugar iron (TSI) agar (Remel, Lenexa, KS), Bile esculin agar (Remel, Lenexa, KS) and with an oxidation/fermentation glucose test (OF glucose) (Remel, Lenexa, KS), all performed according to manufacturer′s recommendations. Each isolate was also evaluated using the Analytical Profile Index (API^®) 20-E Enterobacteriaceae identification kit (bioMerieux, Durham, NC). API values of related taxa were obtained from for Kluyvera, Enterobacter cloacae and Citrobacter koseri (referred to as Citrobacter amalonaticus biogroup 1 by Farmer et al., 1981, 1985). Similar to how this information was used for Escherichia albertii delineation, characters are scored as “+” if > = 85% of the strains are positive, “−” if > = 85% of the strains are negative, “v+” if 50–85% of strains are positive and “v–” if 50-85% of strains are negative (Huys et al., 2003).

Transmission Electron Microscopy

Glycerol stocks of strain BP-1(T) and BP-2 were plated onto MacConkey agar with ceftriaxone (5 μg/ml) overnight. A single colony of each isolate was separately inoculated into 5 ml of LB broth and grown until log phase at 37°C with orbital shaking at 220 rpm.

For analysis of whole bacteria by negative staining, the isolates absorbed onto formvar/carbon-coated copper grids for 2 min. Grids were washed in dH₂O and stained with 1% aqueous uranyl acetate (Ted Pella Inc., Redding CA) for 1 min. Excess liquid was gently wicked off and grids were air dried. Bacteria viewed on a JEOL 1200 EX transmission electron microscope (JEOL USA Inc., Peabody, MA) equipped with an AMT 8-megapixel digital camera (Advanced Microscopy Techniques, Woburn, MA).

For ultrastructural analysis of cross-sections of bacteria, the isolates were fixed in 2% paraformaldehyde/2.5% glutaraldehyde in 100 mM cacodylate buffer, pH 7.2 for 1 h at room temperature. Samples were washed in cacodylate buffer and postfixed in 1% osmium tetroxide (Polysciences Inc.) for 1 h. Samples were then rinsed extensively in dH₂0 prior to en bloc staining with 1% aqueous uranyl acetate (Ted Pella Inc., Redding, CA) for 1 h. Following several rinses in dH₂0, samples were dehydrated in a graded series of ethanol and embedded in Eponate 12 resin (Ted Pella Inc.). Sections of 95 nm were cut with a Leica Ultracut UCT ultramicrotome (Leica Microsystems Inc., Bannockburn, IL), stained with uranyl acetate and lead citrate and viewed on the transmission electron microscope.

Illumina Whole Genome Sequencing

Genomic DNA was extracted using the bacteremia kit (Qiagen, Germantown, MD), from a suspension of ~10 colonies of strain BP-1(T) and BP-2 after overnight growth on blood agar (Hardy Diagnostics, Santa Maria, CA). 0.5 ng of DNA was used as input for constructing Nextera Illumina sequencing libraries (Illumina, San Diego, CA) (Baym et al., 2015). Sample libraries were sequenced on an Illumina NextSeq 550 to obtain ~2.5 million 2 × 150 bp reads. Raw reads had Illumina adapters removed with Trimmomatic and were decontaminated with Deconseq using the commands: “java -Xms1024m -Xmx1024m -jar $TRIMMOMATIC_HOME/trimmomatic-0.36.jar PE -phred33 -trimlog < trimlog.txt> < input forward reads> < input reverse reads> < output paired forward reads> < output unpaired forward reads> < output paired reverse reads> < output unpaired reverse reads> ILLUMINACLIP:/opt/apps/trimmomatic/0.36/adapters/NexteraPE-PE.fa:2:30:10:1:true” and “deconseq.pl -f < paired input reads> -o < paired output reads> -dbs hsref,” respectively (Schmieder and Edwards, 2011; Bolger et al., 2014). Paired forward and reverse reads were used as input for de-novo assembly with SPAdes with the command: “spades.py -k 21,33,55,77 –careful –pe2-1 < input forward reads> –pe2-2 < input reverse reads> -o < output directory>” (Bankevich et al., 2012). Assembly metrics were assessed by using QUAST on all contigs >500 bp in length using the webserver (http://quast.bioinf.spbau.ru/) (Gurevich et al., 2013). Draft whole genome sequences for strain BP-1(T) and BP-2 have been deposited with NCBI BioProject: PRJNA395420.

16S rRNA Similarity

16S rRNA sequences from the BP-1(T) file were identified using Barrnap with the command: “barrnap –quiet < contigs.fasta>” (https://github.com/tseemann/barrnap) and manually retrieved. The sequence is on contig_69 from base pairs 22 to 1559. The complete 16S rRNA sequence was submitted to the EzBioCloud identify service on 08/31/17 (Kim et al., 2014; Yoon et al., 2017). The 16S rRNA sequences for the top 10 most similar taxa were retrieved and aligned using MUSCLE (Edgar, 2004) using the command “muscle -in < input_multifasta> -out < output_aligned_multifasta>.” The aligned 16S rRNA multifasta was made into a phylogenetic tree with 1,000 bootstraps using raxML with the command: “raxmlHPC -s < input_aligned_multifasta> -n output -m GTRGAMMA -p 100 -f a -N 1000 -x 12345.” The resulting newick tree was visualized in FigTree with bootstrap values depicted as node label. The aligned multifasta sequences were uploaded into Jalview v2.0 and colored by nucleotide (Waterhouse et al., 2009). The 16S rRNA sequence has been submitted to Genbank under accession number MG866003.

Core Genome Phylogeny

Fifty-six genomes from type strains in species of Escherichia, Shigella, Citrobacter, Salmonella, Kluyvera, Klebsiella, Kosakonia, Siccibacter and Pasteurella (Table S1) were obtained from NCBI Genomes on 09/01/17. The protein coding sequences of these genomes and strains BP-1(T) and BP-2 were annotated using Prokka with the previously described command (Seemann, 2014). The core genome, representing genes shared by 100% of all isolates in the cohort was constructed and aligned using the PRANK module within roary with the following command: “roary -e –group_limit 700000 -p 8 -i < percent identity> *.gff” (Page et al., 2015). Eighty percent of gene identity was used for core genome alignment of Enterobacteriaceae genomes plus the P. multocida outgroup and for core genome alignment of only the Enterobacteriaceae genomes.

Aligned multifasta files of the 48 genes shared by all Enterobacteriaceae genomes and P. multocida (Table S2) was constructed into a maximum likelihood phylogenetic tree using 1,000 bootstraps in RAxML with the following command: “raxmlHPC -s core_gene_alignment.aln -n output -m GTRGAMMA -o < P. multocida name>-p 100 -f a -N 1000 -x 12345” (Stamatakis, 2014). Trees were visualized in FigTree (http://tree.bio.ed.ac.uk/software/figtree/) with bootstrap values depicted as node labels.

Aligned multifasta files of the 750 genes shared by all Enterobacteriaceae genomes (Table S2) was constructed into a NeighborNet tree using the EqualAngle method in SplitsTree4 (v4.14.5) (Huson and Bryant, 2006).

Average Nucleotide Identity

Whole genome sequences of strains BP-1(T) and BP-2 were compared against each other using the mummer method of average nucleotide identity via the jspecies webserver (http://jspecies.ribohost.com/jspeciesws/) on 09/18/17 (Richter et al., 2016). The contig file for strains BP-1(T), BP-2 and the set of Enterobacteriaceae genomes (Table S1) were compared using the mummer method of average nucleotide identity via the standalone pyANI tool (https://github.com/widdowquinn/pyani) with command “python3 average_nucleotide_identity.py -i < input_fasta_directory> -o < outout_directory> -m ANIm –nucmer_exe < mummer_path> -v –f.” The Hadamard matrix representing the multiplication of the percentage identity and alignment length was converted into a dendrogram using scipy.cluster.hierarchy.dendrogram. within python3 using “euclidian” metric and “weighted” method. The dendrogram was viewed in matplotlib.

Phenotypic Antibiotic Susceptibility Testing

Antimicrobial susceptibility testing was performed on the isolates using Kirby Bauer disk diffusion on Mueller Hinton Agar (Hardy Diagnostics, Santa Maria, CA) according to CLSI standards (CLSI, 2013). Results were interpreted using CLSI M100 criteria for Enterobacteriacace (CLSI, 2013).

Antibiotic Resistance Gene Identification and Visualization

Antibiotic resistance genes (ARGs) from the draft whole genome assemblies of strains BP-1(T) and BP-2 were annotated by submission to the ResFinder webserver (https://cge.cbs.dtu.dk/services/ResFinder/) on 05/18/17 (Zankari et al., 2012). We annotated sequence divergent ARGs using Resfams with the command “annotate_functional_selections.py –Resfams_only -contigs < contig file> -o < output directory>” (Gibson et al., 2015). Putative resistant determinants were then mapped to observed phenotypic resistance in accordance with Resfams and ResFinder classification. Contig_4 in the BP-1(T) assembly was submitted to BLASTN against the nr/nt database on 01/22/18 (Camacho et al., 2009).

Strain BP-1(T) and BP-2 were annotated for protein coding sequences by Prokka using the command: “Prokka < input fasta> < output directory>” (Seemann, 2014). Contigs 54 and 72 from BP-1(T) had BLAST similarity compared against Contig 42 from BP-2 using default settings on WebAct (webact.org/WebACT/generate) (Abbott et al., 2007). BLAST similarity was viewed using EasyFig for sequences with e value 0.001 and length >100 bp (Sullivan et al., 2011).

Results

Colony Appearance and Ultrastructural Analysis

Both strains BP-1(T) and BP-2 appeared identical on agar plates. BP-1(T) formed circular, shiny, mucoid, gray colonies and were non-hemolytic on blood agar plates (Figure 1A). The colonies were small, blue and rigid on HardyCHROM™ ESBL Agar (Figure 1B). BP-1(T) colonies were bright pink on MacConkey agar (i.e., lactose fermenting) with an opaque zone surrounding the edge of the colony (Figure 1C). Coral colored colonies were produced on HE agar (Figure 1D).

FIGURE 1

Figure 1. Colony morphology of strain BP-1(T) on (A) blood, (B) HardyCHROM™ ESBL, (C) MacConkey and (D) Hektoen Enteric agars after overnight growth at 35°C.

Analysis by negative staining and transmission electron microscopy revealed that BP-1(T) and BP-2 are rod shaped and ~2.7 μm in length (Figures 2A,B). Cross sections from resin embedding also show rod shaped bacteria (Figures 2C,D) with two electron dense lipid bilayers of the outer membrane and cytoplasmic membrane (Figure 2E). No evidence of flagella was observed, but pili-like structures were present between a group of isolates (Figures 2F,G).

FIGURE 2

Figure 2. Negative staining of log-phase strains (A) BP-1(T) and (B) BP-2with 1% uranyl acetate show a bacilli shape structure ~2.7 μm in length. Cross sections of strains (C) BP-1(T) and (D) BP-2 from resin embedding identify rod shaped bacteria. The isolates have a double membrane structure characteristic with of Gram-negative bacteria (E). The white arrows indicate the outer membrane, periplasm and cytoplasmic membrane. Negative staining showing several pilus-like structures (black arrows) protruding from strain BP-1(T) isolates (F and enlarged insert shown in G).

Biochemical and Temperature Studies

Strains BP-1(T) and BP-2 had identical results from phenotypic assays (Table 1). Aerobically, BP-1(T) showed no growth at 4°C, scant growth at 25 and 42°C, but robust growth at 35°C (Table 1). BP-1(T) and BP-2 both showed robust growth anaerobically at 35°C using an anaerobic bag (Table 1). The isolates were negative for motility, oxidase, indole, catalase and urease. Both strains BP-1(T) and BP-2 were positive for gas production in TSI and LIA reactions (Table 1). The strains grew on bile esculin agar and can ferment glucose (Table 1).

TABLE 1

Table 1. Phenotypic evaluation of strains BP-1(T) and BP-2.

As evaluated using the API 20E, BP-1(T) and BP-2 were positive for arabinose, amygdalin, melibiose, sacharose, rhamnose, sorbitol and mannitol fermentation, but negative for inositol fermentation (Table 2). BP-1(T) and BP-2 can utilize citrate and decarboxylate ornithine for metabolism in addition to utilizing disaccharides (melibiose) (Table 2). These biochemical test results were unable to assign a high confidence genus identification to strains BP-1(T) and BP-2. The analytic profile index (API) algorithm indicates that strains BP-1(T) and BP-2 most resembles Kluyvera spp. and to a lesser extent E. cloacae and Citrobacter koseri/farmeri (Table 3). Analysis of the biochemical utilization profile from these related taxa indicate that strains BP-1(T)/BP-2 and Kluyvera cryocrescens ATCC 33435 are both negative for arginine dihydrolase activity but E. cloacae and C. koseri are positive for it. Additionally, strains BP-1(T)/BP-2, E. cloacae and C. koseri can ferment D-sorbitol but K. cryocrescens ATCC 33435 cannot. Both BP-1(T) and BP-2 were tested as non-motile but K. cryocrescens ATCC 33435, E. cloacae and C. koseri are positive. Accordingly, BP-1(T)/BP-2 can be uniquely distinguished from its closest neighbors based on being negative for arginine hydrolase activity, positive for D-sorbitol fermentation and non-motile.

TABLE 2

Table 2. Biochemical utilization of strains BP-1(T) and BP-2 with API20E.

TABLE 3

Table 3. API20E interpretation.

Phylogenomic Analysis

We used a phylogenomic approach to better delineate the taxonomic context of strains BP-1(T) and BP-2 within the Enterobacteriaceae. Submission of the 1,512-base pair 16S rRNA locus from strain BP-1(T) to the EzBioCloud taxonomic database indicated that strain BP-1(T) had the highest identity to Citrobacter farmeri CDC 2991-81(T), C. amalonaticus CECT 863(T), C. sedlakii NBRC 105722(T) and C. rodentium NBRC 105723(T) (Table 4). The maximum likelihood tree from alignment of these 16S rRNA sequences placed BP-1(T) in-between a group of Citrobacter species [C. sedlakii NBRC 105722(T), C. rodentium NBRC 105723(T), C. farmeri CDC 2991-81(T) and C. amalonaticus CECT 863(T)] and the other Enterobacteriaceae taxa including Salmonella enterica subsp. Arizonae ATCC 13314(T), Kosakonia sacchari SP1(T), C. europaeus 97/79(T), Enterobacter hormaechei susbp. steigerwaltii DSM 16691(T), an uncharacterized Enterobacter strain termed LGIT s RIT-PI-d and Pluralibacter gergoviae JCM 1234(T) (Figure S1). The greatest amount of variation in 16S rRNA identity occurred at the likely V3 region, around nucleotides 426–450 of the BP-1(T) sequence, using standard E. coli nomenclature (Chakravorty et al., 2007). The 6mer of GAGAAT in the BP-1(T) sequence from sites 426 to 435 is unique in comparison to the other sequences at this position. Although the GAG nucleotides from 426 to 428 are shared by K. sacchari SP1(T)/P. gergoviae JCM 1234(T) and the AT site at 434–435 is shared by C. farmeri CDC 2991-81^T, the A residue at position 429 is discriminatory, as all other taxa contain a thymine at this site (Figure S2). As 98.7–98.65% is a proposed standard for species delineation via rRNA gene sequencing, we endeavored to gain further insight into the taxonomic identity of BP-1(T) and BP-2 based on their whole genome sequence characteristics.

TABLE 4

Table 4. Top hits of strain BP-1(T) 16S rRNA sequence.

To determine the relative taxonomic identity of strains BP-1(T) and BP-2 beyond 16S rRNA similarity, we constructed an equal angle neighbor net tree produced from the aligned 750 core genes shared between type strains in Enterobacter, Siccibacter, Escherichia, Enterobacter, Kosakonia, Kluyvera, Klebsiella, Shigella, Salmonella and Citrobacter at >80% identity (Table S1). Analysis of the network indicates that strains BP-1(T) and BP-2 cluster together and distinctly from these major Enterobacteriaceae taxa (Figure 3A).

FIGURE 3

Figure 3. (A) Unrooted core genome NeighborNet tree from an alignment of the 750 genes shared by Enterobacteriaceae genomes with >80% nucleotide identity. Scale bar represents number of nucleotide substitutions. Strain BP-1(T) and BP-2 are a separate cluster from other Enterobacteriaceae genera. (B) Rooted Core genome phylogeny of Enterobacteriaceae genomes with an outgroup for comparison (Pasteurella multocida), based on 48 genes shared at >80% nucleotide identity by all genomes. Both strains BP-1(T) and BP-2 isolates clustered together in a clade separate from existing Enterobacteriaceae type strains. Bootstrap support values < 95 are depicted as branch labels. Scale bar represents number of nucleotide substitutions per site.

A maximum likelihood tree from the 48 genes shared between the same Enterobacteriaceae genomes and a Pasteurella multocida outgroup with >80% identity again placed strains BP-1(T) and BP-2 in a cluster separate from existing genera. The BP-1(T)/BP-2 branch separates the genomes analyzed into two major clades composed of Escherichia [excluding Escherichia vulneris NBRC 192420(T), which other groups have suggested is misclassified (Walk et al., 2009)], Shigella, Citrobacter, Salmonella in one branch and Klebsiella, Kluyvera, Enterobacter, Siccibacter and Kosakonia on the other branch (Figure 3B).

Pairwise Average Nucleotide Identity (ANI) analysis using MUMer between strains BP-1(T) and BP-2 indicate that they are the same species, as the ANI value of 99.97% is greater than the proposed species cutoff of 95% (Thompson et al., 2013). The dendrogram of a distance matrix representing the product of ANI and percentage of the genome alignment created a clustermap that is similar to the rooted Phylogenetic tree with BP-1(T) and BP-2 separating the Escherichia, Shigella, Citrobacter and Salmonella clade from the Klebsiella, Kluyvera, Enterobacter, Siccibacter and Kosakonia branch (Figure 4). While the ordering of the two Citrobacter clades, Salmonella and Escherichia/Shigella is similar between the ANI dendrogram and the core genome phylogenetic tree, in the ANI dendrogram Kosakonia and not Kluyvera/Klebsiella is the closest genera to BP-1(T)/BP-2.

FIGURE 4

Figure 4. Dendrogram from the distance matrix representing the Hadamard matrix between average nucleotide identity and percent genome aligned. Strain BP-1 and BP-2 are split between two larger clades of Enterobacteriaceae bacteria, similar to Figure 3B.

Phenotypic Antibiotic Susceptibility and Antibiotic Resistance Gene Analysis

Both strains BP-1(T) and BP-2 had phenotypic resistance to all cephalosporins tested except cefepime and cefotetan (Table 5). This susceptibility profile suggests carriage of an extended-spectrum β-lactamase (ESBL). Both isolates were resistant to ampicillin and aztreonam but susceptible to carbapenems (meropenem and imipenem) and β-lactam/β-lactamase inhibitor combinations (pipercillin-tazobactam, ceftazidime-avibactam and ampicillin-sulbactam; Table 5). The isolates were resistant to gentamicin but susceptible to amikacin (Table 5). Each isolate was additionally resistant to trimethoprim-sulfamethoxazole and susceptible to doxycycline, minocycline and tigecycline (Table 5). BP-1(T) was susceptible to fosfomycin but BP-2 had intermediate resistance (Table 5). Both isolates shared the same antibiotic resistance genes. ResFinder identified bla_SHV−12, bla_TEM−1b, qnrB2, aac(6′)IB-cr, strB, aac(6′)-IIC, aacA4, aph(3′)-Ia, strA, sul1, sul2, dfrA18, catA2 and ere(A) (Table 5). Resfams identified the above genes except for sul1, sul2 and dfrA18, but found an additional Class A β-lactamase (Table 5). In BP-1(T) the class A β-lactamase was on contig 4, believed to be chromosomal as the highest blastn match was the complete genome of C. amalonaticus FDAARGOS and had 96% identity with 68% query coverage to a HER family β-lactamase (WP_023281079.1). Preceding it was a tetR response regulator with similarity to a protein from Klebsiella sp. RIT-PI-d (WP_049839935.1).

TABLE 5

Table 5. Phenotypic antibiotic susceptibility and putative resistance determinants.

To gain further insight into the potential for antibiotic resistance gene dissemination from strains BP-1(T) and BP-2, we analyzed the draft genomes of BP-1(T) and BP-2 for contigs that had ARGs co-localized with mobile gene elements. An IS1380 family transposase was found co-localized in both genomes with an erythromycin resistance gene ere(A), a putative rifampin resistance gene that was unidentified by ResFinder or Resfams, a putative aminoglycoside N-acetyltransferase unidentified by Resfinder or Resfams, a hypothetical protein and aac(6′)-IIC (Figure 5).

FIGURE 5

Figure 5. EasyFig construction showing BLAST-based sequence similarity between genetic structure of transposase and putative ARGs on the same contig in BP-1(T) and two different contigs in BP-2. Scale bar represents sequence length.

Discussion

The purpose of our investigation was to determine the taxonomic identity of strains BP-1(T) and BP-2, two isolates recovered from selective culturing of swabs of two Pakistani ICU room surfaces which could not be identified using MALDI-ToF mass spectroscopy. The colony morphology of strains BP-1(T) and BP-2 on blood agar is indistinguishable from non-swarming Enterobacteriaceae, but the results of growth on MacConkey plates and ChromID ESBL agar show phenotypic similarities to lactose fermenting and extended spectrum β-lactamase producing species such as E. coli and K. pneumoniae (Grohs et al., 2013). Ultrastructural analysis showed morphological and size similarity to other Enterobacteriaceae species (Feng et al., 2014). Consistent with the observed lack of a flagellum, strains BP-1(T) and BP-2 were non-motile like Yersinia, Klebsiella and Shigella genera of Enterobacteriaceae (Kramer et al., 2006).

Differential growth assays and API results were not able to assign strains BP-1(T) and BP-2 to an existing genus with high confidence. API testing determined that both isolates share a unique suite of phenotypes that distinguish BP-1(T) and BP-2 from other Enterobacteriaceae genera. The closest identity to strains BP-1(T) and BP-2 based on biochemical utilization and API algorithm predictions were Kluyvera species (41.2%), E. cloacae (30.7%) and Citrobacter koseri/farmeri (8.4%). Strains BP-1(T) and BP-2 differed from these most closely related taxa by API in arginine dihydrolase, indole production, acetoin production, sorbitol fermentation and motility.

To further probe the appropriate classification of BP-1(T) and BP-2, we performed Illumina whole genome sequencing on both isolates. As the BP-1(T) 16S rRNA gene sequence is 98.7% identical to the C. farmeri CDC 2991-81 rRNA gene, we speculated that these isolates may represent a novel species. To gain greater taxonomic resolution beyond the limitations of using just the 16S rRNA sequence, we identified and aligned all shared protein coding sequences between strains BP-1(T), BP-2 and type strains from Enterobacter, Klebsiella, Kluyvera, Citrobacter, Kosakonia, Escherichia, Shigella, Salmonella, and Siccibacter. Although the BP-1^T 16S rRNA sequence had the greatest similarity to several Citrobacter species, the unrooted NeighborNet tree and rooted maximum likelihood tree indicate BP-1(T) and BP-2 form a cluster separate from known Enterobacteriaceae species. A rooted phylogenetic tree placed the taxonomic position of BP-1(T) and BP-2 in between two large clades; one clade contains Escherichia, Shigella, Salmonella and Citrobacter and the other contains Klebsiella, Kluyvera, Enterobacter, Siccibacter and Kosakonia. ANI is the accepted in silico version of a DNA-DNA hybridization assay, the gold standard for demarcation of new bacterial species (Richter and Rossello-Mora, 2009). The 99.97% ANI between strains BP-1(T) and BP-2 indicate that they are the same species and possibly clonal isolates. A dendrogram of the ANI output reflected a tree topology like the rooted core genome tree but had Kosakonia and not Kluyvera/Klebsiella as the closest genera to BP-1(T) and BP-2.

Both strains BP-1(T) and BP-2 are multidrug resistant Enterobacteriaceae which harbor extended spectrum β-lactamases and are therefore classified as a serious antimicrobial threat by the CDC (CDC, 2013; Iredell et al., 2016). The Class A β-lactamase identified by Resfams had high identity to bla_HER family genes, which were previously described in Atlantibacter hermanii but only conferred resistance to penicillins (Beauchef-Havard et al., 2003). As this gene was proximal to a tetR response regulator, it is possible that it may be inducible, like the ampR/ampC system widespread in Enterobacter and Citrobacter (Vadlamani et al., 2015). However further work comparing the expression levels under basal or induced conditions would be needed to ascertain this. TEM-1 β-lactamase genes are widespread in E. coli and K. pneumoniae and can confer resistance to penicillins and first/second generation cephalosporins, but not 3rd generation drugs like ceftazidime (Cantu et al., 1997). bla_TEM−1b, the variant found in strains BP-1(T) and BP-2 is not considered an ESBL (Udomsantisuk et al., 2011). Therefore, a combination of blaTEM-1 and bla_SHV−12 are therefore likely the factor most critical for cephalosporin resistance (Teshager et al., 2000; Newire et al., 2013). However, efflux and porin activity may augment β-lactamase production to achieve clinical resistance, as seen in other Enterobacteriaceae species (Wozniak et al., 2012; Taherpour and Hashemi, 2013). Several predicted antibiotic resistance genes in these genomes were co-localized with IS1247, an IS1380 family transposase. A previously discovered IS1380 family member, ISEcp1 can transfer bla_CTX−M via transposition (Poirel et al., 2005; Toleman et al., 2006).

One limitation of our study is that we have only procured two strains and given the temporal similarity they may be clonally related. Although not ideal from a taxonomic perspective, the Enterobacteriaceae species Klebsiella michiganensis, Enterobacter soli, Enterobacter muelleri and genera Chania and Nissabacter were all established from single strain investigations (Manter et al., 2011; Saha et al., 2013; Kampfer et al., 2015; Ee et al., 2016; Mlaga et al., 2017). This study further demonstrates the utility of whole genome sequencing for pathogen identification in a clinical setting. Similarity, another report found that Vitek2 identified a blood culture derived isolate as E. cloacae but Illumina whole genome sequencing unequivocally determined it to be Kosakonia radicincitans (Bhatti et al., 2017). Further comprehensive work is therefore warranted to investigate clinical outcomes that result from these species misidentification.

Further work using long-read assembly such as PacBio or Oxford Nanopore is warranted to construct a high-quality complete genome. These additional efforts may also unambiguously identify plasmid components that replicate independently of the chromosome but could be collapsed into the same fasta file using our draft genome approach. This may additionally inform the presence of specific ARGs as chromosomal or plasmid-borne. Unfortunately, these techniques were not available to us during the course of this investigation.

Description of Superficieibacter gen. nov.

Superficieibacter (Su.per.fi.ci.e.i.bac′ter. L. fem. n. superficies surface; N.L. masc. n. bacter a rod: N.L. masc. n. Superficieibacter a rod from a surface).

Cells are facultative aerobic Gram-negative bacilli visualized without chains but occasionally in clusters under electron microscopy. Non-motile with no evidence of a flagella. Cells are mesophilic with robust growth at 35°C, scant growth at 25 and 32°C and no growth at 4°C. Cells can grow anaerobically at 35°C. Cells are catalase positive and can produce gas in the TSI and LIA slants. Cells can ferment arabinose, amygdalin, melibiose, sacharose, rhamnose, sorbitol, lactose, glucose and mannitol. Strains can be differed from related genera by being negative for arginine dihydrolase, indole production, acetoin production and motility but positive for sorbitol fermentation. Members belong to class Gammaproteobacteria, order Enterobacteriales and family Enterobacteriaceae. The type species of the genus is Superficieibacter electus.

Description of Superficieibacter electus sp. nov.

Superficieibacter electus (e.lec'tus. L. part. adj. electus chosen).

Main attributes are applicable from genus. Colonies form: circular, shiny, mucoid, non-haemolytic gray colonies on blood agar plates; small, blue and rigid colonies on HardyCHROM™ ESBL Agar; bright pink colonies on MacConkey agar; Coral colonies on HE agar. Cells are rod shaped and ~2.7 μm in length. 16S rRNA sequence of the strain BP-1(T) had highest percentage identity to C. farmeri CDC 2991-81(T) (98.63%) C. amalonaticus CECT 863(T) (98.56%), Citrobacter sedlakii NBRC 105722(T) (97.74%) and Citrobacter rodentium NBRC 105723(T) (97.74%). Strain BP-1(T) and BP-2 form a cluster with each other based on core genome phylogeny with related genera in the family Enterobacteriaceae. GC percentage of strain BP-1(T) and BP-2 is 52.4%. The type strain is BP-1(T) (= ATCC BAA-2937, = NBRC 113412).

Author Contributions

RP and AD wrote the manuscript and performed in silico analysis. MW and AS performed culture work, growth assays and API test. DG isolated the bacteria from hospital surfaces. SP extracted genomic DNA and prepared Illumina sequencing libraries. WB performed electron microscopy. JK, SA, C-AB and GD devised the study.

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

This work was supported in part by a United States Agency for International Development award (3220-29047) to SA, C-AB and GD. RP received support from an NIGMS training grant through award T32 GM007067 (PI: James Skeath). AD received support from The Institutional Program Unifying Population and Laboratory-Based Sciences Burroughs Welcome Fund grant to Washington University. JK received support from the Washington University Institute of Clinical and Translational Sciences grant UL1TR000448, sub-award KL2TR000450 from the National Center for Advancing Translational Sciences of the National Institutes of Health. The content is solely the responsibility of the authors and does not necessarily represent the official views of the funding agencies. We would like to thank Eric Martin, Brian Koebbe, and Jessica Hoisington-López from the Edison Family Center for Genome Sciences & Systems Biology at Washington University School of Medicine for technical support in high-throughput computing and sequencing. We would like to thank George Garrity for recommendations on proper taxonomic analysis. We would like to thank Aharon Oren for extensive advice on bacterial nomenclature.

Supplementary Material

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

References

Abbott, J. C., and Aanensen, D. M. Bentley, S. D. (2007). WebACT: an online genome comparison suite. Methods Mol. Biol. 395, 57–74. doi: 10.1007/978-1-59745-514-5_4

PubMed Abstract | CrossRef Full Text | Google Scholar

Bankevich, A., Nurk, S., Antipov, D., Gurevich, A. A., and Dvorkin, M. Pevzner, P. A. (2012). SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19, 455–477. doi: 10.1089/cmb.2012.0021

PubMed Abstract | CrossRef Full Text | Google Scholar

Baym, M., Kryazhimskiy, S., Lieberman, T. D., Chung, H., and Desai, M. M. Kishony, R. (2015). Inexpensive multiplexed library preparation for megabase-sized genomes. PLoS ONE 10:e0128036. doi: 10.1371/journal.pone.0128036

PubMed Abstract | CrossRef Full Text | Google Scholar

Beauchef-Havard, A., Arlet, G., Gautier, V., Labia, R., and Grimont, P. Philippon, A. (2003). Molecular and biochemical characterization of a novel class A beta-lactamase (HER-1) from Escherichia hermannii. Antimicrob. Agents Chemother. 47, 2669–2673. doi: 10.1128/AAC.47.8.2669-2673.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

Bhatti, M. D., Kalia, A., Sahasrabhojane, P., Kim, J., and Greenberg, D. E. Shelburne, S. A. (2017). Identification and whole genome sequencing of the first case of kosakonia radicincitans causing a human bloodstream infection. Front. Microbiol. 8:62. doi: 10.3389/fmicb.2017.00062

PubMed Abstract | CrossRef Full Text | Google Scholar

Bolger, A. M., and Lohse, M. Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120. doi: 10.1093/bioinformatics/btu170

PubMed Abstract | CrossRef Full Text | Google Scholar

Brady, C., Cleenwerck, I., Venter, S., and Coutinho, T. De Vos, P. (2013). Taxonomic evaluation of the genus Enterobacter based on multilocus sequence analysis (MLSA): proposal to reclassify E. nimipressuralis and E. amnigenus into Lelliottia gen. nov. as Lelliottia nimipressuralis comb. nov. and Lelliottia amnigena comb. nov., respectively, E. gergoviae and E. pyrinus into Pluralibacter gen. nov. as Pluralibacter gergoviae comb. nov. and Pluralibacter pyrinus comb. nov., respectively, E. cowanii, E. radicincitans, E. oryzae and E. arachidis into Kosakonia gen. nov. as Kosakonia cowanii comb. nov., Kosakonia radicincitans comb. nov., Kosakonia oryzae comb. nov. and Kosakonia arachidis comb. nov., respectively and E. turicensis, E. helveticus and E. pulveris into Cronobacter as Cronobacter zurichensis nom. nov., Cronobacter helveticus comb. nov. and Cronobacter pulveris comb. nov., respectively and emended description of the genera Enterobacter and Cronobacter. Syst. Appl. Microbiol. 36, 309–319. doi: 10.1016/j.syapm.2013.03.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Brenner, D. J., Grimont, P. A. D., Steigerwalt, A. G., Fanning, G. R., Ageron, E., and Riddle, C. F. (1993). Classification of citrobacteria by DNA hybridization: designation of Citrobacter farmeri sp. nov., Citrobacter youngae sp. nov., Citrobacter braakii sp. nov., Citrobacter werkmanii sp. nov., Citrobacter sedlakii sp. nov., and three unnamed Citrobacter genomospecies. Int. J. Syst. Bacteriol. 43, 645–658. doi: 10.1099/00207713-43-4-645

PubMed Abstract | CrossRef Full Text | Google Scholar

Camacho, C., Coulouris, G., Avagyan, V., Ma, N., Papadopoulos, J., Bealer, K. L., et al. (2009). BLAST+: architecture and applications. BMC Bioinformatics 10:421. doi: 10.1186/1471-2105-10-421

CrossRef Full Text

Cantu, C. III., and Huang, W. and Palzkill, T. (1997). Cephalosporin substrate specificity determinants of TEM-1 beta-lactamase. J. Biol. Chem. 272, 29144–29150. doi: 10.1074/jbc.272.46.29144

PubMed Abstract | CrossRef Full Text | Google Scholar

CDC (2013). Antibiotic Resistance Threats in the United States. CDC.

PubMed Abstract

Chakravorty, S., Helb, D., Burday, M., and Connell, N. Alland, D. (2007). A detailed analysis of 16S ribosomal RNA gene segments for the diagnosis of pathogenic bacteria. J. Microbiol. Methods 69, 330–339. doi: 10.1016/j.mimet.2007.02.005

PubMed Abstract | CrossRef Full Text | Google Scholar

CLSI (2013). Performance Standards for Antimicrobial Susceptibility Testing: Twenty-Third Informational Supplment M100-S23. Wayne, PA: Clinical and Laboratory Standards Institute.

Edgar, R. C. (2004). MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics 5:113. doi: 10.1186/1471-2105-5-113

PubMed Abstract | CrossRef Full Text | Google Scholar

Ee, R., Madhaiyan, M., Ji, L., Lim, Y. L., Nor, N. M., Tee, K. K., et al. (2016). Chania multitudinisentens gen. nov., sp. nov., an N-acyl-homoserine-lactone-producing bacterium in the family Enterobacteriaceae isolated from landfill site soil. Int. J. Syst. Evol. Microbiol. 66, 2297–2304. doi: 10.1099/ijsem.0.001025

PubMed Abstract | CrossRef Full Text | Google Scholar

Farmer, J. J. III. (1981). “The genus citrobacter,” in The Prokaryotes, eds M. P. Starr, H. Stolp, H. G. Trüper, A. Balows, and H. G. Schlegel (Berlin: Springer-Verlag), 1140–1147.

Google Scholar

Farmer, J. J. III., Davis, B. R., Hickman-Brenner, F. W., McWhorter, A., Huntley-Carter, G. P., Asbury, M. A., et al. (1985). Biochemical identification of new species and biogroups of Enterobacteriaceae isolated from clinical specimens. J. Clin. Microbiol. 21, 46–76.

PubMed Abstract | Google Scholar

Farmer, J. J. III., Fanning, G. R., Huntley-Carter, G. P., Holmes, B., Hickman, F. W., and Richard, C. Brenner, D. J. (1981). Kluyvera, a new (redefined) genus in the family Enterobacteriaceae: identification of Kluyvera ascorbata sp. nov. and Kluyvera cryocrescens sp. nov. in clinical specimens. J. Clin. Microbiol. 13, 919–933.

PubMed Abstract | Google Scholar

Feng, C., Li, J., Qin, D., Chen, L., Zhao, F., Chen, S. P., et al. (2014). Characterization of exoelectrogenic bacteria enterobacter strains isolated from a microbial fuel cell exposed to copper shock load. PLoS ONE 9:e113379. doi: 10.1371/journal.pone.0113379

PubMed Abstract | CrossRef Full Text | Google Scholar

Gibson, M. K., and Forsberg, K. J. Dantas, G. (2015). Improved annotation of antibiotic resistance determinants reveals microbial resistomes cluster by ecology. ISME J. 9, 207–216. doi: 10.1038/ismej.2014.106

PubMed Abstract | CrossRef Full Text | Google Scholar

Grohs, P., Tillecovidin, B., Caumont-Prim, A., Carbonnelle, E., Day, N., Podglajen, I., et al. (2013). Comparison of five media for detection of extended-spectrum Beta-lactamase by use of the wasp instrument for automated specimen processing. J. Clin. Microbiol. 51, 2713–2716. doi: 10.1128/JCM.00077-13

PubMed Abstract | CrossRef Full Text

Gu, C. T., Li, C. Y., Yang, L. J., and Huo, G. C. (2014). Enterobacter xiangfangensis sp. nov., isolated from Chinese traditional sourdough, and reclassification of Enterobacter sacchari Zhu et al. 2013 as Kosakonia sacchari comb. nov. Int. J. Syst. Evol. Microbiol. 64(Pt 8), 2650–2656. doi: 10.1099/ijs.0.064709-0

PubMed Abstract | CrossRef Full Text

Gurevich, A., Saveliev, V., and Vyahhi, N. Tesler, G. (2013). QUAST: quality assessment tool for genome assemblies. Bioinformatics 29, 1072–1075. doi: 10.1093/bioinformatics/btt086

PubMed Abstract | CrossRef Full Text | Google Scholar

Hata, H., Natori, T., Mizuno, T., Kanazawa, I., Eldesouky, I., Hayashi, M., et al. (2016). Phylogenetics of family Enterobacteriaceae and proposal to reclassify Escherichia hermannii and Salmonella subterranea as Atlantibacter hermannii and Atlantibacter subterranea gen. nov., comb. nov. Microbiol. Immunol. 60, 303–311. doi: 10.1111/1348-0421.12374

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoffmann, H., Stindl, S., Ludwig, W., Stumpf, A., Mehlen, A., Monget, D., et al. (2005). Enterobacter hormaechei subsp. oharae subsp. nov., E. hormaechei subsp. Hormaechei comb. nov., and E. hormaechei subsp. Steigerwaltii subsp. nov., three new subspecies of clinical importance. J. Clin. Microbiol. 43, 3297–3303. doi: 10.1128/JCM.43.7.3297-3303.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Holt, K. E., Wertheim, H., Zadoks, R. N., Baker, S., Whitehouse, C. A., Dance, D., et al. (2015). Genomic analysis of diversity, population structure, virulence and antimicrobial resistance in Klebsiella pneumoniae, an urgent threat to public health. Proc. Natl. Acad. Sci. U.S.A. 112, E3574–E3581. doi: 10.1073/pnas.1501049112

PubMed Abstract | CrossRef Full Text | Google Scholar

Horcajada, J. P., Shaw, E., Padilla, B., Pintado, V., Calbo, E., Benito, N., et al. (2013). Healthcare-associated, community-acquired and hospital-acquired bacteraemic urinary tract infections in hospitalized patients: a prospective multicentre cohort study in the era of antimicrobial resistance. Clin. Microbiol. Infect. 19, 962–968. doi: 10.1111/1469-0691.12089

PubMed Abstract | CrossRef Full Text | Google Scholar

Huson, D. H. and Bryant, D. (2006). Application of phylogenetic networks in evolutionary studies. Mol. Biol. Evol. 23, 254–267. doi: 10.1093/molbev/msj030

PubMed Abstract | CrossRef Full Text | Google Scholar

Huys, G., Cnockaert, M., and Janda, J. M. Swings, J. (2003). Escherichia albertii sp. nov., a diarrhoeagenic species isolated from stool specimens of Bangladeshi children. Int. J. Syst. Evol. Microbiol. 53, 807–810. doi: 10.1099/ijs.0.02475-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Iredell, J., and Brown, J. Tagg, K. (2016). Antibiotic resistance in Enterobacteriaceae: mechanisms and clinical implications. BMJ 352:h6420. doi: 10.1136/bmj.h6420

PubMed Abstract | CrossRef Full Text | Google Scholar

Kämpfer, P., and McInroy, J. A. Glaeser, S. P. (2015). Enterobacter muelleri sp. nov., isolated from the rhizosphere of zea mays. Int. J. Syst. Evol. Microbiol. 65, 4093–4099. doi: 10.1099/ijsem.0.000547

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, M., Oh, H. S., and Park, S. C. Chun, J. (2014). Toward a taxonomic coherence between average nucleotide identity and 16S rRNA gene sequence similarity for species demarcation of prokaryotes. Int. J. Syst. Evol. Microbiol. 64, 346–351. doi: 10.1099/ijs.0.059774-0

CrossRef Full Text | Google Scholar

Kramer, A., and Schwebke, I. Kampf, G. (2006). How long do nosocomial pathogens persist on inanimate surfaces? A systematic review. BMC Infect. Dis. 6:130. doi: 10.1186/1471-2334-6-130

CrossRef Full Text

Manter, D. K., and Hunter, W. J. Vivanco, J. M. (2011). Enterobacter soli sp. nov.: a lignin-degrading gamma-proteobacteria isolated from soil. Curr. Microbiol. 62, 1044–1049. doi: 10.1007/s00284-010-9809-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Mlaga, K. D., Lotte, R., Montaudie, H., and Rolain, J. M. Ruimy, R. (2017). ‘Nissabacter archeti' gen. nov., sp. nov., a new member of Enterobacteriaceae family, isolated from human sample at Archet 2 Hospital, Nice, France. New Microbes New Infect. 17, 81–83. doi: 10.1016/j.nmni.2017.02.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Newire, E. A., Ahmed, S. F., House, B., and Valiente, E. Pimentel, G. (2013). Detection of new SHV-12, SHV-5 and SHV-2a variants of extended spectrum beta-lactamase in Klebsiella pneumoniae in Egypt. Ann. Clin. Microbiol. Antimicrob. 12:16. doi: 10.1186/1476-0711-12-16

CrossRef Full Text | Google Scholar

Oberauner, L., Zachow, C., Lackner, S., Hogenauer, C., Smolle, K., and H. Berg, G. (2013). The ignored diversity: complex bacterial communities in intensive care units revealed by 16S pyrosequencing. Sci. Rep. 3:1413. doi: 10.1038/srep01413

PubMed Abstract | CrossRef Full Text | Google Scholar

Page, A. J., Cummins, C. A., Hunt, M., Wong, V. K., Reuter, S., Holden, M. T., et al. (2015). Roary: rapid large-scale prokaryote pan genome analysis. Bioinformatics 31, 3691–3693. doi: 10.1093/bioinformatics/btv421

PubMed Abstract | CrossRef Full Text | Google Scholar

Paradis, S., Boissinot, M., Paquette, N., Belanger, S. D., Martel, E. A., Boudreau, D. K., et al. (2005). Phylogeny of the Enterobacteriaceae based on genes encoding elongation factor Tu and F-ATPase beta-subunit. Int. J. Syst. Evol. Microbiol. 55, 2013–2025. doi: 10.1099/ijs.0.63539-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Poirel, L., Lartigue, M. F., and Decousser, J. W. Nordmann, P. (2005). ISEcp1B-mediated transposition of blaCTX-M in Escherichia coli. Antimicrob. Agents Chemother. 49, 447–450. doi: 10.1128/AAC.49.1.447-450.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Potter, R. F., and D′Souza, A. W. Dantas, G. (2016). The rapid spread of carbapenem-resistant Enterobacteriaceae. Drug Resist. Updat. 29, 30–46. doi: 10.1016/j.drup.2016.09.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Ribeiro, T. G., Clermont, D., Branquinho, R., Machado, E., Peixe, L., and Brisse, S. (2017). Citrobacter europaeus sp. nov., isolated from water and human faecal samples. Int. J. Syst. Evol. Microbiol. 67, 170–173. doi: 10.1099/ijsem.0.001606

PubMed Abstract | CrossRef Full Text | Google Scholar

Richter, M. and Rosselló-Móra, R. (2009). Shifting the genomic gold standard for the prokaryotic species definition. Proc. Natl. Acad. Sci. U.S.A. 106, 19126–19131. doi: 10.1073/pnas.0906412106

PubMed Abstract | CrossRef Full Text | Google Scholar

Richter, M., Rossello-Mora, R., and Oliver Glockner, F. 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

Russotto, V., Cortegiani, A., and Raineri, S. M. Giarratano, A. (2015). Bacterial contamination of inanimate surfaces and equipment in the intensive care unit. J. Intensive Care 3:54. doi: 10.1186/s40560-015-0120-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Saha, R., Farrance, C. E., Verghese, B., and Hong, S. Donofrio, R. S. (2013). Klebsiella michiganensis sp. nov., a new bacterium isolated from a tooth brush holder. Curr. Microbiol. 66, 72–78. doi: 10.1007/s00284-012-0245-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Schauer, D. B., Zabel, B. A., Pedraza, I. F., O'Hara, C. M., Steigerwalt, A. G., and Brenner, D. J. (1995). Genetic and biochemical characterization of Citrobacter rodentium sp. nov. J. Clin. Microbiol. 33, 2064–2068.

PubMed Abstract | Google Scholar

Schmieder, R. and Edwards, R. (2011). Fast identification and removal of sequence contamination from genomic and metagenomic datasets. PLoS ONE 6:e17288. doi: 10.1371/journal.pone.0017288

PubMed Abstract | CrossRef Full Text | Google Scholar

Seemann, T. (2014). Prokka: rapid prokaryotic genome annotation. Bioinformatics 30, 2068–2069. doi: 10.1093/bioinformatics/btu153

PubMed Abstract | CrossRef Full Text | Google Scholar

Sirot, J., Chanal, C., Petit, A., Sirot, D., and Labia, R. Gerbaud, G. (1988). Klebsiella pneumoniae and other Enterobacteriaceae producing novel plasmid-mediated beta-lactamases markedly active against third-generation cephalosporins: epidemiologic studies. Rev. Infect. Dis. 10, 850–859. doi: 10.1093/clinids/10.4.850

PubMed Abstract | CrossRef Full Text | Google Scholar

Stamatakis, A. (2014). RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30, 1312–1313. doi: 10.1093/bioinformatics/btu033

PubMed Abstract | CrossRef Full Text | Google Scholar

Sullivan, M. J., and Petty, N. K. Beatson, S. A. (2011). Easyfig: a genome comparison visualizer. Bioinformatics 27, 1009–1010. doi: 10.1093/bioinformatics/btr039

PubMed Abstract | CrossRef Full Text | Google Scholar

Taherpour, A. and Hashemi, A. (2013). Detection of OqxAB efflux pumps, OmpK35 and OmpK36 porins in extended-spectrum-beta-lactamase-producing Klebsiella pneumoniae isolates from Iran. Hippokratia 17, 355–358.

PubMed Abstract | Google Scholar

Teshager, T., Domínguez, L., Moreno, M. A., Saenz, Y., and Torres, C. Cardenosa, S. (2000). Isolation of an SHV-12 beta-lactamase-producing Escherichia coli strain from a dog with recurrent urinary tract infections. Antimicrob. Agents Chemother. 44, 3483–3484. doi: 10.1128/AAC.44.12.3483-3484.2000

PubMed Abstract | CrossRef Full Text | Google Scholar

Thompson, C. C., Chimetto, L., Edwards, R. A., Swings, J., and Stackebrandt, E. Thompson, F. L. (2013). Microbial genomic taxonomy. BMC Genomics 14:913. doi: 10.1186/1471-2164-14-913

CrossRef Full Text | Google Scholar

Tindall, B. J., Grimont, P. A. D., Garrity, G. M., and Euzéby, J. P. (2005). Nomenclature and taxonomy of the genus Salmonella. Int. J. Syst. Evol. Microbiol. 55(Pt 1), 521–524. doi: 10.1099/ijs.0.63580-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Toleman, M. A., and Bennett, P. M. Walsh, T. R. (2006). ISCR elements: novel gene-capturing systems of the 21st century? Microbiol. Mol. Biol. Rev. 70, 296–316. doi: 10.1128/MMBR.00048-05

PubMed Abstract | CrossRef Full Text | Google Scholar

Udomsantisuk, N., Nunthapisud, P., and Tirawatanapong, T. Dansuputra, M. (2011). Molecular characterization of extended spectrum beta-lactamase among clinical isolates Escherichia coli and Klebsiella pneumoniae. J. Med. Assoc. Thai. 94, 1504–1512.

PubMed Abstract | Google Scholar

Vadlamani, G., Thomas, M. D., Patel, T. R., Donald, L. J., Reeve, T. M., Stetefeld, J., et al. (2015). The beta-lactamase gene regulator AmpR is a tetramer that recognizes and binds the D-Ala-D-Ala motif of its repressor UDP-N-acetylmuramic acid (MurNAc)-pentapeptide. J. Biol. Chem. 290, 2630–2643. doi: 10.1074/jbc.M114.618199

PubMed Abstract | CrossRef Full Text

Walk, S. T., Alm, E. W., Gordon, D. M., Ram, J. L., Toranzos, G. A., Tiedje, J. M., et al. (2009). Cryptic lineages of the genus Escherichia. Appl. Environ. Microbiol. 75, 6534–6544. doi: 10.1128/AEM.01262-09

PubMed Abstract | CrossRef Full Text | Google Scholar

Waterhouse, A. M., Procter, J. B., Martin, D. M., and Clamp, M. Barton, G. J. (2009). Jalview Version 2–a multiple sequence alignment editor and analysis workbench. Bioinformatics 25, 1189–1191. doi: 10.1093/bioinformatics/btp033

PubMed Abstract | CrossRef Full Text | Google Scholar

Wozniak, A., Villagra, N. A., Undabarrena, A., Gallardo, N., Keller, N., Moraga, M., et al. (2012). Porin alterations present in non-carbapenemase-producing Enterobacteriaceae with high and intermediate levels of carbapenem resistance in Chile. J. Med. Microbiol. 61, 1270–1279. doi: 10.1099/jmm.0.045799-0

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

Young, V. M., Kenton, D. M., Hobbs, B. J., and Moody, M. R. (1971). Levinea, a new genus of the family Enterobacteriaceae. Int. J. Syst. Bacteriol. 21, 58–63. doi: 10.1099/00207713-21-1-58

CrossRef Full Text | Google Scholar

Zankari, E., Hasman, H., Cosentino, S., Vestergaard, M., Rasmussen, S., Lund, O., et al. (2012). Identification of acquired antimicrobial resistance genes. J. Antimicrob. Chemother. 67, 2640–2644. doi: 10.1093/jac/dks261

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, B., Zhou, Q., Lin, L., Hu, C., Shen, P., Yang, L., et al. (2013). Enterobacter sacchari sp. nov., a nitrogen-fixing bacterium associated with sugar cane (Saccharum officinarum L.). Int. J. Syst. Evol. Microbiol. 63(Pt 7), 2577–2582. doi: 10.1099/ijs.0.045500-0

PubMed Abstract | CrossRef Full Text | Google Scholar