Global genetic population structure of Bacillus anthracis

Genome sequences of two Bacillus anthracis strains utilized as veterinary vaccines in China

In this report, we present the complete genome sequences of two Bacillus anthracis strains utilized as veterinary vaccines in China. The sequencing was conducted using a hybrid assembly methodology that combined Illumina short reads and PacBio long reads. This approach provides a high-quality representative sequence for the strains mentioned above.

Antimicrobial Resistance Profile, Whole-Genome Sequencing and Core Genome Multilocus Sequence Typing of B. anthracis Isolates in Croatia from 2001 to 2022

Bacillus anthracis, the causative agent of anthrax disease, is a worldwide threat to livestock, wildlife and public health. It is also considered one of the most important pathogens of bioterrorism. Rapid and reliable diagnosis and administration of antimicrobials are essential for effective anthrax treatment. In this study, we determined the in vitro susceptibilities of 40 isolates of B. anthracis isolated in Croatia over the recent two decades to 18 antimicrobials. Whole-genome sequencing was performed, and bioinformatics tools were used to determine virulence factors and antimicrobial resistance genes. Core genome-based multilocus sequence typing was used for isolate comparison and phylogenetic analysis. All isolates were susceptible to all antimicrobials recommended for post-exposure prophylaxis or anthrax therapy. Susceptibility was found to all other tested antimicrobials that are an alternative for primary therapy. We found two beta-lactamase genes, but their expression is not sufficient to confer resistance. In all isolates used in this study, we found 21 virulence genes, 8 of which are responsible for toxin and capsule production. As far as phylogenetic analysis is concerned, the B. anthracis isolates from Croatia are categorised into two clades. The first is clade A, subclade Trans Eurasia, and the other is clade B, subclade B2.

Molecular characterization of an outbreak-involved Bacillus anthracis strain confirms the spillover of anthrax from West Africa

BACKGROUND

Anthrax, a zoonotic disease caused by the spore-forming bacterium Bacillus anthracis, remains a major global public health concern, especially in countries with limited resources. Sierra Leone, a West African country historically plagued by anthrax, has almost been out of report on this disease in recent decades. In this study, we described a large-scale anthrax outbreak affecting both animals and humans and attempted to characterize the pathogen using molecular techniques.

METHODS

The causative agent of the animal outbreak in Port Loko District, Sierra Leone, between March and May 2022 was identified using the nanopore sequencing technique. A nationwide active surveillance was implemented from May 2022 to June 2023 to monitor the occurrence of anthrax-specific symptoms in humans. Suspected cases were subsequently verified using quantitative polymerase chain reaction. Full-genome sequencing was accomplished by combining long-read and short-read sequencing methods. Subsequent phylogenetic analysis was performed based on the full-chromosome single nucleotide polymorphisms.

RESULTS

The outbreak in Port Loko District, Sierra Leone, led to the death of 233 animals between March 26th and May 16th, 2022. We ruled out the initial suspicion of Anaplasma species and successfully identified B. anthracis as the causative agent of the outbreak. As a result of the government's prompt response, out of the 49 suspected human cases identified during the one-year active surveillance, only 6 human cases tested positive, all within the first month after the official declaration of the outbreak. The phylogenetic analysis indicated that the BaSL2022 isolate responsible for the outbreak was positioned in the A.Br.153 clade within the TransEuroAsian group of B. anthracis.

CONCLUSIONS

We successfully identified a large-scale anthrax outbreak in Sierra Leone. The causative isolate of B. anthracis, BaSL2022, phylogenetically bridged other lineages in A.Br.153 clade and neighboring genetic groups, A.Br.144 and A.Br.148, eventually confirming the spillover of anthrax from West Africa. Given the wide dissemination of B. anthracis spores, it is highly advisable to effectively monitor the potential reoccurrence of anthrax outbreaks and to launch campaigns to improve public awareness regarding anthrax in Sierra Leone.

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.

Spatial and phylogenetic patterns reveal hidden infection sources of Bacillus anthracis in an anthrax outbreak in Son La province, Vietnam

Bacillus anthracis, the bacterial cause of anthrax, is a zoonosis affecting livestock and wildlife often spilling over into humans. In Vietnam, anthrax has been nationally reportable since 2015 with cases occurring annually, mostly in the northern provinces. In April 2022, an outbreak was reported in Son La province following the butchering of a water buffalo, Bubalus bubalis. A total of 137 humans from three villages were likely exposed to contaminated meat from the animal. Early epidemiological investigations suggested a single animal was involved in all exposures. Five B. anthracis isolates were recovered from human clinical cases along with one from the buffalo hide, another from associated maggots, and one from soil at the carcass site. The isolates were whole genome sequenced, allowing global, regional, and local molecular epidemiological analyses of the outbreak strains. All recovered B. anthracis belong to the A.Br.001/002 lineage based on canonical single nucleotide polymorphism analysis (canSNP). Although not previously identified in Vietnam, this lineage has been identified in the nearby countries of China, India, Indonesia, Thailand, as well as Australia. A twenty-five marker multi-locus variable number tandem repeat analysis (MLVA-25) was used to investigate the relationship between human, soil, and buffalo strains. Locally, four MLVA-25 genotypes were identified from the eight isolates. This level of genetic diversity is unusual for the limited geography and timing of cases and differs from past literature using MLVA-25. The coupled spatial and phylogenetic data suggest this outbreak originated from multiple, likely undetected, animal sources. These findings were further supported by local news reports that identified at least two additional buffalo deaths beyond the initial animal sampled in response to the human cases. Future outbreak response should include intensive surveillance for additional animal cases and additional molecular epidemiological traceback to identify pathogen sources.

Vibrio vulnificus mutation rate: an in vitro approach

Vibrio vulnificus is a multi-host pathogenic species currently subdivided into five phylogenetic lineages (L) plus one pathovar with the ability to infect fish due to a transmissible virulence plasmid. This plasmid (or a fragment of it) has been transmitted between lineages within the species, contributing to the evolution of V. vulnificus. This study aimed to provide an experimental approximation to the V. vulnificus mutation rate by determining spontaneous mutation rates from bacterial cultures of representants of the different lineages by whole-genome sequencing. To this purpose, synonymous SNP differences, i.e., spontaneous mutation not subjected to the evolutive forces, between initial and final culture after serial growth were evaluated and used for mutation rate calculation.

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.

