#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Comparative genomic analyses reveal diverse virulence factors and antimicrobial resistance mechanisms in clinical Elizabethkingia meningoseptica strains


Authors: Shicheng Chen aff001;  Marty Soehnlen aff002;  Jochen Blom aff003;  Nicolas Terrapon aff004;  Bernard Henrissat aff004;  Edward D. Walker aff001
Authors place of work: Department of Microbiology and Molecular Genetics, Michigan State University, East Lansing, MI, United States of America aff001;  Michigan Department of Health and Human Services, Bureau of Laboratories, Lansing, MI, United States of America aff002;  Bioinformatics and Systems Biology, Justus-Liebig-University, Giessen, Germany aff003;  Architecture et Fonction des Macromolécules Biologiques, Centre National de la Recherche Scientifique (CNRS), Aix-Marseille Université (AMU), UMR 7257, Marseille, France aff004;  Institut National de la Recherche Agronomique (INRA), USC 1408 AFMB, Marseille, France aff005;  Department of Biological Sciences, King Abdulaziz University, Jeddah, Saudi Arabia aff006
Published in the journal: PLoS ONE 14(10)
Category: Research Article
doi: https://doi.org/10.1371/journal.pone.0222648

Summary

Three human clinical isolates of bacteria (designated strains Em1, Em2 and Em3) had high average nucleotide identity (ANI) to Elizabethkingia meningoseptica. Their genome sizes (3.89, 4.04 and 4.04 Mb) were comparable to those of other Elizabethkingia species and strains, and exhibited open pan-genome characteristics, with two strains being nearly identical and the third divergent. These strains were susceptible only to trimethoprim/sulfamethoxazole and ciprofloxacin amongst 16 antibiotics in minimum inhibitory tests. The resistome exhibited a high diversity of resistance genes, including 5 different lactamase- and 18 efflux protein- encoding genes. Forty-four genes encoding virulence factors were conserved among the strains. Sialic acid transporters and curli synthesis genes were well conserved in E. meningoseptica but absent in E. anophelis and E. miricola. E. meningoseptica carried several genes contributing to biofilm formation. 58 glycoside hydrolases (GH) and 25 putative polysaccharide utilization loci (PULs) were found. The strains carried numerous genes encoding two-component system proteins (56), transcription factor proteins (187~191), and DNA-binding proteins (6~7). Several prophages and CRISPR/Cas elements were uniquely present in the genomes.

Keywords:

Genome analysis – DNA-binding proteins – Comparative genomics – Virulence factors – antibiotics – Antimicrobial resistance – Antibiotic resistance – Sialic acids

Introduction

Elizabethkingia meningoseptica is a Gram-negative, non-fermenting, aerobic bacterium occurring in soil, water, plants and animals [1]. E. meningoseptica normally does not cause infection and disease in healthy humans, but it is a serious causative agent of nosocomial pneumonia, neonatal meningitis, bacteremia, and endocarditis in immuno-compromised patient populations [24]. Infections are mostly acquired in hospitals as bacteria are often recovered from medical apparatus and reagents including tap water, disinfection fluid, ventilators, hemodialysis equipment and catheters; however, sporadic, community-acquired infections were also reported [36]. Direct transmission pathways for E. meningoseptica remain largely unknown [7].

Clinical manifestations of E. meningoseptica infections are similar to those of other Elizabethkingia species or other bacteria though some variations exist [8, 9]. Elizabethkingia infections in neonates or immuno-compromised patients are usually associated with a poor outcome [2, 3]. One of the most challenging issues is that routine morphological, biochemical, and molecular tests cannot accurately identify Elizabethkingia species [10]. Sequencing of the 16S rRNA gene often fails to provide sufficient resolution to differentiate various Elizabethkingia species which may show different antibiotic resistance properties [10, 11]. For example, recent studies have demonstrated that E. anophelis or E. meningoseptica were frequently misidentified; it is not surprising that some E. meningoseptica infections were actually reported to be caused by E. anophelis [12]. Thus, accurate and complementary diagnostic methods need to be developed such as MALDI-ToF or genome sequencing.

Elizabethkingia are typically highly resistant to antimicrobials including extended-spectrum beta-lactams, tetracycline, aminoglycosides, and chloramphenicol [13]. Nevertheless, E. meningoseptica-infections are empirically treated with antibiotics used for Gram-positive bacteria [14]. A growing number of studies demonstrate that Elizabethkingia species (even the same species) isolated from different geographical regions have different antibiotic susceptibilities, showing that there is a complex antimicrobial resistance spectrum in Elizabethkingia [3, 13, 15]. Consequently, the high mortality rates (up to 50%) reported for immunosuppressed patients may, at least partially, be the result of a delay in identifying the appropriate antibiotic therapy [3, 10, 13]. Thanks to next generation sequencing technology, the studies on the antimicrobial resistance mechanisms have significantly progressed. Recently, comparative genome analysis in Elizabethkingia has been conducted to discover the breadth of antibiotic susceptibility and epidemiological features [1618]. However, most of the studies are limited to E. anophelis and focused on clinical cases, antimicrobial susceptibility patterns, or characterization of epidemic outbreaks [1618]. Genome analyses and physiological studies have not been conducted in depth for E. meningoseptica.

The aim of this study was to analyze virulence, antibiotic resistance, environmental survival, and adaption mechanisms through genomic, resistome, and antibiotic susceptibility analyses. With the completion of genome sequencing, assembly, and annotation of E. meningoseptica strains, we examined gene repertoire and genetic diversity in comparison with other Elizabethkingia species. Further, the regulatory systems, sugar utilization systems, virulence factors and antibiotic resistance genes were analyzed and compared to those in the selected Elizabethkingia. Collectively, our study offers the opportunity to explore their virulence, antibiotic resistance, environmental survival, and adaptation mechanisms.

Materials and methods

Strain background and culture conditions

E. meningoseptica strains, designated here Em1, Em2 and Em3, were isolated from patients in Michigan between 2015 and 2016 when the Elizabethkingia outbreak occurred in the Midwest regions (Table 1). The three strains were chosen for studying because they were reported to Michigan Department of Health and Human Services. E. meningoseptica Em1 was isolated from tracheal secretions in a female patient on March 6th, 2016. Em2 was isolated from whole blood in a male patient on November 10th, 2015. Em3 was isolated from sputum in a male patient on February 6th, 2016. The three separate patients living in different counties in Michigan visited different clinics during the treatments. The representative E. anophelis and E. miricola strains were selected in this study due to diverse geography and isolation origin. The genome E. meningoseptica strains were grown aerobically in tryptic soy broth (TSB) broth at 30°C. Bacto agar (Difco, Detroit, MI) was added to a final concentration of 20 g/liter. The sheep blood agar (SBA) was purchased from Thermo Scientific (Waltham, MA).

Tab. 1. General features of selected Elizabethkingia genomes.
General features of selected <i>Elizabethkingia</i> genomes.

Antimicrobial susceptibility testing (AST)

Cultures were freshly grown on SBA overnight at 37°C and cell suspensions adjusted in saline solution (0.9%) to a turbidity equivalent of 0.5 McFarland standard. The VITEK 2 microbial ID/AST testing system (version 07.01) with the GNP 70 antibiotic susceptibility cards (bioMerieux) was used to determine the minimum inhibitory concentration (MIC) and classification into resistance phenotypes. MIC results were interpreted according to Clinical and Laboratory Standards Institute criteria (CLSI) [19].

