Differential plague-transmission dynamics determine Yersinia pestis population genetic structure on local, regional, and global scales

The Causes for Genomic Instability and How to Try and Reduce Them Through Rational Design of Synthetic DNA

Genetic engineering has revolutionized our ability to manipulate DNA and engineer organisms for various applications. However, this approach can lead to genomic instability, which can result in unwanted effects such as toxicity, mutagenesis, and reduced productivity. To overcome these challenges, smart design of synthetic DNA has emerged as a promising solution. By taking into consideration the intricate relationships between gene expression and cellular metabolism, researchers can design synthetic constructs that minimize metabolic stress on the host cell, reduce mutagenesis, and increase protein yield. In this chapter, we summarize the main challenges of genomic instability in genetic engineering and address the dangers of unknowingly incorporating genomically unstable sequences in synthetic DNA. We also demonstrate the instability of those sequences by the fact that they are selected against conserved sequences in nature. We highlight the benefits of using ESO, a tool for the rational design of DNA for avoiding genetically unstable sequences, and also summarize the main principles and working parameters of the software that allow maximizing its benefits and impact.

<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.

Epidemiological Characteristics of Human and Animal Plague in Yunnan Province, China, 1950 to 2020

This study analyzed the epidemiological characteristics of 3,464 human plague cases and the distribution pattern of 4,968 Yersinia pestis isolates from humans, hosts, and vector insects from 1950 to 2020 among two natural plague foci in Yunnan Province, China. These foci include the Rattus flavipectus plague focus of the Yunnan, Guangdong, and Fujian provinces and the Apodemus chevrieri-Eothenomys miletus plague focus of the highlands of northwestern Yunnan Province. The case fatality rate for plague in humans was 18.39% (637/3,464), and the total isolation rate of Y. pestis was 0.17% (4,968/2,975,288). Despite that the frequency of human cases declined rapidly, the animal plague fluctuated greatly, alternating between activity and inactivity in these foci. The tendency among human cases can be divided into 4 stages, 1950 to 1955, 1956 to 1989, 1990 to 2005, and 2006 to 2020. Bubonic plague accounted for the majority of cases in Yunnan, where pneumonic and septicemic plague rarely occurred. The natural plague foci have been in a relatively active state due to the stability of local ecology. Dense human population and frequent contact with host animals contribute to the high risk of human infection. This study systematically analyzed the epidemic pattern of human plague and the distribution characteristics of Y. pestis in the natural plague foci in Yunnan, providing a scientific basis for further development and adjustment of plague prevention and control strategies. IMPORTANCE Yunnan is the origin of the third plague pandemic. The analysis of human and animal plague characteristics of plague foci in Yunnan enlightens the prevention and control of the next plague pandemics. The plague characteristics of Yunnan show that human plague occurred when animal plague reached a certain scale, and strengthened surveillance of animal plague and reducing the density of host animals and transmission vectors contribute to the prevention and control of human plague outbreaks. The phenomenon of alternation between the resting period and active period of plague foci in Yunnan further proves the endogenous preservation mechanism of plague.

Evolutionary Stability Optimizer (ESO): A Novel Approach to Identify and Avoid Mutational Hotspots in DNA Sequences While Maintaining High Expression Levels

Modern synthetic biology procedures rely on the ability to generate stable genetic constructs that keep their functionality over long periods of time. However, maintenance of these constructs requires energy from the cell and thus reduces the host’s fitness. Natural selection results in loss-of-functionality mutations that negate the expression of the construct in the population. Current approaches for the prevention of this phenomenon focus on either small-scale, manual design of evolutionary stable constructs or the detection of mutational sites with unstable tendencies. We designed the Evolutionary Stability Optimizer (ESO), a software tool that enables the large-scale automatic design of evolutionarily stable constructs with respect to both mutational and epigenetic hotspots and allows users to define custom hotspots to avoid. Furthermore, our tool takes the expression of the input constructs into account by considering the guanine-cytosine (GC) content and codon usage of the host organism, balancing the trade-off between stability and gene expression, allowing to increase evolutionary stability while maintaining the high expression. In this study, we present the many features of the ESO and show that it accurately predicts the evolutionary stability of endogenous genes. The ESO was created as an easy-to-use, flexible platform based on the notion that directed genetic stability research will continue to evolve and revolutionize current applications of synthetic biology. The ESO is available at the following link: https://www.cs.tau.ac.il/~tamirtul/ESO/.

Burrowing owls as potential phoretic hosts of ground squirrel fleas during a plague epizootic

During the course of a plague epizootic, decimation of rodent host populations may result in the transfer of fleas to alternate or phoretic hosts, including to sympatric raptors that prey on rodents. We studied flea abundance and flea species assemblages on burrowing owls (Athene cunicularia hypugaea) in southwestern Idaho before (2012 - 2014), during (2015 - 2016), and after (2017) an epizootic of plague in Piute ground squirrels (Urocitellis mollis). We examined (1) if a larger proportion of burrowing owl nests contained fleas, (2) the likelihood that owls within a high flea abundance class increased, and (3) if owls harbored ground squirrel fleas during the epizootic. Using a flea abundance index assigned to 1,184 owls from 236 nests, the proportion of nests and the likelihood that owls had high flea abundance decreased rather than increased during epizootic years. Moreover, of 3,538 collected fleas from owls at 143 nests, no fleas were species that Piute ground squirrels typically harbor. Instead, Pulex irritans was the predominant flea collected in all study years (> 99%). Thus, although raptors may play a role in plague, there was no evidence that the die-off of ground squirrels resulted in higher flea intensity in burrowing owls or that they served as frequent accidental or phoretic hosts for ground squirrel fleas that could potentially be infectious with Yersinia pestis.

No evidence for enzootic plague within black-tailed prairie dog (Cynomys ludovicianus) populations

Yersinia pestis, causative agent of plague, occurs throughout the western United States in rodent populations and periodically causes epizootics in susceptible species, including black‐tailed prairie dogs (Cynomys ludovicianus). How Y. pestis persists long‐term in the environment between these epizootics is poorly understood but multiple mechanisms have been proposed, including, among others, a separate enzootic transmission cycle that maintains Y. pestis without involvement of epizootic hosts and persistence of Y. pestis within epizootic host populations without causing high mortality within those populations. We live‐trapped and collected fleas from black‐tailed prairie dogs and other mammal species from sites with and without black‐tailed prairie dogs in 2004 and 2005 and tested all fleas for presence of Y. pestis. Y. pestis was not detected in 2126 fleas collected in 2004 but was detected in 294 fleas collected from multiple sites in 2005, before and during a widespread epizootic that drastically reduced black‐tailed prairie dog populations in the affected colonies. Temporal and spatial patterns of Y. pestis occurrence in fleas and genotyping of Y. pestis present in some infected fleas suggest Y. pestis was introduced multiple times from sources outside the study area and once introduced, was dispersed between several sites. We conclude Y. pestis likely was not present in these black‐tailed prairie dog colonies prior to epizootic activity in these colonies. Although we did not identify likely enzootic hosts, we found evidence that deer mice (Peromyscus maniculatus) may serve as bridging hosts for Y. pestis between unknown enzootic hosts and black‐tailed prairie dogs.

