Aeromonas dhakensis is an emergent human pathogen with medical importance. This study was aimed to determine the sequence types (STs), genetic diversity, and phylogenetic relationships of different clinical sources of 47 A. dhakensis from Malaysia using multilocus sequence typing (MLST), goeBURST, and phylogenetic analyses. The analysis of a concatenated six-gene tree with a nucleotide length of 2994 bp based on six housekeeping genes (gyrB, groL, gltA, metG, ppsA, and recA) and independent analyses of single gene fragments was performed. MLST was able to group 47 A. dhakensis from our collection into 36 STs in which 34 STs are novel STs. The most abundant ST521 consisted of five strains from peritoneal fluid and two strains from stools. Comparison of 62 global A. dhakensis was carried out via goeBURST; 94.4% (34/36) of the identified STs are novel and unique in Malaysia. Two STs (111 and 541) were grouped into clonal complexes among our strains and 32 STs occurred as singletons. Single-gene phylogenetic trees showed varying topologies; groL and rpoD grouped all A. dhakensis into a tight-cluster with bootstrap values of 100% and 99%, respectively. A poor phylogenetic resolution encountered in single-gene analyses was buffered by the multilocus phylogenetic tree that offered high discriminatory power (bootstrap value = 100%) in resolving all A. dhakensis from A. hydrophila and delineating the relationship among other taxa. Genetic diversity analysis showed groL as the most conserved gene and ppsA as the most variable gene. This study revealed novel STs and high genetic diversity among clinical A. dhakensis from Malaysia.
Keywords: Aeromonas dhakensis; Clonal complex; Genetic diversity; Multilocus phylogenetic analysis; Sequence types.