Pre-Columbian origins for North American anthrax

Bacillus anthracis strain differentiation based on SNP and VNTR loci

Bacillus anthracis is the anthrax causative agent. For its epidemiology, it is important not only to identify the etiological agent but also to determine the patterns of its evolution and spread. Modern methods of molecular biology make it possible to detect a number of genetic markers suitable for indicating and differentiating the strains of B. anthracis, including the loci arranged as variable number tandem repeats (VNTRs) and SNPs, one nucleotide-sized differences in the DNA sequence of the loci being compared. The objective of the present study was to examine the effectiveness of SNP analysis and PCR amplif ication of VNTR loci combined with the high-resolution amplicon melting analysis for identif ication and differentiation of the anthrax agent strains. In the study, seven strains of B. anthracis obtained from soil samples and animal carcasses were investigated using vaccine strain STI-1 as a reference. For molecular genetic characterization of these bacteria, analysis of 12 SNPs and variability analysis of eight VNTR loci were carried out. To detect the differences between the strains, their PCR product melting points were measured in the presence of the EvaGreen (Sintol, Russia) intercalating dye. For SNP detection, a PCR assay with double TaqMan probes was applied. It was found that the studied virulent strains, except for B. anthracis No. 1 and 3, could not be attributed to any phylogenetic subgroup of the anthrax agents. The proposed method made it possible to differentiate four out of the seven investigated strains. Strains No. 5-7 had identical SNP and HRM prof iles and, as a result, formed a single cluster. Our investigation has conf irmed that the proposed method can be successfully used for preliminary analysis of an epizootic situation in the case of anthrax.

Bacillus anthracis Phylogeography: New Clues From Kazakhstan, Central Asia

This article describes Bacillus anthracis strains isolated in Kazakhstan since the 1950s until year 2016 from sixty-one independent events associated with anthrax in humans and animals. One hundred and fifty-four strains were first genotyped by Multiple Locus VNTR (variable number of tandem repeats) Analysis (MLVA) using 31 VNTR loci. Thirty-five MLVA31 genotypes were resolved, 28 belong to the A1/TEA group, five to A3/Sterne-Ames group, one to A4/Vollum and one to the B clade. This is the first report of the presence of the B-clade in Kazakhstan. The MLVA31 results and epidemiological data were combined to select a subset of seventy-nine representative strains for draft whole genome sequencing (WGS). Strains from Kazakhstan significantly enrich the known phylogeny of the Ames group polytomy, including the description of a new branch closest to the Texas, United States A.Br.Ames sublineage stricto sensu. Three among the seven currently defined branches in the TEA polytomy are present in Kazakhstan, “Tsiankovskii”, “Heroin”, and “Sanitary Technical Institute (STI)”. In particular, strains from the STI lineage are largely predominant in Kazakhstan and introduce numerous deep branching STI sublineages, demonstrating a high geographic correspondence between “STI” and Kazakhstan, Central Asia. This observation is a strong indication that the TEA polytomy emerged after the last political unification of Asian steppes in the fourteenth century of the Common Era. The phylogenetic analysis of the Kazakhstan data and of currently available WGS data of worldwide origin strengthens our understanding of B. anthracis geographic expansions in the past seven centuries.

Molecular characterization of B. anthracis isolates from the anthrax outbreak among cattle in Karnataka, India

Background

Anthrax, a zoonotic disease is caused by the Gram positive bacterium Bacillus anthracis. During January 2013, an anthrax outbreak among cattle was reported in Gundlupet Taluk, neighboring Bandipur National Park and tiger reserve, India. The present study aims at the molecular identification and characterization of 12 B. anthracis isolates from this outbreak by 16S rRNA gene sequencing, screening B. anthracis specific prophages and chromosomal markers, protective antigen (pag) gene and canonical single nucleotide polymorphism (canSNP) analysis to subtype the isolates into one of the twelve globally identified clonal sub-lineages of B. anthracis.

Results

These isolates had identical 16S rDNA nucleotide sequences with B. anthracis specific dual peaks showing mixed base pair R (G/A) at position 1139 with visual inspection while the automated basecaller software indicated a G. Alternatively the nucleotide A at 1146 position was indicative of the 16S rDNA type 7. Multiple sequence alignment with additional 170 (16S rDNA) sequences of B. cereus sensu lato group from GenBank database revealed 28 new 16S types in addition to eleven 16S types reported earlier. The twelve B. anthracis isolates were found to harbor the four B. anthracis specific prophages (lambdaBa01, lambdaBa02, lambdaBa03, and lambdaBa04) along with its four specific loci markers (dhp 61.183, dhp 77.002, dhp 73.019, and dhp 73.017). The pag gene sequencing identified the isolates as protective antigen (PA) genotype I with phenylalanine-proline-alanine phenotype (FPA phenotype). However, sequence clustering with additional 34 pag sequences from GenBank revealed two additional missense mutations at nucleotide positions 196 bp and 869 bp of the 2294 bp pag sequence among the 5 B. cereus strains with pXO1 like plasmids. The canSNP analysis showed that the isolates belong to A.Br.Aust94 sub-lineage that is distributed geographically in countries of Asia, Africa, Europe and Australia.

Conclusions

The analysis of 16S rDNA sequences reiterated the earlier findings that visual inspection of electropherogram for position 1139 having nucleotide R could be used for B. anthracis identification and not the consensus sequence from base caller. The canSNP results indicated that the anthrax outbreak among cattle was caused by B. anthracis of A.Br.Aust94 sub-lineage.

A classification framework for Bacillus anthracis defined by global genomic structure

Bacillus anthracis, the causative agent of anthrax, is a considerable global health threat affecting wildlife, livestock, and the general public. In this study, whole-genome sequence analysis of over 350 B. anthracis isolates was used to establish a new high-resolution global genotyping framework that is both biogeographically informative and compatible with multiple genomic assays. The data presented in this study shed new light on the diverse global dissemination of this species and indicate that many lineages may be uniquely suited to the geographic regions in which they are found. In addition, we demonstrate that plasmid genomic structure for this species is largely consistent with chromosomal population structure, suggesting vertical inheritance in this bacterium has contributed to its evolutionary persistence. This classification methodology is the first based on population genomic structure for this species and has potential use for local and broader institutions seeking to understand both disease outbreak origins and recent introductions. In addition, we provide access to a newly developed genotyping script as well as the full whole-genome sequence analyses output for this study, allowing future studies to rapidly employ and append their data in the context of this global collection. This framework may act as a powerful tool for public health agencies, wildlife disease laboratories, and researchers seeking to utilize and expand this classification scheme for further investigations into B. anthracis evolution.

