Genetic diversity in a Bacillus anthracis historical collection (1954 to 1988)

Analysis of Whole-Genome Sequences for the Prediction of Penicillin Resistance and β-Lactamase Activity in Bacillus anthracis

Penicillin (PEN) is a low-cost option for anthrax treatment, but naturally occurring resistance has been reported. β-Lactamase expression (bla1, bla2) in Bacillus anthracis is regulated by a sigma factor (SigP) and its cognate anti-sigma factor (RsiP). Mutations leading to truncation of RsiP were previously described as a basis for PEN resistance. Here, we analyze whole-genome sequencing (WGS) data and compare the chromosomal sigP-bla1 regions from 374 B. anthracis strains to determine the frequency of mutations, identify mutations associated with PEN resistance, and evaluate the usefulness of WGS for predicting PEN resistance. Few (3.5%) strains contained at least 1 of 11 different mutations in sigP, rsiP, or bla1. Nine of these mutations have not been previously associated with PEN resistance. Four strains showed PEN resistance (PEN-R) by conventional broth microdilution, including 1 strain with a novel frameshift in rsiP. One strain that carries the same rsiP frameshift mutation as that found previously in a PEN-R strain showed a PEN-susceptible (PEN-S) phenotype and exhibited decreased bla1 and bla2 transcription. An unexpectedly small colony size, a reduced growth rate, and undetectable β-lactamase activity levels (culture supernatant and cell lysate) were observed in this PEN-S strain. Sequence analysis revealed mutations in genes associated with growth defects that may contribute to this phenotype. While B. anthracis rsiP mutations cannot be exclusively used to predict resistance, four of the five strains with rsiP mutations were PEN-R. Therefore, the B. anthracis sigP-bla1 region is a useful locus for WGS-based PEN resistance prediction, but phenotypic testing remains essential. IMPORTANCE Determination of antimicrobial susceptibility of B. anthracis is essential for the appropriate distribution of antimicrobial agents for postexposure prophylaxis (PEP) and treatment of anthrax. Analysis of WGS data allows for the rapid detection of mutations in antimicrobial resistance (AMR) genes in an isolate, but the presence of a mutation in an AMR gene does not always accurately predict resistance. As mutations in the anti-sigma factor RsiP have been previously associated with high-level penicillin resistance in a limited number of strains, we investigated WGS assemblies from 374 strains to determine the frequency of mutations and performed functional antimicrobial susceptibility testing. Of the five strains that contained mutations in rsiP, only four were PEN-R by functional antimicrobial susceptibility testing. We conclude that while sequence analysis of this region is useful for AMR prediction in B. anthracis, genetic analysis should not be used exclusively and phenotypic susceptibility testing remains essential.

Unexpected Relations of Historical Anthrax Strain

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

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

Optical Screening for Rapid Antimicrobial Susceptibility Testing and for Observation of Phenotypic Diversity among Strains of the Genetically Clonal Species Bacillus anthracis

During high-impact events involving Bacillus anthracis, such as the Amerithrax incident of 2001 or the anthrax outbreaks in Russia and Sweden in 2016, critical decisions to reduce morbidity and mortality include rapid selection and distribution of effective antimicrobial agents for treatment and postexposure prophylaxis. Detection of antimicrobial resistance currently relies on a conventional broth microdilution method that requires a 16- to 20-h incubation time for B. anthracis Advances in high-resolution optical screening offer a new technology to more rapidly evaluate antimicrobial susceptibility and to simultaneously assess the growth characteristics of an isolate. Herein, we describe a new method developed and evaluated as a rapid antimicrobial susceptibility test for B. anthracis This method is based on automated digital time-lapse microscopy to observe the growth and morphological effects of relevant antibiotics with an optical screening instrument, the oCelloScope. B. anthracis strains were monitored over time in the presence or absence of penicillin, ciprofloxacin, or doxycycline. Susceptibility to each antibiotic was determined in ≤4 h, 75 to 80% less than the time required for conventional methods. Time-lapse video imaging compiled from the optical screening images revealed unexpected differences in growth characteristics among strains of B. anthracis, which is considered to be a clonal organism. This technology provides a new approach for rapidly detecting phenotypic antimicrobial resistance and for documenting growth attributes that may be beneficial in the further characterization of individual strains.