Metagenomic Insights for Antimicrobial Resistance Surveillance in Soils with Different Land Uses in Brazil

Land-use conversion changes soil properties and their microbial communities, which, combined with the overuse of antibiotics in human and animal health, promotes the expansion of the soil resistome. In this context, we aimed to profile the resistome and the microbiota of soils under different land practices. We collected eight soil samples from different locations in the countryside of São Paulo (Brazil), assessed the community profiles based on 16S rRNA sequencing, and analyzed the soil metagenomes based on shotgun sequencing. We found differences in the communities' structures and their dynamics that were correlated with land practices, such as the dominance of Staphylococcus and Bacillus genera in agriculture fields. Additionally, we surveyed the abundance and diversity of antibiotic resistance genes (ARGs) and virulence factors (VFs) across studied soils, observing a higher presence and homogeneity of the vanRO gene in livestock soils. Moreover, three β-lactamases were identified in orchard and urban square soils. Together, our findings reinforce the importance and urgency of AMR surveillance in the environment, especially in soils undergoing deep land-use transformations, providing an initial exploration under the One Health approach of environmental levels of resistance and profiling soil communities.

Identification of Universally Applicable and Species-Specific Marker Peptides for Bacillus anthracis

Anthrax is a zoonotic infection caused by the bacterium Bacillus anthracis (BA). Specific identification of this pathogen often relies on targeting genes located on two extrachromosomal plasmids, which represent the major pathogenicity factors of BA. However, more recent findings show that these plasmids have also been found in other closely related Bacillus species. In this study, we investigated the possibility of identifying species-specific and universally applicable marker peptides for BA. For this purpose, we applied a high-resolution mass spectrometry-based approach for 42 BA isolates. Along with the genomic sequencing data and by developing a bioinformatics data evaluation pipeline, which uses a database containing most of the publicly available protein sequences worldwide (UniParc), we were able to identify eleven universal marker peptides unique to BA. These markers are located on the chromosome and therefore, might overcome known problems, such as observable loss of plasmids in environmental species, plasmid loss during cultivation in the lab, and the fact that the virulence plasmids are not necessarily a unique feature of BA. The identified chromosomally encoded markers in this study could extend the small panel of already existing chromosomal targets and along with targets for the virulence plasmids, may pave the way to an even more reliable identification of BA using genomics- as well as proteomics-based techniques.

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.

Allelic Polymorphism of Anthrax Pathogenicity Factor Genes as a Means of Estimating Microbiological Risks Associated with Climate Change

Abstract

Climate change brings new risks of emergence of especially dangerous diseases. The paper reports the possibility of assessing the pathogenic potential of bacteria as demonstrated by studying the allelic polymorphism of anthrax bacterium pathogenicity factor genes, which is a prerequisite for assessing the associated microbiological risks. The allelic polymorphism of the capBCADE operon (capB, capC, capA, capD, and capE genes) encoding the capsule biosynthesis proteins of Bacillus anthracis, and the acpA and acpB genes encoding the expression regulators of this operon have been studied for the first time. A number of single nucleotide polymorphisms (SNPs) were described in the strains of the studied sample, including 5 SNPs in the capB gene, 3 in capC, 4 in capA, 14 in capD, 2 in capE, and 15 in acpB, as well as 7 SNPs and one insertion in the acpA gene. As a result, the sample has been divided into sequence types for each gene and 17 genotypes, which are combinations of the identified sequence types. In silico translation of the detected alleles of the studied genes revealed three isoforms of the CapB and CapA proteins, two isoforms of the CapC and CapE proteins, six isoforms of the CapD protein, five isoforms of the AcpA protein, and four isoforms in the AcpB protein. It has been demonstrated that the SNP in the 351A → G position of capC is a marker of A.Br.Aust94 group strains. Based on the results, A.Br.Vollum group strains were divided into two subgroups. The strains in the evolutionary lines B and C differed from the line A strains by the presence of an 853G → A SNP in the acpA gene. In addition, a previously unknown variable number tandem repeat (VNTR), has been found in the acpA gene and the possibility of using it for differentiating and genotyping of B. anthracis strains has been demonstrated.

Rapid Capsular Antigen Immunoassay for Diagnosis of Inhalational Anthrax: Preclinical Studies and Evaluation in a Nonhuman Primate Model

Inhalational anthrax is a fatal infectious disease. Rapid and effective treatment is critically dependent on early and accurate diagnosis. Blood culture followed by identification and confirmation may take days to provide clinically relevant information. In contrast, immunoassay for a shed antigen, the capsular polypeptide gamma-d-polyglutamate (γDPGA), can provide rapid results at the point of care. In this study, a lateral flow immunoassay for γDPGA was evaluated in a robust nonhuman primate model of inhalational anthrax. The results showed that the time to a positive result with the rapid test using either serum or blood as a clinical specimen was similar to the time after infection when a blood culture became positive. In vitro testing showed that the test was equally sensitive with cultures of the three major clades of Bacillus anthracis. Cultures from other Bacillus spp. that are known to produce γDPGA also produced positive results. The test was negative with human sera from 200 normal subjects and 45 subjects with culture-confirmed nonanthrax bacterial or fungal sepsis. Taken together, the results showed that immunoassay for γDPGA is an effective surrogate for blood culture in a relevant cynomolgus monkey model of inhalational anthrax. The test would be a valuable aid in early diagnosis of anthrax, which is critical for rapid intervention and a positive outcome. Use of the test could facilitate triage of patients with signs and symptoms of anthrax in a mass-exposure incident and in low-resource settings where laboratory resources are not readily available. IMPORTANCE Patient outcome in anthrax is critically dependent on early diagnosis followed by effective treatment. We describe a rapid lateral flow immunoassay that detects capsular antigen of Bacillus anthracis that is shed into blood during infection. The test was evaluated in a robust cynomolgus monkey model of inhalational anthrax. Rapid detection of capsular antigen is an effective surrogate for the time-consuming and laboratory-intensive diagnosis by blood culture, direct fluorescent antibody staining, or other molecular testing. The test can be performed at the point of patient contact, is rapid and inexpensive, and can be used by individuals with minimal training.

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.

