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

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

Genomic characterization of Francisella tularensis and other diverse Francisella species from complex samples.

Francisella tularensis, the bacterium that causes the zoonosis tularemia, and its genetic near neighbor species, can be difficult or impossible to cultivate from complex samples. Thus, there is a lack of genomic information for these species that has, among other things, limited the development of robust detection assays for F. tularensis that are both specific and sensitive. The objective of this study was to develop and validate approaches to capture, enrich, sequence, and analyze Francisella DNA present in DNA extracts generated from complex samples. RNA capture probes were designed based upon the known pan genome of F. tularensis and other diverse species in the family Francisellaceae. Probes that targeted genomic regions also present in non-Francisellaceae species were excluded, and probes specific to particular Francisella species or phylogenetic clades were identified. The capture-enrichment system was then applied to diverse, complex DNA extracts containing low-level Francisella DNA, including human clinical tularemia samples, environmental samples (i.e., animal tissue and air filters), and whole ticks/tick cell lines, which was followed by sequencing of the enriched samples. Analysis of the resulting data facilitated rigorous and unambiguous confirmation of the detection of F. tularensis or other Francisella species in complex samples, identification of mixtures of different Francisella species in the same sample, analysis of gene content (e.g., known virulence and antimicrobial resistance loci), and high-resolution whole genome-based genotyping. The benefits of this capture-enrichment system include: even very low target DNA can be amplified; it is culture-independent, reducing exposure for research and/or clinical personnel and allowing genomic information to be obtained from samples that do not yield isolates; and the resulting comprehensive data not only provide robust means to confirm the presence of a target species in a sample, but also can provide data useful for source attribution, which is important from a genomic epidemiology perspective.

Co-Occurrence of Francisella, Spotted Fever Group Rickettsia, and Midichloria in Avian-Associated Hyalomma rufipes.

The migratory behavior of wild birds contributes to the geographical spread of ticks and their microorganisms. In this study, we aimed to investigate the dispersal and co-occurrence of Francisella and spotted fever group Rickettsia (SFGR) in ticks infesting birds migrating northward in the African-Western Palaearctic region (AWPR). Birds were trapped with mist nests across the Mediterranean basin during the 2014 and 2015 spring migration. In total, 575 ticks were collected from 244 birds. We screened the ticks for the species Francisella tularensis, the genus Francisella, and SFGR by microfluidic real-time PCR. Confirmatory analyses and metagenomic sequencing were performed on tick samples that putatively tested positive for F. tularensis during initial screenings. Hyalomma rufipes was the most common tick species and had a high prevalence of Francisella, including co-occurrence of Francisella and SFGR. Metagenomic analysis of total DNA extracted from two H. rufipes confirmed the presence of Francisella, Rickettsia, and Midichloria. Average nucleotide identity and phylogenetic inference indicated the highest identity of the metagenome-assembled genomes to a Francisella-like endosymbiont (FLE), Rickettsia aeschlimannii, and Midichloria mitochondrii. The results of this study suggest that (i) FLE- and SFGR-containing ticks are dispersed by northbound migratory birds in the AWPR, (ii) H. rufipes likely is not involved in transmission of F. tularensis in the AWPR, and (iii) a dual endosymbiosis of FLEs and Midichloria may support some of the nutritional requirements of H. rufipes.

Proteomic Signatures of Antimicrobial Resistance in Yersinia pestis and Francisella tularensis.

Antimicrobial resistance (AMR) is a well-recognized, widespread, and growing issue of concern. With increasing incidence of AMR, the ability to respond quickly to infection with or exposure to an AMR pathogen is critical. Approaches that could accurately and more quickly identify whether a pathogen is AMR also are needed to more rapidly respond to existing and emerging biological threats. We examined proteins associated with paired AMR and antimicrobial susceptible (AMS) strains of Yersinia pestis and Francisella tularensis, causative agents of the diseases plague and tularemia, respectively, to identify whether potential existed to use proteins as signatures of AMR. We found that protein expression was significantly impacted by AMR status. Antimicrobial resistance-conferring proteins were expressed even in the absence of antibiotics in growth media, and the abundance of 10-20% of cellular proteins beyond those that directly confer AMR also were significantly changed in both Y. pestis and F. tularensis. Most strikingly, the abundance of proteins involved in specific metabolic pathways and biological functions was altered in all AMR strains examined, independent of species, resistance mechanism, and affected cellular antimicrobial target. We have identified features that distinguish between AMR and AMS strains, including a subset of features shared across species with different resistance mechanisms, which suggest shared biological signatures of resistance. These features could form the basis of novel approaches to identify AMR phenotypes in unknown strains.

