Yersinia pestis genome sequencing identifies patterns of global phylogenetic diversity

Integrative species delimitation reveals an Idaho-endemic ground squirrel, Urocitellus idahoensis (Merriam 1913)

The "small-eared" species group of Urocitellus ground squirrels (Sciuridae: Xerinae: Marmotini) is endemic to the Great Basin, United States, and surrounding cold desert ecosystems. Most specific and subspecific lineages in this group occupy narrow geographic ranges, and some are of significant conservation concern; despite this, current taxonomy remains largely based on karyotypic or subtle pelage and morphological characteristics. Here, we leverage 2 multilocus DNA sequence data sets and apply formal species delimitation tests alongside morphometric comparisons to demonstrate that the most widespread small-eared species (U. mollis Kennicott, 1863 sensu lato; Piute Ground Squirrel) is comprised of 2 nonsister and deeply divergent lineages. The 2 lineages are geographically separated by the east-west flowing Snake River in southern Idaho, with no sites of sympatry currently known. Based on robust support across the nuclear genome, we elevate populations previously attributed to U. mollis from north of the Snake River to species status under the name Urocitellus idahoensis (Merriam 1913) and propose the common name "Snake River Plains Ground Squirrel" for this taxon. We delimit 2 subspecies within U. idahoensis; U. i. idahoensis (Merriam 1913) in western Idaho and U. i. artemesiae (Merriam 1913) in eastern Idaho. Urocitellus idahoensis is endemic to Idaho and has a maximal range area of roughly 29,700 km2 spanning 22 counties but occurs discontinuously across this area. Our work substantially expands knowledge of ground squirrel diversity in the northern Great Basin and Columbia Plateau and highlights the difficulty in delimiting aridland mammals whose morphological attributes are highly conserved.

Complete genome sequences of Yersinia pestis 6/69 strain isolated from a bubonic plague patient in Madagascar and its isogenic strain cured of pPCP1

We report the complete genome sequences of two valuable strains to investigate plague pathogenesis: (i) Yersinia pestis strain 6/69, which was isolated from a bubonic plague patient in Madagascar and contains pCD1, pMT1, and pPCP1 virulence plasmids, and (ii) the 6/69 strain cured of pPCP1.

Elucidation of peptide screen for targeted identification of Yersinia pestis by nano-liquid chromatography tandem mass spectrometry

Yersinia pestis, a Gram-negative bacterium is the causative agent of the fatal communicable disease plague. The disease had a profound impact on human history. Plague bacteria are usually transmitted to humans through the bite of an infected rat flea. Earlier studies have indicated that Y. pestis can survive in environmental matrices e.g. water and soil. This study aimed to generate a peptide-based screen for identification of Y. pestis particularly from environmental matrices. We employed a shotgun proteomic approach using nano-liquid chromatography-tandem mass spectrometry (nLC-MS/MS) to discover Y. pestis-specific peptides. The pure cultures of Y. pestis and related species were grown, their proteome were delineated and analyzed by in silico tools to discover 61 Y. pestis specific peptides. Additionally, 148 peptides were discovered from proteins of Y. pestis-specific plasmids and chromosomal-associated virulence markers. To validate this screen of 209 peptides, various concentrations of Y. pestis (ranging from 1.3 × 108 to 1.3 × 105 cfu) were spiked into garden soil. Y. pestis could be identified in all samples except un-spiked negative control soil sample. This study offers a valuable method for the identification of Y. pestis, by tandem mass spectrometry which may be used in environmental and clinical matrices.

Hotspots of genetic change in Yersinia pestis

The relative contributions of mutation rate variation, selection, and recombination in shaping genomic variation in bacterial populations remain poorly understood. Here we analyze 3318 Yersinia pestis genomes, spanning nearly a century and including 2336 newly sequenced strains, to shed light on the patterns of genetic diversity and variation distribution at the population level. We identify 45 genomic regions ("hot regions", HRs) that, although comprising a minor fraction of the genome, are hotbeds of genetic variation. These HRs are distributed non-randomly across Y. pestis phylogenetic lineages and are primarily linked to regulatory genes, underscoring their potential functional significance. We explore various factors contributing to the shaping and maintenance of HRs, including genomic context, homologous recombination, mutation rate variation and natural selection. Our findings suggest that positive selection is likely the primary driver behind the emergence of HRs, but not the sole force, as evidenced by the pronounced trend of variation purging within these regions.

Review of genotyping methods for Yersinia pestis in Madagascar

BACKGROUND

Plague, a zoonotic disease caused by Yersinia pestis, was responsible for 3 historical human pandemics that killed millions of people. It remains endemic in rodent populations in Africa, Asia, North America, and South America but human plague is rare in most of these locations. However, human plague is still highly prevalent in Madagascar, which typically records a significant part of all annual global cases. This has afforded an opportunity to study contemporary human plague in detail using various typing methods for Y. pestis.

AIM

This review aims to summarize the methods that have been used to type Y. pestis in Madagascar along with the major discoveries that have been made using these approaches.

METHODS

Pubmed and Google Scholar were used to search for the keywords: "typing Yersinia pestis Madagascar," "evolution Yersinia pestis Madagascar," and "diversity Yersinia pestis Madagascar." Eleven publications were relevant to our topic and further information was retrieved from references cited in those publications.

RESULTS

The history of Y. pestis typing in Madagascar can be divided in 2 periods: the pre-genomics and genomics eras. During the pre-genomics era, ribotyping, direct observation of plasmid content and plasmid restriction fragment length polymorphisms (RFLP) were employed but only revealed a limited amount of diversity among Malagasy Y. pestis strains. Extensive diversity only started to be revealed in the genomics era with the use of clustered regularly interspaced palindromic repeats (CRISPR), multiple-locus variable number tandem repeats (VNTR) analysis (MLVA), and single-nucleotide polymorphisms (SNPs) discovered from whole genome sequences. These higher-resolution genotyping methods have made it possible to highlight the distribution and persistence of genotypes in the different plague foci of Madagascar (Mahajanga and the Central and Northern Highlands) by genotyping strains from the same locations across years, to detect transfers between foci, to date the emergence of genotypes, and even to document the transmission of antimicrobial resistant (AMR) strains during a pneumonic plague outbreak. Despite these discoveries, there still remain topics that deserve to be explored, such as the contribution of horizontal gene transfer to the evolution of Malagasy Y. pestis strains and the evolutionary history of Y. pestis in Madagascar.

CONCLUSIONS

Genotyping of Y. pestis has yielded important insights on plague in Madagascar, particularly since the advent of whole-genome sequencing (WGS). These include a better understanding of plague persistence in the environment, antimicrobial AMR and multi-drug resistance in Y. pestis, and the person-to-person spread of pneumonic plague. Considering that human plague is still a significant public health threat in Madagascar, these insights can be useful for controlling and preventing human plague in Madagascar and elsewhere, and also are relevant for understanding the historical pandemics and the possible use of Y. pestis as a biological weapon.

Phylogenetic Relationships and Evolution of the Genus Eganvirus (186-Type) Yersinia pestis Bacteriophages

Plague is an endemic infectious disease caused by Yersinia pestis. In this study, we isolated fourteen phages with similar sequence arrangements to phage 186; these phages exhibited different lytic abilities in Enterobacteriaceae strains. To illustrate the phylogenetic relationships and evolutionary relationships between previously designated 186-type phages, we analysed the complete sequences and important genes of the phages, including whole-genome average nucleotide identity (ANI) and collinearity comparison, evolutionary analysis of four conserved structural genes (V, T, R, and Q genes), and analysis of the regulatory genes (cI, apl, and cII) and integrase gene (int). Phylogenetic analysis revealed that thirteen of the newly isolated phages belong to the genus Eganvirus and one belongs to the genus Felsduovirus in the family Peduoviridae, and these Eganvirus phages can be roughly clustered into three subgroups. The topological relationships exhibited by the whole-genome and structural genes seemed similar and stable, while the regulatory genes presented different topological relationships with the structural genes, and these results indicated that there was some homologous recombination in the regulatory genes. These newly isolated 186-type phages were mostly isolated from dogs, suggesting that the resistance of Canidae to Y. pestis infection may be related to the wide distribution of phages with lytic capability.

Multiple Introductions of Yersinia pestis during Urban Pneumonic Plague Epidemic, Madagascar, 2017