Genomic DNA preparation, sequencing and assembly

DNA was isolated using a Wizard Genomic DNA Purification Kit (Promega, Madison). The concentration of genomic DNA was measured using a Nanodrop2000 UV-Vis Spectrophotometer (Thermo scientific) and Qubit DNA assay kit. DNA integrity was evaluated by agarose gel assay (1.5%, w/v).

NGS libraries were prepared using the Illumina TruSeq Nano DNA Library Preparation Kit. Completed libraries were evaluated using a combination of Qubit dsDNA HS, Caliper LabChipGX HS DNA and Kapa Illumina Library Quantification qPCR assays. Libraries were combined in a single pool for multiplex sequencing and the pool was loaded on one standard MiSeq flow cell (v2) and sequencing performed in a 2x250 bp, paired end format using a v2, 500 cycle reagent cartridge. Base calling was done by Illumina Real Time Analysis (RTA) v1.18.54 and output of RTA was demultiplexed and converted to FastQ format with Illumina Bcl2fastq v1.8.4. The genomes were assembled into contiguous sequences using SPAdes version 3.9 following the manual as described previously[20], then short and low-coverage contigs were filtered out.

Genome annotation

Annotation of the assembled genome sequences for Em1, Em2 and Em3 was submitted to the NCBI Prokaryotic Genome Automatic Annotation Pipeline (PGAAP). The predicted CDSs were translated and analyzed against the NCBI non-redundant database, Pfam, TIGRfam, InterPro, KEGG and COG. Additional gene prediction and manual revision was performed by using the Integrated Microbial Genomes and Microbiome Samples Expert Review (IMG/MER) platform.

Bioinformatics

Functional categorization and classification for predicted ORFs were performed by RAST server-based SEED viewer [21]. Multi-drug resistance genes were predicted in the Comprehensive Antibiotic Resistance Database [22]. Prophage prediction was done with PHAST [23] and Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) were predicted by CRISPRfinder [24]. For genome similarity assessment, average nucleotide identity (ANI) was computed using ANI calculator (https://www.ezbiocloud.net/tools/ani).

The pan genome, core genome, singletons and specific genes of Em1, Em2 and Em3 were characterized by comparing the representative Elizabethkingia genomes using EDGAR 2.0 [25]. The development of pan genome and core genome sizes was approximated using the methods proposed by Tettelin et al [26]. The core genome calculated by EDGAR 2.0 was the reference used to infer a phylogeny. The 27,144 amino acid sequences (2,088 per genome) of the core genome were aligned set-wise using MUSCLE v3.8.31 [27], resulting in a large multiple alignment with 9,015,188 amino acid residues in total (693,476 per genome). This large alignment was used to construct a phylogenetic tree using the neighbor-joining method as implemented in the PHYLIP package [28].

The regulatory elements were predicted using p2rp with default settings (http://www.p2rp.org). Bacterial protein localization prediction was conducted with tool PSORTb version 3.0.2 (http://www.psort.org/psortb/). Carbohydrate active enzyme families, including enzymes of glycan assembly (glycosyltransferases, GT) and deconstruction (glycoside hydrolases, GH, polysaccharide lyases, PL, carbohydrate esterases, CE), were semi-manually annotated using the Carbohydrate Active Enzyme (CAZy) database curation pipelines [29]. Polysaccharide-utilization loci (PUL) was predicted as described in the PULDB database (www.cazy.org/PULDB/).

Biofilm formation

Biofilm formation was evaluated using a standard assay with crystal violet staining as previously described, [30] with modifications as follows. Em1, Em2 and Em3 were grown in TSB broth to the mid-exponential phase. The cultures were diluted in TSB broth, and 100 μL were deposited in wells of 96-well microtiter polystyrene plates with flat bottoms. The negative controls were TSB medium without inoculum. The plate was incubated at 28°C under static condition for 24 h. The supernatants were discarded, the wells were washed twice with 200 μL of sterile distilled water and then 150 μL of 1% (w/v) crystal violet was added to each well. After 30 min, excess stain was removed by washing the wells four times with 200 μL of sterile distilled water and the stain bound to adherent cells was subsequently released by adding 100 μL of absolute ethanol. The biofilm formation was determined by measuring the OD595 nm using the plate reader.

Accession of the genome sequences

The data from these Whole Genome Shotgun projects have been deposited at DDBJ/ENA/GenBank under accession numbers MCJH00000000, MDTZ00000000 and MDTY00000000 for Em1, Em2 and Em3, respectively. The BioProject designations for this project are PRJNA336273, PRJNA339645 and PRJNA338129, and BioSample accession numbers are SAMN05507161, SAMN05601445 and SAMN05521514 for Em1, Em2 and Em3, respectively.

Statistical analyses

Statistical analyses were performed using SAS (version 9.2; SAS Institute, Cary, NC).

Results

Genome features and phylogenetic inferences

The assembly of strain Em1, Em2 and Em3 contained 10, 14 and 11 contigs with a size of 4.04, 3.89 and 4.04 Mbp, respectively (Table 1), which is congruent with the genome size (ranging from 3.84 to 4.04 Mbp) among the selected E. meningoseptica strains. The Em1, Em2 and Em3 genomes included 3,656, 3,483 and 3,656 coding sequences (CDS) and 55, 53 and 55 RNA genes, respectively. The GC content in the three E. meningoseptica strains (Table 1) was ca. 36.4%, consistent with other E. meningoseptica isolates (36.2~36.6%); collectively, the average GC content in E. meningoseptica genomes (36.4%, n = 26) was slightly higher than that in E. anophelis (average 35.6%, n = 64) or E. miricola (35.9%, n = 17) (S1 Fig). The gene rearrangement comparisions among the Em1, Em2 and Em3 were shown in the synteny plots (S2 Fig).

ANI values indicated that the three isolates belonged to E. meningoseptica species as they were more than 99% identical to the type strain E. meningoseptica ATCC 13253 (Table 2). However, ANI values were low (<80%) in comparisons with E. anophelis or E. miricola (Table 2). The difference in overall genome features between Em1 and Em3 was negligible (Table 1). Phylogenetic trees (S3 Fig) showed that Em1 and Em3 were clustered more closely to E. meningoseptica CSID 30005163359 and G58-80. However, the phylogenetic placement of Em2 was closer to E. meningoseptica CSID 3000516465, which formed a different clade from the other two strains (S3 Fig).

Tab. 2. Average nucleotide identity dendrogram for the selected Elizabethkingia spp*.
Average nucleotide identity dendrogram for the selected <i>Elizabethkingia</i> spp<em class="ref">*</em>.

Gene repertoire of E. meningoseptica

The core genome and pan-genome were sorted and used for gene repertoire analysis in selected E. meningoseptica genomes (Fig 1A and 1B). Core genome analysis showed that the number of shared genes decreased with addition of the input genomes (see Fig 1B). Overall, E. meningoseptica displayed an open pan-genome feature because new genes appeared when more sequenced genomes were added to the analysis (Fig 1A). E. meningoseptica Em1 and Em3 shared at least 3556 genes; only 3 and 5 genes were uniquely present in strains Em1 and Em3, respectively, supporting that Em1 and Em2 are very similar (S4 Fig). Isolate Em2 shared 3280 and 3278 genes with Em1 and Em3, respectively, while there were 125 and 127 unique genes in Em2 (S4 Fig). Moreover, E. meningoseptica Em2 shared 3,184, 3,326, 3,106, 3,322, 3,173 and 3,788 common genes with strains G4076, G4120, NBRC 12535, and 61421 PRCM, respectively (Fig 2A). These genes shared in common accounted for approximately 93.5%, 97.7%, 91.2%, and 97.5% of the encoding genes of Em2, respectively; taken together, Em2 and these four strains shared 3,076 genes (Fig 2A). However, E. meningoseptica Em2 shared far fewer genes with E. anophelis strains (2,772 to 2,823) and E. miricola (2,789 to 2,890), accounting for less than 86% of the total encoding genes (Fig 2B).

Pan, core, and singleton genome evolution according to the number of selected <i>Elizabethkingia</i> genomes.
Fig. 1. Pan, core, and singleton genome evolution according to the number of selected Elizabethkingia genomes.
(A) Number of genes (pan-genome) for a given number of genomes sequentially added. The pan development plot was generated for the following genomes: E. meningoseptica EM2 (NZ_MDTZ01000014), E. meningoseptica NV2016 (NZ_FRFB01000021), E. meningoseptica G4120 (NZ_CP016378), E. meningoseptica CSID3000516359 (NZ_MAHC01000017), E. meningoseptica CCUG214 (NZ_FLSV01000010), E. meningoseptica CSID_3000515919 (NZ_MAGZ01000024), E. meningoseptica CSID_3000516535 (NZ_MAHF01000020), E. meningoseptica ATCC13253 (NBRC_12535), E. meningoseptica EM3 (NZ_MDTY01000011), E. meningoseptica CIP111048 (NZ_FTPF01000022), E. meningoseptica 61421PRCM (NZ_MPOG01000010), E. meningoseptica EM1 (NZ_MCJH01000010), E. meningoseptica 58_80 (NZ_FTRA01000043). (B) Number of shared genes (core genome) as a function of the number of genomes sequentially added. The genomes used for generating the core genome development plot were the same as listed in (A).
Venn diagram of shared and unique genes in selected <i>Elizabethkingia</i>.
Fig. 2. Venn diagram of shared and unique genes in selected Elizabethkingia.
The unique and shared genome among the selected strains was determined by a dual cutoff of 30% or greater amino acid identity and sequence length coverage of more than 70%. EDGAR was used for Venn diagrams. A) 1: E. anophelis CSID_3015183678, 2: E. meningoseptica Em2 and 3: E. miricola BM10. B) 1: E. meningoseptica ATCC 13253, 2: E. meningoseptica 61421 PRCM, 3: E. meningoseptica Em2, E. meningoseptica G4076, and E. meningoseptica G4120.

