Use of Canonical Single Nucleotide Polymorphism (CanSNPs) to characterize Bacillus anthracis outbreak strains in Zambia between 1990 and 2014

One Health Assessment of Bacillus anthracis Incidence and Detection in Anthrax-Endemic Areas of Pakistan

Anthrax, a severe zoonotic disease, is infrequently reported in anthrax-endemic regions of Pakistan. Despite clinical reports indicating its presence, particularly cutaneous anthrax, there is insufficient laboratory evidence regarding disease occurrence and environmental persistence. The present study aimed to confirm Bacillus anthracis presence, accountable for animal mortality and human infection, while exploring environmental transmission factors. Between March 2019 and July 2021, a total of 19 outbreaks were documented. Of these, 11 affected sheep/goats in Zhob district and 8 affected cattle/sheep in Bajour Agency. Clinical signs suggestive of Bacillus anthracis outbreak were observed in 11 animals. Blood and swab samples were collected for confirmation. The study followed a One Health approach, analyzing animal, environmental (soil/plant), and human samples. Of the 19 outbreaks, 11 were confirmed positive for anthrax based on growth characteristics, colony morphology, and PCR. Soil and plant root samples from the outbreak areas were collected and analyzed microscopically and molecularly. Cutaneous anthrax was observed in six humans, and swab samples were taken from the lesions. Human serum samples (n = 156) were tested for IgG antibodies against PA toxin and quantitative analysis of anthrax toxin receptor 1 (ANTXR1). Bacillus anthracis was detected in 65 out of 570 (11.40%) soil samples and 19 out of 190 (10%) plant root samples from the outbreak areas. Four out of six human samples from cutaneous anthrax lesions tested positive for Bacillus anthracis. Human anthrax seroprevalence was found to be 11% and 9% in two districts, with the highest rates among butchers and meat consumers. The highest ANTXR1 levels were observed in butchers, followed by meat consumers, farm employees, meat vendors, veterinarians, and farm owners. These findings highlight the persistence of anthrax in the region and emphasize the potential public health risks.

Population genomics of Bacillus anthracis from an anthrax hyperendemic area reveals transmission processes across spatial scales and unexpected within-host diversity

Genomic sequencing has revolutionized our understanding of bacterial disease epidemiology, but remains underutilized for zoonotic pathogens in remote endemic settings. Anthrax, caused by the spore-forming bacterium Bacillus anthracis, remains a threat to human and animal health and rural livelihoods in low- and middle-income countries. While the global genomic diversity of B. anthracis has been well-characterized, there is limited information on how its populations are genetically structured at the scale at which transmission occurs, critical for understanding the pathogen's evolution and transmission dynamics. Using a uniquely rich dataset, we quantified genome-wide SNPs among 73 B. anthracis isolates derived from 33 livestock carcasses sampled over 1 year throughout the Ngorongoro Conservation Area, Tanzania, a region hyperendemic for anthrax. Genome-wide SNPs distinguished 22 unique B. anthracis genotypes (i.e. SNP profiles) within the study area. However, phylogeographical structure was lacking, as identical SNP profiles were found throughout the study area, likely the result of the long and variable periods of spore dormancy and long-distance livestock movements. Significantly, divergent genotypes were obtained from spatio-temporally linked cases and even individual carcasses. The high number of SNPs distinguishing isolates from the same host is unlikely to have arisen during infection, as supported by our simulation models. This points to an unexpectedly wide transmission bottleneck for B. anthracis, with an inoculum comprising multiple variants being the norm. Our work highlights that inferring transmission patterns of B. anthracis from genomic data will require analytical approaches that account for extended and variable environmental persistence, as well as co-infection.

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.

Whole Genome Sequencing for Studying Bacillus anthracis from an Outbreak in the Abruzzo Region of Italy

Anthrax is a serious infectious disease caused by the gram-positive and spore-forming bacterium Bacillus anthracis. In Italy, anthrax is an endemic disease with sporadic cases each year and few outbreaks, especially in Southern Italy. However, new foci have been discovered in zones without previous history of anthrax. During summer 2016, an outbreak of anthrax caused the death of four goats in the Abruzzo region, where the disease had not been reported before. In order to investigate the outbreak, we sequenced one strain and compared it to 19 Italian B. anthracis genomes. Furthermore, we downloaded 71 whole genome sequences representing the global distribution of canonical SNP lineages and used them to verify the phylogenetic positioning. To this end, we analyzed and compared the genome sequences using canonical SNPs and the whole genome SNP-based analysis. Our results demonstrate that the outbreak strain belonged to the Trans-Eurasian (TEA) group A.Br.011/009, which is the predominant clade in Central-Southern Italy. In conclusion, the high genomic relatedness of the Italian TEA strains suggests their evolution from a common ancestor, while the spread is supposedly driven by trade as well as human and transhumance activities. Here, we demonstrated the capabilities of whole genome sequencing (WGS), which can be used as a tool for outbreak analyses and surveillance activities.

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.