Phylogeographic separation and formation of sexually discrete lineages in a global population of Yersinia pseudotuberculosis

Genome Analysis of Japanese Yersinia pseudotuberculosis Strains Isolated From Kawasaki Disease Patients and Other Sources and Their Phylogenetic Positions in the Global Y. pseudotuberculosis Population

Yersinia pseudotuberculosis (Ypt) is a gram-negative bacterium that infects both humans and animals primarily through fecal‒oral transmission. While Ypt causes acute gastroenteritis in humans, an association with Kawasaki disease (KD), a disease that primarily affects infants and young children and causes multisystemic vasculitis, has also been suspected. Although KD represents a significant health concern worldwide, the highest annual incidence rate is reported in Japan. Previously, a geographical origin-dependent population structure of Ypt comprising the Asian, transitional, and European clades was proposed. However, genomic data on KD-associated Ypt strains is currently unavailable. In this study, to analyze the phylogenetic and genomic features of KD-associated strains, we determined the whole-genome sequences of 35 Japanese Ypt strains, including 11 KD-associated strains, and constructed a genome set (n = 204) representing the global population of Ypt by adding publicly available Ypt genomes. In a phylogenetic analysis, all sequenced Japanese strains, including the KD-associated strains, belonged to the Asian clade, which appeared to be the ancestral clade of Ypt, and the KD-associated strains belonged to multiple lineages in this clade. Strains from patients with Far East scarlet-like fever (FESLF), a KD-related disease, also belonged to the Asian clade. Moreover, no KD strain-specific genes were identified in pan-genome-wide association study analyses. Notably, however, the gene encoding a superantigen called Yersinia pseudotuberculosis-derived mitogen (YPM) showed a distribution pattern highly biased to the Asian clade. Although further studies are needed, our results suggest that Asian clade strains may have a greater potential to trigger KD.

INDEL-Typing of <i>Yersinia pseudotuberculosis</i> Strains

The aim of this study was to develop a new method of intraspecific genetic differentiation of Yersinia pseudotuberculosis, based on the detection of INDEL-markers using PCR. Materials and methods. Analyzed were 308 strains from the NCBI database and 15 strains sequenced within the frames of this study. The nucleotide sequences of the strains were determined using the MiSeq technology platform. The genomes of the strains sequenced in the work, as well as genomes from the NCBI database, were assessed using in silico PCR with 7 pairs of primers designed in the study. As a result of a comparison of genome-wide sequences of 22 Y. pseudotuberculosis strains from the NCBI database, using the author’s software (GenExpert), 7 INDEL-markers were selected that make it possible to effectively distinguish between strains of the causative agent of pseudotuberculosis. Based on these markers, 7 pairs of primers were designed and synthesized for the analysis of different strains using PCR. Analysis of 323 strains in PCR in silico and 70 strains in PCR in vitro allowed for dividing them into 30 genetic groups. Comparison of the results of PCR in silico and in vitro confirmed the possibility of using the proposed primers for intraspecific differentiation of Y. pseudotuberculosis. Based on the data obtained, a dendrogram reflecting the phylogenetic relations of different strains of Y. pseudotuberculosis was constructed. When analyzing the distribution of Y. pseudotuberculosis strains by various clusters and genetic groups, a number of patterns were revealed. Conducted in silico and in vitro PCR show that the proposed method of INDEL-typing can be used for intraspecific genetic differentiation of the causative agent of pseudotuberculosis.

Yersiniosis in New Zealand