Predicting the Geographic Distribution of the Bacillus anthracis A1.a/Western North American Sub-Lineage for the Continental United States: New Outbreaks, New Genotypes, and New Climate Data

Bacillus anthracis, the causative pathogen of anthrax, is a spore-forming, environmentally maintained bacterium that continues to be a veterinary health problem with outbreaks occurring primarily in wildlife and livestock. Globally, the genetic populations of B. anthracis include multiple lineages, and each may have different ecological requirements and geographical distributions. It is, therefore, essential to identify environmental associations within lineages to predict geographical distributions and risk areas with improved accuracy. Here, we model the ecological niche and predict the geography of the most widespread sublineage of B. anthracis in the continental United States using updated MERRA-derived (Modern Era Retrospective analysis for Research and Applications; the NASA atmospheric data reanalysis of satellite information with multiple data products) bioclimate variables (i.e., MERRAclim data) and updated soil variables. We filter the occurrence data associated with the A1.a/Western North American sub-lineage of B. anthracis from historical anthrax outbreaks using the multiple-locus variable-number tandem repeat system. In addition, we also incorporate recent cases associated with B. anthracis A1.a sub-lineage from 2008 to 2012 in Montana, Colorado, and Texas. Our results provide the predicted distribution of the A1.a sub-lineage of B. anthracis for the United States with better predictive accuracy and higher spatial resolution than previous estimates. Our prediction serves as an improved disease risk map to better inform anthrax surveillance and control in the United States, particularly the Dakotas and Montana where this sub-lineage is persistent.

Insights from Bacillus anthracis strains isolated from permafrost in the tundra zone of Russia

This article describes Bacillus anthracis strains isolated during an outbreak of anthrax on the Yamal Peninsula in the summer of 2016 and independently in Yakutia in 2015. A common feature of these strains is their conservation in permafrost, from which they were extracted either due to the thawing of permafrost (Yamal strains) or as the result of paleontological excavations (Yakut strains). All strains isolated on the Yamal share an identical genotype belonging to lineage B.Br.001/002, pointing to a common source of infection in a territory over 250 km in length. In contrast, during the excavations in Yakutia, three genetically different strains were recovered from a single pit. One strain belongs to B.Br.001/002, and whole genome sequence analysis showed that it is most closely related to the Yamal strains in spite of the remoteness of Yamal from Yakutia. The two other strains contribute to two different branches of A.Br.008/011, one of the remarkable polytomies described so far in the B. anthracis species. The geographic distribution of the strains belonging to A.Br.008/011 is suggesting that the polytomy emerged in the thirteenth century, in combination with the constitution of a unified Mongol empire extending from China to Eastern Europe. We propose an evolutionary model for B. anthracis recent evolution in which the B lineage spread throughout Eurasia and was subsequently replaced by the A lineage except in some geographically isolated areas.

Insights from Bacillus anthracis strains isolated from permafrost in the tundra zone of Russia

This article describes Bacillus anthracis strains isolated during an outbreak of anthrax on the Yamal Peninsula in the summer of 2016 and independently in Yakutia in 2015. A common feature of these strains is their conservation in permafrost, from which they were extracted either due to the thawing of permafrost (Yamal strains) or as the result of paleontological excavations (Yakut strains). All strains isolated on the Yamal share an identical genotype belonging to lineage B.Br.001/002, pointing to a common source of infection in a territory over 250 km in length. In contrast, during the excavations in Yakutia, three genetically different strains were recovered from a single pit. One strain belongs to B.Br.001/002, and whole genome sequence analysis showed that it is most closely related to the Yamal strains in spite of the remoteness of Yamal from Yakutia. The two other strains contribute to two different branches of A.Br.008/011, one of the remarkable polytomies described so far in the B. anthracis species. The geographic distribution of the strains belonging to A.Br.008/011 is suggesting that the polytomy emerged in the thirteenth century, in combination with the constitution of a unified Mongol empire extending from China to Eastern Europe. We propose an evolutionary model for B. anthracis recent evolution in which the B lineage spread throughout Eurasia and was subsequently replaced by the A lineage except in some geographically isolated areas.

Comparison of spatiotemporal patterns of historic natural Anthrax outbreaks in Minnesota and Kazakhstan

Disease spread in populations is a consequence of the interaction between host, pathogen, and environment, i.e. the epidemiological triad. Yet the influences of each triad component may vary dramatically for different settings. Comparison of environmental, demographic, socio-economic, and historical backgrounds may support tailoring site-specific control measures. Because of the long-term survival of Bacillus anthracis, Anthrax is a suitable example for studying the influence of triad components in different endemic settings. We compared the spatiotemporal patterns of historic animal Anthrax records in two endemic areas, located at northern latitudes in the western and eastern hemispheres. Our goal was to compare the spatiotemporal patterns in Anthrax progression, intensity, direction, and recurrence (disease hot spots), in relation to epidemiological factors and potential trigger events. Reported animal cases in Minnesota, USA (n = 289 cases between 1912 and 2014) and Kazakhstan (n = 3,997 cases between 1933 and 2014) were analyzed using the spatiotemporal directionality test and the spatial scan statistic. Over the last century Anthrax occurrence in Minnesota was sporadic whereas Kazakhstan experienced a long-term epidemic. Nevertheless, the seasonality was comparable between sites, with a peak in August. Declining number of cases at both sites was attributed to vaccination and control measures. The spatiotemporal directionality test detected a relative northeastern directionality in disease spread for long-term trends in Minnesota, whereas a southwestern directionality was observed in Kazakhstan. In terms of recurrence, the maximum timespans between cases at the same location were 55 and 60 years for Minnesota and Kazakhstan, respectively. Disease hotspots were recognized in both settings, with spatially overlapping clusters years apart. Distribution of the spatiotemporal cluster radii between study sites supported suggestion of site-specific control zones. Spatiotemporal patterns of Anthrax occurrence in both endemic regions were attributed to multiple potential trigger events including major river floods, changes in land use, agriculture, and susceptible livestock populations. Results here help to understand the long-term epidemiological dynamics of Anthrax while providing suggestions to the design and implementation of prevention and control programs, in endemic settings.

