Abstract
Objective:
The objective of this study was to examine the species distribution, genetic relatedness, virulence gene profiles, antimicrobial sensitivities, and resistance gene distribution of clinical Aeromonas strains from Singapore and Malaysia.
Methods:
A total of 210 Aeromonas clinical isolates were investigated: 116 from Singapore General Hospital and 94 archived clinical isolates from University of Malaya Medical Center, Malaysia. The isolates were genetically identified based on the gcat gene screening and the partial sequences of the rpoD housekeeping gene. Genetic relatedness, distribution of 15 virulence genes and 4 beta-lactamase resistance genes, and susceptibility patterns to 11 antimicrobial agents were compared.
Results:
Of the 210 Aeromonas isolates, A. dhakensis–94 (45%) was the dominant species in Singapore and Malaysia. Species composition was similar and enterobacterial repetitive intergenic consensus-PCR did not show genetic relatedness between strains from the two countries. Of the 15 virulence genes, A. dhakensis and A. hydrophila harbored the most compared with other species. Different combinations of 9 virulence genes (exu, fla, lip, eno, alt, dam, hlyA, aexU, and ascV) were present in A. dhakensis, A. hydrophila, and A. veronii from both the countries. Distribution of virulence genes was species and anatomic site related. Majority (>80%) of the strains were susceptible to all antimicrobial agents tested, except amoxicillin and cephalothin. A. dhakensis strains from Malaysia significantly harbored the cphA gene compared with A. dhakensis from Singapore. Multidrug resistance was mostly detected in strains from peritoneal fluids of dialysis patients.
Conclusion:
This study revealed A. dhakensis as the dominant species isolated in both geographic regions, and that it carried a high number of virulence genes. It also highlights the geographic-related differences of virulence gene distribution and antimicrobial resistance profiles of clinical Aeromonas strains from Singapore and Malaysia.
Introduction
A
The taxonomy of Aeromonas is in revision and expansion; the genus currently consists of 36 species. Aeromonads are not difficult to isolate, but identification to species level is challenging due to its phenotypic heterogeneity. 2 Molecular approaches are more useful and definitive in the identification of Aeromonas. 3 The gcat gene was found to be present in practically all Aeromonas strains, and a genus-specific probe of high specificity was designed for Aeromonas detection.4,5 Compared with the use of 16s rRNA gene, housekeeping gene rpoD provided a more definitive identification of the genus.3,6,7
Aeromonas dhakensis (formerly known as A. aquariorum) 8 was found to be the prevalent species among Aeromonas isolates from Malaysia as identified by molecular methods, and these strains displayed various resistance phenotypes and genotypes.7,9 Aeromonas is known to carry numerous virulence determinants such as cell-associated factors, extracellular factors, siderophores, enterotoxins, a variety of proteases, motility factors, and type 3 secretion system (T3SS) effectors. 1
Thus far, there are no reports on the species diversity of Aeromonas from Malaysia's neighboring country, Singapore. Hence, the aim of this study was to (i) molecularly identify clinical isolates of Aeromonas from Singapore, (ii) evaluate the genetic relationship of Aeromonas isolates, (iii) investigate the possible virulence potential of Aeromonas, (iv) study the antimicrobial susceptibility patterns, (v) detect presence of resistance genes, and (vi) perform comparative analysis between isolates from Singapore and Malaysia.
Materials and Methods
Bacterial isolates
Two hundred ten clinical isolates of Aeromonas were investigated: 116 from the Singapore General Hospital and 94 previously characterized clinical isolates from University Malaya Medical Center. 7 These were from stool, peritoneal dialysates, blood, pus, wounds, bile, tissue and body fluids, sputum, bone, lung, stents, and urine. All isolates were cultured and cryopreserved in 20% (v/v) glycerol at −80°C. Working cultures were maintained in Luria-Bertani agar and broth.
Molecular identification and phylogenetic analysis
Isolates were genetically identified using a combination of gcat gene screening and rpoD gene sequencing as described previously. 7 Phylogenetic analysis was performed as described in our previous report 10 using a partial nucleotide sequence (511 bp) of the rpoD gene (GenBank accession numbers: KY198403–KY198518 and JN686647–JN686741) and reference gene sequences as listed in Table 1.
Enterobacterial repetitive intergenic consensus-PCR and fingerprint analyses
Enterobacterial repetitive intergenic consensus-PCR (ERIC-PCR) which amplifies the conserved repetitive DNA regions present throughout bacterial genomeswas performed using the primers and PCR protocol as described previously. 11 Amplification products were electrophoresed on a 1.5% (w/v) agarose gel at 56 V for 5 hours in TBE buffer. Fingerprint profiles were analyzed using software BioNumerics version 7.5 (Applied Maths, Belgium). Similarities between the fingerprints were calculated with Dice coefficient and a dendrogram was constructed using the unweighted pair group method with average linkages.
Virulence gene screening
Detection of 15 virulence genes encoding for DNAse (exu), serine protease (ser), aerolysin/hemolysin (aer), cytotoxic enterotoxin (act), heat-labile cytotonic enterotoxin (alt), heat-stable cytotonic enterotoxin (ast), lipase (lip), flagellin (fla), elastase (ela), ADP-ribosyltransferase toxins (aexT and aexU), DNA adenine methyltransferase (dam), enolase (eno), T3SS membrane component (ascV), and hemolysin (hlyA) was carried out. Specific primers were selected and PCR was performed to screen the 15 virulence genes using previously published protocols. 10 Two-tailed Fisher's exact test was used to determine the presence of combination of virulence genes.
Antimicrobial susceptibility
Susceptibilities to 11 antimicrobial agents of 7 different classes were determined by the disc diffusion method as described by CLSI guideline M45-A212 using Mueller-Hinton agar (Oxoid) and commercial antibiotic discs (Oxoid). The antimicrobial agents were: amoxicillin (10 μg), cefepime (30 μg), cephalothin (30 μg), cefotaxime (30 μg), chloramphenicol (30 μg), tetracycline (30 μg), amikacin (30 μg), gentamicin (10 μg), ciprofloxacin (5 μg), norfloxacin (10 μg), and meropenem (10 μg). Bacterial suspensions of 0.5 McFarland (1.5 × 108 colony-forming unit/ml) were prepared in 0.85% saline solution and inoculated onto Mueller-Hinton agar plates using sterile rayon swabs. Antibiotic discs were dispensed onto the inoculated plates using a dispenser or sterile forceps. The plates were incubated at 35°C for 18 hours and inhibition zones were measured and interpreted based on the CLSI guideline M45-A2. 12 Escherichia coli ATCC 25922 and Pseudomonas aeruginosa ATCC 27853 were used as quality controls.
Detection of resistance genes
All strains were investigated by PCR for genes related to beta-lactam resistance: MOX, VIM, KPC, and cphA using primers as previously published13,14 (Table 2).
Results
Detection of gcat gene and rpoD gene sequencing
Based on gcat gene screening all 210 isolates belonged to the genus Aeromonas and phylogenetic analysis of the partial sequences of the rpoD gene showed distinct clustering of 165 isolates with known Aeromonas species with bootstrap values of 85–99%. Forty-five strains formed a cluster with type strains of A. veronii, A. allosaccharophila, A. finlandiensis and A. aquatilis at a bootstrapping value of 70% and were considered as “A. veronii group” in this study. Based on the phylogenetic tree, 116 strains from Singapore were identified as A. dhakensis–47, A. veronii group–33, A. caviae–22, A. hydrophila–10, A. jandaei–2, A. taiwanensis–1, and A. sanarellii–1 (Supplementary Fig. S1; Supplementary Data are available online at www.liebertpub.com/mdr). The strains from Malaysia were: A. dhakensis–47, A. caviae–18, A. hydrophila–16, A. veronii group–12, and A. trota–1.
Aeromonas dhakensis was the dominant species in both countries and recovered mostly from stool samples (40/94, 43%) (Table 3). Based on pairwise sequence identity matrices, the intraspecies similarities for the strains were 79.6–100% for A. dhakensis, 94.8–100% for A. hydrophila, 87.0–100% for A. caviae, and 78.2–100% for A. veronii group and 94.1–99.1% for A. jandaei.
A. dhakensis from lung, fluid, and stent; A. hydrophila from bone; A. veronii from stent.
S, strains from Singapore; M, strains from Malaysia.
ERIC-PCR and fingerprint analyses
Among the 210 Aeromonas strains, two clusters of 100% similarity level were identified: two isolates of A. veronii group (SC109 and SC110) from Singapore and two isolates of A. hydrophila (C71 and C72) from Malaysia. The isolates within each cluster shared the identical fingerprint and was considered as the same clone (Supplementary Fig. S2). The remaining strains exhibited distinct ERIC fingerprinting patterns and no genetic similarity was observed between strains from Singapore and Malaysia.
Virulence genes
In both regions, A. hydrophila and A. dhakensis harbored more virulence genes than A. veronii group or A. caviae. The average number of virulence genes carried by A. hydrophila, A. dhakensis, A. veronii group, and A. caviae were 11, 10, 6, and 5. One A. dhakensis strain from a stool specimen of a Malaysian patient carried the full complement of all the 15 virulence genes. Two-tailed Fisher's exact test revealed statistically significant (p < 0.05) association of A. dhakensis with alt, dam, and hlyA; A. hydrophila with alt, ast, dam, and hlyA and A. veronii group with act among strains from both regions. Carriage of multiple enterotoxin genes (alt, ast and act) was seen only in A. dhakensis and A. hydrophila with the latter having a higher proportion than A. dhakensis (Table 4). These enterotoxigenic strains were mostly from pus or wound specimens (Table 5).
Significant difference (p < 0.05) between isolates from Singapore and Malaysia.
Very significant difference (p < 0.0001) between isolates from Singapore and Malaysia.
Significant association (p < 0.05) of virulence gene subset with source of isolation.
Very significant association (p < 0.0001) of virulence gene subset with source of isolation.
The frequency of the eno gene was significantly higher among A. dhakensis (p < 0.0001), A. hydrophila (p < 0.05) and A. veronii group (p < 0.05) from Malaysia than those from Singapore. The fla gene was significantly associated with A. dhakensis (p < 0.05) and A. veronii group (p < 0.0001) strains from Singapore. T3SS genes ascV and aexU were almost exclusively (p < 0.0001) detected in A. dhakensis from Malaysia either singly or in combination. Every strain positive for ascV and/or aexU gene carried the aexT gene as well. Strains with a subset of virulence genes aexT/ascV/aexU were the most prevalent (6/14, 43%) from pus associated with chronic osteomyelitis, burn wounds, fractures, and traumatic wounds. Strains from pus/wound specimens harbored the highest number of virulence genes compared with specimens from other anatomic sites (Table 5).
Antimicrobial susceptibility
More than 80% of the strains were susceptible to cefepime, cefotaxime, chloramphenicol, amikacin, gentamicin, ciprofloxacin, norfloxacin, and meropenem. Susceptibility to cephalothin was significantly higher in A. veronii group compared with A. dhakensis, A. caviae, and A. hydrophila (p < 0.0001) for both geographical regions.
Ninety-seven percent of the strains were resistant to amoxicillin. In both countries, A. caviae was significantly more resistant against chloramphenicol than A. dhakensis (p < 0.05), A. veronii group (p < 0.05), and A. jandaei (p < 0.05). Resistance against tetracycline, ciprofloxacin, and norfloxacin was also significantly higher (p < 0.05) in A. caviae than A. dhakensis and A. veronii group.
Resistance to cefotaxime (p < 0.0001), ciprofloxacin (p < 0.05), and norfloxacin (p < 0.05) was more significant in the Singapore strains than those from Malaysia, and vice versa for cephalothin (p < 0.05) and gentamicin (p < 0.0001). At the species level, Malaysian A. dhakensis strains were significantly resistant to gentamicin (p < 0.05), whereas A. dhakensis from Singapore was significantly resistant to cefotaxime (p < 0.05) and tetracycline (p < 0.05) (Table 6).
Significant difference (p < 0.05) between isolates from Singapore and Malaysia.
Very significant difference (p < 0.0001) between isolates from Singapore and Malaysia.
S, strains from Singapore; M, strains from Malaysia.
Multidrug resistance (MDR) is defined as nonsusceptibility to at least one agent in three or more antimicrobial categories: nonsusceptibility refers to resistant, intermediate, and nonsusceptible in in vitro antimicrobial susceptibility testing. 15 MDR patterns were observed in 26% (54/210) of the strains and most prevalent from peritoneal fluids (Table 7). One A. caviae from peritoneal fluid from Singapore was resistant to all classes of antimicrobial agents, except the aminoglycosides. Eleven MDR strains were from peritoneal fluids of renal failure patients, but the MDR trait was not specifically associated with any species or geographic region.
Very significant (p < 0.0001) association with source of isolation.
Detection of resistance genes
The MOX gene was detected in 56 (27%) of the strains and significantly associated with A. caviae regardless of country of origin. The cphA gene was detected in 7 (7%) A. dhakensis, 3 (11%) in A. hydrophila, and 3 (7%) in A. veronii group, but none of the A. caviae strains had the cphA gene (Table 8). All cphA-positive strains were susceptible to meropenem, but the cphA gene was not detected in the two meropenem-resistant strains. There was no significant difference in the distribution of resistance genes between strains of the same species from Singapore and Malaysia, except that the frequency of the MOX gene was significantly higher in A. dhakensis strains from Malaysia than those from Singapore. None of the strains carried the VIM and KPC genes.
Very significant difference (p < 0.0001) between isolates from Singapore and Malaysia.
Discussion
The rpoD gene was used by Puthucheary et al. as a molecular tool for inferring the taxonomy of Aeromonas, 7 but since then 15 new species have been described: A. aquatica, A. aquatilis, A. australiensis, A. cavernicola, A. crassostreae, A. diversa, A. enterica, A. finlandiensis, A. fluvialis, A. intestinalis, A. lacus, A. rivipollensis, A. rivuli, A. sanarellii, and A. taiwanensis.16–24 Using the updated taxonomy we reevaluated the phylogeny of the Malaysian isolates used by Puthucheary et al., 7 and found the identity of each strain to be the same as described previously, except that A. veronii isolates were clustered into A. veronii group consisting of A. veronii, A. allosaccharophila, A. finlandiensis, and A. aquatilis. These four species are closely related, forming a distinct cluster with a bootstrap value of 100% following multilocus phylogenetic analysis and in silico DNA–DNA hybridization. 23 Isolates within the A. veronii group would require additional molecular markers for more refined discrimination.
Aeromonas dhakensis was identified as the predominant species for both countries. Reports of A. dhakensis clinical isolates are scanty probably due to previous identification as A. hydrophila.8,25,26 A. dhakensis is increasingly recognized as a pathogen of considerable virulence with reports of infection in Taiwan and Australia.27,28 Our results confirm a previous report that A. dhakensis and A. hydrophila possess more virulence genes than A. veronii group or A. caviae, 29 but situational expression of virulence factors due to genetic regulation is reported to occur in Aeromonas. 30 Chen et al. found A. dhakensis more pathogenic than A. hydrophila because of higher biofilm formation, higher minimum inhibitory concentration for certain antimicrobial agents, lower survival rate in Caenorhabditis elegans, higher cytotoxicity against human normal skin fibroblast cell lines and mouse fibroblast cell lines, and shorter lifespan of infected mice. 31 Furthermore, A. dhakensis bacteremia contributed to higher mortality rates in patients compared with other species. 28 A. dhakensis is thus probably more virulent than the other species in vitro and in vivo. 25
The Aeromonas species from clinical material were similar in both Singapore and Malaysia, perhaps due to comparable environmental factors, such as climate, temperature and humidity, and dietary patterns with an annual consumption of 20–50 kg of fish per capita, but cooking methods may vary allowing the survival of bacteria. Presence of Aeromonas in different environments has been reported in both countries.10,32,33 Moreover, seafood supply in Singapore relies largely on imports from Malaysia. 34 There are reports of A. hydrophila isolation from aquaculture in Malaysia, and other Asian nations have reported A. dhakensis contamination in marine shrimps cultured in low-salinity ponds in Thailand and farmed eels in Korea35–37 Nevertheless, ERIC fingerprinting demonstrated that strains recovered from patients in Singapore and Malaysia were different clones, suggesting that patients were not infected by clonally related strains. Our data show genetic diversity of the strains and that no predominant pathogenic clones were circulating among patients in both the countries. Similarly, Szczuka and Kaznowski did not find evidence of predominant clones responsible for Aeromonas-associated infections with no genetic relationship between stool isolates from different parts of Europe and Asian countries. 11
Uncommon Aeromonas species detected in this study was A. trota, an ampicillin-susceptible species from a stool specimen of a Malaysian patient with severe gastroenteritis and watery diarrhea, and one strain each of A. taiwanensis and A. sanarellii from peritoneal fluid and foot tissue culture from Singapore.
Different subsets of virulence genes were detected in different Aeromonas species and only one A. dhakensis strain (strain C51) from a stool specimen from Malaysia carried the full complement of all 15 virulence genes, suggesting its high virulence potential. T3SS genes (ascV and aexU) were only present in A. dhakensis from Malaysia either singly or in combination and they were not present in the Singapore strains. These results indicate that different subsets of virulence genes of the T3SS complex are present in subpopulations of Aeromonas, suggesting distinct pathogenic mechanisms.
Aeromonas harboring both alt and ast enterotoxic genes causing watery diarrhea may represent true diarrheal pathogens in Southeast Asia. 38 Ninety-six percent of our A. hydrophila strains carried both the enterotoxin genes, but no significant association with diarrheal cases was recorded (Table 2). Only 12% (10/85) of stool isolates carried both the enterotoxin genes that were reported to be exclusively detected in diarrheal isolates. 38 Our findings imply that enterotoxigenic aeromonads in Singapore and Malaysia were not unequivocally considered diarrheal pathogens as they were frequently isolated from peritoneal fluids and blood cultures. Strains from pus/wound specimens carried the widest array of virulence genes probably indicating their ability to invade the mucosa and deeper tissues.
Cephalothin sensitivity was exclusive to A. veronii group and this feature confirms previous reports that it be used as a phenotypic marker to differentiate A. veronii from other Aeromonas species.39,40 High resistance rates to beta-lactam antibiotics, amoxicillin (97%), was seen. The frequency of the MOX gene–an AmpC beta-lactamase was higher among A. dhakensis strains from Malaysia. Geographic-related variations in antimicrobial susceptibility probably depend on several factors, for example, differences in healthcare systems, guideline recommendations, and compliance of patients. Antimicrobial prescribing policies in Malaysian government hospitals are guided by the national database, 41 whereas healthcare systems in Singapore refer to hospital policies or published guidelines, although antibiotics are not usually prescribed for self-limiting diarrhea.
Resistance to meropenem by the disc diffusion method, but absence of the cphA gene was seen in 1% of our strains, and it may have been due to the presence of other resistance genes, efflux pumps, and porin mutations. 42
Our A. caviae strains were resistant to many antimicrobial agents as 68% of them carried the MOX gene belonging to class C beta lactamase, as AmpC present in Aeromonas can hydrolyze many beta-lactamases. 43 An A. caviae isolate with resistance to all antimicrobial agents was recovered from tropical ornamental fish and carriage water from Singapore, and it is possible that this strain may be disseminated during export of these fish to other countries. 44
MDR was recorded in 13 (68%) of peritoneal strains with 9 of the 13 from renal failure specimens. Colonization and infection rates by MDR strains were reported to be higher in patients on dialysis than patients in other healthcare settings. 45 One A. caviae strain from a postmortem peritoneal specimen from Singapore was susceptible to only the aminoglycosides. The high rate of MDR in peritoneal dialysis patients found in our study and published reports 45 highlights the necessity for increased awareness of adherence to infection control measures in hospitals, and among healthcare workers and caregivers.
Conclusions
The prevalence of A. dhakensis infections in Singapore and Malaysia is an indication of the virulence potential rather than antimicrobial resistance of this species. ERIC fingerprinting suggests that infections in Singapore and Malaysia were not clonally related. The diverse distribution of the virulence genes denotes that different subsets of virulence determinants occur in different subpopulations of Aeromonas species. Definitive species identification using molecular methods is necessary to predict virulence potential and resistant phenotypes. High prevalence of multidrug-resistant strains from peritoneal dialysis fluids highlights the necessity for increased awareness and adherence to good infection control measures.
Footnotes
Acknowledgments
The authors gratefully appreciate and thank members of the participant hospitals for their donation of clinical isolates. This work was funded by the University of Malaya Research Grant RP039A-15HTM and Postgraduate Research Grant PG082-2014B from the University of Malaya.
Disclosure Statement
No competing financial interests exist.
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.