Transmission of Antimicrobial Resistant Yersinia pestis During a Pneumonic Plague Outbreak.

BACKGROUND: Pneumonic plague (PP), caused by Yersinia pestis, is the most feared clinical form of plague due to its rapid lethality and potential to cause outbreaks. PP outbreaks are now rare due to antimicrobial therapy. METHODS: A PP outbreak in Madagascar involving transmission of a Y. pestis strain resistant to streptomycin, the current recommended first-line treatment in Madagascar, was retrospectively characterized using epidemiology, clinical diagnostics, molecular characterization, and animal studies. RESULTS: The outbreak occurred in February 2013 in the Faratsiho district of Madagascar and involved 22 cases, including 3 untreated fatalities. The 19 other cases participated in funeral practices for the fatal cases and fully recovered after combination antimicrobial therapy: intramuscular streptomycin followed by oral co-trimoxazole. The Y. pestis strain that circulated during this outbreak is resistant to streptomycin resulting from a spontaneous point mutation in the 30S ribosomal protein S12 (rpsL) gene. This same mutation causes streptomycin resistance in 2 unrelated Y. pestis strains, one isolated from a fatal PP case in a different region of Madagascar in 1987 and another isolated from a fatal PP case in China in 1996, documenting this mutation has occurred independently at least 3 times in Y. pestis. Laboratory experiments revealed this mutation has no detectable impact on fitness or virulence, and revertants to wild-type are rare in other species containing it, suggesting Y. pestis strains containing it could persist in the environment. CONCLUSIONS: Unique antimicrobial resistant (AMR) strains of Y. pestis continue to arise in Madagascar and can be transmitted during PP outbreaks.

Reorganized Genomic Taxonomy of Francisellaceae Enables Design of Robust Environmental PCR Assays for Detection of Francisella tularensis.

In recent years, an increasing diversity of species has been recognized within the family Francisellaceae. Unfortunately, novel isolates are sometimes misnamed in initial publications or multiple sources propose different nomenclature for genetically highly similar isolates. Thus, unstructured and occasionally incorrect information can lead to confusion in the scientific community. Historically, detection of Francisella tularensis in environmental samples has been challenging due to the considerable and unknown genetic diversity within the family, which can result in false positive results. We have assembled a comprehensive collection of genome sequences representing most known Francisellaceae species/strains and restructured them according to a taxonomy that is based on phylogenetic structure. From this structured dataset, we identified a small number of genomic regions unique to F. tularensis that are putatively suitable for specific detection of this pathogen in environmental samples. We designed and validated specific PCR assays based on these genetic regions that can be used for the detection of F. tularensis in environmental samples, such as water and air filters.

Long-Term Survival of Virulent Tularemia Pathogens outside a Host in Conditions That Mimic Natural Aquatic Environments.

Francisella tularensis, the causative agent of the zoonotic disease tularemia, can cause seasonal outbreaks of acute febrile illness in humans with disease peaks in late summer to autumn. Interestingly, its mechanisms for environmental persistence between outbreaks are poorly understood. One hypothesis is that F. tularensis forms biofilms in aquatic environments. We utilized two fully virulent wild-type strains: FSC200 (Francisella tularensis subsp. holarctica) and Schu S4 (Francisella tularensis subsp. tularensis) and three control strains, the attenuated live vaccine strain (LVS; F. tularensis subsp. holarctica), a Schu S4 ΔwbtI mutant that is documented to form biofilms, and the low-virulence strain U112 of the closely related species Francisella novicida Strains were incubated in saline solution (0.9% NaCl) microcosms for 24 weeks at both 4°C and 20°C, whereupon viability and biofilm formation were measured. These temperatures were selected to approximate winter and summer temperatures of fresh water in Scandinavia, respectively. U112 and Schu S4 ΔwbtI formed biofilms, but F. tularensis strains FSC200 and Schu S4 and the LVS did not. All strains exhibited prolonged viability at 4°C compared to 20°C. U112 and FSC200 displayed remarkable long-term persistence at 4°C, with only 1- and 2-fold log reductions, respectively, of viable cells after 24 weeks. Schu S4 exhibited lower survival, yielding no viable cells by week 20. At 24 weeks, cells from FSC200, but not from Schu S4, were still fully virulent in mice. Taken together, these results demonstrate biofilm-independent, long-term survival of pathogenic F. tularensis subsp. holarctica in conditions that mimic overwinter survival in aquatic environments.IMPORTANCE Tularemia, a disease caused by the environmental bacterium Francisella tularensis, is characterized by acute febrile illness. F. tularensis is highly infectious: as few as 10 organisms can cause human disease. Tularemia is not known to be spread from person to person. Rather, all human infections are independently acquired from the environment via the bite of blood-feeding arthropods, ingestion of infected food or water, or inhalation of aerosolized bacteria. Despite the environmental origins of human disease events, the ecological factors governing the long-term persistence of F. tularensis in nature between seasonal human outbreaks are poorly understood. The significance of our research is in identifying conditions that promote long-term survival of fully virulent F. tularensis outside a mammalian host or insect vector. These conditions are similar to those found in natural aquatic environments in winter and provide important new insights on how F. tularensis may persist long-term in the environment.

