Evidence of Extensive Circulation of Yersinia enterocolitica in Rodents and Shrews in Natural Habitats from Retrospective and Perspective Studies in South Caucasus

Yersinia enterocolitica culture-positive rodents and shrews were reported in different territories across Georgia during 14 of 17 years of investigations conducted for the period of 1981–1997. In total, Y. enterocolitica was isolated from 2052 rodents (15 species) and 33 shrews. Most isolates were obtained from Microtus arvalis, Rattus norvegicus, Mus musculus, and Apodemus spp. During the prospective study (2017−2019), isolates of Yersinia-like bacteria were cultured from 53 rodents collected in four parts of Georgia. All the Yersinia-like isolates were confirmed as Y. enterocolitica based on the API 20E and the BD Phenix50 tests. Whole-genome (WG) sequencing of five rodents and one shrew strain of Y. enterocolitica revealed that they possessed a set of virulence genes characteristic of the potentially pathogenic strains of biogroup 1A. All isolates lacked distinguished virulence determinants for YstA, Ail, TccC, VirF, and virulence plasmid pYV but carried the genes for YstB, YmoA, HemPR-HmuVSTU, YaxAB, PhlA, PldA, ArsCBR, and a flagellar apparatus. One strain contained a gene highly homologous to heat-labile enterotoxin, a chain of E. coli, a function not previously described for Y. enterocolitica. The WG single-nucleotide polymorphism-based typing placed the isolates in four distinct phylogenetic clusters.


Introduction
Yersinia enterocolitica, as a species, belongs to the genus Yersinia within the Enterobacteriaceae family that combines Yersinia pestis, agent of plague, and a large number of other phylogenetically and phenotypically close species and strains [1,2]. Bacteria assigned to Y. enterocolitica represent a heterogeneous and diverse group of strains, but most information about this group arrived from investigations of bacterial cultures obtained from human patients. Because of the origin of these strains, they designated those strains and biogroups as pathogenic. In addition, the pathogenicity of the strains is declared on the presence of some genes (virulence markers), which facilitates colonization of human or animal organisms [3,4]. Since only a limited number of the biogroups obtained from human patients have been investigated, a multitude and diversity of Y. enterocolitica strains circulating in nature remain poorly investigated and could be well beyond the scope of the so-called "pathogenic biogroups". Even among investigated variants, the recognition of pathogenic strains is conditional. For example, strains of biogroup 1A were considered non-pathogenic because of the absence of classical virulence markers until the demonstration of their pathogenic potential [5].
In contrast to Y. pestis that causes plague, and similarly to Y. pseudotuberculosis, some biovars of Y. enterocolitica cause gastrointestinal illness and various clinical manifestations in people, generally termed as 'yersiniosis' [3]. The ecology of Y. enterocolitica is very different from the ecology of Y. pestis. While the plague pathogen, Y. pestis, persists in nature through cycles of flea transmission between specific rodents [6,7], Y. enterocolitica is widespread in the environment (e.g., soil) and is commonly transmitted via contaminated food [8]. At the same time, there are many reports of detection of this bacterium in various animals that led to the recognition of this species as a zoonotic agent [8]. The most well-known zoonotic association of Y. enterocolitica is with domestic pigs since contaminated pork is considered to be the most important source for human yersiniosis. Bacteria of this species have been isolated from other domestic animals (cattle, sheep, and goats) and from numerous species of wild animals (wild boars, deer, bats, foxes, jackals, and raccoons) [9][10][11][12][13].
Much less is known about the distribution of Y. enterocolitica in rodents and shrews. Most reports are limited to culturing these bacteria from commensal rats living in close proximity to pig farms [14]. When identified, the isolates obtained from these rats were undistinguished from strains found in pigs. There are reports of isolation of Y. enterocolitica strains from small free-living rodents from China, Japan, Czechoslovakia, Finland, Poland, Scandinavian countries, and Crimea [15][16][17][18][19][20]. Nine strains of Y. enterocolitica were obtained from common shrews in Finland [18]. The strains of Y. enterocolitica identified in natural sources were antigenically variable, which led to referring to some of them as Y. enterocolitica-like bacteria [1]. Classification of the bacteria that are close but different from well-accepted Y. enterocolitica biovars and await further developments, partially because of the limited representation of Y. enterocolitica-like cultures and genomes obtained from animal sources.
Plague has continuously circulated in the Caucasus region for many centuries [21]. Since 1992, scientists in Georgia have continued the investigation of rodents within natural foci of plague. Occasionally, the bacteriologists working in plague foci reported cultures of Y. enterocolitica from investigated rodents, but those reports were either unpublished or presented in archived documents which were previously classified material and unreachable to outside specialists. An additional source for information about the presence of Y. enterocolitica and related bacteria has emerged after the U.S. Defense Threat Reduction Agency (DTRA) supported a special project to investigate the molecular epidemiology and ecology of Yersinia species in the plague endemic territory in Georgia and neighboring Azerbaijan. This paper summarizes results on the prevalence of Y. enterocolitica in populations of different rodent species based on (1) retrospective analysis of available data on the identification of Y. enterocolitica in rodents during routine plague surveillance conducted across the territory of Georgia and (2) a prospective investigation of rodents captured in a plague focus located in the Eastern part of Georgia.