The rate of yersiniosis in New Zealand (NZ) is high compared with other developed countries, and rates have been increasing over recent years. Typically, >99% of human cases in NZ are attributed to Yersinia enterocolitica (YE), although in 2014, a large outbreak of 220 cases was caused by Yersinia pseudotuberculosis. Up until 2012, the most common NZ strain was YE biotype 4. The emergent strain since this time is YE biotype 2/3 serotype O:9. The pathogenic potential of some YE biotypes remains unclear. Most human cases of yersiniosis are considered sporadic without an identifiable source. Key restrictions in previous investigations included insufficient sensitivity for the isolation of Yersinia spp. from foods, although foodborne transmission is the most likely route of infection. In NZ, YE has been isolated from a variety of sick and healthy domestic and farm animals but the pathways from zoonotic reservoir to human remain unproven. Whole-genome sequencing provides unprecedented discriminatory power for typing Yersinia and is now being applied to NZ epidemiological investigations. A "One-Health" approach is necessary to elucidate the routes of transmission of Yersinia and consequently inform targeted interventions for the prevention and management of yersiniosis in NZ.

Genomic Epidemiology and Phenotyping Reveal on-Farm Persistence and Cold Adaptation of Raw Milk Outbreak-Associated Yersinia pseudotuberculosis

Packaged raw milk contaminated with Yersinia pseudotuberculosis mediated a large yersiniosis outbreak in southern Finland in 2014. The outbreak was traced back to a single dairy farm in southern Finland. Here we explore risk factors leading to the outbreak through epidemiologic investigation of the outbreak farm and through genomic and phenotypic characterization of the farm’s outbreak and non-outbreak associated Y. pseudotuberculosis strains. We show that the outbreak strain persisted on the farm throughout the 7-month study, whereas the non-outbreak strains occurred sporadically. Phylogenomic analysis illustrated that the outbreak strain was related to previously published genomes of wild animal isolates from Finland, implying that wild animals were a potential source of the outbreak strain to the farm. We observed allelic differences between the farm’s outbreak and non-outbreak strains in several genes associated with virulence, stress response and biofilm formation, and found that the outbreak strain formed biofilm in vitro and maintained better growth fitness during cold stress than the non-outbreak strains. Finally, we demonstrate the rapid growth of the outbreak strain in packaged raw milk during refrigerated storage. This study provides insight of the risk factors leading to the Y. pseudotuberculosis outbreak, highlights the importance of pest control to avoid the spread of pathogens from wild to domestic animals, and demonstrates that the cold chain is insufficient as the sole risk management strategy to control Y. pseudotuberculosis risk associated with raw drinking milk.

COMPARATIVE ANALYSIS OF CRISPR-CAS SYSTEM STRUCTURES OF YERSINIA PSEUDOTUBERCULOSIS IP32953 AND IP31758

Background. Pseudotuberculosis is still relevant problem in medical science and public health of Russia and other countries. Typing of Y. рseudotuberculosis strains by their CRISPR systems is a perspective tool for monitoring of Yersinia populations as was shown in Y. pestis.

Aims. Here we describe and compare CRISPR-Cas systems of Yersinia pseudotuberculosis strains IP32953 and IP31758 causing classic pseudotuberculosis and Far-East scarlet-like fever (FESLF) respectively.

Materials and methods. Complete genomes of Y. pseudotuberculosis IP329353 and IP31758 (NC_006155 and NC_009708 respectively) were obtained from NCBI Nucleotide Database. Search; identification; and analysis of CRISPR systems were carried out by online-tools CRISPROne; CRISPRDetect; and CRISPRTarget.

Results and discussion. Analyzed strains have CRISPR-Cas systems that include one set of cas-genes and arrays situated at the long distances from each other. We defined three CRISPR arrays in Y. pseudotuberculosis IP32953 by the combination of program methods. CRISPR-Cas system of this strain consist of array YP1 located near cas-genes; arrays YP2 and YP3. CRISPR-Cas system of Y. pseudotuberculosis IP31758 includes two arrays – YP1 and YP3. CRISPR systems do not share similar spacers. CRISPR systems of the analyzed strains differ in CRISPR loci and cas-protein structures that can be used as specific marks of analyzed strains.

Conclusions. We suggest that acquisition of certain spacers may play a role in evolution and divergence of Y. pseudotuberculosis strains.