Yersinia pestis lineages in Mongolia

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.

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.

Distribution and Characteristics of Human Plague Cases and Yersinia pestis Isolates from 4 Marmota Plague Foci, China, 1950-2019

We analyzed epidemiologic characteristics and distribution of 1,067 human plague cases and 5,958 Yersinia pestis isolates collected from humans, host animals, and insect vectors during 1950–2019 in 4 Marmota plague foci in China. The case-fatality rate for plague in humans was 68.88%; the overall trend slowly decreased over time but fluctuated greatly. Most human cases (98.31%) and isolates (82.06%) identified from any source were from the Marmota himalayana plague focus. The tendency among human cases could be divided into 3 stages: 1950–1969, 1970–2003, and 2004–2019. The Marmota sibirica plague focus has not had identified human cases nor isolates since 1926. However, in the other 3 foci, Y. pestis continues to circulate among animal hosts; ecologic factors might affect local Y. pestis activity. Marmota plague foci are active in China, and the epidemic boundary is constantly expanding, posing a potential threat to domestic and global public health.

Biological activity evaluations of chemical constituents derived from Mongolian medicinal forage plants and their applications in combating infectious diseases and addressing health problems in humans and livestock

Mongolian nomadic people possess traditional knowledge of wild plants that grow in their areas of habitation. Many of these are forage plants in nature and are consumed by livestock. However, these plants are known to have medicinal and/or toxic properties. To establish a scientific understanding of the plants, and in turn, offer sound knowledge on their applications and effective use, it is essential to collect data pertaining to the chemical constituents of each plant. Therefore, the first objective of this study was to identify and determine the structural constituents of the forage plants that were available to our research group. Furthermore, in an attempt to demonstrate the biological activities of the isolated chemical compounds, we focused on solving some of the social issues affecting Mongolian communities, including protozoan diseases affecting livestock, vectors of infectious diseases, and the general health of humans and their livestock. The results of the chemical constituents derived from Mongolian medicinal plants and their biological activities that were studied in the recent decade are also described herein.

Yersinia pestis: the Natural History of Plague

The Gram-negative bacterium Yersinia pestis is responsible for deadly plague, a zoonotic disease established in stable foci in the Americas, Africa, and Eurasia. Its persistence in the environment relies on the subtle balance between Y. pestis-contaminated soils, burrowing and nonburrowing mammals exhibiting variable degrees of plague susceptibility, and their associated fleas. Transmission from one host to another relies mainly on infected flea bites, inducing typical painful, enlarged lymph nodes referred to as buboes, followed by septicemic dissemination of the pathogen. In contrast, droplet inhalation after close contact with infected mammals induces primary pneumonic plague. Finally, the rarely reported consumption of contaminated raw meat causes pharyngeal and gastrointestinal plague. Point-of-care diagnosis, early antibiotic treatment, and confinement measures contribute to outbreak control despite residual mortality. Mandatory primary prevention relies on the active surveillance of established plague foci and ectoparasite control. Plague is acknowledged to have infected human populations for at least 5,000 years in Eurasia. Y. pestis genomes recovered from affected archaeological sites have suggested clonal evolution from a common ancestor shared with the closely related enteric pathogen Yersinia pseudotuberculosis and have indicated that ymt gene acquisition during the Bronze Age conferred Y. pestis with ectoparasite transmissibility while maintaining its enteric transmissibility. Three historic pandemics, starting in 541 AD and continuing until today, have been described. At present, the third pandemic has become largely quiescent, with hundreds of human cases being reported mainly in a few impoverished African countries, where zoonotic plague is mostly transmitted to people by rodent-associated flea bites.

Historical and genomic data reveal the influencing factors on global transmission velocity of plague during the Third Pandemic

Quantitative knowledge about which natural and anthropogenic factors influence the global spread of plague remains sparse. We estimated the worldwide spreading velocity of plague during the Third Pandemic, using more than 200 years of extensive human plague case records and genomic data, and analyzed the association of spatiotemporal environmental factors with spreading velocity. Here, we show that two lineages, 2.MED and 1.ORI3, spread significantly faster than others, possibly reflecting differences among strains in transmission mechanisms and virulence. Plague spread fastest in regions with low population density and high proportion of pasture- or forestland, findings that should be taken into account for effective plague monitoring and control. Temperature exhibited a nonlinear, U-shaped association with spread speed, with a minimum around 20 °C, while precipitation showed a positive association. Our results suggest that global warming may accelerate plague spread in warm, tropical regions and that the projected increased precipitation in the Northern Hemisphere may increase plague spread in relevant regions.

The CRISPR-Cas system in Enterobacteriaceae

In nature, microorganisms are constantly exposed to multiple viral infections and thus have developed many strategies to survive phage attack and invasion by foreign DNA. One of such strategies is the clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated proteins (Cas) bacterial immunological system. This defense mechanism is widespread in prokaryotes including several families such as Enterobacteriaceae. Much knowledge about the CRISPR-Cas system has been generated, including its biological functions, transcriptional regulation, distribution, utility as a molecular marker and as a tool for specific genome editing. This review focuses on these aspects and describes the state of the art of the CRISPR-Cas system in the Enterobacteriaceae bacterial family.

Genotypic Expansion Within the Population Structure of Classical Brucella Species Revealed by MLVA16 Typing of 1404 Brucella Isolates From Different Animal and Geographic Origins, 1974-2006