Ecological Dynamics Impacting Bluetongue Virus Transmission in North America

Bluetongue virus (BTV) is an arbovirus transmitted to domestic and wild ruminants by certain species of Culicoides midges. The disease resulting from infection with BTV is economically important and can influence international trade and movement of livestock, the economics of livestock production, and animal welfare. Recent changes in the epidemiology of Culicoides-transmitted viruses, notably the emergence of exotic BTV genotypes in Europe, have demonstrated the devastating economic consequences of BTV epizootics and the complex nature of transmission across host-vector landscapes. Incursions of novel BTV serotypes into historically enzootic countries or regions, including the southeastern United States (US), Israel, Australia, and South America, have also occurred, suggesting diverse pathways for the transmission of these viruses. The abundance of BTV strains and multiple reassortant viruses circulating in Europe and the US in recent years demonstrates considerable genetic diversity of BTV strains and implies a history of reassortment events within the respective regions. While a great deal of emphasis is rightly placed on understanding the epidemiology and emergence of BTV beyond its natural ecosystem, the ecological contexts in which BTV maintains an enzootic cycle may also be of great significance. This review focuses on describing our current knowledge of ecological factors driving BTV transmission in North America. Information presented in this review can help inform future studies that may elucidate factors that are relevant to longstanding and emerging challenges associated with prevention of this disease.

Prairie Dogs, Persistent Plague, Flocking Fleas, and Pernicious Positive Feedback

Plague (caused by the bacterium Yersinia pestis) is a deadly flea-borne disease that remains a threat to public health nearly worldwide and is particularly disruptive ecologically where it has been introduced. We review hypotheses regarding maintenance and transmission of Y. pestis, emphasizing recent data from North America supporting maintenance by persistent transmission that results in sustained non-epizootic (but variable) rates of mortality in hosts. This maintenance mechanism may facilitate periodic epizootic eruptions "in place" because the need for repeated reinvasion from disjunct sources is eliminated. Resulting explosive outbreaks that spread rapidly in time and space are likely enhanced by synergistic positive feedback (PFB) cycles involving flea vectors, hosts, and the plague bacterium itself. Although PFB has been implied in plague literature for at least 50 years, we propose this mechanism, particularly with regard to flea responses, as central to epizootic plague rather than a phenomenon worthy of just peripheral mention. We also present new data on increases in flea:host ratios resulting from recreational shooting and poisoning as possible triggers for the transition from enzootic maintenance to PFB cycles and epizootic explosions. Although plague outbreaks have received much historic attention, PFB cycles that result in decimation of host populations lead to speculation that epizootic eruptions might not be part of the adaptive evolutionary strategy of Y. pestis but might instead be a tolerated intermittent cost of its modus operandi. We also speculate that there may be mammal communities where epizootics, as we define them, are rare or absent. Absence of plague epizootics might translate into reduced public health risk but does not necessarily equate to inconsequential ecologic impact.

Replacement, Refinement, and Reduction in Animal Studies With Biohazardous Agents

Animal models are critical to the advancement of our knowledge of infectious disease pathogenesis, diagnostics, therapeutics, and prevention strategies. The use of animal models requires thoughtful consideration for their well-being, as infections can significantly impact the general health of an animal and impair their welfare. Application of the 3Rs—replacement, refinement, and reduction—to animal models using biohazardous agents can improve the scientific merit and animal welfare. Replacement of animal models can use in vitro techniques such as cell culture systems, mathematical models, and engineered tissues or invertebrate animal hosts such as amoeba, worms, fruit flies, and cockroaches. Refinements can use a variety of techniques to more closely monitor the course of disease. These include the use of biomarkers, body temperature, behavioral observations, and clinical scoring systems. Reduction is possible using advanced technologies such as in vivo telemetry and imaging, allowing longitudinal assessment of animals during the course of disease. While there is no single method to universally replace, refine, or reduce animal models, the alternatives and techniques discussed are broadly applicable and they should be considered when infectious disease animal models are developed.

A single introduction of Yersinia pestis to Brazil during the 3rd plague pandemic

Yersinia pestis was introduced to Brazil during the third plague pandemic and currently exists in several recognized foci. There is currently limited available phylogeographic data regarding Y. pestis in Brazil. We generated whole genome sequences for 411 Y. pestis strains from six Brazilian foci to investigate the phylogeography of Y. pestis in Brazil; these strains were isolated from 1966 to 1997. All 411 strains were assigned to a single monophyletic clade within the 1.ORI population, indicating a single Y. pestis introduction was responsible for the successful establishment of endemic foci in Brazil. There was a moderate level of genomic diversity but little population structure among the 411 Brazilian Y. pestis strains, consistent with a radial expansion wherein Y. pestis spread rapidly from the coast to the interior of Brazil and became ecologically established. Overall, there were no strong spatial or temporal patterns among the Brazilian strains. However, strains from the same focus tended to be more closely related and strains isolated from foci closer to the coast tended to fall in more basal positions in the whole genome phylogeny than strains from more interior foci. Overall, the patterns observed in Brazil are similar to other locations affected during the 3rd plague pandemic such as in North America and Madagascar.

Genome-based population structure analysis of the strawberry plant pathogen Xanthomonas fragariae reveals two distinct groups that evolved independently before its species description

Xanthomonas fragariae is a quarantine organism in Europe, causing angular leaf spots on strawberry plants. It is spreading worldwide in strawberry-producing regions due to import of plant material through trade and human activities. In order to resolve the population structure at the strain level, we have employed high-resolution molecular typing tools on a comprehensive strain collection representing global and temporal distribution of the pathogen. Clustered regularly interspaced short palindromic repeat regions (CRISPRs) and variable number of tandem repeats (VNTRs) were identified within the reference genome of X. fragariae LMG 25863 as a potential source of variation. Strains from our collection were whole-genome sequenced and used in order to identify variable spacers and repeats for discriminative purpose. CRISPR spacer analysis and multiple-locus VNTR analysis (MLVA) displayed a congruent population structure, in which two major groups and a total of four subgroups were revealed. The two main groups were genetically separated before the first X. fragariae isolate was described and are potentially responsible for the worldwide expansion of the bacterial disease. Three primer sets were designed for discriminating CRISPR-associated markers in order to streamline group determination of novel isolates. Overall, this study describes typing methods to discriminate strains and monitor the pathogen population structure, more especially in the view of a new outbreak of the pathogen.