Antimicrobial susceptibility and antibiotic resistance gene analysis

Susceptibilities of the E. meningoseptica isolates and ATCC13253 with corresponding MICs to tested antibiotics showed very similar antibiotic susceptibility spectra with minor discrepancies across strains (Table 3). The strains were highly resistant to 13 of 16 antimicrobial reagents, showing that they were multi-drug resistance strains (Table 3).

Tab. 3. Antibiotic susceptibility test (AST) in the selected E. menigoseptica.
Antibiotic susceptibility test (AST) in the selected <i>E</i>. <i>menigoseptica</i>.

Twenty-two genes encoding enzymes/proteins conferring resistance to antimicrobial reagents were found by CARD and RAST SEED subsystem (Table 4). Up to 5 lactamase genes encoding β-lactamases (EC 3.5.2.6), metal-dependent hydrolases (superfamily I), class C β-lactamases/penicillin binding protein and other β-lactamases were predicted in Em2, which possibly contributed to the intrinsic resistance to six β-lactam drugs tested in this study (Table 3). One fluoroquinolone resistance gene (gyrB) and three genes (otrA, tetO and tetBP) by RAST analysis possibly involved in tetracycline resistance were discovered. Moreover, 18 multidrug resistance efflux pumps (Table 4) were revealed, which may confer non-specific resistance (Table 3).

Tab. 4. Comparison of similarity (% nucleotide sequence) of antimicrobial resistance genes in various Elizabethkingia*.
Comparison of similarity (% nucleotide sequence) of antimicrobial resistance genes in various <i>Elizabethkingia</i><em class="ref">*</em>.

Virulence factors predicted in E. meningoseptica in comparison to other Elizabethkingia spp.

The pathogenesis mechanisms in Elizabethkingia species remain largely unknown. When PathogenFinder was used to predict the probability of acting as a pathogen for strains Em1, Em2 and Em3, E. meningoseptica showed a “negative” result according to the calculated probability values (~0.162) [31]. Instead, Elizabethkingia genomes only matched 5 non-pathogenic families, which was possibly due to the lack of similar flavobacterial genomes in the PathogenFinder database. However, 766 putative virulence factors in Em2 were successfully predicted by VFDB (cutoff value set to E-10). Forty-four of them in the bacterial virulence database also showed high identity to the selected Elizabethkingia, although there were some minor variations (S1 Table). Many virulence genes were predicted to be involved in CME-1 formation, capsule polysaccharide synthesis, lipooligosaccharide (LOS) synthesis, “attacking” enzymes (proteases), superoxide dismutase, catalases, peroxidase, heat shock protein, two-component regulatory system and many others (see S1 Table).

E. meningoseptica EM2 carried at least 6 genes involved in sialic acid metabolism predicted by RAST while Em1 and Em2 only possessed 4 of them (Table 5). Three genes encoded glucosamine-6-phosphate deaminase (nagB), glucosamine-fructose-6-phosphate aminotransferase (glmS) and phosphoglucosamine mutase (glmB), which were involved in the de novo synthesis pathway of sialic acid. Gene products for sialic acid synthesis were highly conserved (identity > 93%) among the three selected Elizabethkingia spp (Table 5). Two copies of sialic acid transporter genes (neuC1 and neuC2, encoding UDP-N-acetylglucosamine 2-epimerase) were found in strains Em2, G4120 and 61421 PRCM; CSID_3000515919 only carried only one copy (neuC2). An nanH gene, encoding a candidate sialidase, was predicted as the virulence factor (Table 5). Among the selected Elizabethkingia, only E. meningoseptica carries nanH (Table 5). Moreover, neither E. miricola (except strain CSID_3000517120) nor E. anophelis had genes encoding sialidase A (nanH) or sialic acid transporter genes (neu C1 and C2).

Tab. 5. Comparison of selected sialic acid synthesis, transport and utilization genes among Elizabethkingia*.
Comparison of selected sialic acid synthesis, transport and utilization genes among <i>Elizabethkingia</i><em class="ref">*</em>.

The selected E. meningoseptica species produced biofilm on plastic surfaces (S5 Fig). Compared to those in strain Em2, gene clusters for capsular polysaccharide synthesis consisting of cap8E, cap8G, cap8O, cap8D, cap4F, cap8F, cps4D, Cj1137c, capD and cap4D were only conserved in strains G4120, 61421 PRCM as well as E. miricola CSID 3000517120 (S2 Table). Further, the gene products responsive for curli biosynthesis and assembly were highly conserved in E. meningoseptica species while they were absent in other Elizabethkingia species (S2 Table). On the other hand, the flagellin biosynthesis protein FlgD and flagellar motor protein MotB were found in all of Elizabethkingia species (S2 Table). However, Elizabethkingia may not produce functional flagella because they lack of most of the flagellin structure proteins, which is congruent with the fact that bacteria Elizabethkingia are non-motile.

