Abstract
SHV β-lactamases, including SHV extended-spectrum β-lactamases, are widespread throughout the world, and confer a broad spectrum of resistance to antibiotic drugs. Mutations ranging from single base-pair substitutions to small deletions within blaSHV often result in diminished activity and an increased susceptibility to β-lactamase inhibitors. Here, we collected 1,320 clinical isolates from three hospitals in Shanghai. We developed a novel oligonucleotide microarray to detect mutations in the blaSHV gene, and validated the array data by direct sequencing. Sixty-two of the 1,320 isolates carried the blaSHV gene. The genotypes of these 62 isolates were successfully called by the microarray and were consistent with the genotypes identified by bidirectional sequencing. Sixteen different blaSHV alleles were identified. The SHV-1 variant was the most frequent (32.26%), followed by SHV-11 (27.42%) and SHV-12 (25.81%). Of the 62 isolates, 12 contained two different blaSHV alleles. Our microarray significantly facilitated the identification of blaSHV variants, which makes it an attractive option for the detection of SHV variants in clinical laboratories.
Introduction
Conventional laboratory methods, which are phenotype-based antimicrobial susceptibility tests, identify ESBL producers by detecting the presence of the ESBLs themselves. However, some bacteria that are susceptible to β-lactam antibiotics in vitro by routine methods show resistance to them in the clinic. Therefore, the prevalence of SHV-type ESBLs may be underestimated, and tests to identify the presence of an ESBL are not sufficiently reliable to determine ESBL producers definitively.11,36 Furthermore, multiple blaSHV gene copies were observed in most, if not all, strains containing the blaSHV gene16,38; thus, it is important to determine the different alleles of the blaSHV gene within individual isolates. Isoelectric focusing (IEF) can identify the coexistence of SHV and TEM in one isolate. 18 However, many β-lactamases possess an identical pI, and thus the identification of coexisting SHV variants by IEF is far from unambiguous. A high mutation rate allows fast adaptation and enables bacteria to evolve rapidly in one generation. Genes encoding β-lactamases are usually highly polymorphic because of antibiotic selective pressure and high mutation rates. In contrast to phenotype-based tests, variants of the blaSHV gene can be accurately identified by genotyping polymorphic sites in the blaSHV gene.
Various approaches to investigate single-nucleotide polymorphisms (SNPs) have been applied to identify mutations in β-lactamase genes, including single-base extension, 16 polymerase chain reaction (PCR)–restriction fragment length polymorphism, 8 ligase chain reaction, 21 matrix-assisted laser desorption ionization–time-of-flight mass spectrometry, 8 and microarray technology.17,5 Among these methods, only microarray technology provides the ability to genotype SNPs with a high degree of multiplexing. Diagnostic microarrays have been developed for the detection of microorganisms,6,22,23 antibiotic-resistance genes,3,7,13, 18 and drug-metabolism genes. 20 However, an assay that is able to genotype more than 50 mutations of the blaSHV gene simultaneously and discriminate heterozygotes has not yet been described.
In the present study, we developed a microarray to sample the blaSHV gene in clinical samples. This novel oligonucleotide microarray was used to determine 54 mutations in the blaSHV gene. It was able to distinguish 67 variant alleles of the blaSHV gene and discriminate different alleles within an isolate.
Materials and Methods
Bacterial strains and DNA extraction
Strain 46112 and another two K. pneumoniae strains, containing SHV-1, SHV-12, and SHV-18, respectively, were used as controls. These strains were validated by direct sequencing. A total of 1,320 nonreplicate isolates were collected from Ruijin Hospital, Dongfang Hospital, and Shanghai No. 1 Hospital. The isolates were identified using standard biochemical tests. 28 The clones were scraped from the solid culture medium, and plasmid DNA was extracted from them by alkaline lysis.
Primer design
Primers were designed using Primer3 (http://frodo.wi.mit.edu/cgi-bin/primer3/primer3_www_slow.cgi) on the basis of the published sequences of the blaSHV and rrlA genes from the GenBank database and published data on mutations in the blaSHV gene (www.lahey.org/studies/). The primer sequences are listed in Table 1.
Location of primers refers to the ATG translation start site of the blaSHV gene and the GTG translation start site of the rrlA gene.
Sequencing primer.
Universal primer.
Detection of the blaSHV gene
The blaSHV and rrlA genes were co-amplified using the primer pairs SHV-f/SHV-r and 23S-f/23S-r, respectively, in a single duplex reaction by touchdown PCR. 9 Amplification reactions were carried out in a total volume of 30 μl containing 0.3 mM of each deoxynucleoside triphosphate, 10 mM Tris–HCl (pH 8.3), 100 mM KCl, 2 mM MgCl2, 0.5 μM of primers SHV-f and SHV-r, 0.16 μM of primers 23S-f and 23S-r, 60 ng of DNA, and 1.2 U of Taq (Takara). Cycling conditions were as follows: 94°C for 3 min, followed by 10 cycles of 94°C for 30 sec, 63°C for 30 sec with a 0.5°C decrement of the annealing temperature per cycle, and 72°C for 40 sec, followed by 30 cycles of 94°C for 30 sec, 58°C for 30 s, and 72°C for 40 sec, and a final extension of 5 min at 72°C. The PCR products were then separated and detected by 1.8% agarose gel electrophoresis.
Amplification of the full-length blaSHV gene
To identify mutations in the blaSHV gene, a 30-μl PCR reaction was carried out to amplify the full-length blaSHV gene. The reaction mixture contained 0.3 mM of each deoxynucleoside triphosphate, 0.2 μM of primers SHV-Up and SHV-Low, 1× PCR reaction buffer, 2 mM MgCl2, 70 ng of DNA, and 0.8 U of Ex Taq (Takara). The cycling conditions were as follows: 94°C for 3 min; 8 cycles of 94°C for 30 sec, 69°C for 40 sec with a 0.5°C decrement of the annealing temperature per cycle, and 72°C for 55 sec; 30 cycles of 94°C for 30 sec, 65°C for 40 sec, and 72°C for 55 sec; and a final extension of 5 min at 72°C.
Labeling and fragmentation
Fluorescent-target DNA for hybridization was synthesized by amplifying and labeling the full-length blaSHV gene in a 90-μl PCR reaction. The reaction mixture contained 0.3 mM dATP, dGTP, and dTTP; 0.15 mM dCTP and CY3-dCTP (GE Healthcare); 0.2 μM of primers SHV-Up and SHV-Low; 1× PCR reaction buffer; 2 mM MgCl2; 100 ng of DNA; and 2.5 U of Ex Taq (Takara). The cycling conditions were described as above. The full-length amplicons were purified with a QIAquick PCR Purification kit (Qiagen) and were then fragmented for 15 min at 37°C using DNase I (0.0012 U/μg of DNA) in a 20-μl volume containing 1 mM Tris–HCl, 0.25 mM MgCl2, 0.05 mM CaCl2, and 2.5 U of alkaline phosphatase (Fermentas), followed by heat inactivation at 95°C for 6 min.
SHV oligonucleotide microarray
Several allele-specific probes, with a perfectly matched (PM) and a mismatched (MM) probe, were designed to detect each SNP. The MM probe acted as a control for nonspecific hybridization. There were two PM probes for each SNP: one identical to the wild-type allele and one identical to the mutant-type allele. The PM and MM probes for an SNP differed only at the central base within the probe sequence, that is, the SNP site itself, or were offset from the central base by one to two bases in either direction. More than two PM probes were designed for mutations located close together, to avoid an effect of polymorphic sites on hybridization efficiency. Furthermore, to facilitate data analysis, when multiple PM probes were used for a given SNP, each mutant-type PM probe was arranged in a group with the corresponding wild-type PM and MM probe. All probes were 13- to 30-nt long with a melting temperature in the range of 49°C–55°C (Oligo 6 software; Molecular Biology Insights). Each probe had a 15-thymidine nucleotide spacer and was modified by the addition of an amino group to the 5′-end to allow covalent attachment of the probe to the glass slide surface. All 174 amino-modified oligonucleotides with a 15-thymidine spacer were dissolved at 25 μM in 6× SSC buffer (pH 7.0) with 0.05% SDS, and were printed in triplicate as three subarrays on aldehyde-coated glass slides (Cell Associates, Inc.) using a OmniGrid™ 100 microarrayer (GeneMachine). A total of 174 oligonucleotides were printed on each glass slide in a 6.3-mm×8.5-mm grid. The spot spacing was 115 μm, and the spot size was 185 μm. After printing, the SHV microarray was dried at room temperature for 2 hr and then stored desiccated at room temperature.
Hybridization
Fluorescently labeled and fragmented DNAs were added to a hybridization solution containing 6×SSPE, 5% dimethyl sulfoxide, 0.1% Triton X-100, and 1 nM fluorescent control oligonucleotide in a final volume of 25 μl. The hybridization solution was heated to 95°C for 5 min, immediately placed on ice, and 20 μl was transferred to the subarray for hybridization at 50°C for 2 hr. After hybridization, microarrays were washed with 2× SSC and 0.1% SDS at 42°C for 5 min and then washed with 1× SSC and 0.1% SDS at 42°C for 5 min and with 0.5× SSC at room temperature for 5 min.
Data analysis
Microarrays were scanned using a GenePix 4000B scanner (Axon Instruments). The resulting images were processed to obtain hybridization signal intensity values using GenePix Pro 3.0 software (Axon Instruments). We set the mean of the signal intensities of the negative hybridization probe plus three standard deviations to be the cut-off value to block weak or unreliable signals. If the signal intensity value of the probe was less than the cut-off value, its net signal intensity was set as zero. After subtraction of local background signal, the net signal intensity of each probe that was spotted in triplicate was averaged. For each probe set related to one SNP, the average value of one PM was regarded as a specific signal value only if the ratio of the PM minus the MM to the PM plus the MM was higher than 0.3. The genotype could not be assigned if the ratios of all the PMs were lower than or equal to 0.3.
To assign the genotypes, the allelic fractions (AFs) between the signal value from one of the alleles and the sum of the signal values from both possible alleles at each SNP site were calculated from the normalized mean. 42
Sequencing
PCR products were incubated with shrimp alkaline phosphatase and exonuclease I for 1 hr at 37°C followed by 15 min at 80°C. 29 The PCR products were then sequenced in both directions using the primers SHV-s and SHV-r using Big Dye chemistry and run on ABI 3700 (Applied Biosystems). Given that the isolates carried two distinct alleles of the blaSHV gene, we used molecular haplotyping to delineate phase. The full-length PCR product was cloned into the pGEM-T vector (Promega). After transformation, a single colony was chosen for sequencing with T7 and SP6 primers. The difference between the reference sequence (SHV-1) and the cloned sequence provides the phased genotype for the plasmid.
Results
Prevalence of the blaSHV gene in clinical isolates from Shanghai hospitals
We collected 1,320 cultures of bacteria from blood, abscesses, peritoneal fluid, catheter tips, lung, sputum, and throat cultures from patients in three Shanghai hospitals: Ruijin Hospital, Dongfang Hospital, and Shanghai No. 1 Hospital. Some were used for in vitro antimicrobial susceptibility tests. To screen for isolates containing the blaSHV gene, we developed a PCR assay. In this duplex PCR assay, a 494-bp fragment amplified with the primer pair SHV-f and SHV-r indicated the presence of the blaSHV gene; a 340-bp fragment from the rrlA gene served as a positive control. The positive control was amplified from all 1,320 isolates (Supplementary Fig. S1; Supplementary Data are available online at www.liebertpub.com/mdr). The 494-bp fragment, indicative of the presence of the blaSHV gene, was observed in 62 strains, which equates to 4.7% of the clinical samples.
Of the 62 strains, 47 were Enterobacteriaceae, four were Pseudomonas, and five were Staphylococcus. The fact that both Gram-negative and Gram-positive bacteria carried the blaSHV gene implies that this gene is widely spread among divergent bacteria. Table 2 shows that the most common isolates containing the blaSHV gene were from K. pneumoniae (33.87%) followed by Escherichia coli (20.97%), and Enterobacter cloacae (14.52%). The prevalence of the blaSHV gene in Morganella morganii, Pseudomonas fluorescens, and Citrobacter freundii was higher than that in other species tested; however, because of the small sample size for these three species, these findings should be interpreted with caution.
Other species included Chryseobacterium meningosepticum, Enterococcus faecium, Haemophilus influenzae, Stenotrophomonas maltophilia, Streptococcus agalactiae, Proteus mirabilis, Staphylococcus aureus, Serratia marcescens, Staphylococcus warneri, Staphylococcus capitis, and Staphylococcus hominis.
SHV oligonucleotide microarray
We designed 174 oligonucleotide probes to detect five novel mutations and 49 published mutations in the blaSHV gene. The microarray is unique in its ability to discriminate 67 known blaSHV alleles, including five novel alleles.
For assay optimization, we sequenced three K. pneumoniae isolates, provided by Ruijin Hospital, containing SHV-1, SHV-12, and SHV-18 variants, respectively, for reference templates. Equal amounts of PCR products of any two reference templates were mixed and hybridized to the microarray to determine the AF limits. Furthermore, we also set a cut-off signal value and used the MM probes to block weak and nonspecific signals. By monitoring the allele discrimination performance of the microarray, we were able to optimize the hybridization and wash conditions. To assess assay quality, we used the microarray to genotype 62 isolates after complete discrimination of the mutations of the reference strains. All the mutations were successfully called, and the microarray results were in agreement with the genotypes obtained by bidirectional sequencing. These results indicated that our microarray was sufficiently sensitive to detect mutations in the blaSHV gene. A typical microarray image is shown in Supplementary Fig. S2 (Supplementary Data are available online at www.liebertpub.com/mdr).
Allelic distribution of the blaSHV gene in Shanghai
Sequencing avoids misidentification of the β-lactamase genes, and therefore is the standard molecular method to determine the alleles of the blaSHV gene. 1 By bidirectional sequencing, we identified 25 mutations, including five novel nonsynonymous mutations (V75A, A79T, D101H, S203L, and D214) (Table 3) and 10 synonymous mutations. On this basis, 32.26% of blaSHV alleles were SHV-1; 27.42% were SHV-11; and 25.81% were SHV-12 (Table 3). Many other alleles were also identified, including SHV-2A, SHV-27, SHV-28, SHV-32, SHV-33, SHV-36, SHV-71, SHV-77, SHV-93, SHV-94, SHV-95, SHV-96, and SHV-97.
SHV-type ESBLs
Sixteen isolates harbored a known SHV-type ESBL. This corresponds to a frequency of 1.21% among all the isolates, excluding the six isolates carrying novel alleles whose phenotype was not known. The prevalence of ESBL-producing K. pneumoniae and E. coli has been increasing among hospitalized patients. The isolates were K. pneumoniae (three isolates), E. cloacae (six isolates), E. coli (two isolates), M. morganii (one isolate), Acinetobacter baumannii (one isolate), P. aeruginosa (one isolate), P. fluorescens (one isolate), and Staphylococcus epidermidis (one isolate).
Coexistence of multiple blaSHV alleles within individual isolates
We investigated the coexistence of different blaSHV alleles comprehensively by directly sequencing PCR products of the full-length blaSHV gene rather than by clone sequencing. Of the 62 isolates, 12 carried two different blaSHV alleles (19.35%), including six that carried two distinct suballeles. In addition, in three isolates, SHV-12 coexisted with SHV-1, a suballele of SHV-12, and SHV-2a, respectively.
Discussion
We have developed an oligonucleotide microarray that can detect 54 mutations in the blaSHV gene. To our knowledge, this is the first time this technology has been used to detect so many mutations in blaSHV.
Even though many approaches have been developed to detect ESBLs, none is perfect. 4 Therefore, it is necessary to change cephalosporin breakpoints and apply additional ESBL detection methods. 31 Nevertheless, the optimal method to identify ESBL producers is to detect the presence of β-lactamase genes directly, not the β-lactamases they produce. In the present study, we used a positive control PCR to avoid false negatives and ensure detection of the blaSHV gene in all positive samples, not just the ESBL producers.
The coexistence of multiple genes encoding ESBLs within an individual isolate is rare.24–32,41 In contrast, the coexistence of two different alleles of the β-lactamase gene in an individual bacterium is more prevalent. It is of practical significance to distinguish different blaSHV alleles in an isolate, 16 but IEF is unable to discriminate different β-lactamases within an isolate, and current genotyping approaches can only identify a few known mutations.27,37 Because a microarray can simultaneously detect tens of thousands of SNPs, we believe that this is the best method to detect mutations in antibiotic-resistance genes. Microarray analysis also has advantages over conventional molecular approaches and antibiotic-resistance tests. We used AF limits to facilitate the detection of heterozygotes, thus overcoming some of the difficulties in assigning heterozygotes using the microarray technology.10,13
The current disadvantages of microarrays include the inability to detect novel mutations, and the limitations to multiplexing afforded by DNA preparation. 14 When these challenges are overcome, it is possible that a single microarray could detect not only all the variants of the antibiotic-resistance genes but also all the bacteria containing these genes. It is likely that many conventional clinical laboratory tests will be replaced by microarray analysis in the future.
Alleles encoding SHV-type ESBLs included SHV-2a, SHV-12, and SHV-27. Fifteen of the 16 ESBL producers carried SHV-12. Moreover, SHV-12 was found in a broad range of bacterial species, being most prevalent in E. cloacae and P. fluorescens, but also present in P. aeruginosa, M. morganii, S. epidermidis, and A. baumannii. SHV-12, which was originally reported in Switzerland, is the dominant SHV-type ESBL not only in Shanghai but also in other provinces in China15,33,40, as well as in Korea 34 and Italy. 32 SHV-5, previously reported in the Zhejiang province, China, 33 was not found in this study.
In the present study, the occurrence of the blaSHV gene was 14.52% among E. cloacae isolates, but was 33.87% among K. pneumoniae isolates. Because of this high frequency, combined with the high incidence of SHV-type ESBLs in K. pneumoniae, SHV-type ESBLs may be one of the main causes of resistance of E. cloacae to β-lactam antibiotics in Shanghai. In this study, we found that the occurrence of SHV-type ESBLs in E. cloacae (13.63%) was higher than that in K. pneumoniae (2.33%); this differs from previous reports.26,34,35 Our high incidence of SHV-type ESBLs in E. cloacae is consistent with the high prevalence of the blaSHV gene in E. cloacae reported by Hammami et al. 15 A recent survey revealed that CTX-M-type ESBLs are the most common ESBLs in E. coli and K. pneumoniae in China, but the most common type of ESBL differs among the provinces in China, implying a different use of antibiotics and unequal prevalence of ESBL genes across different geographical areas. 33
In conclusion, we developed a microarray to detect mutations in the blaSHV gene and investigated the genetic diversity of the blaSHV gene among clinical isolates from Shanghai by multiplex PCR and sequencing. The allele frequencies of the blaSHV gene in Shanghai are different from those in other regions.24,32,19 By comparison with conventional methods, our microarray can rapidly identify almost all known variants of the blaSHV gene and will be a valuable tool for epidemiological surveillance of the spread of SHV variants.
Footnotes
Acknowledgments
This work was supported by the Shanghai Science and Technology Committee Fund of Shanghai grants 06XD14219, the National High Technology Research and Development Program of China (863 Program) (Grant No. SQ2010AA1000691008), and the National Natural Science Foundation of China (Grant No. 31101091).
Disclosure Statement
The authors have no conflicts of interest to declare.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
