1. Introduction
Bacillus (B.) anthracis, the causative agent of anthrax, belongs to the Bacillus cereus group together with at least six other species (B. cereus, B. thuringiensis, B. mycoides, B. pseudomycoides, B. weihenstephanensis and B. cytotoxicus). They are Gram-positive, spore-forming, aerobic, facultatively anaerobic, rod-shaped bacteria [1]. In contrast to other Bacillus species, it is non-motile, non-haemolytic on sheep blood agar, grows at a temperature of 37 °C, and forms large colonies with irregularly tapered outgrowths. In culture, it tends to form long chains of bacteria. Within a host, however, it appears either as single organisms or as chains consisting of two or three bacilli [2]. Anthrax is primarily a disease of herbivores, although it has also been reported to occur in omnivores, carnivores, and other vertebrates. Herbivorous animals, such as cattle, sheep, and goats, ingest spores in the soil and become infected. Humans can become infected through direct contact with infected animals or tissue from infected animals or through direct exposure to B. anthracis [3]. Based on the initial route of infection, three primary forms of clinical infection in humans are described: inhalational anthrax, gastrointestinal anthrax, and cutaneous anthrax [4] A new form of soft tissue infection associated with injection drug use, termed injectional anthrax, has also been described [5].
B. anthracis is also considered one of the most important pathogens of bioterrorism. Due to the often-fatal outcome of human cases, the rapid administration of clearly effective antimicrobials is critical, whether for prophylaxis, following suspected exposure, or for the treatment of clinical cases. B. anthracis is a bacterium that is sensitive to most antibiotics. However, early treatment is critical to eliminate the bacterium before it releases toxins into the bloodstream [6,7].
When conditions are not favourable for the growth and multiplication of vegetative forms of B. anthracis, they begin to form spores. Spores are highly resistant to biological extremes such as heat, cold, pH, desiccation, chemicals, irradiation, and other unfavourable conditions. The organism can remain in the spore stage for long periods of time, waiting for the moment when conditions are favourable for germination and multiplication [8]. Jacotot and Virat [9] and Wilson and Russell [10] succeeded in reviving spores after more than 60 years, while de Vos [11] recovered anthrax spores from bones after 200 ± 50 years. From the earliest historical records until the development of an effective veterinary vaccine in the middle of the 20th century and the subsequent advent of antibiotics, the disease was a major cause of uncontrolled mortality in cattle, sheep, goats, horses, and pigs worldwide [8]. Anthrax is found all over the world, on every continent except Antarctica. There are endemic areas with more frequent outbreaks, while sporadic outbreaks occur in other areas [8,12,13]. B. anthracis harbours two plasmids, pXO1 (182 kb) and pXO2 (95 kb). The most important virulence factors, the tripartite anthrax toxin genes cya, lef, and pag, are located on the plasmid pXO1 [14,15,16], while the poly-γ-D-glutamic acid capsule genes capABCDE are localised on the plasmid pXO2 [17,18]. Previously, the presence of these two virulence plasmids was considered the main distinguishing feature between B. anthracis and closely related species, particularly Bacillus cereus and Bacillus mycoides [19]. In the meantime, the virulence plasmids pXO1 and pXO2 have also been found in other isolates of the B. cereus group [20,21,22]. B. anthracis is an evolutionarily young pathogen characterised by a lack of molecular variation. Isolates from all over the world show an extremely high genetic homology [23,24]. It is assumed that the dormant state of its spores probably greatly reduces the rate of evolutionary change of B. anthracis and contributes to its extreme homogeneity [25]. Due to this lack of diversity, only modern molecular characterisation techniques have proven effective in distinguishing strains within this highly clonal species. Standard genotyping methods to determine the phylogenetic relationships between B. anthracis isolates include canonical single-nucleotide polymorphisms (canSNPs) and multilocus variable-number tandem repeat analysis (MLVA) [26,27,28]. Recently, core-genome-based multilocus sequence typing (cgMLST) has been increasingly used as a comparative tool using WGS data [29]. Phylogenetic analysis has defined three major clades, with C being basal to clades A and B. Members of clade A are most commonly observed worldwide (~90%), while members of clades B (~10%) and C (<1%) are much less frequent [27]. Clade A is divided into seven subclades: Ancient A, Vollum, V770, Sterne, Ames, Australia 94, and Trans Eurasia [30].
Historical data on anthrax in humans and animals in Croatia can be found in the Supplementary Materials.
Anthrax has occurred only sporadically in Croatia, as immunoprophylaxis was systematically carried out in enzootic areas where anthrax outbreaks had occurred [31]. Since 1981, a total of 59 cases of anthrax in animals have been confirmed at the Croatian Veterinary Institute in Croatia. The course of the disease is characterised by the fact that in some periods the disease does not occur. Within the period of more than four decades, no disease was detected in 20 years (1981–1983, 1992–1996, 1998–2000, 2008–2010, 2013, 2015–2019 and 2021) [31,32]. Only a few major outbreaks have been described in recent decades: 2001/2002 (11 animals–5 isolates) [33], 2006/2007 (11 animals–6 isolates) [34], and the largest outbreak in 2022 (323 animals–25 isolates) [35]. Sporadic cases occur mainly in the southern part of Croatia, where sporadic cases have been recorded over several years (4 isolates tested in this study) [31].
The aim of this study was to describe the reliable and rapid diagnosis of epizootic 2022 and to determine the in vitro antimicrobial resistance (AMR) of 40 isolates of B. anthracis isolated in Croatia in the recent two decades to 18 different antibiotics. In addition, we performed whole-genome sequencing (WGS) on all isolates to determine AMR and virulence genes and used cgMLST as a tool for isolate comparison and phylogenetic analysis.
2. Results
2.1. Isolates
A total of 40 isolates were included in this study. In total, 25 isolates were from the 2022 outbreak, while 15 isolates were obtained from laboratory strain collections from outbreaks in the previous years (2001/2002-5 isolates, 2006/2007-6 isolates) and from areas with sporadic anthrax occurrence (ST) (4 isolates). As shown in Table 1, the presence of the virulent plasmids pXO1 and pXO2 was confirmed in all 40 isolates. Strain Sterne 34F2 lacked the pXO2 plasmid.
2.2. Antimicrobial Susceptibility
The results—expressed as MIC at which 50% and 90% of the tested isolates were inhibited (MIC50, MIC90)—the range (mg/L), the breakpoint, and the percentage of susceptible isolates are shown in Table 2.
All isolates showed susceptibility to quinupristin/dalfopristin, vancomycin, ampicillin, rifampin, levofloxacin, penicillin, oxacillin + 2%NaCl, clindamycin, linezolid, tetracycline, gentamicin, gatifloxacin, and ciprofloxacin. However, 40% of the isolates were intermediate susceptible to ceftriaxone, while 5% were intermediate susceptible to erythromycin. All isolates tested were resistant to daptomycin and trimethoprim/sulfamethoxazole.
2.3. WGS and Bioinformatic Analysis
The genes for antimicrobial resistance and virulence factors, protein function, percentage of overlap, and identities found in 40 B. anthracis strains isolated in Croatia are listed in Table 3.
We identified a total of four AMR genes in all B. anthracis isolates from the 2001/2002 outbreak and isolates from the southern part of Croatia (bla, blaII, fosB, satA), while isolates from the 2006/2007 and 2022 outbreaks lacked the fosB gene.
As far as the virulence genes are concerned, all genes listed in Table 3 were found in all tested isolates without exception.
2.4. Core-Genome-Based Multilocus Sequence Typing (cgMLST)
cgMLST analysis was performed for 3803 core loci and 1263 accessory genomic loci. Forty-seven and forty-five loci were missing in the 2001/2002 outbreak samples and in the samples from southern Croatia, respectively. Complex type 219 was defined for the samples from the 2001/2002 outbreak and 222 for the samples from the southern part of the country.
In total, 145 and 148 loci were missing in the samples from the 2006/2007 and 2022 outbreaks, respectively. Complex type 220 was defined for samples from 2006/2007 and 221 for samples from the 2022 outbreak.
The phylogenetic analysis of the strains using a neighbour-joining tree cluster analysis based on the cgMLST data (Figure 1) shows that the B. anthracis isolates from the outbreaks in Croatia were divided into two clades. Isolates from the 2001/2002 outbreak and isolates from southern part of country were categorised in clade A, subclade Trans Eurasia, while isolates from the 2006/2007 and 2022 outbreaks were categorised in clade B, subclade B2.
3. Discussion
In Europe, anthrax is considered a sporadically occurring disease that almost exclusively affects animals at pastures [36]. In Croatia, anthrax is also rare in animals and humans, but there are sporadic outbreaks, especially in the anthrax districts, which are mainly located in the Lonjsko Polje region and in the southern parts of Croatia [31]. The situation is similar in the neighbouring countries where the occurrence of anthrax has been described: Italy [37], Slovenia [38], Hungary [39], and Bosnia and Herzegovina [40].
Antimicrobial susceptibility testing plays a key role in the treatment of anthrax infections. Rapid diagnosis is also crucial for post-exposure prophylaxis (PEP) or the early treatment of anthrax to eliminate the bacterium before it releases toxins into the bloodstream [6]. Beta-lactam antibiotics such as penicillin are recommended by the WHO and CDC as prophylactic therapy for anthrax, while tetracyclines and fluoroquinolones are recommended by the CDC for PEP and for the treatment of anthrax [8,41].
As shown in this study, a reliable diagnosis of anthrax can be made within 24 h based on colony morphology, haemolysis on blood agar, and the detection of the plasmids pXO1 and pXO2. In this study, a total of 40 B. anthracis strains isolated in Croatia were analysed for their susceptibility to antimicrobial agents. All tested isolates were susceptible to penicillin, tetracycline, and ciprofloxacin, with MIC90 of ≤0.06, ≤2, and ≤0.5, respectively. The first penicillin resistance in a clinical B. anthracis isolate was described in 1976 in a fatal case of anthrax in Northampton, England [42]. Although penicillin resistance of B. anthracis is thought to be at least 10% [43,44], the largest study which included 110 strains from outbreaks in Italy reported a 100% susceptibility to penicillin [45]. In this study, we found two beta-lactamase genes, bla and blaII. It has already been described that B. anthracis possesses these two beta-lactamase genes, bla1 (penicillinase) and bla2 (cephalosporinase), but that their transcription is very low and gene expression is not sufficient to confer resistance to beta-lactam antibiotics [46,47]. Although these genes are present, it has been hypothesised that they are not expressed due to a truncation of the positive regulatory gene, plcR [48]. The transcriptional activator PlcR is responsible for the expression of extracellular proteins such as haemolysins, enterotoxins, and proteases in B. cereus and B. thuringiensis [49,50]. The sequence of the plcR gene of B. anthracis strains indicates that the expression of the gene results in the synthesis of a truncated PlcR protein, in which codon 214 for glutamic acid (GAA or GAG) in B. cereus and B. thuringiensis is converted to a termination codon (TAA) in B. anthracis by a single nucleotide substitution (G→T) [51]. Although the plcR homologue of B. anthracis is truncated, there are 56 putative plcR binding motifs on the chromosome and 2 on pXO2. Extracellular protein genes include phospholipases, enterotoxins, and haemolysins, and the plcR mutation has been shown to be responsible for a drastic reduction in lecithinase, protease, and haemolysin production, giving B. anthracis a different secretome to B. cereus [48,52]. It has been shown that the transcription of beta-lactamases is controlled by an extracytoplasmic sigma factor (ECF) SigP and its anti-sigma factor RsiP [53]. Beta-lactamase expression in B. cereus and B. thuringiensis is inducible by beta-lactam antibiotics, leading to the proteolysis of RsiP, release of SigP, and transcription of the corresponding genes [54]. This beta-lactam feedback system is not functional in prototypic B. anthracis isolates. However, in a penicillin-resistant B. anthracis strains, the rsiP gene was nonsense mutated and the strain constitutively produced beta-lactamases [53].
All B. anthracis isolates in our study were also found to be susceptible to ampicillin and oxacillin +2%NaCl, with MIC90 less or equal to the minimum concentration tested (≤0.12; ≤0.25). This is consistent with previous publications reporting MICs for the other penicillin that closely matched the penicillin MIC [44].
Regarding other fluoroquinolones used in this study, all isolates were found to be susceptible to levofloxacin and gatifloxacin, with MIC90 values ≤0.25 and ≤1, respectively. Ciprofloxacin and tetracycline resistance of B. anthracis strains is very rare and is only described in a few publications [44]. The susceptibility data presented in this study support the current recommendations for the use of penicillin, tetracycline, and ciprofloxacin in PEP and anthrax threats.
Forty percent of the isolates tested in this study were intermediate susceptible to ceftriaxone, while sixty percent were sensitive. This is in partial agreement with Manzulli et al. [45] and Cavallo et al. [43], who categorised all isolates as intermediate susceptible to ceftriaxone. In contrast to these studies, Jones et al. [55] reported all isolates as sensitive to ceftriaxone. Due to this variable susceptibility, non-penicillin beta-lactams are not recommended for the treatment of anthrax.
All B. anthracis isolated in this study were also susceptible to linezolid, rifampin, gentamycin, vankomycin, quinupristin/dalfopristin, and clindamycin. This is consistent with previous published studies [45,56,57].
Linezolid is used for the treatment of a variety Gram-positive infections. It is a bacteriostatic oxazolidinone with a unique mode of action that inhibits bacterial protein synthesis by preventing the formation of the 70S ribosomal complex [58]. For this reason, cross-resistance is unlikely. It is therefore a good alternative for resistant B. anthracis strains. Although all isolates were susceptible to rifampin and clindamycin, they are not considered a good alternative for anthrax treatment as resistance can easily be generated in vitro [59,60]. Nevertheless, rifampin, clyndamicin, and vancomycin together with fluoroquinolones are considered good candidates for the treatment of inhalational anthrax and meningitis [61]. Athamna et al. [62] found that quinupristin/dalfopristin killed B. anthracis the fastest of all antimicrobial agents tested.
Resistance to trimethoprim/sulfamethoxazole was detected in all tested isolates. This result was to be expected as intrinsic resistance to this antimicrobial agent has already been described [63]. All isolates in this study were also resistant to daptomycin with MIC50 and MIC90 of 4 mg/L. This result is consistent with previous publications [56,64].
Another AMR gene found in all isolates was satA, while the fosB gene was only found in the B. anthracis isolates from the 2001/2002 outbreak and isolates from sporadic outbreaks isolated in southern Croatia. SatA belongs to the streptothricin acetyltransferase gene family and confers resistance to nucleoside antibiotics, which are produced by actinomycetes in soil. It has already been described in B. anthracis and B. subtilis [65]. The fosB gene is a chromosomally derived thiol transferase and confers resistance to fosfomycin [48]. In their study, Bruce et al. [66] assumed that the detection of the fosB gene can be a clear prediction for the affiliation of the isolates to clade C. Clade C is largely defined by the absence of the fosB gene, which is found across all other populations. However, this is in partial agreement with the results of our study, in which all isolates that possessed the fosB gene belonged only to clade A, a subclade of Trans Eurasia, but not to clade B. Similar results to ours were also reported by Chiaverini et al. [67].
Pathogenicity-related plasmids, pxO1 and pXO2, were found in all 40 B. anthracis isolates in this study. In addition, all isolates possessed the atxA gene, which is localised on the pathogenicity island pXO1 and is responsible for the expression of the three toxin genes [68,69].
Other virulence genes found in all B. anthracis strains in this study were inhA1, inhA2, entFM, nheABC operon, plcR, sph, plcA, and alo.
InhA1 cleaves multiple host proteins, including collagen, fibrinogen, plasminogen activator inhibitor, and prothrombin [70,71]. It also downregulates BslA, a B. anthracis cell surface protein associated with the adhesion of vegetative bacterial cells to the host endothelium [72]. The B. anthracis strains in this study also contained a gene encoding the metalloprotease InhA2, although it is not known whether this protease is expressed and secreted. This gene is an orthologue of InhA1 described above (68% amino acid identity) [73]. The metalloprotease (inhA) and the phospholipases sph genes are important virulence factors as they enable the bacilli to evade innate and adaptive immune responses during the infective phases [74,75].
The chromosomally encoded nheABC operon was previously described in all B. anthracis isolates [76]. It is responsible for the synthesis of termolabile non-haemolytic enterotoxin, which, together with enterotoxin FM and haemolysin HII, is involved in the diarrhoeic syndrome [77,78,79].
In this study, we used the cgMLST method to genotype B. anthracis isolates from the recent two decades in Croatia. The development of a cgMLST scheme for the entire B. cereus group is timely and ensures high-resolution strain typing of B. cereus group isolates in an open and fully publicly accessible online user interface for the future [80]. The first application of a whole-genome allele typing system was presented by Abdel-Glil et al. [29]. We used cgMLST-searching 3803 loci and accessory cgMLST-searching 1263 loci. After phylogenetic analysis based on the cgMLST data neighbour-joining tree was constructed (Figure 1). The colours represent the clade and subclade, which are also indicated by name, while the node labels correspond to the names of the strains and the country of origin. As shown in Figure 1, the Croatian B. anthracis isolates were categorised into two main clades, depending on the outbreaks: A and B. isolates from the 2006/2007 and 2022 outbreaks were categorised in clade B, subclade B2, while the isolates from the 2001/2002 outbreak and isolates from sporadic cases from the southern part of Croatia (ST2011) were categorised in clade A, subclade Trans Eurasia (TEA). The existence of two different genotypes of B. anthracis isolates in Croatia was also indicated by the finding that the fosB gene was discovered in all isolates from the 2001/2002 outbreak and in isolates from sporadic outbreaks from southern Croatia, while it was not present in isolates from the 2006/2007 and 2022 outbreaks.
It has already been reported that strains belonging to clade A are responsible for the majority of anthrax cases reported worldwide, while the TEA subclade is one of the most common. These isolates are widespread in Europe, the Middle East, Russia, and the Chinese province of Xinjang [27,81]. In Europe, it is the most dominant subclade in Albania, Bulgaria, Hungary, and Italy [27,82]. The MLVA analysis of 234 B. anthracis strains isolated in Italy between 1972 and 2018 revealed the distribution of 55 B. anthracis genotypes, while the CanSNPs analysis categorised 53 of them into a TEA subclade [83]. As can be seen in Figure 1, the Croatian isolates assigned to the TEA subclade showed the greatest homology with the Italian isolate (3IZSLT).
Clade B occurs mainly in South Africa and in the Central Europe, especially in France and Poland. The South African strains mainly belong to the B1 subclade (B.Br.Kruger and B.Br.001/002 groups) [84], while the European strains belong to the B2 subclade (B.Br.CNEVA) [85,86]. Examples of these subclades are rarely found outside these regions [27]. The B2 subclade is widespread in Europe in south-west France, north-east Poland, northern Italy, Switzerland, Germany, and Slovenia, and it is the most dominant subclade in France, Poland, and northern Italy [87,88,89]. The Croatian B. anthracis isolates belonging to subclade B2 showed the most similarities with the Austrian isolates (Tyrol 3520) and, as can be seen in Figure 1, the entire B2 subclade branches from isolate ‘RA3 France’. Similar results were also obtained by Chiaverini et al. in Italy [67].
The TEA subclade has been found to occur in countries where the B2 subclade is dominant and is associated with sporadic cases [90], which is in complete agreement with our study.
These data show that anthrax outbreaks in Croatia are limited to a relatively small area of the country and are caused by two genetically very different B. anthracis genotypes belonging to two different, phylogenetically unrelated clades, indicating a different origin of the isolates.
4. Materials and Methods
4.1. Epidemiological Investigation
At the end of June 2022, cattle breeders noticed an increase in cattle deaths at the Osekovo, Donja Gračenica, and Veliko Svinjičko pastures in the Lonjsko Polje Nature Park area in Croatia (Figure 2).
Considering the fact that it was very close to the border of the anthrax district, the Veterinary and Food Safety Authority authorised the dissection of the two carcasses. The first dissection of non-autolytic cattle was carried out on 15 July 2022 at the rendering plant under controlled conditions and the samples (blood sample, spleen) were delivered to the Laboratory of General Bacteriology and Mycology of the Croatian Veterinary Institute on the same day, as it was suspected that the animals had died of anthrax. Since the beginning of the outbreak of the disease until the end of 2022, 323 animals (246 cattle, 77 horses) have died in the Lonjsko Polje Nature Park area. This number is probably even higher, as we only have data on the number of dead animals brought to the rendering centre.
By the end of 2022, a total of 101 necropsies were performed, mainly on animals that died on the same day or at least the day before. The autopsy of the corpses revealed that there was no rigour mortis and that unclogged blood leaked from all natural orifices. Oedema was found in the subcutaneous tissue. Splenomegaly, haemorrhages in the epicardium and endocardium, and pulmonary oedema with haemorrhages in the parenchyma were noted. The lymph nodes were enlarged and oedematous, and haemorrhages were observed in the parenchyma. Ulcerations were found on the mucosa of the abomasum and small intestine. Petechial haemorrhages were found in the serosa of the intestine.
Anthrax in humans in the Republic of Croatia is characterised by the occurrence of isolated cases. According to data from the Croatian Agency for Agriculture and Food, 16 cases of the disease were recorded in the Republic of Croatia in the period from 2005 to 2016 [91], while, according to data from the Croatian Institute of Public Health [92], two cases of the disease were recorded in humans in the period from 2017 to 2021 (in 2018). Figure 3.
The anthrax epidemic in humans in the Lonjsko Polje Nature Park began on 13 July 2022 and ended on 19 August 2022. Seventeen people fell ill with the cutaneous form of anthrax (Figure 4). Among the infected were 3 children up to the age of 8, 4 aged 15–18, and 10 adults. All had contact with sick cattle and were treated with antibiotics [93].
4.2. Laboratory Diagnosis
Fresh samples (mainly blood, spleen, and ears, with parts of the musculature) were sent to the Laboratory for General Bacteriology and Mycology of the Croatian Veterinary Institute. The spleen samples and ears with muscles were delivered in sterile plastic containers with screw caps, with each container being wrapped in two plastic bags to prevent leakage. The blood samples were delivered in syringes, and each syringe was placed in a plastic container with a screw cap.
Thin smears were made from all samples supplied and a Giemsa stain was used. The samples were inoculated onto blood agar (OXOID blood agar base No. 2 + 5% defibrinated sheep blood). After incubation under aerobic conditions at a temperature of 37 °C for 18–24 h, the plates were examined to determine bacterial growth, morphology, and colour of the bacterial colonies. Bacterial colonies morphologically consistent with B. anthracis were resuspended in distilled water and heated to 95 °C for 20 min. After cooling and centrifugation, the supernatant was used for a PCR reaction to confirm the presence of plasmids pXO1 and pXO2, using the previously described method [94,95].
The presence of plasmids pX01 and pX02 was confirmed in all B. anthracis strains from this epizootic as well as in the other strains from previous outbreaks in Croatia.
4.3. Isolates Used in This Study
We analysed a total of 40 B. anthracis isolates: 5 isolates from the 2001/2002 outbreak, 6 isolates from the 2006/2007 outbreak, 4 isolates from a southern part of Croatia, and 25 isolates from the 2022 outbreak. The molecular part of this study included the B. anthracis Sterne 34F2 (pXO1+, pXO2−) strain.
4.4. Antimicrobial Susceptibility Testing
Antimicrobial susceptibility testing (AST) was performed using the broth microdilution method according to Clinical and Laboratory Standards Institute (CLSI) guidelines [96].
Briefly, isolates were recovered from freezer stocks (tryptic soya broth with 20% glycerol, −80 °C) and incubated overnight on blood agar supplemented with 5% sheep blood in an aerobic atmosphere. AST was performed on GPN3F microtiter plates (Sensititer, Trek Diagnostic Systems Ltd. East Grinstead, West Sussex, UK). Cation-adjusted Mueller–Hinton broth was used to prepare the 0.5 McFarland solution. The microtiter plates were incubated in sealed 96-well microtiter plates at 35 ± 1 °C for 18 ± 2 h.
The susceptibility to erythromycin (ERY; 0.25–4 mg/L), quinupristin/dalfopristin (SYN; 0.12–4 mg/L), vancomycin (VAN; 1–128 mg/L), ampicillin (0.12–16 mg/L), rifampin (0.5–4 mg/L), levofloxacin (0.25–8 mg/L), penicillin (0.06–8 mg/L), trimethoprim/sulfamethoxazole (0.5/9.5–4/76 mg/L), oxacillin + 2%NaCl (0.25–8), ceftriaxone (8–64 mg/L), clindamycin (0.12–2 mg/L), daptomycin (0.25–8 mg/L), linezolid (0.5–8 mg/L), tetracycline (TET 2–16 mg/L), gentamicin (2–16 mg/L), gatifloxacin (1–8 mg/L), ciprofloxacin (0.5–2 mg/L), and streptomycin (1000 mg/L) was tested.
The lowest concentration of an antibiotic that prevented bacterial growth was taken as the minimum inhibitory concentration (MIC).
For penicillin, ampicillin, tetracycline, ciprofloxacin, and levofloxacin, the reported CLSI breakpoints for B. anthracis were used; for ceftriaxone, vancomycin, gentamicin, erythromicin, clindamycin, trimethoprim/sulfamethoxazole, and rifampin, Bacillus spp. breakpoints were used [96]. Interpretation criteria for Staphylococcus sp. were used for quinupristin/dalfopristin, oxacillin + 2%NaCl, linezolid, and gatifloxacin [97].
Streptomycin was not considered further, as there was only one concentration of this antimicrobial agent on the microplate, which was not important for further consideration.
The reference strains Staphylococcus aureus ATCC 29213 and Escherichia coli ATCC 25922 were used for quality control.
4.5. Whole-Genome Sequencing (WGS) and De Novo Assembly
The B. anthracis strains were sequenced using Illumina MiSeq (Illumina, San Diego, CA, USA). DNA was extracted using the Nucleospin Microbial DNA Kit (Macherey-Nagel, Düren, Germany). Paired-end sequencing libraries were prepared using the DNA LP Tagmentation kit (Illumina, USA), and the sequencing was performed using the MiSeq Reagent Kit v2 (300 cycles) (Illumina, USA). The raw sequences were trimmed with Trimmomatic v. 0.38.1. (options SLIDINGWINDOW 4, quality 20; MINLEN 100; LEADING 15; TRAILING 15; HEADCROP 15; AVGQUAL 25) [98]. Genome assembly was performed with Shovill v.1.1.0 (DEPTH 100; GENOME SIZE 5,5 Mb; POSTASSEMBLY CORRECTION 500; MIN COVERAGE 5) [99] for paired-end Illumina data. The assembly statistics were obtained using the program Quast v.5.0.2 [100].
The genome assemblies were deposited in GenBank under the BioProject accession number PRJNA1129309 (SAMN42150189—2001/2002 outbreak; SAMN42150190—2006/2007 outbreak; SAMN42150191—isolates from southern Croatia; SAMN42150192—2022 outbreak).
4.6. Bioinformatic Analysis of Antimicrobial Resistance and Virulence Genes
Virulence was tested with the VFDB database [101]. NCBI AMRFinderPlus v.3.11.26 was used for the identification of the resistance and virulence genes using BLASTX search [102]. MobileElementFinder v.1.0.5 was used for the identification of Mobile Genetic Elements (MGE) [103].
4.7. Core-Genome-Based Multilocus Sequence Typing (cgMLST)
cgMLST typing was performed using Ridom SeqSphere+ v.10.0.0 [104]. Both cgMLST-searching 3803 loci and accessory cgMLST-searching 1263 loci were used. For all analysed genomes, the cgMLST allelic profiles were pairwise-compared (missing values were ignored), and the resultant calculated pairwise distances were used to generate a neighbour-joining tree using SeqSphere+.
5. Conclusions
Anthrax in animals and humans in Croatia is characterised by sporadic outbreaks, which are limited to anthrax districts located mainly in the Lonjsko Polje region and in the southern parts of Croatia.
Rapid and reliable diagnosis and administration of antimicrobials are essential for PEP or effective treatment of anthrax. Molecular diagnosis, based on the presence of pXO1 and pXO2 plasmids and the morphological appearance of colonies, enables reliable anthrax diagnosis within 24 h in material that is not heavily contaminated.
Isolates of B. anthracis isolated in Croatia in the recent two decades were susceptible to all antimicrobials recommended for PEP or anthrax therapy (penicillin, ampicillin, tetracycline, and ciprofloxacin). Susceptibility was observed for all other antimicrobials tested, which represent an alternative for primary therapy. Three or four AMR genes were found in our isolates, depending on the area and year of the outbreak, of which two beta-lactamase genes were not expressed. 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 the outbreaks in Croatia are divided into two clades. The first is clade A, subclade Trans Eurasia, which is predominant in the Mediterranean region and the eastern part of Europe and was responsible for the 2001/2002 outbreak and the outbreaks in the southern part of the country. Isolates from the 2006/2007 and 2022 outbreaks are categorised in clade B, subclade B2, which is dominant in Central European countries. These results indicate that two different B. anthracis clades are circulating in a very small area, which is probably due to the geographical location of Croatia.
Conceptualization, G.K., S.D., S.Š. and B.H.; methodology, G.K., S.D., I.R. and Ž.M.; software, S.D. and I.R.; validation, G.K., S.D., I.R., A.B. and B.H.; formal analysis, G.K., S.D., I.R., A.B., Ž.M., M.B. and L.C.; investigation, G.K., S.D., A.B. and B.H.; data curation, G.K., S.D., I.R. and L.C.; writing—original draft preparation, G.K., S.D. and Ž.M.; writing—review and editing, G.K., S.D., I.R. and A.B.; visualization, G.K. and S.D.; supervision, G.K. All authors have read and agreed to the published version of the manuscript.
Not applicable.
Not applicable.
The data are contained in the article and
The authors declare no conflict of interest.
Footnotes
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Enlarge this image.Figure 1. Phylogenetic analysis of the strains using a neighbour-joining tree based on the cgMLST data. The colours represent the clade and subclade, which are also indicated by the name, while the node labels correspond to the names of the strains and the country of origin. The isolates used in this study are indicated by arrows (red—outbreak 2006/2007; yellow—outbreak 2022; blue—isolates from southern part of Croatia; black—outbreak 2001/2002).
Enlarge this image.Figure 2. Geographical location of the Lonjsko Polje Nature Park. The locations of the 2022 outbreak are marked with blue dots. The red dots represent the anthrax outbreaks of 2001 and 2007. The triangle represents the southern part of Croatia where anthrax occurs sporadically.
PCR confirmation of pXO1 and pXO2 plasmids in 40 strains tested in Croatia.
| Isolates | Number of Isolates | pXO1 | pXO2 |
|---|---|---|---|
| 2001/2002 | 5 | + | + |
| 2006/2007 | 6 | + | + |
| 2022 | 25 | + | + |
| ST | 4 | + | + |
| Sterne 34F2 | 1 | + | - |
Antimicrobial susceptibility of 40 B. anthracis isolates from Croatia.
| Antibiotic | MIC50 | MIC90 | Range (mg/L) | Breakpoint | % S/I/R a | |
|---|---|---|---|---|---|---|
| S (≤) | R (≥) | |||||
| Erythromycin ** | 0.5 | 0.5 | 0.25–4 | 0.5 | 8 | 95/5/0 |
| Quinupristin/dalfopristin *** | 1 | 1 | 0.12–4 | 1 | 4 | 100/0/0 |
| Vancomycin ** | 2 | 2 | 1–128 | 4 | 100/0/0 | |
| Ampicillin * | ≤0.12 | ≤0.12 | 0.12–16 | 0.12 | 0.25 | 100/0/0 |
| Rifampin ** | ≤0.5 | ≤0.5 | 0.5–4 | 1 | 4 | 100/0/0 |
| Levofloxacin * | ≤0.25 | ≤0.25 | 0.25–0.5 | 0.25 | 100/0/0 | |
| Penicillin * | ≤0.06 | ≤0.06 | 0.06–8 | 0.5 | 1 | 100/0/0 |
| Trimethoprim/sulfamethoxazole ** | >4/76 | >4/76 | 0.5/9.5–4/76 | 2/38 | 4/76 | 0/0/100 |
| Oxacillin + 2%NaCl *** | ≤0.25 | ≤0.25 | 0.5–8 | 0.5 | 1 | 100/0/0 |
| Ceftriaxone ** | ≤8 | 16 | 8–64 | 8 | 64 | 60/40/0 |
| Clindamycin ** | 0.25 | 0.25 | 0.12–2 | 0.5 | 4 | 100/0/0 |
| Daptomycin *** | 4 | 4 | 0.25–8 | 1 | 0/0/100 | |
| Linezolid *** | 1 | 1 | 0.5–8 | 4 | 8 | 100/0/0 |
| Tetracycline * | ≤2 | ≤2 | 2–16 | 1 | 100/0/0 | |
| Gentamicin ** | ≤2 | ≤2 | 2–16 | 4 | 16 | 100/0/0 |
| Gatifloxacin *** | ≤1 | ≤1 | 1–8 | 0.5 | 2 | 100/0/0 |
| Ciprofloxacin * | ≤0.5 | ≤0.5 | 0.5–2 | 0.25 | 100/0/0 | |
| Streptomycin | ≤1000 | ≤1000 | 1000 | |||
a S/I/R, susceptible/intermediate/resistant. * CLSI M45 B. anthracis interpretation. ** CLSI M45 Bacillus sp. interpretation. *** CLSI M100 Staphylococcus spp. interpretation.
AMR and virulence factors genes of 40 B. anthracis isolates from Croatia.
| Gene | Name | Type | Aligned Overlap 100%/Identity ≥ 97% | |||
|---|---|---|---|---|---|---|
| 2001/2002 | 2006/2007 | ST | 2022 | |||
| bla | Class A beta-lactamase Bla1 | AMR | + | + | + | + |
| blaII | BcII family subclass B1 metallo-beta-lactamase | AMR | + | + | + | + |
| fosB | fosB/fosD family fosfomycin resistance bacillithiol transferase | AMR | + | - | + | - |
| satA | Streptothricin N- acetyltransferase SatA | AMR | + | + | + | + |
| pagA | anthrax toxin protective antigen | VIR | + | + | + | + |
| lef | anthrax toxin lethal factor | VIR | + | + | + | + |
| cya | anthrax toxin edema factor | VIR | + | + | + | + |
| capA | capsular polyglutamate synthetase CapA | VIR | + | + | + | + |
| capB | capsule biosynthesis protein CapB | VIR | + | + | + | + |
| capC | capsular polyglutamate amide ligase/translocase CapC | VIR | + | + | + | + |
| capD | capsule biosynthesis gamma-glutamyltransferase CapD | VIR | + | + | + | + |
| capE | capsule biosynthesis protein CapE | VIR | + | + | + | + |
| atxA | anthrax toxin expression trans-acting transcriptional regulator AtxA | VIR | + | + | + | + |
| inhA1 | M6 family metalloprotease immune inhibitor InhA1 | VIR | + | + | + | + |
| inhA2 | M6 family metalloprotease immune inhibitor InhA2 | VIR | + | + | + | + |
| entFM | enterotoxin EntFM | VIR | + | + | + | + |
| nheA | non-hemolytic enterotoxin NHE subunit A | VIR | + | + | + | + |
| nheB | non-hemolytic enterotoxin NHE subunit B | VIR | + | + | + | + |
| nheC | non-hemolytic enterotoxin NHE subunit C | VIR | + | + | + | + |
| plcR | transcriptional regulator PlcR | VIR | + | + | + | + |
| sph | sphingomyelinase C | VIR | + | + | + | + |
| hlyIIR | hemolysin II regulator HlyIIR | VIR | + | + | + | + |
| cerA | phospholipase CerA | VIR | + | + | + | + |
| plcA | phosphatidylinositol diacylglycerol-lyase | VIR | + | + | + | + |
| alo | anthrolysin O/cereolysin O family cholesterol-dependent cytolysin Alo | VIR | + | + | + | + |
AMR—antimicrobial resistance gene. VIR—virulence factor gene.
Supplementary Materials
The following supporting information can be downloaded at:
References
1. Liu, Y.; Lai, Q.; Göker, M.; Meier-Kolthoff, J.P.; Wang, M.; Sun, Y.; Wang, L.; Shao, Z. Genomic insights into the taxonomic status of the Bacillus cereus group. Sci. Rep.; 2015; 5, 14082. [DOI: https://dx.doi.org/10.1038/srep14082] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26373441]
2. Swartz, M.N. Recognition and management of anthrax—An update. N. Engl. J. Med.; 2001; 345, pp. 1621-1626. [DOI: https://dx.doi.org/10.1056/NEJMra012892] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/11704686]
3. McKendrick, D.R. Anthrax and its transmission to humans. Cent. Afr. J. Med.; 1980; 26, pp. 126-129. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/7407891]
4. Kiple, K.F. The Cambridge World History of Human Disease; Cambridge University Press: Cambridge, UK, 1993.
5. Ringertz, S.H.; Høiby, E.A.; Jensenius, M.; Maehlen, J.; Caugant, D.A.; Myklebust, A.; Fossum, K. Injectional anthrax in a heroin skin-popper. Lancet; 2000; 356, pp. 1574-1575. [DOI: https://dx.doi.org/10.1016/S0140-6736(00)03133-0] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/11075776]
6. Jernigan, J.A.; Stephens, D.S.; Ashford, D.A.; Omenaca, C.; Topiel, M.S.; Galbraith, M.; Tapper, M.; Fisk, T.L.; Zaki, S.; Popovic, T. et al. Anthrax Bioterrorism Investigation Team. Bioterrorism-related inhalational anthrax: The first 10 cases reported in the United States. Emerg. Infect Dis.; 2001; 7, pp. 933-944. [DOI: https://dx.doi.org/10.3201/eid0706.010604] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/11747719]
7. Fowler, R.A.; Shafazand, S. Anthrax Bioterrorism: Prevention, Diagnosis and Management Strategies. J. Bioterr. Biodef.; 2011; 2, 107. [DOI: https://dx.doi.org/10.4172/2157-2526.1000107]
8. World Health OrganizationFood and Agriculture Organization of the United NationsWorld Organisation for Animal Health. Anthrax in Humans and Animals; WHO: Geneva, Switzerland, 2008; Available online: https://apps.who.int/iris/handle/10665/97503 (accessed on 20 March 2024).
9. Jacotot, H.; Virat, B. Long duration of viability of spores of Bacillus anthracis. Ann. Inst. Pasteur; 1954; 87, pp. 215-217.
10. Wilson, J.B.; Russell, K.E. Isolation of Bacillus anthracis from soil stored 60 years. J. Bacteriol.; 1964; 87, pp. 237-238. [DOI: https://dx.doi.org/10.1128/jb.87.1.237-238.1964] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/14102868]
11. De Vos, V. The ecology of anthrax in the Kruger National Park, South Africa. Salisbury Med. Bull.; 1990; 68S, pp. 19-23.
12. Sterne, M. The effects of different carbon dioxide concentrations on the growth of virulent anthrax strains. Pathogenicity and immunity tests on guinea pigs and sheep with anthrax variants derived from virulent strains. Onderstepoort J. Vet. Sci. Anim. Ind.; 1937; 9, pp. 49-67.
13. Sterne, M. The use of anthrax vaccines prepared from avirulent (uncapsulated) variants of Bacillus anthracis. Onderstepoort J. Vet. Sci. Anim. Ind.; 1939; 13, pp. 307-312.
14. Mikesell, P.; Ivins, B.E.; Ristroph, J.D.; Dreier, T.M. Evidence for plasmid-mediated toxin production in Bacillus anthracis. Infect. Immun.; 1983; 39, pp. 371-376. [DOI: https://dx.doi.org/10.1128/iai.39.1.371-376.1983]
15. Okinaka, R.; Cloud, K.; Hampton, O.; Hoffmaster, A.; Hill, K.; Keim, P.; Koehler, T.; Lamke, G.; Kumano, S.; Manter, D. et al. Sequence, assembly and analysis of pX01 and pX02. J. Appl. Microbiol.; 1999; 87, pp. 261-262. [DOI: https://dx.doi.org/10.1046/j.1365-2672.1999.00883.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/10475962]
16. Young, J.A.; Collier, R.J. Anthrax toxin: Receptor binding, internalization, pore formation, and translocation. Annu. Rev. Biochem.; 2007; 76, pp. 243-265. [DOI: https://dx.doi.org/10.1146/annurev.biochem.75.103004.142728]
17. Candela, T.; Fouet, A. Poly-gamma-glutamate in bacteria. Mol. Microbiol.; 2006; 60, pp. 1091-1098. [DOI: https://dx.doi.org/10.1111/j.1365-2958.2006.05179.x]
18. Green, B.D.; Battisti, L.; Koehler, T.M.; Thorne, C.B.; Ivins, B.E. Demonstration of a capsule plasmid in Bacillus anthracis. Infect. Immun.; 1985; 49, pp. 291-297. [DOI: https://dx.doi.org/10.1128/iai.49.2.291-297.1985]
19. Henderson, I.; Duggleby, C.J.; Turnbull, P.C. Differentiation of Bacillus anthracis from other Bacillus cereus group bacteria with the PCR. Int. J. Syst. Bacteriol.; 1994; 44, pp. 99-105. [DOI: https://dx.doi.org/10.1099/00207713-44-1-99] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/8123566]
20. Klee, S.R.; Ozel, M.; Appel, B.; Boesch, C.; Ellerbrok, H.; Jacob, D.; Holland, G.; Leendertz, F.H.; Pauli, G.; Grunow, R. et al. Characterization of Bacillus anthracis-like bacteria isolated from wild great apes from Cote d’Ivoire and Cameroon. J. Bacteriol.; 2006; 188, pp. 5333-5344. [DOI: https://dx.doi.org/10.1128/JB.00303-06] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16855222]
21. Rasko, D.A.; Rosovitz, M.J.; Økstad, O.A.; Fouts, D.E.; Jiang, L.; Cer, R.Z.; Kolstø, A.B.; Gill, S.R.; Ravel, J. Complete sequence analysis of novel plasmids from emetic and periodontal Bacillus cereus isolates reveals a common evolutionary history among the B. cereus-group plasmids, including Bacillus anthracis pXO1. J. Bacteriol.; 2007; 189, pp. 52-64. [DOI: https://dx.doi.org/10.1128/JB.01313-06]
22. Van der Auwera, G.A.; Andrup, L.; Mahillon, J. Conjugative plasmid pAW63 brings new insights into the genesis of the Bacillus anthracis virulence plasmid pXO2 and of the Bacillus thuringiensis plasmid pBT9727. BMC Genom.; 2005; 6, 103. [DOI: https://dx.doi.org/10.1186/1471-2164-6-103]
23. Keim, P.; Price, L.B.; Klevytska, A.M.; Smith, K.L.; Schupp, J.M.; Okinaka, R.; Jackson, P.J.; Hugh-Jones, M.E. Multiple-locus variable-number tandem repeat analysis reveals genetic relationships within Bacillus anthracis. J. Bacteriol.; 2000; 182, pp. 2928-2936. [DOI: https://dx.doi.org/10.1128/JB.182.10.2928-2936.2000] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/10781564]
24. Keim, P.; Gruendike, J.M.; Klevytska, A.M.; Schupp, J.M.; Challacombe, J.; Okinaka, R. The genome and variation of Bacillus anthracis. Mol. Aspects Med.; 2009; 30, pp. 397-405. [DOI: https://dx.doi.org/10.1016/j.mam.2009.08.005] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19729033]
25. Kolstø, A.B.; Tourasse, N.J.; Økstad, O.A. What sets Bacillus anthracis apart from other Bacillus species?. Annu. Rev. Microbiol.; 2009; 63, pp. 451-476. [DOI: https://dx.doi.org/10.1146/annurev.micro.091208.073255] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19514852]
26. Andersen, G.L.; Simchock, J.M.; Wilson, K.H. Identification of a region of genetic variability among Bacillus anthracis strains and related species. J. Bacteriol.; 1996; 178, pp. 377-384. [DOI: https://dx.doi.org/10.1128/jb.178.2.377-384.1996]
27. Van Ert, M.N.; Easterday, W.R.; Huynh, L.Y.; Okinaka, R.T.; Hugh-Jones, M.E.; Ravel, J.; Zanecki, S.R.; Pearson, T.; Simonson, T.S.; U’Ren, J.M. et al. Global genetic population structure of Bacillus anthracis. PLoS ONE; 2007; 2, e461. [DOI: https://dx.doi.org/10.1371/journal.pone.0000461] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17520020]
28. Thierry, S.; Tourterel, C.; Le Flèche, P.; Derzelle, S.; Dekhil, N.; Mendy, C.; Colaneri, C.; Vergnaud, G.; Madani, N. Genotyping of French Bacillus anthracis strains based on 31-loci multi locus VNTR analysis: Epidemiology, marker evaluation, and update of the internet genotype database. PLoS ONE; 2014; 9, e95131. [DOI: https://dx.doi.org/10.1371/journal.pone.0095131] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24901417]
29. Abdel-Glil, M.Y.; Chiaverini, A.; Garofolo, G.; Fasanella, A.; Parisi, A.; Harmsen, D.; Jolley, K.A.; Elschner, M.C.; Tomaso, H.; Linde, J. et al. A Whole-Genome-Based Gene-by-Gene Typing System for Standardized High-Resolution Strain Typing of Bacillus anthracis. J. Clin. Microbiol.; 2021; 59, e0288920. [DOI: https://dx.doi.org/10.1128/JCM.02889-20] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33827898]
30. Sahl, J.W.; Pearson, T.; Okinaka, R.; Schupp, J.M.; Gillece, J.D.; Heaton, H.; Birdsell, D.; Hepp, C.; Fofanov, V.; Noseda, R. et al. A Bacillus anthracis Genome Sequence from the Sverdlovsk 1979 Autopsy Specimens. mBio; 2016; 7, e01501-16. [DOI: https://dx.doi.org/10.1128/mBio.01501-16] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27677796]
31. Habrun, B.; Račić, I.; Kompes, G.; Špičić, S.; Benić, M.; Mihaljević, Ž.; Cvetnić, Ž. The antimicrobial susceptibility and virulence factors of Bacillus anthracis strains isolated in Croatia. Vet. Med.; 2011; 56, pp. 22-27. [DOI: https://dx.doi.org/10.17221/1570-VETMED]
32. Ministry of Agriculture of the Republic of Croatia. Available online: http://www.veterinarstvo.hr/default.aspx?id=185 (accessed on 7 July 2024).
33. Šoštarić, B.; Habrun, B.; Lipej, Z.; Vicković, I. Prikaz Epizootije Bedrenice u Stadu Ovaca i Koza. Veterinarski Dani, Zagreb, Hrvatska, 17.10.2002–20.10.2002; Hrvatska veterinarska komora: Zagreb, Croatia, 2002; pp. 31-32.
34. Habrun, B.; Kompes, G.; Šeol, B. Anthrax—Epizootic case in Bobovac (Sunja). Vet. Stanica; 2009; 40, pp. 381-388.
35. Miškić, T.; Lohman Janković, I. Epizootiološke Značajke Bedrenice na Području Parka Prirode Lonjsko Polje 2022. Sažetci, BEDRENICA I AFRIČKA SVINJSKA KUGA U HRVATSKOJ, Zagreb, Hrvatska, 13.10.2023; Hrvatska Akademija Znanosti i Umjetnosti: Zagreb, Croatia, 2023; ISBN 978-953-347-371-0 Available online: https://www.info.hazu.hr/wp-content/uploads/2023/10/2023-HAZU-Zbornik-sazetaka-Bedrenica-i-africka-svinjska-kuga.pdf (accessed on 20 March 2024).
36. Schmid, G.; Kaufmann, A. Anthrax in Europe: Its epidemiology, clinical characteristics, and role in bioterrorism. Clin. Microbiol. Infect.; 2002; 8, pp. 479-488. [DOI: https://dx.doi.org/10.1046/j.1469-0691.2002.00500.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/12197870]
37. Nicastri, E.; Vairo, F.; Mencarini, P.; Battisti, A.; Agrati, C.; Cimini, E.; Carrara, S.; D’Arezzo, S.; Adone, R.; Vulcano, A. et al. Unexpected human cases of cutaneous anthrax in Latium region, Italy, August 2017: Integrated human-animal investigation of epidemiological, clinical, microbiological and ecological factors. Euro Surveill.; 2019; 24, 1800685. [DOI: https://dx.doi.org/10.2807/1560-7917.ES.2019.24.24.1800685] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31213220]
38. Hodnik, J.J.; Knific, T.; Starič, J.; Toplak, I.; Ocepek, M.; Hostnik, P.; Ježek, J. Overview of Slovenian Control Programmes for Cattle Diseases Not Regulated by the European Union. Front. Vet. Sci.; 2021; 8, 674515.Erratum in Front. Vet. Sci. 2022, 8, 835395 [DOI: https://dx.doi.org/10.3389/fvets.2021.674515] [DOI: https://dx.doi.org/10.3389/fvets.2021.835395] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34307524]
39. Orlos, Z.; Rakoczi, E.; Misak, O.; Lenart, B.; Ocsai, G.; Kovacs, I.; Gorzsas, S.; Kardos, L.; Lampe, Z.; Szilvassy, Z. et al. Outbreak of anthrax in adults and adolescents: A review of nine cases in a regional teaching hospital in East Hungary. Clin. Microbiol. Infect.; 2017; 23, pp. 52-54. [DOI: https://dx.doi.org/10.1016/j.cmi.2015.11.023] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26794029]
40. Arapović, J.; Skočibusić, S.; Jelavić, B.; Ivanković, H.B.; Jurić, M.; Mamić, D.; Grgić, S.; Lesko, J.; Leventić, M.; Soldo, I. et al. Two cases of human cutaneous anthrax in Bosnia and Herzegovina, September 2014. Euro Surveill.; 2015; 20, pp. 2-4. [DOI: https://dx.doi.org/10.2807/1560-7917.ES2015.20.7.21039] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25719961]
41. Bower, W.A.; Yu, Y.; Person, M.K.; Parker, C.M.; Kennedy, J.L.; Sue, D.; Hesse, E.M.; Cook, R.; Bradley, J.; Bulitta, J.B. et al. CDC Guidelines for the Prevention and Treatment of Anthrax 2023. MMWR Recomm Rep.; 2023; 72, pp. 1-47. [DOI: https://dx.doi.org/10.15585/mmwr.rr7206a1]
42. Severn, M. A fatal case of pulmonary anthrax. Br. Med. J.; 1976; 27, 748. [DOI: https://dx.doi.org/10.1136/bmj.1.6012.748]
43. Cavallo, J.D.; Ramisse, F.; Girardet, M.; Vaissaire, J.; Mock, M.; Hernandez, E. Antibiotic susceptibilities of 96 isolates of Bacillus anthracis isolated in France between 1994 and 2000. Antimicrob. Agents Chemother.; 2002; 46, pp. 2307-2309. [DOI: https://dx.doi.org/10.1128/AAC.46.7.2307-2309.2002]
44. Maxson, T.; Kongphet-Tran, T.; Mongkolrattanothai, T.; Travis, T.; Hendricks, K.; Parker, C.; McLaughlin, H.P.; Bugrysheva, J.; Ambrosio, F.; Michel, P. et al. Systematic Review of In Vitro Antimicrobial Susceptibility Testing for Bacillus anthracis, 1947–2019. Clin. Infect. Dis.; 2022; 75, pp. S373-S378. [DOI: https://dx.doi.org/10.1093/cid/ciac520]
45. Manzulli, V.; Fasanella, A.; Parisi, A.; Serrecchia, L.; Donatiello, A.; Rondinone, V.; Caruso, M.; Zange, S.; Tscherne, A.; Decaro, N. et al. Evaluation of in vitro antimicrobial susceptibility of Bacillus anthracis strains isolated during anthrax outbreaks in Italy from 1984 to 2017. J. Vet. Sci.; 2019; 20, pp. 58-62. [DOI: https://dx.doi.org/10.4142/jvs.2019.20.1.58]
46. Chen, Y.; Succi, J.; Tenover, F.C.; Koehler, T.M. Beta-lactamase genes of the penicillin-susceptible Bacillus anthracis Sterne strain. J. Bacteriol.; 2003; 185, pp. 823-830. [DOI: https://dx.doi.org/10.1128/JB.185.3.823-830.2003] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/12533457]
47. Materon, I.C.; Queenan, A.M.; Koehler, T.M.; Bush, K.; Palzkill, T. Biochemical characterization of beta-lactamases Bla1 and Bla2 from Bacillus anthracis. Antimicrob. Agents Chemother.; 2003; 47, pp. 2040-2042. [DOI: https://dx.doi.org/10.1128/AAC.47.6.2040-2042.2003] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/12760895]
48. Read, T.D.; Peterson, S.N.; Tourasse, N.; Baillie, L.W.; Paulsen, I.T.; Nelson, K.E.; Tettelin, H.; Fouts, D.E.; Eisen, J.A.; Gill, S.R. et al. The genome sequence of Bacillus anthracis Ames and comparison to closely related bacteria. Nature; 2003; 423, pp. 81-86. [DOI: https://dx.doi.org/10.1038/nature01586] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/12721629]
49. Slamti, L.; Perchat, S.; Gominet, M.; Vilas-Bôas, G.; Fouet, A.; Mock, M.; Sanchis, V.; Chaufaux, J.; Gohar, M.; Lereclus, D. Distinct mutations in PlcR explain why some strains of the Bacillus cereus group are nonhemolytic. J. Bacteriol.; 2004; 186, pp. 3531-3538. [DOI: https://dx.doi.org/10.1128/JB.186.11.3531-3538.2004] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15150241]
50. Gohar, M.; Økstad, O.A.; Gilois, N.; Sanchis, V.; Kolstø, A.B.; Lereclus, D. Two-dimensional electrophoresis analysis of the extracellular proteome of Bacillus cereus reveals the importance of the PlcR regulon. Proteomics; 2002; 2, pp. 784-791. [DOI: https://dx.doi.org/10.1002/1615-9861(200206)2:6<784::AID-PROT784>3.0.CO;2-R] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/12112862]
51. Agaisse, H.; Gominet, M.; Okstad, O.A.; Kolstø, A.B.; Lereclus, D. PlcR is a pleiotropic regulator of extracellular virulence factor gene expression in Bacillus thuringiensis. Mol. Microbiol.; 1999; 32, pp. 1043-1053. [DOI: https://dx.doi.org/10.1046/j.1365-2958.1999.01419.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/10361306]
52. Sastalla, I.; Maltese, L.M.; Pomerantseva, O.M.; Pomerantsev, A.P.; Keane-Myers, A.; Leppla, S.H. Activation of the latent PlcR regulon in Bacillus anthracis. Microbiology; 2010; 156, pp. 2982-2993. [DOI: https://dx.doi.org/10.1099/mic.0.041418-0] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20688829]
53. Ross, C.L.; Thomason, K.S.; Koehler, T.M. An extracytoplasmic function sigma factor controls beta-lactamase gene expression in Bacillus anthracis and other Bacillus cereus group species. J. Bacteriol.; 2009; 191, pp. 6683-6693. [DOI: https://dx.doi.org/10.1128/JB.00691-09] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19717606]
54. Ho, T.D.; Nauta, K.M.; Müh, U.; Ellermeier, C.D. Activation of the Extracytoplasmic Function σ Factor σP by β-Lactams in Bacillus thuringiensis Requires the Site-2 Protease RasP. mSphere; 2019; 4, e00511-19. [DOI: https://dx.doi.org/10.1128/mSphere.00511-19]
55. Jones, M.E.; Goguen, J.; Critchley, I.A.; Draghi, D.C.; Karlowsky, J.A.; Sahm, D.F.; Porschen, R.; Patra, G.; DelVecchio, V.G. Antibiotic susceptibility of isolates of Bacillus anthracis, a bacterial pathogen with the potential to be used in biowarfare. Clin. Microbiol. Infect.; 2003; 9, pp. 984-986. [DOI: https://dx.doi.org/10.1046/j.1469-0691.2003.00775.x]
56. Luna, V.A.; King, D.S.; Gulledge, J.; Cannons, A.C.; Amuso, P.T.; Cattani, J. Susceptibility of Bacillus anthracis, Bacillus cereus, Bacillus mycoides, Bacillus pseudomycoides and Bacillus thuringiensis to 24 antimicrobials using Sensititre automated microbroth dilution and Etest agar gradient diffusion methods. J. Antimicrob. Chemother.; 2007; 60, pp. 555-567. [DOI: https://dx.doi.org/10.1093/jac/dkm213]
57. Maho, A.; Rossano, A.; Hächler, H.; Holzer, A.; Schelling, E.; Zinsstag, J.; Hassane, M.H.; Toguebaye, B.S.; Akakpo, A.J.; Van Ert, M. et al. Antibiotic susceptibility and molecular diversity of Bacillus anthracis strains in Chad: Detection of a new phylogenetic subgroup. J. Clin. Microbiol.; 2006; 44, pp. 3422-3425. [DOI: https://dx.doi.org/10.1128/JCM.01269-06]
58. Swaney, S.M.; Aoki, H.; Ganoza, M.C.; Shinabarger, D.L. The oxazolidinone linezolid inhibits initiation of protein synthesis in bacteria. Antimicrob. Agents Chemother.; 1998; 42, pp. 3251-3255. [DOI: https://dx.doi.org/10.1128/AAC.42.12.3251]
59. Pomerantsev, A.P.; Sukovatova, L.V.; Marinin, L.I. Izuchenie svoĭstv Rif-R-populiatsii sibireiazvennogo mikroba [Characterization of a Rif-R population of Bacillus anthracis]. Antibiot. I Khimioterapiia Antibiot. Chemoter.; 1993; 38, pp. 34-38.
60. Kim, H.S.; Choi, E.C.; Kim, B.K. A macrolide-lincosamide-streptogramin B resistance determinant from Bacillus anthracis 590: Cloning and expression of ermJ. J. Gen. Microbiol.; 1993; 139, pp. 601-607. [DOI: https://dx.doi.org/10.1099/00221287-139-3-601]
61. Sejvar, J.J.; Tenover, F.C.; Stephens, D.S. Management of anthrax meningitis. Lancet Infect. Dis.; 2005; 5, pp. 287-295. [DOI: https://dx.doi.org/10.1016/S1473-3099(05)70113-4]
62. Athamna, A.; Massalha, M.; Athamna, M.; Nura, A.; Medlej, B.; Ofek, I.; Bast, D.; Rubinstein, E. In vitro susceptibility of Bacillus anthracis to various antibacterial agents and their time-kill activity. J. Antimicrob. Chemother.; 2004; 53, pp. 247-251. [DOI: https://dx.doi.org/10.1093/jac/dkh016]
63. Odendaal, M.W.; Pieterson, P.M.; de Vos, V.; Botha, A.D. The antibiotic sensitivity patterns of Bacillus anthracis isolated from the Kruger National Park. Onderstepoort J. Vet. Res.; 1991; 58, pp. 17-19.
64. Heine, H.S.; Bassett, J.; Miller, L.; Purcell, B.K.; Byrne, W.R. Efficacy of Daptomycin against Bacillus anthracis in a murine model of anthrax spore inhalation. Antimicrob. Agents Chemother.; 2010; 54, pp. 4471-4473. [DOI: https://dx.doi.org/10.1128/AAC.00210-10]
65. Burckhardt, R.M.; Escalante-Semerena, J.C. Insights into the Function of the N-Acetyltransferase SatA That Detoxifies Streptothricin in Bacillus subtilis and Bacillus anthracis. Appl. Environ. Microbiol.; 2019; 85, e03029-18. [DOI: https://dx.doi.org/10.1128/AEM.03029-18]
66. Bruce, S.A.; Huang, Y.H.; Kamath, P.L.; van Heerden, H.; Turner, W.C. The roles of antimicrobial resistance, phage diversity, isolation source and selection in shaping the genomic architecture of Bacillus anthracis. Microb. Genom.; 2021; 7, 000616. [DOI: https://dx.doi.org/10.1099/mgen.0.000616] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34402777]
67. Chiaverini, A.; Abdel-Glil, M.Y.; Linde, J.; Galante, D.; Rondinone, V.; Fasanella, A.; Cammà, C.; D’Alterio, N.; Garofolo, G.; Tomaso, H. Whole Genome Sequencing for Studying Bacillus anthracis from an Outbreak in the Abruzzo Region of Italy. Microorganisms; 2020; 8, 87. [DOI: https://dx.doi.org/10.3390/microorganisms8010087] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31936409]
68. Uchida, I.; Hornung, J.M.; Thorne, C.B.; Klimpel, K.R.; Leppla, S.H. Cloning and characterization of a gene whose product is a trans-activator of anthrax toxin synthesis. J. Bacteriol.; 1993; 175, pp. 5329-5338. [DOI: https://dx.doi.org/10.1128/jb.175.17.5329-5338.1993] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/8366021]
69. Koehler, T.M.; Dai, Z.; Kaufman-Yarbray, M. Regulation of the Bacillus anthracis protective antigen gene: CO2 and a trans-acting element activate transcription from one of two promoters. J. Bacteriol.; 1994; 176, pp. 586-595. [DOI: https://dx.doi.org/10.1128/jb.176.3.586-595.1994] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/8300513]
70. Chung, M.C.; Popova, T.G.; Millis, B.A.; Mukherjee, D.V.; Zhou, W.; Liotta, L.A.; Petricoin, E.F.; Chandhoke, V.; Bailey, C.; Popov, S.G. Secreted neutral metalloproteases of Bacillus anthracis as candidate pathogenic factors. J. Biol. Chem.; 2006; 281, pp. 31408-31418. [DOI: https://dx.doi.org/10.1016/S0021-9258(19)84053-X] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16926147]
71. Chung, M.C.; Jorgensen, S.C.; Popova, T.G.; Tonry, J.H.; Bailey, C.L.; Popov, S.G. Activation of plasminogen activator inhibitor implicates protease InhA in the acute-phase response to Bacillus anthracis infection. J. Med. Microbiol.; 2009; 58, pp. 737-744. [DOI: https://dx.doi.org/10.1099/jmm.0.007427-0] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19429749]
72. Tonry, J.H.; McNichol, B.A.; Ramarao, N.; Chertow, D.S.; Kim, K.S.; Stibitz, S.; Schneewind, O.; Kashanchi, F.; Bailey, C.L.; Popov, S. et al. Bacillus anthracis protease InhA regulates BslA-mediated adhesion in human endothelial cells. Cell Microbiol.; 2012; 14, pp. 1219-1230. [DOI: https://dx.doi.org/10.1111/j.1462-5822.2012.01791.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22452315]
73. Gomis-Rüth, F.X. Catalytic domain architecture of metzincin metalloproteases. J. Biol. Chem.; 2009; 284, pp. 15353-15357. [DOI: https://dx.doi.org/10.1074/jbc.R800069200]
74. Ramarao, N.; Lereclus, D. The InhA1 metalloprotease allows spores of the B. cereus group to escape macrophages. Cell Microbiol.; 2005; 7, pp. 1357-1364. [DOI: https://dx.doi.org/10.1111/j.1462-5822.2005.00562.x]
75. Oda, M.; Hashimoto, M.; Takahashi, M.; Ohmae, Y.; Seike, S.; Kato, R.; Fujita, A.; Tsuge, H.; Nagahama, M.; Ochi, S. et al. Role of sphingomyelinase in infectious diseases caused by Bacillus cereus. PLoS ONE; 2012; 7, e38054. [DOI: https://dx.doi.org/10.1371/journal.pone.0038054]
76. Cai, Y.; Huang, T.; Xu, Y.; Zhou, G.; Zou, P.; Zeng, G.; Liu, X. Genetic and genomic diversity of NheABC locus from Bacillus strains. Arch. Microbiol.; 2017; 199, pp. 775-785. [DOI: https://dx.doi.org/10.1007/s00203-017-1350-9] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28283680]
77. Asano, S.I.; Nukumizu, Y.; Bando, H.; Iizuka, T.; Yamamoto, T. Cloning of novel enterotoxin genes from Bacillus cereus and Bacillus thuringiensis. Appl. Environ. Microbiol.; 1997; 63, pp. 1054-1057. [DOI: https://dx.doi.org/10.1128/aem.63.3.1054-1057.1997] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/9055420]
78. Boonchai, N.; Asano, S.I.; Bando, H.; Wiwat, C. Study on cytotoxicity and nucleotide sequences of enterotoxin FM of Bacillus cereus isolated from various food sources. J. Med. Assoc. Thai.; 2008; 91, pp. 1425-1432.
79. Castiaux, V.; Laloux, L.; Schneider, Y.J.; Mahillon, J. Screening of Cytotoxic B. cereus on Differentiated Caco-2 Cells and in Co-Culture with Mucus-Secreting (HT29-MTX) Cells. Toxins; 2016; 8, 320. [DOI: https://dx.doi.org/10.3390/toxins8110320] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27827957]
80. Tourasse, N.J.; Jolley, K.A.; Kolstø, A.B.; Økstad, O.A. Core genome multilocus sequence typing scheme for Bacillus cereus group bacteria. Res. Microbiol.; 2023; 174, 104050. [DOI: https://dx.doi.org/10.1016/j.resmic.2023.104050] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36893969]
81. Simonson, T.S.; Okinaka, R.T.; Wang, B.; Easterday, W.R.; Huynh, L.; U’Ren, J.M.; Dukerich, M.; Zanecki, S.R.; Kenefic, L.J.; Beaudry, J. et al. Bacillus anthracis in China and its relationship to worldwide lineages. BMC Microbiol.; 2009; 9, 71. [DOI: https://dx.doi.org/10.1186/1471-2180-9-71] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19368722]
82. Antwerpen, M.; Ilin, D.; Georgieva, E.; Meyer, H.; Savov, E.; Frangoulidis, D. MLVA and SNP analysis identified a unique genetic cluster in Bulgarian Bacillus anthracis strains. Eur. J. Clin. Microbiol. Infect. Dis.; 2011; 30, pp. 923-930. [DOI: https://dx.doi.org/10.1007/s10096-011-1177-2] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21279731]
83. Rondinone, V.; Serrecchia, L.; Parisi, A.; Fasanella, A.; Manzulli, V.; Cipolletta, D.; Galante, D. Genetic characterization of Bacillus anthracis strains circulating in Italy from 1972 to 2018. PLoS ONE; 2020; 15, e0227875. [DOI: https://dx.doi.org/10.1371/journal.pone.0227875] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31931511]
84. Smith, K.L.; DeVos, V.; Bryden, H.; Price, L.B.; Hugh-Jones, M.E.; Keim, P. Bacillus anthracis diversity in Kruger National Park. J. Clin. Microbiol.; 2000; 38, pp. 3780-3784. [DOI: https://dx.doi.org/10.1128/JCM.38.10.3780-3784.2000]
85. Fouet, A.; Smith, K.L.; Keys, C.; Vaissaire, J.; Le Doujet, C.; Lévy, M.; Mock, M.; Keim, P. Diversity among French Bacillus anthracis isolates. J. Clin. Microbiol.; 2002; 40, pp. 4732-4734. [DOI: https://dx.doi.org/10.1128/JCM.40.12.4732-4734.2002]
86. Gierczyński, R.; Kałuzewski, S.; Rakin, A.; Jagielski, M.; Zasada, A.; Jakubczak, A.; Borkowska-Opacka, B.; Rastawicki, W. Intriguing diversity of Bacillus anthracis in eastern Poland--the molecular echoes of the past outbreaks. FEMS Microbiol. Lett.; 2004; 239, pp. 235-240. [DOI: https://dx.doi.org/10.1016/j.femsle.2004.08.038] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15476971]
87. Derzelle, S.; Laroche, S.; Le Flèche, P.; Hauck, Y.; Thierry, S.; Vergnaud, G.; Madani, N. Characterization of genetic diversity of Bacillus anthracis in France by using high-resolution melting assays and multilocus variable-number tandem-repeat analysis. J. Clin. Microbiol.; 2011; 49, pp. 4286-4292. [DOI: https://dx.doi.org/10.1128/JCM.05439-11] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21998431]
88. Garofolo, G.; Serrecchia, L.; Corrò, M.; Fasanella, A. Anthrax phylogenetic structure in Northern Italy. BMC Res. Notes; 2011; 4, 273. [DOI: https://dx.doi.org/10.1186/1756-0500-4-273] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21801397]
89. Girault, G.; Blouin, Y.; Vergnaud, G.; Derzelle, S. High-throughput sequencing of Bacillus anthracis in France: Investigating genome diversity and population structure using whole-genome SNP discovery. BMC Genom.; 2014; 15, 288. [DOI: https://dx.doi.org/10.1186/1471-2164-15-288] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24734872]
90. Derzelle, S.; Thierry, S. Genetic diversity of Bacillus anthracis in Europe: Genotyping methods in forensic and epidemiologic investigations. Biosecur. Bioterror.; 2013; 1, pp. S166-S176. [DOI: https://dx.doi.org/10.1089/bsp.2013.0003]
91. Croatian Food Agency. Available online: https://www.hzjz.hr/wp-content/uploads/2017/11/Godi%C5%A1nje-izvje%C5%A1%C4%87e-o-zoonozama-2015_16.pdf (accessed on 7 July 2024).
92. Croatian Institute for Public Health. Croatian Health Statistics Yearbook. 2022; Available online: https://www.hzjz.hr/wp-content/uploads/2024/05/HZSLj_2022_12-2023.pdf (accessed on 7 July 2024).
93. Pajtlar, S.; Ančić-Birač, S.; Starčević, M. Epidemija Antraksa u Parku Prirode Lonjsko Polje 2022, Sažetci, BEDRENICA I AFRIČKA SVINJSKA KUGA U HRVATSKOJ, Zagreb, Hrvatska, 13.10.2023; Hrvatska Akademija Znanosti I Umjetnosti: Zagreb, Croatia, 2023; ISBN 978-953-347-371-0 Available online: https://www.info.hazu.hr/wp-content/uploads/2023/10/2023-HAZU-Zbornik-sazetaka-Bedrenica-i-africka-svinjska-kuga.pdf (accessed on 20 March 2024).
94. Beyer, W.; Glöckner, P.; Otto, J.; Böhm, R. A nested PCR method for the detection of Bacillus anthracis in environmental samples collected from former tannery sites. Microbiol. Res.; 1995; 150, pp. 179-186. [DOI: https://dx.doi.org/10.1016/S0944-5013(11)80054-6] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/7600011]
95. Hutson, R.A.; Duggleby, C.J.; Lowe, J.R.; Manchee, R.J.; Turnbull, P.C. The development and assessment of DNA and oligonucleotide probes for the specific detection of Bacillus anthracis. J. Appl. Bacteriol.; 1993; 75, pp. 463-472. [DOI: https://dx.doi.org/10.1111/j.1365-2672.1993.tb02803.x]
96. Clinical and Laboratory Standards Institute (CLSI). Methods for Antimicrobial Dilution and Disk Susceptibility Testing of Infrequently Isolated or Fastidious Bacteria M45; 3rd ed. CLSI: Wayne, PA, USA, 2016.
97. Clinical and Laboratory Standards Institute (CLSI). Performance Standards for Antimicrobial Susceptibility Testing, M100; 32nd ed. CLSI: Wayne, PA, USA, 2022.
98. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics; 2014; 30, pp. 2114-2120. [DOI: https://dx.doi.org/10.1093/bioinformatics/btu170] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24695404]
99. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D. et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol.; 2012; 19, pp. 455-477. [DOI: https://dx.doi.org/10.1089/cmb.2012.0021]
100. Gurevich, A.; Saveliev, V.; Vyahhi, N.; Tesler, G. QUAST: Quality assessment tool for genome assemblies. Bioinformatics; 2013; 29, pp. 1072-1075. [DOI: https://dx.doi.org/10.1093/bioinformatics/btt086]
101. Chen, L.; Zheng, D.; Liu, B.; Yang, J.; Jin, Q. VFDB 2016: Hierarchical and refined dataset for big data analysis—10 years on. Nucleic Acids Res.; 2016; 44, pp. D694-D697. [DOI: https://dx.doi.org/10.1093/nar/gkv1239] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26578559]
102. Feldgarden, M.; Brover, V.; Haft, D.H.; Prasad, A.B.; Slotta, D.J.; Tolstoy, I.; Tyson, G.H.; Zhao, S.; Hsu, C.H.; McDermott, P.F. et al. Validating the AMRFinder Tool and Resistance Gene Database by Using Antimicrobial Resistance Genotype-Phenotype Correlations in a Collection of Isolates. Antimicrob. Agents Chemother.; 2019; 63, e00483-19.Erratum in Antimicrob. Agents Chemother. 2020, 64, e00361-20 [DOI: https://dx.doi.org/10.1128/AAC.00483-19] [DOI: https://dx.doi.org/10.1128/AAC.00361-20] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31427293]
103. Johansson, M.H.K.; Bortolaia, V.; Tansirichaiya, S.; Aarestrup, F.M.; Roberts, A.P.; Petersen, T.N. Detection of mobile genetic elements associated with antibiotic resistance in Salmonella enterica using a newly developed web tool: MobileElementFinder. J. Antimicrob. Chemother.; 2021; 76, pp. 101-109. [DOI: https://dx.doi.org/10.1093/jac/dkaa390] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33009809]
104. Jünemann, S.; Sedlazeck, F.J.; Prior, K.; Albersmeier, A.; John, U.; Kalinowski, J.; Mellmann, A.; Goesmann, A.; von Haeseler, A.; Stoye, J. et al. Updating benchtop sequencing performance comparison. Nat. Biotechnol.; 2013; 31, pp. 294-296. [DOI: https://dx.doi.org/10.1038/nbt.2522] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23563421]
105. Maestro, M. Bedrenica Domaćih Životinja u BIVŠOJ SAVSKOJ BANOVINI od 1931–1940 Godine. Ph.D. Thesis; Veterinarski fakultet Hrvatskog sveučilišta u Zagrebu: Zagreb, Croatia, 1946.
106. Anon. Zakona o uređenju veterinarstva u Kraljevinah Hrvatskoj i Slavoniji. Nar. Nov. (Tisk. Zavod Zagreb); 1888; 27, kolovoza
107. Mlinac, F. Pojavljivanje Bedrenice na Području Bivše Hrvatske i Slavonije. Ph.D. Thesis; Veterinarski fakultet Kraljevine Jugoslavije: Sveučiliše u Zagrebu, Zagreb, Croatia, 1930.
108. Zaharija, I. Zarazne Bolesti Domaćih Životinja; Školska knjiga: Zagreb, Croatia, 1978; pp. 111-112.
109. Fališevac, J.; Hellenbach, H. Antraks kao profesionalno oboljenje sa specijalnim osvrtom na prilike u Jugoslaviji. Arh. Hig. Rada Toksikol.; 1951; 2, pp. 428-464. Available online: https://hrcak.srce.hr/clanak/277177 (accessed on 15 September 2019).
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Details
; Duvnjak, Sanja 2
; Reil, Irena 2
; Mihaljević, Željko 3
; Habrun, Boris 1
; Benić, Miroslav 4 ; Cvetnić, Luka 4 ; Špičić, Silvio 2
; Bagarić, Antonela 1 1 Laboratory for General Bacteriology and Mycology, Department for Bacteriology and Parasitology, Croatian Veterinary Institute, 10000 Zagreb, Croatia;
2 Laboratory for Bacterial Zoonoses and Molecular Diagnostics of Bacterial Diseases, Department for Bacteriology and Parasitology, Croatian Veterinary Institute, 10000 Zagreb, Croatia;
3 Laboratory for Pathology, Department for Pathology, Croatian Veterinary Institute, 10000 Zagreb, Croatia;
4 Laboratory for Mastitis and Raw Milk Quality, Department for Bacteriology and Parasitology, Croatian Veterinary Institute, 10000 Zagreb, Croatia;