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

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.

Evolutionary selection of biofilm-mediated extended phenotypes in Yersinia pestis in response to a fluctuating environment

Yersinia pestis is transmitted from fleas to rodents when the bacterium develops an extensive biofilm in the foregut of a flea, starving it into a feeding frenzy, or, alternatively, during a brief period directly after feeding on a bacteremic host. These two transmission modes are in a trade-off regulated by the amount of biofilm produced by the bacterium. Here by investigating 446 global isolated Y. pestis genomes, including 78 newly sequenced isolates sampled over 40 years from a plague focus in China, we provide evidence for strong selection pressures on the RNA polymerase ω-subunit encoding gene rpoZ. We demonstrate that rpoZ variants have an increased rate of biofilm production in vitro, and that they evolve in the ecosystem during colder and drier periods. Our results support the notion that the bacterium is constantly adapting-through extended phenotype changes in the fleas-in response to climate-driven changes in the niche.

Natural Mega-Focus of Yersinia pestis Main Subspecies, Antique Biovar, Phylogenetic Line 4.ANT in Gorny Altai

Objective of this study was to investigate the areal of Yersinia pestis strains of antique biovar, main subspecies, phylogenetic line 4.ANT, as well as to establish the borders and spatial structure of mega-focus in the territory of Gorny Altai. Materials and methods. Complex comprehensive analysis of the properties in Yersinia pestis strains of the main subspecies, circulating in Gorny Altai has been conducted. 33 out of them, isolated predominantly in 2012–2017 have been sequenced. Whole genome SNP analysis and search of marker SNPs was performed using Wombac 2.0 software package. Tree diagram was built applying Maximum Likelihood algorithm, using PhyML 3.1 software on the basis of HKY85 model. Results and conclusions. Based on the results of whole genome SNPs analysis of 33 endemic strains and creation of the tree diagram of Y. pestis strains, the presence of natural mega-focus of Y. pestis belonging to the main subspecies, antique biovar, phylogenetic line 4.ANT has been substantiated. Epizootic manifestations on multiple local areas, characterized by persistent autonomous nature of plague manifestations, are registered on an annual basis. Within the boundaries of the areal of Yersinia pestis main subspecies, antique biovar, phylogenetic line 4.ANT, existence of joint natural foci of Yersinia pestis belonging to non-main subspecies ssp. altaica and ulegeica is established. Location of natural foci of the main and non-main subspecies of Yersinia pestis in different altitudinal belts of the Altai Mountains Range on the whole provides for observed multi-host and multi-vector feature of epizootic manifestations. For the first time ever, the data on the areal of the main subspecies of plague microbe are used for setting the boundaries of its natural focus.

Genetic diversity and spatial-temporal distribution of Yersinia pestis in Qinghai Plateau, China

BACKGROUND

Plague, caused by the bacterium Yersinia pestis, is a highly infectious, zoonotic disease. Hundreds of human plague cases are reported across the world annually. Qinghai Plateau is one of the most severely affected plague regions in China, with more than 240 fatal cases of Y. pestis in the last 60 years. Conventional epidemiologic analysis has effectively guided the prevention and control of local plague transmission; however, molecular genetic analysis is more effective for investigating population diversity and transmission. In this report, we employed different genetic markers to analyze the population structure of Y. pestis in Qinghai Plateau.

METHODOLOGY/PRINCIPAL FINDING

We employed a two-step hierarchical strategy to analyze the phylogeny of 102 Qinghai Plateau isolates of Y. pestis, collected between 1954 and 2011. First, we defined the genealogy of Y. pestis by constructed minimum spanning tree based on 25 key SNPs. Seven groups were identified, with group 1.IN2 being identified as the dominant population. Second, two methods, MLVA and CRISPR, were applied to examine the phylogenetic detail of group 1.IN2, which was further divided into three subgroups. Subgroups of 1.IN2 revealed a clear geographic cluster, possibly associated with interaction between bacteriophage and Y. pestis. More recently, Y. pestis populations appear to have shifted from the east toward the center and west of Qinghai Plateau. This shift could be related to destruction of the local niche of the original plague focus through human activities. Additionally, we found that the abundance and relative proportion of 1.IN2 subgroups varied by decade and might be responsible for the fluctuations of plague epidemics in Qinghai Plateau.