Phylogenetic Analysis of Francisella tularensis Group A.II Isolates from 5 Patients with Tularemia, Arizona, USA, 2015-2017.

We examined 5 tularemia cases in Arizona, USA, during 2015-2017. All were caused by Francisella tularensis group A.II. Genetically similar isolates were found across large spatial and temporal distances, suggesting that group A.II strains are dispersed across long distances by wind and exhibit low replication rates in the environment.

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.

Phenotypic characterization and whole genome analysis of extended-spectrum beta-lactamase-producing bacteria isolated from dogs in Germany.

Asymptomatic colonization with extended-spectrum beta-lactamase (ESBL) producing Enterobacteriaceae has been described for humans, various mammal species, and birds. Here, antimicrobial resistant bacteria were recovered from dog feces originating in Germany, Kosovo, Afghanistan, Croatia, and Ukraine, with a subset of mostly E. coli isolates obtained from a longitudinal collection over twelve months. In vitro antimicrobial resistance testing revealed various patterns of resistance against single or all investigated beta-lactam antibiotics, with none of the 101 isolates resistant against two tested carbapenem antibiotics. Whole genome sequence analysis revealed bacteria species-specific patterns for 23 antimicrobial resistance coding DNA sequences (CDS) that were unapparent from the in vitro analysis alone. Phylogenetic analysis of single nucleotide polymorphisms (SNP) revealed clonal bacterial isolates originating from different dogs, suggesting transmission between dogs in the same community. However, individual resistant E. coli clones were not detected over a period longer than seven days. Multi locus sequence typing (MLST) of 85 E. coli isolates revealed 31 different sequence types (ST) with an accumulation of ST744 (n = 9), ST10 (n = 8), and ST648 (n = 6), although the world-wide hospital-associated CTX-M beta-lactamase producing ST131 was not detected. Neither the antimicrobial resistance CDSs patterns nor the phylogenetic analysis revealed an epidemiological correlation among the longitudinal isolates collected from a period longer than seven days. No genetic linkage could be associated with the geographic origin of isolates. In conclusion, healthy dogs frequently carry ESBL-producing bacteria, independent to prior treatment, which may be transmitted between individual dogs of the same community. Otherwise, these antimicrobial resistant bacteria share few commonalities, making their presence eerily unpredictable.

Evaluating the long-term persistence of Bacillus spores on common surfaces.

Bacillus spores resist inactivation, but the extent of their persistence on common surfaces is unclear. This work addresses knowledge gaps regarding biothreat agents in the environment to reduce uncertainty in risk assessment models. Studies were conducted to investigate the long-term inactivation of Bacillus anthracis and three commonly used surrogate organisms - B. cereus, B. atrophaeus and B. thuringiensis on three materials: laminate countertop, stainless steel and polystyrene Petri dishes. Viable spores were measured at 1, 30, 90, 196, 304 and 1038 days. Twelve different persistence models were fit to the data using maximum likelihood estimation and compared. The study found that (1) spore inactivation was not log-linear, as commonly modelled; (2) B. thuringiensis counts increased at 24 h on all materials, followed by a subsequent decline; (3) several experiments showed evidence of a 'U' shape, with spore counts apparently decreasing and then increasing between 1 and 304 days; (4) spores on polystyrene showed little inactivation; and (5) the maximum inactivation of 56% was observed for B. atrophaeus spores on steel at 196 days. Over the range of surfaces, time durations and conditions (humidity controlled vs. uncontrolled) examined, B. thuringiensis most closely matched the behaviour of B. anthracis.

Low cost, low tech SNP genotyping tools for resource-limited areas: Plague in Madagascar as a model.

