Comparative Genomic Analysis of 19 Clinical Isolates of Tigecycline-Resistant Acinetobacter baumannii

Abstract

To assess the genomic profiles of tigecycline (Tgc)-resistant Acinetobacter baumannii, including antibiotic resistance (AR) genes and virulence factors (VF), whole-genome shotgun sequencing was performed on 19 Tgc-resistant (TgcR) A. baumannii strains collected in a tertiary hospital during the early phase of the clinical introduction of Tgc in China from late 2012 to mid-2014. The major sample types containing TgcR strains were sputum and drain fluid. Data from an average of 624 Mbp of sequence was generated on each bacterial genome, with Q30 quality of 90%, and an average coverage of 96.6%. TCDC-AB0715 was used as a reference genome. The genome sequences were annotated for functional elements including AR genes, VFs, genome islands, and inserted sequences before they were comparatively analyzed. The antibiotic susceptibility phenotypes of the strains were examined by a broth microdilution method to determine the minimal inhibitory concentration (MIC) of strains against major clinical antibiotics. The AR genes (ARGs) were annotated using the Comprehensive Antibiotic Resistance Database (CARD). Thirty-three ARGs were shared by all 19 TgcR strains, and 24 ARGs were distributed differently among strains. A total of 391 VFs were found to be diversely distributed in all TgcR strains. Based on ARG number distribution, the 19 TgcR strains were divided into several groups. Highly differentiated genes included gpi, mphG, armA, msrE, adec, catB8, aadA, sul1, blaOXA–435, aph3i, and blaTEM–1, which may represent gene markers for TgcR A. baumannii sub-types. In addition, when compared with Tgc-sensitive (TgcS) strains collected during the same period, TgcR strains featured enrichment of ARGs including aph6id, aph3ib, and teta. Compared with 26 other whole-genome sequences of A. baumannii deposited in GeneBank, TgcR strains in this study commonly lacked the EF-Tu mutation for elfamycin resistance. Previous investigation of three A. baumannii strains isolated from one patient indicated genomic exchange and a homologous recombination event associated with generation of tigecycline resistance. This study further analyzed additional TgcR strains. Phylogenetic analysis revealed a close evolutionary relationship between 19 TgcR strains and to isolates in East and Northeast China. In short, the comprehensive functional and comparative genomic analysis of 19 clinical TgcR A. baumannii strains isolated in the early stage of Tgc usage in China revealed their close phylogenetic relationship yet variable genetic background involving multiple resistance mechanisms. Using a simple ARG or VF gene number diversity method and marker genes, TgcR strain sub-types can be identified. The distinct characteristics of TgcR A. baumannii strains with versatile genomic resistance and regulation patterns raise concern regarding prediction and control of Tgc resistance in the clinic.