Phylogeographical Landscape of Citrobacter portucalensis Carrying Clinically Relevant Resistomes

ABSTRACT During a surveillance study conducted to assess the occurrence and genomic landscape of critical priority pathogens circulating at the human-animal-environment interface in Brazil, as part of the Grand Challenges Explorations-New Approaches to Characterize the Global Burden of Antimicrobial Resistance program, two multidrug-resistant (MDR) Citrobacter portucalensis carrying blaCTX-M-15 extended-spectrum β-lactamase (ESBL) genes, isolated from green sea turtles, were characterized. Genomic and phylogeographical analysis of C. portucalensis genomes available in public databases revealed the intercontinental dissemination of clades carrying different arrays of clinically relevant genes conferring resistance to carbapenems, broad-spectrum cephalosporins, cephamycins, aminoglycosides and fluoroquinolones, disinfectants, and heavy metals. Our observations suggest that C. portucalensis could be emerging as critical priority bacteria of both public and One Health importance worldwide. IMPORTANCE The global spread of antibiotic-resistant priority pathogens beyond the hospital setting is a critical issue within a One Health context that integrates the human-animal-environment interfaces. On the other hand, next-generation sequencing technologies along with user-friendly and high-quality bioinformatics tools have improved the identification of bacterial species, and bacterial resistance surveillance. The novel Citrobacter portucalensis species was proposed in 2017 after taxonomic reclassification and definition of the strain A60T isolated in 2008. Here, we presented genomic data showing the occurrence of multidrug-resistant C. portucalensis isolates carrying blaCTX-M-15 ESBL genes in South America. Additionally, we observed the intercontinental dissemination of clades harboring a broad resistome to clinically relevant antibiotics. Therefore, these findings highlight that C. portucalensis is a global MDR bacteria that carries intrinsic blaCMY- and qnrB-type genes and has become a critical priority pathogen due to the acquisition of clinically relevant resistance determinants, such as ESBL and carbapenemase-encoding genes.

