The evolutionary trend and genomic features of an emerging lineage of elizabethkingia anophelis strains in Taiwan

[[abstract]]The incidence of Elizabethkingia anophelis bacteremia increased significantly in a tertiary hospital, Changhua Christian Hospital (CCH) since 2013. The infection density was 1.3 and 8.1 cases per 100,000 patient-days between 2005 and 2012 and 2013 and 2020, respectively (P < 0.05). During an outbreak investigation, a specific lineage of E. anophelis strains was identified by the pulsed-field gel electrophoresis analysis. To evaluate the evolution of the specific E. anophelis lineage, whole-genome sequencing was performed, and unique genomic features (GRs) were determined by comparative genomic analysis. The specific E. anophelis lineage was novel compared to worldwide strains ever reported by cg-MLST phylogenic and whole-genome comparative analysis. Multiplex PCR using primers designed from unique GRs were performed for prevalence screening among isolates from the CCH and nationwide isolates from the Taiwan surveillance of Antimicrobial Resistance (TSAR) Program. The proportion of the specific E. anophelis lineage increased from 7.9% (3/38) during 2005-2012 to 89.2% (223/250) during 2013-2020 (P < 0.05). Although E. anophelis usually confers resistance to multiple antibiotics with limited therapeutic options, the E. anophelis strains in the specific lineage had higher ciprofloxacin resistance (100% [226/226] versus 27.4% [17/62], P < 0.05) and was associated with a higher 14-day mortality rates (33.2% [37/226] versus 16.1% [10/62], P < 0.05) than other strains at CCH. A similarly increasing trend was also found in the national TSAR program during 2002-2018 (p for trend