In Silico Identification of Three Types of Integrative and Conjugative Elements in Elizabethkingia anophelis Strains Isolated from around the World

De novo circular genome assembly of Elizabethkingia anophelis found in the mosquito Aedes aegypti from an Australian colony

We report the complete circular genome assembly of Elizabethkingia anophelis (Flavobacteriales) generated with the ONT and Illumina sequences from a laboratory-reared Aedes aegypti mosquito. This genome sequence does not belong to the lineage of known isolates from Anopheles mosquitoes, indicating that E. anophelis is genomically diverse across mosquito disease vectors.

Phylogenomics analysis of multidrug-resistant Elizabethkingia anophelis in industrial wastewater treatment plant

AIMS

This study investigated the phylogenetic relatedness of multidrug-resistant Elizabethkingia anophelis recovered from an industrial wastewater treatment plant (WWTPi).

METHODS AND RESULTS

The wastewater samples were plated in brain heart infusion agar (4 mg/L ceftazidime, 8 mg/L meropenem, and 2 mg/L polimixin). Four isolates recovered from four stages of WWTPi (influent, aeration, decantation, and treated effluent) were identified and evaluated of susceptibility profiles in the VITEK 2 system. These strains identified as E. meningoseptica were confirmed to be E. anophelis by whole genomic sequencing (Miseq-Illumina) and showed antimicrobial resistance genes of β-lactams, aminoglycosides, and tetracycline's classes. The ribosomal multilocus sequence typing showed that they belong to the rST 65620 together with clinical strains. The phylogenomic tree revealed the similarity of our strains to those belonging to sublineage 11 and the single nucleotide polymorphism analysis confirmed that they belong to a single clade.

CONCLUSIONS

To the best of our knowledge, this is the first study reporting the persistence of multidrug-resistant E. anophelis sublineage 11 along the wastewater treatment.

Population genomics of emerging Elizabethkingia anophelis pathogens reveals potential outbreak and rapid global dissemination

Elizabethkingia anophelis is an emerging species and have increasingly been reported to cause life-threatening infections and even outbreaks in humans. Nevertheless, there is little data regarding the E. anophelis geographical distribution, phylogenetic structure, and transmission across the globe, especially in Asia. We utilize whole genome sequencing (WGS) data to define a global population framework, phylogenetic structure, geographical distribution, and transmission evaluation of E. anophelis pathogens. The geographical distribution diagram revealed the emerging pathogenic bacteria already distributed in various countries worldwide, especially in the USA and China. Strikingly, phylogenetic analysis showed a part of our China original E. anophelis shared the same ancestor with the USA outbreak strain, which implies the possibility of localized outbreaks and global spread. These closer related strains also contained ICEEaI, which might insert into a disrupted DNA repair mutY gene and made the strain more liable to mutation and outbreak infection. BEAST analysis showed that the most recent common ancestor for ICEEaI E. anophelis was dated twelve years ago, and China might be the most likely recent source of this bacteria. Our study sheds light on the potential possibility of E. anophelis causing the large-scale outbreak and rapid global dissemination. Continued genomic surveillance of the dynamics of E. anophelis populations will generate further knowledge for optimizing future prevent global outbreak infections.

The Evolutionary Trend and Genomic Features of an Emerging Lineage of Elizabethkingia anophelis Strains in Taiwan

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 <0.05). We concluded that a novel lineage of E. anophelis strains has emerged dominantly in Taiwan. The genomic features are important for further investigations of epidemiology, resistance, virulence, and appropriate treatment.

IMPORTANCEElizabethkingia anophelis is an emerging multidrug resistant pathogen caused several global outbreaks recently. E. anophelis was frequently misidentified as E. meningoseptica in the past by conventional culture methods; therefore, the prevalence was often underestimated. Through revised identification, an increasing trend of E. anophelis infection was noted in a tertiary hospital and a dominant lineage of strains was recognized by genotyping. To our best knowledge, the dominant lineage of E. anophelis is novel in comparison to other worldwide strains by whole-genome comparative analysis and several unique genomic regions were found. The whole-genome sequencing data also demonstrated multiple putative virulence factors and genes associated with multidrug resistance. In our study, we identified a specially evolved E. anophelis in Taiwan with increasing nationwide dominance. This study will assist in further epidemiology surveillance and developing corresponsive infection control policies to restrain it potential of global dissemination.

