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

New Research on the Bacillus anthracis Genetic Diversity in Siberia

Anthrax is a particularly dangerous infection of humans and ungulates caused by the Gram-positive spore-forming bacterium Bacillus anthracis. The highly monomorphic and clonal species B. anthracis is commonly divided into three main lineages, A, B, and C, which in turn are divided into several canSNP groups. We report here a phylogenetic analysis based on the whole-genome sequence (WGS) data of fifteen strains isolated predominantly in Siberia or Central and Southern Russia. We confirm the wide distribution of the cluster of strains of the B.Br.001/002 group, endemic to the Russian Arctic, which is also present in the steppe zone of Southern Siberia. We characterize additional branches within the major A.Br.001/002 polytomy comprising the A.Br.Ames and A.Br.Sterne lineages, one of which is identified in the Arctic.

Ancient Rapid Radiation Explains Most Conflicts Among Gene Trees and Well-Supported Phylogenomic Trees of Nostocalean Cyanobacteria

Prokaryotic genomes are often considered to be mosaics of genes that do not necessarily share the same evolutionary history due to widespread horizontal gene transfers (HGTs). Consequently, representing evolutionary relationships of prokaryotes as bifurcating trees has long been controversial. However, studies reporting conflicts among gene trees derived from phylogenomic data sets have shown that these conflicts can be the result of artifacts or evolutionary processes other than HGT, such as incomplete lineage sorting, low phylogenetic signal, and systematic errors due to substitution model misspecification. Here, we present the results of an extensive exploration of phylogenetic conflicts in the cyanobacterial order Nostocales, for which previous studies have inferred strongly supported conflicting relationships when using different concatenated phylogenomic data sets. We found that most of these conflicts are concentrated in deep clusters of short internodes of the Nostocales phylogeny, where the great majority of individual genes have low resolving power. We then inferred phylogenetic networks to detect HGT events while also accounting for incomplete lineage sorting. Our results indicate that most conflicts among gene trees are likely due to incomplete lineage sorting linked to an ancient rapid radiation, rather than to HGTs. Moreover, the short internodes of this radiation fit the expectations of the anomaly zone, i.e., a region of the tree parameter space where a species tree is discordant with its most likely gene tree. We demonstrated that concatenation of different sets of loci can recover up to 17 distinct and well-supported relationships within the putative anomaly zone of Nostocales, corresponding to the observed conflicts among well-supported trees based on concatenated data sets from previous studies. Our findings highlight the important role of rapid radiations as a potential cause of strongly conflicting phylogenetic relationships when using phylogenomic data sets of bacteria. We propose that polytomies may be the most appropriate phylogenetic representation of these rapid radiations that are part of anomaly zones, especially when all possible genomic markers have been considered to infer these phylogenies. [Anomaly zone; bacteria; horizontal gene transfer; incomplete lineage sorting; Nostocales; phylogenomic conflict; rapid radiation; Rhizonema.].

Spanish Outbreak Isolates Bridge Phylogenies of European and American Bacillus anthracis

The geographical origin of a major present-day phylogenetic group (A branch WNA; A.Br.WNA) of American Bacillus anthracis is controversial. One hypothesis postulated that the anthrax pathogen reached North America via a then-existing land bridge from northeastern Asia thousands of years ago. A competing hypothesis suggested that B. anthracis was introduced to America a couple of hundred years ago, related to European colonization. The latter view is strongly supported by genomic analysis of a group of French B. anthracis isolates that are phylogenetically closely related to the North American strains of the A branch A.Br.WNA clade. In addition, three West African strains also belong to this relationship group. Recently, we have added a Spanish strain to these close relatives of the WNA lineage of American B. anthracis. Nevertheless, the diversity of Spanish B. anthracis remains largely unexplored, and phylogenetic links to European or American relatives are not well resolved. Here, we genome sequenced and characterized 29 new B. anthracis isolates (yielding 18 unique genotypes) from outbreaks in west central and central Spain in 2021. Applying comparative chromosomal analysis, we placed the chromosomes of these isolates within the established phylogeny of the A.Br.008/009 (A.Br.TEA) canonical SNP group. From this analysis, a new sub-clade, named A.Br.11/ESPc, emerged that constitutes a sister group of American A.Br.WNA.

