Xi Guo1,2,3†
Min Wang1,2,3†Lu Wang3†Yao Wang3Tingting Chen3Pan Wu3Min Chen4
Bin Liu1,2,3*Lu Feng1,2,3*- 1Key Laboratory of Molecular Microbiology and Technology, Ministry of Education–Tianjin Economic-Technological Development Area, Tianjin, China
- 2Tianjin Key Laboratory of Microbial Functional Genomics, Tianjin Economic-Technological Development Area, Tianjin, China
- 3TEDA Institute of Biological Sciences and Biotechnology, Nankai University–Tianjin Economic-Technological Development Area, Tianjin, China
- 4Shanghai Municipal Center for Disease Control and Prevention, Shanghai, China
Serotyping based on surface polysaccharide antigens is important for the clinical detection and epidemiological surveillance of pathogens. Polysaccharide gene clusters (PSgcs) are typically responsible for the diversity of bacterial surface polysaccharides. Through whole-genome sequencing and analysis, eight putative PSgc types were identified in 23 Enterobacter aerogenes strains from several geographic areas, allowing us to present the first molecular serotyping system for E. aerogenes. A conventional antigenic scheme was also established and correlated well with the molecular serotyping system that was based on PSgc genetic variation, indicating that PSgc-based molecular typing and immunological serology provide equally valid results. Further, a multiplex Luminex-based array was developed, and a double-blind test was conducted with 97 clinical specimens from Shanghai, China, to validate our array. The results of these analyses indicated that strains containing PSgc4 and PSgc7 comprised the predominant groups. We then examined 86 publicly available E. aerogenes strain genomes and identified an additional seven novel PSgc types, with PSgc10 being the most abundant type. In total, our study identified 15 PSgc types in E. aerogenes, providing the basis for a molecular serotyping scheme. From these results, differing epidemic patterns were identified between strains that were predominant in different regions. Our study highlights the feasibility and reliability of a serotyping system based on PSgc diversity, and for the first time, presents a molecular serotyping system, as well as an antigenic scheme for E. aerogenes, providing the basis for molecular diagnostics and epidemiological surveillance of this important emerging pathogen.
Introduction
Enterobacter aerogenes is a Gram-negative bacterium that is ubiquitous in the human gastrointestinal tract and various other environments and is generally non-pathogenic to healthy humans (Chen et al., 2015). However, E. aerogenes has become an important opportunistic pathogen since the early 1990s and has frequently been isolated from respiratory, urinary, blood, and gastrointestinal tract infections (Langley et al., 2001). Moreover, E. aerogenes has been reported to readily cause septic shock in patients, thus leading to high mortality rates (Lavigne et al., 2012). E. aerogenes strains that are isolated from hospitalized patients generally exhibit high resistance to commonly used broad-spectrum antibiotics (De Gheldre et al., 1997; Chen et al., 2008; Lavigne et al., 2013). Consequently, this species has been considered an important emerging multidrug resistant (MDR) pathogen over the past two decades (Chevalier et al., 2008).
Surface polysaccharides, particularly the heat-stable somatic (O) and/or capsular (K) antigens, are major targets for both host immune systems and bacteriophages and therefore are some of the most variable constituents of the cell. The genes involved in surface polysaccharide synthesis are generally arranged in a gene cluster, and most of the variation among polysaccharides is due to genetic variation in the polysaccharide gene clusters (PSgcs), thereby providing the molecular basis for serotyping. These genes are commonly classified into three main classes: (i) nucleotide sugar precursor synthesis genes for sugars that are specific to the particular polysaccharide; (ii) sugar transferase genes that are associated with the O- or K-unit assemblies, and which are specific for the donor and acceptor sugars and generate a specific linkage between them; and (iii) genes (wzx/wzy or wzm/wzt) for O- or K-unit translocation and polymerization (Liu et al., 2014). In addition to the three gene classes above, a wzi-wza-wab-wzc gene set that is highly conserved among different serotypes exists in the locus of most Klebsiella capsular types and the group 1 capsule of E. coli (Whitfield, 2006; Pan et al., 2015). The protein complex Wza/b/c accepts the polymer made by Wzy and transports it to the cell surface, forming the capsule (K) antigen (Whitfield, 2006). The remainder, Wzi, is proposed to play a role in capsule assembly, but the mechanistic details of this process remain unclear (Rahn et al., 2003). The complete O-antigen is ligated to the lipid A/core by WaaL, which is encoded by the waaL gene in the core gene cluster, thus forming the LPS structure. The K-repeat polymer can also be added to the lipid A/core, and this is known as the KLPS, which does not require the Wza/b/c complex (Hu et al., 2013).
The recent development of molecular detection methods, including pulsed-field gel electrophoresis (PFGE), multilocus sequence typing (MLST), plasmid fingerprinting (PF), restriction fragment length polymorphism (RFLP) and gene-specific PCR, have advanced the ability to identify pathogens. Ribotyping, arbitrarily primed PCR and PFGE assays have been reported to be used for epidemiological study of nosocomial E. aerogenes isolates (Grattard et al., 1995; Jalaluddin et al., 1998). However, serotyping based on the diversity of surface polysaccharides that underlie antigenic schemes remains the ‘gold standard’ for detecting pathogenic strains in clinical specimens and environmental samples, as well as for epidemiological surveillance and tracing (Wang et al., 2010). In contrast to widely studied pathogens, including Salmonella spp. (Guibourdenche et al., 2010), E. coli (Orskov et al., 1977), S. pneumoniae (Konradsen and Pneumococcus Reference laboratories in Europe, 2005), and Vibrio parahaemolyticus (Chen et al., 2012) and other major pathogens, antigenic schemes for some important emerging pathogens have not been established. Furthermore, even the genetic and structural basis for serological diversity has not yet been elucidated in several of these emerging pathogens. The development of next-generation sequencing (NGS) has made it possible to rapidly detect almost all genetic features of numerous pathogens, including elucidation of putative PSgcs in novel, emerging pathogenic bacteria. Consequently, a molecular serotyping scheme could be developed based on the sequence diversity of PSgcs that is theoretically consistent with the traditional antigenic system based on the structural diversity of surface polysaccharide antigens.
Here, the genomes of 23 E. aerogenes strains were sequenced, and eight forms of putative PSgcs (1–8) were obtained and analyzed, providing confirmation of the molecular serotyping scheme for E. aerogenes. Antisera against type strains that represented each PSgc type were then produced and showed good absorption specificity to their homologous strains, indicating a high level of accuracy for the molecular serotyping scheme that we presented. Furthermore, a Luminex bead-based suspension array targeting these eight PSgc types was developed and evaluated for specificity and sensitivity in double-blind tests. We then examined another 86 E. aerogenes strains with genome sequences available in the GenBank database and identified an additional seven new gene cluster types (PSgcs 9–15).
Materials and Methods
Sequencing and Bioinformatic Analyses
Bacterial cultivation and DNA extraction were conducted in a BS-2 lab.
Whole-genome sequencing of 23 E. aerogenes strains (Supplementary Table 1) was performed on the Solexa paired-end sequencing platform. Genomic DNA was sheared, polished, and prepared using the Illumina Sample Preparation Kit according to the manufacturer’s protocols. Genomic libraries were constructed that contained 500-bp paired-end inserts, and sequencing was performed with Solexa sequencing technologies (Illumina, Inc.) to produce ∼100-fold coverage for each genome. Sequence reads were assembled using the de novo genome-assembly program Velvet to generate multi-contig draft genomes. Gaps within the putative PSgc were closed by directed PCR, and the products were sequenced using BigDye terminator chemistry on ABI 3730 capillary sequencers.
Artemis (Rutherford et al., 2000) was used to annotate genes, and the lockMaker program (Henikoff et al., 1995) was used to identify conserved motifs. BLAST and PSI-BLAST (Altschul et al., 1997) were used to search genes and proteins against available databases, including GenBank1 and the Pfam protein motif databases2. The TMHMM analysis program v2.03 was used to identify potential transmembrane domains within protein sequences.
Identification of Putative PSgcs From Genomes
Eighty-six E. aerogenes genome sequences were downloaded from the GenBank database (Supplementary Table 1) and combined with the 23 genome sequences obtained in this study to assess PSgc distributions. Putative PSgc sequences that were located between the two housekeeping genes galF and gnd were retrieved from the genomes. PSgc sequences sharing high DNA or protein level identity (>97%) and possessing the same gene organization were considered identical serotypes.
Preparation of Antigens for Immunization
All strains that were used for immunization were cultivated in 400 ml of media and were harvested by centrifugation at 5,000 × g for 20 min, washed in 20 ml of 0.85% NaCl, and suspended in 20 ml of 0.85% NaCl (approximately 1011 cells/ml). The cell suspension was subsequently treated with 0.5% formaldehyde overnight. Only strains that did not exhibit autoagglutination were used for immunization. The resultant antigens were stored at 4°C.
Production of Antisera
Adult New Zealand white rabbits (2.5 kg weight) were injected intravenously with the prepared antigens. Three rabbits were used for each Enterobacter strain. At 3-day intervals, injections were given at doses of 0.5, 1, 2, and 4 ml. Rabbits were exsanguinated 8 days after the final injection, and the separated antisera were stored at 4°C.
Tube Agglutination
Serotyping and titer determination were performed by tube agglutination. For antigen preparation, cultures were taken from agar plates after incubation at 37°C overnight, autoclaved, and then washed in 20 ml of 0.85% NaCl, and suspended in 20 ml of 0.85% NaCl. Tube agglutination was performed in 96-well U-bottom microtiter plates with 25 μl of antigen and 25 μl of twofold diluted antiserum in PBS (50 mM phosphate buffer, 150 mM NaCl, pH 7.2). Microtiter plates were incubated at 37°C for 4 h and then stored at 4°C overnight. The titer was determined based on the most diluted concentration of a serum that gave a positive reaction compared to the PBS negative control.
Antiserum Absorption
To assess the absorption of agglutinins, cell suspensions were prepared by inoculating moist, thickly poured (∼1 ml medium per plate) 90-mm infusion agar plates. Plates were incubated in an upright position for 18–24 h at 37°C. Cultures from the plates were then suspended in 0.85% NaCl and treated with 0.5% formaldehyde overnight and then centrifuged at 5,000 × g for 20 min. The cell pellets were then gently washed with 0.85% NaCl three times and resuspended in 2 ml of antisera. The mixtures were incubated at room temperature for 2 h and subsequently centrifuged at 10,000 × g for 15 min, and the supernatants were collected. The absorbed antisera were then tested against all antigens that reacted with the unabsorbed antisera. This process was repeated until cross-reactions no longer occurred.
Design of Primers and Probes
PSgc1-8 specific primers and probes were designed based on the DNA sequences that were obtained in this study using Primer Premier v5.0 software (Premier Biosoft International, Palo Alto, CA, United States). The specificity of each individual primer pair was confirmed using BLAST program and was subsequently validated by a single PCR amplification. For Luminex-based detection, each probe was synthesized with a 5′ amino modifier C12 that enabled coupling to the carboxyl group located on the microsphere, and each reverse primer was labeled with biotin at its 5′ end.
Multiplex PCR Amplification
Multiplex PCR amplification was performed in a 50 μl reaction mixture consisting of 100 ng of genomic DNA, 1× Goldstar PCR buffer, 20 μM of each deoxynucleoside triphosphate (dNTP), 2.5 units of Goldstar DNA polymerase, 0.5 μM of each forward primer, and 2 μM of each reverse primer. The reaction parameters were as follows: hot start at 95°C for 5 min; amplification stage of 30 cycles at 95°C for 30 s, 60°C for 30 s, and 72°C for 1 min, followed by a final extension at 72°C for 5 min.
Hybridization and Luminex Analysis
Each PSgc probe was bound to a different carboxylated microsphere as described previously (Chen et al., 2010). The working microspheres consisted of eight types of beads, and each bead was coupled to a PSgc-specific probe. Hybridization was performed in a 50 μl mixture that included 17 μl of biotin-labeled PCR product and 33 μl of working microspheres, with denaturation at 95°C for 10 min and incubation at 55°C for 17 min in a thermal cycler. The hybridization product was then transferred to a filter plate and washed three times with 1× TMAC buffer at 1,000 rpm for 1 min. For detection, 80 μl of streptavidin-R-phycoerythrin in 1× TMAC buffer was added to each well, followed by incubation at 53°C for 20 min. Finally, the signals for each set of beads were measured using a Bio-plex 100 reader (Bio-Rad) according to the manufacturer’s protocol. Data were analyzed using Bio-plex Manager 4.0, and the data are presented as the median fluorescence intensity (MFI). The cut-off value for a positive result was defined as three times greater than the mean MFI value of the background.
Ethics Statement
All animal experiments were performed in accordance with the standards established in the Guide for the Care and Use of Laboratory Animals published by the Institute of Laboratory Animal Resources of the National Research Council (United States). The animal research procedures were approved by the Institutional Animal Care Committee at Nankai University and Tianjin Institute of Pharmaceutical Research New Drug Evaluation Co., Ltd. (IACUC number: 2016032102), Tianjin, China. Efforts were made to minimize animal suffering and to reduce the number of animals used.
Results
General Features of Putative PSgcs in E. aerogenes
In this study, genomes were sequenced from two reference E. aerogenes strains from the American Type Culture Collection (ATCC), two from the German Collection of Microorganisms and Cultures (DSMZ), and 19 clinical strains that were isolated from the Shanghai Municipal Center for Disease Control and Prevention (SCDC) during 2012–2015 (Supplementary Table 1). Putative PSgc regions of E. aerogenes were first located by us, and a total of eight different putative PSgc types were found. All eight putative PSgcs mapped between two housekeeping genes, galF and gnd, and ranged in size from 11,194 bp to 23,479 bp, with all genes transcribed from galF and gnd. The average % GC content of all putative PSgcs is ∼40%, which is significantly lower than for the E. aerogenes genome as a whole (55%). The PSgc allocation and accession numbers for all type strains are summarized in Supplementary Table 2.
In each PSgc, the conserved gene set cpsACP-wzi-wza-wzb-wzc was located at the 5′ end. Further, an initial glycosyltransferase gene (wbaP or wcaJ) was present, whose products have been characterized as initial glycosyl transferases (ITs) that transfer galactose-1-phosphate or glucose-1-phosphate, respectively, to undecaprenol phosphate to initiate polysaccharide synthesis (Liu et al., 1993; Patel et al., 2012). In several PSgc types, manCB, which is responsible for the formation of GDP-D-Man (Liu et al., 2008; Perepelov et al., 2015), rmlAB/vioAB for dDTP-D-Qui4NAc (Wang et al., 2007), and glf for UDP-Galf (Nassau et al., 1996), were also presented. Biosynthetic pathways for all putative rare sugars present in E. aerogenes surface polysaccharides are indicated in Supplementary Figure 1. Meanwhile, each PSgc possessed non-initial glycosyltransferase genes (GTs), and oligosaccharide unit processing genes were observed in most PSgc types (Figure 1). Characteristics of open reading frames (ORFs) for all putative PSgc types are summarized in Supplementary Table 3.
FIGURE 1. Schematic diagram of the eight putative PSgc types identified from the 23 Enterobacter aerogenes strain genomes sequenced in this study. Genes are represented by arrows and colored according to the gene key at the bottom with gene names indicated above each arrow.
A conserved 39 bp JUMPstar sequence has been reported to be present in the non-coding regions upstream of PSgcs in several bacteria, and it was proposed to play two possible roles, regulating the transcription of PSgcs and involvement in the recombination of sugar pathway and assembly genes between loci (Hu et al., 2013). We aligned the non-coding regions of each PSgc between cpsACP and wzi using MEGA 4 software and identified homologous JUMPstar sequences that were also present in these gene clusters (Supplementary Figure 2).
The results outlined above indicate that the genetic region between galF and gnd should be a PSgc, and a potential molecular serotyping system could be presented based on the genetic diversity of PSgc types.
The Conventional Antigenic Scheme Correlated Well With the Molecular Serotyping System Based on PSgc Genetic Variation
Antisera were prepared using each type strain that represented PSgc types 1–8 and tested against all isolates. The homologous and heterologous titers of these antisera are summarized in Table 1. Generally, homologous titers were high, varying from 320 to 1280. However, none of the antisera were specific for only their homologous strains, with the exception of PSgc6 antiserum, and each antiserum resulted in unique cross-agglutination with one or more of the other antigens. Thus, the specific antisera needed to be absorbed with the corresponding heterologous antigens. Each antiserum agglutinated with only the corresponding homologous antigen after absorption, with serum titers ranging from 160 to 640 (Table 2). The eight absorbed antisera were then tested against all available strains, and each produced specific agglutination only against the strains that carried the same PSgc. These results suggested that the traditional immunoassay system was perfectly consistent with the molecular serotyping scheme that we developed based on PSgc diversity.
TABLE 1. Homologous and heterologous agglutinin titers of E. aerogenes antisera.
TABLE 2. Agglutination of absorbed E. aerogenes antisera.
Development of a Multiplexed Luminex-Based Array
Compared with nucleotide sugar precursor synthesis genes, glycosyltransferase genes and processing genes are very heterogeneous among serotypes owing to the wide range of possible sugar linkages, and they thereby provide the potential as genetic targets for serotypes (Li and Reeves, 2000). Furthermore, in contrast to the glycosyltransferase gene, the wzy gene encoding polymerase and the wzx gene encoding flippase are much more highly serotype determinative (Ballmer et al., 2007). Consequently, all primers and probes, except those for PSgc4, were designed based on wzy. For PSgc4, as it processes no wzy or wzx gene and as all three glycosyltransferase genes give no identity to any sequences via BLAST search, we selected GT2 as the specific gene (Table 3).
TABLE 3. Primers and probes used in the multiplex Luminex-based assay.
Twenty-three E. aerogenes strains carrying PSgc1-8 and other pathogenic bacteria, including Escherichia coli (n = 2), Salmonella spp. (n = 2), Shigella spp. (n = 2), Klebsiella pneumoniae (n = 2), Klebsiella oxytoca (n = 2), Enterobacter cloacae (n = 1), Enterobacter sakazakii (n = 1), Vibrio cholerae (n = 2), Vibrio parahaemolyticus (n = 1), and Citrobacter freundii (n = 1), were used to assess the specificity of our multiplexed Luminex-based array. Multiplex PCR produced only a single band of the expected size for each type strain, and no non-specific amplicons were detected (Figure 2A). The subsequent Luminex-based array analyses indicated that each PSgc-specific probe detected the homologous strains correctly, and signals corresponding to heterologies or other pathogens were not observed (Figure 2B). A double-blind test with 97 clinical isolates that were obtained from the SCDC (i.e., hereafter referred to as Shanghai isolates) was performed. The results from these tests indicated that each isolate belonged to one of the eight PSgc types (PSgc1-8). Among these isolates, 32 (33%) contained PSgc4, 16 (16%) contained PSgc7, 12 (12%) contained PSgc 6, and 11 (11%) contained either PSgc2 or PSgc8, with other PSgc types each corresponding to <10% of the isolates. This result was confirmed to be correct by ABI sequencing of each single PCR amplicon. The limit of detection for each of the target strains was determined by examining six serial 10-fold dilutions (10 ng to 0.1 pg) of each type strain’s genomic DNA. As a result, the sensitivity of our assay using genomic DNA was 0.1 ng for each strain.
FIGURE 2. Results of multiplex Luminex-based array. (A) Multiplex PCR of type strains representing each PSgc type. M, DNA Marker I (600, 500, 400, 300, 200, and 100 bp); 1, G2351; 2, G5305; 3, G5306; 4, G5307; 5, G5308; 6, G5310; 7, G5313; 8, G5319. (B) Specific detection based on the Luminex-based array. Biotin-labeled PCR products were separated by probe-coupled beads, and the hybridization signals are presented in terms of median fluorescence intensity (MFI) on the y-axis. Each sample representing the corresponding type strain is indicated on the x-axis. Data are combined from three independent experiments with error bars representing the standard deviation.
Molecular Serotyping Allocation of E. aerogenes Strains Based on Publicly Available Genome Sequences
Putative PSgc regions were extracted and identified from 86 E. aerogenes genomes that were obtained from the GenBank database. Thirty-two strains contained one of the seven PSgc types described above (PSgc3 was only found in G5306). In addition, another seven PSgc types (PSgc9-15) were distributed among the remaining 54 isolates (Figure 3 and Supplementary Table 2). Details for all ORFs corresponding to these seven types are provided in Supplementary Table 3. Thus, our combined analyses identified a total of 15 PSgc types within the E. aerogenes molecular serotyping scheme. Among the 86 isolates, 37 (43%) contained PSgc10, followed by PSgc4 (12 of 86) and PSgc7 (7 of 86). Since the majority (74%) of the 86 strains were isolated from the United States and the remainder originated from other or unknown regions, the strains are tentatively referred to as “United States isolates” in the following discussion.
FIGURE 3. Schematic diagram of the seven novel putative PSgc types that were identified in 54 strains whose genomes were publicly available. Genes are represented by arrows and colored according to the gene key at the bottom with gene names indicated above each arrow.
Within the 15 PSgcs, several anomalies were observed, including the following: (i) neither wzx nor wzy was annotated in PSgc4, PSgc14, or PSgc15, there was no wzx gene in PSgc7, and there was no wzy gene in PSgc2 or PSgc11. This atypical feature has been reported in Klebsiella capsular K11 and K34 without the wzx gene and in K29 and K50 without the wzy gene (Pan et al., 2015). In Salmonella serogroups A, B, and D1, the wzy genes also map far from the O-antigen gene cluster (Wang et al., 2002). Thus, it is likely that in the isolates with the above PSgc types, the functional wzx/wzy gene(s) are located outside of the PSgc region. However, they are difficult to annotate because of their low similarities to their counterparts, and further studies should be performed to confirm our hypothesis. (ii) In PSgc14, only manB was annotated. To our knowledge, this is the first evidence that manC and manB do not coexist in a gene cluster. (iii) It is noteworthy that the ugd and gnd genes are transcribed in the opposite direction in PSgc14 and that ugd is located at the 3′ end, meaning that there is probably a promoter at the end of the gene cluster.
Discussion
In the present study, we developed a molecular serotyping system for Enterobacter aerogenes that is based on the PSgc region, which correlated perfectly with the conventional antigenic scheme. Hence, from the serotyping perspective, molecular PSgc-based typing and immunological serology are equally valid.
Differences in serotype distribution typically occur between geographical regions. For instance, among the genus Shigella, S. flexneri is globally and traditionally isolated frequently in resource-poor countries, with predominant serotypes including S. flexneri 2a, 3a, and 6, whereas S. sonnei is detected most commonly in high-income regions (Thompson et al., 2015). Using the molecular serotyping scheme presented here, we revealed the occurrence of diverse clusters of E. aerogenes serotypes between far-distant regions (Supplementary Figure 3). Strains containing the PSgc4 and PSgc7 types were the dominant serotypes among Shanghai isolates, while those with PSgc10 comprised the majority of United States isolates. It has been proposed that as one of the important components exposed on the cell surface, the bacterial surface polysaccharide is highly immunogenic and therefore subject to intense selection by the host immune system (Reeves and Wang, 2002). This selection probably accounts for the maintenance of surface polysaccharide diversity and offers certain serotype(s) selective advantages in specific niches of different hosts. Overall, our findings suggest an important distinction for clinical microbiologists, infectious disease specialists, and infection control agencies, although more data should be gathered to confirm these results.
In this study, we prepared antisera for each PSgc type strain and performed agglutination tests, providing evidence of the reliability of our molecular serotyping scheme. However, several disadvantages exist in the development and application of serotyping antisera and should be noted. First, the development of an antigenic scheme is sophisticated and time consuming (Orskov et al., 1977; Pan et al., 2015). Second, serotyping is labor-intensive, and only a few reference laboratories have the capacity to perform the analyses necessary to identify serotypes (Ballmer et al., 2007; Lacher et al., 2014). In addition, conventional serological assays that are based on agglutination reactions are also limited by the risks of the high prevalence of non-typeable isolates due to the expression of novel polysaccharide forms, which is common in clinical isolates (Jenney et al., 2006; Thrane et al., 2016; Guo et al., 2017). The development of NGS technologies and automated data analysis pipelines has led to a number of proof-of-concept studies. These studies have suggested the possibility of increased resolution of the existing typing methods and potential advantages in public health microbiology for outbreak investigation and surveillance (Dallman et al., 2015; Joensen et al., 2015; Zhang et al., 2015; Thrane et al., 2016). Moreover, a few practical genomics-based approaches have been reported for the typing of bacteria and pathogen surveillance (den Bakker et al., 2014; Dallman et al., 2015; Kwong et al., 2016). However, a few shortcomings of NGS seem insurmountable for the routine detection and identification of pathogens over a short period of time including (i) the need for pure isolate cultures to obtain enough quality genomic DNA (usually 50 ng to 1 μg depending on different sequencing platforms); (ii) the lengthy turnaround time for library preparation and NGS platform sequencing runs; and (iii) the relatively high cost per isolate genome sequence. Consequently, typing methods that are based on PSgc-specific gene(s) are likely to be widely applied and economically viable alternatives for the foreseeable future. These methodologies, including RFLP, multiplex PCR, high throughput real-time PCR, and microarrays, have been widely utilized to identify serotypes based on sero-specific genes in major pathogens, including E. coli (Coimbra et al., 2000; Botkin et al., 2012), Streptococcus pneumoniae (Azzari et al., 2010), Salmonella spp. (Guo et al., 2013) and Streptococcus suis (Bai et al., 2015). The Luminex-based array system is a multiplexed microsphere-based suspension system that offers a molecular diagnostic platform and provides an open and attractive approach for simultaneous, high-throughput and multiplex detection of up to 100 targets in protein and nucleic acid studies (Yan et al., 2017). It has been approved by the US FDA for clinical diagnosis and has been used in various applications (Dunbar, 2006; Liu et al., 2011; Glushakova et al., 2015; Silbereisen et al., 2015). Although the assay presented here only targeted E. aerogenes strains of types PSgc1-8, which appear to be more common in Shanghai and its surrounding environs, thus highlighting a limitation of the analysis (i.e., strain availability), the detection range of our assay could easily be expanded, as PSgc data can be obtained without limitation. Indeed, evaluation and revisiting of the serotyping system presented here should be performed as more isolates become available in the future. Regardless, our method shows significant potential for the characterization and epidemiological surveillance of this emerging pathogen.
We propose that all of the putative PSgc types identified in this study are responsible for the biosynthesis of capsular (K) antigens of E. aerogenes. The reasons for this supposition are as follows: (i) a conserved gene set, cpsACP-wzi-wza-wzb-wzc, is located at the 5′ end of each PSgc, which shares the typical features of most K antigen gene clusters in the genus Klebsiella, which is closely related to E. aerogenes, and the group 1 capsule gene clusters of E. coli (Whitfield, 2006; Pan et al., 2015). (ii) In the genome of E. aerogenes CAV1320, which contained PSgc5, another genetic region (loci: ABY61_23370 to ABY61_23395) containing six genes (wzm, wzt, glf, and GT1-3) was found upstream of hisI, which resembled an O antigen gene cluster (Supplementary Figure 4). This region has been purported to have been introduced into an E. coli strain by an IS element, forming the O62 serotype in our recent study (Hou et al., 2017). We also examined the putative O antigen gene cluster among the strains for which whole genomes were available. Ten types of putative O-antigen gene clusters were characterized and distributed among 105 genomes, except for another four incomplete genomes due to poor sequencing (data not shown).
Seven of eight antisera gave slight serological cross-reactivity with heterologous antigens before absorption. It is hard to predict whether this is attributable to the similar polysaccharides, as no chemical structure has ever been reported. However, the PSgc types of each strain are quite different, meaning probably that no similar structures are shared among them. We therefore propose that the weak cross-reaction is due to the presence of common epitopes of other cell surface components. More detailed studies of each PSgc and the corresponding polysaccharide structures should be performed, and the correlation between them will be necessary to completely elucidate the heredity and evolution of E. aerogenes. In this regard, the results presented here provide a valuable framework from which to further assess the evolution of E. aerogenes surface polysaccharide structures in response to host and bacteriophage interactions.
Author Contributions
XG, BL, and LF conceived the project. XG and MW prepared the strain samples, preformed the sequence analyses, and developed the molecular serotyping system. LW conducted the bioinformatics analyses. YW and TC developed the conventional antigenic scheme. PW developed the multiplexed Luminex-based array. MC performed the double-blind test. XG, BL, and LF prepared the manuscript. All authors read and approved the final manuscript.
Funding
This work was supported by National Key Programs for Infectious Diseases of China (2017ZX10303405-001 and 2017ZX10104002-001-006), a Tianjin Municipal Natural Science Foundation Grant (17JCYBJC24300), a National Natural Science Foundation of China (NSFC) Key Program Grant (31530083), NSFC General Program Grants (81471904, 81772148, and 31470194), and a National Special Project on Research and Development of Key Biosafety Technologies Grant.
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.00501/full#supplementary-material
Footnotes
References
Altschul, S. F., Madden, T. L., Schäffer, A. A., Zhang, J., Zhang, Z., Miller, W., et al. (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25, 3398–3402. doi: 10.1093/nar/25.17.3389
Azzari, C., Moriondo, M., Indolfi, G., Cortimiglia, M., Canessa, C., Becciolini, L., et al. (2010). Realtime PCR is more sensitive than multiplex PCR for diagnosis and serotyping in children with culture negative pneumococcal invasive disease. PLoS One 5:e9282. doi: 10.1371/journal.pone.0009282
Bai, X., Liu, Z., Ji, S., Gottschalk, M., Zheng, H., and Xu, J. (2015). Simultaneous detection of 33 Streptococcus suis serotypes using the luminex xTAG® assayTM. J. Microbiol. Methods 117, 95–99. doi: 10.1016/j.mimet.2015.07.018
Ballmer, K., Korczak, B. M., Kuhnert, P., Slickers, P., Ehricht, R., and Hächler, H. (2007). Fast DNA serotyping of Escherichia coli by use of an oligonucleotide microarray. J. Clin. Microbiol. 45, 370–379. doi: 10.1128/JCM.01361-06
Botkin, D. J., Galli, L., Sankarapani, V., Soler, M., Rivas, M., and Torres, A. G. (2012). Development of a multiplex PCR assay for detection of Shiga toxin-producing Escherichia coli, enterohemorrhagic E. coli, and enteropathogenic E. coli strains. Front. Cell. Infect. Microbiol. 2:8. doi: 10.3389/fcimb.2012.00008
Chen, M., Guo, D., Wong, H. C., Zhang, X., Liu, F., Chen, H., et al. (2012). Development of O-serogroup specific PCR assay for detection and identification of Vibrio parahaemolyticus. Int. J. Food Microbiol. 159, 122–129. doi: 10.1016/j.ijfoodmicro.2012.08.012
Chen, R., Bi, Y., Yang, G., Liu, Z., Liu, Z., Zeng, B., et al. (2010). Development of a fluorescent microsphere-based multiplex assay for simultaneous rapid detection of Mycobacterium tuberculosis complex and differentiation of M. tuberculosis and M. bovis in clinical samples. Diagn. Mol. Pathol. 19, 172–179. doi: 10.1097/PDM.0b013e3181d8c241
Chen, Y. G., Zhang, Y., Yu, Y. S., Qu, T. T., Wei, Z. Q., Shen, P., et al. (2008). In vivo development of carbapenem resistance in clinical isolates of Enterobacter aerogenes producing multiple beta-lactamases. Int. J. Antimicrob. Agents 32, 302–307. doi: 10.1016/j.ijantimicag.2008.02.014
Chen, Z., Li, H., Feng, J., Li, Y., Chen, X., Guo, X., et al. (2015). NDM-1 encoded by a pNDM-BJ01-like plasmid p3SP-NDM in clinical Enterobacter aerogenes. Front. Microbiol. 14:294. doi: 10.3389/fmicb.2015.00294
Chevalier, J., Mulfinger, C., Garnotel, E., Nicolas, P., Davin-Régli, A., and Pagès, J. M. (2008). Identification and evolution of drug efflux pump in clinical Enterobacter aerogenes strains isolated in 1995 and 2003. PLoS One 3:e3203. doi: 10.1371/journal.pone.0003203
Coimbra, R. S., Grimont, F., Lenormand, P., Burguière, P., Beutin, L., and Grimont, P. A. (2000). Identification of Escherichia coli O-serogroups by restriction of the amplified O-antigen gene cluster (rfb-RFLP). Res. Microbiol. 151, 639–654. doi: 10.1016/S0923-2508(00)00134-0
Dallman, T. J., Byrne, L., Ashton, P. M., Cowley, L. A., Perry, N. T., Adak, G., et al. (2015). Whole-genome sequencing for national surveillance of Shiga toxin-producing Escherichia coli O157. Clin. Infect. Dis. 61, 305–312. doi: 10.1093/cid/civ318
De Gheldre, Y., Maes, N., Rost, F., De Ryck, R., Clevenbergh, P., Vincent, J. L., et al. (1997). Molecular epidemiology of an outbreak of multidrug-resistant Enterobacter aerogenes infections and in vivo emergence of imipenem resistance. J. Clin. Microbiol. 35, 152–160.
den Bakker, H. C., Allard, M. W., Bopp, D., Brown, E. W., Fontana, J., Iqbal, Z., et al. (2014). Rapid whole-genome sequencing for surveillance of Salmonella enterica serovar enteritidis. Emerg. Infect. Dis. 20, 1306–1314. doi: 10.3201/eid2008.131399
Dunbar, S. A. (2006). Applications of Luminex xMAP technology for rapid, high-throughput multiplexed nucleic acid detection. Clin. Chim. Acta 363, 71–82. doi: 10.1016/cccn.2005.06.023
Glushakova, L. G., Bradley, A., Bradley, K. M., Alto, B. W., Hoshika, S., Hutter, D., et al. (2015). High-throughput multiplexed xMAP Luminex array panel for detection of twenty two medically important mosquito-borne arboviruses based on innovations in synthetic biology. J. Virol. Methods 214, 60–74. doi: 10.1016/j.jviromet.2015.01.003
Grattard, F., Pozzetto, B., Tabard, L., Petit, M., Ros, A., and Gaudin, O. G. (1995). Characterization of nosocomial strains of Enterobacter aerogenes by arbitrarily primed-PCR analysis and ribotyping. Infect. Control Hosp. Epidemiol. 16, 224–230. doi: 10.2307/30140982
Guibourdenche, M., Roggentin, P., Mikoleit, M., Fields, P. I., Bockemuhl, J., Grimont, P. A., et al. (2010). Supplement 2003-2007 (no. 47) to the White-Kauffmann-Le Minor scheme. Res. Microbiol. 161, 26–29. doi: 10.1016/j.resmic.2009.10.002
Guo, D., Liu, B., Liu, F., Cao, B., Chen, M., Hao, X., et al. (2013). Development of a DNA microarray for molecular identification of all 46 Salmonella O serogroups. Appl. Environ. Microbiol. 79, 3392–3399. doi: 10.1128/AEM.00225-13
Guo, X., Liu, B., Chen, M., Wang, Y., Wang, L., Chen, H., et al. (2017). Genetic and serological identification of three Vibrio parahaemolyticus strains as candidates for novel provisional O serotypes. Int. J. Food Microbiol. 245, 53–58. doi: 10.1016/j.ijfoodmicro.2017.01.010
Henikoff, S., Henikoff, J. G., Alford, W. J., and Pietrokovski, S. (1995). Automated construction and graphical presentation of protein blocks from unaligned sequences. Gene 163, GC17–GC26. doi: 10.1016/0378-1119(95)00486-P
Hou, X., Perepelov, A. V., Guo, X., Senchenkova, S. N., Shashkov, A. S., Liu, B., et al. (2017). A gene cluster at an unusual chromosomal location responsible for the novel O-antigen synthesis in Escherichia coli O62 by the ABC transporter-dependent pathway. Glycobiology 27, 669–676. doi: 10.1093/glycob/cwx030
Hu, D., Liu, B., Dijkshoorn, L., Wang, L., and Reeves, P. R. (2013). Diversity in the major polysaccharide antigen of Acinetobacter baumannii assessed by DNA sequencing, and development of a molecular serotyping scheme. PLoS One 8:e70329. doi: 10.1371/journal.pone.0070329
Jalaluddin, S., Devaster, J. M., Scheen, R., Gerard, M., and Butzler, J. P. (1998). Molecular epidemiological study of nosocomial Enterobacter aerogenes isolates in a Belgian hospital. J. Clin. Microbiol. 36, 1846–1852.
Jenney, A. W., Clements, A., Farn, J. L., Wijburg, O. L., McGlinchey, A., Spelman, D. W., et al. (2006). Seroepidemiology of Klebsiella pneumoniae in an Australian tertiary hospital and its implications for vaccine development. J. Clin. Microbiol. 44, 102–107. doi: 10.1128/JCM.44.1.102-107.2006
Joensen, K. G., Tetzschner, A. M., Iguchi, A., Aarestrup, F. M., and Scheutz, F. (2015). Rapid and easy in silico serotyping of Escherichia coli isolates by use of whole-genome sequencing data. J. Clin. Microbiol. 53, 2410–2426. doi: 10.1128/JCM.00008-15
Konradsen, H. B., and Pneumococcus Reference laboratories in Europe (2005). Validation of serotyping of Streptococcus pneumoniae in Europe. Vaccine 23, 1368–1373. doi: 10.1016/j.vaccine.2004.09.011
Kwong, J. C., Mercoulia, K., Tomita, T., Easton, M., Li, H. Y., Bulach, D. M., et al. (2016). Prospective whole-genome sequencing enhances national surveillance of Listeria monocytogenes. J. Clin. Microbiol. 54, 333–342. doi: 10.1128/JCM.02344-15
Lacher, D. W., Gangiredla, J., Jackson, S. A., Elkins, C. A., and Feng, P. C. (2014). Novel microarray design for molecular serotyping of shiga toxin producing Escherichia coli strains isolated from fresh produce. Appl. Environ. Microbiol. 80, 4677–4682. doi: 10.1128/AEM.01049-14
Langley, J. M., Hanakowski, M., and Leblanc, J. C. (2001). Unique epidemiology of nosocomial urinary tract infection in children. Am. J. Infect. Control 29, 94–98. doi: 10.1067/mic.2001.111537
Lavigne, J. P., Sotto, A., Nicolas-Chanoine, M. H., Bouziges, N., Bourg, G., Davin-Regli, A., et al. (2012). Membrane permeability, a pivotal function involved in antibiotic resistance and virulence in Enterobacter aerogenes clinical isolates. Clin. Microbiol. Infect. 18, 539–545. doi: 10.1111/j.1469-0691.2011.03607.x
Lavigne, J. P., Sotto, A., Nicolas-Chanoine, M. H., Bouziges, N., Pages, J. M., and Davin-Regli, A. (2013). An adaptive response of Enterobacter aerogenes to imipenem: regulation of porin balance in clinical isolates. Int. J. Antimicrob. Agents 41, 130–136. doi: 10.1016/j.ijantimicag.2012.10.010
Li, Q., and Reeves, P. R. (2000). Genetic variation of dTDP-L-rhamnose pathway genes in Salmonella enterica. Microbiology 136, 2291–2307. doi: 10.1099/00221287-146-9-2291
Liu, B., Knirel, Y. A., Feng, L., Perepelov, A. V., Senchenkova, S. N., Reeves, P. R., et al. (2014). Structural diversity in Salmonella O antigens and its genetic basis. FEMS Microbiol. Rev. 38, 56–89. doi: 10.1111/1574-6976.12034
Liu, B., Knirel, Y. A., Feng, L., Perepelov, A. V., Senchenkova, S. N., Wang, Q., et al. (2008). Structure and genetics of Shigella O antigens. FEMS Microbiol. Rev. 32, 627–653. doi: 10.1111/j.1574-6976.2008.00114.x
Liu, D., Haase, A. M., Lindgvist, L., Lindberg, A. A., and Reeves, P. R. (1993). Glycosyl transferases of O-antigen biosynthesis in Salmonella enterica: identification and characterization of transferase genes of groups B, C2, and E1. J. Bacteriol. 175, 3408–3413. doi: 10.1128/jb.175.11.3408-3413.1993
Liu, J., Kibiki, G., Maro, V., Maro, A., Kumburu, H., Swai, N., et al. (2011). Multiplex reverse transcription PCR Luminex assay for detection and quantitation of viral agents of gastroenteritis. J. Clin. Virol. 50, 308–313. doi: 10.1016/j.jcv.2010.12.009
Nassau, P. M., Martin, S. L., Brown, R. E., Weston, A., Monsey, D., McNeil, M. R., et al. (1996). Galactofuranose biosynthesis in Escherichia coli K-12: identification and cloning of UDP-galactopyranose mutase. J. Bacteriol. 178, 1047–1052. doi: 10.1128/jb.178.4.1047-1052.1996
Orskov, I., Orskov, F., Jann, B., and Jann, K. (1977). Serology, chemistry, and genetics of O and K antigens of Escherichia coli. Bacteriol. Rev. 41, 667–710.
Pan, Y. J., Lin, T. L., Chen, C. T., Chen, Y. Y., Hsieh, P. F., Hsu, C. R., et al. (2015). Genetic analysis of capsular polysaccharide synthesis gene clusters in 79 capsular types of Klebsiella spp. Sci. Rep. 5:15573. doi: 10.1038/srep15573
Patel, K. B., Toh, E., Fernandez, X. B., Hanuszkiewicz, A., Hardy, G. G., Brun, Y. V., et al. (2012). Functional characterization of UDP-glucose: undecaprenyl-phosphate glucose-1-phosphate transferases of Escherichia coli and Caulobacter crescentus. J. Bacteriol. 194, 2646–2657. doi: 10.1128/JB.06052-11
Perepelov, A. V., Shashkov, A. S., Guo, X., Filatov, A. V., Weintraub, A., Widmalm, G., et al. (2015). Structure and genetics of the O-antigen of Escherichia coli O169 related to the O-antigen of Shigella boydii type 6. Carbohydr. Res. 23, 46–50. doi: 10.1016/j.carres.2015.05.016
Rahn, A., Beis, K., Naismith, J. H., and Whitfield, C. (2003). A novel outer membrane protein, Wzi, is involved in surface assembly of the Escherichia coli K30 group 1 capsule. J. Bacteriol. 185, 5882–5890. doi: 10.1128/JB.185.19.5882-5890.2003
Reeves, P. P., and Wang, L. (2002). Genomic organization of LPS-specific loci. Curr. Top. Microbiol. Immunol. 264, 109–135.
Rutherford, K., Parkhill, J., Crook, J., Horsnell, T., Rice, P., Rajandream, M. A., et al. (2000). Artemis: sequence visualisation and annotation. Bioinformatics 16, 944–945. doi: 10.1093/bioinformatics/16.10.944
Silbereisen, A., Tamborrini, M., Wittwer, M., Schürch, N., and Pluschke, G. (2015). Development of a bead-based Luminex assay using lipopolysaccharide specific monoclonal antibodies to detect biological threats from Brucella species. BMC Microbiol. 15:198. doi: 10.1186/s12866-015-0534-1
Thompson, C. N., Duy, P. T., and Baker, S. (2015). The rising dominance of Shigella sonnei: an intercontinental shift in the etiology of bacillary dysentery. PLoS Negl. Trop. Dis. 9:e0003708. doi: 10.1371/journal.pntd.0003708
Thrane, S. W., Taylor, V. L., Lund, O., Lam, J. S., and Jelsbak, L. (2016). Application of WGS data for O-specific antigen analysis and in silico serotyping of Pseudomonas aeruginosa isolates. J. Clin. Microbiol. 54, 1782–1788. doi: 10.1128/JCM.00349-16
Wang, L., Andrianopoulos, K., Liu, D., Popoff, M. Y., and Reeves, P. R. (2002). Extensive variation in the O-antigen gene cluster within one Salmonella enterica serogroup reveals an unexpected complex history. J. Bacteriol. 184, 1669–1677. doi: 10.1128/JB.184.6.1669-1677.2002
Wang, Q., Ruan, X., Wei, D., Hu, Z., Wu, L., Yu, T., et al. (2010). Development of a serogroup-specific multiplex PCR assay to detect a set of Escherichia coli serogroups based on the identification of their O-antigen gene clusters. Mol. Cell. Probes 24, 286–290. doi: 10.1016/j.mcp.2010.06.002
Wang, Y., Xu, Y., Perepelov, A. V., Qi, Y., Kinrel, Y. A., Wang, L., et al. (2007). Biochemical characterization of dTDP-D-Qui4N and dTDP-D-Qui4NAc biosynthetic pathways in Shigella dysenteriae type 7 and Escherichia coli O7. J. Baceriol. 189, 8626–8635. doi: 10.1128/JB.00777-07
Whitfield, C. (2006). Biosynthesis and Assembly of capsular polysaccharides in Escherichia coli. Annu. Rev. Biochem. 75, 39–68. doi: 10.1146/annurev.biochem.75.103004.142545
Yan, Y., Luo, J. Y., Chen, Y., Wang, H. H., Zhu, G. Y., He, P. Y., et al. (2017). A multiplex liquid-chip assay based on Luminex xMAP technology for simultaneous detection of six common respiratory virus. Oncotarget 8, 96913–96923. doi: 10.18632/oncotarget.18533
Keywords: Enterobacter aerogenes, polysaccharide gene cluster, molecular serotyping system, antigenic scheme, Luminex-based array
Citation: Guo X, Wang M, Wang L, Wang Y, Chen T, Wu P, Chen M, Liu B and Feng L (2018) Establishment of a Molecular Serotyping Scheme and a Multiplexed Luminex-Based Array for Enterobacter aerogenes. Front. Microbiol. 9:501. doi: 10.3389/fmicb.2018.00501
Received: 23 November 2017; Accepted: 05 March 2018;
Published: 19 March 2018.
Edited by:
Dongsheng Zhou, Beijing Institute of Microbiology and Epidemiology, ChinaReviewed by:
Guang Zhao, Qingdao Institute of Bioenergy and Bioprocess Technology (CAS), ChinaHan Zheng, Chinese Center for Disease Control and Prevention, China
Quan Wang, Tianjin Medical University, China
Xin Deng, City University of Hong Kong, Hong Kong
Copyright © 2018 Guo, Wang, Wang, Wang, Chen, Wu, Chen, Liu and Feng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Bin Liu, liubin1981@nankai.edu.cn Lu Feng, fenglu63@nankai.edu.cn
†These authors have contributed equally to this work.