Burkholderia pseudomallei distribution in Australasia is linked to paleogeographic and anthropogenic history

Burkholderia pseudomallei is the environmental bacillus that causes melioidosis; a disease clinically significant in Australia and Southeast Asia but emerging in tropical and sub-tropical regions around the globe. Previous studies have placed the ancestral population of the organism in Australia with a single lineage disseminated to Southeast Asia. We have previously characterized B. pseudomallei isolates from New Guinea and the Torres Strait archipelago; remote regions that share paleogeographic ties with Australia. These studies identified regional biogeographical boundaries. In this study, we utilize whole-genome sequencing to reconstruct ancient evolutionary relationships and ascertain correlations between paleogeography and present-day distribution of this bacterium in Australasia. Our results indicate that B. pseudomallei from New Guinea fall into a single clade within the Australian population. Furthermore, clades from New Guinea are region-specific; an observation possibly linked to limited recent anthropogenic influence in comparison to mainland Australia and Southeast Asia. Isolates from the Torres Strait archipelago were distinct yet scattered among those from mainland Australia. These results provide evidence that the New Guinean and Torres Strait lineages may be remnants of an ancient portion of the Australian population. Rising sea levels isolated New Guinea and the Torres Strait Islands from each other and the Australian mainland, and may have allowed long-term isolated evolution of these lineages, providing support for a theory of microbial biogeography congruent with that of macro flora and fauna. Moreover, these findings indicate that contemporary microbial biogeography theories should consider recent and ongoing impacts of globalisation and human activity.

Phylogenetic conservatism of thermal traits explains dispersal limitation and genomic differentiation of Streptomyces sister-taxa

The latitudinal diversity gradient is a pattern of biogeography observed broadly in plants and animals but largely undocumented in terrestrial microbial systems. Although patterns of microbial biogeography across broad taxonomic scales have been described in a range of contexts, the mechanisms that generate biogeographic patterns between closely related taxa remain incompletely characterized. Adaptive processes are a major driver of microbial biogeography, but there is less understanding of how microbial biogeography and diversification are shaped by dispersal limitation and drift. We recently described a latitudinal diversity gradient of species richness and intraspecific genetic diversity in Streptomyces by using a geographically explicit culture collection. Within this geographically explicit culture collection, we have identified Streptomyces sister-taxa whose geographic distribution is delimited by latitude. These sister-taxa differ in geographic distribution, genomic diversity, and ecological traits despite having nearly identical SSU rRNA gene sequences. Comparative genomic analysis reveals genomic differentiation of these sister-taxa consistent with restricted gene flow across latitude. Furthermore, we show phylogenetic conservatism of thermal traits between the sister-taxa suggesting that thermal trait adaptation limits dispersal and gene flow across climate regimes as defined by latitude. Such phylogenetic conservatism of thermal traits is commonly associated with latitudinal diversity gradients for plants and animals. These data provide further support for the hypothesis that the Streptomyces latitudinal diversity gradient was formed as a result of historical demographic processes defined by dispersal limitation and driven by paleoclimate dynamics.

Unexpected Relations of Historical Anthrax Strain

In 1998, it was claimed that an 80-year-old glass tube intentionally filled with Bacillus anthracis and embedded in a sugar lump as a WWI biological weapon still contained viable spores. Today, genome sequencing of three colonies isolated in 1998 and subjected to phylogenetic analysis surprisingly identified a well-known B. anthracis reference strain isolated in the United States in 1981, pointing to accidental laboratory contamination.

Next-generation sequencing and subsequent phylogenetic analyses are useful and reliable tools for the classification of recent and historical samples. The reliability of sequences obtained and bioinformatic algorithms has increased in recent years, and research has uncovered the identity of a presumed bioweapon agent as a contaminant.

Genome Surfing As Driver of Microbial Genomic Diversity

Historical changes in population size, such as those caused by demographic range expansions, can produce nonadaptive changes in genomic diversity through mechanisms such as gene surfing. We propose that demographic range expansion of a microbial population capable of horizontal gene exchange can result in genome surfing, a mechanism that can cause widespread increase in the pan-genome frequency of genes acquired by horizontal gene exchange. We explain that patterns of genetic diversity within Streptomyces are consistent with genome surfing, and we describe several predictions for testing this hypothesis both in Streptomyces and in other microorganisms.

Genome Sequence of the Soviet/Russian Bacillus anthracis Vaccine Strain 55-VNIIVViM

Bacillus anthracis strain 55-VNIIVViM is a live-attenuated nonencapsulated Soviet/Russian veterinary anthrax vaccine strain. We report here the genome of 55-VNIIVViM and confirm its phylogenetic placement in the global population structure of B. anthracis.

Comparative genomics of Bacillus anthracis from the wool industry highlights polymorphisms of lineage A.Br.Vollum

BACKGROUND

With the advent of affordable next-generation sequencing (NGS) technologies, major progress has been made in the understanding of the population structure and evolution of the B. anthracis species. Here we report the use of whole genome sequencing and computer-based comparative analyses to characterize six strains belonging to the A.Br.Vollum lineage. These strains were isolated in Switzerland, in 1981, during iterative cases of anthrax involving workers in a textile plant processing cashmere wool from the Indian subcontinent.

RESULTS

We took advantage of the hundreds of currently available B. anthracis genomes in public databases, to investigate the genetic diversity existing within the A.Br.Vollum lineage and to position the six Swiss isolates into the worldwide B. anthracis phylogeny. Thirty additional genomes related to the A.Br.Vollum group were identified by whole-genome single nucleotide polymorphism (SNP) analysis, including two strains forming a new evolutionary branch at the basis of the A.Br.Vollum lineage. This new phylogenetic lineage (termed A.Br.H9401) splits off the branch leading to the A.Br.Vollum group soon after its divergence to the other lineages of the major A clade (i.e. 6 SNPs). The available dataset of A.Br.Vollum genomes were resolved into 2 distinct groups. Isolates from the Swiss wool processing facility clustered together with two strains from Pakistan and one strain of unknown origin isolated from yarn. They were clearly differentiated (69 SNPs) from the twenty-five other A.Br.Vollum strains located on the branch leading to the terminal reference strain A0488 of the lineage. Novel analytic assays specific to these new subgroups were developed for the purpose of rapid molecular epidemiology.

CONCLUSIONS