Pneumonic plague (PP) is characterized by high infection rate, person-to-person transmission, and rapid progression to severe disease. In 2017, a PP epidemic occurred in 2 Madagascar urban areas, Antananarivo and Toamasina. We used epidemiologic data and Yersinia pestis genomic characterization to determine the sources of this epidemic. Human plague emerged independently from environmental reservoirs in rural endemic foci >20 times during August-November 2017. Confirmed cases from 5 emergences, including 4 PP cases, were documented in urban areas. Epidemiologic and genetic analyses of cases associated with the first emergence event to reach urban areas confirmed that transmission started in August; spread to Antananarivo, Toamasina, and other locations; and persisted in Antananarivo until at least mid-November. Two other Y. pestis lineages may have caused persistent PP transmission chains in Antananarivo. Multiple Y. pestis lineages were independently introduced to urban areas from several rural foci via travel of infected persons during the epidemic.

The Natural and Clinical History of Plague: From the Ancient Pandemics to Modern Insights

The human pathogen Yersinia pestis is responsible for bubonic, septicemic, and pneumonic plague. A deeply comprehensive overview of its historical context, bacteriological characteristics, genomic analysis based on ancient DNA (aDNA) and modern strains, and its impact on historical and actual human populations, is explored. The results from multiple studies have been synthesized to investigate the origins of plague, its transmission, and effects on different populations. Additionally, molecular interactions of Y. pestis, from its evolutionary origins to its adaptation to flea-born transmission, and its impact on human and wild populations are considered. The characteristic combinations of aDNA patterns, which plays a decisive role in the reconstruction and analysis of ancient genomes, are reviewed. Bioinformatics is fundamental in identifying specific Y. pestis lineages, and automated pipelines are among the valuable tools in implementing such studies. Plague, which remains among human history's most lethal infectious diseases, but also other zoonotic diseases, requires the continuous investigation of plague topics. This can be achieved by improving molecular and genetic screening of animal populations, identifying ecological and social determinants of outbreaks, increasing interdisciplinary collaborations among scientists and public healthcare providers, and continued research into the characterization, diagnosis, and treatment of these diseases.

Population structure, dispersion patterns and genetic diversity of two major invasive and commensal zoonotic disease hosts (Rattus norvegicus and Rattus tanezumi) from the southeastern coast of China

Background: The invasive brownrat (Rattus norvegicus) and the Oriental rats (Rattus tanezumi) are common commensal murid that are important hosts for rodent-borne diseases in southeast Asia. Understanding their population structure and genetic diversity is essential to uncover their invasion biology and distribution dynamics that are essential for controlling rodent-borne diseases. Methods: TA total of 103 R. norvegicus and 85 R. tanezumi were collected from 13 to 9 coastal areas of six provincial monitoring sentinel sites, respectivelyto assess patterns in their microsatellite loci and their mitochondrial coxl gene region. Results: Eleven sampled populations of R. norvegicus were divided into two major clusters by region. The observed heterozygosity values of all regional populations were smaller than expected genetic diversity heterozygosity values and deviated from Hardy-Weinberg equilibrium Nine sample populations of R. tanezumi were divided into three clusters; two that included sample from Hainan and Fujian provinces, and one that included samples from the other provinces and cities. The genetic diversity of R. tanezumi was highest in samples from Jiangsu and Guangdong provinces. Conclusion: The data in this paper confirm the two invasive rodent species from the southeastern coastal region of China may have relied on maritime transport to spread from the southern region of China to the Yangtze River basin. R. tanezumi may then hanve migrated unidirectionally, along the southeastern provinces of China towards the north, while R. norvegicus spread in a complex and multidirectional manner in Hainan, Fujian, Zhejiang and Jiangsu Provinces of the country.

Comparative sequence analysis elucidates the evolutionary patterns of Yersinia pestis in New Mexico over thirty-two years

Background

Yersinia pestis, a Gram-negative bacterium, is the causative agent of plague. Y. pestis is a zoonotic pathogen that occasionally infects humans and became endemic in the western United States after spreading from California in 1899.

Methods

To better understand evolutionary patterns in Y. pestis from the southwestern United States, we sequenced and analyzed 22 novel genomes from New Mexico. Analytical methods included, assembly, multiple sequences alignment, phylogenetic tree reconstruction, genotype-phenotype correlation, and selection pressure.

Results

We identified four genes, including Yscp and locus tag YPO3944, which contained codons undergoing negative selection. We also observed 42 nucleotide sites displaying a statistically significant skew in the observed residue distribution based on the year of isolation. Overall, the three genes with the most statistically significant variations that associated with metadata for these isolates were sapA, fliC, and argD. Phylogenetic analyses point to a single introduction of Y. pestis into the United States with two subsequent, independent movements into New Mexico. Taken together, these analyses shed light on the evolutionary history of this pathogen in the southwestern US over a focused time range and confirm a single origin and introduction into North America.

Genomic diversity of Yersinia pestis from Yunnan Province, China, implies a potential common ancestor as the source of two plague epidemics

Plague, caused by Yersinia pestis, is a zoonotic disease that can reemerge and cause outbreaks following decades of latency in natural plague foci. However, the genetic diversity and spread pattern of Y. pestis during these epidemic-silent cycles remain unclear. In this study, we analyze 356 Y. pestis genomes isolated between 1952 and 2016 in the Yunnan Rattus tanezumi plague focus, China, covering two epidemic-silent cycles. Through high-resolution genomic epidemiological analysis, we find that 96% of Y. pestis genomes belong to phylogroup 1.ORI2 and are subdivided into two sister clades (Sublineage1 and Sublineage2) characterized by different temporal-spatial distributions and genetic diversity. Most of the Sublineage1 strains are isolated from the first epidemic-silent cycle, while Sublineage2 strains are predominantly from the second cycle and revealing a west to east spread. The two sister clades evolved in parallel from a common ancestor and independently lead to two separate epidemics, confirming that the pathogen responsible for the second epidemic following the silent interval is not a descendant of the causative strain of the first epidemic. Our results provide a mechanism for defining epidemic-silent cycles in natural plague foci, which is valuable in the prevention and control of future plague outbreaks.

Genomic epidemiological analysis of county-scale Yersinia pestis spread pattern over 50 years in a Southwest Chinese prefecture

Plague, one of the most devastating infectious diseases in human history, is caused by the bacterium Yersinia pestis. Since the 1950s, the Dehong Dai-Jingpo Autonomous Prefecture (DH) in Yunnan Province, China, has recorded plague outbreaks that have resulted in 1,153 human cases and 379 deaths. The genetic diversity and transmission characteristics of Y. pestis strains in this region remain unknown. Here, we performed high-resolution genomic epidemiological analysis of 175 Y. pestis strains isolated from five counties and 19 towns in DH between 1953 and 2007. Phylogenetic analysis revealed that most DH strains were located in lineage 1.ORI2, which could be further subdivided into seven sub-phylogroups (SPG1-SPG7). The dominant sub-phylogroups of Y. pestis in DH varied during different periods and presented a population shift. Genomic evidence showed that plague might have emerged from the southwest of DH (e.g., Longchuan or Ruili counties) or its bordering countries, and subsequently spread to the northeast in multiple waves between 1982 and 2007. Our study infers a fine-scale phylogeny and spread pattern of the DH Y. pestis population, which extends our knowledge regarding its genetic diversity and provides clues for the future prevention and control of plague in this region.

The complex genomic diversity of Yersinia pestis on the long-term plague foci in Qinghai-Tibet plateau

Plague is a typical natural focus disease that circulates in different ecology of vectors and reservoir hosts. We conducted genomic population and phylogenetic analyses of the Yersinia pestis collected from the 12 natural plague foci in China with more than 20 kinds of hosts and vectors. Different ecological landscapes with specific hosts, vectors, and habitat which shape various niches for Y. pestis. The phylogeographic diversity of Y. pestis in different kinds plague foci in China showed host niches adaptation. Most natural plague foci strains are region-and focus-specific, with one predominant subpopulation; but the isolates from the Qinghai-Tibet plateau harbor a higher genetic diversity than other foci. The Y. pestis from Marmota himalayana plague foci are defined as the ancestors of different populations at the root of the evolutionary tree, suggesting several different evolutionary paths to other foci. It has the largest pan-genome and widest SNP distances with most accessory genes enriched in mobilome functions (prophages, transposons). Geological barriers play an important role in the maintenance of local Y. pestis species and block the introduction of non-native strains. This study provides new insights into the control of plague outbreaks and epidemics, deepened the understanding of the evolutionary history of MHPF (M. himalayana plague focus) in China. The population structure and identify clades among different natural foci of China renewed the space cognition of the plague.

Knockdown resistance mutations are common and widely distributed in Xenopsylla cheopis fleas that transmit plague in Madagascar

BACKGROUND