Target-enrichment sequencing yields valuable genomic data for challenging-to-culture bacteria of public health importance

Genomic data contribute invaluable information to the epidemiological investigation of pathogens of public health importance. However, whole-genome sequencing (WGS) of bacteria typically relies on culture, which represents a major hurdle for generating such data for a wide range of species for which culture is challenging. In this study, we assessed the use of culture-free target-enrichment sequencing as a method for generating genomic data for two bacterial species: (1) Bacillus anthracis, which causes anthrax in both people and animals and whose culture requires high-level containment facilities; and (2) Mycoplasma amphoriforme, a fastidious emerging human respiratory pathogen. We obtained high-quality genomic data for both species directly from clinical samples, with sufficient coverage (>15×) for confident variant calling over at least 80% of the baited genomes for over two thirds of the samples tested. Higher qPCR cycle threshold (Ct) values (indicative of lower pathogen concentrations in the samples), pooling libraries prior to capture, and lower captured library concentration were all statistically associated with lower capture efficiency. The Ct value had the highest predictive value, explaining 52 % of the variation in capture efficiency. Samples with Ct values ≤30 were over six times more likely to achieve the threshold coverage than those with a Ct > 30. We conclude that target-enrichment sequencing provides a valuable alternative to standard WGS following bacterial culture and creates opportunities for an improved understanding of the epidemiology and evolution of many clinically important pathogens for which culture is challenging.

Identification of the molecular characteristics of Bacillus anthracis (1982-2020) isolates in East Indonesia using multilocus variable-number tandem repeat analysis

Background and Aim:

Anthrax is one of the endemic strategic diseases in East Indonesia, particularly in the provinces of South Sulawesi, West Sulawesi, Gorontalo, East Nusa Tenggara, and West Nusa Tenggara. Anthrax is an important disease due to its zoonotic and economic impact on the livestock industry. This study aimed to identify the molecular characteristics of Bacillus anthracis in East Indonesia using multilocus variable-number tandem repeat (VNTR) analysis (MLVA).

Materials and Methods:

Isolates were obtained from an investigation of anthrax outbreaks in five provinces of East Indonesia from 1982 to 2020. Conventional polymerase chain reaction for B. anthracis was used to identify MLVA-8. Deoxyribonucleic acid sequencing analysis was based on MLVA-8 primers for VNTR identification of the phylogenetic relationship among 24 isolates of B. anthracis obtained from 17 distinct districts/cities in East Indonesia. Tandem Repeats Finder was used for VNTR identification, and Molecular Evolutionary Genetics Analysis X was used to construct phylogenetic analysis.

Results:

In this study, 24 isolates were classified as genotype or lineage A. There were four subgroups of B. anthracis circulating in East Indonesia based on eight molecular marker loci sequence results.

Conclusion:

The findings of this study show that MLVA-8 typing might be useful as a subtyping tool for the epidemiological investigation of identical genotypes and low genetic diversity of B. anthracis. No other lineage of B. anthracis was circulating in East Indonesia. Other molecular methods are needed, such as extended MLVA, whole-genome sequencing, and canonical single-nucleotide polymorphism, for a more precise study of B. anthracis genetic diversity.

Reoccurring bovine anthrax in Germany on the same pasture after 12 years

The zoonotic disease anthrax, caused by the endospore-forming bacterium Bacillus anthracis, is very rare in Germany. In the state of Bavaria, the last case occurred in July of 2009, resulting in four dead cows. In August of 2021, the disease reemerged after heavy rains, killing one gestating cow. Notably, both outbreaks affected the same pasture, suggesting a close epidemiological connection. B. anthracis could be grown from blood culture, and the presence of both virulence plasmids (pXO1 and pXO2) was confirmed by PCR. Also, recently developed diagnostic tools enabled rapid detection of B. anthracis cells and nucleic acids directly in clinical samples. The complete genome of the strain isolated from blood, designated BF-5, was DNA sequenced and phylogenetically grouped within the B.Br.CNEVA clade, which is typical for European B. anthracis strains. The genome was almost identical to BF-1, the isolate from 2009, separated only by three single nucleotide polymorphisms (SNPs) on the chromosome, one on plasmid pXO2 and three indel regions. Further, B. anthracis DNA was detected by PCR from soil samples taken from spots in the pasture where the cow had fallen. New tools based on phage receptor-binding proteins enabled the microscopic detection and isolation of B. anthracis directly from soil samples. These environmental isolates were genotyped and found to be identical to BF-5 in terms of SNPs. Therefore, it seems that the BF-5 genotype is currently the prevalent one at the affected premises. The area contaminated by the cadaver was subsequently disinfected with formaldehyde.

Molecular Typing of Ukrainian Bacillus anthracis Strains by Combining Whole-Genome Sequencing Techniques

Anthrax is a recurrent zoonosis in the Ukraine with outbreaks occurring repeatedly in certain areas. For determining whether several Bacillus anthracis genotypes are circulating in this region, four strains from various sources isolated from different regions of the Ukraine were investigated. By combining long- and short-read next-generation sequencing techniques, highly accurate genomes were reconstructed, enabling detailed in silico genotyping. Thus, the strains could be assigned to the Tsiankovskii subgroup of the “TransEurAsia” clade, which is commonly found in this region. Their high genetic similarity suggests that the four strains are members of the endemic population whose progenitor was once introduced in the Ukraine and bordering regions. This study provides information on B. anthracis strains from a region where there is little knowledge of the local population, thereby adding to the picture of global B. anthracis genotype distribution. We also emphasize the importance of surveillance and prevention methods regarding anthrax outbreaks, as other studies predicted a higher number of cases in the future due to global warming.

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.

In-Depth Analysis of Bacillus anthracis 16S rRNA Genes and Transcripts Reveals Intra- and Intergenomic Diversity and Facilitates Anthrax Detection