Whole genome SNP surveys greatly expand upon our knowledge on the sub-structure of the A.Br.Vollum lineage. Possible origin and route of spread of this lineage worldwide are discussed.

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.

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.

Latitude delineates patterns of biogeography in terrestrial Streptomyces

The biogeography of Streptomyces was examined at regional spatial scales to identify factors that govern patterns of microbial diversity. Streptomyces are spore forming filamentous bacteria which are widespread in soil. Streptomyces strains were isolated from perennial grass habitats sampled across a spatial scale of more than 6000 km. Previous analysis of this geographically explicit culture collection provided evidence for a latitudinal diversity gradient in Streptomyces species. Here the hypothesis that this latitudinal diversity gradient is a result of evolutionary dynamics associated with historical demographic processes was evaluated. Historical demographic phenomena have genetic consequences that can be evaluated through analysis of population genetics. Population genetic approaches were applied to analyze population structure in six of the most numerically abundant and geographically widespread Streptomyces phylogroups from our culture collection. Streptomyces population structure varied at regional spatial scales, and allelic diversity correlated with geographic distance. In addition, allelic diversity and gene flow are partitioned by latitude. Finally, it was found that nucleotide diversity within phylogroups was negatively correlated with latitude. These results indicate that phylogroup diversification is constrained by dispersal limitation at regional spatial scales, and they are consistent with the hypothesis that historical demographic processes have influenced the contemporary biogeography of Streptomyces.

A Latitudinal Diversity Gradient in Terrestrial Bacteria of the Genus Streptomyces

We show that Streptomyces biogeography in soils across North America is influenced by the regional diversification of microorganisms due to dispersal limitation and genetic drift. Streptomyces spp. form desiccation-resistant spores, which can be dispersed on the wind, allowing for a strong test of whether dispersal limitation governs patterns of terrestrial microbial diversity. We employed an approach that has high sensitivity for determining the effects of genetic drift. Specifically, we examined the genetic diversity and phylogeography of physiologically similar Streptomyces strains isolated from geographically distributed yet ecologically similar habitats. We found that Streptomyces beta diversity scales with geographic distance and both beta diversity and phylogenetic diversity manifest in a latitudinal diversity gradient. This pattern of Streptomyces biogeography resembles patterns seen for diverse species of plants and animals, and we therefore evaluated these data in the context of ecological and evolutionary hypotheses proposed to explain latitudinal diversity gradients. The data are consistent with the hypothesis that niche conservatism limits dispersal, and historical patterns of glaciation have limited the time for speciation in higher-latitude sites. Most notably, higher-latitude sites have lower phylogenetic diversity, higher phylogenetic clustering, and evidence of range expansion from lower latitudes. In addition, patterns of beta diversity partition with respect to the glacial history of sites. Hence, the data support the hypothesis that extant patterns of Streptomyces biogeography have been driven by historical patterns of glaciation and are the result of demographic range expansion, dispersal limitation, and regional diversification due to drift.

Biogeographic patterns provide insight into the evolutionary and ecological processes that govern biodiversity. However, the evolutionary and ecological processes that govern terrestrial microbial diversity remain poorly characterized. We evaluated the biogeography of the genus Streptomyces to show that the diversity of terrestrial bacteria is governed by many of the same processes that govern the diversity of many plant and animal species. While bacteria of the genus Streptomyces are a preeminent source of antibiotics, their evolutionary history, biogeography, and biodiversity remain poorly characterized. The observations we describe provide insight into the drivers of Streptomyces biodiversity and the processes that underlie microbial diversification in terrestrial habitats.

Comparison of French and Worldwide Bacillus anthracis Strains Favors a Recent, Post-Columbian Origin of the Predominant North-American Clade

Background

Bacillus anthracis, the highly dangerous zoonotic bacterial pathogen species is currently composed of three genetic groups, called A, B and C. Group A is represented worldwide whereas group B is present essentially in Western Europe and Southern Africa. Only three strains from group C have been reported. This knowledge is derived from the genotyping of more than 2000 strains collected worldwide. Strains from both group A and group B are present in France. Previous investigations showed that the majority of sporadic French strains belong to the so-called A.Br.011/009 group A clade and define a very remarkable polytomy with six branches. Here we explore the significance of this polytomy by comparing the French B. anthracis lineages to worldwide lineages. We take advantage of whole genome sequence data previously determined for 122 French strains and 45 strains of various origins.

Results

A total of 6690 SNPs was identified among the available dataset and used to draw the phylogeny. The phylogeny of the French B group strains which belongs to B.Br.CNEVA indicates an expansion from the south-east part of France (the Alps) towards the south-west (Massif-Central and Pyrenees). The relatively small group A strains belonging to A.Br.001/002 results from at least two independent introductions. Strikingly, the data clearly demonstrates that the currently predominant B. anthracis lineage in North America, called WNA for Western North American, is derived from one branch of the A.Br.011/009 polytomy predominant in France.

Conclusions/Significance

The present work extends the range of observed substitution rate heterogeneity within B. anthracis, in agreement with its ecology and in contrast with some other pathogens. The population structure of the six branches A.Br.011/009 polytomy identified in France, diversity of branch length, and comparison with the WNA lineage, suggests that WNA is of post-Columbian and west European origin, with France as a likely source. Furthermore, it is tempting to speculate that the polytomy’s most recent common ancestor -MRCA- dates back to the Hundred Years' war between France and England started in the mid-fourteenth century. These events were associated in France with deadly epidemics and major economic and social changes.

Whole Genome Analysis of Injectional Anthrax Identifies Two Disease Clusters Spanning More Than 13 Years

Background

Anthrax is a rare disease in humans but elicits great public fear because of its past use as an agent of bioterrorism. Injectional anthrax has been occurring sporadically for more than ten years in heroin consumers across multiple European countries and this outbreak has been difficult to trace back to a source.

Methods

We took a molecular epidemiological approach in understanding this disease outbreak, including whole genome sequencing of Bacillus anthracis isolates from the anthrax victims. We also screened two large strain repositories for closely related strains to provide context to the outbreak.

Findings

Analyzing 60 Bacillus anthracis isolates associated with injectional anthrax cases and closely related reference strains, we identified 1071 Single Nucleotide Polymorphisms (SNPs). The synapomorphic SNPs (350) were used to reconstruct phylogenetic relationships, infer likely epidemiological sources and explore the dynamics of evolving pathogen populations. Injectional anthrax genomes separated into two tight clusters: one group was exclusively associated with the 2009–10 outbreak and located primarily in Scotland, whereas the second comprised more recent (2012–13) cases but also a single Norwegian case from 2000.