Plague, caused by the bacterium Yersinia pestis, remains an important disease in Madagascar, where the oriental rat flea, Xenopsylla cheopis, is a primary vector. To control fleas, synthetic pyrethroids (SPs) have been used for >20 years, resulting in resistance in many X. cheopis populations. The most common mechanisms of SP resistance are target site mutations in the voltage-gated sodium channel (VGSC) gene.

METHODOLOGY/PRINCIPAL FINDINGS

We obtained 25 collections of X. cheopis from 22 locations across Madagascar and performed phenotypic tests to determine resistance to deltamethrin, permethrin, and/or dichlorodiphenyltrichloroethane (DDT). Most populations were resistant to all these insecticides. We sequenced a 535 bp segment of the VGSC gene and identified two different mutations encoding distinct substitutions at amino acid position 1014, which is associated with knockdown resistance (kdr) to SPs in insects. Kdr mutation L1014F occurred in all 25 collections; a rarer mutation, L1014H, was found in 12 collections. There was a significant positive relationship between the frequency of kdr alleles and the proportion of individuals surviving exposure to deltamethrin. Phylogenetic comparisons of 12 VGSC alleles in Madagascar suggested resistant alleles arose from susceptible lineages at least three times. Because genotype can reasonably predict resistance phenotype, we developed a TaqMan PCR assay for the rapid detection of kdr resistance alleles.

CONCLUSIONS/SIGNIFICANCE

Our study provides new insights into VGSC mutations in Malagasy populations of X. cheopis and is the first to report a positive correlation between VGSC genotypes and SP resistance phenotypes in fleas. Widespread occurrence of these two SP resistance mutations in X. cheopis populations in Madagascar reduces the viability of these insecticides for flea control. However, the TaqMan assay described here facilitates rapid detection of kdr mutations to inform when use of these insecticides is still warranted to reduce transmission of plague.

Parallelism in Speciation and Intraspecific Diversification of the Plague Microbe Yersinia pestis

Modern "molecular genetic (MG) phylogenies" of the plague microbe Yersinia pestis, built on models of neutral evolution using statistical methods of phylogenetic analysis, contradict numerous obvious environmental (ECO) patterns and are not consistent with the concept of adaptatiogenesis. The reason for the discrepancy between MG and ECO phylogenies is seen in the underestimation by the MG approach of parallelisms in the processes of speciation and intraspecific diversification of the plague microbe. ECO methods showed the parallel tritope (almost) simultaneous speciation of three primary genovariants (populations, subspecies) Y. pestis 2.ANT3, 3.ANT2, and 4.ANT1 in three geographical populations of the Mongolian marmot (Marmota sibirica), which in the MG approach is mistaken for polytomy ("Big Bang"), caused by unknown natural phenomena on the eve of the first pandemic (Justinian's plague, 6th-8th centuries AD). The discrepancy between the MG and ECO interpretations of the evolution of intraspecifically-derived phylogenetic subbranches 0.PE and 2.MED is also associated with parallel evolutionary processes in independent lines, based on genovariants 2.ANT3, 3.ANT2, and 4.ANT1. The independence of these phylogenetic lines and parallelisms of sub-branches 0.PE and 2.MED associated with them are not taken into account in the MG approach. The prospect of creating a real phylogenetic tree for Y. pestis depends on a creative synthesis of the two approaches-MG and ECO.

Phylogenetic analysis of the origin and spread of plague in Madagascar

BACKGROUND

Plague is a zoonotic disease caused by the bacterium Yersinia pestis, highly prevalent in the Central Highlands, a mountainous region in the center of Madagascar. After a plague-free period of over 60 years in the northwestern coast city of Mahajanga, the disease reappeared in 1991 and caused several outbreaks until 1999. Previous research indicates that the disease was reintroduced to the city of Mahajanga from the Central Highlands instead of reemerging from a local reservoir. However, it is not clear how many reintroductions occurred and when they took place.

METHODOLOGY/PRINCIPAL FINDINGS

In this study we applied a Bayesian phylogeographic model to detect and date migrations of Y. pestis between the two locations that could be linked to the re-emergence of plague in Mahajanga. Genome sequences of 300 Y. pestis strains sampled between 1964 and 2012 were analyzed. Four migrations from the Central Highlands to Mahajanga were detected. Two resulted in persistent transmission in humans, one was responsible for most of the human cases recorded between 1995 and 1999, while the other produced plague cases in 1991 and 1992. We dated the emergence of the Y. pestis sub-branch 1.ORI3, which is only present in Madagascar and Turkey, to the beginning of the 20th century, using a Bayesian molecular dating analysis. The split between 1.ORI3 and its ancestor lineage 1.ORI2 was dated to the second half of the 19th century.

CONCLUSIONS/SIGNIFICANCE

Our results indicate that two independent migrations from the Central Highlands caused the plague outbreaks in Mahajanga during the 1990s, with both introductions occurring during the early 1980s. They happened over a decade before the detection of human cases, thus the pathogen likely survived in wild reservoirs until the spillover to humans was possible. This study demonstrates the value of Bayesian phylogenetics in elucidating the re-emergence of infectious diseases.

Yersinia pestis Δail Mutants Are Not Susceptible to Human Complement Bactericidal Activity in the Flea

Ail confers serum resistance in humans and is a critical virulence factor of Y. pestis, the causative agent of plague. Here, the contribution of Ail for Y. pestis survival in the flea vector was examined. Rat or human but not mouse sera were bactericidal against a Y. pestis Δail mutant at 28°C in vitro. Complement components deposited rapidly on the Y. pestis surface as measured by immunofluorescent microscopy. Ail reduced the amount of active C3b on the Y. pestis surface. Human sera retained bactericidal activity against a Y. pestis Δail mutant in the presence of mouse sera. However, in the flea vector, the serum protective properties of Ail were not required. Flea colonization studies using murine sera and Y. pestis KIM6+ wild type, a Δail mutant, and the Δail/ail+ control showed no differences in bacterial prevalence or numbers during the early stage of flea colonization. Similarly, flea studies with human blood showed Ail was not required for serum resistance. Finally, a variant of Ail (AilF100V E108_S109insS) from a human serum-sensitive Y. pestis subsp. microtus bv. Caucasica 1146 conferred resistance to human complement when expressed in the Y. pestis KIM6+ Δail mutant. This indicated that Ail activity was somehow blocked, most likely by lipooligosaccharide, in this serum sensitive strain. IMPORTANCE This work contributes to our understanding of how highly virulent Y. pestis evolved from its innocuous enteric predecessor. Among identified virulence factors is the attachment invasion locus protein, Ail, that is required to protect Y. pestis from serum complement in all mammals tested except mice. Murine sera is not bactericidal. In this study, we asked, is bactericidal sera from humans active in Y. pestis colonized fleas? We found it was not. The importance of this observation is that it identifies a protective niche for the growth of serum sensitive and nonsensitive Y. pestis strains.

Yersiniomics, a Multi-Omics Interactive Database for Yersinia Species

The genus Yersinia includes a large variety of nonpathogenic and life-threatening pathogenic bacteria, which cause a broad spectrum of diseases in humans and animals, such as plague, enteritis, Far East scarlet-like fever (FESLF), and enteric redmouth disease. Like most clinically relevant microorganisms, Yersinia spp. are currently subjected to intense multi-omics investigations whose numbers have increased extensively in recent years, generating massive amounts of data useful for diagnostic and therapeutic developments. The lack of a simple and centralized way to exploit these data led us to design Yersiniomics, a web-based platform allowing straightforward analysis of Yersinia omics data. Yersiniomics contains a curated multi-omics database at its core, gathering 200 genomic, 317 transcriptomic, and 62 proteomic data sets for Yersinia species. It integrates genomic, transcriptomic, and proteomic browsers, a genome viewer, and a heatmap viewer to navigate within genomes and experimental conditions. For streamlined access to structural and functional properties, it directly links each gene to GenBank, the Kyoto Encyclopedia of Genes and Genomes (KEGG), UniProt, InterPro, IntAct, and the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) and each experiment to Gene Expression Omnibus (GEO), the European Nucleotide Archive (ENA), or the Proteomics Identifications Database (PRIDE). Yersiniomics provides a powerful tool for microbiologists to assist with investigations ranging from specific gene studies to systems biology studies. IMPORTANCE The expanding genus Yersinia is composed of multiple nonpathogenic species and a few pathogenic species, including the deadly etiologic agent of plague, Yersinia pestis. In 2 decades, the number of genomic, transcriptomic, and proteomic studies on Yersinia grew massively, delivering a wealth of data. We developed Yersiniomics, an interactive web-based platform, to centralize and analyze omics data sets on Yersinia species. The platform allows user-friendly navigation between genomic data, expression data, and experimental conditions. Yersiniomics will be a valuable tool to microbiologists.