Previous studies have shown the usefulness of MLVA16 as a rapid molecular identification and classification method for Brucella species and biovars including recently described novel Brucella species from wildlife. Most studies were conducted on a limited number of strains from limited geographic/host origins. The objective of this study was to assess genetic diversity of Brucella spp. by MLVA16 on a larger scale. Thus, 1404 animal or human isolates collected from all parts of the world over a period of 32 years (1974-2006) were investigated. Selection of the 1404 strains was done among the approximately 4000 strains collection of the BCCN (Brucella Culture Collection Nouzilly), based on classical biotyping and on the animal/human/geographic origin over the time period considered. MLVA16 was performed on extracted DNAs using high throughput capillary electrophoresis. The 16 loci were amplified in four multiplex PCR reactions. This large scale study firstly confirmed the accuracy of MLVA16 typing for Brucella species and biovar identification and its congruence with the recently described Extended Multilocus Sequence Analysis. In addition, it allowed identifying novel MLVA11 (based upon 11 slowly evolving VNTRs) genotypes representing an increase of 15% relative to the previously known Brucella MLVA11 genotypes. Cluster analysis showed that among the MLVA16 genotypes some were genetically more distant from the major classical clades. For example new major clusters of B. abortus biovar 3 isolated from cattle in Sub-Saharan Africa were identified. For other classical species and biovars this study indicated also genotypic expansion within the population structure of classical Brucella species. MLVA proves to be a powerful tool to rapidly assess genetic diversity of bacterial populations on a large scale, as here on a large collection of strains of the genomically homogeneous genus Brucella. The highly discriminatory power of MLVA appears of particular interest as a first step for selection of Brucella strains for whole-genome sequencing. The MLVA data of this study were added to the public Brucella MLVA database at http://microbesgenotyping.i2bc.paris-saclay.fr. Current version Brucella_4_3 comprises typing data from more than 5000 strains including in silico data analysis of public whole genome sequence datasets.

Investigation of Yersinia pestis and Yersinia pseudotuberculosis strains from Georgia and neighboring countries in the Caucasus by high-density SNP microarray

Yersinia pestis, the causative agent of plague, is a recently evolved clone of the enteropathogenic bacterium Yersinia pseudotuberculosis. Y. pestis has been extensively studied for decades; however, there are insufficient data about the intra-species diversity of this microorganism in certain parts of the world, including the Caucasus region. Using a high-density single-nucleotide polymorphism (SNP) microarray, we genotyped a total of 46 Y. pestis isolates from two plague foci in Georgia and neighboring Caucasus countries together with 12 Y. pseudotuberculosis isolates from Georgia. The genotyping microarray captured a total of 13,525 SNP positions across the Y. pestis and Y. pseudotuberculosis genomes and plasmids with high-throughput capability and superior reproducibility. From this analysis, we confirmed the presence of two independent and relatively distant phylogenetic groups of Y. pestis in the Caucasus region. The signature SNP patterns identified from this study will allow assay development for plague surveillance and pseudotuberculosis diagnostics.

Reemergence of human plague in Yunnan, China in 2016

The third plague pandemic originated from Yunnan Province, China in the middle of the 19th century. The last human plague epidemic in Yunnan occurred from 1986-2005. On June 6, 2016, a case of human plague was reported in the Xishuangbanna Prefecture, Yunnan. The patient suffered from primary septicemic plague after exposure to a dead house rat (Rattus flavipectus), which has been identified as the main plague reservoir in the local epizootic area. Moreover, a retrospective investigation identified another bubonic plague case in this area. Based on these data, human plague reemerged after a silent period of ten years. In this study, three molecular typing methods, including a clustered regularly interspaced short palindromic repeats (CRISPR) analysis, different region analysis (DFR), and multiple-locus variable number of tandem repeats analysis (MLVA), were used to illustrate the molecular characteristics of Yersinia pestis (Y. pestis) strains isolated in Yunnan. The DFR profiles of the strains isolated in Yunnan in 2016 were the same as the strains that had previously been isolated in this Rattus flavipectus plague focus. The c3 spacer present in the previously isolated strains was absent in the spacer arrays of the Ypc CRISPR loci of the strains isolated in 2016. The MLVA analysis using MLVA (14+12) showed that the strains isolated from the human plague case and host animal plague infection in 2016 in Yunnan displayed different molecular patterns than the strains that had previously been isolated from Yunnan and adjacent provinces.

Ten years of surveillance of the Yulong plague focus in China and the molecular typing and source tracing of the isolates

Plague, caused by Yersinia pestis, was classified as a reemerging infectious disease by the World Health Organization. The five human pneumonic plague cases in Yulong County in 2005 gave rise to the discovery of a Yulong plague focus in Yunnan province, China. Thereafter, continuous wild rodent plague (sylvatic plague) was identified as the main plague reservoir of this focus. In this study, the epizootics in Yulong focus were described, and three molecular typing methods, including the different region (DFR) analysis, clustered regularly interspaced short palindromic repeats (CRISPRs), and the multiple-locus variable number of tandem repeats (VNTR) analysis (MLVA) (14+12), were used for the molecular typing and source tracing of Y. pestis isolates in the Yulong plague focus. Simultaneously, several isolates from the vicinity of Yunnan were used as controls. The results showed that during the 10-year period from 2006 to 2016, an animal plague epidemic occurred in 6 of those years, and 5 villages underwent an animal plague epidemic within a 30-km² area of the Yulong plague focus. Searching for dead mice was the most effective monitoring method in this plague focus. No positive sample has been found in 6937 captured live rodents thus far, suggesting that the virulence of strains in the Yulong plague focus is stronger and the survival time of mice is shorter after infection. Strains from Lijiang, Sichuan and Tibet were of the same complex based on a typing analysis of DFR and CRISPR. The genetic relationship of Y. pestis illustrated by MLVA “14+12” demonstrates that Tibet and Sichuan strains evolved from the strains 1.IN2 (Qinghai, 1970 and Tibet, 1976), and Lijiang strains are closer to Batang strains (Batang County in Sichuan province, 2011, Himalaya marmot plague foci) in terms of genetic or phylogenic relationships. In conclusion, we have a deeper understanding of this new plague focus throughout this study, which provides a basis for effective prevention and control.

