- Research article
- Open Access
SNR analysis: molecular investigation of an anthrax epidemic
BMC Veterinary Researchvolume 6, Article number: 11 (2010)
In Italy, anthrax is endemic but occurs sporadically. During the summer of 2004, in the Pollino National Park, Basilicata, Southern Italy, an anthrax epidemic consisting of 41 outbreaks occurred; it claimed the lives of 124 animals belonging to different mammal species. This study is a retrospective molecular epidemiological investigation carried out on 53 isolates collected during the epidemic. A 25-loci Multiple Locus VNTR Analysis (MLVA) MLVA was initially performed to define genetic relationships, followed by an investigation of genetic diversity between epidemic strains through Single Nucleotide Repeat (SNR) analysis.
53 Bacillus anthracis strains were isolated. The 25-loci MLVA analysis identified all of them as belonging to a single genotype, while the SNR analysis was able to detect the existence of five subgenotypes (SGTs), allowing a detailed epidemic investigation. SGT-1 was the most frequent (46/53); SGTs 2 (4/53), 3 (1/53) 4 (1/53) and 5 (1/53) were detected in the remaining seven isolates.
The analysis revealed the prevalent spread, during this epidemic, of a single anthrax clone. SGT-1 - widely distributed across the epidemic area and present throughout the period in question - may, thus, be the ancestral form. SGTs 2, 3 and 4 differed from SGT-1 at only one locus, suggesting that they could have evolved directly from the latter during the course of this epidemic. SGT-5 differed from the other SGTs at 2-3 loci. This isolate, thus, appears to be more distantly related to SGT-1 and may not be a direct descendant of the lineage responsible for the majority of cases in this epidemic. These data confirm the importance of molecular typing and subtyping methods for in-depth epidemiological analyses of anthrax epidemics.
In the region of Basilicata, Southern Italy, anthrax outbreaks are typically isolated, self containing, and involve unvaccinated herbivores. Epidemics are rare, and often occur when a rainy spring is followed by a dry summer [1–5]. During the spring and summer of 2004, as a result of such weather conditions in the Pollino National Park, an anthrax epidemic occurred. The affected area included 13 towns and involved 41 farms over an area of about 900 Km2, with a livestock population numbering about 7,000 cattle and 33,000 between sheep and goats. In 40 days, 81 cattle, 15 sheep, nine goats, eleven horses and eight red deer died (Figure 1). The anthrax epidemic evolved in three different phases. The first, counting 26 outbreaks, was the most critical. The second and third phases, with eight and six outbreaks, respectively, were less severe.
An additional outbreak preceded the epidemic by about one month . Several epidemiological factors may have contributed to this phenomenon. In this, as in other endemic areas, spores resulting from previous outbreaks may remain in the soil, thus facilitating the spread of anthrax among livestock through grazing [7, 8]. In addition, anthrax infected carcasses are seldom removed. These carcasses are not isolated from the wild animals populating the Pollino National Park (deer, wild boars), resulting in a persistent source of infection in the environment . Furthermore, the abundance, at this time of year, of both biting (e.g. tabanid) and non biting flies, which may act as mechanical vectors, could also have contributed to the persistence of anthrax [9–14].
Genetically, B. anthracis is a relatively homogeneous bacteria species. Not surprisingly, then, discriminating between strains isolated from epidemiologically linked outbreaks is not an easy task . Different studies typed and differentiated B. anthracis isolates using Single Nucleotide Polymorphisms (SNP) analysis and Multiple Locus VNTR analysis (MLVA) [16–22]. In an epidemic, however, these methodologies are not likely to find genetic variation. The Single Nucleotide Repeats (SNR) analysis described by Stratilo et al. increases the likelihood of differentiating closely related isolates . Unfortunately, due to the presence of poly-A sequences, such polymorphisms are difficult to detect both with electorphoretic fragment analysis and with direct sequencing.
In this study, a retrospective molecular epidemiological investigation was performed, comparing 25-loci MLVA and two SNR analyses. We applied the modified SNR technique described by Kenefic et al. (KEN-MTD) as well as Stratilo's original SNR method (STR-MTD), selecting the four loci with the highest diversity indices (D = 0.57-0.90; where D = 1-Σ [allele frequency]2) [23, 24]. The SNR primer panels used have two loci in common (CL33, CL12) and two distinct loci (STR-MTD: CL1, CL37) (KEN-MTD: CL10, CL35). Two different genetic analyzers (DNA sequencers) were used to verify results. This was done to address the technical difficulty in correct allele assignment.
SNR analysis identified five SGTs. Of 53 isolates, 46 were classified as SGT-1, compared to which four isolates, classified as SGT-2, had a single-base pair deletion corresponding to locus CL33; SGT-3 and SGT-4, each with a single isolate, exhibited insertions in loci CL33 and CL12, respectively, while SGT-5, again with a single isolate, differed from SGT-1 at three loci, with a five-base pair deletion at CL33, an insertion into the CL12 locus, and a deletion at CL10 (Table 1; figure 2).
Due to homoplasy, the high mutation rate of SNR loci, estimated at 10-4 per generation in B. anthracis, does not allow a correct definition of phylogenetic relationships between different isolates . SNR is, however, able to detect genetic diversity between closely related strains as would occur in an epidemic . Our method of analysis involved the initial use of 25-loci MLVA to define genetic relationships, and a subsequent investigation of genetic diversity between epidemic strains through SNR analysis.
The 25-loci MLVA assigned all strains to a single genotype within the cluster A1.a, evidence of their autochthonous origin. This particular genotype is frequently found in Basilicata. The two SNR analyses, STR-MTD and KEN-MTD, yielded the same result, identifying five SGTs. The difference between these methods in terms of the amplicon sizes obtained for the common loci is a consequence of the use of different primer pairs. We found the KEN-MTD more useful, as it allows for a multiplexed reaction and permits faster and more reliable detection of fragment sizes, owing to smaller amplicons sizes. Moreover, the KEN-MTD revealed an additional difference between SGT-5 and SGT-1, an allele in the CL10 (Table 1). Polymorphisms were only discovered among the loci with the highest diversity indices (CL33, CL12, CL10) .
SGT-1, the most common, was distributed across the epidemic area and present throughout the period under study. SGT-1 was also detected in soil from a grave site, and in the feces of wild boar collected in the same area, suggesting it is the probable ancestral strain. SGTs 2-4 were rarer (6/53) and not present in all phases of the epidemic. These SGTs could represent intra epidemic mutational steps descended from SGT-1. SGT-5, with three loci exhibiting fragment length polymorphism as compared to SGT-1, was the first to be isolated in the affected area, but may have descended from a different, closely related previous outbreak (Table 2).
The epidemic under study was characterized by a single anthrax clone. Although the epidemic spread over a large area, it involved extremely closely related isolates, 86.7% of which (SGT-1) were identical across the highly discriminating SNR markers. This epidemic was probably exacerbated by the absence of a careful monitoring of animals, making it possible for spores from a recent victim to spread in the environment. The mutational steps found in SGT-2, 3 and 4 may be associated with infective cycles subsequent to the first, while the presence of a more distantly related strain, SGT-5, testifies to the evolutionary history of this lineage in this region. This divergent mutational clone may have originated from different outbreaks in the past. The high throughput genotyping system used in this study proved to be a useful tool for the study of closely related B. anthracis strains, and is therefore potentially valuable not only for the study of epidemics, but also for other contexts requiring the characterization of closely related strains such as the study of biological contaminations through environmental isolates or the forensic investigation of bioterrorist events.
Bacillus anthracis isolates
In this study, we analyzed 53 B. anthracis strains associated with a single anthrax epidemic.
Each B. anthracis strain was streaked onto 5% sheep blood agar plates and then incubated at + 37°C for 24 hours. After heat inactivation (98°C for 20 min.), microbial DNA was extracted using DNAeasy Blood and Tissue kit (Qiagen), following the protocol for Gram positive bacteria.
25-loci MLVA and SNR analyses
We utilized 5' fluorescent-labelled oligos, deprotected and desalted, specifically selected for the VNTRs and for the SNRs used.
The 25 specific primer pairs for the MLVA were selected as described by Lista et al. . The eight specific primer pairs for SNR reactions were selected following Stratilo et al. and Kenefic et al. [23, 24] (Table 3).
MLVA PCRs were performed in four multiplex reactions in a final volume of 15 μl. The reaction mixture contained: 1× PCR reaction buffer (Roche), 1 U of Taq DNA polymerase (Roche), dNTPs (0.2 mM each), and appropriate concentrations of each primer as reported in Table 3. The thermocycling conditions were as follows: 96°C for 3 min; 36 cycles at 95°C for 20 s, at 60°C for 30 s, and at 72°C for 1 min; and finally, 72°C for 10 min.
The KEN-MTD PCR was performed in a multiplex reaction in a final volume of 25 μl containing 1× AmpliTaq Gold PCR buffer and 0.5 U of AmpliTaq Gold DNA polymerase (Applied Biosystems Inc.), 3,5 mM MgCl2, dNTPs (0.2 mM each), and appropriate concentrations of forward and reverse primers as reported in Table 3.
STR-MTD PCRs were performed in four singleplex reactions in a final volume of 25 μl containing 1× AmpliTaq Gold PCR buffer and 0.5 U of AmpliTaq Gold DNA polymerase (Applied Biosystems Inc.), 4 mM MgCl2, dNTPs (0.2 mM each), and appropriate concentrations of forward and reverse primers as reported in Table 3. The thermocycling conditions were as follows: 95°C for 5 min; 35 cycles at 94°C for 30 s, 60°C for 30 s, and 72°C for 30 s; and finally, 72°C for 7 min.
Automated genotype analysis
MLVA PCR products were diluted 1:5. Five μl of solution were added to a mix containing 40 μl of Sample Loading Solution (Beckman Coulter) and 0.5 μl of MapMarker 1000 size marker (BioVentures Inc.). Amplicons were separated by electrophoresis on a CEQ 8000 automated DNA Analysis System (Beckman Coulter) and sized by CEQ Fragment Analysis System software.
Amplified SNR PCR products were diluted 1:80 and subjected to capillary electrophoresis on ABI Prism 3130 genetic analyzer (Applied Biosystems) with 0.25 μl of GeneScan 120 Liz and 500 Liz size standards for KEN-MTD and STR-MTD amplicons, respectively, and sized by GeneMapper 4.0 (Applied Biosystems Inc.). DNA extracted from each sample was tested by two different laboratories: the Anthrax Reference Institute of Italy, and the Army Medical and Veterinary Research Institute.
Fasanella A, Van Ert M, Altamura SA, Garofolo G, Buonavoglia C, Leori G, Huynh L, Zanecki S, Keim P: Molecular diversity of Bacillus anthracis in Italy. J Clin Microbiol. 2005, 43: 3398-3401. 10.1128/JCM.43.7.3398-3401.2005.
Hugh-Jones ME, de Vos V: Anthrax and wildlife. Rev Sci Tech. 2002, 2: 359-389.
Hugh-Jones ME: 1996-97 Global Anthrax report. J Appl Microbiol. 1999, 87: 189-191. 10.1046/j.1365-2672.1999.00867.x.
Turner AJ, Galvin JW, Rubira RJ, Condron RJ, Brandley T: Experiences with vaccination and epidemiological investigations on an anthrax outbreak in Australia in 1997. J of Appl Microbiol. 1999, 87: 294-298. 10.1046/j.1365-2672.1999.00894.x.
Turnbull PC, Bell RH, Saigawa K, Munyenyembe FE, Mulenga CK: Anthrax in wildlife in Luangwa Valley, Zambia. Vet Rec. 1991, 128: 399-403.
Fasanella A, Palazzo L, Petrella A, Quaranta V, Romanelli B, Garofolo G: Anthrax in red deer (Cervus elaphus), Italy. Emerg Infect Dis. 2007, 13 (7): 1118-9.
Dragon DC, Rennie RP: ecology of anthrax spores: tough but not invincible. Can Vet J. 1999, 36 (5): The295-301.
Dragon DC, Bader DE, Mitchell J, Woollen N: Natural dissemination of Bacillus anthracis spores in Northern Canada. Appl Environ Microbiol. 2005, 71 (3): 1610-1615. 10.1128/AEM.71.3.1610-1615.2005.
Mitzmain MB: Experimental insect transmission of anthrax. Public Health Reports. 1914, 29: 75-7716.
Kraneveld FC, Djaenoedin R: Test on the dissemination of anthrax by Tabanus rubidus in horses and buffalo. Overgedrukt uit de Nederlands-Indische Bladen Voor Diergeneeskunde. 1940, 52: 339-80.
Rao NS, Mohiyudeen S: Tabanus flies as transmitters of anthrax - a field experience. Indian Veterinary Journal. 1958, 35: 348-353.
Davies JC: A major epidemic of anthrax in Zimbabwe. Part II. Cent Afr J Med. 1983, 29 (1): 8-12.
Turrel MJ, Knudson GB: Mechanical transmission of Bacillus anthracis by stable flies (Stomoxys calcitrans) and mosquitoes (Aedes aegypti and Aedes taeniorhynchus). Infection and immunity. 1987, 55: 1859-1961.
Hugh-Jones M, Blackburn J: The ecology of Bacillus anthracis. Molecular Aspects of Medicine. 2009, 30 (6): 356-67. 10.1016/j.mam.2009.08.003.
Keim P, van Ert MN, Pearson T, Vogler AJ, Huynh LY, Wagner DM: Anthrax molecular epidemiology and forensics: using the appropriate marker for different evolutionary scales. Infect Genet Evol. 2004, 4: 205-213. 10.1016/j.meegid.2004.02.005.
Keim P, Price LB, Klevytska AM, Smith KL, Schupp JM, Okinaka R, Jackson P, Hugh Jones ME: Multiple-Locus Variable-Number Tandem Repeat Analysis Reveals Genetic Relationships within Bacillus anthracis. J Bacteriology. 2000, 182: 2928-2936. 10.1128/JB.182.10.2928-2936.2000.
Lista F, Faggioni G, Valjevac S, Ciammaruconi A, Vaissaire J, le Doujet C, Gorge O, de Santis R, Carattoli A, Ciervo A, Fasanella A, Orsini F, D'Amelio R, Pourcel C, Cassone A, Vergnaud G: Genotyping of Bacillus anthracis strains based on automated capillary. BMC Microbiology. 2006, 6: 6-33. 10.1186/1471-2180-6-33.
Fouet A, Smith KL, 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: 4732-4734. 10.1128/JCM.40.12.4732-4734.2002.
Gierczynski R, Jakubczak A, Jagielski M: Extended multiple-locus variable tandem-repeat analysis of Bacillus anthracis isolated in Poland. Pol J Microbiol. 2009, 58 (1): 3-7.
Van Ert M, Easterday WR, Huynh LY, Okinaka RT, Hugh-Jones ME, Ravel J, Zanecki SR, Pearson T, Simonson TS, U'Ren JM, Kachur SM, Leadem-Dougherty RR, Rhoton SD, Zinser G, Farlow J, Coker PR, Smith KL, Wang B, Kenefic LJ, Fraser-Liggett CM, Wagner DM, Keim P: Global Genetic Population Structure of Bacillus anthracis. PlosOne. 2007, 5: 461-471.
Pearson T, Busch JD, Ravel J, Read TD, Rhoton SD, U'Ren JM, Simonson TS, Kachur SM, Leadem RR, Cardon ML, Van Ert MN, Huynh LY, Fraser CM, Keim P: Phylogenetic discovery bias in Bacillus anthracis using single-nucleotide polymorphisms from whole-genome sequencing. Procl Nat Acad Sci USA. 2004, 101 (37): 13536-41. 10.1073/pnas.0403844101.
Read TD, Salzberg SL, Pop M, Shumway M, Umayam L, Jiang L, Holtzapple E, Busch JD, Smith KL, Schupp JM, Solomon D, Keim P, Fraser CM: Comparative genome sequencing for discovery of novel polymorphisms in. Bacillus anthracis. 2002, 296 (5575): 2028-33.
Stratilo CW, Lewis CT, Bryden L, Mulvey MR, Bader D: Single-nucleotide repeat analysis for subtyping Bacillus anthracis. J Clin Microbiol. 2006, 44 (3): 777-8. 10.1128/JCM.44.3.777-782.2006.
Kenefic LJ, Beaudry J, Trim C, Huynh L, Zanecki S, Matthews M, Schupp J, Van Ert M, Keim P: A high resolution four-locus multiplex single nucleotide repeat (SNR) genotyping system in Bacillus anthracis. J Microbiol Methods. 2008, 73 (3): 269-72. 10.1016/j.mimet.2007.11.014.
This work was financed by the Italian Ministry of Health (DIAGNOVA funds, Ricerca Corrente 2006). We thank Angela Aceti, Giuseppe Stramaglia, Nicola Nigro and Rosa d'Errico for their excellent technical support.
GG conceived the study, participated in its design, coordinated the molecular tests and drafted the manuscript. AC participated in performing the 25-loci MLVA and the SNR analyses and helped to draft the manuscript. AF designed the study and coordinated the collection of samples. SS performed the SNR analysis and helped draft the manuscript. VP participated in performing the 25-loci MLVA. RA helped design the study. FL coordinated the work of the two laboratories and revised the final version of the manuscript. All authors read and approved the final manuscript.