<i>Yersinia pestis</i> ssp. <i>pestis</i> Spatial MLVA25 Genotypic Structure in the Transboundary Saylyugem Natural Plague Focus

Advanced molecular-genetic methods for the diagnosis and typing of Yersinia pestis ssp. pestis in the field and clinical material are used for epidemiological surveillance of plague in the Saylyugem natural focus. The aim of the work was to study the spatial genotypic structure of Y. pestis ssp. pestis in the transboundary Saylyugem natural plague focus using MLVA25 typing. Materials and methods. The MLVA25 typing of 160 strains of Y. pestis ssp. Pestis isolated in the Saylyugem natural plague focus in 2012–2021 was carried out. Phylogenetic tree construction was performed with the help of UPGMA and MST methods. Results and discussion. The Y. pestis ssp. pestis strains isolated from the Saylyugem natural plague focus were differentiated into 15 MLVA types by the 25 VNTR loci cluster analysis. The studied strains form a homogeneous complex of MLVA25 types without marked geographical distribution across seven spatial groups. The analysis of the frequency of occurrence of the tandem repeats number for three variable loci of Y. pestis ssp. pestis strains shows the significant differences between the samples from the Mongolian and Russian parts of the Saylyugem natural plague focus. The most pronounced differences in spatial genotypic structure are traced through the yp4280ms62 locus.

Historical and Modern Classifications of the Plague Agent

The review presents the data on domestic and foreign phenotypic classifications of Yersinia pestis strains developed in the XX century; genetic classifications of the XXI century; as well as on the genealogy of ancient strains of the plague microbe, reconstructed using paleogenomic technologies. Since the discovery of the plague agent in 1894, many classifications were created that corresponded to the level of development of microbiology at that time. The intraspecific classification schemes of the XX century were based on three principles: phenotypic differences between strains, features of the species composition of carriers, and geographical affiliation. With the development of molecular microbiology early on in the XXI century, a genetic nomenclature of the branches of the pathogen evolution was developed and a number of classifications based on the analysis of the population structure of Y. pestis were created. Through the prism of the genetic diversity of Y. pestis strains from natural plague foci in Russia, near and far abroad countries, an improved classification with a division into seven subspecies has been developed: pestis, tibetica, caucasica, qinghaica, angolica, central asiatica, ulegeica, which allocates the subspecies according to the phylogenetic principle and epidemic significance. With the advancements in paleomicrobiology, prehistoric lineages of evolution have been included in the genealogy of Y. pestis, which expand the data on the intraspecific diversity of the plague microbe.

Plagued by a cryptic clock: insight and issues from the global phylogeny of Yersinia pestis

Plague has an enigmatic history as a zoonotic pathogen. This infectious disease will unexpectedly appear in human populations and disappear just as suddenly. As a result, a long-standing line of inquiry has been to estimate when and where plague appeared in the past. However, there have been significant disparities between phylogenetic studies of the causative bacterium, Yersinia pestis, regarding the timing and geographic origins of its reemergence. Here, we curate and contextualize an updated phylogeny of Y. pestis using 601 genome sequences sampled globally. Through a detailed Bayesian evaluation of temporal signal in subsets of these data we demonstrate that a Y. pestis-wide molecular clock is unstable. To resolve this, we developed a new approach in which each Y. pestis population was assessed independently, enabling us to recover substantial temporal signal in five populations, including the ancient pandemic lineages which we now estimate may have emerged decades, or even centuries, before a pandemic was historically documented from European sources. Despite this methodological advancement, we only obtain robust divergence dates from populations sampled over a period of at least 90 years, indicating that genetic evidence alone is insufficient for accurately reconstructing the timing and spread of short-term plague epidemics.

Yersinia pestis and Plague: some knowns and unknowns

Since its first identification in 1894 during the third pandemic in Hong Kong, there has been significant progress of understanding the lifestyle of Yersinia pestis, the pathogen that is responsible for plague. Although we now have some understanding of the pathogen's physiology, genetics, genomics, evolution, gene regulation, pathogenesis and immunity, there are many unknown aspects of the pathogen and its disease development. Here, we focus on some of the knowns and unknowns relating to Y. pestis and plague. We notably focus on some key Y. pestis physiological and virulence traits that are important for its mammal-flea-mammal life cycle but also its emergence from the enteropathogen Yersinia pseudotuberculosis. Some aspects of the genetic diversity of Y. pestis, the distribution and ecology of plague as well as the medical countermeasures to protect our population are also provided. Lastly, we present some biosafety and biosecurity information related to Y. pestis and plague.

Retrospective analysis of dissemination of the 2.MED1 phylogenetic branch of Yersinia pestis in the Caucasus

The 2.MED1 phylogenetic branch of Yersinia pestis of the medieval biovar became widespread in the Caspian Sea region, the Caucasus, and the Northern Aral Sea region in the 20th century, causing outbreaks and epizootics of plague there. Some of the formed natural foci of 2.MED1 still show epizootic activity and retain their epidemic potential. In this work, we carried out a phylogenetic analysis of 46 Y. pestis strains of the medieval biovar isolated in the Caucasus, the Caspian Sea, and the Northern Aral Sea regions during epidemic outbreaks and epizootics from 1922-2014. The obtained phylogenetic data, together with epidemiological and epizootological data accumulated over a period of about a hundred years, indicate the presence of two waves of penetration of the 2.MED1 branch into the Caucasus. The first occurred, apparently, in the first half of the 20th century as a result of the penetration of 2.MED1 from the foci of the Northern and North-Western Caspian Sea. The second wave was caused by the spread of 2.MED1 from the Northern Aral to the foci of the North-Western, Northern and Eastern Caspian Sea regions at the beginning of the second half of the 20th century, followed by introduction into the Pre-Caucasus and Transcaucasia. The rapid spread of 2.MED1 could be associated with the transfer of the pathogen by land and sea transport in the process of economic activity of the population.

Two Novel Yersinia pestis Bacteriophages with a Broad Host Range: Potential as Biocontrol Agents in Plague Natural Foci

Bacteriophages (phages) have been successfully used as disinfectors to kill bacteria in food and the environment and have been used medically for curing human diseases. The objective of this research was to elucidate the morphological and genomic characteristics of two novel Yersinia pestis phages, vB_YpeM_ MHS112 (MHS112) and vB_YpeM_GMS130 (GMS130), belonging to the genus Gaprivervirus, subfamily Tevenvirinae, family Myoviridae. Genome sequencing showed that the sizes of MHS112 and GMS130 were 170507 and 168552 bp, respectively. A total of 303 and 292 open reading frames with 2 tRNA and 3 tRNA were predicted in MHS112 and GMS130, respectively. The phylogenetic relationships were analysed among the two novel Y. pestis phages, phages in the genus Gaprivervirus, and several T4-like phages infecting the Yersinia genus. The bacteriophage MHS112 and GMS130 exhibited a wider lytic host spectrum and exhibited comparative temperature and pH stability. Such features signify that these phages do not need to rely on Y. pestis as their host bacteria in the ecological environment, while they could be based on more massive Enterobacteriales species to propagate and form ecological barriers against Y. pestis pathogens colonised in plague foci. Such characteristics indicated that the two phages have potential as biocontrol agents for eliminating the endemics of animal plague in natural plague foci.

Ecological scenario of the plague microbe <i>Yersinia pestis</i> speciation underlying adequate molecular evolutionary model

It is known that the psychrophilic pseudotuberculosis microbe serotype 1 (Y. pseudotuberculosis 0:1b) causing Far East scarlet-like fever (FESLF) an intestinal infection found in a wide range of invertebrates and vertebrates inhabiting cold regions in the Northern and Central Asia as well as Far East is direct ancestor of the plague causative agent Yersinia pestis. However, the mechanism of Y. pestis speciation remains poorly elucidated. Numerous Y. pestis phylogenies created by using molecular genetic (MG) technologies are largely contradictory, being not in line with reliable data obtained by natural science approaches (e.g., ecology, epizootology, biogeography, and paleontology), which disagree with current evolutionary doctrine (synthetic theory of evolution). The MG approach provides no definitive answer to the questions of where, when, how, and under what circumstances the species Y. pestis arose. One of the reasons for such situation might be due to inadequacy of using the molecular evolutionary model for Y. pestis phylogenetics. Knowledge of the life cycles for the ancestral pseudotuberculosis and derivative plague microbes as well as related unique environmental features allows to create a reliable ecological model for the plague microbe evolution to be further used for assessing patterns of molecular variability and building proper molecular model that might be accepted for MG-reconstruction of plague microbe history. According to the ecological model, the species Y. pestis was formed in a tritopic manner (almost) simultaneously from FESLF clones (populations) in the three geographical populations of the Mongolian marmot-tarbagan (Marmota sibirica) and the flea Oropsylla silantiewi parasitizing on it. The inducer of speciation was coupled to the last maximum (Sartan) cooling in Central Asia occurred 2215 thousand years ago. Soil cooling and deep freezing resulted in altered behavior of the marmot flea larvae with emergence of facultative hematophagy, which, in turn, led to a unique traumatic (compared to routine alimentary) infection route of sleeping marmots with FESLF and, as a result, a unique way of Y. pestis speciation. The molecular model should predict a Y. pestis peripatric tritope speciation, existing numerous parallelisms in intraspecific variability associated with tritope speciation, and the quantum principle of speciation in the highly variable heterothermic (heteroimmune) stressful marmot-flea (Marmota sibirica Oropsylla silantiewi) host-vector environment involving stress-induced mutagenesis. Such molecular model of evolution may be useful for improving molecular methodology of phylogenetic constructions for a wide range of parasitic microorganisms.