CONCLUSION/SIGNIFICANCE

Molecular genotyping methods provided us with detailed information on population diversity and the spatial-temporal distribution of dominant populations of Y. pestis, which will facilitate future surveillance, prevention, and control of plague in Qinghai Plateau.

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.

Variant site strain typer (VaST): efficient strain typing using a minimal number of variant genomic sites

BACKGROUND

Targeted PCR amplicon sequencing (TAS) techniques provide a sensitive, scalable, and cost-effective way to query and identify closely related bacterial species and strains. Typically, this is accomplished by targeting housekeeping genes that provide resolution down to the family, genera, and sometimes species level. Unfortunately, this level of resolution is not sufficient in many applications where strain-level identification of bacteria is required (biodefense, forensics, clinical diagnostics, and outbreak investigations). Adding more genomic targets will increase the resolution, but the challenge is identifying the appropriate targets. VaST was developed to address this challenge by finding the minimum number of targets that, in combination, achieve maximum strain-level resolution for any strain complex. The final combination of target regions identified by the algorithm produce a unique haplotype for each strain which can be used as a fingerprint for identifying unknown samples in a TAS assay. VaST ensures that the targets have conserved primer regions so that the targets can be amplified in all of the known strains and it also favors the inclusion of targets with basal variants which makes the set more robust when identifying previously unseen strains.

RESULTS

We analyzed VaST's performance using a number of different pathogenic species that are relevant to human disease outbreaks and biodefense. The number of targets required to achieve full resolution ranged from 20 to 88% fewer sites than what would be required in the worst case and most of the resolution is achieved within the first 20 targets. We computationally and experimentally validated one of the VaST panels and found that the targets led to accurate phylogenetic placement of strains, even when the strains were not a part of the original panel design.

CONCLUSIONS

VaST is an open source software that, when provided a set of variant sites, can find the minimum number of sites that will provide maximum resolution of a strain complex, and it has many different run-time options that can accommodate a wide range of applications. VaST can be an effective tool in the design of strain identification panels that, when combined with TAS technologies, offer an efficient and inexpensive strain typing protocol.

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.

Phenotypic and Molecular Genetic Characteristics of Yersinia pestis at an Emerging Natural Plague Focus, Junggar Basin, China

The 15th natural plague focus in China, the Junggar Basin plague focus, is located near an important communication route connecting China and Central Asia and was discovered after 2005. To characterize the phenotypic and genetic diversity of the Yersinia pestis population in this newly established focus, we collected 25 Y. pestis strains from six counties across Junggar Basin in 2005-2006, and determined their biochemical features and genotypes based on multiple-locus variable number of tandem repeats analysis and clustered regularly interspaced short palindromic repeats analysis. We inferred the phylogenetic positions and possible sources of the Junggar strains by comparing their genotypes with the genetic diversity for known representative Y. pestis strains. Our results indicate that the major genotype of Junggar strains belongs to 2.MED1, a lineage of biovar Medievalis with identical biochemical characters and high virulence in mice. Although share a similar ecology, the 2.MED1 in Junggar Basin are not descended from known strains in the neighboring Central Asian Desert plague foci. Therefore, the emergence of the Junggar Basin plague focus is not attributable to the recent clonal spread of Y. pestis from Central Asia. We also identified two distinct minor genotypes in Junggar Basin, one of which clusters genetically with the 0.ANT1 strains of the Tianshan Mountain natural plague focus and another belongs to a 1.IN lineage not previously reported. Our study clarifies the phenotypic and genetic characters of Junggar Y. pestis strains. These findings extend our knowledge of the population diversity of Y. pestis and will facilitate future plague surveillance and prevention in Junggar Basin and adjacent regions.