Polysaccharide utilization loci and carbohydrate active enzymes

E. meningoseptica carried 107 (Em2) or 109 (Em1, Em3) CAZyme-encoding genes occupying ca. 3% of the bacterial genome (S3 Table). The predicted CAZyme repertoires in the three E. meningoseptica genomes were similar, with 58 glycoside hydrolases (GHs) and two polysaccharide lyases (PLs) distributed across the same families. The main distinction between Em2 and Em1/Em3 was absence of GT2 and GT32 sequences (S3 Table). The CAZyme repertoires of Em1, Em2 and Em3 strains were similar to other E. meningoseptica strains and showed the same features which distinguished E. meningoseptica from other Elizabethkingia species (S3 Table). Notably, E. meningoseptica strains did not encode single copies of GH1 (β-glycosidase), GH5 (subfamily 46) and CBM6 (b-glucan binding) genes systematically found in other species, and had a reduced number of genes encoding enzymes from families GH95 and GH130 (S3 Table). However, E. meningoseptica strains specifically encoded a GH33 (sialidase) and an extra GH13 protein. Analysis revealed that E. meningoseptica utilized a battery of carbon sources, including D-maltose, D-trehalose, D-gentibiose, D- melibiose, D-glucose, D-mammose, D-fructose, D-fucose, D-mannitol and D-glycerol[11]. Strains Em1, Em2 and Em3 contained 25 loci encoding SusC-SusD homologs in their respective genomes [32], with only four PULs containing GHs. These numbers were similar to those found in other E. meningoseptica (23–24 SusCD loci; 5 PULs with GHs) and E. anophelis (25 SusCD loci and 4–5 PULs with GHs per genome; n = 17) strains. The three clinical strains did not exhibit the PUL that contains the gene encoding the GH13_5 enzyme (α-amylase; conserved but translocated elsewhere in the genome), that was conserved in all Elizabethkingia genomes analyzed, including other E. meningoseptica strains. However, all six E. meningoseptica strains exhibited a specific PUL encoding a GH30_3 enzyme (β-1,6-glucanase) not present in other Elizabethkingia genomes (S4 Table). They also harbor an additional GH13 gene in the Elizabethkingia-conserved PUL with GH63-GH97 which likely targets storage polysaccharides like glycogen and starch, as in the prototypic starch utilization system [33], and an additional GH16 gene in the Elizabethkingia-conserved PUL that encodes a GH29 appended to a CBM32 (S4 Table). Interestingly, the only GH-containing PUL conserved in all Elizabethkingia genomes was also widely distributed across the Bacteroidetes phylum, and encoded enzymes from families GH3 (various b-glycosidases) and GH144 (β-1,2-glucanase) which suggests action on bacterial β-1,2-glucans.

Regulatory systems in Elizabethkingia spp.

The Em1 genome encoded 56 two-component system proteins, 191 transcription factor proteins, and 7 other DNA-binding proteins (Table 6). Em3 had the same regulatory systems as those in Em1 (Table 6). The Em2 genome encoded 56 predicted two-component system proteins, 187 transcription factor proteins, and 6 other DNA-binding proteins. It seemed that strain Em1 or Em3 had more regulatory capacity than strain Em2 primarily due to the number of transcriptional regulators (TRs), one-component systems (OCS), and DNA-binding proteins (ODP) proteins. Most E. meningoseptica strains possessed similar numbers of total regulatory protein (ranging from 246 to 254) and components (Table 6). Comparable regulatory systems occurred in E. anophelis, although there were more ODPs in E. anophelis CSID_3015183678 and Ag1. Remarkably, the total number of TRs (>140) in E. miricola was much higher than those in E. meningoseptica and E. anophelis (Table 6). Particularly, the TRs mostly accounted for the higher regulatory systems due to increasing AraC family transcription factors (Table 6). AraC family regulators control a variety of cellular processes including carbon metabolism, stress responses and virulence. Moreover, the two component systems and DNA-binding proteins were more abundant in the three selected E. miricola strains than those in E. meningoseptica.

Tab. 6. Predicted regulatory systems in Elizabethkingia*.
Predicted regulatory systems in <i>Elizabethkingia</i><em class="ref">*</em>.

CRISPR-Cas systems and prophages

The clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated protein (CRISPR/Cas) system was predicted in the Em1 and Em3 genome; no such sequences were found in the Em2 chromosome and the other 5 E. meningoseptica strains (S5 Table). The analyses unveiled that Em1 had a directed repeat (47-bp, GTTGTGCAGTATCACAAATATACTGTAAAATGAAAGCAGTTCACAAC) and had 40 spacers (S5 Table). CRISPR sequences were completely conserved in Em3 (coverage 100%, identity 100%). By scanning the other selected E. meningoseptica genomes in the NCBI database, E. meningoseptica strain 4076 (CP016376) and NBRC 12535 had a predicted CRISPR (DR length: 47, number of spacers: 21). Em1and Em3 genomes carried only one cas gene (type II CRISPR RNA-guided endonuclease Cas9) flanking the CRISPR regions. The amino acid sequences of Cas9 in Em1 showed the highest identity (68%) to that in Capnocytophaga spp. Only incomplete prophage elements were found in Em1 (2), Em2 (1), and Em3 (2). Moreover, most of these E. meningoseptica strains did not have the complete prophages except strain CSID_3000515919 (S6 Table).

Discussion

E. meningoseptica isolated from hospitalized patients in Michigan carried many virulence factors participating in biofilm formation, proteases, lipooligosaccharide (LOS) synthesis, iron uptake and transportation, heat shock proteins, and capsule formation, highlighting that they have great potential to invade and colonize animal hosts. Further, our results also showed that these E. meningoseptica strains had several uncommon features including multi-drug resistance, sialic acid synthesis and transportation, and curli formation. Comparative genome analysis also showed that E. meningoseptica differed from E. anophelis or E. miricola in many ways. For example, the average GC content in E. meningoseptica was higher than in E. anophelis or E. miricola (Table 1), highlighting that E. meningoseptica evolves differently from E. anophelis or E. miricola. The predicted functional genes indicated that representative E. meningoseptica (Em2) shared less than 86% of the total encoding genes with E. anophelis strains and E. miricola. In particular, the divergence of E. meningoseptica from the other Elizabethkingia species was revealed by the different prophages, virulence factors and CRISPR-Cas systems.

By studying how pan-genome size increases as a function of the number of genomes sampled, we can gain insight into a species’ genetic repertoire [25]. Regression analysis of the pan genomes of E. meningoseptica showed that the pan-genome of E. meningoseptica is evolving through the loss or gain of a range of genes since their divergence as those reported in many pathogens including E. anophelis, Flavobacterium psychrophilum and F. spartansii [16, 3436]. It is not surprising that E. meningoseptica has an open pan-genome because it lives in diverse ecological niches (both aquatic and terrestrial environments), and colonizes and multiplies in animal and plant hosts [1, 37, 38].