Polygenic plague resistance in the great gerbil uncovered by population sequencing

Pathogens can elicit high selective pressure on hosts, potentially altering genetic diversity over short evolutionary timescales. Intraspecific variation in immune response is observable as variable survivability from specific infections. The great gerbil (Rhombomys opimus) is a rodent plague host with a heterogenic but highly resistant phenotype. Here, we investigate the genomic basis for plague-resistant phenotypes by exposing wild-caught great gerbils to plague (Yersinia pestis). Whole genome sequencing of 10 survivors and 10 moribund individuals revealed a subset of genomic regions showing elevated differentiation. Gene ontology analysis of candidate genes in these regions demonstrated enrichment of genes directly involved in immune functions, cellular metabolism and the regulation of apoptosis as well as pathways involved in transcription, translation, and gene regulation. Transcriptomic analysis revealed that the early activated great gerbil immune response to plague consisted of classical components of the innate immune system. Our approach combining challenge experiments with transcriptomics and population level sequencing, provides new insight into the genetic background of plague-resistance and confirms its complex nature, most likely involving multiple genes and pathways of both the immune system and regulation of basic cellular functions.

<i>Yersinia pestis</i> Strains of the 1.ORI Line as Etiological Agent of the Plague Pandemic III

Yersinia pestis strains of the 1.ORI lineage originate from China as a result of evolution of the 1.ANT phylogenetic branch. Strains of the biovar orientalis are divided into three major lines of evolution: 1.ORI1, 1.ORI2, 1.ORI3. Lines 1.ORI1 and 1.ORI2 originated in China and then spread across the east and west coasts of India, respectively. Strains of the biovar orientalis have widely spread throughout the world, mainly as a result of introduction by sea. This way, the 1.ORI1 line was imported onto the territory of North America. 1.ORI2 line has spread to Southeast Asia, Africa, Europe, and South America. In addition, the strains of the biovar orientalis were brought to the territory of Australia, however, the formation of natural foci did not occur. The spread of strains to new territories during the third plague pandemic, as a rule, took place with the participation of one strain, which caused epizootics among synanthropic rodents. After that, outbreaks were recorded among the population of port cities, followed by drifting into the countryside and the formation of natural foci under suitable natural conditions. In the absence of such, the plague pathogen was eliminated from natural biotopes, and the formation of a natural focus did not occur. In recent decades, most cases of human plague in the world have been caused by strains of the biovar orientalis (1.ORI). However, the emergence and spread of the evolutionary line “1” is insufficiently studied. Currently, there is a lack of both historical data and strains that are ancestors of modern strains in many countries to clarify the details of the irradiation of strains of the biovar orientalis. As a result, the concepts of dissemination of many evolution branches of the strains, biovar orientalis are in the form of hypotheses to date. In this work, the collection and analysis of literature data on the history and epidemiology of plague over the third pandemic, a search for a connection between epidemic manifestations and the appurtenance of the strains that caused them to certain phylogenetic lineages was carried out.

EnteroBase: hierarchical clustering of 100 000s of bacterial genomes into species/subspecies and populations