Some Peculiarities of Anthrax Epidemiology in Herbivorous and Carnivorous Animals

Anthrax is an especially dangerous zooanthroponosis caused by the Gram-positive spore-forming bacterium Bacillus anthracis. A notable feature of this disease is the difference in susceptibility to it among different groups of animals. Anthrax primarily affects herbivorous ungulate mammals; they are easily infected, and their disease often leads to rapid, even sudden, death. However, predators and scavengers are extremely resistant to anthrax, and if they become infected, they usually become mildly ill. As the result of the increased sensitivity of ungulates to anthrax and the possibility of disease transmission from them to humans, most studies of anthrax have focused on the diagnosis, prevention, and treatment of infection in farm animals and humans. The issues of anthrax in other animals, such as predators, and the peculiarities of anthrax epidemiology in wild ungulates have not been sufficiently detailed in the literature. In this article, we provide a review of literature sources that describe the differential susceptibility to infection of various groups of animals to anthrax and some epidemiological features of anthrax in animals that are not the main hosts of B. anthracis.

High-Resolution Melting PCR as Rapid Genotyping Tool for Brucella Species

Brucella sp. are the causative agents of brucellosis. One of the main characteristics of the Brucella genus concerns its very high genetic homogeneity. To date, classical bacteriology typing is still considered as the gold standard assay for direct diagnosis of Brucella. Molecular approaches are routinely used for the identification of Brucella at the genus level. However, genotyping is more complex, and to date, no method exists to quickly assign a strain into species and biovar levels, and new approaches are required. Next generation sequencing (NGS) opened a new era into the diagnosis of bacterial diseases. In this study, we designed a high-resolution melting (HRM) method for the rapid screening of DNA and direct assignment into one of the 12 species of the Brucella genus. This method is based on 17 relevant single nucleotide polymorphisms (SNPs), identified and selected from a whole genome SNP (wgSNP) analysis based on 988 genomes (complete and drafts). These markers were tested against the collection of the European Reference Laboratory (EU-RL) for brucellosis (1440 DNAs extracted from Brucella strains). The results confirmed the reliability of the panel of 17 SNP markers, allowing the differentiation of each species of Brucella together with biovars 1, 2, and 3 of B. suis and vaccine strain Rev1 (B. melitensis) within 3 h, which is a considerable gain of time for brucellosis diagnosis. Therefore, this genotyping tool provides a new and quick alternative for Brucella identification based on SNPs with the HRM-PCR assay.

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.

Sequence Variability of pXO1-Located Pathogenicity Genes of Bacillus anthracis Natural Strains of Different Geographic Origin

The main pathogenic factor of Bacillus anthracis is a three-component toxin encoded by the pagA, lef, and cya genes, which are located on the pXO1 plasmid. The atxA gene, which encodes the primary regulator of pathogenicity factor expression, is located on the same plasmid. In this work, we evaluated the polymorphism of the pagA, lef, cya, and atxA genes for 85 B. anthracis strains from different evolutionary lineages and canSNP groups. We have found a strong correlation of 19 genotypes with the main evolutionary lineages, but the correlation with the canSNP group of the strain was not as strong. We have detected several genetic markers indicating the geographical origin of the strains, for example, their source from the steppe zone of the former USSR. We also found that strains of the B.Br.001/002 group caused an anthrax epidemic in Russia in 2016 and strains isolated during paleontological excavations in the Russian Arctic have the same genotype as the strains of the B.Br.CNEVA group circulating in Central Europe. This data could testify in favor of the genetic relationship of these two groups of strains and hypothesize the ways of distribution of their ancestral forms between Europe and the Arctic.

A Whole-Genome-Based Gene-by-Gene Typing System for Standardized High-Resolution Strain Typing of Bacillus anthracis