Yersinia pestis Survival and Replication in Potential Ameba Reservoir

Plague ecology is characterized by sporadic epizootics, then periods of dormancy. Building evidence suggests environmentally ubiquitous amebae act as feral macrophages and hosts to many intracellular pathogens. We conducted environmental genetic surveys and laboratory co-culture infection experiments to assess whether plague bacteria were resistant to digestion by 5 environmental ameba species. First, we demonstrated that Yersinia pestis is resistant or transiently resistant to various ameba species. Second, we showed that Y. pestis survives and replicates intracellularly within Dictyostelium discoideum amebae for ˃48 hours postinfection, whereas control bacteria were destroyed in <1 hour. Finally, we found that Y. pestis resides within ameba structures synonymous with those found in infected human macrophages, for which Y. pestis is a competent pathogen. Evidence supporting amebae as potential plague reservoirs stresses the importance of recognizing pathogen-harboring amebae as threats to public health, agriculture, conservation, and biodefense.

Temporal phylogeography of Yersinia pestis in Madagascar: Insights into the long-term maintenance of plague

BACKGROUND

Yersinia pestis appears to be maintained in multiple, geographically separate, and phylogenetically distinct subpopulations within the highlands of Madagascar. However, the dynamics of these locally differentiated subpopulations through time are mostly unknown. To address that gap and further inform our understanding of plague epidemiology, we investigated the phylogeography of Y. pestis in Madagascar over an 18 year period.

METHODOLOGY/PRINCIPAL FINDINGS

We generated whole genome sequences for 31 strains and discovered new SNPs that we used in conjunction with previously identified SNPs and variable-number tandem repeats (VNTRs) to genotype 773 Malagasy Y. pestis samples from 1995 to 2012. We mapped the locations where samples were obtained on a fine geographic scale to examine phylogeographic patterns through time. We identified 18 geographically separate and phylogenetically distinct subpopulations that display spatial and temporal stability, persisting in the same locations over a period of almost two decades. We found that geographic areas with higher levels of topographical relief are associated with greater levels of phylogenetic diversity and that sampling frequency can vary considerably among subpopulations and from year to year. We also found evidence of various Y. pestis dispersal events, including over long distances, but no evidence that any dispersal events resulted in successful establishment of a transferred genotype in a new location during the examined time period.

CONCLUSIONS/SIGNIFICANCE

Our analysis suggests that persistent endemic cycles of Y. pestis transmission within local areas are responsible for the long term maintenance of plague in Madagascar, rather than repeated episodes of wide scale epidemic spread. Landscape likely plays a role in maintaining Y. pestis subpopulations in Madagascar, with increased topographical relief associated with increased levels of localized differentiation. Local ecological factors likely affect the dynamics of individual subpopulations and the associated likelihood of observing human plague cases in a given year in a particular location.

Protein abundances can distinguish between naturally-occurring and laboratory strains of Yersinia pestis, the causative agent of plague

The rapid pace of bacterial evolution enables organisms to adapt to the laboratory environment with repeated passage and thus diverge from naturally-occurring environmental ("wild") strains. Distinguishing wild and laboratory strains is clearly important for biodefense and bioforensics; however, DNA sequence data alone has thus far not provided a clear signature, perhaps due to lack of understanding of how diverse genome changes lead to convergent phenotypes, difficulty in detecting certain types of mutations, or perhaps because some adaptive modifications are epigenetic. Monitoring protein abundance, a molecular measure of phenotype, can overcome some of these difficulties. We have assembled a collection of Yersinia pestis proteomics datasets from our own published and unpublished work, and from a proteomics data archive, and demonstrated that protein abundance data can clearly distinguish laboratory-adapted from wild. We developed a lasso logistic regression classifier that uses binary (presence/absence) or quantitative protein abundance measures to predict whether a sample is laboratory-adapted or wild that proved to be ~98% accurate, as judged by replicated 10-fold cross-validation. Protein features selected by the classifier accord well with our previous study of laboratory adaptation in Y. pestis. The input data was derived from a variety of unrelated experiments and contained significant confounding variables. We show that the classifier is robust with respect to these variables. The methodology is able to discover signatures for laboratory facility and culture medium that are largely independent of the signature of laboratory adaptation. Going beyond our previous laboratory evolution study, this work suggests that proteomic differences between laboratory-adapted and wild Y. pestis are general, potentially pointing to a process that could apply to other species as well. Additionally, we show that proteomics datasets (even archived data collected for different purposes) contain the information necessary to distinguish wild and laboratory samples. This work has clear applications in biomarker detection as well as biodefense.

Multiple Locus Variable number of tandem repeat Analysis: A molecular genotyping tool for Paenibacillus larvae

American Foulbrood, caused by Paenibacillus larvae, is the most severe bacterial disease of honey bees (Apis mellifera). To perform genotyping of P. larvae in an epidemiological context, there is a need of a fast and cheap method with a high resolution. Here, we propose Multiple Locus Variable number of tandem repeat Analysis (MLVA). MLVA has been used for typing a collection of 209 P. larvae strains from which 23 different MLVA types could be identified. Moreover, the developed methodology not only permits the identification of the four Enterobacterial Repetitive Intergenic Consensus (ERIC) genotypes, but allows also a discriminatory subdivision of the most dominant ERIC type I and ERIC type II genotypes. A biogeographical study has been conducted showing a significant correlation between MLVA genotype and the geographical region where it was isolated.

VNTR diversity in Yersinia pestis isolates from an animal challenge study reveals the potential for in vitro mutations during laboratory cultivation

Underlying mutation rates and other evolutionary forces shape the population structure of bacteria in nature. Although easily overlooked, similar forces are at work in the laboratory and may influence observed mutations. Here, we investigated tissue samples and Yersinia pestis isolates from a rodent laboratory challenge with strain CO92 using whole genome sequencing and multi-locus variable-number tandem repeat (VNTR) analysis (MLVA). We identified six VNTR mutations that were found to have occurred in vitro during laboratory cultivation rather than in vivo during the rodent challenge. In contrast, no single nucleotide polymorphism (SNP) mutations were observed, either in vivo or in vitro. These results were consistent with previously published mutation rates and the calculated number of Y. pestis generations that occurred during the in vitro versus the in vivo portions of the experiment. When genotyping disease outbreaks, the potential for in vitro mutations should be considered, particularly when highly variable genetic markers such as VNTRs are used.