Genomic analysis of Elizabethkingia species from aquatic environments: Evidence for potential clinical transmission

•

Identification of closely related (< 50 SNV) clinical and environmental aquatic Elizabethkingia anophelis isolates.

•

Identification of a provisional novel species Elizabethkingia umaracha.

•

Novel bla_GOB and bla_B carbapenemases and extended spectrum β-lactamase bla_CME alleles identified in Elizabethkingia spp.

•

Analysis of the global phylogeny and pangenome of Elizabethkingia spp.

•

Identification of novel ICE elements carrying uncharacterised genetic cargo in 67 / 94 (71.3%) of the aquatic environments Elizabethkingia spp.

Elizabethkingia species are ubiquitous in aquatic environments, colonize water systems in healthcare settings and are emerging opportunistic pathogens with reports surfacing in 25 countries across six continents. Elizabethkingia infections are challenging to treat, and case fatality rates are high. Chromosomal bla_B, bla_GOB and bla_CME genes encoding carbapenemases and cephalosporinases are unique to Elizabethkingia spp. and reports of concomitant resistance to aminoglycosides, fluoroquinolones and sulfamethoxazole-trimethoprim are known. Here, we characterized whole-genome sequences of 94 Elizabethkingia isolates carrying multiple wide-spectrum metallo-β-lactamase (bla_Band bla_GOB) and extended-spectrum serine‑β-lactamase (bla_CME) genes from Australian aquatic environments and performed comparative phylogenomic analyses against national clinical and international strains. qPCR was performed to quantify the levels of Elizabethkingia species in the source environments. Antibiotic MIC testing revealed significant resistance to carbapenems and cephalosporins but susceptibility to fluoroquinolones, tetracyclines and trimethoprim-sulfamethoxazole. Phylogenetics show that three environmental E. anophelis isolates are closely related to E. anophelis from Australian clinical isolates (∼36 SNPs), and a new species, E. umeracha sp. novel, was discovered. Genomic signatures provide insight into potentially shared origins and a capacity to transfer mobile genetic elements with both national and international isolates.

The Integrative and Conjugative Element ICECspPOL2 Contributes to the Outbreak of Multi-Antibiotic-Resistant Bacteria for Chryseobacterium Spp. and Elizabethkingia Spp

Antibiotic resistance genes (ARGs) and horizontal transfer of ARGs among bacterial species in the environment can have serious clinical implications as such transfers can lead to disease outbreaks from multidrug-resistant (MDR) bacteria. Infections due to antibiotic-resistant Chryseobacterium and Elizabethkingia in intensive care units have been increasing in recent years. In this study, the multi-antibiotic-resistant strain Chryseobacterium sp. POL2 was isolated from the wastewater of a livestock farm. Whole-genome sequencing and annotation revealed that the POL2 genome encodes dozens of ARGs. The integrative and conjugative element (ICE) ICECspPOL2, which encodes ARGs associated with four types of antibiotics, including carbapenem, was identified in the POL2 genome, and phylogenetic affiliation analysis suggested that ICECspPOL2 evolved from related ICEEas of Elizabethkingia spp. Conjugation assays verified that ICECspPOL2 can horizontally transfer to Elizabethkingia species, suggesting that ICECspPOL2 contributes to the dissemination of multiple ARGs among Chryseobacterium spp. and Elizabethkingia spp. Because Elizabethkingia spp. is associated with clinically significant infections and high mortality, there would be challenges to clinical treatment if these bacteria acquire ICECspPOL2 with its multiple ARGs, especially the carbapenem resistance gene. Therefore, the results of this study support the need for monitoring the dissemination of this type of ICE in Chryseobacterium and Elizabethkingia strains to prevent further outbreaks of MDR bacteria. IMPORTANCE Infections with multiple antibiotic-resistant Chryseobacterium and Elizabethkingia in intensive care units have been increasing in recent years. In this study, the mobile integrative and conjugative element ICECspPOL2, which was associated with the transmission of a carbapenem resistance gene, was identified in the genome of the multi-antibiotic-resistant strain Chryseobacterium sp. POL2. ICECspPOL2 is closely related to the ICEEas from Elizabethkingia species, and ICECspPOL2 can horizontally transfer to Elizabethkingia species with the tRNA-Glu-TTC gene as the insertion site. Because Elizabethkingia species are associated with clinically significant infections and high mortality, the ability of ICECspPOL2 to transfer carbapenem resistance from environmental strains of Chryseobacterium to Elizabethkingia is of clinical concern.