Archived Documents
We have carefully analyzed annual reports of the Republican Anti-plague Station for the period of 1960-1997 archived at the National Center for Disease Control & Public Health (NCDC). Information on the detection of Y. enterocolitica in rodents and shrews was found in reports starting in 1981. The reported analysis covers the period from 1981 to 1997. Only common names of rodents have been provided in the reports, and where possible, the names follow current rodent taxonomy with the use of Latin names. Generally, we referred to the local Caucasian guide [22] and the checklist in [23], with some changes according to modern issues [24]. By doing this, we wish to avoid any controversy and debates between zoologists regarding taxonomic revisions. A selection of some Latin names (e.g., Apodemus versus Sylvaemus for wood mice) was made purely for a conventional purpose and does not reflect our preference for any particular taxonomic scheme. All mice belonging to the genus Apodemus were combined together because of the impossibility of species name verification.

Field Investigations
In total, 573 rodents and shrews were collected during 2017-2019 for projects conducted by NCDC. All procedures for capturing and processing small rodents were approved by the Institutional Animal Care and Use Committee at the Georgian Association for Laboratory Animal Science (Protocol No. PA 005/17). A combination of two kinds of the Sherman live traps (Sherman Traps Inc., Tallahassee, FL, USA) were used to capture rodents: 7.6 cm × 8.9 cm × 22.9 cm aluminum folding and non-folding traps collecting for mice and voles and 12.7 cm × 12.7 cm × 38.1 cm Large Non-Folding Heavy-Duty traps to target rats and jirds. The traps were set out with intervals of approximately 4-6 m during daylight hours, usually late afternoon, over a 1-2-day period for each trap line. The small mammals were collected from four sites located in different regions of Georgia: (1) Dedoplistskaro in the Kakheti region of eastern Georgia close to Azerbaijan; (2) Aspindza in Samtskhe-Javakheti region of southern Georgia; (3) Marneuli in the Kvemo Kartli region in south-eastern Georgia; and (4) near Batumi City located on the Black Sea coast in Ajaria of southwestern Georgia.
Prior to blood and ectoparasite collection, the captured small mammals were anesthetized using open-drop exposure to Isoflurane (Aesica, Queenborough, UK) using a 20% volume/volume mixture for mice and voles and a 30% volume/volume mixture for rats and jirds. Blood samples from mice, voles, and small jirds were collected by bleeding from the retro-orbital plexus and by femoral venipunctures from medium-size animals (large jirds and rats). Rodents were then euthanized via either exsanguination or exposure to a lethal dose of Isoflurane. Tissue samples (blood, spleen, liver, and intestine) were sterilely collected from each rodent and placed in separate pre-labeled cryogenic vials.