Microbial Genomics of a Host-Associated Commensal Bacterium in Fragmented Populations of Endangered Takahe

Isolation of wildlife into fragmented populations as a consequence of anthropogenic-mediated environmental change may alter host-pathogen relationships. Our understanding of some of the epidemiological features of infectious disease in vulnerable populations can be enhanced by the use of commensal bacteria as a proxy for invasive pathogens in natural ecosystems. The distinctive population structure of a well-described meta-population of a New Zealand endangered flightless bird, the takahe (Porphyrio hochstetteri), provided a unique opportunity to investigate the influence of host isolation on enteric microbial diversity. The genomic epidemiology of a prevalent rail-associated endemic commensal bacterium was explored using core genome and ribosomal multilocus sequence typing (rMLST) of 70 Campylobacter sp. nova 1 isolated from one third of the takahe population resident in multiple locations. While there was evidence of recombination between lineages, bacterial divergence appears to have occurred and multivariate analysis of 52 rMLST genes revealed location-associated differentiation of C. sp. nova 1 sequence types. Our results indicate that fragmentation and anthropogenic manipulation of populations can influence host-microbial relationships, with potential implications for niche adaptation and the evolution of micro-organisms in remote environments. This study provides a novel framework in which to explore the complex genomic epidemiology of micro-organisms in wildlife populations.

Genetic variation at the MHC DRB1 locus is similar across Gunnison's prairie dog (Cynomys gunnisoni) colonies regardless of plague history

Yersinia pestis was introduced to North America around 1900 and leads to nearly 100% mortality in prairie dog (Cynomys spp.) colonies during epizootic events, which suggests this pathogen may exert a strong selective force. We characterized genetic diversity at an MHC class II locus (DRB1) in Gunnison's prairie dog (C. gunnisoni) and quantified population genetic structure at the DRB1 versus 12 microsatellite loci in three large Arizona colonies. Two colonies, Seligman (SE) and Espee Ranch (ES), have experienced multiple plague‐related die‐offs in recent years, whereas plague has never been documented at Aubrey Valley (AV). We found fairly low allelic diversity at the DRB1 locus, with one allele (DRB1*01) at high frequency (0.67–0.87) in all colonies. Two other DRB1 alleles appear to be trans‐species polymorphisms shared with the black‐tailed prairie dog (C. ludovicianus), indicating that these alleles have been maintained across evolutionary time frames. Estimates of genetic differentiation were generally lower at the MHC locus (F_ST = 0.033) than at microsatellite markers (F_ST = 0.098). The reduced differentiation at DRB1 may indicate that selection has been important for shaping variation at MHC loci, regardless of the presence or absence of plague in recent decades. However, genetic drift has probably also influenced the DRB1 locus because its level of differentiation was not different from that of microsatellites in an F_ST outlier analysis. We then compared specific MHC alleles to plague survivorship in 60 C. gunnisoni that had been experimentally infected with Y. pestis. We found that survival was greater in individuals that carried at least one copy of the most common allele (DRB1*01) compared to those that did not (60% vs. 20%). Although the sample sizes of these two groups were unbalanced, this result suggests the possibility that this MHC class II locus, or a nearby linked gene, could play a role in plague survival.

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.

Ecological Traits Driving the Outbreaks and Emergence of Zoonotic Pathogens

Infectious diseases that are transmitted from wildlife hosts to humans, such as the Ebola virus and MERS virus, can be difficult to understand because the pathogens emerge from complex multifaceted ecological interactions. We use a wildlife–pathogen system—prairie dogs (Cynomys ludovicianus) and the plague bacterium (Yersinia pestis)—to describe aspects of disease ecology that apply to many cases of emerging infectious disease. We show that the monitoring and surveillance of hosts and vectors during the buildup to disease outbreaks are crucial for understanding pathogen-transmission dynamics and that a community-ecology framework is important to identify reservoir hosts. Incorporating multidisciplinary approaches and frameworks may improve wildlife–pathogen surveillance and our understanding of seemingly sporadic and rare pathogen outbreaks.

Genome and Evolution of Yersinia pestis

This chapter summarizes researches on genome and evolution features of Yersinia pestis, the young pathogen that evolved from Y. pseudotuberculosis at least 5000 years ago. Y. pestis is a highly clonal bacterial species with closed pan-genome. Comparative genomic analysis revealed that genome of Y. pestis experienced highly frequent rearrangement and genome decay events during the evolution. The genealogy of Y. pestis includes five major branches, and four of them seemed raised from a "big bang" node that is associated with the Black Death. Although whole genome-wide variation of Y. pestis reflected a neutral evolutionary process, the branch length in the genealogical tree revealed over dispersion, which was supposedly caused by varied historical molecular clock that is associated with demographical effect by alternate cycles of enzootic disease and epizootic disease in sylvatic plague foci. In recent years, palaeomicrobiology researches on victims of the Black Death, and Justinian's plague verified that two historical pandemics were indeed caused by Y. pestis, but the etiological lineages might be extinct today.

Investigation of Yersinia pestis Laboratory Adaptation through a Combined Genomics and Proteomics Approach

The bacterial pathogen Yersinia pestis, the cause of plague in humans and animals, normally has a sylvatic lifestyle, cycling between fleas and mammals. In contrast, laboratory-grown Y. pestis experiences a more constant environment and conditions that it would not normally encounter. The transition from the natural environment to the laboratory results in a vastly different set of selective pressures, and represents what could be considered domestication. Understanding the kinds of adaptations Y. pestis undergoes as it becomes domesticated will contribute to understanding the basic biology of this important pathogen.

In this study, we performed a parallel serial passage experiment (PSPE) to explore the mechanisms by which Y. pestis adapts to laboratory conditions, hypothesizing that cells would undergo significant changes in virulence and nutrient acquisition systems. Two wild strains were serially passaged in 12 independent populations each for ~750 generations, after which each population was analyzed using whole-genome sequencing, LC-MS/MS proteomic analysis, and GC/MS metabolomics. We observed considerable parallel evolution in the endpoint populations, detecting multiple independent mutations in ail, pepA, and zwf, suggesting that specific selective pressures are shaping evolutionary responses. Complementary LC-MS/MS proteomic data provide physiological context to the observed mutations, and reveal regulatory changes not necessarily associated with specific mutations, including changes in amino acid metabolism and cell envelope biogenesis. Proteomic data support hypotheses generated by genomic data in addition to suggesting future mechanistic studies, indicating that future whole-genome sequencing studies be designed to leverage proteomics as a critical complement.