E. meningoseptica harbored an extensive array of specialized CAZymes (up to 117) for the metabolism of glucans with up to 58 GHs, congruent with metabolic capacity to utilize various carbon sources. Some of the GH genes are organized like the typical PULs that are widely found in Bacteroidetes [3941]. Differing from aquatic bacteria, terrestrial flavobacteria always carry diverse CAZymes with a large potential for the breakdown of polysaccharides from plants or animals [4143]. Instead, less CAZyme genes are expected in the genomes of marine flavobacteria because they usually utilize peptides (rhodopsins) and/or harvesting light under the nutrient stressing environment [41, 44, 45]. The CAZyme components and assembly patterns are important properties that reveal Elizabethkingia may utilize a variety of carbon sources [4, 4649]. These results further support the notion that Elizabethkingia can adapt to environmental variations between the terrestrial and aquatic native habitats they colonize [16, 18, 46]. Besides soil and water, Elizabethkingia are often associated with the amphibious animals such as frogs and insects (i.e. mosquitoes) living in both aquatic and terrestrial environments [49, 50].

The detailed mechanisms for Elizabethkingia transmission from the environment to healthcare facilities and to patients remain unclear [18]. However, infection by Elizabethkingia can be mediated by multiple pathways of exposure [18]. Regardless of the transmission pathways, the genes in E. meningoseptica strains involved in biofilm on the abiotic material are one of the most important virulence factors [1]. The biofilm allows bacteria to persist on the surface of the medical utilities and resist against disinfection reagents [1]. Moreover, attachment/adhesion to the external surfaces and/or tissues of animal hosts is critical in the course of E. meningoseptica infection [1, 36, 51]. Many virulence factor genes including liposaccharide, hemagglutinin, capsule, and curli formation are conserved in E. meningoseptica (S1 Table). Capsule components are well known to play a vital role in the adhesion and/or biofilm formation [16]. A previous study has shown that expression of hemagglutinin adhesins allowed bacterial attachment and subsequent cell accumulation on target substrates [52]. Curli fibers not only participate in biofilm formation but also contribute to bacterial adhesion to animal cells [52]. E. meningoseptica Em1, Em2 and Em3 had hydrophilic cell surfaces. However, it remains unclear if such features affect bacterial adhesion ability [1]. Jacobs and Chenia (2011) investigated biofilm formation and adherence characteristics of E. meningoseptica isolated from freshwater tilapia [1]. The results showed that E. meningoseptica displayed better biofilm formation in nutrient-rich medium than that in nutrient-limited one at both low and high temperature [1].

Sialic acid is generally found at the terminal position(s) within glycan molecules covering animal cell surfaces [53]. Sialic acid is involved in various cellular processes including intercellular adhesion, cell signaling and immune system invasion [54]. In bacteria, sialic acid molecules participate in many pathogenesis processes [5456]. For instance, it can be integrated into cell components (i.e. the membrane lipopolysaccharide and capsule) that mimic the surface molecules of the host cell (called molecular mimicry), and thereby assist pathogens to escape the host innate immune response [56]. Moreover, sialic acid is also a good carbon as well as nitrogen source when environmental nutrients are limited [57]. Sialic acid enters bacteria through transporters and is next converted into fructose-6-phosphate; thus, it enters the central metabolism pathway [57, 58]. Additionally, sialic acid-rich glycolipids, glycoproteins and proteoglycans maintain water molecules at the bacterial surface, contributing to the uptake of polar molecules [56]. The detailed physiological roles of sialic acid remain unknown in flavobacteria. However, some E. meningoseptica strains (G4120, PRCM and CSID_3000515919) potentially acquire sialic acid by either uptake or the de novo synthesis pathway (Table 5). Other Elizabethkingia may only utilize sialic acid by de novo synthesis because they lack the exosialidase and sialic acid uptake machinery (Table 5). G. vaginalis carrying the putative sialidase A gene (nanH) was associated with the presence of vaginal biofilms, indicating that sialidase A is an important virulence factor for bacterial vaginosis (BV) [59]. nanH mutant failed to form biofilms in T. forsythia [60]. However, it restored the biofilm when purified NanH protein was added, showing that sialidase A was critical for pathogens to acquire nutrient source [60]. The same study further suggested that sialidase inhibitors might be useful adjuncts in periodontal therapy. In this study, genes encoding sialidase A (nanH) or sialic acid transporter genes (neuC1 and C2) were absent from E. miricola (except strain CSID_3000517120) and E. anophelis, indicating that different Elizabethkingia may display different pathogenesis mechanisms.

Antimicrobial susceptibility pattern for Elizabethkingia is often controversial [3, 18, 61]. Therefore, the selection of appropriate empiric therapy early in the course of bacterial infection can be very challenging. For example, Han et al [13] reported that 100% of E. meningoseptica isolates from South Korea were susceptible to piperacillin-tazobactam and less susceptible to ciprofloxacin (23% of isolates). However, E. meningoseptica clinical isolates from Taiwan were resistant to piperacillin/tazobactam while 74.4% of the strains were susceptible to trimethoprim/sulfamethoxaz [7]. The antibiotic susceptibility patterns in Michigan isolates mimic those isolates from Taiwan. Such discrepancies indicate that E. meningoseptica from different geographical regions may evolve different antibiotic resistance mechanisms [62]. Thus, the investigation of the molecular mechanisms involving in antibiotic-resistance is required, which could directly contribute to the empirical treatment of E. meningoseptica infections [62, 63]. Remarkably, the genomes of E. meningoseptica carried intrinsic class A extended-spectrum β-lactamases (ESBLs) and inherent class B metallo-β-lactamases (MBLs) [15]. These genes are likely to confer resistance to these β-lactam antibiotics including the combination drug piperacillin/tazobactam (Table 3). E. meningoseptica genomes shared many antibiotic-resistance genes with those from other Elizabethkingia species while there was minor differences [14, 64]. Besides those antibiotic-inactivating enzymes/proteins, Elizabethkingia may utilize the multidrug efflux pumps to excrete a wide range of antibiotics [65], including abes S, msrB, macB, ceoB, adeG and many others (Table 4). For example, CeoB, a component of RND multidrug efflux system, pumps out chloramphenicol and ciprofloxacin [65, 66]. MsrB, as an ABC-efflux pump, was reported in Staphylococcus species to confer resistance to erythromycin and streptogramin B antibiotics [67]. Genes involving in tetracycline and vancomycin resistance were detected in Elizabethkingia in our study and others [64], which may explain why these drugs are not very effective in eliminating Elizabethkingia infections. Analysis of the antibiotic resistance in E. meningoseptica showed that they carried genes conferring to drug resistance including aminoglycosides, isoniazid, streptogramin, aminosalicylate, fluoroquinolone, macrolide, tetracycline, diaminopyrimidine, phenico, glycopeptide and rifamycin. Instead, we did not detect genes involving in resistance of sulfonamide (folate pathway inhibitors) in the genomes of Michigan isolates, which agrees with their susceptibility to trimethoprim/sulfamethoxazole. Mutation of the gyrB gene encoding the DNA gyrase subunit B in strains Em1, Em2 and Em3 was predicted to contribute to fluoroquinolones’ resistance (such as ciprofloxacin) in this study.