Structural characterization of a Type B chloramphenicol acetyltransferase from the emerging pathogen Elizabethkingia anophelis NUHP1

Elizabethkingia anophelis is an emerging multidrug resistant pathogen that has caused several global outbreaks. E. anophelis belongs to the large family of Flavobacteriaceae, which contains many bacteria that are plant, bird, fish, and human pathogens. Several antibiotic resistance genes are found within the E. anophelis genome, including a chloramphenicol acetyltransferase (CAT). CATs play important roles in antibiotic resistance and can be transferred in genetic mobile elements. They catalyse the acetylation of the antibiotic chloramphenicol, thereby reducing its effectiveness as a viable drug for therapy. Here, we determined the high-resolution crystal structure of a CAT protein from the E. anophelis NUHP1 strain that caused a Singaporean outbreak. Its structure does not resemble that of the classical Type A CATs but rather exhibits significant similarity to other previously characterized Type B (CatB) proteins from Pseudomonas aeruginosa, Vibrio cholerae and Vibrio vulnificus, which adopt a hexapeptide repeat fold. Moreover, the CAT protein from E. anophelis displayed high sequence similarity to other clinically validated chloramphenicol resistance genes, indicating it may also play a role in resistance to this antibiotic. Our work expands the very limited structural and functional coverage of proteins from Flavobacteriaceae pathogens which are becoming increasingly more problematic.

A dominant strain of Elizabethkingia anophelis emerged from a hospital water system to cause a three-year outbreak in a respiratory care center

BACKGROUND

Elizabethkingia species are ubiquitous bacteria but uncommonly cause human infection. An outbreak of Elizabethkingia anophelis bacteraemia was observed in a respiratory care center of a tertiary hospital in Taiwan from 2015 to 2018.

METHODS

Clinical and environmental isolates were collected for the outbreak investigation. Pulsed-field gel electrophoresis (PFGE) and complete-genome sequencing were conducted to elucidate the mechanism of transmission.

FINDINGS

The three-year outbreak involved 26 patients with E. anophelis bacteraemia and the incidence significantly increased during the outbreak period compared with that observed from 2010 to 2014 (P<0.05). All 26 clinical isolates during the outbreak period belonged to a cluster by PFGE analysis. In contrast, the PFGE pattern was heterogeneous among comparative historical strains. Hospital tap water was highly contaminated by Elizabethkingia species (18/34, 52.9%); among that, five E. anophelis belonged to the outbreak cluster (5/18, 27.8%). As for the inanimate surface survey, 3.4% sites (4/117) revealed positive growth of E. anophelis including two from feeding tubes/bags and two from sputum suction regulators. All four isolates belonged to the outbreak clone. The outbreak strain had no apparent relationship to currently known E. anophelis strains worldwide through complete-genome sequencing analysis. Specific infection control strategies aimed at water source control and environmental disinfection were implemented subsequently and the outbreak ended in mid-2018.

CONCLUSIONS

A specific E. anophelis strain was identified from a three-year outbreak. The elucidation of the mechanism of dominance and intra-hospital transmission is crucial for development of corresponsive infection control policies and outbreak control.

Comparative Genomics of 86 Whole-Genome Sequences in the Six Species of the Elizabethkingia Genus Reveals Intraspecific and Interspecific Divergence