Whole-genome sequencing (WGS) has been established for bacterial subtyping and is regularly used to study pathogen transmission, to investigate outbreaks, and to perform routine surveillance. Core-genome multilocus sequence typing (cgMLST) is a bacterial subtyping method that uses WGS data to provide a high-resolution strain characterization. This study aimed at developing a novel cgMLST scheme for Bacillus anthracis, a notorious pathogen that causes anthrax in livestock and humans worldwide. The scheme comprises 3,803 genes that were conserved in 57 B. anthracis genomes spanning the whole phylogeny. The scheme has been evaluated and applied to 584 genomes from 50 countries. On average, 99.5% of the cgMLST targets were detected. The cgMLST results confirmed the classical canonical single-nucleotide-polymorphism (SNP) grouping of B. anthracis into major clades and subclades. Genetic distances calculated based on cgMLST were comparable to distances from whole-genome-based SNP analysis with similar phylogenetic topology and comparable discriminatory power. Additionally, the application of the cgMLST scheme to anthrax outbreaks from Germany and Italy led to a definition of a cutoff threshold of five allele differences to trace epidemiologically linked strains for cluster typing and transmission analysis. Finally, the association of two clusters of B. anthracis with human cases of injectional anthrax in four European countries was confirmed using cgMLST. In summary, this study presents a novel cgMLST scheme that provides high-resolution strain genotyping for B. anthracis. This scheme can be used in parallel with SNP typing methods to facilitate rapid and harmonized interlaboratory comparisons, essential for global surveillance and outbreak analysis. The scheme is publicly available for application by users, including those with little bioinformatics knowledge.

Phylogenetics of Bacillus anthracis isolates from Russia and bordering countries

Anthrax is a concern for public health and veterinary medicine in Russia. The available phylogenetic data on isolates from Russia and neighboring CIS countries are clearly not enough to gain a better understanding of their position in the global phylogenetic population structure of this pathogen. In this study, we analyzed the genomes of 66 Bacillus anthracis strains, which were isolated between 1935 and 2019 from different sources in Russia, as well as in Ukraine, Azerbaijan, Georgia, Armenia and Moldova. Whole genome SNP analysis of genomes of 66 strains obtained in this study along with 222 B. anthracis genomes available in the GenBank database revealed 7242 SNPs used to construct a phylogenetic reconstruction with the method of Maximum Likelihood. Studied strains belong to 6 different genetic groups: A.Br.008(A.Br.008/009), A.Br.081(Ames), A.Br.014(A.Br.Aust94), A.Br.082(A.Br.001/002), A.Br.034(A.Br.005/006, Ancient A) and B.Br.002 (B.Br.001/002). Within the group A.Br.014(A.Br.Aust94) a subcluster A.Br.029 of strains isolated in Georgia, Armenia, Azerbaijan, Russia (Republic of Dagestan) and Turkey, named Caucasus-East Anatolia (CEA), was identified. In the subgroup A.Br.105(Tsiankovskii) the cluster A.Br.117 of strains from Russia, Ukraine and Slovakia are assigned, in the subgroup A.Br118 (STI) - cluster A.Br.123 with strains from Russia and Georgia and cluster A.Br.125 with strains from Republic of Dagestan. New subclusters B.Br.017("EUROPE") were identified in the B.Br.002(B.Br.001/002) cluster, represented by strains from the European part of Russia, as well as from South Korea and Finland. For 8 clusters and subclusters, the SNP markers were identified. The study confirmed a significant genetic diversity of the strains isolated in Russia and border countries and clarified their position in the phylogenetic structure of the global B. anthracis population. New genetic clusters A.Br.029 (CEA), A.Br.117, A.Br.123, A.Br.125, and B.Br.017 («EUROPE») were defined. 96 marker SNPs specific for these clusters were identified.

Proposal of a Taxonomic Nomenclature for the Bacillus cereus Group Which Reconciles Genomic Definitions of Bacterial Species with Clinical and Industrial Phenotypes