Prophages are known to modulate virulence and antibiotic resistance gene expression, which alters the production and/or secretion of toxins during the infection course [68, 69]. Further, they are one of the important contributors to genetic diversification by transferring the functional genes among the different strains [70]. The rare occurrence of complete prophages in E. meningoseptica remains unclear. Our analysis revealed only few of the selected genomes (Em1 and Em3) had “confirmed” CRISPRs (Supplementary S5 Table). The similar observations were reported in E. anophelis genomes [16]. The defense of the invasions of foreign genetic elements such as plasmids, transposons or phages may require both restriction modification systems (RMs) and Clustered Regularly Interspaced Short Palindromic repeat sequences (CRISPRs) in Elizabethkingia [71]. However, the detailed mechanisms need to be further investigated.

Supporting information

S1 Fig [docx]
Comparisions of the GC contents in . , . and . .

S2 Fig [docx]
Synteny plots show the conservation of gene order among the selected genomes.

S3 Fig [docx]
Phylogenetic relationship among the selected .

S4 Fig [docx]
The shared genes among the selected . Em1, Em2 and Em3.

S5 Fig [docx]
biofilm assay in the selected . .

S1 Table [docx]
Comparations of the predicted virulence factors among selected spp.

S2 Table [docx]
Comparative analysis of the putative biofilm formation in spp.

S3 Table [xlsx]
Carbohydrate-Active Enzymes (CAZyme) analysis in selected spp.

S4 Table [xlsx]
Polysaccharide Utilization Loci (PULs) in the selected spp.

S5 Table [xlsx]
Predicted CRISPR-Cas systems in spp.

S6 Table [xlsx]
PHAST in spp.


Zdroje

1. Jacobs A, Chenia HY. Biofilm formation and adherence characteristics of an Elizabethkingia meningoseptica isolate from Oreochromis mossambicus. Ann Clin Microbiol Antimicrob. 2011;10:16–. doi: 10.1186/1476-0711-10-16 PMC3112384. 21545730

2. Luke SPM, Daniel SO, Annette J, Jane FT, Simon A, Hugo D, et al. Waterborne Elizabethkingia meningoseptica in adult critical care. Emerg Infect Diseases. 2016;22(1):9. doi: 10.3201/eid2201.150139 26690562

3. Hsu M-S, Liao C-H, Huang Y-T, Liu C-Y, Yang C-J, Kao K-L, et al. Clinical features, antimicrobial susceptibilities, and outcomes of Elizabethkingia meningoseptica (Chryseobacterium meningosepticum) bacteremia at a medical center in Taiwan, 1999–2006. Eur J Clin Microbiol Infect Dis 2011;30(10):1271–8. doi: 10.1007/s10096-011-1223-0 21461847

4. Jean SS, Lee WS, Chen FL, Ou TY, Hsueh PR. Elizabethkingia meningoseptica: an important emerging pathogen causing healthcare-associated infections. J Hosp Infect. 2014;86(4):244–9. doi: 10.1016/j.jhin.2014.01.009 24680187

5. Balm MND, Salmon S, Jureen R, Teo C, Mahdi R, Seetoh T, et al. Bad design, bad practices, bad bugs: frustrations in controlling an outbreak of Elizabethkingia meningoseptica in intensive care units. J Hosp Infect 2013;85(2):134–40. doi: 10.1016/j.jhin.2013.05.012 23958153

6. Chawla K, Gopinathan A, Varma M, Mukhopadhyay C. Elizabethkingia meningoseptica outbreak in intensive care unit. J Global Infect Dis. 2015;7(1):43–4. doi: 10.4103/0974-777X.150890 25722622.

7. Chang Y-C, Lo H-H, Hsieh H-Y, Chang S-M. Identification and epidemiological relatedness of clinical Elizabethkingia meningoseptica isolates from central Taiwan. J Microbiol Immunol Infect 2014;47(4):318–23. doi: 10.1016/j.jmii.2013.03.007 23726463

8. de Carvalho Filho ÉB, Marson FAL, Levy CE. Challenges in the identification of Chryseobacterium indologenes and Elizabethkingia meningoseptica in cases of nosocomial infections and patients with cystic fibrosis. New Microbes New Infect. 2017;20:27–33. doi: 10.1016/j.nmni.2017.09.002 29062487.

9. Cenk MH, Özlem T, Serpil Ö, Burçin Ş, Dilşa TG, Uğur Ö, et al. Clinical strains of Chryseobacterium and Elizabethkingia spp. isolated from pediatric patients in a university hospital: performance of MALDI-TOF MS-based identification, antimicrobial susceptibilities, and baseline patient characteristics. Microbial Drug Resistance. 2018;24(6):816–21. doi: 10.1089/mdr.2017.0206 29227188.

10. Lau SKP, Chow W-N, Foo C-H, Curreem SOT, Lo GC-S, Teng JLL, et al. Elizabethkingia anophelis bacteremia is associated with clinically significant infections and high mortality. Sci Rep. 2016;6:26045. doi: 10.1038/srep26045 PMC4868968. 27185741

11. Kim KK, Kim MK, Lim JH, Park HY, Lee S-T. Transfer of Chryseobacterium meningosepticum and Chryseobacterium miricola to Elizabethkingia gen. nov. as Elizabethkingia meningoseptica comb. nov. and Elizabethkingia miricola comb. nov. Int J Syst Evol Microbiol. 2005;55(3):1287–93. doi: 10.1099/ijs.0.63541–0

12. Chew KL, Cheng B, Lin RTP, Teo JWP. Elizabethkingia anophelis is the dominant Elizabethkingia species found in blood cultures in Singapore. J Clin Microbiol. 2018;56(3):e01445–17. doi: 10.1128/JCM.01445-17 29237782.

13. Han M-S, Kim H, Lee Y, Kim M, Ku NS, Choi JY, et al. Relative prevalence and antimicrobial susceptibility of clinical isolates of Elizabethkingia species based on 16S rRNA gene sequencing. J Clin Microbiol. 2017;55(1):274–80. doi: 10.1128/JCM.01637-16 27847376

14. Opota O, Diene SM, Bertelli C, Prod'hom G, Eckert P, Greub G. Genome of the carbapenemase-producing clinical isolate Elizabethkingia miricola EM_CHUV and comparative genomics with Elizabethkingia meningoseptica and Elizabethkingia anophelis: evidence for intrinsic multidrug resistance trait of emerging pathogens. Int J Antimicrob Agents. 2017;49(1):93–7. doi: 10.1016/j.ijantimicag.2016.09.031 27913093

15. González LJ, Vila AJ. Carbapenem resistance in Elizabethkingia meningoseptica is mediated by metallo-β-lactamase BlaB. Antimicrob Agents and Chemother. 2012;56(4):1686–92. doi: 10.1128/aac.05835-11 22290979

16. Perrin A, Larsonneur E, Nicholson AC, Edwards DJ, Gundlach KM, Whitney AM, et al. Evolutionary dynamics and genomic features of the Elizabethkingia anophelis 2015 to 2016 Wisconsin outbreak strain. Nat Commun. 2017;8:15483. doi: 10.1038/ncomms15483 https://www.nature.com/articles/ncomms15483-supplementary-information. 28537263

17. Lin J-N, Lai C-H, Yang C-H, Huang Y-H, Lin H-H. Genomic features, phylogenetic relationships, and comparative genomics of Elizabethkingia anophelis strain EM361-97 isolated in Taiwan. Sci Rep. 2017;7(1):14317. doi: 10.1038/s41598-017-14841-8 29085032