The definition of bacterial species is traditionally a taxonomic issue while bacterial populations are identified by population genetics. These assignments are species specific, and depend on the practitioner. Legacy multilocus sequence typing is commonly used to identify sequence types (STs) and clusters (ST Complexes). However, these approaches are not adequate for the millions of genomic sequences from bacterial pathogens that have been generated since 2012. EnteroBase (http://enterobase.warwick.ac.uk) automatically clusters core genome MLST allelic profiles into hierarchical clusters (HierCC) after assembling annotated draft genomes from short-read sequences. HierCC clusters span core sequence diversity from the species level down to individual transmission chains. Here we evaluate HierCC's ability to correctly assign 100 000s of genomes to the species/subspecies and population levels for Salmonella, Escherichia, Clostridoides, Yersinia, Vibrio and Streptococcus. HierCC assignments were more consistent with maximum-likelihood super-trees of core SNPs or presence/absence of accessory genes than classical taxonomic assignments or 95% ANI. However, neither HierCC nor ANI were uniformly consistent with classical taxonomy of Streptococcus. HierCC was also consistent with legacy eBGs/ST Complexes in Salmonella or Escherichia and with O serogroups in Salmonella. Thus, EnteroBase HierCC supports the automated identification of and assignment to species/subspecies and populations for multiple genera. This article is part of a discussion meeting issue 'Genomic population structures of microbial pathogens'.

Subversion of GBP-mediated host defense by E3 ligases acquired during Yersinia pestis evolution

Plague has caused three worldwide pandemics in history, including the Black Death in medieval ages. Yersinia pestis, the etiological agent of plague, has evolved a powerful arsenal to disrupt host immune defenses during evolution from enteropathogenic Y. pseudotuberculosis. Here, we find that two functionally redundant E3 ligase of Y. pestis, YspE1 and YspE2, can be delivered via type III secretion injectisome into host cytosol where they ubiquitinate multiple guanylate-binding proteins (GBPs) for proteasomal degradation. However, Y. pseudotuberculosis has no such capability due to lacking functional YspE1/2 homologs. YspE1/2-mediated GBP degradations significantly promote the survival of Y. pestis in macrophages and strongly inhibit inflammasome activation. By contrast, Gbp^{chr3−/−, chr5−/−} macrophages exhibit much lowered inflammasome activation independent of YspE1/2, accompanied with an enhanced replication of Y. pestis. Accordingly, Gbp^{chr3−/−, chr5−/−} mice are more susceptible to Y. pestis. We demonstrate that Y. pestis utilizes E3 ligases to subvert GBP-mediated host defense, which appears to be newly acquired by Y. pestis during evolution.

Guanylate-binding proteins (GBPs) recognize pathogen containing vacuoles, leading to lysis of this intracellular niche and induction of inflammasomes. Here, Cao et al. show that Y. pestis, the causative agent of plague, secret two functionally redundant E3 ligase, YspE1 and YspE2, into the host’s cytosol to ubiquitinate multiple GBPs for proteasomal degradation to subvert host immune defense. This capability appears to be newly acquired by Y. pestis during evolution, since its closely related progenitor Y. pseudotuberculosis is unable to do so.

Ancient Yersinia pestis and Salmonella enterica genomes from Bronze Age Crete

During the late 3^rd millennium BCE, the Eastern Mediterranean and Near East witnessed societal changes in many regions, which are usually explained with a combination of social and climatic factors.^1-4 However, recent archaeogenetic research forces us to rethink models regarding the role of infectious diseases in past societal trajectories.⁵ The plague bacterium Yersinia pestis, which was involved in some of the most destructive historical pandemics,^5-8 circulated across Eurasia at least from the onset of the 3^rd millennium BCE,^9-13 but the challenging preservation of ancient DNA in warmer climates has restricted the identification of Y.pestis from this period to temperate climatic regions. As such, evidence from culturally prominent regions such as the Eastern Mediterranean is currently lacking. Here, we present genetic evidence for the presence of Y. pestis and Salmonella enterica, the causative agent of typhoid/enteric fever, from this period of transformation in Crete, detected at the cave site Hagios Charalambos. We reconstructed one Y. pestis genome that forms part of a now-extinct lineage of Y. pestis strains from the Late Neolithic and Bronze Age that were likely not yet adapted for transmission via fleas. Furthermore, we reconstructed two ancient S. enterica genomes from the Para C lineage, which cluster with contemporary strains that were likely not yet fully host adapted to humans. The occurrence of these two virulent pathogens at the end of the Early Minoan period in Crete emphasizes the necessity to re-introduce infectious diseases as an additional factor possibly contributing to the transformation of early complex societies in the Aegean and beyond.

Identifying treponemal disease in early East Asia

Objectives

Historic records suggest that a virulent form of treponematosis, sexually transmitted syphilis was introduced to Asia from Europe by the da Gama crew, who landed in India in 1498. Our objective is to assess the gross pathology of human skeletal remains from the Tang dynasty of China to test the presence of treponemal infection in East Asia before 1498. We interpret this paleopathological evidence in the context of site ecology and sociocultural changes during the Tang dynasty.

Materials and methods

We examined the gross pathology of 1598 human skeletons from Xingfulindai (AD 618 to AD 1279) archeological site located on the Central Plain of China. Using the modified diagnostic criteria defined by Hackett's classical work, we classify the pathology as consistent, strongly suggestive, or pathognomonic for treponemal infection.

Results

Twelve adult individuals from Xingfulindai had bone lesions suggestive of systemic pathology. Two of these individuals displayed a combination of lesion patterns pathognomonic of treponemal disease and one had lesions consistent with treponematosis. The radiocarbon dates for the bone samples from these skeletons place them before AD 1200.

Conclusions

The location of Xingfulindai in a continental climatic zone is not typical for yaws and bejel ecology, because these strains occur in the tropics, or in hot, dry environments, respectively. The urban setting, where there is documented evidence for increased interaction between multiple ethnic groups and a developed institution of courtesans during the Tang dynasty, favors sexually transmitted syphilis as the more likely diagnosis. This study supports an earlier spread of syphilis to China than 1498.

The source of the Black Death in fourteenth-century central Eurasia

The origin of the medieval Black Death pandemic (AD 1346-1353) has been a topic of continuous investigation because of the pandemic's extensive demographic impact and long-lasting consequences1,2. Until now, the most debated archaeological evidence potentially associated with the pandemic's initiation derives from cemeteries located near Lake Issyk-Kul of modern-day Kyrgyzstan1,3-9. These sites are thought to have housed victims of a fourteenth-century epidemic as tombstone inscriptions directly dated to 1338-1339 state 'pestilence' as the cause of death for the buried individuals9. Here we report ancient DNA data from seven individuals exhumed from two of these cemeteries, Kara-Djigach and Burana. Our synthesis of archaeological, historical and ancient genomic data shows a clear involvement of the plague bacterium Yersinia pestis in this epidemic event. Two reconstructed ancient Y. pestis genomes represent a single strain and are identified as the most recent common ancestor of a major diversification commonly associated with the pandemic's emergence, here dated to the first half of the fourteenth century. Comparisons with present-day diversity from Y. pestis reservoirs in the extended Tian Shan region support a local emergence of the recovered ancient strain. Through multiple lines of evidence, our data support an early fourteenth-century source of the second plague pandemic in central Eurasia.

Pathogenomes and variations in Shiga toxin production among geographically distinct clones of Escherichia coli O113:H21

Infections with globally disseminated Shiga toxin-producing Escherichia coli (STEC) of the O113:H21 serotype can progress to severe clinical complications, such as hemolytic uremic syndrome (HUS). Two phylogeographically distinct clonal complexes have been established by multi locus sequence typing (MLST). Infections with ST-820 isolates circulating exclusively in Australia have caused severe human disease, such as HUS. Conversely, ST-223 isolates prevalent in the US and outside Australia seem to rarely cause severe human disease but are frequent contaminants. Following a genomic epidemiology approach, we wanted to gain insights into the underlying cause for this disparity. We examined the plasticity in the genome make-up and Shiga toxin production in a collection of 20 ST-820 and ST-223 strains isolated from produce, the bovine reservoir, and clinical cases. STEC are notorious for assembly into fragmented draft sequences when using short-read sequencing technologies due to the extensive and partly homologous phage complement. The application of long-read technology (LRT) sequencing yielded closed reference chromosomes and plasmids for two representative ST-820 and ST-223 strains. The established high-resolution framework, based on whole genome alignments, single nucleotide polymorphism (SNP)-typing and MLST, includes the chromosomes and plasmids of other publicly available O113:H21 sequences and allowed us to refine the phylogeographical boundaries of ST-820 and ST-223 complex isolates and to further identify a historic non-shigatoxigenic strain from Mexico as a quasi-intermediate. Plasmid comparison revealed strong correlations between the strains' featured pO113 plasmid genotypes and chromosomally inferred ST, which suggests coevolution of the chromosome and virulence plasmids. Our pathogenicity assessment revealed statistically significant differences in the Stx_2a-production capabilities of ST-820 as compared to ST-223 strains under RecA-induced Stx phage mobilization, a condition that mimics Stx-phage induction. These observations suggest that ST-820 strains may confer an increased pathogenic potential in line with the strain-associated epidemiological metadata. Still, some of the tested ST-223 cultures sourced from contaminated produce or the bovine reservoir also produced Stx at levels comparable to those of ST-820 isolates, which calls for awareness and for continued surveillance of this lineage.

Knowledge of Emerging and Reemerging Infectious Diseases in the Public of Guangzhou, Southern China

Objective

The objective of this study is to get the overall picture about the knowledge of emerging and reemerging infectious diseases in public in Guangzhou and provide a scientific basis for developing health information strategies.

Methods

We used the structured questionnaire to interview 1,000 Guangzhou residents by health enquiry hotline. Descriptive analysis was presented to evaluate the knowledge of the participants. Multiple logistic regression model was performed to determine the influence factors for knowledge of emerging and reemerging infectious diseases

Results

A total of 801 individuals completed the survey. About one-third had heard of Middle East respiratory syndrome (MERS) and Zika, whereas Ebola and plague about 50%. A total of 32.08% participants had never heard of any of the four diseases. Only 2.08% knew the sexual transmission of Zika and 90.17% had no idea about the epidemic region of plague. No more than 15% knew they should check their health status after returning from the epidemic region. Education level and income were the key factors that influenced knowledge rate.

Conclusions

The low-level knowledge called for the improvement in health information to the public, especially those with low level of education and income. Effective and precise health information was urged to carry out to improve the prevention for the emerging and reemerging infectious diseases.

Sialoglycan binding patterns of bacterial AB5 toxin B subunits correlate with host range and toxicity, indicating evolution independent of A subunits

Many pathogenic bacteria secrete AB₅ toxins that can be virulence factors. Cytotoxic A subunits are delivered to the cytosol following B subunit binding to specific host cell surface glycans. Some B subunits are not associated with A subunits, for example, YpeB of Yersinia pestis, the etiologic agent of plague. Plague cannot be eradicated because of Y. pestis' adaptability to numerous hosts. We previously showed selective binding of other B₅ pentamers to a sialoglycan microarray, with sialic acid (Sia) preferences corresponding to those prominently expressed by various hosts, for example, N-acetylneuraminic acid (Neu5Ac; prominent in humans) or N-glycolylneuraminic acid (Neu5Gc; prominent in ruminant mammals and rodents). Here, we report that A subunit phylogeny evolved independently of B subunits and suggest a future B subunit nomenclature based on bacterial species names. We also found via phylogenetic analysis of B subunits, which bind Sias, that homologous molecules show poor correlation with species phylogeny. These data indicate ongoing lateral gene transfers between species, including mixing of A and B subunits. Consistent with much broader host range of Y. pestis, we show that YpeB recognizes all mammalian Sia types, except for 4-O-acetylated ones. Notably, YpeB alone causes dose-dependent cytotoxicity, which is abolished by a mutation (Y77F) eliminating Sia recognition, suggesting that cell proliferation and death are promoted via lectin-like crosslinking of cell surface sialoglycoconjugates. These findings help explain the host range of Y. pestis and could be important for pathogenesis. Overall, our data indicate ongoing rapid evolution of both host Sias and pathogen toxin-binding properties.

PgtE Enzyme of Salmonella enterica Shares the Similar Biological Roles to Plasminogen Activator (Pla) in Interacting With DEC-205 (CD205), and Enhancing Host Dissemination and Infectivity by Yersinia pestis

Yersinia pestis, the cause of plague, is a newly evolved Gram-negative bacterium. Through the acquisition of the plasminogen activator (Pla), Y. pestis gained the means to rapidly disseminate throughout its mammalian hosts. It was suggested that Y. pestis utilizes Pla to interact with the DEC-205 (CD205) receptor on antigen-presenting cells (APCs) to initiate host dissemination and infection. However, the evolutionary origin of Pla has not been fully elucidated. The PgtE enzyme of Salmonella enterica, involved in host dissemination, shows sequence similarity with the Y. pestis Pla. In this study, we demonstrated that both Escherichia coli K-12 and Y. pestis bacteria expressing the PgtE-protein were able to interact with primary alveolar macrophages and DEC-205-transfected CHO cells. The interaction between PgtE-expressing bacteria and DEC-205-expressing transfectants could be inhibited by the application of an anti-DEC-205 antibody. Moreover, PgtE-expressing Y. pestis partially re-gained the ability to promote host dissemination and infection. In conclusion, the DEC-205-PgtE interaction plays a role in promoting the dissemination and infection of Y. pestis, suggesting that Pla and the PgtE of S. enterica might share a common evolutionary origin.

Plague risk in the western United States over seven decades of environmental change

After several pandemics over the last two millennia, the wildlife reservoirs of plague (Yersinia pestis) now persist around the world, including in the western United States. Routine surveillance in this region has generated comprehensive records of human cases and animal seroprevalence, creating a unique opportunity to test how plague reservoirs are responding to environmental change. Here, we test whether animal and human data suggest that plague reservoirs and spillover risk have shifted since 1950. To do so, we develop a new method for detecting the impact of climate change on infectious disease distributions, capable of disentangling long-term trends (signal) and interannual variation in both weather and sampling (noise). We find that plague foci are associated with high-elevation rodent communities, and soil biochemistry may play a key role in the geography of long-term persistence. In addition, we find that human cases are concentrated only in a small subset of endemic areas, and that spillover events are driven by higher rodent species richness (the amplification hypothesis) and climatic anomalies (the trophic cascade hypothesis). Using our detection model, we find that due to the changing climate, rodent communities at high elevations have become more conducive to the establishment of plague reservoirs-with suitability increasing up to 40% in some places-and that spillover risk to humans at mid-elevations has increased as well, although more gradually. These results highlight opportunities for deeper investigation of plague ecology, the value of integrative surveillance for infectious disease geography, and the need for further research into ongoing climate change impacts.

Local-scale diversity of Yersinia pestis: A case study from Ambohitromby, Ankazobe District, Madagascar

Plague is a re-emerging zoonotic disease and a major public health concern in several portions of the world, especially in Madagascar. We report on the presence of different subtypes of Yersinia pestis co-occurring in the same locality. After confirmation of a human plague case in Ambohitromby Commune (Ankazobe District) via isolation of Y. pestis, we undertook small mammal trapping to identify the circulation of Y. pestis amongst rodents in this locality; blood samples were collected from rodents for seroprevalence analysis. Of the 60 individuals of Rattus rattus captured, one yielded an isolate of Y. pestis, 13 others were positive for F1 antigen of Y. pestis using a rapid diagnostic test, and 4 were PCR positive targeting the caf1 and pla genes; 28/60 (46.7%) of the captured R. rattus were seropositive for Y. pestis. Whole-genome SNP analyses revealed that the two isolates obtained from the human case, and the R. rattus belonged to two different subtypes of Y. pestis (s05 and s13, respectively) that were circulating concurrently in Ambohitromby in 2016. Three Y. pestis subtypes (s03, s05 and s13) have now been isolated from Ambohitromby. Subtype s05 had been persisting there for >10 years but one or both of the other subtypes may have been introduced from the Central Highlands region as they were not observed in previous years (s13) or only observed once previously (s03). High seroprevalence against Y. pestis in R. rattus suggests that a portion of the local murine population may have acquired resistance to Y. pestis. Future research should focus on genomically characterizing Y. pestis strains circulating in Ankazobe District and other plague-endemic regions of Madagascar to better understand the overall phylogeography of Y. pestis.

Analysis of Genomic DNA from Medieval Plague Victims Suggests Long-Term Effect of Yersinia pestis on Human Immunity Genes

Pathogens and associated outbreaks of infectious disease exert selective pressure on human populations, and any changes in allele frequencies that result may be especially evident for genes involved in immunity. In this regard, the 1346-1353 Yersinia pestis-caused Black Death pandemic, with continued plague outbreaks spanning several hundred years, is one of the most devastating recorded in human history. To investigate the potential impact of Y. pestis on human immunity genes, we extracted DNA from 36 plague victims buried in a mass grave in Ellwangen, Germany in the 16th century. We targeted 488 immune-related genes, including HLA, using a novel in-solution hybridization capture approach. In comparison with 50 modern native inhabitants of Ellwangen, we find differences in allele frequencies for variants of the innate immunity proteins Ficolin-2 and NLRP14 at sites involved in determining specificity. We also observed that HLA-DRB1*13 is more than twice as frequent in the modern population, whereas HLA-B alleles encoding an isoleucine at position 80 (I-80+), HLA C*06:02 and HLA-DPB1 alleles encoding histidine at position 9 are half as frequent in the modern population. Simulations show that natural selection has likely driven these allele frequency changes. Thus, our data suggest that allele frequencies of HLA genes involved in innate and adaptive immunity responsible for extracellular and intracellular responses to pathogenic bacteria, such as Y. pestis, could have been affected by the historical epidemics that occurred in Europe.

The emergence and transmission of COVID-19 in European countries, 2019-2020: a comprehensive review of timelines, cases and containment

When it emerged in late 2019, COVID-19 was carried via travelers to Germany, France and Italy, where freedom of movement accelerated its transmission throughout Europe. However, effective non-pharmaceutical interventions introduced by European governments led to containment of the rapid increase in cases within European nations. Electronic searches were performed to obtain the number of confirmed cases, incident rates and non-pharmaceutical government measures for each European country. The spread and impact of non-pharmaceutical interventions throughout Europe were assessed and visualized. Specifically, heatmaps were used to represent the number of confirmed cases and incident rates for each of the countries over time. In addition, maps were created showing the number of confirmed cases and incident rates in Europe on three different dates (15 March, 15 April and 15 May 2020), which allowed us to assess the geographic and temporal patterns of the disease.

Genetic source tracking of human plague cases in Inner Mongolia-Beijing, 2019

On 12 November 2019, one couple from the Sonid Left Qi (County) in the Inner Mongolia Autonomous Region was diagnosed with pneumonic plague in Beijing. The wife acquired the infection from her husband. Thereafter, two bubonic plague cases were identified in Inner Mongolia on November 16th and 24th. In this study, genome-wide single nucleotide polymorphism (SNP) analysis was used to identify the phylogenetic relationship of Yersinia pestis strains isolated in Inner Mongolia. Strains isolated from reservoirs in 2018 and 2019 in Inner Mongolia, together with the strain isolated from Patient C, were further clustered into 2.MED3m, and two novel lineages (2.MED3q, 2.MED3r) in the 2.MED3 population. According to the analysis of PCR-based molecular subtyping methods, such as the MLVA 14 scheme and seven SNP allele sequencing, Patients A/B and D were classified as 2.MED3m. In addition, strains from rodents living near the patients' residences were clustered into the same lineage as patients. Such observations indicated that human plague cases originated from local reservoirs. Corresponding phylogenetic analysis also indicated that rodent plague strains in different areas in Inner Mongolia belong to different epizootics rather than being caused by spreading from the same epizootic in Meriones unguiculatus in 2019.

A novel mechanism of streptomycin resistance in Yersinia pestis: Mutation in the rpsL gene

Streptomycin is considered to be one of the effective antibiotics for the treatment of plague. In order to investigate the streptomycin resistance of Y. pestis in China, we evaluated streptomycin susceptibility of 536 Y. pestis strains in China in vitro using the minimal inhibitory concentration (MIC) and screened streptomycin resistance-associated genes (strA and strB) by PCR method. A clinical Y. pestis isolate (S19960127) exhibited high-level resistance to streptomycin (the MIC was 4,096 mg/L). The strain (biovar antiqua) was isolated from a pneumonic plague outbreak in 1996 in Tibet Autonomous Region, China, belonging to the Marmota himalayana Qinghai–Tibet Plateau plague focus. In contrast to previously reported streptomycin resistance mediated by conjugative plasmids, the genome sequencing and allelic replacement experiments demonstrated that an rpsL gene (ribosomal protein S12) mutation with substitution of amino-acid 43 (K43R) was responsible for the high-level resistance to streptomycin in strain S19960127, which is consistent with the mutation reported in some streptomycin-resistant Mycobacterium tuberculosis strains. Streptomycin is used as the first-line treatment against plague in many countries. The emergence of streptomycin resistance in Y. pestis represents a critical public health problem. So streptomycin susceptibility monitoring of Y. pestis isolates should not only include plasmid-mediated resistance but also include the ribosomal protein S12 gene (rpsL) mutation, especially when treatment failure is suspected due to antibiotic resistance.

The plague natural foci are widely distributed in the world, and correspondingly, the plague still poses a significant threat to human health in some countries with endemic plague foci. Streptomycin is used as the first-line treatment against plague in many countries for the antibiotic is considered to be one of the effective antibiotics, particularly for the treatment of pneumonic plague. The resistance to streptomycin had been reported in Y. pestis strains from Madagascar in previous studies. In this study, we reported the high-level resistance to streptomycin in a clinical isolate of Y. pestis from a pneumonic patient in Tibet Autonomous Region, China, and a novel mechanism of streptomycin resistance, i.e. mutation in the rpsL gene were identified. The knowledge acquired about streptomycin resistance in Y. pestis will remain of great practical value. For the emergence of resistance to streptomycin in Y. pestis would render the treatment failure, thus corresponding antibiotic monitoring should be routinely carried out in countries threatened by plague. In addition, based on our further understanding about streptomycin resistance of Y. pestis isolates, such monitoring should not only include plasmid-mediated resistance but also include the ribosomal protein S12 gene (rpsL) mutation in Y. pestis isolates.

Yersinia pestis strains isolated in natural plague foci of Caucasus and Transcaucasia in the context of the global evolution of species

BACKGROUND

Plague is a highly dangerous vector-borne infectious disease that has left a significant mark on history of humankind. There are 13 natural plague foci in the Caucasus, located on the territory of the Russian Federation, Azerbaijan, Armenia and Georgia. We performed whole-genome sequencing of Y. pestis strains, isolated in the natural foci of the Caucasus and Transcaucasia. Using the data of whole-genome SNP analysis and Bayesian phylogeny methods, we carried out an evolutionary-phylogeographic analysis of modern population of the plague pathogen in order to determine the phylogenetic relationships of Y. pestis strains from the Caucasus with the strains from other countries.

RESULTS

We used 345 Y. pestis genomes to construct a global evolutionary phylogenetic reconstruction of species based on whole-genome SNP analysis. The genomes of 16 isolates were sequenced in this study, the remaining 329 genomes were obtained from the GenBank database. Analysis of the core genome revealed 3315 SNPs that allow differentiation of strains. The evolutionary phylogeographic analysis showed that the studied Y. pestis strains belong to the genetic lineages 0.PE2, 2.MED0, and 2.MED1. It was shown that the Y. pestis strains isolated on the territory of the East Caucasian high-mountain, the Transcaucasian high-mountain and the Priaraksinsky low-mountain plague foci belong to the most ancient of all existing genetic lineages - 0.PE2.

CONCLUSIONS

On the basis of the whole-genome SNP analysis of 345 Y. pestis strains, we describe the modern population structure of the plague pathogen and specify the place of the strains isolated in the natural foci of the Caucasus and Transcaucasia in the structure of the global population of Y. pestis. As a result of the retrospective evolutionary-phylogeographic analysis of the current population of the pathogen, we determined the probable time frame of the divergence of the genetic lineages of Y. pestis, as well as suggested the possible paths of the historical spread of the plague pathogen.

New Genotype of Yersinia pestis Found in Live Rodents in Yunnan Province, China

Yunnan Province, China is thought to be the original source of biovar Orientalis of Yersinia pestis, the causative agent of the third plague pandemic that has spread globally since the end of the 19th century. Although encompassing a large area of natural plague foci, Y. pestis strains have rarely been found in live rodents during surveillance in Yunnan, and most isolates are from rodent corpses and their fleas. In 2017, 10 Y. pestis strains were isolated from seven live rodents and three fleas in Heqing County of Yunnan. These strains were supposed to have low virulence to local rodents Eothenomys miletus and Apodemus chevrieri because the rodents were healthy and no dead animals were found in surrounding areas, as had occurred in previous epizootic disease. We performed microscopic and biochemical examinations of the isolates, and compared their whole-genome sequences and transcriptome with those of 10 high virulence Y. pestis strains that were isolated from nine rodents and one parasitic flea in adjacent city (Lijiang). We analyzed the phenotypic, genomic, and transcriptomic characteristics of live rodent isolates. The isolates formed a previously undefined monophyletic branch of Y. pestis that was named 1.IN5. Six SNPs, two indels, and one copy number variation were detected between live rodent isolates and the high virulence neighbors. No obvious functional consequence of these variations was found according to the known annotation information. Among genes which expression differential in the live rodent isolates compared to their high virulent neighbors, we found five iron transfer related ones that were significant up-regulated (| log₂ (FC) | > 1, p.adjust < 0.05), indicating these genes may be related to the low-virulence phenotype. The novel genotype of Y. pestis reported here provides further insights into the evolution and spread of plague as well as clues that may help to decipher the virulence mechanism of this notorious pathogen.

HierCC: A multi-level clustering scheme for population assignments based on core genome MLST

Motivation

Routine infectious disease surveillance is increasingly based on large-scale whole-genome sequencing databases. Real-time surveillance would benefit from immediate assignments of each genome assembly to hierarchical population structures. Here we present pHierCC, a pipeline that defines a scalable clustering scheme, HierCC, based on core genome multi-locus typing that allows incremental, static, multi-level cluster assignments of genomes. We also present HCCeval, which identifies optimal thresholds for assigning genomes to cohesive HierCC clusters. HierCC was implemented in EnteroBase in 2018 and has since genotyped >530 000 genomes from Salmonella, Escherichia/Shigella, Streptococcus, Clostridioides, Vibrio and Yersinia.

Availability and implementation

https://enterobase.warwick.ac.uk/ and Source code and instructions: https://github.com/zheminzhou/pHierCC

Supplementary information

Supplementary data are available at Bioinformatics online.

Molecular characterisation of Mycobacterium avium subsp. paratuberculosis in Australia

BACKGROUND

Mycobacterium avium subsp. paratuberculosis (Map) causes Johne's disease (JD), a chronic enteritis widespread in ruminants, resulting in substantial economic losses, especially to the dairy industry. Understanding the genetic diversity of Map in Australia will assist epidemiological studies for tracking disease transmission and identify subtype characteristics for use in development of improved diagnostic typing methods. Here we investigated the phylogenetic relationships of 351 Map isolates and compared different subtyping methods to assess their suitability for use in diagnostics and accuracy.

RESULTS

SNP-based phylogenetic analysis of 228 Australian isolates and 123 publicly available international isolates grouped Type S and Type C strains into two distinct lineages. Type C strains were highly monomorphic with only 20 SNP differences separating them. Type S strains, when aligned separately to the Telford strain, fell into two distinct clades: The first clade contained seven international isolates while the second clade contained one international isolate from Scotland and all 59 Australian isolates. The Australian Type B strain clustered with US bison strains. IS1311 PCR and Restriction Enzyme Analysis (REA) intermittently generated incorrect results when compared to Long Sequence Polymorphism (LSP) analysis, whole genome SNP-based phylogenetic analysis, IS1311 sequence alignment and average nucleotide identity (ANI). These alternative methods generated consistent Map typing results. A published SNP based assay for genotyping Map was found to be unsuitable for differentiating between Australian and international strain types of Map.

CONCLUSION

This is the first phylogenetic analysis of Australian Map isolates. The Type C lineage was highly monomorphic, and the Type S lineage clustered all Australian isolates into one clade with a single Scottish sheep strain. The Australian isolate classified as Type B by IS1311 PCR and REA is likely to be descended from bison and most closely related to US bison strains. Limitations of the current typing methods were identified in this study.

Evolution and circulation of Yersinia pestis in the Northern Caspian and Northern Aral Sea regions in the 20th-21st centuries

According to the whole genome SNP analysis of 38 Yersinia pestis strains isolated in the foci of the Northern Caspian and Northern Aral Sea regions in the 20th–early 21st centuries, between 1912 and 2015, the spatial and temporal structure of the 2.MED population of a medieval biovar in this region was determined. A phylogenetic branch 2.MED4 was identified which preceded the 2.MED1 branch that diverged later. 2.MED1 strains became the etiological agent of high-mortality plague outbreaks that occurred in the Northern Caspian region at the beginning of the 20th century. Later in the 20th century, the 2.MED1 branch became widespread in the Caspian Sea region, Caucasus, and vast areas of Central Asia. Based on the data of phylogenetic analysis, as well as epidemiological and epizootiological data, we reconstructed the paths of spread of the 2.MED1 branch in the Northern Caspian Sea region and in the Northern subzone of the Central Asian deserts. It is shown, that the reason for the activation of plague foci in the Northern Caspian region in the second half of the 20th century after a long inter-epizootic period caused by cyclical climate warming was the return of 2.MED1 from the foci of the Northern Aral Sea region. This led to the formation of stable plague foci in the Northern Caspian Sea region and Pre-Caucasus, which manifested epizootic activity in the second half of the 20th and early 21st centuries.