The Bacillus cereus group comprises numerous closely related species, including bioterrorism agent B. anthracis, foodborne pathogen B. cereus, and biopesticide B. thuringiensis Differentiating organisms capable of causing illness or death from those used in industry is essential for risk assessment and outbreak preparedness. However, current species definitions facilitate species-phenotype incongruences, particularly when horizontally acquired genes are responsible for a phenotype. Using all publicly available B. cereus group genomes (n = 2,231), we show that current species definitions lead to overlapping genomospecies clusters, in which 66.2% of genomes belong to multiple genomospecies at a conventional 95 average nucleotide identity (ANI) genomospecies threshold. A genomospecies threshold of ≈92.5 ANI is shown to reflect a natural gap in genome similarity for the B. cereus group, and medoid genomes identified at this threshold are shown to yield resolvable genomospecies clusters with minimal overlap (six of 2,231 genomes assigned to multiple genomospecies; 0.269%). We thus propose a nomenclatural framework for the B. cereus group which accounts for (i) genomospecies using resolvable genomospecies clusters obtained at ≈92.5 ANI, (ii) established lineages of medical importance using a formal collection of subspecies names, and (iii) heterogeneity of clinically and industrially important phenotypes using a formalized and extended collection of biovar terms. We anticipate that the proposed nomenclature will remain interpretable to clinicians, without sacrificing genomic species definitions, which can in turn aid in pathogen surveillance; early detection of emerging, high-risk genotypes; and outbreak preparedness.IMPORTANCE Historical species definitions for many prokaryotes, including pathogens, have relied on phenotypic characteristics that are inconsistent with genome evolution. This scenario forces microbiologists and clinicians to face a tradeoff between taxonomic rigor and clinical interpretability. Using the Bacillus cereus group as a model, a conceptual framework for the taxonomic delineation of prokaryotes which reconciles genomic definitions of species with clinically and industrially relevant phenotypes is presented. The nomenclatural framework outlined here serves as a model for genomics-based bacterial taxonomy that moves beyond arbitrarily set genomospecies thresholds while maintaining congruence with phenotypes and historically important species names.

A Rare Glimpse into the Past of the Anthrax Pathogen Bacillus anthracis

The bacterium Bacillus anthracis is the causative agent of the zoonotic disease anthrax. While genomics of extant B. anthracis isolates established in-depth phylogenomic relationships, there is scarce information on the historic genomics of the pathogen. Here, we characterized the oldest documented B. anthracis specimen. The inactive 142-year-old material originated from a bovine diseased in Chemnitz (Germany) in 1878 and is contemporary with the seminal studies of Robert Koch on B. anthracis. A specifically developed isolation method yielded high-quality DNA from this specimen for genomic sequencing. The bacterial chromosome featuring 242 unique base-characters placed it into a major phylogenetic clade of B. anthracis (B.Branch CNEVA), which is typical for central Europe today. Our results support the notion that the CNEVA-clade represents part of the indigenous genetic lineage of B. anthracis in this part of Europe. This work emphasizes the value of historic specimens as precious resources for reconstructing the past phylogeny of the anthrax pathogen.

Epidemiological and comparative genomic analysis of Bacillus anthracis isolated from northern Vietnam

To understand the epidemiological and genetic background of anthrax cases occurring in Vietnam from 2011 to 2015, we surveilled and genetically analyzed Bacillus anthracis isolated in the north of the country. Epidemiological surveillance showed that most human cutaneous anthrax cases occurred in association with animal dissection. Whole-genome sequences were obtained from six B. anthracis strains from human patients with cutaneous anthrax in the endemic area. Comparative genomic analysis showed that the genetic homogeneity among Vietnamese B. anthracis strains was very high. All Vietnamese B. anthracis strains belonged to the canSNP lineage of A.Br.011/009, which mostly consists of strains of the trans-Eurasian (TEA) group, including the most closely related strain, Carbosap. To clarify the genetic diversity of Vietnamese strains and strains belonging to A.Br.011/009 and A.Br.008/011 canSNP lineages, we applied a reference genome-based single-nucleotide polymorphism (SNP) and gene-by-gene genomic analysis (whole-genome MLST) strategy. The phylogeny from core genome SNPs revealed that the Vietnamese strains were positioned close to each other; moreover, several SNPs specific to Vietnamese B. anthracis were identified. Whole-genome MLST analysis revealed the differences in the number of SNPs between Vietnamese strains, which could enable discrimination at the strain level.