Analysis of 16S rRNA (rRNA) genes provides a central means of taxonomic classification of bacterial species. Based on presumed sequence identity among species of the Bacillus cereus sensu lato group, the 16S rRNA genes of B. anthracis have been considered unsuitable for diagnosis of the anthrax pathogen. With the recent identification of a single nucleotide polymorphism in some 16S rRNA gene copies, specific identification of B. anthracis becomes feasible. Here, we designed and evaluated a set of in situ, in vitro, and in silico assays to assess the unknown 16S state of B. anthracis from different perspectives. Using a combination of digital PCR, fluorescence in situ hybridization, long-read genome sequencing, and bioinformatics, we were able to detect and quantify a unique 16S rRNA gene allele of B. anthracis (16S-BA-allele). This allele was found in all available B. anthracis genomes and may facilitate differentiation of the pathogen from any close relative. Bioinformatics analysis of 959 B. anthracis SRA data sets inferred that abundances and genomic arrangements of the 16S-BA-allele and the entire rRNA operon copy numbers differ considerably between strains. Expression ratios of 16S-BA-alleles were proportional to the respective genomic allele copy numbers. The findings and experimental tools presented here provide detailed insights into the intra- and intergenomic diversity of 16S rRNA genes and may pave the way for improved identification of B. anthracis and other pathogens with diverse rRNA operons.

IMPORTANCE For severe infectious diseases, precise pathogen detection is crucial for antibiotic therapy and patient survival. Identification of Bacillus anthracis, the causative agent of the zoonosis anthrax, can be challenging when querying specific nucleotide sequences such as in small subunit rRNA (16S rRNA) genes, which are commonly used for typing of bacteria. This study analyzed on a broad genomic scale a cryptic and hitherto underappreciated allelic variant of the bacterium’s 16S rRNA genes and their transcripts using a set of in situ, in vitro, and in silico assays and found significant intra- and intergenomic heterogeneity in the distribution of the allele and overall rRNA operon copy numbers. This allelic variation was uniquely species specific, which enabled sensitive pathogen detection on both DNA and transcript levels. The methodology used here is likely also applicable to other pathogens that are otherwise difficult to discriminate from their less harmful relatives.

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.

Multiple stages of evolutionary change in anthrax toxin receptor expression in humans

The advent of animal husbandry and hunting increased human exposure to zoonotic pathogens. To understand how a zoonotic disease may have influenced human evolution, we study changes in human expression of anthrax toxin receptor 2 (ANTXR2), which encodes a cell surface protein necessary for Bacillus anthracis virulence toxins to cause anthrax disease. In immune cells, ANTXR2 is 8-fold down-regulated in all available human samples compared to non-human primates, indicating regulatory changes early in the evolution of modern humans. We also observe multiple genetic signatures consistent with recent positive selection driving a European-specific decrease in ANTXR2 expression in multiple tissues affected by anthrax toxins. Our observations fit a model in which humans adapted to anthrax disease following early ecological changes associated with hunting and scavenging, as well as a second period of adaptation after the rise of modern agriculture.

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.

The roles of environmental variation and parasite survival in virulence-transmission relationships

Disease outbreaks are a consequence of interactions among the three components of a host-parasite system: the infectious agent, the host and the environment. While virulence and transmission are widely investigated, most studies of parasite life-history trade-offs are conducted with theoretical models or tractable experimental systems where transmission is standardized and the environment controlled. Yet, biotic and abiotic environmental factors can strongly affect disease dynamics, and ultimately, host-parasite coevolution. Here, we review research on how environmental context alters virulence-transmission relationships, focusing on the off-host portion of the parasite life cycle, and how variation in parasite survival affects the evolution of virulence and transmission. We review three inter-related 'approaches' that have dominated the study of the evolution of virulence and transmission for different host-parasite systems: (i) evolutionary trade-off theory, (ii) parasite local adaptation and (iii) parasite phylodynamics. These approaches consider the role of the environment in virulence and transmission evolution from different angles, which entail different advantages and potential biases. We suggest improvements to how to investigate virulence-transmission relationships, through conceptual and methodological developments and taking environmental context into consideration. By combining developments in life-history evolution, phylogenetics, adaptive dynamics and comparative genomics, we can improve our understanding of virulence-transmission relationships across a diversity of host-parasite systems that have eluded experimental study of parasite life history.

Genetic diversity of Francisella tularensis subsp. holarctica in Kazakhstan

Tularemia is a highly dangerous zoonotic infection due to the bacteria Francisella tularensis. Low genetic diversity promoted the use of polymorphic tandem repeats (MLVA) as first-line assay for genetic description. Whole genome sequencing (WGS) is becoming increasingly accessible, opening the perspective of a time when WGS might become the universal genotyping assay. The main goal of this study was to describe F. tularensis strains circulating in Kazakhstan based on WGS data and develop a MLVA assay compatible with in vitro and in silico analysis. In vitro MLVA genotyping and WGS were performed for the vaccine strain and for 38 strains isolated in Kazakhstan from natural water bodies, ticks, rodents, carnivores, and from one migratory bird, an Isabellina wheatear captured in a rodent burrow. The two genotyping approaches were congruent and allowed to attribute all strains to two F. tularensis holarctica lineages, B.4 and B.12. The seven tandem repeats polymorphic in the investigated strain collection could be typed in a single multiplex PCR assay. Identical MLVA genotypes were produced by in vitro and in silico analysis, demonstrating full compatibility between the two approaches. The strains from Kazakhstan were compared to all publicly available WGS data of worldwide origin by whole genome SNP (wgSNP) analysis. Genotypes differing at a single SNP position were collected within a time interval of more than fifty years, from locations separated from each other by more than one thousand kilometers, supporting a role for migratory birds in the worldwide spread of the bacteria.

Genotyping of Francisella tularensis has become a routine practice in epidemiology. Despite rapidly accumulating knowledge, the phylogeography of the pathogen is still poorly understood and discussions about geographic and temporal origins continue. One important reason is the poor characterization of the pathogen in many tularemia-endemic countries. This article describes the genetic diversity of Francisella tularensis subsp. holarctica in Kazakhstan using tandem repeat polymorphisms as well as whole genome sequencing. Thirty-nine strains were analyzed and two lineages were identified, namely B.4 and B.12. The study demonstrates a wider distribution of genotype B.4 in Asia, and identified a more basal branching point in this subclade. The obtained data support the Asian origin hypothesis for F. tularensis. The finding of identical genotypes in strains separated in time by decades and a thousand-kilometers geographic distance, confirms the ability of the bacteria for long-term preservation and fast long distances spread. The isolation of F. tularensis subsp. holarctica from the bird species Isabellina wheatear allows speculating about a major contribution of birds to the phylogeography of the pathogen. A genotyping protocol was developed utilizing seven polymorphic tandem repeats, two of which were identified within the framework of this work. The in vitro and in silico results are identical when using sequencing reads of 300 base-pairs or more.

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.