Plague is a type of zoonosis that is highly lethal to humans. The surveillance of animal hosts is critical for the prevention and control of plague. The Yulong plague focus is a newly discovered plague focus in China in recent years. The plague outbreak had attracted widespread attention because 5 people were infected in 2005, 2 of whom died. We have monitored the plague focus for a decade, and isolated strains and DNAs of Yersinia pestis were studied. The structure, origin and evolutionary trend of the Yulong plague focus were clarified, which provides a scientific basis for the effective prevention and control of human plague. This article also provides a set of paradigms for the systematic study of new plague foci, which is a perfect combination of traditional monitoring methods and modern research methods.

Genome rearrangements and phylogeny reconstruction in Yersinia pestis

Genome rearrangements have played an important role in the evolution of Yersinia pestis from its progenitor Yersinia pseudotuberculosis. Traditional phylogenetic trees for Y. pestis based on sequence comparison have short internal branches and low bootstrap supports as only a small number of nucleotide substitutions have occurred. On the other hand, even a small number of genome rearrangements may resolve topological ambiguities in a phylogenetic tree. We reconstructed phylogenetic trees based on genome rearrangements using several popular approaches such as Maximum likelihood for Gene Order and the Bayesian model of genome rearrangements by inversions. We also reconciled phylogenetic trees for each of the three CRISPR loci to obtain an integrated scenario of the CRISPR cassette evolution. Analysis of contradictions between the obtained evolutionary trees yielded numerous parallel inversions and gain/loss events. Our data indicate that an integrated analysis of sequence-based and inversion-based trees enhances the resolution of phylogenetic reconstruction. In contrast, reconstructions of strain relationships based on solely CRISPR loci may not be reliable, as the history is obscured by large deletions, obliterating the order of spacer gains. Similarly, numerous parallel gene losses preclude reconstruction of phylogeny based on gene content.

How old are bacterial pathogens?

Only few molecular studies have addressed the age of bacterial pathogens that infected humans before the beginnings of medical bacteriology, but these have provided dramatic insights. The global genetic diversity of Helicobacter pylori, which infects human stomachs, parallels that of its human host. The time to the most recent common ancestor (tMRCA) of these bacteria approximates that of anatomically modern humans, i.e. at least 100 000 years, after calibrating the evolutionary divergence within H. pylori against major ancient human migrations. Similarly, genomic reconstructions of Mycobacterium tuberculosis, the cause of tuberculosis, from ancient skeletons in South America and mummies in Hungary support estimates of less than 6000 years for the tMRCA of M. tuberculosis Finally, modern global patterns of genetic diversity and ancient DNA studies indicate that during the last 5000 years plague caused by Yersinia pestis has spread globally on multiple occasions from China and Central Asia. Such tMRCA estimates provide only lower bounds on the ages of bacterial pathogens, and additional studies are needed for realistic upper bounds on how long humans and animals have suffered from bacterial diseases.

Russian isolates enlarge the known geographic diversity of Francisella tularensis subsp. mediasiatica

Francisella tularensis, a small Gram-negative bacterium, is capable of infecting a wide range of animals, including humans, and causes a plague-like disease called tularemia—a highly contagious disease with a high mortality rate. Because of these characteristics, F. tularensis is considered a potential agent of biological terrorism. Currently, F. tularensis is divided into four subspecies, which differ in their virulence and geographic distribution. Two of them, subsp. tularensis (primarily found in North America) and subsp. holarctica (widespread across the Northern Hemisphere), are responsible for tularemia in humans. Subsp. novicida is almost avirulent in humans. The fourth subspecies, subsp. mediasiatica, is the least studied because of its limited distribution and impact in human health. It is found only in sparsely populated regions of Central Asia. In this report, we describe the first focus of naturally circulating F. tularensis subsp. mediasiatica in Russia. We isolated and characterized 18 strains of this subspecies in the Altai region. All strains were highly virulent in mice. The virulence of subsp. mediasiatica in a vaccinated mouse model is intermediate between that of subsp. tularensis and subsp. holarctica. Based on a multiple-locus variable number tandem repeat analysis (MLVA), we show that the Altaic population of F. tularensis subsp. mediasiatica is genetically distinct from the classical Central Asian population, and probably is endemic to Southern Siberia. We propose to subdivide the mediasiatica subspecies into three phylogeographic groups, M.I, M.II and M.III.

[An history of the CRISPR-Cas systems discovery]

From 1987 and during the following 20 years, a few research teams exploring bacteria and archea genome sequences uncover the prokaryotic adaptative immune system made of the CRISPR sequence and associated cas genes. First believed to be similar to the eukaryote RNA interference system, CRISPR-Cas turned out to be unique and of an amazing genetic complexity. The comparative studies of CRISPR arrays and of cas, and later of microbiotes metagenomes allowed to propose an evolution scenario for these systems. The results demonstrate the importance of a naturalistic approach, without a priori, for the understanding of living organisms.

Two Isoforms of Yersinia pestis Plasminogen Activator Pla: Intraspecies Distribution, Intrinsic Disorder Propensity, and Contribution to Virulence