Population genetic structure of a common host predicts the spread of white-nose syndrome, an emerging infectious disease in bats

Landscape complexity influences patterns of animal dispersal, which in turn may affect both gene flow and the spread of pathogens. White-nose syndrome (WNS) is an introduced fungal disease that has spread rapidly throughout eastern North America, causing massive mortality in bat populations. We tested for a relationship between the population genetic structure of the most common host, the little brown myotis (Myotis lucifugus), and the geographic spread of WNS to date by evaluating logistic regression models of WNS risk among hibernating colonies in eastern North America. We hypothesized that risk of WNS to susceptible host colonies should increase with both geographic proximity and genetic similarity, reflecting historical connectivity, to infected colonies. Consistent with this hypothesis, inclusion of genetic distance between infected and susceptible colonies significantly improved models of disease spread, capturing heterogeneity in the spatial expansion of WNS despite low levels of genetic differentiation among eastern populations. Expanding our genetic analysis to the continental range of little brown myotis reveals strongly contrasting patterns of population structure between eastern and western North America. Genetic structure increases markedly moving westward into the northern Great Plains, beyond the current distribution of WNS. In western North America, genetic differentiation of geographically proximate populations often exceeds levels observed across the entire eastern region, suggesting infrequent and/or locally restricted dispersal, and thus relatively limited opportunities for pathogen introduction in western North America. Taken together, our analyses suggest a possibly slower future rate of spread of the WNS pathogen, at least as mediated by little brown myotis.

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.

Predicting the Genetic Stability of Engineered DNA Sequences with the EFM Calculator

Unwanted evolution can rapidly degrade the performance of genetically engineered circuits and metabolic pathways installed in living organisms. We created the Evolutionary Failure Mode (EFM) Calculator to computationally detect common sources of genetic instability in an input DNA sequence. It predicts two types of mutational hotspots: deletions mediated by homologous recombination and indels caused by replication slippage on simple sequence repeats. We tested the performance of our algorithm on genetic circuits that were previously redesigned for greater evolutionary reliability and analyzed the stability of sequences in the iGEM Registry of Standard Biological Parts. More than half of the parts in the Registry are predicted to experience >100-fold elevated mutation rates due to the inclusion of unstable sequence configurations. We anticipate that the EFM Calculator will be a useful negative design tool for avoiding volatile DNA encodings, thereby increasing the evolutionary lifetimes of synthetic biology devices.

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.

Age at Vaccination May Influence Response to Sylvatic Plague Vaccine (SPV) in Gunnison's Prairie Dogs (Cynomys gunnisoni)

Gunnison's prairie dogs (Cynomys gunnisoni) have been considered at greater risk from Yersinia pestis (plague) infection in the montane portion of their range compared to populations at lower elevations, possibly due to factors related to flea transmission of the bacteria or greater host susceptibility. To test the latter hypothesis and determine whether vaccination against plague with an oral sylvatic plague vaccine (SPV) improved survival, we captured prairie dogs from a C. g. gunnisoni or "montane" population and a C. g. zuniensis or "prairie" population for vaccine efficacy and challenge studies. No differences (P = 0.63) were found in plague susceptibility in non-vaccinated animals between these two populations; however, vaccinates from the prairie population survived plague challenge at significantly higher rates (P < 0.01) than those from the montane population. Upon further analysis, we determined that response to immunization was most likely associated with differences in age, as the prairie group was much younger on average than the montane group. Vaccinates that were juveniles or young adults survived plague challenge at a much higher rate than adults (P < 0.01 and P = 0.02, respectively), but no difference (P = 0.83) was detected in survival rates between control animals of different ages. These results suggest that host susceptibility is probably not related to the assumed greater risk from plague in the C. g. gunnisoni or "montane" populations of Gunnison's prairie dogs, and that SPV could be a useful plague management tool for this species, particularly if targeted at younger cohorts.

Single-Nucleotide Polymorphisms Reveal Spatial Diversity Among Clones of Yersinia pestis During Plague Outbreaks in Colorado and the Western United States

BACKGROUND

In western North America, plague epizootics caused by Yersinia pestis appear to sweep across landscapes, primarily infecting and killing rodents, especially ground squirrels and prairie dogs. During these epizootics, the risk of Y. pestis transmission to humans is highest. While empirical models that include climatic conditions and densities of rodent hosts and fleas can predict when epizootics are triggered, bacterial transmission patterns across landscapes, and the scale at which Y. pestis is maintained in nature during inter-epizootic periods, are poorly defined. Elucidating the spatial extent of Y. pestis clones during epizootics can determine whether bacteria are propagated across landscapes or arise independently from local inter-epizootic maintenance reservoirs.

MATERIAL AND METHODS

We used DNA microarray technology to identify single-nucleotide polymorphisms (SNPs) in 34 Y. pestis isolates collected in the western United States from 1980 to 2006, 21 of which were collected during plague epizootics in Colorado. Phylogenetic comparisons were used to elucidate the hypothesized spread of Y. pestis between the mountainous Front Range and the eastern plains of northern Colorado during epizootics. Isolates collected from across the western United States were included for regional comparisons.

RESULTS

By identifying SNPs that mark individual clones, our results strongly suggest that Y. pestis is maintained locally and that widespread epizootic activity is caused by multiple clones arising independently at small geographic scales. This is in contrast to propagation of individual clones being transported widely across landscapes. Regionally, our data are consistent with the notion that Y. pestis diversifies at relatively local scales following long-range translocation events. We recommend that surveillance and prediction by public health and wildlife management professionals focus more on models of local or regional weather patterns and ecological factors that may increase risk of widespread epizootics, rather than predicting or attempting to explain epizootics on the basis of movement of host species that may transport plague.

Analysis of Salmonella enterica serovar Typhimurium variable-number tandem-repeat data for public health investigation based on measured mutation rates and whole-genome sequence comparisons

Variable-number tandem repeats (VNTRs) mutate rapidly and can be useful markers for genotyping. While multilocus VNTR analysis (MLVA) is increasingly used in the detection and investigation of food-borne outbreaks caused by Salmonella enterica serovar Typhimurium (S. Typhimurium) and other bacterial pathogens, MLVA data analysis usually relies on simple clustering approaches that may lead to incorrect interpretations. Here, we estimated the rates of copy number change at each of the five loci commonly used for S. Typhimurium MLVA, during in vitro and in vivo passage. We found that loci STTR5, STTR6, and STTR10 changed during passage but STTR3 and STTR9 did not. Relative rates of change were consistent across in vitro and in vivo growth and could be accurately estimated from diversity measures of natural variation observed during large outbreaks. Using a set of 203 isolates from a series of linked outbreaks and whole-genome sequencing of 12 representative isolates, we assessed the accuracy and utility of several alternative methods for analyzing and interpreting S. Typhimurium MLVA data. We show that eBURST analysis was accurate and informative. For construction of MLVA-based trees, a novel distance metric, based on the geometric model of VNTR evolution coupled with locus-specific weights, performed better than the commonly used simple or categorical distance metrics. The data suggest that, for the purpose of identifying potential transmission clusters for further investigation, isolates whose profiles differ at one of the rapidly changing STTR5, STTR6, and STTR10 loci should be collapsed into the same cluster.

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.