Investigation on Anthrax in Bangladesh during the Outbreaks of 2011 and Definition of the Epidemiological Correlations

In 2011, in Bangladesh, 11 anthrax outbreaks occurred in six districts of the country. Different types of samples were collected from May to September in the six districts where anthrax had occurred in order to detect and type Bacillus anthracis (B. anthracis) strains. Anthrax was detected in 46.6% of the samples analysed, in particular in soils, but also in bone samples, water, animal feed, and rumen ingesta of dead animals. Canonical single nucleotide polymorphisms (CanSNPs) analysis showed that all the isolates belonged to the major lineage A, sublineage A.Br.001/002 of China and Southeast Asia while the multi-locus variable number of tandem repeats (VNTRs) analysis (MLVA) with 15 VNTRs demonstrated the presence of five genotypes, of which two resulted to be new genotypes. The single nucleotide repeats (SNRs) analysis showed 13 SNR types; nevertheless, due to its higher discriminatory power, the presence of two isolates with different SNR-type polymorphisms was detected within two MLVA genotypes. This study assumes that soil is not the only reason for the spread of the disease in Bangladesh; contaminated feed and water can also play an important role in the epidemiology of anthrax. Possible explanations for these epidemiological relationships are discussed.

Hidden Diversity within Common Protozoan Parasites as Revealed by a Novel Genomotyping Scheme

Giardia duodenalis (syn. Giardia lamblia, Giardia intestinalis) is the causative agent of giardiasis, one of the most common diarrheal infections in humans. Evolutionary relationships among G. duodenalis genotypes (or subtypes) of assemblage B, one of two genetic assemblages causing the majority of human infections, remain unclear due to poor phylogenetic resolution of current typing methods. In this study, we devised a methodology to identify new markers for a streamlined multilocus sequence typing (MLST) scheme based on comparisons of all core genes against the phylogeny of whole-genome sequences (WGS). Our analysis identified three markers with resolution comparable to that of WGS data. Using newly designed PCR primers for our novel MLST loci, we typed an additional 68 strains of assemblage B. Analyses of these strains and previously determined genome sequences showed that genomes of this assemblage can be assigned to 16 clonal complexes, each with unique gene content that is apparently tuned to differential virulence and ecology. Obtaining new genomes of Giardia spp. and other eukaryotic microbial pathogens remains challenging due to difficulties in culturing the parasites in the laboratory. Hence, the methods described here are expected to be widely applicable to other pathogens of interest and advance our understanding of their ecology and evolution.IMPORTANCE Giardia duodenalis assemblage B is a major waterborne pathogen and the most commonly identified genotype causing human giardiasis worldwide. The lack of morphological characters for classification requires the use of molecular techniques for strain differentiation; however, the absence of scalable and affordable next-generation sequencing (NGS)-based typing methods has prevented meaningful advancements in high-resolution molecular typing for further understanding of the evolution and epidemiology of assemblage B. Prior studies have reported high sequence diversity but low phylogenetic resolution at standard loci in assemblage B, highlighting the necessity of identifying new markers for accurate and robust molecular typing. Data from comparative analyses of available genomes in this study identified three loci that together form a novel high-resolution typing scheme with high concordance to whole-genome-based phylogenomics and which should aid in future public health endeavors related to this parasite. In addition, data from newly characterized strains suggest evidence of biogeographic and ecologic endemism.

TaqMan Assays for Simultaneous Detection of Bacillus anthracis and Bacillus cereus biovar anthracis

Anthrax is a worldwide zoonotic disease caused by the spore-forming bacterium Bacillus anthracis. Primarily a disease of herbivores, human infections often result from direct contact with contaminated animal products (cutaneous and inhalational anthrax) or through consumption of infected meat (gastrointestinal anthrax). The genetic near neighbor, Bacillus cereus biovar anthracis (Bcbva), causes an anthrax-like illness in the wildlife and livestock of west and central Africa due to the presence and expression of B. anthracis-specific virulence factors in this background. While Bcbva infections have not been reported in humans, a recent seroprevalence study detected Bcbva antibodies in the rural population around Taï National Park. This work describes the development of new TaqMan multiplex PCRs for the simultaneous detection of B. anthracis and Bcbva. The assays are designed to amplify Ba-1, capB, and lef markers in B. anthracis and genomic island IV (GI4), capB, and lef in Bcbva. Our assays allow for the rapid discrimination of B. anthracis and Bcbva and will provide insights into the molecular epidemiology of these two important pathogens that share an overlapping geographical range in west and central Africa.

Convergent evolution of diverse Bacillus anthracis outbreak strains toward altered surface oligosaccharides that modulate anthrax pathogenesis

Bacillus anthracis, a spore-forming gram-positive bacterium, causes anthrax. The external surface of the exosporium is coated with glycosylated proteins. The sugar additions are capped with the unique monosaccharide anthrose. The West African Group (WAG) B. anthracis have mutations rendering them anthrose deficient. Through genome sequencing, we identified 2 different large chromosomal deletions within the anthrose biosynthetic operon of B. anthracis strains from Chile and Poland. In silico analysis identified an anthrose-deficient strain in the anthrax outbreak among European heroin users. Anthrose-deficient strains are no longer restricted to West Africa so the role of anthrose in physiology and pathogenesis was investigated in B. anthracis Sterne. Loss of anthrose delayed spore germination and enhanced sporulation. Spores without anthrose were phagocytized at higher rates than spores with anthrose, indicating that anthrose may serve an antiphagocytic function on the spore surface. The anthrose mutant had half the LD50 and decreased time to death (TTD) of wild type and complement B. anthracis Sterne in the A/J mouse model. Following infection, anthrose mutant bacteria were more abundant in the spleen, indicating enhanced dissemination of Sterne anthrose mutant. At low sample sizes in the A/J mouse model, the mortality of ΔantC-infected mice challenged by intranasal or subcutaneous routes was 20% greater than wild type. Competitive index (CI) studies indicated that spores without anthrose disseminated to organs more extensively than a complemented mutant. Death process modeling using mouse mortality dynamics suggested that larger sample sizes would lead to significantly higher deaths in anthrose-negative infected animals. The model was tested by infecting Galleria mellonella with spores and confirmed the anthrose mutant was significantly more lethal. Vaccination studies in the A/J mouse model showed that the human vaccine protected against high-dose challenges of the nonencapsulated Sterne-based anthrose mutant. This work begins to identify the physiologic and pathogenic consequences of convergent anthrose mutations in B. anthracis.