It has been shown previously that several endemic Y. pestis isolates with limited virulence contained the I259 isoform of the outer membrane protease Pla, while the epidemic highly virulent strains possessed only the T259 Pla isoform. Our sequence analysis of the pla gene from 118 Y. pestis subsp. microtus strains revealed that the I259 isoform was present exclusively in the endemic strains providing a convictive evidence of more ancestral origin of this isoform. Analysis of the effects of the I259T polymorphism on the intrinsic disorder propensity of Pla revealed that the I259T mutation slightly increases the intrinsic disorder propensity of the C-terminal tail of Pla and makes this protein slightly more prone for disorder-based protein-protein interactions, suggesting that the T259 Pla could be functionally more active than the I259 Pla. This assumption was proven experimentally by assessing the coagulase and fibrinolytic activities of the two Pla isoforms in human plasma, as well as in a direct fluorometric assay with the Pla peptide substrate. The virulence testing of Pla-negative or expressing the I259 and T259 Pla isoforms Y. pestis subsp. microtus and subsp. pestis strains did not reveal any significant difference in LD₅₀ values and dose-dependent survival assays between them by using a subcutaneous route of challenge of mice and guinea pigs or intradermal challenge of mice. However, a significant decrease in time-to-death was observed in animals infected with the epidemic T259 Pla-producing strains as compared to the parent Pla-negative variants. Survival curves of the endemic I259 Pla⁺ strains fit between them, but significant difference in mean time to death post infection between the Pla⁻strains and their I259 Pla⁺ variants could be seen only in the isogenic set of subsp. pestis strains. These findings suggest an essential role for the outer membrane protease Pla evolution in Y. pestis bubonic infection exacerbation that is necessary for intensification of epidemic process from endemic natural focality with sporadic cases in men to rapidly expanding epizootics followed by human epidemic outbreaks, local epidemics or even pandemics.

Yersinia pestis Caf1 Protein: Effect of Sequence Polymorphism on Intrinsic Disorder Propensity, Serological Cross-Reactivity and Cross-Protectivity of Isoforms

Yersinia pestis Caf1 is a multifunctional protein responsible for antiphagocytic activity and is a key protective antigen. It is generally conserved between globally distributed Y. pestis strains, but Y. pestis subsp. microtus biovar caucasica strains circulating within populations of common voles in Georgia and Armenia were reported to carry a single substitution of alanine to serine. We investigated polymorphism of the Caf1 sequences among other Y. pestis subsp. microtus strains, which have a limited virulence in guinea pigs and in humans. Sequencing of caf1 genes from 119 Y. pestis strains belonging to different biovars within subsp. microtus showed that the Caf1 proteins exist in three isoforms, the global type Caf1NT1 (Ala48 Phe117), type Caf1NT2 (Ser48 Phe117) found in Transcaucasian-highland and Pre-Araks natural plague foci #4-7, and a novel Caf1NT3 type (Ala48 Val117) endemic in Dagestan-highland natural plague focus #39. Both minor types are the progenies of the global isoform. In this report, Caf1 polymorphism was analyzed by comparing predicted intrinsic disorder propensities and potential protein-protein interactivities of the three Caf1 isoforms. The analysis revealed that these properties of Caf1 protein are minimally affected by its polymorphism. All protein isoforms could be equally detected by an immunochromatography test for plague at the lowest protein concentration tested (1.0 ng/mL), which is the detection limit. When compared to the classic Caf1NT1 isoform, the endemic Caf1NT2 or Caf1NT3 had lower immunoreactivity in ELISA and lower indices of self- and cross-protection. Despite a visible reduction in cross-protection between all Caf1 isoforms, our data suggest that polymorphism in the caf1 gene may not allow the carriers of Caf1NT2 or Caf1NT3 variants escaping from the Caf1NT1-mediated immunity to plague in the case of a low-dose flea-borne infection.

CRISPR-Based Typing and Next-Generation Tracking Technologies

Bacteria occur ubiquitously in nature and are broadly relevant throughout the food supply chain, with diverse and variable tolerance levels depending on their origin, biological role, and impact on the quality and safety of the product as well as on the health of the consumer. With increasing knowledge of and accessibility to the microbial composition of our environments, food supply, and host-associated microbiota, our understanding of and appreciation for the ratio of beneficial to undesirable bacteria are rapidly evolving. Therefore, there is a need for tools and technologies that allow definite, accurate, and high-resolution identification and typing of various groups of bacteria that include beneficial microbes such as starter cultures and probiotics, innocuous commensals, and undesirable pathogens and spoilage organisms. During the transition from the current molecular biology-based PFGE (pulsed-field gel electrophoresis) gold standard to the increasingly accessible omics-level whole-genome sequencing (WGS) N-gen standard, high-resolution technologies such as CRISPR-based genotyping constitute practical and powerful alternatives that provide valuable insights into genome microevolution and evolutionary trajectories. Indeed, several studies have shown potential for CRISPR-based typing of industrial starter cultures, health-promoting probiotic strains, animal commensal species, and problematic pathogens. Emerging CRISPR-based typing methods open new avenues for high-resolution typing of a broad range of bacteria and constitute a practical means for rapid tracking of a diversity of food-associated microbes.

Taxonomy of Yersinia pestis

