Anna Kopf
Boyke Bunk
Sina M. Coldewey
Florian Gunzer6
Thomas Riedel
Percy Schröttner- 1Medical Microbiology and Virology, University Hospital Carl Gustav Carus, Dresden, Germany
- 2Clinic for Hematology and Oncology, Carl-Thiem-Klinikum, Cottbus, Germany
- 3German Collection of Microorganisms and Cell Cultures GmbH, Leibniz Institute DSMZ, Braunschweig, Germany
- 4Department of Anesthesiology and Intensive Care Medicine, Jena University Hospital, Jena, Germany
- 5Septomics Research Center, Jena University Hospital, Jena, Germany
- 6Department of Hospital Infection Control, University Hospital Carl Gustav Carus, Dresden, Germany
- 7German Center for Infection Research (DZIF), Partner Site Hannover-Braunschweig, Braunschweig, Germany
Recent studies suggest that Wohlfahrtiimonas chitiniclastica may be the cause of several diseases in humans including sepsis and bacteremia making the bacterium as a previously underappreciated human pathogen. However, very little is known about the pathogenicity and genetic potential of W. chitiniclastica; therefore, it is necessary to conduct systematic studies to gain a deeper understanding of its virulence characteristics and treatment options. In this study, the entire genetic repertoire of all publicly available W. chitiniclastica genomes was examined including in silico characterization of bacteriophage content, antibiotic resistome, and putative virulence profile. The pan-genome of W. chitiniclastica comprises 3819 genes with 1622 core genes (43%) indicating a putative metabolic conserved species. Furthermore, in silico analysis indicated presumed resistome expansion as defined by the presence of genome-encoded transposons and bacteriophages. While macrolide resistance genes macA and macB are located within the core genome, additional antimicrobial resistance genotypes for tetracycline (tetH, tetB, and tetD), aminoglycosides (ant(2’’)-Ia, aac(6’)-Ia,aph(3’’)-Ib, aph(3’)-Ia, and aph(6)-Id)), sulfonamide (sul2), streptomycin (strA), chloramphenicol (cat3), and beta-lactamase (blaVEB) are distributed among the accessory genome. Notably, our data indicate that the type strain DSM 18708T does not encode any additional clinically relevant antibiotic resistance genes, whereas drug resistance is increasing within the W. chitiniclastica clade. This trend should be monitored with caution. To the best of our knowledge, this is the first comprehensive genome analysis of this species, providing new insights into the genome of this opportunistic human pathogen.
Introduction
Wohlfahrtiimonas chitiniclastica was first isolated from the larvae of Wohlfahrtia magnifica (Tóth et al., 2008), an obligate parasitic fly that causes myiasis by depositing eggs and larvae in mammalian wounds in both animals and humans (Robbins and Khachemoune, 2010). W. chitiniclastica are Gram-negative, strictly aerobic, non-motile rods. A key feature is strong chitinase activity, which may be an indicator of a symbiotic relationship with its host fly while also playing an important role in metamorphosis (Tóth et al., 2008; Schröttner et al., 2017; Snyder et al., 2020). Apart from the close association between W. chitiniclastica and certain flies (Tóth et al., 2008; Cao et al., 2013; Maleki-Ravasan et al., 2020), the bacteria have recently been detected in various habitats around the world such as arsenic-affected soils from Bangladesh (Sanyal et al., 2016), chicken meat purchased in Brazilian supermarkets (Matos et al., 2019), the pancreas of a Zebra in China (Zhou et al., 2016), and human soft tissue infection in Estonia (Kõljalg et al., 2015), to name but a few. Finally, yet importantly, recent studies indicate that W. chitiniclastica can be the cause of several diseases in animals (Thaiwong et al., 2014; Diaz-Delgado et al., 2015; Qi et al., 2016) and humans, including sepsis and bacteremia (Almuzara et al., 2011; Kõljalg et al., 2015; Campisi et al., 2015; Suryalatha et al., 2015), making the bacterium a previously underestimated human pathogen (Schröttner et al., 2017). Although the pathogenesis of W. chitiniclastica has not been fully elucidated, the bacterium is thought to invade traumatic skin lesions via fly larvae, resulting in severe myiasis and/or wound contamination (Robbins and Khachemoune, 2010; Thaiwong et al., 2014; Schröttner et al., 2017). However, since most clinicians are unfamiliar with this species and conventional methods often lead to misidentification (Kõljalg et al., 2015; Kopf et al., 2021), it can be assumed that W. chitiniclastica has been poorly recognized as a possible cause and is even more common than originally thought (Kopf et al., 2021).
To date, NCBI lists 26 genomes of W. chitiniclastica strains, 22 of which have been isolated in the course of human disease. In addition, three draft genomes of strains isolated from an animal source have been published (Cao et al., 2013; Zhou et al., 2016; Matos et al., 2019), and annotations revealed genes encoding for macrolide-specific efflux pumps (macA and macB) (Matos et al., 2019), a blaVEB-1 gene cassette, which confers resistance to ceftazidime, ampicillin, and tetracycline (Zhou et al., 2016), and a genome-encoded 25.9kb intact phage (Matos et al., 2019). However, apart from these preliminary genomic studies, very little is known about the genetic potential of W. chitiniclastica, making it necessary to initiate systematic studies in order to gain more insight into its virulence characteristics as well as treatment options. Therefore, in this study, the entire genetic repertoire of all publicly available W. chitiniclastica genomes was investigated; including in silico characterization of the antibiotic resistome, prophage content, and virulence potential. In addition, we performed a pan-genome analysis to elucidate the major genome features and genetic variability of W. chitiniclastica. To the best of our knowledge, this is the first comprehensive genome analysis of W. chitiniclastica and allows us to better understand this previously underestimated human pathogen.
Materials and Methods
Genomic Strain Collection
All publicly available W. chitiniclastica genomes were included in this study (n = 26). These include a total of 14 W. chitiniclastica isolates from Dresden (Germany) that have been collected in routine diagnostics over a period of six years (Schröttner et al., 2017; Kopf et al., 2021). These isolates were recovered exclusively from diagnostic cultures analyzed at the Institute for Medical Microbiology and Virology, University Hospital Carl Gustav Carus (Dresden, Germany). Whole-genome sequences were submitted to NCBI GenBank under Acc. Nos JAGIBR000000000-JAGICE000000000, applying the NCBI Prokaryotic Annotation Pipeline PGAP (Tatusova et al., 2016) as previously reported (Kopf et al., 2021). In addition, the remaining publicly available W. chitiniclastica genomes at NCBI as of April 2021 (n = 12) were included in this study. The corresponding datasets were retrieved in preassembled nucleotide FASTA files and GenBank files. These include the type strain DSM 18708T (AQXD01000000) (Tóth et al., 2008), BM-Y (LVXD00000000) (Zhou et al., 2016), Strain 20 (LWST01000000) (Matos et al., 2019), SH04 (AOBV01000000) (Cao et al., 2013), and 8 genomes, which were submitted by the Centers for Disease Control and Prevention of the United States but are not associated to a citable publication as far as we know. The latter comprise of ATCC 51249 (NEFL01000000), F6512 (NEFK01000000), F6513 (NEFJ01000000), F6514 (NEFI01000000), F6515 (NEFH01000000), F6516 (NEFC01000000), F9188 (NEFE01000000), and G9145 (NEFC01000000). Isolation source, host, year of isolation, and geographical origins of isolates were taken from published research papers, otherwise estimated using dates/locations on public databases as indicated. Results are displayed in Table 1; Supplementary Table S1.
Table 1 Overview of the general genome features of the W. chitiniclastica genomes analyzed in this study.
In Silico Genome Analysis
Functional genome analysis was performed using the freely available computational tools with default parameters from January till June 2021, unless indicated otherwise. Preassembled FASTA files were annotated using Prokka (Galaxy version 1.14.6+galaxy0) (Seemann, 2014) and strarmar (Galaxy Version 0.7.2+galaxy0) (Petkau, 2018). Results are displayed in Figure 3; Supplementary Tables S6, S7. Search for antimicrobial resistance profiling was extended by using the comprehensive antibiotic resistance database CARD (https://card.mcmaster.ca/) (Alcock et al., 2020) retaining “Perfect hit and Strict hit only” and “High-quality/coverage”. Results are displayed in Figure 3; Supplementary Table S8. Phage analysis was performed using PHASTER (PHAge Search Tool Enhanced Release) (https://phaster.ca/) (Arndt et al., 2016). Results are displayed in Table 1; Supplementary Table S2. Analysis of CRISPR (clustered regularly interspaced short palindromic repeats) and their associated (Cas) proteins was done using the CRISPRCasFinder (Couvin et al., 2018) (https://crisprcas.i2bc.paris-saclay.fr/). Only results with evidence levels 3 and 4 were included in the analysis. Results are displayed in Table 1; Supplementary Table S1. AcrFinder (http://bcb.unl.edu/AcrFinder/index.php) (Yi et al., 2020) was used for the detection of Anti-CRISPR (Acr) proteins. Results are displayed in Table 1; Supplementary Table S3.
Phylogenetic Identification
For phylogenomic identification; genomic contigs were submitted to the Type Strain Genome Server at https://tygs.dsmz.de/ (Meier-Kolthoff and Göker, 2019). Gene sequence of the 16S rRNA gene and rpoB gene were retrieved from the previous results using Prokka (Galaxy version 1.14.6+galaxy0), and BLAST analysis for the homology of the 16S rRNA gene and the rpoB gene for the identification of W. chitiniclastica was performed on https://blast.ncbi.nlm.nih.gov (Altschul et al., 1990).
Pan-Genome Assembly and Visualization
Preassembled GenBank files were converted to GFF3 using the ‘Genbank to GFF3’ converter (Galaxy Version 1.1) (Stajich et al., 2002). Annotated GFF3 files of 26 W. chitiniclastica genomes were submitted to Roary (Galaxy Version 3.13.0+galaxy1) (Page et al., 2015) for pan genome analysis using default parameters. A gene-absence-presence data matrix was derived and visualized with Phandango (Hadfield et al., 2018). Results are displayed in Figure 1; Supplementary Tables S4, S5.
Figure 1 Linearized visualization of W. chitiniclastica pan-genome along with the phylogenetic tree. The pangenome was visualized based on the software Phandango (Hadfield et al., 2018).
Results
Characterization of Central Genome Features
Over a period of six years, a total of 14 W. chitiniclastica strains were recovered from wound swabs in routine diagnostics (Table 1) (Kopf et al., 2021). The genome size of these W. chitiniclastica isolates, linked in terms of geographic location and isolation source, ranged from 2.07 to 2.14 million bases with an average of 1941 predicted genes (Table 1). To broaden the picture, we extended our study to include all W. chitiniclastica genomes, which were publicly available at NCBI as of April 2021. Thus, we arrived at 26 genomes comprising 22 isolates from human sources and 4 strains from animal sources (Table 1). Genome size ranged from 1.99 to 2.18 million bases with an average of 1967 predicted genes, 51 tRNAs, and five rRNAs. Notably, all isolates from Dresden and Strain 20 have three rRNAs, whereas the other genomes contain an average of seven rRNAs. Interestingly, BM-Y (n = 4) and Strain 20 (n = 2) contain multiple copies of the 16S rRNA gene (Supplementary Table S1).
The CRISPRCasFinder (Couvin et al., 2018) was used for the identification of potential genes encoding for CRISPR (clustered regularly interspaced short palindromic repeats) arrays and their associated (Cas) proteins. Because CRISPR arrays of evidence levels 1 and 2 are potentially invalid (Couvin et al., 2018), we focused on the results of evidence levels 3 and 4, which are considered as highly likely candidates. All strains contain CRISPR repeats, spacers, and the cas cluster CAS-Type IF (Table 1; Supplementary Table S1). On average, the isolates contain two CRISPR and 85 spacers. DSM 108048 stands out with four CRISPR and 130 spacers, while DSM 100375 contains only one potential CRISPR sequence and eight spacers. The highest number of spacers was detected in Strain 20 isolated from an animal source. Interestingly, all W. chitiniclastica isolates also contain genomically encoded anti-CRISPR (acr) genes (Table 1). DSM 108048, SH04, ATCC 51249 and G9145 stand out with 10 genomically encoded Acr proteins, whereas the type strain DSM 18708, Strain 20, and DSM 105712 contain only two sequences. Although some acr genes were labeled ‘low confidence’ (Supplementary Table S3), the majority are homologs of known Acr proteins, making the actual presence of potential anti-CRISPR proteins highly likely.
Phylogenetic Identification
Genome-based-taxonomy analysis of all strains revealed correct assignment to W. chitiniclastica. Hereby, digital DNA:DNA hybridization (dDDH) values of 74.0-75.2% were computed against the type strain DSM 18708T and therefore fulfilling the criteria for correct bacterial species identification (Meier-Kolthoff and Göker, 2019) (Supplementary Table S1). Construction of a phylogenomic tree based on whole-genome sequences revealed that all strains cluster in one subclade with the type strain DSM 18708T, and 25 strains form a subspecies using a 79% dDDH threshold (Meier-Kolthoff et al., 2014) (Figure 2). Surprisingly, the type strain does not belong to the subspecies. Furthermore, 16S rRNA gene and rpoB gene sequences were compared for sequence similarity to other sequences using the BLAST database. All isolates were identified with a minimum of 98% out of 100% sequence identity to corresponding 16S rRNA and rpoB reference genes, respectively, leading to the correct assignment of W. chitiniclastica (Supplementary Table S1). Notably, the additional 16S rRNA gene copies of BM-Y (n = 4) and Strain 20 (n = 2) revealed correct taxonomic assignment and share high homology, ranging from 99.7%-99.9%.
Figure 2 Phylogenomic tree of the W. chitiniclastica species and subspecies delineation based on the GBDP phylogenetic analyses retrieved and modified based on the Type (Strain) Genome Server (TYGS). The tree was inferred with FastME 2.1.6.1 (Lefort et al., 2015) from GBDP distances calculated from genome sequences and was subjected to a clustering using established thresholds for delineating species (DDH > 70%) (Meier-Kolthoff et al., 2013) as well as subspecies (DDH > 79%) (Meier-Kolthoff et al., 2014). The branch lengths are scaled in terms of GBDP distance formula d5 and are represented in red numbers. The numbers above branches are GBDP pseudo-bootstrap support values >60% from 100 replications, with an average branch support of 83.0%.
Pan-Genome Construction
We constructed the pan-genome of W. chitiniclastica encompassing 26 genomes, which represents the first and largest analysis of this type to date. Roary (Page et al., 2015) was used to cluster the genes encoding complete protein sequences into core (hard core and soft core) and accessory (shell and cloud) genomes. The core genome is shared by every genome tested in this study and its genes are often related to housekeeping functions. It can further be divided into hard core genes, which are defined to be present in >99% of the genomes, and soft core genes, which are present in 95–99%. The accessory genome is shared by a subset of the genomes and is associated with, but not limited to, pathogenicity or environmental adaptation. It is subdivided into shell genes, which are present in 15–95%, and cloud genes, which are found in less than 15% of genomes. The latter include singletons or genes found in only one of the genomes.
The pan-genome of all 26 W. chitiniclastica isolates comprises 3819 genes; 1622 core genes (43%), and 2197 accessory genes (57%) (Figure 1) with 1117 unique genes (29%) defined as genes only present in one strain. We also observed a remarkable abundance of 1240 (32%) genes with an unknown function. The core genome divides into 1175 hard (31%) and 447 soft core (12%) genes. Notably, 92% of the core genes (1494/1622) encode a known function, while only 8% (128/1622) are assigned to hypothetical proteins. The accessory genome comprises 635 shell (17%) and 1562 cloud (41%) genes. Remarkably, 51% (1122/2197) of the accessory genes code for hypothetical proteins, a fact that highlights the limited characterization of the W. chitiniclastica genome.
The majority of core genes encodes for protein families associated with housekeeping functions such as amino acid metabolism, energy production and translation, to name but a few (Supplementary Table S4). Genes associated with defense mechanisms such as antimicrobial resistance (AMR) genes are mainly present in the accessory genome (Supplementary Table S5). Notably, major multidrug efflux systems, on the other hand, are encoded within the core genome (Supplementary Table S4). These include the adenosine triphosphate (ATP)-binding cassette (ABC) superfamily, resistance nodulation-division (RND) family, major facilitator superfamily (MFS), small multidrug resistance (SMR) family, multidrug and toxic compound extrusion (MATE) family, proteobacterial antimicrobial compound efflux (PACE) family, and p-aminobenzoyl-glutamate transporter (AbgT) family. In addition, we detected several protein families encoding TRAP transporters and TonB dependent transport systems.
The generated phylogenetic tree of the pan-genome shows three main clades (Figure 1). DSM 108048 and ATCC 51249 are clustered in a single lineage (clade 1), and clade 2 is composed of six isolates from Dresden. The remaining isolates are grouped in clade 3, suggesting potential spread and transmission between hosts without any clear geographical links or host association.
Prediction of Arsenic Resistance
Pan-genome analysis revealed two protein families encoding for arsenic resistance proteins within the core genome (Supplementary Table S4). Interestingly, genes belonging to the ubiquitous ars operon (arsRDABC) such as arsC, arsD, arsA, and arsB (Carlin et al., 1995) were only detected in DSM 100375, DSM 105712, DSM 105838, F6512, F6513, and F6514 (Figure 3). In addition, these six isolates contain the gene for the inorganic arsenic efflux pump acr3 (Fekih et al., 2018).
Figure 3 Heat-map visualizing genes derived from mobile genetic elements (red), and genes associated with putative virulence traits (blue), arsenic resistance (green) and antimicrobial resistance (yellow). Cell colors indicate the presence of genes: grey (absence); red, blue, green and yellow (presence).
Prediction of Bacteriophages and Other Mobile Genetic Elements
Prediction of prophage sequences within the W. chitiniclastica genomes revealed a total of 75 prophages (Table 1), 18 of which are classified as intact (Supplementary Table S2). The latter are most likely to be complete and functional, and were found in 16 different genomes all isolated from a human source. The most common intact phage was identified either as Mannheimia (n = 12) or Enterobacteria phage (n = 4) (Supplementary Table S2). The remaining prophages were classified as “questionable” (n = 4) and “incomplete” (n = 53). Notably, none of the strains isolated from an animal source contained intact phages.
We then extended our search for the presence of other mobile genetic elements as it can improve the fitness and enables bacteria to acquire new AMR genes (Kottara et al., 2018; Stalder et al., 2019). Each genome contained genes homologous to Tn7 transposon genes such as tnsA, tnsB, tnsC, and tnsE (Peters, 2014; Babakhani and Oloomi, 2018) (Figure 3). This is consistent with pan-genome analysis, which revealed a ubiquitous presence of a transposase encoding protein family within the core genome (Supplementary Table S4). Interestingly, the accessory genome showed numerous additional transposases as well as integrases and phage proteins (Supplementary Table S5). Detailed genome analysis based on Prokka confirmed these results (Supplementary Table S6). For example, F9188 and G9145 harbor Tn10 encoded tetracycline resistance genes tetR (Smith and Bertrand, 1988) and tetC (Babakhani and Oloomi, 2018) (Figure 3). Prophage integrase intA gene, which is necessary for the integration of the phage into the host genome by site-specific recombination (Kirby et al., 1994), was present in all genomes except DSM 100375. Integrase gene intS (Panis et al., 2010), on the other hand, was only detected in seven genomes including DSM 100375. Notably, BM-Y contains five intS copies (Supplementary Table S6). Additional genome-integrated phage genes such as yfdh, yfdg (Rudd, 1999) and gpFI (Popovic et al., 2012) were distributed among the W. chitiniclastica genomes (Figure 3).
In Silico Annotation of Potential Virulence Traits
We analyzed all genomes including the type strain DSM 18708T with respect to potential virulence traits (Supplementary Table S6). To begin with, all isolates harbor the conserved virulence factor B (cvfB) (Matsumoto et al., 2007) and genes encoding for active multidrug efflux systems. The latter are one of the major mechanisms of bacterial resistance to drugs. For example, we detected the multidrug resistance transporter gene mdtA (Nagakubo et al., 2002), the MATE transporter genes norM (Nishino et al., 2006) and mepA (Kaatz et al., 2006), the TonB dependent transport genes tdhA (Thomas et al., 1998), exbD, and exbB (Bosch et al., 2002; Holden et al., 2012), and multidrug efflux pumps genes acrB and arcA (Eicher et al., 2009), to name but a few. Secretion systems are another key element in the pathogenesis of bacterial infections, and all genomes harbor copies of the type II secretion (TS2) genes xcpT and epsF (Sandkvist, 2001).
Other meaningful elements involved in pathogenicity are toxins. Although W. chitiniclastica has not been reported to be a prolific toxin-producing organism, all genomes were manually searched for potential candidates. Based on our in silico results, 14 isolates contain the toxin relG gene (Korch et al., 2015). Notably, F9188, G9145 and SH04 harbor two copies of relG.
Toxin-antitoxin (TA) modules are ubiquitous among bacteria, and are involved in diverse physiological processes including virulence (Lobato-Márquez et al., 2016). Interestingly, the core genome harbors a protein family coding for the type II toxin-antitoxin system RatA (Supplementary Table S4), which is described as antisense RNA that blocks the accumulation of the mRNA for the TxpA toxin (Silvaggi et al., 2005). However, the corresponding txpA gene appears to be missing. Complete TA modules, on the contrary, can be found in the accessory genome (Supplementary Table S5). These include the RelE/ParE family, which encompasses several smaller toxin families including YoeB (Anantharaman and Aravind, 2003). In-depth genome analysis using Prokka (Seemann, 2014) confirmed these results by revealing genes encoding for the TA type II system YefM-YoeB (Norton and Mulvey, 2012) in 14 W. chitiniclastica strains (Figure 3). Notably, SH04 harbors two copies of yefM and yoeB. The PasTI is another known TA module (Norton and Mulvey, 2012), and F6512 contains both genes, pasT and pasI, respectively. The remaining 25 isolates only harbor the pasT gene.
In Silico Analysis of Antimicrobial Resistance Genes
In order to perform an thorough and comprehensive search for AMR genes, we extended our previous analysis using the comprehensive antibiotic resistance database CARD (Alcock et al., 2020) retaining ‘Perfect’ and ‘Strict’ hits only while excluding ‘Loose’ hits (Supplementary Table S8). Noteworthy, the ‘Strict’ algorithm represents a flexible sequence variation but lies within the curated BLAST bit score cut-off (Alcock et al., 2020), and by that, making a correct identification highly feasible. The ‘Loose’ algorithm, on the other hand, works outside of the detection model cut-off to provide detection of new and more distant homologs of AMR genes. Although it could help to identify potential resistance genes and/or shed light on new unknown modifications, in silico results based on loose hits should always be taken with caution and require further research. In order to gain high-quality results we decided to restrict the CARD search to ‘Perfect’ and ‘Strict’ hits, and combine the outcome with results gained by Prokka (Seemann, 2014) and strarmar (Petkau, 2018). Summarized results are displayed in Figure 3, and detailed information can be retrieved from Supplementary Tables S6, S7.
Genes coding for macrolide-specific efflux pumps (macA and macB) (Kobayashi et al., 2000; Yum et al., 2009), and tellurite resistance methyltransferase (tehB) (Choudhury et al., 2011) were detected in all 26 genomes. Tetracycline resistant efflux protein tetH gene (Roberts, 2005) was found in 7 genomes; notably, BM-Y harbors double tetH genes. In addition, DSM 100676, DSM 100917, DSM 105708, and BM-Y contained the tetracycline repressor protein tetD gene (Roberts, 2005). G9145 and F9188, on the other hand, harbor the tetracycline efflux protein tetB gene, which confers resistance to tetracycline, doxycycline, and minocycline, but not tigecycline (Roberts, 2005). BM-Y contains two aminoglycoside resistance genes encoding for the adenylyltransferase ant(2’’)-Ia (Cox et al., 2015) and acetyltransferase aac(6’)-Ia (Parent and Roy, 1992). The remaining aminoglycoside resistance genes such as aph(3’’)-Ib, aph(3’)-Ia (Tauch et al., 2000), and aph-ld (Cao et al., 2013)-Id (Chiou and Jones, 1995) were detected in F6513, F6515, G9188, and F9188, respectively. Sulfonamide resistance gene sul2 and streptomyicn resistance gene strA, which have been described to be encoded within the same resistance gene cassette (Anantham and Hall, 2012), were detected in F6513, F6514, G9145, and F9188. In addition, BM-Y harbored a gene encoding for the VEB-1 beta-lactamase (Nordmann and Naas, 1994), and G9145 contains the chloramphenicol acetyltransferase cat3 gene (Vassort-Bruneau et al., 1996). Finally, yet importantly, DSM 108045 contains the antibiotic efflux pump encoding abaF gene, which has been reported, when expressed, to increase fosfomycin resistance (Sharma et al., 2017).
Discussion
In the past 12 years, several case reports have shown that W. chitiniclastica is capable of causing sepsis and bacteremia in humans labeling this organism as a newly underestimated pathogen. However, little is known about its pathogenicity and genome content. Our current study, analyzing all publicly available W. chitiniclastica genomes to date, highlights significant genomic characteristics including potential virulence factors, and AMR genes. Moreover, we provide the first pan genome analysis and shed light on the core features of W. chitiniclastica.
Distinct Genomic Characteristics of W. Chitiniclastica and Their Effect on the Assessment of Microbial Diversity
Genome size ranged from 1.99 to 2.18 million bases with an average of 1967 predicted genes, 51 tRNAs, and 5 rRNAs. Interestingly, all isolates from Dresden and Strain 20 harbor 3 rRNAs, while the remaining genomes contained an average of 7 rRNAs. It has recently been proposed that multiple rRNA operons confer a selective advantage to respond quickly and grow rapidly in environments characterized by fluctuations in resource availability (Stevenson and Schmidt, 2004). Based on this hypothesis, with 12 and 9 rRNAs, respectively, SH04 and BM-Y should have a fitness advantage when compared to the other isolates. In addition, BM-Y and Strain 20 harbor multiple 16S rRNA gene copies; a fact that has not been reported for any W. chitiniclastica strain yet. Although most bacterial genomes exhibit only one or two 16S rRNA genes (Armougom, 2009), some microbes contain multiple or varying numbers (Rainey et al., 1996; Acinas et al., 2004). For example, Bacillus subtilis has 10 copies (Stewart et al., 1982) and Clostridium paradoxum has up to 15 copies with heterogeneous intervening sequences (Rainey et al., 1996). Those multiple copies are often associated with nucleotide sequence variability (Rainey et al., 1996; Acinas et al., 2004), and/or provide insufficient taxonomic resolution power at the species level or with closely related species (Janda and Abbott, 2007). Our study indicates that the phylogenetic identification of those additional 16S rRNA gene copies of BM-Y and Strain 20 provides correct classification; however it can lead to an overestimation in terms of abundance and diversity composition (Armougom, 2009), a fact that one has to keep in mind within the scope of a thorough microbial community profiling. 16S rRNA-based identification generally faces a number of challenges, including sequencing error (Schirmer et al., 2015), primer bias (Klindworth et al., 2012), and varying discrimination power between variable regions (Sune et al., 2020) and certain bacterial genera (Roux et al., 2011). In light of these drawbacks, an alternative, or at least complementary taxonomic markers, should be considered. Our study indicates that both dDDH and rpoB based analysis have proven a worthy identification method for W. chitiniclastica. While dDDH is a very costly and time-consuming technique, the rpoB gene has emerged as a new marker gene candidate for phylogenetic analyses and identification of bacteria (Adékambi et al., 2009). Although the rpoB gene does not have a database currently as comprehensive as that of the 16S rRNA gene (Bondoso et al., 2013), the approach has a number of advantages. For example, rpoB excels with an increased phylogenetic resolution on the genus level or lower, and it is a single-copy protein-encoding gene enabling a phylogenetic analysis on amino acid and nucleotide level (Case et al., 2007; Adékambi et al., 2009). However, since neither the dDDH nor the rpoB approach is widely established in clinical routine diagnostics yet, 16S rRNA-based identification most likely remains by far the most frequently used method. Nevertheless, in case of doubtful results, additional methods should be considered as an alternative or complement.
Mobile Genetic Elements and Their Putative Effect in Shaping the Genetic Diversity of W. Chitiniclastica
Genome scanning for phages and MGE revealed the ubiquitous presence of transposons and bacteriophages within all W. chitiniclastica isolates tested in this study, and we believe that they are some of the key elements in shaping the genetic diversity within the W. chitiniclastica clade. Interestingly, all genomes harbor CRISPR-Cas elements, which are described to constitute the adaptive immune system in prokaryotes in order to provide resistance against invasive genetic elements including viruses, plasmids, and transposons (Barrangou et al., 2007). Therefore, in theory, W. chitiniclastica should be well equipped against its invasion. However, the ubiquitous distribution of MGEs among prokaryotes suggests that CRISPR systems are not always functional (Grissa et al., 2007), participate in other processes, such as signal transduction and gene regulation (Westra et al., 2014; Faure et al., 2019), and/or have a yet undiscovered function (Aydin et al., 2017). Moreover, many phages have evolved an anti-CRISPR system that inhibits the CRISPR immune response of their host (Landsberger et al., 2018). Although Acr proteins were first discovered in Pseudomonas phages and other prophages (Bondy-Denomy et al., 2013), they have also been detected in other prokaryotes such as Moraxella bovoculi (Marino et al., 2018). In fact, >30% of P. aeruginosa strains carrying a CRISPR-Cas system also encode one or more cognate acr genes (van Belkum et al., 2015); therefore, the ubiquitous presence of genome-encoded Acr proteins among W. chitiniclastica may partly explain the widespread distribution of MGE and phage-related genes. Although the conditions and extent to which these immunosuppressive genes allow bacteriophages to persist in their bacterial host remain unclear (Landsberger et al., 2018), the balance between CRISPR-Cas immunity and Acr activities may be a central element in shaping the genetic diversity of W. chitiniclastica including antimicrobial resistome expansion.
Pan-Genome Composition of W. Chitiniclastica
To begin with, the pan-genomic phylogenetic tree indicated no clear host or geographical clustering suggesting a potential spread and transmission. Although six strains isolated from Dresden clustered within a subclade, the analysis might be limited in terms of capturing total diversity due to the fact that the majority are from the same location, and only four are associated with an animal source. With hopefully increasing numbers of available genomes from various locations, we recommend repeating the phylogenetic pan-genome analysis to see whether a certain niche specificity emerges.
The composition of the pan-genome revealed a core genome of 43%, which appears to be conserved when compared to other reported species such as Clostridium perfringens (12.6%) (Kiu et al., 2017), Aliarcobacter butzleri (22%) (Buzzanca et al., 2021), Staphylococcus aureus (32%), Pseudomonas aeruginosa (26%) (Costa et al., 2020), Klebsiella pneumoniae (26%) and Salmonella enterica (16%) (McInerney et al., 2017), to name but a few. Microbes with large and diverse accessory genomes, on the other hand, are considered metabolically versatile species with the ability to migrate to new niches, and to adapt to changing environmental conditions (Rouli et al., 2015; Costa et al., 2020). Accessory genes are often acquired by Horizontal gene transfer (HGT) (Costa et al., 2020), and are related to virulence, antimicrobial defense or confer a fitness advantage (Jordan et al., 2001). Therefore, species with a large and conserved core genome often lack a diverse pool of virulence and AMR factors. This is congruent with recent studies regarding Bordetella pertussis. Members of this species have a large core genome (59%) suggesting that due to this low genomic diversity, antibiotics and vaccines are quite effective against this species (Tettelin et al., 2008; Carbonetti, 2016; Costa et al., 2020). This is consistent with recent observations regarding W. chitiniclastica isolates, which are described to be susceptible to the majority of known antibiotics with the exception of fosfomycin (Schröttner et al., 2017; Matos et al., 2019), and by that, underlines the assumption that members of this species appear to be metabolically conserved when compared to others. However, as mentioned above, with increasing numbers of available genomes a large-scale pan-genome analysis is recommended to confirm and/or reevaluate the findings.
Further Evidence of a Previously Newly Described Subspecies of W. Chitiniclastica
Originally isolated from Wohlfahrtia magnifica larvae in Hungary (Tóth et al., 2008), there is increasing evidence that W. chitiniclastica inhabits diverse niches and environmental habitats such as soil (Sanyal et al., 2016), humans (Schröttner et al., 2017), other mammals (Thaiwong et al., 2014; Diaz-Delgado et al., 2015), fish (Reddy and Mastan, 2013), and diverse flies (Lysaght et al., 2018; Maleki-Ravasan et al., 2020). In most of these published studies, identification was based on 16S rRNA gene sequence, which is known to lack sufficient resolution to distinguish between closely related species (Roux et al., 2011). The dDDH analysis performed during this work surprised with further evidence of a previously newly described subspecies of W. chitiniclastica (Kopf et al., 2021) (Figure 2). Originally thought to be the adaptation to a human environment and geographic location (Kopf et al., 2021), this work rather suggest a broad host and environmental range of W. chitiniclastica. This observation is also reflected in the pangenomic phylogenetic tree (Figure 1). The fact that W. chitiniclastica can colonize different species should be considered an advantage for the bacterium. However, on the other hand, it also poses an increased risk for zoonotic transmission, whose dynamic interactions between humans, animals, and pathogens should be considered in the context of the “One Health” approach (Welch et al., 2007; Cantas and Suer, 2014).
Prediction of Arsenic Resistance
W. chitiniclastica has recently been detected in arsenic-affected soils from Bangladesh (Sanyal et al., 2016) indicating its ability to inhabit different habitats including soil, humans, and animals (Tóth et al., 2008; Thaiwong et al., 2014; Diaz-Delgado et al., 2015; Sanyal et al., 2016; Schröttner et al., 2017; Matos et al., 2019). Arsenic occurs naturally in aquatic and terrestrial environments, and despite its relatively low abundance, the high toxicity of arsenic derivatives is considered a severe problem of public health worldwide (Fekih et al., 2018). In this process, microorganisms are known to play a crucial role in global arsenic geocycles (Zhu et al., 2014) subsequently leading to the ubiquitous presence of arsenic resistance genes among microbes (Fekih et al., 2018). With the presence of arsenic-resistant protein families within the core genome, W. chitiniclastica is no exception. However, genes of the common arsRDABC operon were only detected in six genomes. Notably, the regulator protein-encoding arsR gene, which acts as a repressor of the arsRDABC operon in the absence of arsenic (Cai et al., 2009), appears to be missing. These findings indicate the development of a yet unknown regulation and/or arsenic tolerance mechanism distributed among W. chitiniclastica. Moreover, these 6 isolates harbor the arsenic efflux pump acr3 gene (Fekih et al., 2018), which is surprising as most prokaryotic species are described to have either an arsB or acr3 gene (Yang et al., 2012). Notably, arsenic-resistant bacterial isolates from highly arsenic-contaminated soils showed a predominance of acr3 genes over arsB genes (Cai et al., 2009) suggesting that DSM 100375, DSM 105712, DSM 105838, F6512, F6513, and F6514 have a severe fitness advantage in highly enriched arsenic habitats.
In Silico Profiling of Potential Virulence Traits
To detect possible genomic signatures linked to virulence, all genomes were manually searched for genes putatively associated with host-pathogen interaction. The ubiquitous presence of diverse multidrug efflux systems emphasizes a central role in the pathogenesis of W. chitiniclastica. For example, the core genome harbors a PACE efflux transporter, which is described to confer resistance to a wide range of biocides used as disinfectants and antiseptics, and are encoded by many Gram-negative human pathogens (Hassan et al., 2018). Other potential virulence components are represented by membrane-associated proteins, like TonB dependent transport systems, which are involved the in virulence of Shigella dysenteriae, Haemophilus influenzae and E. coli, to name but a few (Reeves et al., 2000; Torres et al., 2001; Morton et al., 2012). Moreover, the ubiquitous presence of the conserved virulence factor B (cvfB) suggests a central role in the virulence of W. chitiniclastica. Recent studies showed that deletion of CvfB results in reduced virulence in S. aureus and decreased production of hemolysin, DNase, and protease (Matsumoto et al., 2007), which further emphasized its importance for pathogenicity. However, in-depth research including target specific manipulations is required to unravel its function in W. chitiniclastica.
TA modules are involved in diverse physiological processes providing bacteria with pronounced fitness advantages dependent on toxin expression levels and the specific environmental niche occupied (Ma et al., 2021). These include bacterial adaptation to hostile environments, mediating stress response, stabilization of chromosomal regions, and bacterial survival during infection (Lobato-Márquez et al., 2016). Our in silico analysis indicates the presence of the TA system YefM-YoeB (Norton and Mulvey, 2012) in 14 W. chitiniclastica strains. This TA module has been described to be involved in the niche-specific colonization, stress resistance, and survival inside the host (Norton and Mulvey, 2012). Although the natural habitat of W. chitiniclastica is not well investigated, we could envision that TA systems might be involved in invading different habitats including persistence as part of a polymicrobial infection. Moreover, F6512 harbors a second TA system comprised of the toxin PasT and the antitoxin PasI. The PasTI module assures cell formation in the presence of antibiotics and increases pathogen resistance to nutrient limitation as well as oxidative and nitrosative stresses (Norton and Mulvey, 2012). Notably, the remaining 25 isolates only harbor the pasT gene, whose function has recently been reannotated based on new experimental evidence. While it was shown that PasT sustains antibiotic tolerance, and is critical for the formation or survival of ciprofloxacin-tolerant cells, the function of PasTI as a TA system could not be confirmed (Fino et al., 2020). Instead, the supposed toxin PasT is a bacterial homolog of mitochondrial protein Coq10 suggesting a central role in respiratory electron transport by acting as an important accessory factor in the ubiquinone-dependent electron transport chain (Fino et al., 2020). This leaves us to speculate, whether the pasT gene of W. chitiniclastica is primarily involved in virulence and/or energy production.
The secretion of proteins and toxins have a major role in the pathogenesis of bacterial infections, and several highly specialized pathways have evolved for this purpose such as the T2S system. The latter has been widely discovered in a number of bacterial species including several human pathogens like Chlamydia trachomatis, Escherichia coli, Klebsiella pneumoniae, Legionella pneumophila, Vibrio cholerae, and P. aerunginosa, to name but a few (Sandkvist, 2001; Durand et al., 2003; Cianciotto and White, 2017). In general, proteins secreted by T2S systems are associated with the destruction of various tissues, cell damage and diseases such as proteases, cellulases, pectinases, phospholipases, lipases, and toxins (Sandkvist, 2001). Based on our in silico analysis, W. chitiniclastica harbors, for example, a copy of the xcpT gene, which encodes for pseudopilin XcpT of the T2S machinery of P. aeruginosa (Sandkvist, 2001; Durand et al., 2003). Notably, assembly of the type II pseudopilus also confers increased bacterial adhesive capabilities (Durand et al., 2003), which raises the question, whether its function in W. chitiniclastica is primarily involved in adhesion and/or endotoxin secretion. Unfortunately, the analysis regarding the toxin profile of W. chitiniclastica provided limited information. Although some isolates harbor the toxin-encoding gene relG, which is described to inhibit mycobacterial growth when expressed independently (Korch et al., 2015), other exotoxin encoding genes appear to be missing or are yet unknown. Based on this observation, W. chitiniclastica seems to contain a limited toxin profile when compared to other prolific toxin-producing organisms like C. perfringens (Kiu et al., 2017), making the participation of XcpT in bacterial cell adhesion feasible. However, T2S modules are not restricted to exotoxin secretion; in fact, they can export a wide range of substances. For example, in V. cholerae the T2S system supports the secretion of cholera toxin, hemagglutinin-protease, and chitinase (Connell et al., 1998; Davis et al., 2000; Sandkvist, 2001). Notably, W. chitiniclastica is known to have strong chitinase activity (Tóth et al., 2008) suggesting that T2S systems might be involved in its secretion.
In Silico Analysis of Antimicrobial Resistance Genes
Previous studies have reported W. chitiniclastica to be susceptible to the majority of known antibiotics with the exception of fosfomycin (Schröttner et al., 2017; Matos et al., 2019; Kopf et al., 2021). This is in line with our in silico analysis, which showed that the majority lacks essential AMR genes suggesting a broad susceptibility against several clinical important antibiotics including β-lactamases and fluoroquinolones. This is congruent with recent case reports, where infections caused by W. chitiniclastica were successfully treated with levofloxacin (Schröttner et al., 2017; Bueide et al., 2021) and cephalosporins (Rebaudet et al., 2009; Campisi et al., 2015; Suryalatha et al., 2015; Snyder et al., 2020; Bueide et al., 2021), respectively. However, it should be noted that BM-Y carries a blaVEB-1 gene cassette, thus conferring resistance to ceftazidime and ampicillin as previously reported (Zhou et al., 2016).
Surprisingly, the previously reported fosfomycin resistance (Schröttner et al., 2017; Matos et al., 2019; Kopf et al., 2021) is not reflected within the core resistome profile. DSM 108045 harbors the MFS transporter gene abaF, which is described to confer resistance to fosfomycin (Sharma et al., 2017), but apart from several hits for multidrug efflux proteins within the core genome, we did not detect any known fosfomycin resistance genes such fosA, fosC, or fomB (Silver, 2017), to name but a few. Our results rather indicate a natural fosfomycin resistance most likely based on a yet unknown resistance mechanism as previously anticipated (Kopf et al., 2021). Macrolide resistance genes macA and macB, on the other hand, are encoded in the core genome. Unfortunately, there are no case studies available, that either support or deny our observation. Nevertheless, based on our in silico analysis macrolide resistance appears to be feasible, although further research is still required to uncover the macrolide resistance profile fully. Moreover, W. chitiniclastica appears to be resistant to tellurite, which is not surprising, since potassium tellurite was used intensively as an antimicrobial agent in the past, and as a consequence, many Gram-positive and Gram-negative bacteria developed resistance (Valková et al., 2007). Notably, tellurite resistance genes have also been reported to increase oxidative stress resistance in bacteria (Valková et al., 2007), which is another explanation for their core genome presence in W. chitiniclastica.
Additional genes involved in antimicrobial defense are distributed within the accessory genome indicating resistome expansion due to enormous selective pressure. For example, in silico analysis with respect to aminoglycosides revealed putative resistance in BM-Y, F6513, F6514, F9188 and G9145. Moreover, four isolates harbor sulfonamide resistance genes, while the remaining 22 strains appear to be susceptible. This is in line with recent case reports, in which W. chitiniclastica was susceptible to diverse antibiotics including trimethoprim/sulfamethoxazole (Chavez et al., 2017; Katanami et al., 2018; Connelly et al., 2019; Snyder et al., 2020; Bueide et al., 2021), while the first reported case in South Africa surprised with trimethoprim/sulfamethoxazole resistance (Hoffmann et al., 2016). A similar picture was observed for tetracycline resistance, which showed a rather diverse distribution among the isolates. This observation is also reflected in recent case and research studies, in which some isolates were susceptible to tetracyclines (Almuzara et al., 2011; Nogi et al., 2016), and some resistant (Snyder et al., 2020; Kopf et al., 2021). The presence of Tn10 encoded tetracycline resistance genes tetR and tetC (Babakhani and Oloomi, 2018) further emphasizes the assumption that the majority of resistance genes within W. chitinclastica genomes were required via HGT. Notably, our data indicate that the type strain DSM 18708T does not encode any additional clinical relevant AMR genes, while other strains harbor comparatively more. In particular BM-Y, F6513, F6514, F9188 and G9145 acquired an extended AMR profile; however, this is still limited when compared to other Gram-negative pathogens such as Acinetobacter lwoffii (Hu et al., 2011). Nevertheless, there is an increasing incidence of drug resistance within the W. chitiniclastica clade, whose development should be observed with caution.
Conclusion
The present study provides novel insights on the genetic diversity and pan-genome composition of W. chitiniclastica a rare but potential new emerging human pathogen. Our analysis of all publicity available strains indicate a surprisingly conserved pan-genome without clear host or geographical clustering suggesting a potential spread and transmission. However, with an increasing number of available genomes, reanalysis is strongly recommended to confirm and/or reevaluate the findings. In silico genome studies revealed first insights into genomic features including putative virulence factors and AMR genes, that potentially influence pathogenicity. Interestingly, no clear toxin profile could be determined suggesting an alternative virulence profile. Our results could offer advantages in order to identify potential candidates for target specific manipulations and experimental studies to gain deeper insight into the pathogenic lifestyle of this emerging pathogen. With regard to empirical antibiotic therapy, no general validity for the species W. chitiniclastica can yet be derived from this study. With increasing numbers of available strains, preferably from different regions, but especially with a clear medical history, the analysis should be repeated and even extended to confirm and/or re-evaluate the results. Overall, our results provide the first overview of the genetic mechanisms and AMR profile of W. chitiniclastica that has never been presented in this form, laying the foundation for the best possible therapy.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics Statement
The study was approved by the Ethics Committee at the Technical University of Dresden (EK 61022019). Written informed consent was not obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author Contributions
PS had the idea and the concept for the study. AK analyzed the data, and wrote the first version of the manuscript. BB provided the bioinformatic data from the whole genome sequences. BB, TR, SC, FG, and PS contributed text passages for the manuscript. All authors contributed to the revision of the manuscript and approved the present version. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Federal Ministry of Education and Research, Germany (BMBF; ZIK Septomics Research Centre, Translational Septomics, award no. 03Z22JN12 to SC).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors thank Franziska Klann and Stefan Tiede for excellent technical assistance.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2022.912427/full#supplementary-material
References
Acinas, S. G., Marcelino, L. A., Klepac-Ceraj, V., Polz, M. F. (2004). Divergence and Redundancy of 16S rRNA Sequences in Genomes With Multiple Rrn Operons. J. Bacteriol. 186 (9), 2629–2635. doi: 10.1128/JB.186.9.2629-2635.2004
Adékambi, T., Drancourt, M., Raoult, D. (2009). The rpoB Gene as a Tool for Clinical Microbiologists. Trends Microbiol 17(1):37–5. doi: 10.1016/j.tim.2008.09.008
Alcock, B. P., Raphenya, A. R., Lau, T. T. Y., Tsang, K. K., Bouchard, M., Edalatmand, A., et al. (2020). CARD 2020: Antibiotic Resistome Surveillance With the Comprehensive Antibiotic Resistance Database. Nucleic Acids Res. 48 (D1), D517–D525. doi: 10.1093/nar/gkz935
Almuzara, M. N., Palombarani, S., Tuduri, A., Figueroa, S., Gianecini, A., Sabater, L., et al. (2011). First Case of Fulminant Sepsis Due to Wohlfahrtiimonas Chitiniclastica. J. Clin. Microbiol. 49 (6), 2333–2335. doi: 10.1128/JCM.00001-11
Altschul, S. F., Gish, W., Miller, W., Myers, E. W., Lipman, D. J. (1990). Basic Local Alignment Search Tool. J. Mol. Biol. [Internet]. 215 (3), 403–410. doi: 10.1016/S0022-2836(05)80360-2
Anantham, S., Hall, R. M. (2012). PCERC1, a Small, Globally Disseminated Plasmid Carrying the Dfra14 Cassette in the strA Gene of the Sul2-strA-strB Gene Cluster. Microb. Drug Resist. 18 (4), 364–371. doi: 10.1089/mdr.2012.0008
Anantharaman, V., Aravind, L. (2003). New Connections in the Prokaryotic Toxin-Antitoxin Network : Relationship With the Eukaryotic Nonsense-Mediated RNA Decay System. Genome Biol. 4. doi: 10.1186/gb-2003-4-12-r81
Armougom, F. (2009). Exploring Microbial Diversity Using 16s rRNA High-Throughput Methods. J. Comput. Sci. Syst. Biol. 02 (01), 74–92. doi: 10.4172/jcsb.1000019
Arndt, D., Grant, J. R., Marcu, A., Sajed, T., Pon, A., Liang, Y., et al. (2016). PHASTER: A Better, Faster Version of the PHAST Phage Search Tool. Nucleic Acids Res. 44 (W1), W16–W21. doi: 10.1093/nar/gkw387
Aydin, S., Personne, Y., Newire, E., Laverick, R., Russell, O., Roberts, A. P., et al. (2017). Presence of Type I-F CRISPR/Cas Systems Is Associated With Antimicrobial Susceptibility in Escherichia Coli. J. Antimicrob. Chemother. 72 (8), 2213–2218. doi: 10.1093/jac/dkx137
Babakhani, S., Oloomi, M. (2018). Transposons: The Agents of Antibiotic Resistance in Bacteria. J. Basic Microbiol. 58 (11), 905–917. doi: 10.1002/jobm.201800204
Barrangou, R., Fremaux, C., Deveau, H., Richards, M., Boyaval, P., Moineau, S., et al. (2007). CRISPR Provides Acquired Resistance Against Viruses in Prokaryotes. Science 315 (5819), 1709–1712. doi: 10.1126/science.1138140
Bondoso, J., Harder, J., Lage, O. M. (2013). RpoB Gene as a Novel Molecular Marker to Infer Phylogeny in Planctomycetales. Antonie van Leeuwenhoek Int. J. Gen. Mol. Microbiol. 104 (4), 477–488. doi: 10.1007/s10482-013-9980-7
Bondy-Denomy, J., Pawluk, A., Maxwell, K. L., Davidson, A. R. (2013). Bacteriophage Genes That Inactivate the CRISPR/Cas Bacterial Immune System. Nat. [Internet]. 493 (7432), 429–432. doi: 10.1038/nature11723
Bosch, M., Garrido, E., Llagostera, M., Pérez de Rozas, A. M., Badiola, I., Barbé, J. (2002). Pasteurella Multocida Exbb, exbD and tonB Genes Are Physically Linked But Independently Transcribed. FEMS Microbiol. Lett. 210 (2), 201–208. doi: 10.1111/j.1574-6968.2002.tb11181.x
Bueide, P., Hunt, J., Bande, D., Guerrero, D. M. (2021). Maggot Wound Therapy Associated With Wohlfahrtiimonas Chitiniclastica Blood Infection. Cureus. 13 (1), 10–13. doi: 10.7759/cureus.12471
Buzzanca, D., Botta, C., Ferrocino, I., Alessandria, V., Houf, K., Rantsiou, K. (2021). Genomics Functional Pangenome Analysis Reveals High Virulence Plasticity of Aliarcobacter Butzleri and Affinity to Human Mucus. Genomics [Internet]. 113 (4), 2065–2076. doi: 10.1016/j.ygeno.2021.05.001
Cai, L., Liu, G., Rensing, C., Wang, G. (2009). Genes Involved in Arsenic Transformation and Resistance Associated With Different Levels of Arsenic-Contaminated Soils. BMC Microbiol. 9, 1–11. doi: 10.1186/1471-2180-9-4
Campisi, L., Mahobia, N., Clayton, J. J. (2015). Wohlfahrtiimonas Chitiniclastica Bacteremia Associated With Myiasis, United Kingdom. Emerg. Infect. Dis. 21 (6), 1068–1069. doi: 10.3201/eid2106.140007
Cantas, L., Suer, K. (2014). Review: The Important Bacterial Zoonoses in “One Health” Concept. Front. Public Heal. ;2 (OCT), 1–8. doi: 10.3389/fpubh.2014.00144
Cao, XMX-M, Chen, T., Xu, L-ZLZ, Yao, L-SLS, Qi, J., Zhang, XLX-L, et al. (2013). Complete Genome Sequence of Wohlfahrtiimonas Chitiniclastica Strain SH04, Isolated From Chrysomya Megacephala Collected From Pudong International Airport in China. Genome Announc. 1 (2), 4–5. doi: 10.1128/genomeA.00119-13
Carbonetti, N. H. (2016). Pertussis Leukocytosis: Mechanisms, Clinical Relevance and Treatment. Pathog. Dis. 74 (7), 1–8. doi: 10.1093/femspd/ftw087
Carlin, A., Shi, W., Dey, S., Rosen, B. P. (1995). The Ars Operon of Escherichia Coli Confers Arsenical and Antimonial Resistance. J. Bacteriol. 177 (4), 981–986. doi: 10.1128/jb.177.4.981-986.1995
Case, R. J., Boucher, Y., Dahllöf, I., Holmström, C., Doolittle, W. F., Kjelleberg, S. (2007). Use of 16S rRNA and rpoB Genes as Molecular Markers for Microbial Ecology Studies. Appl. Environ. Microbiol. 73 (1), 278–288. doi: 10.1128/AEM.01177-06
Chavez, J. A., Alexander, A. J., Balada-Llasat, J. M., Pancholi, P. (2017). A Case of Wohlfahrtiimonas Chitiniclastica Bacteremia in Continental United States. JMM Case Rep. 4 (12), 10–12. doi: 10.1099/jmmcr.0.005134
Chiou, C.-S., Jones, A. L. (1995). Expression and Identification of the strA-strB Gene Pair From Streptomycin-Resistant Erwinia Amylovora. Gene [Internet]. 152 (1), 47–51. doi: 10.1016/0378-1119(94)00721-4
Choudhury, H. G., Cameron, A. D., Iwata, S., Beis, K. (2011). Structure and Mechanism of the Chalcogen-Detoxifying Protein TehB From Escherichia Coli. Biochem. J. 435 (1), 85–91. doi: 10.1042/BJ20102014
Cianciotto, N. P., White, R. C. (2017). Expanding Role of Type II Secretion in Bacterial Pathogenesis and Beyond. Infect. Immun. 85 (5), 11–16. doi: 10.1128/IAI.00014-17
Connell, T. D., Metzger, D. J., Lynch, J., Folster, J. P. (1998). Endochitinase is Transported to the Extracellular Milieu by the Eps- Encoded General Secretory Pathway of Vibrio Cholerae. J. Bacteriol. 180 (21), 5591–5600. doi: 10.1128/JB.180.21.5591-5600.1998
Connelly, K., Freeman, E., Smibert, O., Lin, B. (2019). Wohlfahrtiimonas Chitiniclastica Bloodstream Infection Due to a Maggot-Infested Wound in a 54-Year-Old Male. J. Glob Infect. Dis. [Internet]. 11 (3), 125–126. doi: 10.4103/jgid.jgid_58_18
Costa, S. S., Guimarães, L. C., Silva, A., Soares, S. C., Baraúna, R. A. (2020). First Steps in the Analysis of Prokaryotic Pan-Genomes. Bioinform. Biol. Insights 14. doi: 10.1177/1177932220938064
Couvin, D., Bernheim, A., Toffano-Nioche, C., Touchon, M., Michalik, J., Néron, B., et al. (2018). CRISPRCasFinder, an Update of CRISRFinder, Includes a Portable Version, Enhanced Performance and Integrates Search for Cas Proteins. Nucleic Acids Res. 46 (W1), W246–W251. doi: 10.1093/nar/gky425
Cox, G., Stogios, P. J., Savchenko, A., Wright, G. D. (2015). Structural and Molecular Basis for Resistance to Aminoglycoside Antibiotics by the Adenylyltransferase ANT(2″)-Ia. MBio. 6 (1), 1–9. doi: 10.1128/mBio.02180-14
Davis, B. M., Lawson, E. H., Sandkvist, M., Ali, A., Sozhamannan, S., Waldor, M. K. (2000). Convergence of the Secretory Pathways for Cholera Toxin and the Filamentous Phage, CTXphi. Science. 288 (5464), 333–335. doi: 10.1126/science.288.5464.333
Diaz-Delgado, J., Eva, S., Isabel, V. A., Lucas, D., Marisa, A., Manuel, A., et al. (2015). Endocarditis Associated With Wohlfahrtiimonas Chitiniclastica in a Short-Beaked Common Dolphin (Delphinus Delphis). J. Wildl Dis. 51 (1), 283–286. doi: 10.7589/2014-03-072
Durand, É, Bernadac, A., Ball, G., Lazdunski, A., Sturgis, J. N., Filloux, A. (2003). Type II Protein Secretion in Pseudomonas Aeruginosa: The Pseudopilus Is a Multifibrillar and Adhesive Structure. J. Bacteriol. 185 (9), 2749–2758. doi: 10.1128/JB.185.9.2749-2758.2003
Eicher, T., Brandstätter, L., Pos, K. M. (2009). Structural and Functional Aspects of the Multidrug Efflux Pump AcrB. Biol. Chem. 390 (8), 693–699. doi: 10.1515/BC.2009.090
Faure, G., Makarova, K. S., Koonin, E. V. (2019). CRISPR–Cas: Complex Functional Networks and Multiple Roles Beyond Adaptive Immunity. J. Mol. Biol. [Internet]. 431 (1), 3–20. doi: 10.1016/j.jmb.2018.08.030
Fekih, I. B., Zhang, C., Li, Y. P., Zhao, Y., Alwathnani, H. A., Saquib, Q., et al. (2018). Distribution of Arsenic Resistance Genes in Prokaryotes. Front. Microbiol. 9 (OCT), 1–11. doi: 10.3389/fmicb.2018.02473
Fino, C., Vestergaard, M., Ingmer, H., Pierrel, F., Gerdes, K., Harms, A. (2020). PasT of Escherichia Coli Sustains Antibiotic Tolerance and Aerobic Respiration as a Bacterial Homolog of Mitochondrial Coq10. Microbiologyopen. 9 (8), 1–36. doi: 10.1002/mbo3.1064
Grissa, I., Vergnaud, G., Pourcel, C. (2007). The CRISPRdb Database and Tools to Display CRISPRs and to Generate Dictionaries of Spacers and Repeats. BMC Bioinf. 8, 1–10. doi: 10.1186/1471-2105-8-172
Hadfield, J., Croucher, N. J., Goater, R. J., Abudahab, K., Aanensen, D. M., Harris, S. R. (2018). Phandango: An Interactive Viewer for Bacterial Population Genomics. Bioinformatics. 34 (2), 292–293. doi: 10.1093/bioinformatics/btx610
Hassan, K. A., Liu, Q., Elbourne, L. D. H., Ahmad, I., Sharples, D., Naidu, V., et al. (2018). Pacing Across the Membrane: The Novel PACE Family of Efflux Pumps Is Widespread in Gram-Negative Pathogens. Res. Microbiol. 169 (7–8), 450–454. doi: 10.1016/j.resmic.2018.01.001
Hoffmann, R., Fortuin, F., Newton-Foot, M., Singh, S. (2016). First Report of Wohlfahrtiimonas Chitiniclastica Bacteraemia in South Africa. SAMJ South Afr. Med. J. [Internet]. 106, 1062. doi: 10.7196/SAMJ.2016.v106i11.11449
Holden, K. M., Browning, G. F., Noormohammadi, A. H., Markham, P. F., Marenda, M. S. (2012). TonB Is Essential for Virulence in Avian Pathogenic Escherichia Coli. Comp. Immunol. Microbiol. Infect. Dis. [Internet]. 35 (2), 129–138. doi: 10.1016/j.cimid.2011.12.004
Hu, Y., Zhang, W., Liang, H., Liu, L., Peng, G., Pan, Y., et al. (2011). Whole-Genome Sequence of a Multidrug-Resistant Clinical Isolate of Acinetobacter Lwoffii. J. Bacteriol. 193 (19), 5549–5550. doi: 10.1128/JB.05617-11
Janda, J. M., Abbott, S. L. (2007). 16s rRNA Gene Sequencing for Bacterial Identification in the Diagnostic Laboratory: Pluses, Perils, and Pitfalls. J. Clin. Microbiol. 45 (9), 2761–2764. doi: 10.1128/JCM.01228-07
Jordan, I. K., Makarova, K. S., Spouge, J. L., Wolf, Y. I., Koonin, E. V. (2001). Lineage-Specific Gene Expansions in Bacterial and Archaeal Genomes. Genome Res. 11 (4), 555–565. doi: 10.1101/gr.166001
Kaatz, G. W., DeMarco, C. E., Seo, S. M. (2006). MepR, a Represser of the Staphylococcus Aureus MATE Family Multidrug Efflux Pump MepA, Is a Substrate-Responsive Regulatory Protein. Antimicrob. Agents Chemother. 50 (4), 1276–1281. doi: 10.1128/AAC.50.4.1276-1281.2006
Katanami, Y., Kutsuna, S., Nagashima, M., Takaya, S., Yamamoto, K., Takeshita, N., et al. (2018). Wohlfahrtiimonas Chitiniclastica Bacteremia Hospitalized Homeless Man With Squamous Cell Carcinoma. Emerg. Infect. Dis. [Internet]. 24 (9), 1746–1748. doi: 10.3201/eid2409.170080
Kirby, J. E., Trempy, J. E., Gottesman, S. (1994). Excision of a P4-Like Cryptic Prophage Leads to Alp Protease Expression in Escherichia Coli. J. Bacteriol. 176 (7), 2068–2081. doi: 10.1128/jb.176.7.2068-2081.1994
Kiu, R., Caim, S., Alexander, S., Pachori, P., Hall, L. J. (2017). Probing Genomic Aspects of the Multi-Host Pathogen Clostridium Perfringens Reveals Significant Pangenome Diversity, and a Diverse Array of Virulence Factors. Front. Microbiol. 8 (DEC). doi: 10.3389/fmicb.2017.02485
Klindworth, A., Peplies, J., Pruesse, E., Schweer, T., Glöckner, F. O., Quast, C., et al. (2012). Evaluation of General 16s Ribosomal RNA Gene PCR Primers for Classical and Next-Generation Sequencing-Based Diversity Studies. Nucleic Acids Res. 41 (1), e1–e1. doi: 10.1093/nar/gks808
Kobayashi, S., Kuzuyama, T., Seto, H. (2000). Characterization of the fomA and fomB Gene Products From Streptomyces Wedmorensis, Which Confer Fosfomycin Resistance on Escherichia Coli. Antimicrob. Agents Chemother. 44 (3), 647–650. doi: 10.1128/AAC.44.3.647-650.2000
Kõljalg, S., Telling, K., Huik, K., Murruste, M., Saarevet, V., Pauskar, M., et al. (2015). First Report of Wohlfahrtiimonas Chitiniclastica From Soft Tissue and Bone Infection at an Unusually High Northern Latitude. Folia Microbiol. (Praha) [Internet]. 60 (2), 155–158. doi: 10.1007/s12223-014-0355-x
Kopf, A., Bunk, B., Coldewey, S. M., Gunzer, F., Riedel, T., Schröttner, P. (2021). Identification and Antibiotic Profiling of Wohlfahrtiimonas Chitiniclastica, an Underestimated Human Pathogen. Front. Microbiol. 12 (September). doi: 10.3389/fmicb.2021.712775
Korch, S. B., Malhotra, V., Contreras, H., Clark-Curtiss, J. E. (2015). The Mycobacterium Tuberculosis relBE Toxin:Antitoxin Genes Are Stress-Responsive Modules That Regulate Growth Through Translation Inhibition. J. Microbiol. 53 (11), 783–795. doi: 10.1007/s12275-015-5333-8
Kottara, A., Hall, J. P. J., Harrison, E., Brockhurst, M. A. (2018). Variable Plasmid Fitness Effects and Mobile Genetic Element Dynamics Across Pseudomonas Species. FEMS Microbiol. Ecol. 94 (1), 1–7. doi: 10.1093/femsec/fix172
Landsberger, M., Gandon, S., Meaden, S., Rollie, C., Chevallereau, A., Chabas, H., et al. (2018). Anti-CRISPR Phages Cooperate to Overcome CRISPR-Cas Immunity. Cell. 174 (4), 908–916.e12. doi: 10.1016/j.cell.2018.05.058
Lefort, V., Desper, R., Gascuel, O. (2015). FastME 2.0: A Comprehensive, Accurate, and Fast Distance-Based Phylogeny Inference Program. Mol. Biol. Evol. 32 (10), 2798–2800. doi: 10.1093/molbev/msv150
Lobato-Márquez, D., Díaz-Orejas, R., García-del Portillo, F. (2016). Toxin-Antitoxins and Bacterial Virulencea. FEMS Microbiol. Rev. 40 (5), 592–609. doi: 10.1093/femsre/fuw022
Lysaght, T. B., Wooster, M. E., Jenkins, P. C., Koniaris, L. G. (2018). Myiasis-Induced Sepsis: A Rare Case Report of Wohlfahrtiimonas Chitiniclastica and Ignatzschineria Indica Bacteremia in the Continental United States. Med. (Baltimore). 97 (52), e13627. doi: 10.1097/MD.0000000000013627
Ma, D., Gu, H., Shi, Y., Huang, H., Sun, D., Hu, Y. (2021). Edwardsiella Piscicida YefM-YoeB: A Type II Toxin-Antitoxin System That Is Related to Antibiotic Resistance, Biofilm Formation, Serum Survival, and Host Infection. Front. Microbiol. 12 (March), 1–15. doi: 10.3389/fmicb.2021.646299
Maleki-Ravasan, N., Ahmadi, N., Soroushzadeh, Z., Raz, A. A., Zakeri, S., Dinparast Djadid, N. (2020). New Insights Into Culturable and Unculturable Bacteria Across the Life History of Medicinal Maggots Lucilia Sericata (Meigen) (Diptera: Calliphoridae). Front. Microbiol. 11 (April), 1–17. doi: 10.3389/fmicb.2020.00505
Marino, N. D., Zhang, J. Y., Rorges, A. L., Sousa, A. A., Leon, L. M., Rauch, B. J., et al. (2018). Discovery of Widespread Type I and Type V CRISPR-Cas Inhibitors. Sci. (80- ). 07415 (October), 240–242. doi: 10.1126/science.aau5174
Matos, J., Faria, A. R., Carvalho Assef, A. P. D., de Freitas-Almeida, ÂC, Albano, R. M., Queiroz, M. L. P. (2019). Draft Genome Sequence of a Wohlfahrtiimonas Chitiniclastica Strain Isolated From Frozen Chicken in Rio De Janeiro, Brazil. Microbiol. Resour Announc. 8 (49), 1–2. doi: 10.1128/MRA.00352-19
Matsumoto, Y., Kaito, C., Morishita, D., Kurokawa, K., Sekimizu, K. (2007). Regulation of Exoprotein Gene Expression by the Staphylococcus Aureus cvfB Gene. Infect. Immun. 75 (4), 1964–1972. doi: 10.1128/IAI.01552-06
McInerney, J. O., McNally, A., O’Connell, M. J. (2017). Why Prokaryotes Have Pangenomes. Nat. Microbiol. 2 (March). doi: 10.1038/nmicrobiol.2017.40
Meier-Kolthoff, J. P., Auch, A. F., Klenk, H.-P. P., Göker, M. (2013). Genome Sequence-Based Species Delimitation With Confidence Intervals and Improved Distance Functions. BMC Bioinf. [Internet]. 14 (1), 60. doi: 10.1186/1471-2105-14-60
Meier-Kolthoff, J. P., Göker, M. (2019). TYGS is an Automated High-Throughput Platform for State-Of-the-Art Genome-Based Taxonomy. Nat. Commun. [Internet]. 10 (1), 2182. doi: 10.1038/s41467-019-10210-3
Meier-Kolthoff, J. P., Hahnke, R. L., Petersen, J., Scheuner, C., Michael, V., Fiebig, A., et al. (2014). Complete Genome Sequence of DSM 30083T, the Type Strain (U5/41T) of Escherichia Coli, and a Proposal for Delineating Subspecies in Microbial Taxonomy. Stand Genomic Sci. [Internet]. 9 (1), 2. doi: 10.1186/1944-3277-9-2
Morton, D. J., Hempel, R. J., Seale, T. W., Whitby, P. W., Stull, T. L. (2012). A Functional tonB Gene Is Required for Both Virulence and Competitive Fitness in a Chinchilla Model of Haemophilus Influenzae Otitis Media. BMC Res. Notes. 5. doi: 10.1186/1756-0500-5-327
Nagakubo, S., Nishino, K., Hirata, T., Yamaguchi, A. (2002). The Putative Response Regulator BaeR Stimulates Multidrug Resistance of Escherichia Coli via a Novel Multidrug Exporter System, MdtABC. J. Bacteriol. 184 (15), 4161–4167. doi: 10.1128/JB.184.15.4161-4167.2002
Nishino, K., Latifi, T., Groisman, E. A. (2006). Virulence and Drug Resistance Roles of Multidrug Efflux Systems of Salmonella Enterica Serovar Typhimurium. Mol. Microbiol. 59 (1), 126–141. doi: 10.1111/j.1365-2958.2005.04940.x
Nogi, M., Bankowski, M. J., Pien, F. D. (2016). Wohlfahrtiimonas Chitiniclastica Infections in 2 Elderly Patients, Hawaii, USA. Emerg. Infect. Dis. 22 (3), 567–568. doi: 10.3201/eid2203.151701
Nordmann, P., Naas, T. (1994). Sequence Analysis of PER-1 Extended-Spectrum β-Lactamase From Pseudomonas Aeruginosa and Comparison With Class A β-Lactamases. Antimicrob. Agents Chemother. 38 (1), 104–114. doi: 10.1128/AAC.38.1.104
Norton, J. P., Mulvey, M. A. (2012). Toxin-Antitoxin Systems Are Important for Niche-Specific Colonization and Stress Resistance of Uropathogenic Escherichia Coli. PLoS Pathog. 8 (10). doi: 10.1371/journal.ppat.1002954
Page, A. J., Cummins, C. A., Hunt, M., Wong, V. K., Reuter, S., Holden, M. T. G., et al. (2015). Roary: Rapid Large-Scale Prokaryote Pan Genome Analysis. Bioinformatics. 31 (22), 3691–3693. doi: 10.1093/bioinformatics/btv421
Panis, G., Duverger, Y., Courvoisier-Dezord, E., Champ, S., Talla, E., Ansaldi, M. (2010). Tight Regulation of the intS Gene of the KplE1 Prophage: A New Paradigm for Integrase Gene Regulation. PLoS Genet. [Internet]. 6 (10), e1001149. doi: 10.1371/journal.pgen.1001149
Parent, R., Roy, P. H. (1992). The Chloramphenicol Acetyltransferase Gene of Tn2424: A New Breed of Cat. J. Bacteriol. 174 (9), 2891–2897. doi: 10.1128/jb.174.9.2891-2897.1992
Petkau, A. (2018). Staramr (GitHub). Available at: https://github.com/phac-nml/staramr.
Popovic, A., Wu, B., Arrowsmith, C. H., Edwards, A. M., Davidson, A. R., Maxwell, K. L. (2012). Structural and Biochemical Characterization of Phage λ FI Protein (gpFI) Reveals a Novel Mechanism of DNA Packaging Chaperone Activity. J. Biol. Chem. 287 (38), 32085–32095. doi: 10.1074/jbc.M112.378349
Qi, J., Gao, Y., Wang, G., Li, L., Li, L., Zhao, X., et al. (2016). Identification of Wohlfahrtiimonas Chitiniclastica Isolated From an Infected Cow With Hoof Fetlow, China. Infect. Genet. Evol. [Internet]. 41, 174–176. doi: 10.1016/j.meegid.2016.04.008
Rainey, F. A., Ward-Rainey, N. L., Janssen, P. H., Hippe, H., Stackebrandt, E. (1996). Clostridium Paradoxum DSM 7308t Contains Multiple 16s rRNA Genes With Heterogeneous Intervening Sequences. Microbiology. 142, 2087–2095. doi: 10.1099/13500872-142-8-2087
Rebaudet, S., Genot, S., Renvoise, A., Fournier, P. E., Stein, A. (2009). Wohlfahrtiimonas Chitiniclastica Bacteremia in Homeless Woman. Emerg. Infect. Dis. [Internet]. 15 (6), 985–987. doi: 10.3201/eid1506.080232
Reddy, M. R. K., Mastan, S. A. (2013). Wohlfahrtiimonas Chitiniclastica Fulminant Sepsis in Pangasius Sutchi-First Report. Turkish J. Fish Aquat Sci. 13 (4), 753–758. doi: 10.4194/1303-2712-v13_4_21
Reeves, S. A., Torres, A. G., Payne, S. M. (2000). TonB is Required for Intracellular Growth and Virulence of Shigella Dysenteriae. Infect. Immun. 68 (11), 6329–6336. doi: 10.1128/IAI.68.11.6329-6336.2000
Robbins, K., Khachemoune, A. (2010). Cutaneous Myiasis: A Review of the Common Types of Myiasis. Int. J. Dermatol. [Internet]. 49 (10), 1092–1098. doi: 10.1111/j.1365-4632.2010.04577.x
Roberts, M. C. (2005). Update on Acquired Tetracycline Resistance Genes. FEMS Microbiol. Lett. 245 (2), 195–203. doi: 10.1016/j.femsle.2005.02.034
Rouli, L., Merhej, V., Fournier, P. E., Raoult, D. (2015). The Bacterial Pangenome as a New Tool for Analysing Pathogenic Bacteria. N. Microbes N. Infect. [Internet]. 7, 72–85. doi: 10.1016/j.nmni.2015.06.005
Roux, S., Enault, F., le Bronner, G., Debroas, D. (2011). Comparison of 16S rRNA and Protein-Coding Genes as Molecular Markers for Assessing Microbial Diversity (Bacteria and Archaea) in Ecosystems. FEMS Microbiol. Ecol. 78 (3), 617–628. doi: 10.1111/j.1574-6941.2011.01190.x
Rudd, K. E. (1999). Novel Intergenic Repeats of Escherichia Coli K-12. Res. Microbiol. 150 (9–10), 653–664. doi: 10.1016/S0923-2508(99)00126-6
Sandkvist, M. (2001). Type II Secretion and Pathogenesis. Infect. Immun. 69 (6), 3523–3535. doi: 10.1128/IAI.69.6.3523-3535.2001
Sanyal, S. K., Mou, T. J., Chakrabarty, R. P., Hoque, S., Hossain, M. A., Sultana, M. (2016). Diversity of Arsenite Oxidase Gene and Arsenotrophic Bacteria in Arsenic Affected Bangladesh Soils. AMB Express. 6 (1). doi: 10.1186/s13568-016-0193-0
Schirmer, M., Ijaz, U. Z., D’Amore, R., Hall, N., Sloan, W. T., Quince, C. (2015). Insight Into Biases and Sequencing Errors for Amplicon Sequencing With the Illumina MiSeq Platform. Nucleic Acids Res. 43 (6). doi: 10.1093/nar/gku1341
Schröttner, P., Rudolph, W. W., Damme, U., Lotz, C., Jacobs, E., Gunzer, S. (2017). Wohlfahrtiimonas Chitiniclastica: Current Insights Into an Emerging Human Pathogen. Epidemiol. Infect. 145 (7), 1292–1303. doi: 10.1017/S0950268816003411
Seemann, T. (2014). Prokka: Rapid Prokaryotic Genome Annotation. Bioinformatics. 30 (14), 2068–2069. doi: 10.1093/bioinformatics/btu153
Sharma, A., Sharma, R., Bhattacharyya, T., Bhando, T., Pathania, R. (2017). Fosfomycin Resistance in Acinetobacter Baumannii Is Mediated by Efflux Through a Major Facilitator Superfamily (MFS) Transporter-AbaF. J. Antimicrob. Chemother. 72 (1), 68–74. doi: 10.1093/jac/dkw382
Silvaggi, J. M., Perkins, J. B., Losick, R. (2005). Small Untranslated RNA Antitoxin in Bacillus Subtilis. J. Bacteriol. 187 (19), 6641–6650. doi: 10.1128/JB.187.19.6641-6650.2005
Silver, L. L. (2017). Fosfomycin: Mechanism and Resistance. Cold Spring Harb. Perspect. Med. 7 (2), 1–11. doi: 10.1101/cshperspect.a025262
Smith, L. D., Bertrand, K. P. (1988). Mutations in the Tn10 Tet Repressor That Interfere With Induction. Location of the Tetracycline-Binding Domain. J. Mol. Biol. 203 (4), 949–959. doi: 10.1016/0022-2836(88)90120-9
Snyder, S., Singh, P., Goldman, J. (2020). Emerging Pathogens: A Case of Wohlfahrtiimonas Chitiniclastica and Ignatzschineria Indica Bacteremia. IDCases [Internet]. 19, e00723. doi: 10.1016/j.idcr.2020.e00723
Stajich, J. E., Block, D., Boulez, K., Brenner, S. E., Chervitz, S. A., Dagdigian, C., et al. (2002). The Bioperl Toolkit: Perl Modules for the Life Sciences. Genome Res. 12 (10), 1611–1618. doi: 10.1101/gr.361602
Stalder, T., Press, M. O., Sullivan, S., Liachko, I., Top, E. M. (2019). Linking the Resistome and Plasmidome to the Microbiome. ISME J. [Internet]. 13 (10), 2437–2446. doi: 10.1038/s41396-019-0446-4
Stevenson, B. S., Schmidt, T. M. (2004). Life History Implications of rRNA Gene Copy Number in Escherichia Coli. Appl. Environ. Microbiol. 70 (11), 6670–6677. doi: 10.1128/AEM.70.11.6670-6677.2004
Stewart, G. C., Wilson, F. E., Bott, K. F. (1982). Detailed Physical Mapping of the Ribosomal RNA Genes of Bacillus Subtilis. Gene. 19 (2), 153–162. doi: 10.1016/0378-1119(82)90001-4
Sune, D., Rydberg, H., Augustinsson, ÅN, Serrander, L., Jungeström, M. B. (2020). Optimization of 16S rRNA Gene Analysis for Use in the Diagnostic Clinical Microbiology Service. J. Microbiol. Methods 170 (December 2019), 105854. doi: 10.1016/j.mimet.2020.105854
Suryalatha, K., John, J., Thomas, S. (2015). Wohlfahrtiimonas Chitiniclastica-Associated Osteomyelitis: A Rare Case Report. Future Microbiol. [Internet]. 10 (7), 1107–1109. doi: 10.2217/fmb.15.44
Tatusova, T., DiCuccio, M., Badretdin, A., Chetvernin, V., Nawrocki, E. P., Zaslavsky, L., et al. (2016). NCBI Prokaryotic Genome Annotation Pipeline. Nucleic Acids Res. 44 (14), 6614–6624. doi: 10.1093/nar/gkw569
Tauch, A., Krieft, S., Kalinowski, J., Pühler, A. (2000). The 51,409-Bp R-Plasmid Ptp10 From the Multiresistant Clinical Isolate Corynebacterium Striatum M82B Is Composed of DNA Segments Initially Identified in Soil Bacteria and in Plant, Animal, and Human Pathogens. Mol. Gen. Genet. 263 (1), 1–11. doi: 10.1007/PL00008668
Tettelin, H., Riley, D., Cattuto, C., Medini, D. (2008). Comparative Genomics: The Bacterial Pan-Genome. Curr. Opin. Microbiol. 11 (5), 472–477. doi: 10.1016/j.mib.2008.09.006
Thaiwong, T., Kettler, N. M., Lim, A., Dirkse, H., Kiupel, M. (2014). First Report of Emerging Zoonotic Pathogen Wohlfahrtiimonas Chitiniclastica in the United States. J. Clin. Microbiol. 52 (6), 2245–2247. doi: 10.1128/JCM.00382-14
Thomas, C. E., Olsen, B., Elkins, C. (1998). Cloning and Characterization of Tdha, a Locus Encoding a TonB-Dependent Heme Receptor From Haemophilus Ducreyi. Infect. Immun. 66 (9), 4254–4262. doi: 10.1128/IAI.66.9.4254-4262.1998
Torres, A. G., Redford, P., Welch, R. A., Payne, S. M. (2001). TonB-Dependent Systems of Uropathogenic Escherichia Coli: Aerobactin and Heme Transport and TonB are Required for Virulence in the Mouse. Infect. Immun. [Internet]. 69 (10), 6179–6185. doi: 10.1128/IAI.69.10.6179-6185.2001
Tóth, E. M., Schumann, P., Borsodi, A. K., Kéki, Z., Kovács, A. L., Márialigeti, K. (2008). Wohlfahrtiimonas Chitiniclastica Gen. Nov., Sp. Nov., a New Gammaproteobacterium Isolated From Wohlfahrtia Magnifica (Diptera: Sarcophagidae). Int. J. Syst. Evol. Microbiol. 58 (4), 976–981. doi: 10.1099/ijs.0.65324-0
Valková, D., Valkovičová, L., Vávrová, S., Kováčová, E., Mravec, J., Turňa, J. (2007). The Contribution of Tellurite Resistance Genes to the Fitness of Escherichia Coli Uropathogenic Strains. Cent Eur. J. Biol. 2 (2), 182–191. doi: 10.2478/s11535-007-0019-9
van Belkum, A., Soriaga, L. B., Lafave, M. C., Akella, S., Veyrieras, J., Barbu, E. M., et al. (2015). Phylogenetic Distribution of CRISPR-Cas Systems in Antibiotic-Resistant Pseudomonas aeruginosa. mBio. 6, 6, 1–13. doi: 10.1128/mBio.01796-15
Vassort-Bruneau, C., Lesage-Descauses, M. C., Martel, J. L., Lafont, J. P., Chaslus-Dancla, E. (1996). CAT III Chloramphenicol Resistance in Pasteurella Haemolytica and Pasteurella Multocida Isolated From Calves. J. Antimicrob. Chemother. 38 (2), 205–213. doi: 10.1093/jac/38.2.205
Welch, T. J., Fricke, W. F., McDermott, P. F., White, D. G., Rosso, M.-L., Rasko, D. A., et al. (2007). Multiple Antimicrobial Resistance in Plague: An Emerging Public Health Risk. PLoS One 2 (3), e309. doi: 10.1371/journal.pone.0000309
Westra, E. R., Buckling, A., Fineran, P. C. (2014). CRISPR-Cas Systems: Beyond Adaptive Immunity. Nat. Rev. Microbiol. 12 (5), 317–326. doi: 10.1038/nrmicro3241
Yang, H.-C., Fu, H.-L., Lin, Y.-F., Rosen, B. P., Lutsenko SBT-CT in M (2012). “Chapter Twelve - Pathways of Arsenic Uptake and Efflux,” in Metal Transporters. Ed. Argüello, J. M. (Academic Press), 325–358.
Yi, H., Huang, L., Yang, B., Gomez, J., Zhang, H., Yin, Y. (2020). AcrFinder: Genome Mining Anti-CRISPR Operons in Prokaryotes and Their Viruses. Nucleic Acids Res. 48 (W1), W358–W365. doi: 10.1093/nar/gkaa351
Yum, S., Xu, Y., Piao, S., Sim, S.-H., Kim, H.-M., Jo, W.-S., et al. (2009). Crystal Structure of the Periplasmic Component of a Tripartite Macrolide-Specific Efflux Pump. J. Mol. Biol. 387 (5), 1286–1297. doi: 10.1016/j.jmb.2009.02.048
Zhou, W., Li, M., Zhu, L., Hua, F., Ji, X., Sun, Y., et al. (2016). Complete Genome Sequence of Wohlfahrtiimonas Chitiniclastica Strain BM-Y, Isolated From the Pancreas of a Zebra in China. Genome Announc. 4 (3), 2015–2016. doi: 10.1128/genomeA.00643-16
Keywords: W. chitiniclastica, pan-genome, antimicrobial resistance, CRISPR, Acr (anti-CRISPR), rpoB, arsenic resistance
Citation: Kopf A, Bunk B, Coldewey SM, Gunzer F, Riedel T and Schröttner P (2022) Comparative Genomic Analysis of the Human Pathogen Wohlfahrtiimonas Chitiniclastica Provides Insight Into the Identification of Antimicrobial Resistance Genotypes and Potential Virulence Traits. Front. Cell. Infect. Microbiol. 12:912427. doi: 10.3389/fcimb.2022.912427
Received: 04 April 2022; Accepted: 02 June 2022;
Published: 07 July 2022.
Edited by:
Yang Zhang, University of Pennsylvania, United StatesReviewed by:
Yu Zhou, Institut Pasteur of Shanghai, Chinese Academy of Sciences (CAS), ChinaQiaosi Tang, University of Pennsylvania, United States
Copyright © 2022 Kopf, Bunk, Coldewey, Gunzer, Riedel and Schröttner. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Percy Schröttner, cGVyY3kuc2Nocm9ldHRuZXJAdHUtZHJlc2Rlbi5kZQ==