18. Michael Janda J., Lopez DL. Mini review: New pathogen profiles: Elizabethkingia anophelis. Diagn Microbiol Infect Dis 2017;88(2):201–5. doi: 10.1016/j.diagmicrobio.2017.03.007 28342565

19. CLSI. Performance Standards for Antimicrobial Susceptibility Testing. 26th ed. CLSI supplement M100S. Wayne PCaLSI.

20. Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19(5):455–77. doi: 10.1089/cmb.2012.0021 PMC3342519. 22506599

21. Overbeek R, Olson R, Pusch GD, Olsen GJ, Davis JJ, Disz T, et al. The SEED and the rapid annotation of microbial genomes using subsystems technology (RAST). Nucleic Acids Res. 2014;42(Database issue):D206–D14. doi: 10.1093/nar/gkt1226 PMC3965101. 24293654

22. Montaner B, Navarro S, Piqué M, Vilaseca M, Martinell M, Giralt E, et al. Prodigiosin from the supernatant of Serratia marcescens induces apoptosis in haematopoietic cancer cell lines. Br J Pharmacol. 2000;131(3):585–93. doi: 10.1038/sj.bjp.0703614 PMC1572367. 11015311

23. Zhou Y, Liang Y, Lynch KH, Dennis JJ, Wishart DS. PHAST: A Fast Phage Search Tool. Nucleic Acids Research. 2011. doi: 10.1093/nar/gkr485 21672955

24. Grissa I, Vergnaud G, Pourcel C. CRISPRFinder: a web tool to identify clustered regularly interspaced short palindromic repeats. Nucleic Acids Res. 2007;35(suppl 2):W52–W7. doi: 10.1093/nar/gkm360 17537822

25. Blom J, Kreis J, Spänig S, Juhre T, Bertelli C, Ernst C, et al. EDGAR 2.0: an enhanced software platform for comparative gene content analyses. Nucleic Acids Res. 2016:W22–W8. doi: 10.1093/nar/gkw255 27098043

26. Tettelin H, Masignani V, Cieslewicz MJ, Donati C, Medini D, Ward NL, et al. Genome analysis of multiple pathogenic isolates of Streptococcus agalactiae: Implications for the microbial “pan-genome”. PNAS USA. 2005;102(39):13950–5. doi: 10.1073/pnas.0506758102 16172379

27. Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Research. 2004;32(5):1792–7. doi: 10.1093/nar/gkh340 15034147

28. Felsenstein J. PHYLIP—phylogeny inference package (Version 3.2). Cladistics. 1989;5:164–6.

29. Lombard V, Golaconda Ramulu H, Drula E, Coutinho PM, Henrissat B. The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res. 2014;42(D1):D490–D5. doi: 10.1093/nar/gkt1178 24270786

30. Chen S, Blom J, Walker ED. Genomic, physiologic, and symbiotic characterization of Serratia marcescens strains isolated from the mosquito Anopheles stephensi. Front Microbiol. 2017;8(1483). doi: 10.3389/fmicb.2017.01483 28861046

31. Cosentino S, Voldby Larsen M, Møller Aarestrup F, Lund O. PathogenFinder—distinguishing friend from foe using bacterial whole genome sequence data. PLOS ONE. 2013;8(10):e77302. doi: 10.1371/journal.pone.0077302 24204795

32. Terrapon N, Lombard V, Drula É, Lapébie P, Al-Masaudi S, Gilbert HJ, et al. PULDB: the expanded database of polysaccharide utilization loci. Nucleic Acids Res. 2018;46(D1):D677–D83. doi: 10.1093/nar/gkx1022 29088389

33. Foley MH, Cockburn DW, Koropatkin NM. The Sus operon: a model system for starch uptake by the human gut Bacteroidetes. Cellular and molecular life sciences: CMLS. 2016;73(14):2603–17. doi: 10.1007/s00018-016-2242-x 27137179.

34. Duchaud E, Rochat T, Habib C, Barbier P, Loux V, Guérin C, et al. Genomic Diversity and Evolution of the Fish Pathogen Flavobacterium psychrophilum. Frontiers in Microbiology. 2018;9(138). doi: 10.3389/fmicb.2018.00138 29467746

35. Chen S, Blom J, Loch TP, Faisal M, Walker ED. The emerging fish pathogen Flavobacterium spartansii isolated from chinook salmon: comparative genome analysis and molecular manipulation. Front Microbiol. 2017;8(2339). doi: 10.3389/fmicb.2017.02339 29250046

36. Teo J, Tan SY-Y, Liu Y, Tay M, Ding Y, Li Y, et al. Comparative genomic analysis of malaria mosquito vector associated novel pathogen Elizabethkingia anophelis. Genome Biol Evol 2014: doi: 10.1093/gbe/evu094 24803570

37. Chen S, Zhao J, Joshi D, Xi Z, Norman B, Walker ED. Persistent infection by Wolbachia wAlbB has no effect on composition of the gut microbiota in adult female Anopheles stephensi. Front Microbiol. 2016;7(1485). doi: 10.3389/fmicb.2016.01485 27708633

38. Akhouayri IG, Habtewold T, Christophides GK. Melanotic pathology and vertical transmission of the gut commensal Elizabethkingia meningoseptica in the major malaria vector Anopheles gambiae. PLoS ONE. 2013;8(10):e77619. doi: 10.1371/journal.pone.0077619 24098592

39. Dodd D, Moon YH, Swaminathan K, Mackie RI, Cann IK. Transcriptomic analyses of xylan degradation by Prevotella bryantii and insights into energy acquisition by xylanolytic Bacteroidetes. J Biol Chem. 2010;285(39):30261–73. Epub 2010/07/14. doi: 10.1074/jbc.M110.141788 20622018; PubMed Central PMCID: PMC2943253.

40. Flint HJ, Scott KP, Duncan SH, Louis P, Forano E. Microbial degradation of complex carbohydrates in the gut. Gut Microbes. 2012;3(4):289–306. doi: 10.4161/gmic.19897 PMC3463488. 22572875

41. Kolton M, Sela N, Elad Y, Cytryn E. Comparative genomic analysis indicates that niche adaptation of terrestrial Flavobacteria is strongly linked to plant glycan metabolism. PLoS ONE. 2013;8(9):e76704. doi: 10.1371/journal.pone.0076704 PMC3784431. 24086761

42. Chen S, Kaufman MG, Miazgowicz KL, Bagdasarian M, Walker ED. Molecular characterization of a cold-active recombinant xylanase from Flavobacterium johnsoniae and its applicability in xylan hydrolysis. Bioresour Technol. 2013;128(0):145–55. http://dx.doi.org/10.1016/j.biortech.2012.10.087.

43. Nagal S, Jain PC. Production of feather hydrolysate by Elizabethkingia meningoseptica KB042 (MTCC 8360) in submerged fermentation. ndian J Microbiol. 2010;50(1):41–5. doi: 10.1007/s12088-010-0014-0 22815570

44. Tekedar HC, Karsi A, Reddy JS, Nho SW, Kalindamar S, Lawrence ML. Comparative genomics and transcriptional analysis of Flavobacterium columnare strain ATCC 49512. Front Microbiol. 2017;8(588). doi: 10.3389/fmicb.2017.00588 28469601

45. Yoshizawa S, Kumagai Y, Kim H, Ogura Y, Hayashi T, Iwasaki W, et al. Functional characterization of flavobacteria rhodopsins reveals a unique class of light-driven chloride pump in bacteria. Proc Natl Acad Sci U S A 2014;111(18):6732–7. doi: 10.1073/pnas.1403051111 24706784