This chapter summarized the taxonomy and typing works of Yersinia pestis since it's firstly identified in Hong Kong in 1894. Phenotyping methods that based on phenotypic characteristics, including biotyping, serotyping, antibiogram analysis, bacteriocin typing, phage typing, and plasmid typing, were firstly applied in classification of Y. pestis in subspecies level. And then, with the advancement of molecular biological technology, the methods based on outer membrane protein profiles, fatty acid composition, and bacterial mass fingerprinting were also used to identify the populations within Y. pestis. However, Y. pestis is a highly homogenous species; therefore, the above typing methods could only provide low resolution, e.g., only one serotype and one phage type were observed for the whole species. Since the 1990s, molecular typing based on DNA variations, including single-nucleotide polymorphism, gene gain/loss, variable-number tandem repeats, clustered regularly interspaced short palindromic repeat, etc., was introduced and improved the resolution and robust of typing result. Especially in recent years, genotyping-based whole-genome-wide variations were successfully employed in Y. pestis, which built the "gold standard" of typing scheme of the species and could distinguish the samples under the strain level. The taxonomy and typing works leaved us enormous polymorphism data; therefore, a comprehensive fingerprint database of Y. pestis was needed to collect and standardize these data, for facilitating future works on evolution, plague surveillance and control, anti-bioterrorism, and microbial forensic researches.

Whole genome multilocus sequence typing as an epidemiologic tool for Yersinia pestis

Human plague is a severe and often fatal zoonotic disease caused by Yersinia pestis. For public health investigations of human cases, nonintensive whole genome molecular typing tools, capable of defining epidemiologic relationships, are advantageous. Whole genome multilocus sequence typing (wgMLST) is a recently developed methodology that simplifies genomic analyses by transforming millions of base pairs of sequence into character data for each gene. We sequenced 13 US Y. pestis isolates with known epidemiologic relationships. Sequences were assembled de novo, and multilocus sequence typing alleles were assigned by comparison against 3979 open reading frames from the reference strain CO92. Allele-based cluster analysis accurately grouped the 13 isolates, as well as 9 publicly available Y. pestis isolates, by their epidemiologic relationships. Our findings indicate wgMLST is a simplified, sensitive, and scalable tool for epidemiologic analysis of Y. pestis strains.

Nineteen Whole-Genome Assemblies of Yersinia pestis subsp. microtus, Including Representatives of Biovars caucasica, talassica, hissarica, altaica, xilingolensis, and ulegeica

The etiologic agent of plague, Yersinia pestis, includes two subspecies, of which Y. pestis subsp. microtus contains the strains that cause only occasional diseases in humans that are not accompanied by human-to-human transmission. Here, we report the draft genome sequences of 19 Y. pestis strains (across 6 biovars of Y. pestis subsp. microtus).

A review of methods for subtyping Yersinia pestis: From phenotypes to whole genome sequencing

Numerous subtyping methods have been applied to Yersinia pestis with varying success. Here, we review the various subtyping methods that have been applied to Y. pestis and their capacity for answering questions regarding the population genetics, phylogeography, and molecular epidemiology of this important human pathogen. Methods are evaluated in terms of expense, difficulty, transferability among laboratories, discriminatory power, usefulness for different study questions, and current applicability in light of the advent of whole genome sequencing.

Pneumonic plague outbreak, Northern Madagascar, 2011

Yersinia pestis, the causative agent of plague, is endemic to Madagascar, particularly to the central highlands. Although plague has not been previously reported in northern Madagascar, an outbreak of pneumonic plague occurred in this remote area in 2011. Over a 27-day period, 17 suspected, 2 presumptive, and 3 confirmed human cases were identified, and all 15 untreated 20 patients died. Molecular typing of Y. pestis isolated from 2 survivors and 5 Rattus rattus rat samples identified the Madagascar-specific 1.ORI3-k single-nucleotide polymorphism genotype and 4 clustered regularly interspaced short palindromic repeat patterns. This outbreak had a case-fatality rate of 100% for nontreated patients. The Y. pestis 1.ORI3-k single-nucleotide polymorphism genotype might cause larger epidemics. Multidrug-resistant strains and persistence of the pathogen in natural foci near human settlements pose severe risks to populations in plague-endemic regions and require outbreak response strategies.

Diverse Genotypes of Yersinia pestis Caused Plague in Madagascar in 2007

BACKGROUND

Yersinia pestis is the causative agent of human plague and is endemic in various African, Asian and American countries. In Madagascar, the disease represents a significant public health problem with hundreds of human cases a year. Unfortunately, poor infrastructure makes outbreak investigations challenging.

METHODOLOGY/PRINCIPAL FINDINGS

DNA was extracted directly from 93 clinical samples from patients with a clinical diagnosis of plague in Madagascar in 2007. The extracted DNAs were then genotyped using three molecular genotyping methods, including, single nucleotide polymorphism (SNP) typing, multi-locus variable-number tandem repeat analysis (MLVA), and Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) analysis. These methods provided increasing resolution, respectively. The results of these analyses revealed that, in 2007, ten molecular groups, two newly described here and eight previously identified, were responsible for causing human plague in geographically distinct areas of Madagascar.

CONCLUSIONS/SIGNIFICANCE

Plague in Madagascar is caused by numerous distinct types of Y. pestis. Genotyping method choice should be based upon the discriminatory power needed, expense, and available data for any desired comparisons. We conclude that genotyping should be a standard tool used in epidemiological investigations of plague outbreaks.

A novel multiple locus variable number of tandem repeat (VNTR) analysis (MLVA) method for Propionibacterium acnes