BACKGROUND: Genetic analysis of pathogenic organisms is a useful tool for linking human cases together and/or to potential environmental sources. The resulting data can also provide information on evolutionary patterns within a targeted species and phenotypic traits. However, the instruments often used to generate genotyping data, such as single nucleotide polymorphisms (SNPs), can be expensive and sometimes require advanced technologies to implement. This places many genotyping tools out of reach for laboratories that do not specialize in genetic studies and/or lack the requisite financial and technological resources. To address this issue, we developed a low cost and low tech genotyping system, termed agarose-MAMA, which combines traditional PCR and agarose gel electrophoresis to target phylogenetically informative SNPs. METHODOLOGY/PRINCIPAL FINDINGS: To demonstrate the utility of this approach for generating genotype data in a resource-constrained area (Madagascar), we designed an agarose-MAMA system targeting previously characterized SNPs within Yersinia pestis, the causative agent of plague. We then used this system to genetically type pathogenic strains of Y. pestis in a Malagasy laboratory not specialized in genetic studies, the Institut Pasteur de Madagascar (IPM). We conducted rigorous assay performance validations to assess potential variation introduced by differing research facilities, reagents, and personnel and found no difference in SNP genotyping results. These agarose-MAMA PCR assays are currently employed as an investigative tool at IPM, providing Malagasy researchers a means to improve the value of their plague epidemiological investigations by linking outbreaks to potential sources through genetic characterization of isolates and to improve understanding of disease ecology that may contribute to a long-term control effort. CONCLUSIONS: The success of our study demonstrates that the SNP-based genotyping capacity of laboratories in developing countries can be expanded with manageable financial cost for resource constraint laboratories. This is a practical formula that reduces resource-driven limitations to genetic research and promises to advance global collective knowledge of infectious diseases emanating from resource limited regions of the world.

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.

Pneumonic Plague Transmission, Moramanga, Madagascar, 2015.

During a pneumonic plague outbreak in Moramanga, Madagascar, we identified 4 confirmed, 1 presumptive, and 9 suspected plague case-patients. Human-to-human transmission among close contacts was high (reproductive number 1.44) and the case fatality rate was 71%. Phylogenetic analysis showed that the Yersinia pestis isolates belonged to group q3, different from the previous outbreak.

A Bacillus anthracis Genome Sequence from the Sverdlovsk 1979 Autopsy Specimens.

Anthrax is a zoonotic disease that occurs naturally in wild and domestic animals but has been used by both state-sponsored programs and terrorists as a biological weapon. A Soviet industrial production facility in Sverdlovsk, USSR, proved deficient in 1979 when a plume of spores was accidentally released and resulted in one of the largest known human anthrax outbreaks. In order to understand this outbreak and others, we generated a Bacillus anthracis population genetic database based upon whole-genome analysis to identify all single-nucleotide polymorphisms (SNPs) across a reference genome. Phylogenetic analysis has defined three major clades (A, B, and C), B and C being relatively rare compared to A. The A clade has numerous subclades, including a major polytomy named the trans-Eurasian (TEA) group. The TEA radiation is a dominant evolutionary feature of B. anthracis, with many contemporary populations having resulted from a large spatial dispersal of spores from a single source. Two autopsy specimens from the Sverdlovsk outbreak were deep sequenced to produce draft B. anthracis genomes. This allowed the phylogenetic placement of the Sverdlovsk strain into a clade with two Asian live vaccine strains, including the Russian Tsiankovskii strain. The genome was examined for evidence of drug resistance manipulation or other genetic engineering, but none was found. The Soviet Sverdlovsk strain genome is consistent with a wild-type strain from Russia that had no evidence of genetic manipulation during its industrial production. This work provides insights into the world's largest biological weapons program and provides an extensive B. anthracis phylogenetic reference. IMPORTANCE: The 1979 Russian anthrax outbreak resulted from an industrial accident at the Soviet anthrax spore production facility in the city of Sverdlovsk. Deep genomic sequencing of two autopsy specimens generated a draft genome and phylogenetic placement of the Soviet Sverdlovsk anthrax strain. While it is known that Soviet scientists had genetically manipulated Bacillus anthracis with the potential to evade vaccine prophylaxis and antibiotic therapeutics, there was no genomic evidence of this from the Sverdlovsk production strain genome. The whole-genome SNP genotype of the Sverdlovsk strain was used to precisely identify it and its close relatives in the context of an extensive global B. anthracis strain collection. This genomic identity can now be used for forensic tracking of this weapons material on a global scale and for future anthrax investigations.

High Leptospira Diversity in Animals and Humans Complicates the Search for Common Reservoirs of Human Disease in Rural Ecuador.