Interpretation

Genome-based differentiation of injectional anthrax isolates argues for at least two separate disease events spanning > 12 years. The genomic similarity of the two clusters makes it likely that they are caused by separate contamination events originating from the same geographic region and perhaps the same site of drug manufacturing or processing. Pathogen diversity within single patients challenges assumptions concerning population dynamics of infecting B. anthracis and host defensive barriers for injectional anthrax.

Funding

This work was supported by the United States Department of Homeland Security grant no. HSHQDC-10-C-00,139 and via a binational cooperative agreement between the United States Government and the Government of Germany. This work was supported by funds from the German Ministry of Defense (Sonderforschungsprojekt 25Z1-S-431,214). Support for sequencing was also obtained from Illumina, Inc. These sources had no role in the data generation or interpretation, and had not role in the manuscript preparation.

Panel 1: Research in Context Systematic Review

We searched PubMed for any article published before Jun. 17, 2015, with the terms “Bacillus anthracis” and “heroin”, or “injectional anthrax”. Other than our previously published work (Price et al., 2012), we found no other relevant studies on elucidating the global phylogenetic relationships of B. anthracis strains associated with injectional anthrax caused by recreational heroin consumption of spore-contaminated drug. There were, however, publically available genome sequences of two strains involved (Price et al., 2012, Grunow et al., 2013) and the draft genome sequence of Bacillus anthracis UR-1, isolated from a German heroin user (Ruckert et al., 2012) with only limited information on the genotyping of closely related strains (Price et al., 2012, Grunow et al., 2013).

Lay Person Interpretation

Injectional anthrax has been plaguing heroin drug users across Europe for more than 10 years. In order to better understand this outbreak, we assessed genomic relationships of all available injectional anthrax strains from four countries spanning a > 12 year period. Very few differences were identified using genome-based analysis, but these differentiated the isolates into two distinct clusters. This strongly supports a hypothesis of at least two separate anthrax spore contamination events perhaps during the drug production processes. Identification of two events would not have been possible from standard epidemiological analysis. These comprehensive data will be invaluable for classifying future injectional anthrax isolates and for future geographic attribution.

•

Whole genome sequences of injectional anthrax B. anthracis isolates fall in two tight but distinct genomic clusters.

•

The distinct genomic clusters are consistent with two or more disease events that overlap in time and space.

•

Defining pathogen clusters will lead to better public health responses to difficult to track disease outbreaks.

Whole Genome-Sequencing and Phylogenetic Analysis of a Historical Collection of Bacillus anthracis Strains from Danish Cattle

Bacillus anthracis, the causative agent of anthrax, is known as one of the most genetically monomorphic species. Canonical single-nucleotide polymorphism (SNP) typing and whole-genome sequencing were used to investigate the molecular diversity of eleven B. anthracis strains isolated from cattle in Denmark between 1935 and 1988. Danish strains were assigned into five canSNP groups or lineages, i.e. A.Br.001/002 (n = 4), A.Br.Ames (n = 2), A.Br.008/011 (n = 2), A.Br.005/006 (n = 2) and A.Br.Aust94 (n = 1). The match with the A.Br.Ames lineage is of particular interest as the occurrence of such lineage in Europe is demonstrated for the first time, filling an historical gap within the phylogeography of the lineage. Comparative genome analyses of these strains with 41 isolates from other parts of the world revealed that the two Danish A.Br.008/011 strains were related to the heroin-associated strains responsible for outbreaks of injection anthrax in drug users in Europe. Eight novel diagnostic SNPs that specifically discriminate the different sub-groups of Danish strains were identified and developed into PCR-based genotyping assays.

Phylogeography of Bacillus anthracis in the country of Georgia shows evidence of population structuring and is dissimilar to other regional genotypes

Sequence analyses and subtyping of Bacillus anthracis strains from Georgia reveal a single distinct lineage (Aust94) that is ecologically established. Phylogeographic analysis and comparisons to a global collection reveals a clade that is mostly restricted to Georgia. Within this clade, many groups are found around the country, however at least one subclade is only found in the eastern part. This pattern suggests that dispersal into and out of Georgia has been rare and despite historical dispersion within the country, for at least for one lineage, current spread is limited.

Genetic diversity of Bacillus anthracis in Europe: genotyping methods in forensic and epidemiologic investigations

Bacillus anthracis, the etiological agent of anthrax, a zoonosis relatively common throughout the world, can be used as an agent of bioterrorism. In naturally occurring outbreaks and in criminal release of this pathogen, a fast and accurate diagnosis is crucial to an effective response. Microbiological forensics and epidemiologic investigations increasingly rely on molecular markers, such as polymorphisms in DNA sequence, to obtain reliable information regarding the identification or source of a suspicious strain. Over the past decade, significant research efforts have been undertaken to develop genotyping methods with increased power to differentiate B. anthracis strains. A growing number of DNA signatures have been identified and used to survey B. anthracis diversity in nature, leading to rapid advances in our understanding of the global population of this pathogen. This article provides an overview of the different phylogenetic subgroups distributed across the world, with a particular focus on Europe. Updated information on the anthrax situation in Europe is reported. A brief description of some of the work in progress in the work package 5.1 of the AniBioThreat project is also presented, including (1) the development of a robust typing tool based on a suspension array technology and multiplexed single nucleotide polymorphisms scoring and (2) the typing of a collection of DNA from European isolates exchanged between the partners of the project. The know-how acquired will contribute to improving the EU's ability to react rapidly when the identity and real origin of a strain need to be established.

Genotyping of French Bacillus anthracis strains based on 31-loci multi locus VNTR analysis: epidemiology, marker evaluation, and update of the internet genotype database

Background

Bacillus anthracis is known to have low genetic variability. In spite of this lack of diversity, multiple-locus variable-number tandem repeat (VNTR) analysis (MLVA) and single nucleotide polymorphisms (SNPs) including the canonical SNPs assay (canSNP) have proved to be highly effective to differentiate strains. Five different MLVA schemes based on a collection of 31 VNTR loci (MLVA8, MLVA15, MLVA20, MLVA25 and MLVA31) with increased resolving power have been described.

Results