Microbiology
The laboratory procedures for culturing Yersinia bacteria were conducted at the Richard Lugar Center for Public Health Research, Tbilisi, Georgia, following the standard operating procedures approved by NCDC's QMS committee (SOP-GBG -116-Identification of Yersinia species _V3.0, 2020). Suspensions of small intestine samples were inoculated into 1% peptone water followed by growth for 28 days at 4 • C. The inoculum was plated every 5th day on Endo and CIN agar (bioMérieux, Marcy-l'Étoile, France), and plates were incubated at 28 • C for 7 days. Bacterial colonies with a morphology resembling Yersinia (round shape, 0.5-1 mm diameter size, and grayish-pinkish color; a "bulls-eye" appearance on CIN agar) were sub-cultured on Endo agar plates at 28 • C for 24-48 h. The colonies obtained were analyzed microscopically after Gram staining. BD Phoenix (Becton & Dickinson, Eysins, Switzerland) and API (Analytical Profile Index, bioMérieux, France) tests were used according to the manufacturers' guidance for further analysis. Identification of Y. enterocolitica was conducted with the use of the API20E test system (Analytical Profile Index, bioMérieux, France) and the BD Phenix50 Automatic System using a Gram-Negative ID Panel (Becton & Dickinson, Switzerland) according to the manufacturers' guidance for further analysis.

Determination of Plasmid Content
The plasmid profile of Y. enterocolitica isolates was determined by the phenol method [25] as described by us previously [26]. The vaccine strain of Y. pestis, EV76, containing three plasmids of known molecular weight (101 kb, 70.5 kb, and 9.5 kb), was used as the plasmid size standard.

Selection of Y. enterocolitica Isolates for Sequencing
The isolates chosen for the WG sequencing are listed in Table 1. They were selected from different geographical locations across the country ( Figure 1) and were isolated from different rodent species and a shrew during the same year. Three of these strains were isolated from rodents collected in Dedoplistskaro, Kakheti Region, during the investigation of the active plague epizootic.

Genome Sequencing and Assembly
DNA for WG sequencing of six selected Y. enterocolitica strains were prepared using the QIAGEN Genomic DNA Kit (QIAGEN, Valencia, CA, USA). Extracted DNAs were quantified using a Qubit DNA HS Quantitation kit (Invitrogen, Thermo Fisher Scientific, Waltham, MA, USA) on a Qubit 3.0 fluorometer. Sequencing libraries were prepared using the NEBNext ® Ultra DNA library prep kit for Illumina (New England Biolabs, San Diego, CA, USA) according to the manufacturer's instructions. DNA (1 g) was sheared by sonication using a Covaris M220 sonication system (Covaris, Inc., Woburn, MA, USA) to an average of 500 bp fragment size. For library quantification, a Qubit dsDNA HS kit (Invitrogen) was used. Fragment size distribution was verified with a high sensitivity DNA chip on a Bioanalyzer (Agilent Technologies, Inc. Germany). Libraries were pooled and sequenced on an Illumina MiSeq instrument in paired-end mode 2 × 251 nt (500 cycles) using a Miseq reagent kit V2.