Phylogeography and population structure of the biologically invasive phytopathogen Erwinia amylovora inferred using minisatellites

Erwinia amylovora causes a major disease of pome fruit trees worldwide, and is regulated as a quarantine organism in many countries. While some diversity of isolates has been observed, molecular epidemiology of this bacterium is hindered by a lack of simple molecular typing techniques with sufficiently high resolution. We report a molecular typing system of E. amylovora based on variable number of tandem repeats (VNTR) analysis. Repeats in the E. amylovora genome were identified with comparative genomic tools, and VNTR markers were developed and validated. A Multiple-Locus VNTR Analysis (MLVA) was applied to E. amylovora isolates from bacterial collections representing global and regional distribution of the pathogen. Based on six repeats, MLVA allowed the distinction of 227 haplotypes among a collection of 833 isolates of worldwide origin. Three geographically separated groups were recognized among global isolates using Bayesian clustering methods. Analysis of regional outbreaks confirmed presence of diverse haplotypes but also high representation of certain haplotypes during outbreaks. MLVA analysis is a practical method for epidemiological studies of E. amylovora, identifying previously unresolved population structure within outbreaks. Knowledge of such structure can increase our understanding on how plant diseases emerge and spread over a given geographical region.

Dynamic of mutational events in variable number tandem repeats of Escherichia coli O157:H7

VNTRs regions have been successfully used for bacterial subtyping; however, the hypervariability in VNTR loci is problematic when trying to predict the relationships among isolates. Since few studies have examined the mutation rate of these markers, our aim was to estimate mutation rates of VNTRs specific for verotoxigenic E. coli O157:H7. The knowledge of VNTR mutational rates and the factors affecting them would make MLVA more effective for epidemiological or microbial forensic investigations. For this purpose, we analyzed nine loci performing parallel, serial passage experiments (PSPEs) on 9 O157:H7 strains. The combined 9 PSPE population rates for the 8 mutating loci ranged from 4.4 × 10⁻⁰⁵ to 1.8 × 10⁻⁰³ mutations/generation, and the combined 8-loci mutation rate was of 2.5 × 10⁻⁰³ mutations/generation. Mutations involved complete repeat units, with only one point mutation detected. A similar proportion between single and multiple repeat changes was detected. Of the 56 repeat mutations, 59% were insertions and 41% were deletions, and 72% of the mutation events corresponded to O157-10 locus. For alleles with up to 13 UR, a constant and low mutation rate was observed; meanwhile longer alleles were associated with higher and variable mutation rates. Our results are useful to interpret data from microevolution and population epidemiology studies and particularly point out that the inclusion or not of O157-10 locus or, alternatively, a differential weighting data according to the mutation rates of loci must be evaluated in relation with the objectives of the proposed study.

Duration of plague (Yersinia pestis) outbreaks in black-tailed prairie dog (Cynomys ludovicianus) colonies of northern Colorado

Plague, caused by the bacterium Yersinia pestis, triggers die-offs in colonies of black-tailed prairie dogs (Cynomys ludovicianus), but the time-frame of plague activity is not well understood. We document plague activity in fleas from prairie dogs and their burrows on three prairie dog colonies that suffered die-offs. We demonstrate that Y. pestis transmission occurs over periods from several months to over a year in prairie dog populations before observed die-offs.

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.

Plague circulation and population genetics of the reservoir Rattus rattus: the influence of topographic relief on the distribution of the disease within the Madagascan focus

BACKGROUND

Landscape may affect the distribution of infectious diseases by influencing the population density and dispersal of hosts and vectors. Plague (Yersinia pestis infection) is a highly virulent, re-emerging disease, the ecology of which has been scarcely studied in Africa. Human seroprevalence data for the major plague focus of Madagascar suggest that plague spreads heterogeneously across the landscape as a function of the relief. Plague is primarily a disease of rodents. We therefore investigated the relationship between disease distribution and the population genetic structure of the black rat, Rattus rattus, the main reservoir of plague in Madagascar.

METHODOLOGY/PRINCIPAL FINDINGS

We conducted a comparative study of plague seroprevalence and genetic structure (15 microsatellite markers) in rat populations from four geographic areas differing in topology, each covering about 150-200 km(2) within the Madagascan plague focus. The seroprevalence levels in the rat populations mimicked those previously reported for humans. As expected, rat populations clearly displayed a more marked genetic structure with increasing relief. However, the relationship between seroprevalence data and genetic structure differs between areas, suggesting that plague distribution is not related everywhere to the effective dispersal of rats.

CONCLUSIONS/SIGNIFICANCE

Genetic diversity estimates suggested that plague epizootics had only a weak impact on rat population sizes. In the highlands of Madagascar, plague dissemination cannot be accounted for solely by the effective dispersal of the reservoir. Human social activities may also be involved in spreading the disease in rat and human populations.

Host resistance, population structure and the long-term persistence of bubonic plague: contributions of a modelling approach in the Malagasy focus

Although bubonic plague is an endemic zoonosis in many countries around the world, the factors responsible for the persistence of this highly virulent disease remain poorly known. Classically, the endemic persistence of plague is suspected to be due to the coexistence of plague resistant and plague susceptible rodents in natural foci, and/or to a metapopulation structure of reservoirs. Here, we test separately the effect of each of these factors on the long-term persistence of plague. We analyse the dynamics and equilibria of a model of plague propagation, consistent with plague ecology in Madagascar, a major focus where this disease is endemic since the 1920s in central highlands. By combining deterministic and stochastic analyses of this model, and including sensitivity analyses, we show that (i) endemicity is favoured by intermediate host population sizes, (ii) in large host populations, the presence of resistant rats is sufficient to explain long-term persistence of plague, and (iii) the metapopulation structure of susceptible host populations alone can also account for plague endemicity, thanks to both subdivision and the subsequent reduction in the size of subpopulations, and extinction-recolonization dynamics of the disease. In the light of these results, we suggest scenarios to explain the localized presence of plague in Madagascar.