Bacteria of the genus Elizabethkingia are emerging infectious agents that can cause infection in humans. The number of published whole-genome sequences of Elizabethkingia is rapidly increasing. In this study, we used comparative genomics to investigate the genomes of the six species in the Elizabethkingia genus, namely E. meningoseptica, E. anophelis, E. miricola, E. bruuniana, E. ursingii, and E. occulta. In silico DNA–DNA hybridization, whole-genome sequence-based phylogeny, pan genome analysis, and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed, and clusters of orthologous groups were evaluated. Of the 86 whole-genome sequences available in GenBank, 21 were complete genome sequences and 65 were shotgun sequences. In silico DNA–DNA hybridization clearly delineated the six Elizabethkingia species. Phylogenetic analysis confirmed that E. bruuniana, E. ursingii, and E. occulta were closer to E. miricola than to E. meningoseptica and E. anophelis. A total of 2,609 clusters of orthologous groups were identified among the six type strains of the Elizabethkingia genus. Metabolism-related clusters of orthologous groups accounted for the majority of gene families in KEGG analysis. New genes were identified that substantially increased the total repertoire of the pan genome after the addition of 86 Elizabethkingia genomes, which suggests that Elizabethkingia has shown adaptive evolution to environmental change. This study presents a comparative genomic analysis of Elizabethkingia, and the results of this study provide knowledge that facilitates a better understanding of this microorganism.

ICEKp2: description of an integrative and conjugative element in Klebsiella pneumoniae, co-occurring and interacting with ICEKp1

Klebsiella pneumoniae is a human pathogen, prominent in antimicrobial-resistant and nosocomial infection. The integrative and conjugative element ICEKp1 is present in a third of clinical isolates and more prevalent in invasive disease; it provides genetic diversity and enables the spread of virulence-associated genes. We report a second integrative conjugative element that can co-occur with ICEKp1 in K. pneumoniae. This element, ICEKp2, is similar to the Pseudomonas aeruginosa pathogenicity island PAPI. We identified ICEKp2 in K. pneumoniae sequence types ST11, ST258 and ST512, which are associated with carbapenem-resistant outbreaks in China and the US, including isolates with and without ICEKp1. ICEKp2 was competent for excision, but self-mobilisation to recipient Escherichia coli was not detected. In an isolate with both elements, ICEKp2 positively influenced the efficiency of plasmid mobilisation driven by ICEKp1. We propose a putative mechanism, in which a Mob2 ATPase of ICEKp2 may contribute to the ICEKp1 conjugation machinery. Supporting this mechanism, mob2, but not a variant with mutations in the ATPase motif, restored transfer efficiency to an ICEKp2 knockout. This is the first demonstration of the interaction between integrative and conjugative genetic elements in a single Gram-negative bacterium with implications for understanding evolution by horizontal gene transfer.

The antibiotic resistance and pathogenicity of a multidrug-resistant Elizabethkingia anophelis isolate

Elizabethkingia anophelis 12012‐2 PRCM was isolated from a patient with multiple organ dysfunction syndrome and lower respiratory tract infection in China. Minimum inhibitory concentration (MIC) analysis demonstrated that it was resistant to 20 antibiotics including trimethoprim/sulfamethoxazole and ciprofloxacin, which were effective for the elimination of other Elizabethkingia infections. To investigate multidrug resistance and pathogenicity mechanisms, we analyzed genome features of 12012‐2 PRCM and compared them to the other Elizabethkingia species. The draft genome size was 4.02 Mb with a GC content of 32%, comparable to that of other E. anophelis strains. Phylogenetic analysis showed that E. anophelis 12012‐2 PRCM formed a sister group with E. anophelis 502, distinct from clades formed by other clinical and environmental E. anophelis isolates. E. anophelis 12012‐2 PRCM contained multiple copies of β‐lactamase genes as well as genes predicted to function in antimicrobial efflux. It also contained 92 genes that were potentially involved in virulence, disease, and defense, and were associated with resistance and pathogenicity. Comparative genomic analysis showed high homology among three clinical and two environmental E. anophelis strains having a variety of similar antibiotic resistance and virulence factor genes, and similar genomic structure. Applications of this analysis will contribute to understanding the antibiotic resistance and pathogenic mechanisms of E. anophelis infections, which will assist in the management of infections as it increases in prevalence.