Propionibacterium acnes plays a central role in the pathogenesis of acne and is responsible for severe opportunistic infections. Numerous typing schemes have been developed that allow the identification of phylotypes, but they are often insufficient to differentiate subtypes. To better understand the genetic diversity of this species and to perform epidemiological analyses, high throughput discriminant genotyping techniques are needed. Here we describe the development of a multiple locus variable number of tandem repeats (VNTR) analysis (MLVA) method. Thirteen VNTRs were identified in the genome of P. acnes and were used to genotype a collection of clinical isolates. In addition, publically available sequencing data for 102 genomes were analyzed in silico, providing an MLVA genotype. The clustering of MLVA data was in perfect congruence with whole genome based clustering. Analysis of the clustered regularly interspaced short palindromic repeat (CRISPR) element uncovered new spacers, a supplementary source of genotypic information. The present MLVA13 scheme and associated internet database represents a first line genotyping assay to investigate large number of isolates. Particular strains may then be submitted to full genome sequencing in order to better analyze their pathogenic potential.

Characterization and evolution of Salmonella CRISPR-Cas systems

Prokaryotic CRISPR-Cas (clustered regularly interspaced short palindromic repeats and CRISPR-associated genes) systems provide adaptive immunity from invasive genetic elements and encompass three essential features: (i) cas genes, (ii) a CRISPR array composed of spacers and direct repeats and (iii) an AT-rich leader sequence upstream of the array. We performed in-depth sequence analysis of the CRISPR-Cas systems in >600 Salmonella, representing four clinically prevalent serovars. Each CRISPR-Cas feature is extremely conserved in the Salmonella, and the CRISPR1 locus is more highly conserved than CRISPR2. Array composition is serovar-specific, although no convincing evidence of recent spacer acquisition against exogenous nucleic acids exists. Only 12 % of spacers match phage and plasmid sequences and self-targeting spacers are associated with direct repeat variants. High nucleotide identity (>99.9 %) exists across the cas operon among isolates of a single serovar and in some cases this conservation extends across divergent serovars. These observations reflect historical CRISPR-Cas immune activity, showing that this locus has ceased undergoing adaptive events. Intriguingly, the high level of conservation across divergent serovars shows that the genetic integrity of these inactive loci is maintained over time, contrasting with the canonical view that inactive CRISPR loci degenerate over time. This thorough characterization of Salmonella CRISPR-Cas systems presents new insights into Salmonella CRISPR evolution, particularly with respect to cas gene conservation, leader sequences, organization of direct repeats and protospacer matches. Collectively, our data suggest that Salmonella CRISPR-Cas systems are no longer immunogenic; rather, their impressive conservation indicates they may have an alternative function in Salmonella.

Staphylococcus aureus from 152 cases of bovine, ovine and caprine mastitis investigated by Multiple-locus variable number of tandem repeat analysis (MLVA)

Staphylococcus aureus is one of the main etiological agents of mastitis in ruminants. In the present retrospective study, we evaluated the potential interest of a previously described automated multiple loci Variable Number of Tandem Repeats (VNTR) Assay (MLVA) comprising 16 loci as a first line tool to investigate the population structure of S. aureus from mastitis. We determined the genetic diversity of S. aureus strains from cases of clinical and subclinical mastitis in dairy cattle (n = 118, of which 16 were methicillin-resistant), sheep (n = 18) and goats (n = 16). The 152 strains could be subdivided into 115 MLVA genotypes (including 14 genotypes for the ovine strains and 15 genotypes for the caprine strains). This corresponds to a discriminatory index (D) value of 0.9936. Comparison with published MLVA data obtained using the same protocol applied to strains from diverse human and animal origins revealed a low number (8.5%) of human-related MLVA genotypes among the present collection. Eighteen percent of the S. aureus mastitis collection belonged to clonal complexes apparently not associated with other pathological conditions. Some of them displayed a relatively low level of diversity in agreement with a restricted ecological niche. These findings provide arguments suggesting that specific S. aureus lineages particularly adapted to ruminant mammary glands have emerged and that MLVA is a convenient tool to provide a broad overview of the population, owing to the availability via internet of databases compiling published MLVA genotypes.

Dynamics of CRISPR loci in microevolutionary process of Yersinia pestis strains

The potential use of CRISPR loci genotyping to elucidate population dynamics and microevolution of 146 Yersinia pestis strains from different biovars and locations was investigated in this work. The majority of strains from the Orientalis biovar presented specific spacer arrays, allowing for the establishment of a CRISPR signature for their respective isolates. Twenty-one new spacers were found in the Y. pestis strains from plague foci in Brazil. Ninety-three (64%) strains were grouped in the G1 genotype, whereas the others were distributed in 35 genotypes. This study allowed observing a microevolutionary process in a group of Y. pestis isolated from Brazil. We also identified specific genotypes of Y. pestis that were important for the establishment of the bacteria in plague foci in Brazil. The data have provided supporting evidence for the diversity and dynamics of CRISPR loci present in the genome of Y. pestis strains from plague foci in Brazil.

Evolutionary causes and consequences of diversified CRISPR immune profiles in natural populations

Host-pathogen co-evolution is a significant force which shapes the ecology and evolution of all types of organisms, and such interactions are driven by resistance and immunity mechanisms of the host. Diversity of resistance and immunity can affect the co-evolutionary trajectory of both host and pathogen. The microbial CRISPR (clustered regularly interspaced short palindromic repeats)-Cas (CRISPR-associated) system is one host immunity mechanism which offers a tractable model for examining the dynamics of diversity in an immune system. In the present article, we review CRISPR variation observed in a variety of natural populations, examine the forces which can push CRISPRs towards high or low diversity, and investigate the consequences of various levels of diversity on microbial populations.

Yersinia pseudotuberculosis ST42 (O:1) Strain Misidentified as Yersinia pestis by Mass Spectrometry Analysis