Capillary electrophoresis based on nucleic acid detection for diagnosing human infectious disease

Rapid transmission, high morbidity, and mortality are the features of human infectious diseases caused by microorganisms, such as bacteria, fungi, and viruses. These diseases may lead within a short period of time to great personal and property losses, especially in regions where sanitation is poor. Thus, rapid diagnoses are vital for the prevention and therapeutic intervention of human infectious diseases. Several conventional methods are often used to diagnose infectious diseases, e.g. methods based on cultures or morphology, or biochemical tests based on metabonomics. Although traditional methods are considered gold standards and are used most frequently, they are laborious, time consuming, and tedious and cannot meet the demand for rapid diagnoses. Disease diagnosis using capillary electrophoresis methods has the advantages of high efficiency, high throughput, and high speed, and coupled with the different nucleic acid detection strategies overcomes the drawbacks of traditional identification methods, precluding many types of false positive and negative results. Therefore, this review focuses on the application of capillary electrophoresis based on nucleic detection to the diagnosis of human infectious diseases, and offers an introduction to the limitations, advantages, and future developments of this approach.

Reannotation of Yersinia pestis Strain 91001 Based on Omics Data

Yersinia pestis is among the most dangerous human pathogens, and systematic research of this pathogen is important in bacterial pathogenomics research. To fully interpret the biological functions, physiological characteristics, and pathogenesis of Y. pestis, a comprehensive annotation of its entire genome is necessary. The emergence of omics-based research has brought new opportunities to better annotate the genome of this pathogen. Here, the complete genome of Y. pestis strain 91001 was reannotated using genomics and proteogenomics data. One hundred and thirty-seven unreliable coding sequences were removed, and 41 homologous genes were relocated with their translational initiation sites, while the functions of seven pseudogenes and 392 hypothetical genes were revised. Moreover, annotations of noncoding RNAs, repeat sequences, and transposable elements have also been incorporated. The reannotated results are freely available at http://tody.bmi.ac.cn.

Two Distinct Yersinia pestis Populations Causing Plague among Humans in the West Nile Region of Uganda

Background

Plague is a life-threatening disease caused by the bacterium, Yersinia pestis. Since the 1990s, Africa has accounted for the majority of reported human cases. In Uganda, plague cases occur in the West Nile region, near the border with Democratic Republic of Congo. Despite the ongoing risk of contracting plague in this region, little is known about Y. pestis genotypes causing human disease.

Methodology/Principal Findings

During January 2004–December 2012, 1,092 suspect human plague cases were recorded in the West Nile region of Uganda. Sixty-one cases were culture-confirmed. Recovered Y. pestis isolates were analyzed using three typing methods, single nucleotide polymorphisms (SNPs), pulsed field gel electrophoresis (PFGE), and multiple variable number of tandem repeat analysis (MLVA) and subpopulations analyzed in the context of associated geographic, temporal, and clinical data for source patients. All three methods separated the 61 isolates into two distinct 1.ANT lineages, which persisted throughout the 9 year period and were associated with differences in elevation and geographic distribution.

Conclusions/Significance

We demonstrate that human cases of plague in the West Nile region of Uganda are caused by two distinct 1.ANT genetic subpopulations. Notably, all three typing methods used, SNPs, PFGE, and MLVA, identified the two genetic subpopulations, despite recognizing different mutation types in the Y. pestis genome. The geographic and elevation differences between the two subpopulations is suggestive of their maintenance in highly localized enzootic cycles, potentially with differing vector-host community composition. This improved understanding of Y. pestis subpopulations in the West Nile region will be useful for identifying ecologic and environmental factors associated with elevated plague risk.