MLVA31 was applied to characterize the French National Reference Laboratory collection. The total collection of 130 strains is resolved in 35 genotypes. The 119 veterinary and environmental strains collection in France were resolved into 26 genotypes belonging to three canSNP lineages and four MLVA clonal complexes (CCs) with particular geographical clustering. A subset of seven loci (MLVA7) is proposed to constitute a first line assay. The loci are compatible with moderate resolution equipment such as agarose gel electrophoresis and show a good congruence value with MLVA31. The associated MLVA and SNP data was imported together with published genotyping data by taking advantage of major enhancements to the MLVAbank software and web site.

Conclusions

The present report provides a wide coverage of the genetic diversity of naturally occurring B. anthracis strains in France as can be revealed by MLVA. The data obtained suggests that once such coverage is achieved, it becomes possible to devise optimized first-line MLVA assays comprising a sufficiently low number of loci to be typed either in one multiplex PCR or on agarose gels. Such a selection of seven loci is proposed here, and future similar investigations in additional countries will indicate to which extend the same selection can be used worldwide as a common minimum set. It is hoped that this approach will contribute to an efficient and low-cost routine surveillance of important pathogens for biosecurity such as B. anthracis.

High-throughput sequencing of Bacillus anthracis in France: investigating genome diversity and population structure using whole-genome SNP discovery

Background

Single nucleotide polymorphisms (SNPs) are ideal signatures for subtyping monomorphic pathogens such as Bacillus anthracis. Here we report the use of next-generation sequencing technology to investigate the historical, geographic and genetic diversity of Bacillus anthracis in France. 122 strains isolated over a 60-years period throughout the country were whole-genome sequenced and comparative analyses were carried out with a focus on SNPs discovery to discriminate regional sub-groups of strains.

Results

A total of 1581 chromosomal SNPs precisely establish the phylogenetic relationships existing between the French strains. Phylogeography patterns within the three canSNP sub-lineages present in France (i.e. B.Br.CNEVA, A.Br.011/009 and A.Br.001/002) were observed. One of the more remarkable findings was the identification of a variety of genotypes within the A.Br.011/009 sub-group that are persisting in the different regions of France. The 560 SNPs defining the A.Br.011/009- affiliated French strains split the Trans-Eurasian sub-group into six distinct branches without any intermediate nodes. Distinct sub-branches, with some geographic clustering, were resolved. The 345 SNPs defining the major B.Br CNEVA sub-lineage clustered three main phylogeographic clades, the Alps, the Pyrenees, and the Massif Central, with a small Saône-et-Loire sub-cluster nested within the latter group. The French strains affiliated to the minor A.Br.001/002 group were characterized by 226 SNPs. All recent isolates collected from the Doubs department were closely related. Identification of SNPs from whole-genome sequences facilitates high-resolution strain tracking and provides the level of discrimination required for outbreak investigations. Eight diagnostic SNPs, representative of the main French-specific phylogeographic clusters, were therefore selected and developed into high-resolution melting SNP discriminative assays.

Conclusions

This work has established one of the most accurate phylogenetic reconstruction of B. anthracis population structure in a country. An extensive next-generation sequencing (NGS) dataset of 122 French strains have been created that allowed the identification of novel diagnostic SNPs useful to rapidly determine the geographic origin of any strain found in France.

When outgroups fail; phylogenomics of rooting the emerging pathogen, Coxiella burnetii

Rooting phylogenies is critical for understanding evolution, yet the importance, intricacies and difficulties of rooting are often overlooked. For rooting, polymorphic characters among the group of interest (ingroup) must be compared to those of a relative (outgroup) that diverged before the last common ancestor (LCA) of the ingroup. Problems arise if an outgroup does not exist, is unknown, or is so distant that few characters are shared, in which case duplicated genes originating before the LCA can be used as proxy outgroups to root diverse phylogenies. Here, we describe a genome-wide expansion of this technique that can be used to solve problems at the other end of the evolutionary scale: where ingroup individuals are all very closely related to each other, but the next closest relative is very distant. We used shared orthologous single nucleotide polymorphisms (SNPs) from 10 whole genome sequences of Coxiella burnetii, the causative agent of Q fever in humans, to create a robust, but unrooted phylogeny. To maximize the number of characters informative about the rooting, we searched entire genomes for polymorphic duplicated regions where orthologs of each paralog could be identified so that the paralogs could be used to root the tree. Recent radiations, such as those of emerging pathogens, often pose rooting challenges due to a lack of ingroup variation and large genomic differences with known outgroups. Using a phylogenomic approach, we created a robust, rooted phylogeny for C. burnetii. [Coxiella burnetii; paralog SNPs; pathogen evolution; phylogeny; recent radiation; root; rooting using duplicated genes.].

Ecological niche modeling of Bacillus anthracis on three continents: evidence for genetic-ecological divergence?

We modeled the ecological niche of a globally successful Bacillus anthracis sublineage in the United States, Italy and Kazakhstan to better understand the geographic distribution of anthrax and potential associations between regional populations and ecology. Country-specific ecological-niche models were developed and reciprocally transferred to the other countries to determine if pathogen presence could be accurately predicted on novel landscapes. Native models accurately predicted endemic areas within each country, but transferred models failed to predict known occurrences in the outside countries. While the effects of variable selection and limitations of the genetic data should be considered, results suggest differing ecological associations for the B. anthracis populations within each country and may reflect niche specialization within the sublineage. Our findings provide guidance for developing accurate ecological niche models for this pathogen; models should be developed regionally, on the native landscape, and with consideration to population genetics. Further genomic analysis will improve our understanding of the genetic-ecological dynamics of B. anthracis across these countries and may lead to more refined predictive models for surveillance and proactive vaccination programs. Further studies should evaluate the impact of variable selection of native and transferred models.

Draft Genome Sequences of Two Bulgarian Bacillus anthracis Strains

Bacillus anthracis strains previously isolated from Bulgaria form a unique subcluster within the A1.a cluster that is typical for isolates from southeastern Europe. Here, we report the draft genome sequences of two Bulgarian B. anthracis strains belonging to the A branch (A.Br.)008/009 canonical single nucleotide polymorphism (SNP) group of the major A branch.

Grappling with Proteus: population level approaches to understanding microbial diversity