Bubonic plague, known to have marked human history by three deadly pandemics, is an infectious disease which mainly circulates in wild rodent populations and is transmitted by fleas. Although this disease can be quickly lethal to its host, it has persisted on long-term in many rodent populations around the world. The reasons for this persistence remain poorly known. Two mechanisms have been invoked, but not yet explicitly and independently tested: first, the spatial structure of rodent populations (subdivision into several subpopulations) and secondly, the presence of, not only plague-susceptible rodents, but also plague-resistant ones. To gain insight into the role of the above two factors in plague persistence, we analysed a mathematical model of plague propagation. We applied our analyses to the case of Madagascar, where plague has persisted on central highlands since the 1920s and is responsible for about 30% of the human cases worldwide. We found that the long-term persistence of plague can be explained by the presence of any of the above two factors. These results allowed us to propose scenarios to explain the localized presence of plague in the Malagasy highlands, and help understand the persistence of plague in many wild foci.

Within-host evolution of Brucella canis during a canine brucellosis outbreak in a kennel

Background

Little is currently known about Brucella evolution within the host during infection. The current study is the first to employ fine-scale genotyping on an isolate collection derived from a Brucella canis outbreak. Eight isolates of B. canis, cultured from different tissues of three dogs (female, stud dog, puppy of another female) from a single kennel over three months were genetically characterized with a 15-marker multi-locus, variable-number tandem repeat (VNTR) analysis (MLVA) to assess the genetic relatedness of isolates and potential rapid mutational changes.

Results

MLVA discriminated among the otherwise indistinguishable isolates from different animals and from isolates collected at different time points within each host, with different VNTR alleles being detected at multiple dates and tissue sites. We suspect that all isolates cultured from the female, puppy, and stud dogs originated from the same strain, with subsequent rapid in vivo mutations. However, high mutation rates and apparent in several of the loci prevented making definitive epidemiological relationships among isolates.

Conclusions

This investigation highlights the rapid in vivo genetic mutations of several VNTRs of B. canis over a short time period in the host and the emergence of alternate alleles. However, this work also suggests the challenges of using highly mutable VNTRs to infer epidemiological relationships of strains within a short duration outbreak.

Rapid focused sequencing: a multiplexed assay for simultaneous detection and strain typing of Bacillus anthracis, Francisella tularensis, and Yersinia pestis

Background

The intentional release of Bacillus anthracis in the United States in 2001 has heightened concern about the use of pathogenic microorganisms in bioterrorism attacks. Many of the deadliest bacteria, including the Class A Select Agents Bacillus anthracis, Francisella tularensis, and Yersinia pestis, are highly infectious via the pulmonary route when released in aerosolized form. Hence, rapid, sensitive, and reliable methods for detection of these biothreats and characterization of their potential impact on the exposed population are of critical importance to initiate and support rapid military, public health, and clinical responses.

Methodology/Principal Findings

We have developed microfluidic multiplexed PCR and sequencing assays based on the simultaneous interrogation of three pathogens per assay and ten loci per pathogen. Microfluidic separation of amplified fluorescently labeled fragments generated characteristic electrophoretic signatures for identification of each agent. The three sets of primers allowed significant strain typing and discrimination from non-pathogenic closely-related species and environmental background strains based on amplicon sizes alone. Furthermore, sequencing of the 10 amplicons per pathogen, termed “Rapid Focused Sequencing,” allowed an even greater degree of strain discrimination and, in some cases, can be used to determine virulence. Both amplification and sequencing assays were performed in microfluidic biochips developed for fast thermal cycling and requiring 7 µL per reaction. The 30-plex sequencing assay resulted in genotypic resolution of 84 representative strains belonging to each of the three biothreat species.

Conclusions/Significance

The microfluidic multiplexed assays allowed identification and strain differentiation of the biothreat agents Bacillus anthracis, Francisella tularensis, and Yersinia pestis and clear discrimination from closely-related species and several environmental background strains. The assays may be extended to detect a large number of pathogens, are applicable to the evaluation of both environmental and clinical samples, and have the potential to be applied in military, public health, and clinical diagnostic settings.

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.

Historical variations in mutation rate in an epidemic pathogen, Yersinia pestis

The genetic diversity of Yersinia pestis, the etiologic agent of plague, is extremely limited because of its recent origin coupled with a slow clock rate. Here we identified 2,326 SNPs from 133 genomes of Y. pestis strains that were isolated in China and elsewhere. These SNPs define the genealogy of Y. pestis since its most recent common ancestor. All but 28 of these SNPs represented mutations that happened only once within the genealogy, and they were distributed essentially at random among individual genes. Only seven genes contained a significant excess of nonsynonymous SNP, suggesting that the fixation of SNPs mainly arises via neutral processes, such as genetic drift, rather than Darwinian selection. However, the rate of fixation varies dramatically over the genealogy: the number of SNPs accumulated by different lineages was highly variable and the genealogy contains multiple polytomies, one of which resulted in four branches near the time of the Black Death. We suggest that demographic changes can affect the speed of evolution in epidemic pathogens even in the absence of natural selection, and hypothesize that neutral SNPs are fixed rapidly during intermittent epidemics and outbreaks.

Yersinia pestis lineages in Mongolia

Background

Whole genome sequencing allowed the development of a number of high resolution sequence based typing tools for Yersinia (Y.) pestis. The application of these methods on isolates from most known foci worldwide and in particular from China and the Former Soviet Union has dramatically improved our understanding of the population structure of this species. In the current view, Y. pestis including the non or moderate human pathogen Y. pestis subspecies microtus emerged from Yersinia pseudotuberculosis about 2,600 to 28,600 years ago in central Asia. The majority of central Asia natural foci have been investigated. However these investigations included only few strains from Mongolia.

Methodology/Principal Findings

Clustered Regularly Interspaced Short Prokaryotic Repeats (CRISPR) analysis and Multiple-locus variable number of tandem repeats (VNTR) analysis (MLVA) with 25 loci was performed on 100 Y. pestis strains, isolated from 37 sampling areas in Mongolia. The resulting data were compared with previously published data from more than 500 plague strains, 130 of which had also been previously genotyped by single nucleotide polymorphism (SNP) analysis. The comparison revealed six main clusters including the three microtus biovars Ulegeica, Altaica, and Xilingolensis. The largest cluster comprises 78 isolates, with unique and new genotypes seen so far in Mongolia only. Typing of selected isolates by key SNPs was used to robustly assign the corresponding clusters to previously defined SNP branches.

Conclusions/Significance