Genome Sequence Processing and Phylogenetic Analyses
Obtained raw sequence reads were processed using the CLC Genomics Workbench 12.0.3 (CLC Bio, Aarhus, Denmark). After quality trimming, the sequence reads were de novo assembled (k-mer size = 20). Contigs generated in CLC were uploaded into the EDGE bioinformatics platform version 2.3.1 and processed using the following analytical workflows: reference-based analysis, annotation (Prokka), taxonomy classification, and gene family analysis. For comparison against known references, RefSeq complete genomes of Y. entercolitica were selected from the EDGE drop-down list. The virulence genes were identified using ShortBRED with a database generated by the developers of EDGE using VFDB.
Illumina sequence reads were also uploaded into the integrated software environment EnteroBase (http://enterobase.warnick.ac.uk) under the Yersinia database. Sequence reads were processed with the EnteroBase pipeline, including the sequence read quality check, assembly (SPAdes), and annotation. The six uploaded Y. enterocolitica isolates were compared to a publicly available whole-genome sequence of selected 125 Y. enterocolitica strains in the database using an in silico SNP-typing scheme module. Maximum likelihood phylogenetic trees of core, non-repetitive SNPs were called against a reference draft genome (strain 1127) and constructed using the inbuilt dendrogram option of SNP Projects.

Prevalence of Y. entercolitica in Rodents and Shrews Based on Archived Records
During 1981−1997, Y. enterocolitica culture-positive rodents and shrews were reported during 14 of 17 years of the investigations. In total, Y. enterocolitica bacteria were isolated from 2052 rodents of 14 species (Table 2). Of those, most isolates were obtained from common voles Microtus arvalis (n = 1386), Norway rats Rattus norvegicus (n = 198), wood mice Apodemus spp. (n = 180), and house mice Mus musculus (n = 131) ( Table 2). Infected rodents and shrews were identified in different territories across Georgia, but rodent species carrying Y. enterocolitica varied. Infected commensal rodents, such as R. rattus, R. norvegicus, and M. musculus, were found in Tbilisi (capital of Georgia) and in seaports on the Black Sea coast (Batumi, Poti, and Sukhumi). Most infected common voles were found in southern and southwestern regions of Georgia, while infected social voles (Microtus socialis) and Libyan jirds (Meriones libycus) were captured in the steppe foothills of Eastern Georgia (Figure 2). Interestingly, there was a report of isolation of Y. enterocolitica from an Oriental rat flea Xenopsilla cheopis collected from R. norvegicus from the city of Tbilisi in 1987.

Isolation of Y. enterocolitica and Y. enterocolitica-Like Organisms during the Perspective Study
Growth of bacterial colonies morphologically suspected as Yersinia was observed after five days of incubation of 53 out of 572 freshly collected intestines from rodents and shrews randomly selected for this study. The characteristics of the colonies on CIN agar included a deep red center with transparent margins (a "bulls-eye" appearance). Gramstained microscopy showed gram-negative bacillar-shaped microorganisms. All Yersinialike isolates obtained from rodents were identified as Y. enterocolitica or Y. enterocolitica-like bacteria based on the API20E tests. Positive samples were obtained from R. norvegicus (7/63), M. musculus (10/36), Macedonian mice, M. macedonicus (4/85), Microtus socialis (5/39), M. arvalis (2/10), Apodemus spp. (21/187), bicolored shrews Crocidura leucodon (4/41), and from an unidentified, white-toothed shrew of the genus Crocidura (presumably, a lesser shrew C. suaveolens based on the habitat characteristics, but a complete identification cannot be confirmed).

Phenotypical Characterization of Freshly Isolated Y. enterocolitica Strains
Gram-stained microscopy showed gram-negative bacillar-shaped microorganisms. A bacterial suspension was used to rehydrate each of the wells, and the strips were incubated. All positive and negative API test results were compiled to obtain a profile number, which was then compared with profile numbers in a commercial codebook. All 53 Yersinia-like isolates were identified as Y. enterocolitica based on the API 20E tests and the BD Phenix50 Automatic System. The phenotypic API20E characteristics of 53 Y. enterocolitica isolates are shown in Supplementary Table S1.
Out of six isolates selected for sequencing, three had identical biochemical properties, strains 18ENCDC AsAtc Rint 11-12, 18E17 065 D I Rint-91, and 18E17 065 D I Rint-95. The latter two were isolated from the plague focus in Dedoplistskaro from different mouse species (Table 1). Another strain from this region isolated from a vole, 18E17 065 D I-33, differed in one position. Two other isolates, LcENCDC Rint 5 and B18 Rint-231 differed in two non-identical positions.  Table 3.
Both present and absent virulence functions pointed out that these isolates belong to the potential pathogenicity Y. enterocolitica biogroup A1. Nevertheless, the isolates differed from each other in a number of virulence proteins. Notably, only three strains contained classical invasin Inv, typically located within the flagellar cluster in Y. enterocolitica. Nevertheless, all strains, including those missing the Inv function, had putative invasin-like proteins with low homology to Inv. Three strains isolated in the same region (Dedoplistskaro) did not contain MyfA fimbrial protein, which is an adhesin essential for virulence and in vivo-expressed subtilisin/kexin-like protease HreP, RtxA-like putative leukotoxin, and cytolethal distending toxin subunit B-like protein CdtB. The strains isolated in other regions all contained these four virulence factors. Moreover, the Dedoplistskaro isolates contained tRNA(fMet)-specific endonuclease VapC, which was absent in the genomes of all other strains.
Two strains from Dedoplistskaro, 18E17 065 D I Rint-91 and 18E17 065 D I Rint-95, isolated in close proximity, although from different animal species (Table 1, Figure 1), were identical with regards to the possession of virulence genes. In contrast, the third strain from this region (18E17 065 D I-33) isolated distantly from the former two strains lacked ADP-rybosyltransferase Pertussis-like toxin YtxA, the function identified in all our strains. Finally, this strain contained a distinct virulence protein, YltA, heat-labile enterotoxin, a chain with significant homology to LTA toxin of enterotoxigenic E. coli and other microorganisms, such as Providencia alcalifaciens, Cronobacter, and Citrobacter sp. S-77 (Supplementary Figure S1). The unique feature of 18ENCDC AsAtc Rint-11-12 isolate was an absence of VapBC type II toxin-antitoxin system, and the B18 Rint-231 isolate lacked critical components of chromosomal T3SS; both functions were found in all other strains.

Plasmids
The screening of Y. enterocolitica isolates, from small mammals, by gel electrophoresis revealed the presence of plasmids in all three isolates from Dedoplistskaro (Supplementary Figure S2), while other isolates had no plasmids. Both strains 18E17 065 D I Rint-91 and 18E17 065 D I Rint-95 that had identical sets of virulence genes and phenotypic properties by API20E showed the same plasmid profile. We were not able to identify plasmid contigs in the sequenced genomes since these plasmids were large and probably had low copy numbers. The crude estimate of their molecular weight, based on the standards from Y. pestis vaccine strain EV76, suggested that the sizes of the large plasmids were about 200 kb, 90 kb, and 60 kb.
Another Dedoplistskaro isolate, 18E17 065 D I -33, likely had similar 200 kb and 60 kb plasmids and lacked the 90 kb plasmid. Nevertheless, we were able to identify contigs associated with the small multi-copy number plasmids in 18E17 065 D I Rint-91 and 18E17 065 D I Rint-95 isolates. The plasmids from both strains were identical except for a small deletion of 29 bp found in 18E17 065 D I Rint-91. The size of the plasmid was 4.5 kb, and it contained replication initiation protein repE_1, two short hypothetical proteins, and to our surprise, the type 4 secretion system (T4SS) VirB6 conjugal transfer protein. This protein has significant homology with other VirB6 proteins found in Yersinia pseudotuberculosis, Yersinia kristensenii, Yersinia massiliensis, Klebsiella pneumoniae NCTC9632, E. coli, and Salmonella enterica PNUSAS047895 (Supplementary Figure S3).

Phylogenetic Relationship
To define genetic relationships among the small mammal Y. enterocolitica isolates from Georgia and worldwide isolates of Y. enterocolitica strains, 125 core genome SNP-typing was performed using the EnteroBase platform. We included the isolates of Y. enterocolitica obtained from bats in Georgia [26] in this analysis. The results of the SNP-typing are shown in Supplementary Figure S4. As expected from the genomic analysis of the presence of virulence genes, all the isolates obtained from Georgian small mammals belonged to the biogroup 1A.
The three strains isolated from Dedoplistskaro grouped together and formed a separate cluster, which was quite distant from the rest of Y. enterocolitica isolates in the entire database. The strains 18E17 065 D I Rint-91 and 18E17 065 D I Rint-95 had identical SNP-types, despite being isolated from rodents captured in close proximity, though from different rodent species (M. macedonicus and Apodemus sp.). This suggests a potential transmission of this variant between mice of different species. Interestingly, the infected mice were identified within a plague-endemic territory where Y. pestis was recently detected in wild rodents and their fleas (unpublished data). The third strain (18E17 065 D I-33) from the same territory (Dedoplistskaro) was isolated from a vole M. socialis and was phylogenetically distinct from the above-mentioned two strains. There was only one strain in the entire EnteroBase, which belonged to the same cluster as the Dedoplistskaro strains, the strain Y18EE025 of an undetermined biogroup and serotype. Interestingly, the SNP-type of the strain 18ENCDC AsAtc Rint-11-12 isolated from R. norvegicus was identical to the Y. enterocolitica strains obtained from dead bats described in Georgia [26].
This finding suggests a wide distribution of this genotype of Y. enterocolitica in this part of Georgia. The SNP-type of the strain (B18 Rint -231) isolated from a shrew (Adilia, suburb of Batumi, Black Sea port) was also unique in the entire database. The only strain grouped together on tree 18005024 of biogroup 1Awas an undetermined human isolate serotype from Luxemburg, 2018. Finally, the LcENCDC Rint 5 strain isolated from R. norvegicus in another region in Southern Georgia also had a distinct SNP profile and grouped together with the strain YE13/03 of biogroup 1A, serotype O:6,30 isolated in the United Kingdom from human feces in 2003 (BioSample: SAMEA980154; SRA: ERS008611).

Discussions
Until recently, frequent reports of isolation of strains, which were identified as Y. enterocolitica from wild rodents collected in plague endemic areas in the Caucasus and Central Asia, remained undocumented in English language literature. The Soviet Union's Anti-plague system covered a vast territory of these regions. Zoologists and bacteriologists working for the Anti-plague stations captured and processed an enormous number of rodents in search of plague and other rodent-borne diseases. According to the laboratory procedure for the identification of plague pathogens, bacteriophage typing was routinely applied for bacteria morphologically identified as Yersinia. Bacteria of Y. pestis are sensitive to specific bacteriophage lysis at 20 • C, whereas bacteria of Y. enterocolitica do not lyse at temperatures below 28 • C [27].
Isolation of numerous cultures of Y. enterocolitica from various species of wild mammals in Georgia strongly suggests a wide circulation of this species of Yersinia in natural habitats across all landscape zones of the Caucasian region. This discovery was unexpected because previously, Y. enterocolitica has been primarily known as a human enteric pathogen with some involvement of domestic animals (pigs) with rare reports from commensal rats living closely to pig farms. Therefore, it was not entirely surprising that several specimens of R. rattus, and house mice M. musculus, were found culture-positive for Y. enterocolitica. Unexpected was a high occurrence of this bacterium in wild rodents inhabiting various natural habitats across the region.
During the perspective study organized specifically to identify Yersinia species in 14 rodents of different species, two species of voles (M. socialis and M. majori) and two species of wild mice (A. uralensis and A. witherbyi) were found to be infected with this bacterium. These rodents commonly occupy natural habitats far from human settlements and pig farms. Noticeably, four Caucasus field mice (A. ponticus), which were caught in a forest near a small pig farm where pigs reportedly grazed under the forest cover, were tested negative for Y. enterocolitica.
Another important question has arisen about the interaction between bacteria of two species belonging to the genus Yersinia (Y. pestis and Y. enterocolitica), which circulate in the same populations of rodents, which requires future investigations. Interestingly, Fukushima et al. (2001) proposed a role of Y. enterocolitica as a possible barrier against Y. pestis based on the investigation of rodents in a natural plague focus in Ningxia, China. The co-circulation of bacteria belonging to the genus of Yersinia in the same rodent populations within plague endemic territories of Georgia presents an excellent opportunity to address this question.
The results from this investigation and few recent publications [28,29] have challenged the dominant perception of Y. enterocolitica as a pathogen of people and domestic animals. Cultures obtained provide material for further genetic and bacteriological characterization of isolates of Y. enterocolitica and to compare them with human isolates. Such a comparative analysis can help to evaluate the risk for public health of the strains of Y. enterocolitica carried by wild rodents.
The phylogenetic analysis revealed that the small mammal Y. enterocolitica strains from Georgia belonged to the potentially pathogenic biogroup 1A [3,5,30]. The six isolates investigated here were quite heterogeneous with regard to the presence of different virulence genes, except two mouse strains 18E17 065 D I Rint-95 and 18E17 065 D I Rint-91 from the Dedoplistskaro region, which also contained identical plasmid and biochemical profiles (Supplementary Table S1 and Figure S2). These two strains were isolated in close proximity from two different mouse species, A. witherbyi, and M. macedonicus, respectively, and were generally indistinguishable by WG sequencing. The vole (M. socialis) strain from this area 18E17 065 D I-33 had an almost identical set of virulence determinants, only missing YtxA Pertussis-like toxin and two of four plasmids characteristic to the former strains.
The interesting feature of all three Dedoplistskaro isolates was a lack of fimbrial adhesin MyfA, suggesting a reduced intestinal colonization potential [5]. Nevertheless, the vole strain 18E17 065 D I-33 contained a gene encoding the A chain of heat-labile enterotoxin (LTA), a function that was not previously seen in Y. enterocolitica and rarely found in pathogens other than enterotoxigenic E. coli (Supplementary Figure S1). Since LTA was shown to enhance enteric pathogen adherence and intestinal colonization [31], this may compensate for the lack of the MyfA adhesin in 18E17 065 D I-33 isolate. However, a hallmark of the heat-labile enterotoxin in Escherichia coli is its remarkable ability to contribute to diarrhea in travelers, children, and young animals [32]. It remains to be seen whether this Y. enterocolitica isolate can produce a holotoxin, where the LTA is assembled with the receptor binding B-subunits (LTB) to promote diarrhea in animal models. The 18E17 065 D I-33 genome search did not reveal obvious candidates for the gene encoding the LTB, which may have a low homology with the LTB from E. coli and a distinct receptor binding specificity.
Overall, all six small mammal bacteria strains (Y. enterocolitica) from Georgia, as other typical biogroup A1 strains, were deficient in essential virulence factors, such as heatstable enterotoxin A YstA, adhesin Ail, insecticidal toxin TccC, T3SS from the virulence plasmid pYV, and virulence regulator VirF (Supplementary Table S1). Nevertheless, all strains carried the genes for heat-stable enterotoxin B YstB, a modulator of expression for virulence functions YmoA, ferric enterobactin transport system FepBDGC, direct heme uptake system HemPR-HmuVSTU, pore-forming cytotoxin YaxAB, phospholipases PhlA and PldA, arcenic cluster ArsCBR, flagellar apparatus, and a type two secretion system. The isolates varied in the presence of invasin Inv, adhesin MyfA, subtilisin/kexin-like protease HreP, RtxA-like putative leukotoxin, endonuclease VapC, Pertussis-like toxin YtxA and cytolethal distending toxin subunit B-like protein CdtB. The differences in the virulence gene composition suggested a marked diversity of small mammal Y. enterocolitica infections in Georgia, although clonal domination can be possible as indicated by the similarity of the strains from Dedoplistskaro, which were isolated from different rodent species.
The diversity of small mammal Y. enterocolitica infections in Georgia is clearly seen on the whole-genome SNP-based phylogenetic tree (Supplementary Figure S4). The strains belonged to four distant clusters with little or no other Y. enterocolitica strains from the Enterobase present in such clusters. The exception was the strain LcENCDC Rint 5 isolated from R. norvegicus in South Georgia, which grouped with several strains of biogroup 1A. This Georgian isolate was phylogenetically close to the strain YE13/03, serotype O:6,30 isolated in the United Kingdom from human feces in 2003. Interestingly, the strain 18ENCDC AsAtc Rint-11-12 isolated from Rattus norvegicus formed a distinct and unique cluster with Y. enterocolitica isolates from dead bats in Georgia [26]. Despite the distance between the two sites of isolation being only about 100 kilometers (Figure 1), the two places are separated by mountains, likely preventing a direct communication of the rat and bat populations between the areas. Nevertheless, this might be another example of a wide clonal distribution of Y. enterocolitica variants in certain parts of Georgia.
Finally, the most intriguing finding was the observation that Y. enterocolitica isolated from different animal species in the natural plague focus of the Dedoplistskaro region were grouped together and might form a different clade or even biogroup in the future. The cluster containing these isolates did not possess Y. enterocolitica of a known biogroup, although the set of virulence genes suggests that these isolates are closely related to biogroup 1A. The role of ongoing outbreaks of Y. pestis in this area for the selection of these unique Y. enterocolitica variants will need further investigation. Overall, the reported observations contribute to a better understanding of the potential public health importance of Yersinia species of small mammals beyond plague [10,18,26,32].
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/pathogens10080939/s1, Figure S1: Alignment of heat-labile enterotoxin, A chain; Figure S2: Plasmid profiles of the strains obtained by Kado & Liu method; Figure S3: Alignment of ORF1 of small plasmid; Figure S4: Maximum likelihood-phylogenetic SNP-based tree constructed in Enterobase showing the position of sequenced Y. enterocolitica isolates from Georgia in the context of worldwide Y. enterocolitica biogroup distribution, Table S1: Yersinia enterocolitica strains phenotypic characterization by API20E; Table S2: Presence of the classical virulence-associated genes of Yersinia enterocolitica in the strains from Georgia.