46. Kukutla P, Lindberg BG, Pei D, Rayl M, Yu W, Steritz M, et al. Insights from the genome annotation of Elizabethkingia anophelis from the malaria vector Anopheles gambiae. PLoS ONE. 2014;9(5):10.1371/journal.pone.0097715. doi: 10.1371/journal.pone.0097715 24842809

47. Kämpfer P, Busse H-J, McInroy JA, Glaeser SP. Elizabethkingia endophytica sp. nov., isolated from Zea mays and emended description of Elizabethkingia anophelis Kämpfer et al. 2011. Int J Syst Evol Microbiol. 2015;65(7):2187–93. doi: 10.1099/ijs.0.000236 25858248

48. Kämpfer P, Matthews H, Glaeser SP, Martin K, Lodders N, Faye I. Elizabethkingia anophelis sp. nov., isolated from the midgut of the mosquito Anopheles gambiae. Int J Syst Evol Microbiol. 2011;61(11):2670–5. doi: 10.1099/ijs.0.026393–0

49. Chen S, Bagdasarian M, Walker ED. Elizabethkingia anophelis: molecular manipulation and interactions with mosquito hosts. Appl Environ Microbiol. 2015. doi: 10.1128/aem.03733-14 25595771

50. Ruixue H, Junfa Y, Yin M, Zhe W, Zemao G. Pathogenic Elizabethkingia miricola infection in cultured black-spotted frogs, China, 2016. Emerg Infect Diseases. 2017;23(12):2055. doi: 10.3201/eid2312.170942 29148374

51. Chen S, Soehnlen M, Downes FP, Walker ED. Insights from the draft genome into the pathogenicity of a clinical isolate of Elizabethkingia meningoseptica Em3. Stand Genomic Sci. 2017;12:56. doi: 10.1186/s40793-017-0269-8 PMC5602931. 28932346

52. Darvish Alipour Astaneh S, Rasooli I, Mousavi Gargari SL. The role of filamentous hemagglutinin adhesin in adherence and biofilm formation in Acinetobacter baumannii ATCC19606(T). Microb Pathog. 2014;74:42–9. doi: 10.1016/j.micpath.2014.07.007 25086432.

53. Cioffi DL, Pandey S, Alvarez DF, Cioffi EA. Terminal sialic acids are an important determinant of pulmonary endothelial barrier integrity. Am J Physiol Lung Cell Mol Physiol. 2012;302(10):L1067–L77. doi: 10.1152/ajplung.00190.2011 PMC3362258. 22387293

54. Reutter W, Stäsche R, Stehling P, Baum O. The biology of sialic acids: insights into their structure, metabolism and function in particular during viral infection. Glycosciences.

55. Opal SM. Significance of sialic acid in Klebsiella pneumoniae K1 capsules. Virulence. 2014;5(6):648–9. doi: 10.4161/viru.34349 PMC4139404. 25105481

56. Anderson GG, Goller CC, Justice S, Hultgren SJ, Seed PC. Polysaccharide capsule and sialic acid-mediated regulation promote biofilm-like intracellular bacterial communities during cystitis. Infect Immun. 2010;78(3):963–75. doi: 10.1128/IAI.00925-09 20086090

57. McDonald ND, Lubin J-B, Chowdhury N, Boyd EF. Host-derived sialic acids are an important nutrient source required for optimal bacterial fitness in vivo. mBio. 2016;7(2). doi: 10.1128/mBio.02237-15 27073099

58. Severi E, Hood DW, Thomas GH. Sialic acid utilization by bacterial pathogens. Microbiol. 2007;153(9):2817–22. doi: 10.1099/mic.0.2007/009480-0 17768226

59. Hardy L, Jespers V, Van den Bulck M, Buyze J, Mwambarangwe L, Musengamana V, et al. The presence of the putative Gardnerella vaginalis sialidase A gene in vaginal specimens is associated with bacterial vaginosis biofilm. PLOS ONE. 2017;12(2):e0172522. doi: 10.1371/journal.pone.0172522 28241058

60. Honma K, Mishima E, Sharma A. Role of Tannerella forsythia NanH sialidase in epithelial cell attachment. Infect Immun. 2011;79(1):393–401. doi: 10.1128/IAI.00629-10 21078857

61. Leclercq R, Cantón R, Brown DFJ, Giske CG, Heisig P, MacGowan AP, et al. EUCAST expert rules in antimicrobial susceptibility testing. Clin Microbiol Infect. 2013;19(2):141–60. doi: 10.1111/j.1469-0691.2011.03703.x 22117544

62. Munita JM, Arias CA. Mechanisms of antibiotic resistance. Microbiol Spectr. 2016;4(2):10.1128/microbiolspec.VMBF-0016-2015. doi: 10.1128/microbiolspec.VMBF-0016-2015 PMC4888801. 27227291

63. Thabit AK, Crandon JL, Nicolau DP. Antimicrobial resistance: impact on clinical and economic outcomes and the need for new antimicrobials. Expert Opin Pharmacother. 2015;16(2):159–77. doi: 10.1517/14656566.2015.993381 25496207

64. Lin J-N, Lai C-H, Yang C-H, Huang Y-H, Lin H-H. Genomic features, phylogenetic relationships, and comparative genomics of Elizabethkingia anophelis strain EM361-97 isolated in Taiwan. Scientific Reports. 2017;7(1):14317. doi: 10.1038/s41598-017-14841-8 29085032

65. Li X-Z, Plésiat P, Nikaido H. The challenge of efflux-mediated antibiotic resistance in Gram-negative bacteria. Clin Microbiol Rev 2015;28(2):337–418. doi: 10.1128/CMR.00117-14 25788514

66. Burns JL, Hedin LA, Lien DM. Chloramphenicol resistance in Pseudomonas cepacia because of decreased permeability. Antimicrob Agents Chemother. 1989;33(2):136–41. doi: 10.1128/aac.33.2.136 2719457

67. Leclercq R. Mechanisms of resistance to macrolides and lincosamides: nature of the resistance elements and their clinical implications. Clin Infect Dis. 2002;34(4):482–92. doi: 10.1086/324626 11797175

68. Casas V, Miyake J, Balsley H, Roark J, Telles S, Leeds S, et al. Widespread occurrence of phage-encoded exotoxin genes in terrestrial and aquatic environments in Southern California. FEMS Microbiol Lett. 2006;261(1):141–9. doi: 10.1111/j.1574-6968.2006.00345.x 16842371

69. Wagner PL, Waldor MK. Bacteriophage control of bacterial virulence. Infect Immun. 2002;70(8):3985–93. doi: 10.1128/IAI.70.8.3985-3993.2002 PMC128183. 12117903

70. Waldor MK, Friedman DI. Phage regulatory circuits and virulence gene expression. Curr Opin Microbiol. 2005;8(4):459–65. doi: 10.1016/j.mib.2005.06.001 15979389

71. Dupuis M-È, Villion M, Magadán AH, Moineau S. CRISPR-Cas and restriction–modification systems are compatible and increase phage resistance. Nat Commun. 2013;4:2087. doi: 10.1038/ncomms3087 23820428


Článek vyšel v časopise

PLOS One


2019 Číslo 10
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy Podcasty Doporučená témata Časopisy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#