Recent developments in the understanding and use of anthrax vaccine adsorbed: achieving more with less

Anthrax Vaccine Adsorbed (AVA, BioThrax™) is the only Food and Drug Administration (FDA) approved vaccine for the prevention of anthrax in humans. Recent improvements in pre-exposure prophylaxis (PrEP) use of AVA include intramuscular (IM) administration and simplification of the priming series to three doses over 6 months. Administration IM markedly reduced the frequency, severity and duration of injection site reactions. Refinement of animal models for inhalation anthrax, identification of immune correlates of protection and cross-species modeling have created opportunities for reductions in the PrEP booster schedule and were pivotal in FDA approval of a post-exposure prophylaxis (PEP) indication. Clinical and nonclinical studies of accelerated PEP schedules and divided doses may provide prospects for shortening the PEP antimicrobial treatment period. These data may assist in determining feasibility of expanded coverage in a large-scale emergency when vaccine demand may exceed availability. Enhancements to the AVA formulation may broaden the vaccine's PEP application.

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

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

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

Background

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

Results

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

Conclusions

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

Antimicrobial susceptibility and molecular subtyping of 55 Turkish Bacillus anthracis strains using 25-loci multiple-locus VNTR analysis

Anthrax, which is caused by the bacterium Bacillus anthracis, is one of the oldest documented infectious diseases in both livestock and humans. The differentiation of B. anthracis strains is difficult because of their highly homogeneous genomes. We used multiple-locus variable-number tandem repeat analysis (MLVA) with 25 markers to genotype 55 B. anthracis isolates from 16 distinct regions of Turkey. The antimicrobial susceptibility of the isolates was investigated using the agar dilution method. An eight-loci MLVA assay revealed six unique genotypes (G(K)13, G(K)27, G(K)35, G(K)43, G(K)44, and G(K)61). However, the 25-loci MLVA was more discriminatory, revealing the presence of ten genotypes instead of six. The additional genotypes resulted from the split of four subtypes: G(K)35 (b and c), G(K)43 (a and f), G(K)44 (d and e), and G(K)61 (i and j). All of the Turkish B. anthracis isolates were susceptible to ciprofloxacin, levofloxacin, tigecycline, linezolid, and vancomycin. One isolate was resistant to penicillin and to doxycycline. A total of 34 isolates were susceptible, 20 isolates were partially susceptible, and one isolate was resistant to erythromycin. None of the isolates exhibited susceptibility to cefotaxime. A total of 53 isolates were susceptible to gentamicin, and two were resistant. The genotypes G(K)35 (n=24), G(K)44 (n=13), and G(K)43 (n=10) were the most prevalent in 10, 6, and 5 regions, respectively, of the total 16 provinces. The B. anthracis isolates collected from these regions implied that the movement of B. anthracis is a result of the increased transportation of livestock and the resultant cross contamination.

Distribution and molecular evolution of bacillus anthracis genotypes in Namibia

The recent development of genetic markers for Bacillus anthracis has made it possible to monitor the spread and distribution of this pathogen during and between anthrax outbreaks. In Namibia, anthrax outbreaks occur annually in the Etosha National Park (ENP) and on private game and livestock farms. We genotyped 384 B. anthracis isolates collected between 1983-2010 to identify the possible epidemiological correlations of anthrax outbreaks within and outside the ENP and to analyze genetic relationships between isolates from domestic and wild animals. The isolates came from 20 animal species and from the environment and were genotyped using a 31-marker multi-locus-VNTR-analysis (MLVA) and, in part, by twelve single nucleotide polymorphism (SNP) markers and four single nucleotide repeat (SNR) markers. A total of 37 genotypes (GT) were identified by MLVA, belonging to four SNP-groups. All GTs belonged to the A-branch in the cluster- and SNP-analyses. Thirteen GTs were found only outside the ENP, 18 only within the ENP and 6 both inside and outside. Genetic distances between isolates increased with increasing time between isolations. However, genetic distance between isolates at the beginning and end of the study period was relatively small, indicating that while the majority of GTs were only found sporadically, three genetically close GTs, accounting for more than four fifths of all the ENP isolates, appeared dominant throughout the study period. Genetic distances among isolates were significantly greater for isolates from different host species, but this effect was small, suggesting that while species-specific ecological factors may affect exposure processes, transmission cycles in different host species are still highly interrelated. The MLVA data were further used to establish a model of the probable evolution of GTs within the endemic region of the ENP. SNR-analysis was helpful in correlating an isolate with its source but did not elucidate epidemiological relationships.

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

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