We report here the draft sequence of strain CEB14_0017, alias HIAD_DUP, recovered from a human patient and initially identified as Yersinia pestis by mass spectrometry analysis. Genotyping based on tandem repeat polymorphism assigned the strain to Yersinia pseudotuberculosis sequence type 42 (ST42). The total assembly length is 4,894,739 bp.

The role of CRISPR-Cas systems in virulence of pathogenic bacteria

Clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated (Cas) genes are present in many bacterial and archaeal genomes. Since the discovery of the typical CRISPR loci in the 1980s, well before their physiological role was revealed, their variable sequences have been used as a complementary typing tool in diagnostic, epidemiologic, and evolutionary analyses of prokaryotic strains. The discovery that CRISPR spacers are often identical to sequence fragments of mobile genetic elements was a major breakthrough that eventually led to the elucidation of CRISPR-Cas as an adaptive immunity system. Key elements of this unique prokaryotic defense system are small CRISPR RNAs that guide nucleases to complementary target nucleic acids of invading viruses and plasmids, generally followed by the degradation of the invader. In addition, several recent studies have pointed at direct links of CRISPR-Cas to regulation of a range of stress-related phenomena. An interesting example concerns a pathogenic bacterium that possesses a CRISPR-associated ribonucleoprotein complex that may play a dual role in defense and/or virulence. In this review, we describe recently reported cases of potential involvement of CRISPR-Cas systems in bacterial stress responses in general and bacterial virulence in particular.

CRISPRs: molecular signatures used for pathogen subtyping

Rapid and accurate strain identification is paramount in the battle against microbial outbreaks, and several subtyping approaches have been developed. One such method uses clustered regular interspaced short palindromic repeats (CRISPRs), DNA repeat elements that are present in approximately half of all bacteria. Though their signature function is as an adaptive immune system against invading DNA such as bacteriophages and plasmids, CRISPRs also provide an excellent framework for pathogen tracking and evolutionary studies. Analysis of the spacer DNA sequences that reside between the repeats has been tremendously useful for bacterial subtyping during molecular epidemiological investigations. Subtyping, or strain identification, using CRISPRs has been employed in diverse Gram-positive and Gram-negative bacteria, including Mycobacterium tuberculosis, Salmonella enterica, and the plant pathogen Erwinia amylovora. This review discusses the several ways in which CRISPR sequences are exploited for subtyping. This includes the well-established spoligotyping methodologies that have been used for 2 decades to type Mycobacterium species, as well as in-depth consideration of newer, higher-throughput CRISPR-based protocols.

Features of Variable Number of Tandem Repeats in Yersinia pestis and the Development of a Hierarchical Genotyping Scheme

Background

Variable number of tandem repeats (VNTRs) that are widely distributed in the genome of Yersinia pestis proved to be useful markers for the genotyping and source-tracing of this notorious pathogen. In this study, we probed into the features of VNTRs in the Y. pestis genome and developed a simple hierarchical genotyping system based on optimized VNTR loci.

Methodology/Principal Findings

Capillary electrophoresis was used in this study for multi-locus VNTR analysis (MLVA) in 956 Y. pestis strains. The general features and genetic diversities of 88 VNTR loci in Y. pestis were analyzed with BioNumerics, and a “14+12” loci-based hierarchical genotyping system, which is compatible with single nucleotide polymorphism-based phylogenic analysis, was established.

Conclusions/Significance

Appropriate selection of target loci reduces the impact of homoplasies caused by the rapid mutation rates of VNTR loci. The optimized “14+12” loci are highly discriminative in genotyping and source-tracing Y. pestis for molecular epidemiological or microbial forensic investigations with less time and lower cost. An MLVA genotyping datasets of representative strains will improve future research on the source-tracing and microevolution of Y. pestis.

Yersinia pestis DNA from skeletal remains from the 6(th) century AD reveals insights into Justinianic Plague

Yersinia pestis, the etiologic agent of the disease plague, has been implicated in three historical pandemics. These include the third pandemic of the 19^th and 20^th centuries, during which plague was spread around the world, and the second pandemic of the 14^th–17^th centuries, which included the infamous epidemic known as the Black Death. Previous studies have confirmed that Y. pestis caused these two more recent pandemics. However, a highly spirited debate still continues as to whether Y. pestis caused the so-called Justinianic Plague of the 6^th–8^th centuries AD. By analyzing ancient DNA in two independent ancient DNA laboratories, we confirmed unambiguously the presence of Y. pestis DNA in human skeletal remains from an Early Medieval cemetery. In addition, we narrowed the phylogenetic position of the responsible strain down to major branch 0 on the Y. pestis phylogeny, specifically between nodes N03 and N05. Our findings confirm that Y. pestis was responsible for the Justinianic Plague, which should end the controversy regarding the etiology of this pandemic. The first genotype of a Y. pestis strain that caused the Late Antique plague provides important information about the history of the plague bacillus and suggests that the first pandemic also originated in Asia, similar to the other two plague pandemics.

Plague is a notorious and fatal human disease caused by the bacterium Yersinia pestis that is endemic in many countries around the world. Three of the most devastating pandemics in human history have been associated with plague. The second pandemic originated in central Asia and peaked in Europe between 1348 and 1350 (a period of time known as the Black Death). The third pandemic began in the Yunnan province of China in the mid-1850s and subsequently spread to Africa, the Americas, Australia, Europe, and other parts of Asia. The second and third pandemics are well documented and scientifically proven. However, the first pandemic, which began in the 6^th century and is also known as Justinianic Plague, is still a matter of controversy. Recently it has been suggested that Justinian's plague was not caused by Y. pestis. We detected Y. pestis DNA in samples obtained from multiple 6^th century skeletons from Germany. This confirms that Justinianic Plague crossed the Alps and affected local populations there, including current day Bavaria. Furthermore, we used DNA fingerprinting approaches to determine Asia as the likely geographic origin for these strains.