IMPORTANCE The global spread of antibiotic-resistant priority pathogens beyond the hospital setting is a critical issue within a One Health context that integrates the human-animal-environment interfaces. On the other hand, next-generation sequencing technologies along with user-friendly and high-quality bioinformatics tools have improved the identification of bacterial species, and bacterial resistance surveillance. The novel Citrobacter portucalensis species was proposed in 2017 after taxonomic reclassification and definition of the strain A60 T isolated in 2008. Here, we presented genomic data showing the occurrence of multidrug-resistant C. portucalensis isolates carrying bla CTX-M-15 ESBL genes in South America. Additionally, we observed the intercontinental dissemination of clades harboring a broad resistome to clinically relevant antibiotics. Therefore, these findings highlight that C. portucalensis is a global MDR bacteria that carries intrinsic bla CMY -and qnrB-type genes and has become a critical priority pathogen due to the acquisition of clinically relevant resistance determinants, such as ESBL and carbapenemaseencoding genes.
M embers of the genus Citrobacter are part of the normal intestinal flora of humans and animals and have been isolated from a variety of environmental sources, including soil, water, and food; being considered opportunistic pathogens for humans (1,2). The clinical significance of Citrobacter is due to its ability to easily acquire multiple resistance genes to antibiotics used in therapy. In this regard, some species of the genus carry intrinsic AmpC-type b-lactamase genes (i.e., bla CMY -type), which can be related to mobile genetic elements (1,3). Additionally, ESBLs and carbapenemases have been identified in Citrobacter spp., expanding the resistance to broad-spectrum cephalosporins and carbapenems. Currently, this genus comprises 15 species (https:// lpsn.dsmz.de/genus/citrobacter) with the most recently described Citrobacter species being C. portucalensis, C. europaeus, and C. cronae (1,4,5). Although the novel C. portucalensis species was proposed in 2017, the strain A60 T was isolated in 2008 from a water well sample collected in Cantanhede city, Centre region of Portugal (1,2). Consequently, by using genome data of C. portucalensis A60 T , the average nucleotide identity (ANI) results have led to change originally genomes submitted as C. freundii and Citrobacter spp., to C. portucalensis since 2018.
Antibiotic-resistant priority pathogens have globally disseminated beyond the hospital environment, threatening wildlife and ecosystems within a One Health context that integrates the human-animal-environment interface (6)(7)(8). During a nationwide surveillance study (One Health Brazilian Resistance Program, OneBR), conducted in Brazil to assess the occurrence and genomic features of World Health Organization (WHO) priority pathogens circulating at the human-animal-environment interface (6), two ceftriaxone-resistant Gram-negative bacilli isolated from cloacal swab samples collected from 2 green sea turtles affected by fibropapillomatosis (a debilitating neoplastic disease with immunosuppressive effects), were sequenced. Both green sea turtles were admitted to a rescue and rehabilitation center in Southeastern Brazil (23°27'08.0"S, 45°04'13.6"W). Isolates were recovered from MacConkey agar plates supplemented with ceftriaxone (2 mg/mL) and incubated at 37°C overnight. Bacterial identification and antimicrobial susceptibility testing were performed by Vitek 2 system (bioMérieux, USA), disc diffusion, Etest, and/or agar dilution methods (9). Colistin MIC was determined by broth microdilution method, according to EUCAST (http://www.eucast.org). Extended-spectrum b-lactamase (ESBL) production was screened by the double-disk synergy test (DDST) (10).
Results from WGS analysis are summarized in Table S1. The resistome of both C. portucalensis strains was composed of genes conferring resistance to b-lactams (bla CTX-M-15 , bla CMY-129 -like, bla TEM-1B , bla OXA-1 ), aminoglycosides, quinolones, sulfonamides, and trimethoprim. C. portucalensis TV06 harbored additional resistance genes to phenicols, fosfomycin, and tetracyclines. On the other hand, genes conferring resistance to heavy metals (arsenic, mercury, silver, and tellurium), disinfectants (quaternary ammonium compounds [QACs]), and pesticides (glyphosate) were detected in both genomes (Table S1). Plasmid analysis revealed the presence of IncHI2 in C. portucalensis TV06, and IncHI1, IncFII, and IncFIB plasmids in the strain TV13. In this regard, in silico analysis mapping the contigs against paired-end short reads along to BLASTn analysis revealed that bla CTX-M-15 genes were carried by IncHI and IncFII plasmids in C. portucalensis TV06 and TV13, respectively (Table S1). Unfortunately, because we used short-read sequencing technology, it was not possible to assemble the complete plasmids.
Although, MLST analysis revealed that the TV13 strain belonged to ST63, it was not possible to predict the ST of the TV06 strain. SNP-based phylogenetic analysis grouped CTX-M-15-positive TV06 and TV13 strains in two distinct clades, both related to C. portucalensis strains isolated from human and nonhuman sources in distinct geographic regions (Fig. 1). The matrix of SNP-based phylogeny analysis is quoted in Table S2. Specifically, the TV06 strain clustered (14549 to 15041 SNP differences) with C. portucalensis strains MBTC-1222 (isolated from food in Nigeria; assembly accession number: GCA_002843195.2), 1001302B_160321_C4 (isolated from a human in the USA; assembly accession number: GCA_015668205.1), and EC_47 and EC_49 (isolated from environmental sources in Nigeria; Assembly accession numbers: GCA_019584605.1 and GCA_019584545.1) (Fig. 1,  Table S2). On the other hand, TV13 was clustered (21254 to 21829 SNP differences) with the human strains AM17-37 and sc18407397, identified in China and Germany, respectively (assembly accession numbers: GCA_003471655.   (Table S1). HM, heavy metal resistome. DSF, disinfectant resistome. NI, not informed.
After having your updated name formally recognized, C. portucalensis has been further isolated from uziza leaves in Nigeria (15), nasal and urine samples in dogs from Japan (3), urinary tract infection in a cat from Austria (16), poultry dropping samples from Bangladesh (17), soil from USA (18), and human stool and sputum samples in China (19,20). Moreover, from genome data deposited in the NCBI database, vegetables, humans, animals, and the environment have been registered as sources of isolation of C. portucalensis in Bangladesh, China, Colombia, Denmark, Germany, Japan, Nepal, Nigeria, Pakistan, Portugal, Rwanda, South Korea, Spain, Thailand, United Kingdom and USA (Table S1).
In this study, we report the emergence of CTX-M-15-producing C. portucalensis in marine wildlife, in South America. In this regard, green turtles (C. mydas) are endangered marine animals that have a global distribution, occurring throughout tropical and subtropical waters (https://www.iucnredlist.org/species/4615/11037468). Due to their highly migratory behavior, green turtles have been investigated as bioindicators of polluted ecosystems contaminated with MDR bacteria (21). Therefore, the occurrence of this sort of bacteria in natural environments can have serious ecological implications (22), adding another layer of difficulty to marine wildlife conservation strategies. On the other hand, considering that C. portucalensis is closely related to C. freundii (an important nosocomial pathogen) (1), continuous surveillance is required to provide more information about the global dissemination of this emergent species and its adaptation to different hosts.
In summary, our findings highlight that C. portucalensis is a global MDR pathogen with intrinsic genes encoding AmpC enzymes active on cephamycin (CMY), which has become a critical priority specie due to the acquisition of clinically significant resistance determinants, such as plasmid-mediated bla CTX-M -type ESBL and carbapenemaseencoding genes (i.e., bla GES , bla VIM , bla NDM , bla IMP , and bla KPC ). Finally, because marine animals have been overlooked in the epidemiology of critically important pathogens, more research focusing on the transmission pathways of ESBL-producing bacteria in marine environments is required to understand the clinical and epidemiological impacts on their populations.

SUPPLEMENTAL MATERIAL
Supplemental material is available online only. We thank Cefar Diagnóstica Ltda. (São Paulo, Brazil) and CEFAP-GENIAL facility for kindly supplying antibiotic discs for susceptibility testing and Illumina sequencing, respectively.
We declare no conflict of interest.