The emerging fields of microbial population genetics and genomics provide an avenue to study the ecological rules that govern how communities form, function, and evolve. Our struggle to understand the causes and consequences of microbial diversity stems from our inability to define ecologically and evolutionarily meaningful units of diversity. The 16S rRNA-based tools that have been so useful in charting microbial diversity may lack sufficient sensitivity to answer many questions about the ecology and evolution of microbes. Examining genetic diversity with increased resolution is vital to understanding the forces shaping community structure. Population genetic analyses enabled by whole genome sequencing, multilocus sequence analyses, or single-nucleotide polymorphism analyses permit the testing of hypotheses pertaining to the geographic distribution, migration, and habitat preference of specific microbial lineages. Furthermore, these approaches can reveal patterns of gene exchange within and between populations and communities. Tools from microbial population genetics and population genomics can be used to increase the resolution with which we measure microbial diversity, enabling a focus on the scale of genetic diversity at which ecological processes impact evolutionary events. This tighter focus promises to improve our understanding of the causes and consequences of microbial community structure.

Molecular epidemiologic investigation of an anthrax outbreak among heroin users, Europe

In December 2009, two unusual cases of anthrax were diagnosed in heroin users in Scotland. A subsequent anthrax outbreak in heroin users emerged throughout Scotland and expanded into England and Germany, sparking concern of nefarious introduction of anthrax spores into the heroin supply. To better understand the outbreak origin, we used established genetic signatures that provided insights about strain origin. Next, we sequenced the whole genome of a representative Bacillus anthracis strain from a heroin user (Ba4599), developed Ba4599-specific single-nucleotide polymorphism assays, and genotyped all available material from other heroin users with anthrax. Of 34 case-patients with B. anthracis-positive PCR results, all shared the Ba4599 single-nucleotide polymorphism genotype. Phylogeographic analysis demonstrated that Ba4599 was closely related to strains from Turkey and not to previously identified isolates from Scotland or Afghanistan, the presumed origin of the heroin. Our results suggest accidental contamination along the drug trafficking route through a cutting agent or animal hides used to smuggle heroin into Europe.

Genomic characterization of the Bacillus cereus sensu lato species: backdrop to the evolution of Bacillus anthracis

The key genes required for Bacillus anthracis to cause anthrax have been acquired recently by horizontal gene transfer. To understand the genetic background for the evolution of B. anthracis virulence, we obtained high-redundancy genome sequences of 45 strains of the Bacillus cereus sensu lato (s.l.) species that were chosen for their genetic diversity within the species based on the existing multilocus sequence typing scheme. From the resulting data, we called more than 324,000 new genes representing more than 12,333 new gene families for this group. The core genome size for the B. cereus s.l. group was ∼1750 genes, with another 2150 genes found in almost every genome constituting the extended core. There was a paucity of genes specific and conserved in any clade. We found no evidence of recent large-scale gene loss in B. anthracis or for unusual accumulation of nonsynonymous DNA substitutions in the chromosome; however, several B. cereus genomes isolated from soil and not previously associated with human disease were degraded to various degrees. Although B. anthracis has undergone an ecological shift within the species, its chromosome does not appear to be exceptional on a macroscopic scale compared with close relatives.

Characterization of genetic diversity of Bacillus anthracis in France by using high-resolution melting assays and multilocus variable-number tandem-repeat analysis

Using high-resolution melting (HRM) analysis, we developed a cost-effective method to genotype a set of 13 phylogenetically informative single-nucleotide polymorphisms (SNPs) within the genome of Bacillus anthracis. SNP discrimination assays were performed in monoplex or duplex and applied to 100 B. anthracis isolates collected in France from 1953 to 2009 and a few reference strains. HRM provided a reliable and cheap alternative to subtype B. anthracis into one of the 12 major sublineages or subgroups. All strains could be correctly positioned on the canonical SNP (canSNP) phylogenetic tree, except the divergent Pasteur vaccine strain ATCC 4229. We detected the cooccurrence of three canSNP subgroups in France. The dominant B.Br.CNEVA sublineage was found to be prevalent in the Alps, the Pyrenees, the Auvergne region, and the Saône-et-Loire department. Strains affiliated with the A.Br.008/009 subgroup were observed throughout most of the country. The minor A.Br.001/002 subgroup was restricted to northeastern France. Multiple-locus variable-number tandem-repeat analysis using 24 markers further resolved French strains into 60 unique profiles and identified some regional patterns. Diversity found within the A.Br.008/009 and B.Br.CNEVA subgroups suggests that these represent old, ecologically established clades in France. Phylogenetic relationships with strains from other parts of the world are discussed.

New insights into the biological effects of anthrax toxins: linking cellular to organismal responses

The anthrax toxins lethal toxin (LT) and edema toxin (ET) are essential virulence factors produced by Bacillus anthracis. These toxins act during two distinct phases of anthrax infection. During the first, prodromal phase, which is often asymptomatic, anthrax toxins act on cells of the immune system to help the pathogen establish infection. Then, during the rapidly progressing (or fulminant) stage of the disease bacteria disseminate via a hematological route to various target tissues and organs, which are typically highly vascularized. As bacteria proliferate in the bloodstream, LT and ET begin to accumulate rapidly reaching a critical threshold level that will cause death even when the bacterial proliferation is curtailed by antibiotics. During this final phase of infection the toxins cause an increase in vascular permeability and a decrease in function of target organs including the heart, spleen, kidney, adrenal gland, and brain. In this review, we examine the various biological effects of anthrax toxins, focusing on the fulminant stage of the disease and on mechanisms by which the two toxins may collaborate to cause cardiovascular collapse. We discuss normal mechanisms involved in maintaining vascular integrity and based on recent studies indicating that LT and ET cooperatively inhibit membrane trafficking to cell-cell junctions we explore several potential mechanisms by which the toxins may achieve their lethal effects. We also summarize the effects of other potential virulence factors secreted by B. anthracis and consider the role of toxic factors in the evolutionarily recent emergence of this devastating disease.

Bovine Bacillus anthracis in Cameroon

Bovine Bacillus anthracis isolates from Cameroon were genetically characterized. They showed a strong homogeneity, and they belong, together with strains from Chad, to cluster Aβ, which appears to be predominant in western Africa. However, one strain that belongs to a newly defined clade (D) and cluster (D1) is penicillin resistant and shows certain phenotypes typical of Bacillus cereus.

Cost-effective interrogation of single nucleotide polymorphisms using the mismatch amplification mutation assay and capillary electrophoresis