Phylogenomic structure of Bacillus anthracis isolates in the Northern Cape Province, South Africa revealed novel single nucleotide polymorphisms

Bacillus anthracis, the aetiological agent of anthrax, is regarded as a highly monomorphic pathogen that presents a low genetic diversity using standard molecular techniques. Whole genome sequencing and single nucleotide polymorphisms (SNPs) are definitive signatures for subtyping of B. anthracis. Here we employed whole genome single nucleotide polymorphism (wgSNP) analysis to investigate the genetic diversity of B. anthracis in the historically endemic region of Northern Cape Province (NCP), South Africa. Twenty-six isolates from anthrax outbreaks that occurred between 1998 and 2008/9 in NCP as well as from Namibia-South Africa Transfontier Conservation area and Botswana were compared to global B. anthracis genomes. Most NCP B. anthracis strains (n = 22) clustered in the A.Br.003/004 (A.Br.101) branch and are closely related to the Zimbabwe and Mozambique strains (A.Br.102 branch). A total of 4923 parsimony informative-SNPs accurately established the A.Br.003/004 phylogenetic relationships of the NCP isolates into two distinct sub-clades and SNP markers designated as A.Br.172 and A.Br.173 were developed. Other NCP strains (n = 2) grouped in the A.Br.001/002 (Sterne) branch while strains (n = 2) from the Namibia-South Africa Transfontier Conservation area and Botswana clustered in A.Br.005/006 (Ancient A) branch. The sequenced B. anthracis strains (A0094, A0096 and A0097) that clustered in the A.Br.064 (V770) clade were isolated from Vaalbos National Park and similar strains have not been isolated. The B. anthracis A0088 strain cluster with the NCP strains in the A.Br.003/004 (A.Br.172) SNP branch which has been isolated in NCP, South Africa. This study highlights the phylogenetic structure of NCP B. anthracis strains with distinctive SNP branches important for forensic tracing and novel SNP discovery purposes. The sequenced strains will serve as a means to further trace the dissemination of B. anthracis outbreaks in NCP, South Africa, and on the continent, as well as for forensic tracking on a global scale.

Phylogenetic Placement of Isolates Within the Trans-Eurasian Clade A.Br.008/009 of Bacillus anthracis

The largest phylogenetic lineage known to date of the anthrax pathogen Bacillus anthracis is the wide-spread, so-called Trans-Eurasian clade systematically categorized as the A.Br.008/009 group sharing two defining canonical single-nucleotide polymorphisms (canSNP). In this study, we genome-sequenced a collection of 35 B. anthracis strains of this clade, derived from human infections, animal outbreaks or soil, mostly from European countries isolated between 1936 and 2008. The new data were subjected to comparative chromosomal analysis, together with 75 B. anthracis genomes available in public databases, and the relative placements of these isolates were determined within the global phylogeny of the A.Br.008/009 canSNP group. From this analysis, we have detected 3754 chromosomal SNPs, allowing the assignation of the new chromosomal sequences to established sub-clades, to define new sub-clades, such as two new Spanish, one Bulgarian or one German group(s), or to introduce orphan lineages. SNP-based results were compared with that of a multilocus variable number of tandem repeat analysis (MLVA). This analysis indicated that MLVA typing might provide additional information in cases when genomics yields identical genotypes or shows only minor differences. Introducing the delayed mismatch amplification assay (DMAA) PCR-analysis, we developed a cost-effective method to interrogate for a set of ten phylogenetically informative SNPs within genomes of A.Br.008/009 canSNP clade strains of B. anthracis. By this approach, additional 32 strains could be assigned to five of ten defined clades.

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.

Comparative genomics and phylogenetic analysis of Bacillus anthracis strains isolated from domestic animals in Japan