BACKGROUND: Leptospirosis is a zoonotic disease responsible for high morbidity around the world, especially in tropical and low income countries. Rats are thought to be the main vector of human leptospirosis in urban settings. However, differences between urban and low-income rural communities provide additional insights into the epidemiology of the disease. METHODOLOGY/PRINCIPAL FINDINGS: Our study was conducted in two low-income rural communities near the coast of Ecuador. We detected and characterized infectious leptospira DNA in a wide variety of samples using new real time quantitative PCR assays and amplicon sequencing. We detected infectious leptospira in a high percentage of febrile patients (14.7%). In contrast to previous studies on leptospirosis risk factors, higher positivity was not found in rats (3.0%) but rather in cows (35.8%) and pigs (21.1%). Six leptospira species were identified (L. borgpetersenii, L kirschnerii, L santarosai, L. interrogans, L noguchii, and an intermediate species within the L. licerasiae and L. wolffii clade) and no significant differences in the species of leptospira present in each animal species was detected (χ2 = 9.89, adj.p-value = 0.27). CONCLUSIONS/SIGNIFICANCE: A large portion of the world's human population lives in low-income, rural communities, however, there is limited information about leptospirosis transmission dynamics in these settings. In these areas, exposure to peridomestic livestock is particularly common and high prevalence of infectious leptospira in cows and pigs suggest that they may be the most important reservoir for human transmission. Genotyping clinical samples show that multiple species of leptospira are involved in human disease. As these genotypes were also detected in samples from a variety of animals, genotype data must be used in conjunction with epidemiological data to provide evidence of transmission and the importance of different potential leptospirosis reservoirs.

Clonality of erythromycin resistance in Francisella tularensis.

OBJECTIVES: We analysed diverse strains of Francisella tularensis subsp. holarctica to assess if its division into biovars I and II is associated with specific mutations previously linked to erythromycin resistance and to determine the distribution of this resistance trait across this subspecies. METHODS: Three-hundred and fourteen F. tularensis subsp. holarctica strains were tested for erythromycin susceptibility and whole-genome sequences for these strains were examined for SNPs in genes previously associated with erythromycin resistance. Each strain was assigned to a global phylogenetic framework using genome-wide canonical SNPs. The contribution of a specific SNP to erythromycin resistance was examined using allelic exchange. The geographical distribution of erythromycin-resistant F. tularensis strains was further investigated by literature search. RESULTS: There was a perfect correlation between biovar II strains (erythromycin resistance) and the phylogenetic group B.12. Only B.12 strains had an A → C SNP at position 2059 in the three copies of the rrl gene. Introducing 2059C into an rrl gene of an erythromycin-susceptible F. tularensis strain resulted in resistance. An additional 1144 erythromycin-resistant strains were identified from the scientific literature, all of them from Eurasia. CONCLUSIONS: Erythromycin resistance in F. tularensis is caused by an A2059C rrl gene mutation, which exhibits a strictly clonal inheritance pattern found only in phylogenetic group B.12. This group is an extremely successful clone, representing the most common type of F. tularensis throughout Eurasia.

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.

Long-range dispersal moved Francisella tularensis into Western Europe from the East.

For many infections transmitting to humans from reservoirs in nature, disease dispersal patterns over space and time are largely unknown. Here, a reversed genomics approach helped us understand disease dispersal and yielded insight into evolution and biological properties of Francisella tularensis, the bacterium causing tularemia. We whole-genome sequenced 67 strains and characterized by single-nucleotide polymorphism assays 138 strains, collected from individuals infected 1947-2012 across Western Europe. We used the data for phylogenetic, population genetic and geographical network analyses. All strains (n=205) belonged to a monophyletic population of recent ancestry not found outside Western Europe. Most strains (n=195) throughout the study area were assigned to a star-like phylogenetic pattern indicating that colonization of Western Europe occurred via clonal expansion. In the East of the study area, strains were more diverse, consistent with a founder population spreading from east to west. The relationship of genetic and geographic distance within the F. tularensis population was complex and indicated multiple long-distance dispersal events. Mutation rate estimates based on year of isolation indicated null rates; in outbreak hotspots only, there was a rate of 0.4 mutations/genome/year. Patterns of nucleotide substitution showed marked AT mutational bias suggestive of genetic drift. These results demonstrate that tularemia has moved from east to west in Europe and that F. tularensis has a biology characterized by long-range geographical dispersal events and mostly slow, but variable, replication rates. The results indicate that mutation-driven evolution, a resting survival phase, genetic drift and long-distance geographical dispersal events have interacted to generate genetic diversity within this species.