Changed epidemiology of anthrax and molecular characteristics of Bacillus anthracis in Inner Mongolia Autonomous Region, China

Anthrax is a natural foci disease in Inner Mongolia, which poses a severe threat to public health. In this study, the incidence number, rate and constituent ratio were used to describe the epidemiological characteristics of anthrax in the region from 1956-2018. The molecular correlation and genetic characteristics of the strains were investigated using canonical single nucleotide polymorphisms (CanSNP), multiple-locus variable-number tandem repeat analysis (MLVA-15) and whole genome sequencing (WGS). The epidemiological characteristics of anthrax in Inner Mongolia have altered significantly. The incidence of anthrax has decreased annually without vaccination, and the regional distribution of anthrax gradually transferred from central and western regions to the eastern. Moreover, the occupation distribution evolved from multiple early occupations to predominated by farmers and herdsmen. This change is closely related to policy factors and to changes in the means of production and the living habits of the local population. This indicates that reformulating the control and prevention strategies is essential. Both A. Br. Ames and A. Br. 001/002 subgroups were the predominant CanSNP genotypes of Bacillus anthracis in Inner Mongolia. A total of 36 strains constituted six shared MLVA-15 genotypes, suggesting an epidemiological link between the strains of each shared genotype. The six shared genotypes ([GT1, 9, 11 and 15] and [GT8 and 12]) consisting of 2-7 strains confirmed the occurrence of multiple point outbreaks and cross-regional transmission caused by multiple common sources of infection. Phylogenetic analysis based on the WGS core genome showed that strains from this study formed an independent clade (C.V.), and they were positioned close to each other, suggesting a common origin. Further comparison analysis should be performed to ascertain the geographic origin of these strains.

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.

Rapid Microscopic Detection of Bacillus anthracis by Fluorescent Receptor Binding Proteins of Bacteriophages

Bacillus anthracis, the etiological agent of anthrax disease, is typically diagnosed by immunological and molecular methods such as polymerase chain reaction (PCR). Alternatively, mass spectrometry techniques may aid in confirming the presence of the pathogen or its toxins. However, because of the close genetic relationship between B. anthracis and other members of the Bacillus cereus sensu lato group (such as Bacillus cereus or Bacillus thuringiensis) mis- or questionable identification occurs frequently. Also, bacteriophages such as phage gamma (which is highly specific for B. anthracis) have been in use for anthrax diagnostics for many decades. Here we employed host cell-specific receptor binding proteins (RBP) of (pro)-phages, also known as tail or head fibers, to develop a microscopy-based approach for the facile, rapid and unambiguous detection of B. anthracis cells. For this, the genes of (putative) RBP from Bacillus phages gamma, Wip1, AP50c and from lambdoid prophage 03 located on the chromosome of B. anthracis were selected. Respective phage genes were heterologously expressed in Escherichia coli and purified as fusions with fluorescent proteins. B. anthracis cells incubated with either of the reporter fusion proteins were successfully surface-labeled. Binding specificity was confirmed as RBP fusion proteins did not bind to most isolates of a panel of other B. cereus s.l. species or to more distantly related bacteria. Remarkably, RBP fusions detected encapsulated B. anthracis cells, thus RBP were able to penetrate the poly-γ-d-glutamate capsule of B. anthracis. From these results we anticipate this RBP-reporter assay may be useful for rapid confirmative identification of B. anthracis.

Genetic Diversity of Australian Bacillus anthracis Isolates Revealed by Multiple-Locus Variable-Number Tandem Repeat Analysis

Outbreaks of anthrax occur sporadically in Australia and most commonly in the "anthrax belt", a region which extends from southern Queensland through the centre of New South Wales and into northern Victoria. Little is known about the epidemiological links between Bacillus anthracis isolates taken from different outbreaks and the diversity of strains within Australia. We used multiple-locus variable-number tandem repeat analysis employing 25 markers (MLVA25) to genotype 99 B. anthracis isolates from an archival collection of Australian isolates. MLVA25 genotyping revealed eight unique genotypes which clustered within the previously defined A3 genotype of B. anthracis. Genotyping of B. anthracis strains from outbreaks of disease in Victoria identified the presence of multiple genotypes associated with these outbreaks. The geographical distribution of genotypes within Australia suggests that a single genotype was introduced into the eastern states of Australia, followed by the spread and localised differentiation of the pathogen (MLVA25 genotypes MG1-MG6) throughout the anthrax belt. In contrast, unexplained occurrences of disease in areas outside of this anthrax belt which are associated with different genotypes, (MLVA25 genotypes MG7 and MG8) indicate separate introductions of B. anthracis into Australia.

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.

Epidemiological Investigation and Etiological Analysis of a Cutaneous Anthrax Epidemic Caused by Butchering Sick Cattle in Guizhou, China

A suspected human cutaneous anthrax epidemic caused by butchering sick cattle occurred in Zhijin County of Guizhou Province, Southwest of China, in 2016. Epidemiological investigation and etiological analysis were performed to provide a scientific basis for the source tracking of the epidemic. The epidemic was epidemiologically investigated, and skin blister samples collected from patients and soil samples collected from the butchering spots were used for Bacillus anthracis isolation. The suspicious B. anthracis isolates were identified using conventional methods and PCR, followed by genotyping using multiple-locus variable-number tandem repeats (VNTRs) analysis (MLVA-15) and canonical single-nucleotide polymorphism (canSNP). The genetic relationship of epidemic strains and isolates collected from other regions was analyzed. Epidemiological investigation results showed that the patients may be infected by B. anthracis during butchering sick cattle. Two suspected B. anthracis strains were isolated from blood samples and blister fluids, respectively. Conventional methods identified the two suspected isolates as B. anthracis, while PCR results showed that anti-protective antigen (PA) and capsule (CAP) gene were positive in the two isolates. MLVA-15 showed that the MLVA profiles of the two isolates were 9-20-12-53-16-2-8-8-8-4-4-4-4-10-4, which is different from the MLVA profiles of representative strains from other regions. CanSNP analysis showed that the two strains belonged to cluster A.Br.001/002. Clustering analysis and minimum spanning tree (MST) demonstrated that the two isolates were clustered with strains previously isolated from Guizhou Province. The results indicated that B. anthracis was the pathogen for this epidemic, and the patients were infected during butchering the sick. The genetic characteristics and the relationship of the B. anthracis isolates to strains from other regions indicated that the epidemic was a local occurrence.