Anthrax, caused by Bacillus anthracis, is a severe zoonosis with a great impact on both human and animal health. In the present study, we identified the phylogenetic relationships among 16 Japanese strains of B. anthracis, including eight bovine strains, two equine strains, five swine strains, and one former vaccine strain, using in silico canonical single nucleotide polymorphism (canSNP) and core genome SNP analyses. The results of our in silico canSNP analysis suggest that these 16 Japanese strains could be divided into four lineages: i) one equine strain in A.Br.Ames, ii) one equine and six bovine strains in A.Br.001/002, iii) five swine and one bovine strain in A.Br.Aust94, and iv) one bovine and one vaccine strain in A.Br.008/011. A comparison with non-Japanese B. anthracis strains revealed a total of 3787 SNPs identified from the whole genome sequences of the Japanese strains; these SNP data were subjected to a phylogenetic analysis using the maximum parsimony (MP) method. Our core genome SNP analysis was also able to detect differences of a few chromosomal SNPs across clonal strains from the same cases that had different storage and passage histories. Additionally, our whole genome SNP analysis clearly indicated that the Japanese swine anthrax cases of 1982 were caused by at least three independent strains; however, their phylogeny revealed no clear relationship with swine strains from other countries. The bovine strain belonging to the A.Br.008/011 lineage differed from a former Japanese vaccine strain by only 12 SNPs. Together with the phylogenic results and epidemiological circumstances, the diversity of strains reveals that the B. anthracis available in Japan probably resulted from multiple relatively recent import events, rather than reflecting the persistence of a more ancient ecologically established group.

Loss of Bacitracin Resistance Due to a Large Genomic Deletion among Bacillus anthracis Strains

Anthrax is caused by Bacillus anthracis, an endospore-forming soil bacterium. The genetic diversity of B. anthracis is known to be low compared with that of Bacillus species. In this study, we performed whole-genome sequencing of Zambian isolates of B. anthracis to understand the genetic diversity between closely related strains. Comparison of genomic sequences revealed that closely related strains were separated into three groups based on single nucleotide polymorphisms distributed throughout the genome. A large genomic deletion was detected in the region containing a bacitracin resistance gene cluster flanked by rRNA operons, resulting in the loss of bacitracin resistance. The structure of the deleted region, which was also conserved among species of the Bacillus cereus group, has the potential for both deletion and amplification and thus might be enabling the species to flexibly control the level of bacitracin resistance for adaptive evolution.

Bacillus anthracis is a Gram-positive endospore-forming bacterial species that causes anthrax in both humans and animals. In Zambia, anthrax cases are frequently reported in both livestock and wildlife, with occasional transmission to humans, causing serious public health problems in the country. To understand the genetic diversity of B. anthracis strains in Zambia, we sequenced and compared the genomic DNA of B. anthracis strains isolated across the country. Single nucleotide polymorphisms clustered these strains into three groups. Genome sequence comparisons revealed a large deletion in strains belonging to one of the groups, possibly due to unequal crossing over between a pair of rRNA operons. The deleted genomic region included genes conferring resistance to bacitracin, and the strains with the deletion were confirmed with loss of bacitracin resistance. Similar deletions between rRNA operons were also observed in a few B. anthracis strains phylogenetically distant from Zambian strains. The structure of bacitracin resistance genes flanked by rRNA operons was conserved only in members of the Bacillus cereus group. The diversity and genomic characteristics of B. anthracis strains determined in this study would help in the development of genetic markers and treatment of anthrax in Zambia.

IMPORTANCE Anthrax is caused by Bacillus anthracis, an endospore-forming soil bacterium. The genetic diversity of B. anthracis is known to be low compared with that of Bacillus species. In this study, we performed whole-genome sequencing of Zambian isolates of B. anthracis to understand the genetic diversity between closely related strains. Comparison of genomic sequences revealed that closely related strains were separated into three groups based on single nucleotide polymorphisms distributed throughout the genome. A large genomic deletion was detected in the region containing a bacitracin resistance gene cluster flanked by rRNA operons, resulting in the loss of bacitracin resistance. The structure of the deleted region, which was also conserved among species of the Bacillus cereus group, has the potential for both deletion and amplification and thus might be enabling the species to flexibly control the level of bacitracin resistance for adaptive evolution.

Genotyping and phylogenetic placement of Bacillus anthracis isolates from Finland, a country with rare anthrax cases

