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

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.

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.

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.

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.

Pathogenicity, population genetics and dissemination of Bacillus anthracis

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

HRM and SNaPshot as alternative forensic SNP genotyping methods

Single nucleotide polymorphisms (SNPs) have been widely used in forensics for prediction of identity, biogeographical ancestry (BGA) and externally visible characteristics (EVCs). Single base extension (SBE) assays, most notably SNaPshot® (Thermo Fisher Scientific), are commonly used for forensic SNP genotyping as they can be employed on standard instrumentation in forensic laboratories (e.g. capillary electrophoresis). High resolution melt (HRM) analysis is an alternative method and is a simple, fast, single tube assay for low throughput SNP typing. This study compares HRM and SNaPshot®. HRM produced reproducible and concordant genotypes at 500 pg, however, difficulties were encountered when genotyping SNPs with high GC content in flanking regions and differentiating variants of symmetrical SNPs. SNaPshot® was reproducible at 100 pg and is less dependent on SNP choice. HRM has a shorter processing time in comparison to SNaPshot®, avoids post PCR contamination risk and has potential as a screening tool for many forensic applications.

Genome Sequence of Historical Bacillus anthracis Strain Tyrol 4675 Isolated from a Bovine Anthrax Case in Austria

In 1988, an outbreak of anthrax occurred among cattle in the Austrian state of Tyrol. Since then, Austria has been declared anthrax-free. Here, we report the draft genome sequence of one of these last outbreak strains, Bacillus anthracis Tyrol 4675, isolated from a diseased cow.

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

Background

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

Results

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

Conclusions/Significance

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

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

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

Investigation of Anthrax Cases in North-East China, 2010-2014

We determined the genotypes of seven Bacillus anthracis strains that were recovered from nine anthrax outbreaks in North-East China from 2010 to 2014, and two approved vaccine strains that are currently in use in China. The causes of these cases were partly due to local farmers being unaware of the presence of anthrax, and butchers with open wounds having direct contact with anthrax-contaminated meat products. The genotype of five of the seven recovered strains was A.Br.001/002 sub-lineage, which was concordant with previously published research. The remaining two cases belongs to the A.Br.Ames sub-lineage. Both of these strains displayed an identical SNR pattern, which was the first time that this genotype was identified in North-East China. Strengthening education in remote villages of rural China is an important activity aimed at fostering attempts to prevent and control anthrax. The genotype of the vaccine strain Anthrax Spore Vaccine No.II was A.Br.008/009 and A.Br.001/002 for the vaccine strain Anthrax Spore Vaccine Non-capsulated. Further studies of their characteristics are clearly warranted.

Genetic Characterization of Bacillus anthracis 17 JB strain

BACKGROUND AND OBJECTIVES

Bacillus anthracis is one of the most homogenous bacteria ever described. Some level of diversity. Bacillus anthracis 17JB is a laboratory strain It is broadly used as a challenge strain in guinea pigs for potency test of anthrax vaccine.

MATERIAL AND METHODS

This work describes genetic characterization of B. anthracis 17 JB strain using the SNPs and MLVA genotyping.

RESULTS AND CONCLUSION

In SNPs typing, the originally French 17JB strain represented the A.Br. 008/009 subgroup. In Levy's genotyping method, 843, 451 and 864 bp long fragments were identified at AA03, AJ03 and AA07 loci, respectively. In the vaccine manufacturer perspective these findings are much valuable on their own account, but similar research is required to extend molecular knowledge of B. anthracis epidemiology in Persia.

Draft Genomes of Three Strains Representative of the Bacillus anthracis Diversity Found in France

We report here the draft genomes of three Bacillus anthracis strains isolated in France: 08-8_20 (A.Br.001/002), 99-100 (A.Br.011/009), and 00-82 (B.Br CNEVA). The total lengths of assemblies are 5,440,708 bp, 5,446,472 bp, and 5,436,014 bp for 08-8_20, 99-100, and 00-82, respectively.

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

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

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

Background

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

Results

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

Conclusions

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

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

Background

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

Results

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

Conclusions

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

Rapid identification of bio-molecules applied for detection of biosecurity agents using rolling circle amplification

Detection and identification of pathogens in environmental samples for biosecurity applications are challenging due to the strict requirements on specificity, sensitivity and time. We have developed a concept for quick, specific and sensitive pathogen identification in environmental samples. Target identification is realized by padlock- and proximity probing, and reacted probes are amplified by RCA (rolling-circle amplification). The individual RCA products are labeled by fluorescence and enumerated by an instrument, developed for sensitive and rapid digital analysis. The concept is demonstrated by identification of simili biowarfare agents for bacteria (Escherichia coli and Pantoea agglomerans) and spores (Bacillus atrophaeus) released in field.