We show that Mongolia hosts the most recent microtus clade (Ulegeica). Interestingly no representatives of the ancestral Y. pestis subspecies pestis nodes previously identified in North-western China were identified in this study. This observation suggests that the subsequent evolution steps within Y. pestis pestis did not occur in Mongolia. Rather, Mongolia was most likely re-colonized by more recent clades coming back from China contemporary of the black death pandemic, or more recently in the past 600 years.

Comparative genomics of 2009 seasonal plague (Yersinia pestis) in New Mexico

Plague disease caused by the Gram-negative bacterium Yersinia pestis routinely affects animals and occasionally humans, in the western United States. The strains native to the North American continent are thought to be derived from a single introduction in the late 19^th century. The degree to which these isolates have diverged genetically since their introduction is not clear, and new genomic markers to assay the diversity of North American plague are highly desired. To assay genetic diversity of plague isolates within confined geographic areas, draft genome sequences were generated by 454 pyrosequencing from nine environmental and clinical plague isolates. In silico assemblies of Variable Number Tandem Repeat (VNTR) loci were compared to laboratory-generated profiles for seven markers. High-confidence SNPs and small Insertion/Deletions (Indels) were compared to previously sequenced Y. pestis isolates. The resulting panel of mutations allowed clustering of the strains and tracing of the most likely evolutionary trajectory of the plague strains. The sequences also allowed the identification of new putative SNPs that differentiate the 2009 isolates from previously sequenced plague strains and from each other. In addition, new insertion points for the abundant insertion sequences (IS) of Y. pestis are present that allow additional discrimination of strains; several of these new insertions potentially inactivate genes implicated in virulence. These sequences enable whole-genome phylogenetic analysis and allow the unbiased comparison of closely related isolates of a genetically monomorphic pathogen.

Microevolution of pandemic Vibrio parahaemolyticus assessed by the number of repeat units in short sequence tandem repeat regions

The emergence of the pandemic strain Vibrio parahaemolyticus O3:K6 in 1996 caused a large increase of diarrhea outbreaks related to seafood consumption in Southeast Asia, and later worldwide. Isolates of this strain constitutes a clonal complex, and their effectual differentiation is possible by comparison of their variable number tandem repeats (VNTRs). The differentiation of the isolates by the differences in VNTRs will allow inferring the population dynamics and microevolution of this strain but this requires knowing the rate and mechanism of VNTRs' variation. Our study of mutants obtained after serial cultivation of clones showed that mutation rates of the six VNTRs examined are on the order of 10(-4) mutant per generation and that difference increases by stepwise addition of single mutations. The single stepwise mutation (SSM) was deduced because mutants with 1, 2, 3, or more repeat unit deletions or insertions follow a geometric distribution. Plausible phylogenetic trees are obtained when, according to SSM, the genetic distance between clusters with different number of repeats is assessed by the absolute differences in repeats. Using this approach, mutants originated from different isolates of pandemic V. parahaemolyticus after serial cultivation are clustered with their parental isolates. Additionally, isolates of pandemic V. parahaemolyticus from Southeast Asia, Tokyo, and northern and southern Chile are clustered according their geographical origin. The deepest split in these four populations is observed between the Tokyo and southern Chile populations. We conclude that proper phylogenetic relations and successful tracing of pandemic V. parahaemolyticus requires measuring the differences between isolates by the absolute number of repeats in the VNTRs considered.

Phylogeography and molecular epidemiology of Yersinia pestis in Madagascar

Background

Plague was introduced to Madagascar in 1898 and continues to be a significant human health problem. It exists mainly in the central highlands, but in the 1990s was reintroduced to the port city of Mahajanga, where it caused extensive human outbreaks. Despite its prevalence, the phylogeography and molecular epidemiology of Y. pestis in Madagascar has been difficult to study due to the great genetic similarity among isolates. We examine island-wide geographic-genetic patterns based upon whole-genome discovery of SNPs, SNP genotyping and hypervariable variable-number tandem repeat (VNTR) loci to gain insight into the maintenance and spread of Y. pestis in Madagascar.

Methodology/Principal Findings

We analyzed a set of 262 Malagasy isolates using a set of 56 SNPs and a 43-locus multi-locus VNTR analysis (MLVA) system. We then analyzed the geographic distribution of the subclades and identified patterns related to the maintenance and spread of plague in Madagascar. We find relatively high levels of VNTR diversity in addition to several SNP differences. We identify two major groups, Groups I and II, which are subsequently divided into 11 and 4 subclades, respectively. Y. pestis appears to be maintained in several geographically separate subpopulations. There is also evidence for multiple long distance transfers of Y. pestis, likely human mediated. Such transfers have resulted in the reintroduction and establishment of plague in the port city of Mahajanga, where there is evidence for multiple transfers both from and to the central highlands.

Conclusions/Significance

The maintenance and spread of Y. pestis in Madagascar is a dynamic and highly active process that relies on the natural cycle between the primary host, the black rat, and its flea vectors as well as human activity.

Plague, caused by the bacterium Yersinia pestis, has been a problem in Madagascar since it was introduced in 1898. It mainly affects the central highlands, but also has caused several large outbreaks in the port city of Mahajanga, after it was reintroduced there in the 1990s. Despite its prevalence, the genetic diversity and related geographic distribution of different genetic groups of Y. pestis in Madagascar has been difficult to study due to the great genetic similarity among isolates. We subtyped a set of Malagasy isolates and identified two major genetic groups that were subsequently divided into 11 and 4 subgroups, respectively. Y. pestis appears to be maintained in several geographically separate subpopulations. There is also evidence for multiple long distance transfers of Y. pestis, likely human mediated. Such transfers have resulted in the reintroduction and establishment of plague in the port city of Mahajanga where there is evidence for multiple transfers both from and to the central highlands. The maintenance and spread of Y. pestis in Madagascar is a dynamic and highly active process that relies on the natural cycle between the primary host, the black rat, and its flea vectors as well as human activity.

Phylogeography of Francisella tularensis ssp. holarctica in France

AIM

To investigate the phylogeography of French Francisella tularensis ssp. holarctica isolates.

METHODS AND RESULTS

  Canonical SNPs and MLVA were used to genotype 103 French F. tularensis ssp. holarctica isolates. We confirmed the presence of one subclade, the central and western European group (B.Br.FTNF002-00), and identified four major MLVA genotypes with no obvious geographical differentiation.

CONCLUSIONS

  The lack of geographical resolution among MLVA genotypes suggests rapid dispersal, convergent evolution or a combination of the two.

SIGNIFICANCE AND IMPACT OF THE STUDY

  This study expands knowledge of the phylogeography of one of the two dominant European F. tularensis ssp. holarctica subclades and illustrates the need for additional SNP discovery within this subclade.