Background

Anthrax, the zoonotic disease caused by the gram-positive bacterium Bacillus anthracis, is nowadays rare in northern parts of Europe including Finland and Scandinavia. Only two minor outbreaks of anthrax in 1988 and in 2004 and one sporadic infection in 2008 have been detected in animals in Finland since the 1970’s. Here, we report on two Finnish B. anthracis strains that were isolated from spleen and liver of a diseased calf related to the outbreak in 1988 (strain HKI4363/88) and from a local scrotum and testicle infection of a bull in 2008 (strain BA2968). These infections occurred in two rural Finnish regions, i.e., Ostrobothnia in western Finland and Päijänne Tavastia in southern Finland, respectively.

Results

The isolates were genetically characterized by PCR-based methods such as multilocus variable number of tandem repeat analysis (MLVA) and whole genome-sequence analysis (WGS). Phylogenetic comparison of the two strains HKI4363/88 and BA2968 by chromosomal single nucleotide polymorphism (SNP) analysis grouped these organisms within their relatives of the minor canonical A-branch canSNP-group A.Br.003/004 (A.Br.V770) or canonical B-branch B.Br.001/002, respectively. Strain HKI4363/88 clustered relatively closely with other members of the A.Br.003/004 lineage from Europe, South Africa, and South America. In contrast, strain BA2968 clearly constituted a new sublineage within B.Br.001/002 with its closest relative being HYO01 from South Korea.

Conclusions

Our results suggest that Finland harbors both unique (autochthonous) and more widely distributed, common clades of B. anthracis. We suspect that members of the common clades such as strains HKI4363/88 have been introduced only recently by anthropogenic activities involving importation of contaminated animal products. On the other hand, autochthonous strains such as isolate BA2968 probably have an older history of their introduction into Finland as evidenced by a high number of single nucleotide variant sites in their genomes.

Electronic supplementary material

The online version of this article (10.1186/s12866-018-1250-4) contains supplementary material, which is available to authorized users.

Pathogenicity, population genetics and dissemination of Bacillus anthracis

Bacillus anthracis, the etiological agent of anthrax, procures its particular virulence by a capsule and two AB type toxins: the lethal factor LF and the edema factor EF. These toxins primarily disable immune cells. Both toxins are translocated to the host cell by the adhesin-internalin subunit called protective antigen PA. PA enables LF to reach intra-luminal vesicles, where it remains active for long periods. Subsequently, LF translocates to non-infected cells, leading to inefficient late therapy of anthrax. B. anthracis undergoes slow evolution because it alternates between vegetative and long spore phases. Full genome sequence analysis of a large number of worldwide strains resulted in a robust evolutionary reconstruction of this bacterium, showing that B. anthracis is split in three main clades: A, B and C. Clade A efficiently disseminated worldwide underpinned by human activities including heavy intercontinental trade of goat and sheep hair. Subclade A.Br.WNA, which is widespread in the Northern American continent, is estimated to have split from clade A reaching the Northern American continent in the late Pleistocene epoch via the former Bering Land Bridge and further spread from Northwest southwards. An alternative hypothesis is that subclade A.Br.WNA. evolved from clade A.Br.TEA tracing it back to strains from Northern France that were assumingly dispatched by European explorers that settled along the St. Lawrence River. Clade B established mostly in Europe along the alpine axis where it evolved in association with local cattle breeds and hence displays specific geographic subclusters. Sequencing technologies are also used for forensic applications to trace unintended or criminal acts of release of B. anthracis. Under natural conditions, B. anthracis generally affects domesticated and wild ruminants in arid ecosystems. The more recently discovered B. cereus biovar anthracis spreads in tropical forests, where it threatens particularly endangered primate populations.

Best practice assessment of disease modelling for infectious disease outbreaks