A rapid antimicrobial susceptibility test for Bacillus anthracis

An effective public health response to a deliberate release of Bacillus anthracis will require a rapid distribution of antimicrobial agents for postexposure prophylaxis and treatment. However, conventional antimicrobial susceptibility testing for B. anthracis requires a 16- to 20-h incubation period. To reduce this time, we have combined a modified broth microdilution (BMD) susceptibility testing method with real-time quantitative PCR (qPCR). The growth or inhibition of growth of B. anthracis cells incubated in 2-fold dilutions of ciprofloxacin (CIP) (0.015 to 16 microg/ml) or doxycycline (DOX) (0.06 to 64 microg/ml) was determined by comparing the fluorescence threshold cycle (C(T)) generated by target amplification from cells incubated with each drug concentration with the C(T) of the no-drug (positive growth) control. This DeltaC(T) readily differentiated susceptible and nonsusceptible strains. Among susceptible strains, the median DeltaC(T) values were > or = 7.51 cycles for CIP and > or = 7.08 cycles for DOX when drug concentrations were at or above the CLSI breakpoint for susceptibility. For CIP- and DOX-nonsusceptible strains, the DeltaC(T) was < 1.0 cycle at the breakpoint for susceptibility. When evaluated with 14 genetically and geographically diverse strains of B. anthracis, the rapid method provided the same susceptibility results as conventional methods but required less than 6 h, significantly decreasing the time required for the selection and distribution of appropriate medical countermeasures.

Discriminating virulence mechanisms among Bacillus anthracis strains by using a murine subcutaneous infection model

Bacillus anthracis strains harboring virulence plasmid pXO1 that encodes the toxin protein protective antigen (PA), lethal factor, and edema factor and virulence plasmid pXO2 that encodes capsule biosynthetic enzymes exhibit different levels of virulence in certain animal models. In the murine model of pulmonary infection, B. anthracis virulence was capsule dependent but toxin independent. We examined the role of toxins in subcutaneous (s.c.) infections using two different genetically complete (pXO1(+) pXO2(+)) strains of B. anthracis, strains Ames and UT500. Similar to findings for the pulmonary model, toxin was not required for infection by the Ames strain, because the 50% lethal dose (LD(50)) of a PA-deficient (PA(-)) Ames mutant was identical to that of the parent Ames strain. However, PA was required for efficient s.c. infection by the UT500 strain, because the s.c. LD(50) of a UT500 PA(-) mutant was 10,000-fold higher than the LD(50) of the parent UT500 strain. This difference between the Ames strain and the UT500 strain could not be attributed to differences in spore coat properties or the rate of germination, because s.c. inoculation with the capsulated bacillus forms also required toxin synthesis by the UT500 strain to cause lethal infection. The toxin-dependent phenotype of the UT500 strain was host phagocyte dependent, because eliminating Gr-1(+) phagocytes restored virulence to the UT500 PA(-) mutant. These experiments demonstrate that the dominant virulence factors used to establish infection by B. anthracis depend on the route of inoculation and the bacterial strain.

A simple method for the rapid removal of Bacillus anthracis spores from DNA preparations

This study establishes a filtration method for the safe removal of Bacillus anthracis spores which may contaminate DNA preparations. Centrifugal filtration with 0.1-microm filter units can be used following extraction of DNA from B. anthracis spores to render samples safe without compromising the sensitivity of diagnostic real-time PCR assays for B. anthracis.

Single nucleotide polymorphism typing of Bacillus anthracis from Sverdlovsk tissue

A small number of conserved canonical single nucleotide polymorphisms (canSNP) that define major phylogenetic branches for Bacillus anthracis were used to place a Sverdlovsk patient’s B. anthracis genotype into 1 of 12 subgroups. Reconstruction of the pagA gene also showed a unique SNP that defines a new lineage for B. anthracis.