Coalescence modeling of intrainfection Bacillus anthracis populations allows estimation of infection parameters in wild populations

Bacillus anthracis, the etiological agent of anthrax, is a well-established model organism. For B. anthracis and most other infectious diseases, knowledge regarding transmission and infection parameters in natural systems, in large part, comprises data gathered from closely controlled laboratory experiments. Fatal, natural anthrax infections transmit the bacterium through new host-pathogen contacts at carcass sites, which can occur years after death of the previous host. For the period between contact and death, all of our knowledge is based upon experimental data from domestic livestock and laboratory animals. Here we use a noninvasive method to explore the dynamics of anthrax infections, by evaluating the terminal diversity of B. anthracis in anthrax carcasses. We present an application of population genetics theory, specifically, coalescence modeling, to intrainfection populations of B. anthracis to derive estimates for the duration of the acute phase of the infection and effective population size converted to the number of colony-forming units establishing infection in wild plains zebra (Equus quagga). Founding populations are small, a few colony-forming units, and infections are rapid, lasting roughly between 1 d and 3 d in the wild. Our results closely reflect experimental data, showing that small founding populations progress acutely, killing the host within days. We believe this method is amendable to other bacterial diseases from wild, domestic, and human systems.

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.

Genetic diversity of Bacillus anthracis Ames lineage strains in China

BACKGROUND

Anthrax is an endemic disease that persists in the rural regions of China. The global genetic population structure of B.anthracis has also been defined by the canonical single-nucleotide polymorphisms (canSNP) and multiple-locus variable-number tandem repeat analysis (MLVA). Five canSNP lineages were found in China, and the A.Br.Ames lineage has been the second predominant group in recent years in China. The objective of this study was to reveal genetic diversity of the Ames lineage strains by MLVA.

METHODS

Two molecular typing methods, canSNP and MLVA with 15markers were used to study the genetic relationship among the Ames lineage strains. The outbreak information associated with these strains was also collected and investigated.

RESULTS

From 2007 to 2018, a total of 21 human anthrax infection outbreaks (68 patients) associated with B. anthracis Ames lineage strains were reported in China. Ames lineage strain-associated human anthrax is mainly distributed in the northern part of China, including the provinces of Inner Mongolia, Liaoning, Gansu, and Xinjiang. In the study, a total of 30 Ames lineage strains were included and 10 MLVA15 genotypes were identified. These strains were mainly found in northeast China, Inner Mongolia and Liaoning. In recent years, the Ames lineage strains were isolated in the two provinces every year. The 18 Ames lineage strains isolated from Inner Mongolia were divided into eight MLVA15 genotypes. From 2010 to 2015, there were continuous reports of outbreaks in Keyouzhongqi County, Inner Mongolia, and the strains that were isolated annually in succession belonged to the MLVA15-30 genotype.

CONCLUSIONS

The Ames lineage strains are widely distributed in northern China. Their genetic diversity can be illustrated by the results of the MLVA. The genetic characteristics of the Ames lineage strains from outbreaks in different provinces varied. In some areas, human anthrax outbreaks occurred annually in succession, and these related strains grouped together. These observations indicate that the local environment was persistently contaminated with B. anthracis spores, vaccination of livestock should become the fundamental control measure in the areas.

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.

Laboratory strains of Bacillus anthracis lose their ability to rapidly grow and sporulate compared to wildlife outbreak strains

Bacillus anthracis is the causative agent of anthrax in animals and humans. The organism lies in a dormant state in the soil until introduced into an animal via, ingestion, cutaneous inoculation or inhalation. Once in the host, spores germinate into rapidly growing vegetative cells elaborating toxins. When animals die of anthrax, vegetative bacteria sporulate upon nutrient limitation in the carcass or soil while in the presence of air. After release into the soil environment, spores form a localized infectious zone (LIZ) at and around the carcass. Laboratory strains of B. anthracis produce fewer proteins associated with growth and sporulation compared to wild strains isolated from recent zoonotic disease events. We verified wild strains grow more rapidly than lab strains demonstrating a greater responsiveness to nutrient availability. Sporulation was significantly more rapid in these wild strains compared to lab strains, indicating wild strains are able to sporulate faster due to nutrient limitation while laboratory strains have a decrease in the speed at which they utilize nutrients and an increase in time to sporulation. These findings have implications for disease control at the LIZ as well as on the infectious cycle of this dangerous zoonotic pathogen.

Genotyping and population diversity of Bacillus anthracis in China based on MLVA and canSNP analysis

In pastoral parts of China, anthrax still presents a major risk to livestock and threatens the health of local human populations. Currently, whole-genome-based molecular markers, such as single-nucleotide polymorphisms (SNPs) and variable number tandem repeats (VNTRs), are the most effective tools for genotyping Bacillus anthracis. In this study, 191 isolates were selected to assess the diversity of B. anthracis in China. Five isolates were confirmed not to be B. anthracis by clustered regularly interspaced short palindromic repeat analysis, while the remaining 186 isolates were typed using canonical SNP (canSNP) and VNTR analyses. Five sublineages/subgroups, A.Br.001/002, A.Br.Vollum, A.Br.Aust.94, A.Br.Ames, and A.Br.008/009, were detected based on 13 canSNP sites. The 186 isolates were further assigned 114 sequence types based on 27 VNTR loci, with major branches correlating with the canSNP analysis. We then used a simplified multiple-locus variable number tandem repeat analysis (MLVA) protocol (MLVAmin) based on eight high-resolution VNTR sites to analyze the Chinese isolates, with the resulting phylogeny again agreeing with the canSNP analysis. We also developed two schemes, MLVAc and MLVAp, using various numbers of VNTRs to analyze different canSNP sublineages to increase the typing resolution of the canSNP protocol. The results showed a highly imbalanced geographical distribution of the B. anthracis population, with four different sublineages observed in Xinjiang Province, while only one sublineage, A.Br.001/002, was found in the other six provinces, except for three A.Br.Ames strains isolated from Inner Mongolia. Based on the MLVA and canSNP analysis, the spread of B. anthracis appears to have occurred from west to east via three independent routes.