During emerging disease outbreaks, public health, emergency management officials and decision-makers increasingly rely on epidemiological models to forecast outbreak progression and determine the best response to health crisis needs. Outbreak response strategies derived from such modelling may include pharmaceutical distribution, immunisation campaigns, social distancing, prophylactic pharmaceuticals, medical care, bed surge, security and other requirements. Infectious disease modelling estimates are unavoidably subject to multiple interpretations, and full understanding of a model's limitations may be lost when provided from the disease modeller to public health practitioner to government policymaker. We review epidemiological models created for diseases which are of greatest concern for public health protection. Such diseases, whether transmitted from person-to-person (Ebola, influenza, smallpox), via direct exposure (anthrax), or food and waterborne exposure (cholera, typhoid) may cause severe illness and death in a large population. We examine disease-specific models to determine best practices characterising infectious disease outbreaks and facilitating emergency response and implementation of public health policy and disease control measures.

The reporting of a Bacillus anthracis B-clade strain in South Africa after more than 20 years

OBJECTIVES

Anthrax is a disease with an age old history in Africa caused by the Gram-positive endospore forming soil bacterium Bacillus anthracis. Epizootics of wild ungulates occur annually in the enzootic region of Pafuri, Kruger National Park (KNP) in the Limpopo Province of South Africa. Rigorous routine surveillance and diagnostics in KNP, has not revealed these rare isolates since the 1990s, despite unabated annual outbreaks. In 2011 a cheetah was diagnosed as anthrax positive from a private game reserve in Limpopo Province and reported to State Veterinary Services for further investigation. Isolation, molecular diagnostics, whole genome sequencing and comparative genomics were carried out for B. anthracis KC2011.

RESULTS

Bacteriological and molecular diagnostics confirmed the isolate as B. anthracis. Subsequent typing and whole genome single nucleotide polymorphisms analysis indicated it clustered alongside B. anthracis SA A0091 in the B.Br.010 SNP branch. Unlike B. anthracis KrugerB strain, KC2011 strain has unique SNPs and represents a new branch in the B-clade. The isolation and genotypic characterisation of KC2011 demonstrates a gap in the reporting of anthrax outbreaks in the greater Limpopo province area. The identification of vulnerable and susceptible cheetah mortalities due to this strain has implications for conservation measures and disease control.

Persistent anthrax as a major driver of wildlife mortality in a tropical rainforest

Anthrax is a globally important animal disease and zoonosis. Despite this, our current knowledge of anthrax ecology is largely limited to arid ecosystems, where outbreaks are most commonly reported. Here we show that the dynamics of an anthrax-causing agent, Bacillus cereus biovar anthracis, in a tropical rainforest have severe consequences for local wildlife communities. Using data and samples collected over three decades, we show that rainforest anthrax is a persistent and widespread cause of death for a broad range of mammalian hosts. We predict that this pathogen will accelerate the decline and possibly result in the extirpation of local chimpanzee (Pan troglodytes verus) populations. We present the epidemiology of a cryptic pathogen and show that its presence has important implications for conservation.

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

[This corrects the article DOI: 10.1371/journal.pone.0146216.].

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.

Unexpected genomic relationships between Bacillus anthracis strains from Bangladesh and Central Europe

The zoonosis anthrax caused by the bacterium Bacillus anthracis has a broad geographical distribution. Active enzootic areas are typically located away from central and northern Europe where cases of the disease occur only sporadically and in limited numbers. In contrast, a few out of the 64 districts of Bangladesh are hyper-endemic for anthrax and there the disease causes major losses in live-stock. In this study we genotyped eight strains of B. anthracis collected from the districts of Sirajganj and Tangail in 2013. All these strains belonged to canSNP group A.Br.001/002 Sterne differing only in a few of 31 tandem-repeat (MLVA)-markers. Whole genome sequences were obtained from five of these strains and compared with genomic information of B. anthracis strains originating from various geographical locations. Characteristic signatures were detected defining two "Bangladesh" clusters potentially useful for rapid molecular epidemiology. From this data high-resolution PCR assays were developed and subsequently tested on additional isolates from Bangladesh and Central Europe. Remarkably, this comparative genomic analysis focusing on SNP-discovery revealed a close genetic relationship between these strains from Bangladesh and historic strains collected between 1991 and 2008 in The Netherlands and Germany, respectively. Possible explanations for these phylogenetic relationships are discussed.