Plague, a severe and often fatal zoonotic disease, is caused by the bacterium Yersinia pestis. Currently, the majority of human cases have been reported from resource limited areas of Africa, where the proximity to commensal rats and other small mammals increases the likelihood for human contact with infected animals or their fleas. Over a 9 year time period, >1000 suspect cases were recorded in the West Nile region of Uganda within the districts of Arua and Zombo. Culture-confirmed cases were shown by three independent typing methods to be due to two distinct 1.ANT genetic subpopulations of Y. pestis. The two genetic subpopulations persisted throughout the 9 year time period, consistent with their ongoing maintenance in local enzootic cycles. Additionally, the two subpopulations were found to differ with respect to geographic location and elevation, with SNP Group 1 strains being found further north and at lower elevations as compared to SNP Group 2. The relative independence of the two Y. pestis subpopulations is suggestive of their maintenance in distinct foci involving enzootic cycles with differing vector-host community composition.

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.

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.

Strain discrimination of Yersinia pestis using a SERS-based electrochemically driven melting curve analysis of variable number tandem repeat sequences

Strain discrimination within genetically highly similar bacteria is critical for epidemiological studies and forensic applications. An electrochemically driven melting curve analysis monitored by SERS has been utilised to reliably discriminate strains of the bacterial pathogen Yersinia pestis, the causative agent of plague. DNA amplicons containing Variable Number Tandem Repeats (VNTRs) were generated from three strains of Y. pestis: CO92, Harbin 35 and Kim. These amplicons contained a 10 base pair VNTR repeated 6, 5, and 4 times in CO92, Harbin 35 and Kim respectively. The assay also included a blocker oligonucleotide comprising 3 repeats of the 10-mer VNTR sequence. The use of the blocker reduced the effective length of the target sequence available to bind to the surface bound probe and significantly improved the sensitivity of the discrimination. The results were consistent during three replicates that were carried out on different days, using different batches of PCR product and different SERS sphere segment void (SSV) substrate. This methodology which combines low cost, speed and sensitivity is a promising alternative to the time consuming current electrophoretic methods.

Genetic variations of live attenuated plague vaccine strains (Yersinia pestis EV76 lineage) during laboratory passages in different countries

Plague, one of the most devastating infectious diseases in human history, is caused by the bacterial species Yersinia pestis. A live attenuated Y. pestis strain (EV76) has been widely used as a plague vaccine in various countries around the world. Here we compared the whole genome sequence of an EV76 strain used in China (EV76-CN) with the genomes of Y. pestis wild isolates to identify genetic variations specific to the EV76 lineage. We identified 6 SNPs and 6 Indels (insertions and deletions) differentiating EV76-CN from its counterparts. Then, we screened these polymorphic sites in 28 other strains of EV76 lineage that were stored in different countries. Based on the profiles of SNPs and Indels, we reconstructed the parsimonious dissemination history of EV76 lineage. This analysis revealed that there have been at least three independent imports of EV76 strains into China. Additionally, we observed that the pyrE gene is a mutation hotspot in EV76 lineages. The fine comparison results based on whole genome sequence in this study provide better understanding of the effects of laboratory passages on the accumulation of genetic polymorphisms in plague vaccine strains. These variations identified here will also be helpful in discriminating different EV76 derivatives.

Two-step source tracing strategy of Yersinia pestis and its historical epidemiology in a specific region

Source tracing of pathogens is critical for the control and prevention of infectious diseases. Genome sequencing by high throughput technologies is currently feasible and popular, leading to the burst of deciphered bacterial genome sequences. Utilizing the flooding genomic data for source tracing of pathogens in outbreaks is promising, and challenging as well. Here, we employed Yersinia pestis genomes from a plague outbreak at Xinghai county of China in 2009 as an example, to develop a simple two-step strategy for rapid source tracing of the outbreak. The first step was to define the phylogenetic position of the outbreak strains in a whole species tree, and the next step was to provide a detailed relationship across the outbreak strains and their suspected relatives. Through this strategy, we observed that the Xinghai plague outbreak was caused by Y. pestis that circulated in the local plague focus, where the majority of historical plague epidemics in the Qinghai-Tibet Plateau may originate from. The analytical strategy developed here will be of great help in fighting against the outbreaks of emerging infectious diseases, by pinpointing the source of pathogens rapidly with genomic epidemiological data and microbial forensics information.