The ability to characterize SNPs is an important aspect of many clinical diagnostic, genetic and evolutionary studies. Here, we designed a multiplexed SNP genotyping method to survey a large number of phylogenetically informative SNPs within the genome of the bacterium Bacillus anthracis. This novel method, CE universal tail mismatch amplification mutation assay (CUMA), allows for PCR multiplexing and automatic scoring of SNP genotypes, thus providing a rapid, economical and higher throughput alternative to more expensive SNP genotyping techniques. CUMA delivered accurate B. anthracis SNP genotyping results and, when multiplexed, saved reagent costs by more than 80% compared with TaqMan real-time PCR. When real-time PCR technology and instrumentation is unavailable or the reagents are cost-prohibitive, CUMA is a powerful alternative for SNP genotyping.

MLVA and SNP analysis identified a unique genetic cluster in Bulgarian Bacillus anthracis strains

A collection of 40 Bacillus anthracis strains mostly isolated from soil in Bulgaria between 1960 and 1980 were investigated. All strains were proven to be B. anthracis by culture and amplification of a B. anthracis-specific chromosomal marker. PCR demonstrated that in nine strains both virulence plasmids (pX01+/pX02+) and in four strains only one plasmid (pX02+) were present, whereas the majority of strains (n = 27) lacked both plasmids (pX01-/pX02-). Multi-locus-variable number of tandem repeat-analysis (MLVA) using 15 markers differentiated three genotypes. Comparison with typing data of more than 1,000 different B. anthracis strains revealed that Bulgarian genotypes affiliated with the A1.a cluster and form their own unique cluster different from clusters containing strains isolated in geographical proximity, e.g., Turkey, Georgia, Hungary, Albania or Italy. In addition, a new allele of one marker (vrrC2) was identified. Canonical single nucleotide polymorphisms analysis allocated 31 Bulgarian strains into the A.Br.008/009 and nine strains into the A.Br.WNA group, which is the first description of B. anthracis strains of the A.Br.WNA group on the Eurasian continent.

Anthrax outbreak in a Swedish beef cattle herd--1st case in 27 years: Case report

After 27 years with no detected cases, an outbreak of anthrax occurred in a beef cattle herd in the south of Sweden. The outbreak was unusual as it occurred in winter, in animals not exposed to meat-and-bone meal, in a non-endemic country. The affected herd consisted of 90 animals, including calves and young stock. The animals were kept in a barn on deep straw bedding and fed only roughage. Seven animals died during 10 days, with no typical previous clinical signs except fever. The carcasses were reportedly normal in appearance, particularly as regards rigor mortis, bleeding and coagulation of the blood. Subsequently, three more animals died and anthrax was suspected at necropsy and confirmed by culture and PCR on blood samples. The isolated strain was susceptible to tetracycline, ciprofloxacin and ampicillin. Subtyping by MLVA showed the strain to cluster with isolates in the A lineage of Bacillus anthracis. Environmental samples from the holding were all negative except for two soil samples taken from a spot where infected carcasses had been kept until they were picked up for transport. The most likely source of the infection was concluded to be contaminated roughage, although this could not be substantiated by laboratory analysis. The suspected feed was mixed with soil and dust and originated from fields where flooding occurred the previous year, followed by a dry summer with a very low water level in the river allowing for the harvesting on soil usually not exposed. In the early 1900s, animal carcasses are said to have been dumped in this river during anthrax outbreaks and it is most likely that some anthrax spores could remain in the area. The case indicates that untypical cases in non-endemic areas may be missed to a larger extent than previously thought. Field tests allowing a preliminary risk assessment of animal carcasses would be helpful for increased sensitivity of detection and prevention of further exposure to the causative agent.

Humans and evolutionary and ecological forces shaped the phylogeography of recently emerged diseases

The development of human civilizations and global commerce has led to the emergence and worldwide circulation of many infectious diseases. Anthrax, plague and tularaemia are three zoonotic diseases that have been intensely studied through genome characterization of the causative species and phylogeographical analyses. A few highly fit genotypes in each species represent the causative agents for most of the observed disease cases. Together, mutational and selective forces create highly adapted pathogens, but this must be coupled with ecological opportunities for global expansion. This Review describes the distributions of the bacteria that cause anthrax, plague and tularaemia and investigates the forces that created clonal structures in these species.

Bacillus anthracis in China and its relationship to worldwide lineages

Background

The global pattern of distribution of 1033 B. anthracis isolates has previously been defined by a set of 12 conserved canonical single nucleotide polymorphisms (canSNP). These studies reinforced the presence of three major lineages and 12 sub-lineages and sub-groups of this anthrax-causing pathogen. Isolates that form the A lineage (unlike the B and C lineages) have become widely dispersed throughout the world and form the basis for the geographical disposition of "modern" anthrax. An archival collection of 191 different B. anthracis isolates from China provides a glimpse into the possible role of Chinese trade and commerce in the spread of certain sub-lineages of this pathogen. Canonical single nucleotide polymorphism (canSNP) and multiple locus VNTR analysis (MLVA) typing has been used to examine this archival collection of isolates.

Results

The canSNP study indicates that there are 5 different sub-lineages/sub-groups in China out of 12 previously described world-wide canSNP genotypes. Three of these canSNP genotypes were only found in the western-most province of China, Xinjiang. These genotypes were A.Br.008/009, a sub-group that is spread across most of Europe and Asia; A.Br.Aust 94, a sub-lineage that is present in Europe and India, and A.Br.Vollum, a lineage that is also present in Europe. The remaining two canSNP genotypes are spread across the whole of China and belong to sub-group A.Br.001/002 and the A.Br.Ames sub-lineage, two closely related genotypes. MLVA typing adds resolution to the isolates in each canSNP genotype and diversity indices for the A.Br.008/009 and A.Br.001/002 sub-groups suggest that these represent older and established clades in China.

Conclusion

B. anthracis isolates were recovered from three canSNP sub-groups (A.Br.008/009, A.Br.Aust94, and A.Br.Vollum) in the western most portion of the large Chinese province of Xinjiang. The city of Kashi in this province appears to have served as a crossroads for not only trade but the movement of diseases such as anthrax along the ancient "silk road". Phylogenetic inference also suggests that the A.Br.Ames sub-lineage, first identified in the original Ames strain isolated from Jim Hogg County, TX, is descended from the A.Br.001/002 sub-group that has a major presence in most of China. These results suggest a genetic discontinuity between the younger Ames sub-lineage in Texas and the large Western North American sub-lineage spread across central Canada and the Dakotas.