Draft Genome Sequences of Five Yersinia pseudotuberculosis ST19 Isolates and One Isolate Variant

We report the first draft genome sequences of five Yersinia pseudotuberculosis isolates of sequence type (ST) 19 and of a variant from one of the five isolates. The total length of assemblies ranged from 4,226,485 bp to 4,274,148 bp, including between 3,808 and 3,843 predicted coding sequences.

A decade of plague in Mahajanga, Madagascar: insights into the global maritime spread of pandemic plague

A cluster of human plague cases occurred in the seaport city of Mahajanga, Madagascar, from 1991 to 1999 following 62 years with no evidence of plague, which offered insights into plague pathogen dynamics in an urban environment. We analyzed a set of 44 Mahajanga isolates from this 9-year outbreak, as well as an additional 218 Malagasy isolates from the highland foci. We sequenced the genomes of four Mahajanga strains, performed whole-genome sequence single-nucleotide polymorphism (SNP) discovery on those strains, screened the discovered SNPs, and performed a high-resolution 43-locus multilocus variable-number tandem-repeat analysis of the isolate panel. Twenty-two new SNPs were identified and defined a new phylogenetic lineage among the Malagasy isolates. Phylogeographic analysis suggests that the Mahajanga lineage likely originated in the Ambositra district in the highlands, spread throughout the northern central highlands, and was then introduced into and became transiently established in Mahajanga. Although multiple transfers between the central highlands and Mahajanga occurred, there was a locally differentiating and dominant subpopulation that was primarily responsible for the 1991-to-1999 Mahajanga outbreaks. Phylotemporal analysis of this Mahajanga subpopulation revealed a cycling pattern of diversity generation and loss that occurred during and after each outbreak. This pattern is consistent with severe interseasonal genetic bottlenecks along with large seasonal population expansions. The ultimate extinction of plague pathogens in Mahajanga suggests that, in this environment, the plague pathogen niche is tenuous at best. However, the temporary large pathogen population expansion provides the means for plague pathogens to disperse and become ecologically established in more suitable nonurban environments.

Maritime spread of plague led to the global dissemination of this disease and affected the course of human history. Multiple historical plague waves resulted in massive human mortalities in three classical plague pandemics: Justinian (6th and 7th centuries), Middle Ages (14th to 17th centuries), and third (mid-1800s to the present). Key to these events was the pathogen’s entry into new lands by “plague ships” via seaport cities. Although initial disease outbreaks in ports were common, they were almost never sustained for long and plague pathogens survived only if they could become established in ecologically suitable habitats. Although plague pathogens’ ability to invade port cities has been essential for intercontinental spread, these regions have not proven to be a suitable long-term niche. The disease dynamics in port cities such as Mahajanga are thus critical to plague pathogen amplification and dispersal into new suitable ecological niches for the observed global long-term maintenance of plague pathogens.

The combination of CRISPR-MVLST and PFGE provides increased discriminatory power for differentiating human clinical isolates of Salmonella enterica subsp. enterica serovar Enteritidis

Salmonella enterica subsp. enterica serovar Enteritidis (S. Enteritidis) is a major cause of foodborne salmonellosis. Rapid, efficient and accurate methods for identification are required to track specific strains of S. Enteritidis during outbreaks of human salmonellosis. By exploiting the hypervariable nature of virulence genes and Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs), we previously developed a powerful sequence-based subtyping approach, designated CRISPR-MVLST. To substantiate the applicability of CRISPR-MVLST, we analyzed a broad set of S. Enteritidis isolates collected over a six-year period. Among 141 isolates we defined 22 Enteritidis Sequence Types (ESTs), the majority of which were novel. Notably, strains exhibiting the common PFGE pattern, JEGX01.0004 (characteristic of ∼40% of S. Enteritidis isolates in the United States), were separated into twelve distinct sequence types. Conversely, isolates of EST4, the most predominant EST we observed, comprised eight different PFGE patterns. Importantly, we showed that some genotypes that were previously associated with the food supply chain at the farm level have now been identified in clinical samples. CRISPR sequence data shows subtle but distinct differences among different alleles of S. Enteritidis, suggesting that evolution of these loci occurs vertically, as opposed to previously reported evolution by spacer acquisition in other bacteria.

Lipopolysaccharide of Yersinia pestis, the Cause of Plague: Structure, Genetics, Biological Properties

The present review summarizes data pertaining to the composition and structure of the carbohydrate moiety (core oligosaccharide) and lipid component (lipid A) of the various forms of lipopolysaccharide (LPS), one of the major pathogenicity factors ofYersinia pestis, the cause of plague. The review addresses the functions and the biological significance of genes for the biosynthesis of LPS, as well as the biological properties of LPS in strains from various intraspecies groups ofY. pestis and their mutants, including the contribution of LPS to the resistance of bacteria to factors of the innate immunity of both insect-vectors and mammal-hosts. Special attention is paid to temperature-dependent variations in the LPS structure, their genetic control and roles in the pathogenesis of plague. The evolutionary aspect is considered based on a comparison of the structure and genetics of the LPS ofY. pestis and other enteric bacteria, including otherYersinia species. The prospects of development of live plague vaccines created on the basis ofY. pestis strains with the genetically modified LPS are discussed.