Genetic characterization of Bacillus anthracis strains circulating in Italy from 1972 to 2018

In Italy anthrax is an endemic disease, with a few outbreaks occurring almost every year. We surveyed 234 B. anthracis strains from animals (n = 196), humans (n = 3) and the environment (n = 35) isolated during Italian outbreaks in the years 1972-2018. Despite the considerable genetic homogeneity of B. anthracis, the strains were effectively differentiated using canonical single nucleotide polymorphisms (CanSNPs) assay and multiple-locus variable-number tandem repeat analysis (MLVA). The phylogenetic identity was determined through the characterization of 14 CanSNPs. In addition, a subsequent 31-loci MLVA assay was also used to further discriminate B. anthracis genotypes into subgroups. The analysis of 14 CanSNPs allowed for the identification of four main lineages: A.Br.011/009, A.Br.008/011 (respectively belonging to A.Br.008/009 sublineage, also known Trans-Eurasian or TEA group), A.Br.005/006 and B.Br.CNEVA. A.Br.011/009, the most common subgroup of lineage A, is the major genotype of B. anthracis in Italy. The MLVA analysis revealed the presence of 55 different genotypes in Italy. Most of the genotypes are genetically very similar, supporting the hypothesis that all strains evolved from a local common ancestral strain, except for two genotypes representing the branch A.Br.005/006 and B.Br.CNEVA. The genotyping analysis applied in this study remains a very valuable tool for studying the diversity, evolution, and molecular epidemiology of B. anthracis.

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.

Nucleotide polymorphism assay for the identification of west African group Bacillus anthracis: a lineage lacking anthrose

Background

The exosporium of the anthrax-causing Bacillus anthracis endospores display a tetrasaccharide composed of three rhamnose residues and an unusual sugar termed anthrose. Anthrose is a proposed potential target for immunotherapy and for specific detection of B. anthracis. Although originally thought to be ubiquitous in B. anthracis, previous work identified an anthrose negative strain from a West African lineage isolated from cattle that could represent a vaccine escape mutant. These strains carry genes required for expression of the anthrose operon but premature stop codons resulting from an 8-bp insertion in BAS3320 (an amino-transferase) and a C/T substitution at position 892 of the BAS3321 (a glycosyltransferase) gene prevent anthrose expression. Various other single nucleotide polymorphisms (SNPs) have been identified throughout the operon and could be the basis for detection of anthrose-deficient strains.

Results

In this study, we evaluated rhAmp genotypic assays based on SNPs at positions 892 and 1352 of BAS3321 for detection and differentiation of anthrose negative (Ant⁻) West African strains. Discrimination of anthrose negative West African isolates was achieved with as low as 100 fg of DNA, whereas consistent genotyping of Sterne necessitated at least 1 pg of DNA.

Conclusions

Screening of a global panel of B. anthracis isolates showed anthrose-expressing alleles are prevalent worldwide whereas the anthrose-deficient phenotype is to date limited to West Africa. Our work also revealed a third, previously unreported anthrose genotype in which the operon is altogether missing from a Polish B. anthracis isolate.

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.

Biological Properties and Molecular Genetic Characteristics of Bacillus Anthracis Strains Isolated During Anthrax Outbreaks in Kazakhstan in 2016

Today, anthrax is recorded in many countries around the world, in Kazakhstan – in the form of sporadic cases or small outbreaks. Despite the endemicity of anthrax in Kazakhstan, the genetic diversity strains is not well described. At the moment, MLVA-25, which is sufficient for molecular and epidemiological monitoring, is recognized as one of the most discriminatory methods of genotyping.Objective: to determine the cultural and morphological properties, the genomic characteristics of the strains of the anthrax pathogen, their geographical distribution in the territory of Kazakhstan during the outbreaks in 2016, a comparative analysis with collection strains isolated since 1962.Methods: microbiological, genetic research methods were used in the work.Results. We investigated 11 strains of B. anthracis, which were isolated in 2016 in Kazakhstan. For comparison of genotypes, 26 strains were used from the collection of pathogenic microorganisms of the A.M. Aykimbaev’s Kazakh Scientific Center for Quarantine and Zoonotic Diseases. Phylogenetic analysis clustered 37 strains of B. anthracis into three clusters and 23 genotypes.The results of studying the phenotypic properties of anthrax strains by the main identification tests showed that all the strains studied in the experiment and isolated from 1961 to 2016 had biological properties characteristic of typical of B. anthracis strains.The anthrax strains isolated in 2016 are characterized as significant variability, as well as the circulation of the same genotypes and clusters in different areas of Kazakhstan. MLVA-profiles of analyzed Kazakhstan strains are unique and do not fully coincide with any studied strain from MLVAbank. On the MST-tree, Kazakhstan’s strains are located in three clusters, as on the phylogenetic tree.Conclusion: Molecular genetic analysis of B. anthracis strains enhances the ability of epidemiologists to track the sources and pathways of infection.It is necessary to improve the tracking system for strains of especially dangerous infections in  Kazakhstan using modern molecular genetic methods.

The signatures of microorganisms and of human and environmental biomes can now be used to provide evidence in legal cases

The microorganisms with which we share our world go largely unnoticed. We are, however, beginning to be able to exploit their apparently silent presence as witnesses to events that are of legal concern. This information can be used to link forensic samples to criminal events and even perpetrators. Once dead, our bodies are rapidly colonised, internally and externally. The progress of these events can be charted to inform how long and even by what means a person has died. A small number of microbial species could actually be the cause of such deaths as a result of biocrime or bioterrorism. The procedures and techniques to respond to such attacks have matured in the last 20 years. The capability now exists to identify malicious intent, characterise the threat agent to isolate level and potentially link it to perpetrators with a high level of confidence.