Application of High-Density DNA Resequencing
Microarray for Detection and Characterization of
Botulinum Neurotoxin-Producing Clostridia

Jessica Vanhomwegen, Nicolas Berthet, Christelle Mazuet, Ghislaine Guigon, Tatiana Vallaeys, Rayna Stamboliyska, Philippe Dubois, Giulia Kennedy, Stewart Cole, Valerie Caro, et al.

► To cite this version:

Jessica Vanhomwegen, Nicolas Berthet, Christelle Mazuet, Ghislaine Guigon, Tatiana Vallaeys, et al. Application of High-Density DNA Resequencing Microarray for Detection and Characterization of Botulinum Neurotoxin-Producing Clostridia. PLoS ONE, 2013, 8 (6), pp.e67510. 10.1371/journal.pone.0067510. pasteur-01765377

HAL Id: pasteur-01765377

https://pasteur.hal.science/pasteur-01765377v1

Submitted on 12 Apr 2018

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# Application of High-Density DNA Resequencing Microarray for Detection and Characterization of Botulinum Neurotoxin-Producing Clostridia

Jessica Vanhomwegen¹*, Nicolas Berthet²,³, Christelle Mazuet⁴, Ghislaine Guigon⁵, Tatiana Vallaeys⁶, Rayna Stamboliyska⁷, Philippe Dubois¹, Giulia C. Kennedy⁸, Stewart T. Cole⁹, Valérie Caro⁵, Jean-Claude Manuguerra¹, Michel-Robert Popoff⁴

¹Institut Pasteur, Laboratory for Urgent Responses to Biological Threats, Paris, France, ²Institut Pasteur, Epidemiology and Pathophysiology of Oncogenic Viruses, Paris, France, ³CNRS, UMR3569, Paris, France, ⁴Institut Pasteur, Anaerobic Bacteria and Toxins Unit, Paris, France, ⁵Institut Pasteur, Genotyping of Pathogens and Public Health, Paris, France, ⁶CNRS - CC093 Université Montpellier II, UMR5119 Ecosystèmes lagunaires, Montpellier, France, ⁷Department of Biology II, University of Munich (LMU), Planegg-Martinsried, Germany, ⁸Department of Research and Development, Veracyte, Inc., South San Francisco, California, United States of America, ⁹Global Health Institute, Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland

## Abstract

**Background:** Clostridium botulinum and related clostridia express extremely potent toxins known as botulinum neurotoxins (BoNTs) that cause severe, potentially lethal intoxications in humans. These BoNT-producing bacteria are categorized in seven major toxinotypes (A through G) and several subtypes. The high diversity in nucleotide sequence and genetic organization of the gene cluster encoding the BoNT components poses a great challenge for the screening and characterization of BoNT-producing strains.

**Methodology/Principal Findings:** In the present study, we designed and evaluated the performances of a resequencing microarray (RMA), the PathogenId v2.0, combined with an automated data approach for the simultaneous detection and characterization of BoNT-producing clostridia. The unique design of the PathogenID v2.0 array allows the simultaneous detection and characterization of 48 sequences targeting the BoNT gene cluster components. This approach allowed successful identification and typing of representative strains of the different toxinotypes and subtypes, as well as the neurotoxin-producing C. botulinum strain in a naturally contaminated food sample. Moreover, the method allowed fine characterization of the different neurotoxin gene cluster components of all studied strains, including genomic regions exhibiting up to 24.65% divergence with the sequences tiled on the arrays.

**Conclusions/Significance:** The severity of the disease demands rapid and accurate means for performing risk assessments of BoNT-producing clostridia and for tracing potentials sources of contamination in outbreak situations. The RMA approach constitutes an essential higher echelon component in a diagnostics and surveillance pipeline. In addition, it is an important asset to characterise potential outbreak related strains, but also environment isolates, in order to obtain a better picture of the molecular epidemiology of BoNT-producing clostridia.

**Citation:** Vanhomwegen J, Berthet N, Mazuet C, Guigon G, Vallaey s T, et al. (2013) Application of High-Density DNA Resequencing Microarray for Detection and Characterization of Botulinum Neurotoxin-Producing Clostridia. PLoS ONE 8(6): e67510. doi:10.1371/journal.pone.0067510

**Editor:** Yung-Fu Chang, Cornell University, United States of America

**Received** May 21, 2012; **Accepted** May 22, 2013; **Published** June 20, 2013

**Copyright:** © 2013 Vanhomwegen et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** This program was supported by grant UC1 AI062613 (to GC Kennedy) from the U.S. National Institute of Allergy and Infectious Diseases, National Institute of Health; by the "Programme Transversal de Recherche" (DEVA No. 246), financed by the Institut Pasteur (Paris, France), sponsorship by the Fondation Total-Institut Pasteur and the Elisabeth Taub award from the Académie Nationale de Medecine. The funders had no role in study design, data analysis or preparation of the manuscript.

**Competing Interests:** One of the authors, Dr. Giulia C. Kennedy, is Chief Scientific Officer at Veracyte Inc, a commercial company. This does not alter the authors' adherence to all the PLOS ONE policies on sharing data and materials.

* E-mail: jvhomweg@pasteur.fr

## Introduction

Botulism is a severe neuroparalytic disease caused by botulinum toxin, and characterized by acute descending flaccid paralysis. The disease is caused by consumption of food contaminated with pre-formed toxin (foodborne botulism), or by absorption of toxin produced in *situ* in wounds (wound botulism) or colonized intestinal tracts (infant/intestinal adult botulism) [1]. Botulinum neurotoxins (BoNT) are the most potent toxins known and are

considered as one of the six highest risk threat agents of bioterrorism [2],[3].

BoNT are produced by six physiologically and genetically distinct bacteria, namely Clostridium botulinum Groups I to IV, and occasionally strains of C. butyricum and C. barati [4,5]. These neurotoxin-producing bacteria can be further categorized in seven major toxinotypes (A, B, C, D, E, F and G) based on the antigenic properties of the toxins they produce [6],[7]. Toxinotypes A, B, E and more rarely F are responsible for human botulism cases, while infections by C and D toxinotypes are observed mainly in animals.
Although most strains produce only one toxin, bivalent strains producing two different toxins (Ab, Af, Ba, and Bf) have also been reported [8–11]. In addition, *C. botulinum* strains producing a single toxin but carrying a silent gene for another (A(B)) [10], and strains producing a chimeric neurotoxin (C/D or D/C) have been described [12].

The genes encoding the neurotoxins (*bont*) can be found on the chromosome, a plasmid or a bacteriophage, and are located in a cluster with other genes encoding associated non-toxic proteins (ANTPs) forming the botulinum neurotoxin complexes [13]. The recent sequencing of the *bont* genes has revealed significant sequence variations within the toxinotypes and has led to the recognition of several subtypes, with differences ranging between 2–32% at the amino acid level [14]. In addition, the BoNT gene cluster varies in structure and organization [5], leading to the differentiation of two different conserved cluster types, termed the “ha cluster” and the “orf-X cluster” (Figure 1). The “ha cluster” consists of a set of hemaglutinin (HA) genes while the “orf-X cluster” consists of a set of genes of unknown function called “orf-Xs” associated to the *bont* gene. The diversity in nucleotide sequence and genetic cluster organization can be explained by the occurrence of many recombination and insertion events [15], as well as *bont* gene transfer between some toxinotypes of *C. botulinum*, and with non-botulinum species [7].

This variability not only has implications for antibody neutralization strategies [14], but also poses a great challenge for screening and characterization of BoNT-producing strains. Sensitive laboratory detection techniques are necessary for establishing monitoring programs to track *C. botulinum* contamination in animals, food and environmental samples. In addition, characterization methods are essential for tracing potentials sources of contamination in outbreak situations or forensic investigations, and for performing molecular risk assessments of BoNT-producing clostridia. Moreover, characterization of genetic differences allows the development of improved diagnostics and therapeutic agents for the treatment of botulism [14]. In consequence, efforts have been centered on the development of molecular techniques allowing simultaneous detection and typing of *C. botulinum* toxinotypes, such as multiplex real time PCR and focused DNA microarray [16,17][17–19]. These assays are highly sensitive but rely on hybridization of sequence-specific primers for

sufficient DNA amplification and/or on a limited number of probes for detection of target genomic regions, and may therefore be suboptimal when bacterial sequence divergence is high. In addition, although some multiplex assays allow simultaneous detection and differentiation of toxinotype A subtypes and strains [16,20], they do not allow typing further than the toxin type or in-depth genetic characterization of strains belonging to other known toxinotypes. Several genotyping methods such as pulse-field gel electrophoresis, randomly amplified polymorphic DNA, multiplex variable number tandem repeat analysis and comparative genomic hybridization microarray have been applied to further differentiate strains within a toxinotype. However, such methods are often technically demanding, difficult to standardize between laboratories and sometimes produce poorly interpretable results, as subtype and strain discrimination relies on differences in observed electrophoretic or hybridization patterns [18,19,21–24]. In addition, these approaches often require high amounts of genomic DNA obtained from pure bacterial cultures. The use of high-yield random DNA amplification methods, combined to sequence-based technologies targeting multiple molecular markers, is indispensable to overcome these limitations and obtain a better overall resolution.

To overcome the problem of limiting amount of available DNA in biological samples, whole genome amplification (WGA) methods have been developed and successfully applied for a number of genotyping assays [25–29]. In addition, high-density resequencing microarrays (RMA) have emerged as a rapid detection and molecular characterization tool for a broad range of bacterial and viral agents [30],[31]. This technology has demonstrated reliable detection of sequences differing up to 10–15% from the prototype sequences on the array, enhancing the spectrum of detection [31],[32]. In order to determine a target nucleotide sequence, resequencing microarrays use closely overlapping (“tiled”) probe sets of 25 mers, which contain one perfectly matched and three mismatched probes per base for both strands of the target genes [33]. In the present study, we evaluated the PathogenID v2.0 resequencing microarray, containing probes for the detection and characterization of neurotoxin-producing *Clostridium* species.

« ha - Cluster » : A1¹, A5, B, C², D², G³

<figure><img src="image_2.png" /></figure>

« OrfX - Cluster » : A1¹, A2, A3, A4, F, E⁴

<figure><img src="image_3.png" /></figure>

**Figure 1. Neurotoxin gene cluster organization.** In the gene cluster encoding the botulinum neurotoxin complex, the “ha cluster” appears to be associated with type A1, A5, B, C, D and G *bont* genes, and the “orf-X cluster” with type A1, A2, A3, A4, F and E *bont* genes.¹Although the type A1 *bont* gene is mostly found in a ha cluster, it has been found associated with an orf-X cluster in some single toxin strains and in all bivalent strains [18],[17].²The *botR* gene is in front of the ha genes in all toxinotype C and D strains.³The *Ha33* gene is absent in all toxinotype G strains.⁴The *botR* gene is absent in all toxinotype E strains
doi:10.1371/journal.pone.0067510.g001
## Results

### Design of PathogenID microarray for *Clostridium* detection and characterization

Detailed description of prototype sequences selected for *Clostridium* detection and characterization is shown in table 1. The array contains 4 tiled sequences targeting housekeeping genes of *Clostridium* species: the complete *rrs* gene of *C. botulinum* toxinotype A and partial *rpoB* genes of *C. difficile*, *C. perfringens* and *C. botulinum* toxinotype A. In addition, the RMA includes 48 tiled sequences ranging from 149 to 200 nt in length, targeting representative botulinum neurotoxin gene cluster components (*bont*, *ntnh*, *botR*, *ha17*, *ha70*, *ha33*, *orf-X1*, *orf-X2*, *orf-X3* and *p47*) of toxinotypes A through G, as well as three toxinotype A subtypes (A1, A2 and A3). Finally, it contains 19 tiled regions of 200 nt in length for the detection and characterization of other specific toxins produced by *C. tetani*, *C. difficile*, *C. perfringens*, *C. sordellii*, *C. septicum* and *C. oedematiens*.

### Resequencing of *Clostridium* with the PathogenID v2.0 microarray

The ability of the PathogenID v2.0 prototype sequences to detect and subtype neurotoxin-producing *C. botulinum* strains was assessed by performing a blind analysis of genomic DNA from eleven well-characterized strains representative of the main toxinotypes and subtypes. The observed call-rate for each tiled sequence is detailed in table 2. The call-rates marked in bold indicate sequences retained after filtering using the defined threshold. The average base call rate was 91.6% (range 52.5–100%) for strains for which genome sequence was identical to the tiled sequences. The analysis of more distant genomic sequences showed that these were successfully retrieved by RMA (Figure 2), and retained after filtering when they presented up to 24.7% divergence. The mean resequencing accuracy (correctly re-sequenced bases in comparison to reference sequencing results or database sequences) was 98.9% (range 87.1–100.0%). For all strains analysed by the RMA and the automated filtering approach, at least one sequence targeting a housekeeping gene and one sequence targeting a neurotoxin gene cluster component were successfully retrieved. For all bivalent strains analysed by RMA, two distinct *bont* gene prototype sequences could be retrieved: *bont*/B and *F* sequences for strain 168.08 (toxinotype *Bf2*), and *bont*/*A* and *B* sequences for strains NCTC 2916 [toxinotype *A1*(*B*)] and 1430-11 [toxinotype *A5*(*B*')]. Additionally, sequences of two complete sets of *ha* and *orf-X* cluster components were retrieved for strains 168.08 and NCTC 2916, but as expected only sequences of *ha* cluster components were recovered for strain 1430-11. Finally, the RMA did not only allow the retrieval of sequences specifically targeting BoNT gene cluster components *antp* genes from the toxinotype of the strain analysed, but also sequences targeting *antp* genes from other toxinotypes belonging to the same taxonomic group. In particular, toxinotype *B* prototype sequences (e.g. *botR*/*B*, *ha17*/*B*) were successfully retrieved for strain Hall (toxinotype *A*), although these sequences show up to 11% divergence with the corresponding regions of the Hall strain genome. Similarly, prototype sequences targeting toxinotype *A1* *antp* genes were retrieved for the BL6 strain (toxinotype *B*), and prototype sequences targeting toxinotypes *A2* and *A3* *antp* genes were retrieved for the NCIB10658 strain (toxinotype *F*). This was also observed for strains belonging to the *C. botulinum* taxonomic group III strains, as both toxinotype *C* and toxinotype *D* *antp* sequences were retrieved for strains 468 and 1873.

### Automated filtering and taxonomic identification of resequencing results

Automated filtering by TaxFinder using the defined parameters retained between 7 and 27 *C. botulinum* sequences for each strain (Table 2). Of the retained sequences, 98.8% (168/170) allowed retrieval of at least one BLAST hit, and 99,4% (167/168) of the retrieved BLAST hits lead to correct taxonomic identification at least to the toxinotype level. Moreover, blasting of all remaining unfiltered *C. botulinum* prototype sequences (n = 358) did not return any BLAST hit, indicating that no useful sequence was excluded. For all reference strains the systematic BLAST strategy of retained sequences obtained by the RMA allowed successful identification of the clostridial species and taxonomic group using housekeeping genes prototype sequences (Table 3). The strategy also permitted in-depth characterization of all tested strains using the BoNT gene cluster prototype sequences (Table 3). Moreover, BLAST analysis allowed simultaneous typing of the neurotoxin up to the subtype level, as for each strain a unique best hit or multiple best bacterial hits were retrieved and corresponded to the correct subtype. Finally, the turnaround time of this automated approach was inferior to 2 hours, extensively decreasing the time needed for filtering of sequences, blasting of retained sequences and result analysis.

### Identification, typing and characterization of a foodborne botulism outbreak strain

To demonstrate the capacity of the RMA approach to detect *C. botulinum* target sequences in the context of a foodborne outbreak, DNA extracts from naturally contaminated food samples were analysed. Food specimens were collected during the investigation of a *C. botulinum* family outbreak in Corsica [34]. RMA was performed on the DNA extract from a contaminated salad sample, as well as the enrichment culture of the salad sample. For DNA obtained directly from the salad specimen, a total of 286 bacterial prototype sequences were retained, including 3 sequences targeting *C. botulinum* genes (*rrs*, *rpoB* and *ha33*). However, none led to correct taxonomic identification of the *C. botulinum* infecting strain. For the DNA extract obtained from the enrichment culture of the contaminated salad, a total of 15 *C. botulinum* prototype sequences were retained after filtering, with call-rates ranging between 46.7 to 97.7%. In addition, sequences targeting partial *C. perfringens*-specific housekeeping genes (*rpoB*, *gyrA*, and *parC*) and toxin genes (*cpa*, *cpb2* and *pfoR*) were retrieved (call-rates ranging between 70.1 and 98.3%). All filtered sequences led to the retrieval of at least one BLAST hit and corresponding taxonomic identity. The highest scoring BLAST results obtained for *C.perfringens* prototype sequences retrieved by RMA all designated *C. perfringens* strains as taxonomic hit (data not shown). The results of the highest scoring BLAST alignment for each retained *C. botulinum* sequence (total score, coverage and taxonomic identity) are presented in figure 3. Strains belonging to the toxin subtype *A2* were recognized as best hit for most of the retrieved sequences of the neurotoxin gene cluster components. In addition, a gene cluster type *ha-*/orfX+ was identified. Phylogenetic analysis of the retrieved neurotoxin sequence with the corresponding sequences of five strains representative of the different toxin A subtypes (A1 to A5), demonstrated an overall identity of 100% with the neurotoxin sequence of the toxin subtype *A2* reference strain. The only mismatches in the multiple alignment were due to unidentified bases (Figure 4).
<table><caption>**Table 1.** Description of the sequences tiled on the PathogenID v2.0 microarray for the detection and characterization of BoNT-producing clostridia.</caption><thead><tr><th>Tiled sequence description</th><th>Species</th><th>Type</th><th>Strain</th><th>Genbank accession no.</th><th>Position (nt)</th><th>length (nt)</th></tr></thead><tbody><tr><td colspan="7"><strong>Housekeeping genes</strong></td></tr><tr><td>rrs</td><td>16S rRNA</td><td>C. botulinum</td><td>A1</td><td>DSM1734</td><td>X73442.1</td><td>7..1513</td><td>1506</td></tr><tr><td>rpoB</td><td>RNA polymerase (subunit β)</td><td>C. botulinum</td><td>A1</td><td>NCTC7272</td><td>Y16466.1</td><td>1249..1755</td><td>506</td></tr><tr><td>rpoB</td><td>RNA polymerase (subunit β)</td><td>C. perfringens</td><td></td><td>ATCC13124</td><td>CP000246.1</td><td>2986532..2986026</td><td>506</td></tr><tr><td>rpoB</td><td>RNA polymerase (subunit β)</td><td>C. difficile</td><td></td><td>630</td><td>AM180355.1</td><td>91744..92253</td><td>509</td></tr><tr><td colspan="7"><strong>Neurotoxin complex genes</strong></td></tr><tr><td>bont/A1</td><td>Botulinum neurotoxin type A1</td><td>C. botulinum</td><td>A1</td><td>Hall</td><td>AF461540.1</td><td>11250..11450</td><td>200</td></tr><tr><td>bont/A2</td><td>Botulinum neurotoxin type A2</td><td>C. botulinum</td><td>A2</td><td>Kyoto</td><td>X73423.1</td><td>1771..1971</td><td>200</td></tr><tr><td>bont/B1</td><td>Botulinum neurotoxin type B1</td><td>C. botulinum</td><td>B1</td><td>Okra</td><td>AB232927.1</td><td>14210..14410</td><td>200</td></tr><tr><td>bont/Bnp</td><td>Botulinum neurotoxin type Bnp</td><td>C. botulinum</td><td>Bnp</td><td>ATCC25765</td><td>X71343.1</td><td>3586..3786</td><td>200</td></tr><tr><td>bont/C1</td><td>Botulinum neurotoxin type C1</td><td>C. botulinum</td><td>C1</td><td>468</td><td>X53751.1</td><td>1876..2076</td><td>200</td></tr><tr><td>bont/D</td><td>Botulinum neurotoxin type D</td><td>C. botulinum</td><td>D</td><td>1873</td><td>X54254.1</td><td>1697..1897</td><td>200</td></tr><tr><td>bont/E</td><td>Botulinum neurotoxin type E</td><td>C. botulinum</td><td>E1</td><td>Beluga</td><td>X62089.1</td><td>1213..1413</td><td>200</td></tr><tr><td>bont/F</td><td>Botulinum neurotoxin type F</td><td>C. botulinum</td><td>F6</td><td>202</td><td>M92906.1</td><td>606..806</td><td>200</td></tr><tr><td>bont/G</td><td>Botulinum neurotoxin type G</td><td>C. botulinum</td><td>G</td><td>NCFB3012</td><td>X74162.1</td><td>1651..1851</td><td>200</td></tr><tr><td>ntnh/A1</td><td>Non-toxigenic non-hemagglutinin</td><td>C. botulinum</td><td>A1</td><td>Hall</td><td>AF461540</td><td>8501..8695</td><td>194</td></tr><tr><td>ntnh/A2</td><td>Non-toxigenic non-hemagglutinin</td><td>C. botulinum</td><td>A2</td><td>Kyoto</td><td>X87974</td><td>2453..2647</td><td>194</td></tr><tr><td>ntnh/A3</td><td>Non-toxigenic non-hemagglutinin</td><td>C. botulinum</td><td>A3</td><td>Mascarpone</td><td>DQ310546</td><td>5690..5884</td><td>194</td></tr><tr><td>ntnh/Bp</td><td>Non-toxigenic non-hemagglutinin</td><td>C. botulinum</td><td>B1</td><td>Okra</td><td>AB232927</td><td>9649..9843</td><td>194</td></tr><tr><td>ntnh/C</td><td>Non-toxigenic non-hemagglutinin</td><td>C. botulinum</td><td>C1</td><td>468</td><td>X72793</td><td>4757..4966</td><td>209</td></tr><tr><td>ntnh/D</td><td>Non-toxigenic non-hemagglutinin</td><td>C. botulinum</td><td>D</td><td>1873</td><td>AB012112</td><td>4483..4692</td><td>209</td></tr><tr><td>ntnh/E</td><td>Non-toxigenic non-hemagglutinin</td><td>C. botulinum</td><td>E</td><td>Mashike</td><td>D12697</td><td>2644..2838</td><td>194</td></tr><tr><td>ntnh/Ebut</td><td>Non-toxigenic non-hemagglutinin</td><td>C. butyricum</td><td>E4</td><td>BL6340</td><td>D12739</td><td>2644..2838</td><td>194</td></tr><tr><td>ntnh/F</td><td>Non-toxigenic non-hemagglutinin</td><td>C. botulinum</td><td>F6</td><td>202</td><td>S73676</td><td>2664..2858</td><td>194</td></tr><tr><td>ntnh/G</td><td>Non-toxigenic non-hemagglutinin</td><td>C. botulinum</td><td>G</td><td>ATCC27322</td><td>X87972</td><td>3390..3599</td><td>209</td></tr><tr><td>ha17/A</td><td>Hemagglutinin 17</td><td>C. botulinum</td><td>A1</td><td>Hall</td><td>AF461540</td><td>4022..3822</td><td>200</td></tr><tr><td>ha17/Bp</td><td>Hemagglutinin 17</td><td>C. botulinum</td><td>B1</td><td>Okra</td><td>AB232927</td><td>5170..4970</td><td>200</td></tr><tr><td>ha17/C</td><td>Hemagglutinin 17</td><td>C. botulinum</td><td>C1</td><td>468</td><td>X72793</td><td>3181..2981</td><td>200</td></tr><tr><td>ha17/D</td><td>Hemagglutinin 17</td><td>C. botulinum</td><td>D</td><td>1873</td><td>AB012112</td><td>2907..2707</td><td>200</td></tr><tr><td>ha17/G</td><td>Hemagglutinin 17</td><td>C. botulinum</td><td>G</td><td>ATCC27322</td><td>X87972</td><td>1749..1600</td><td>149</td></tr><tr><td>ha70/A</td><td>Hemagglutinin 70</td><td>C. botulinum</td><td>A</td><td>Hall</td><td>AF461540</td><td>3166..3005</td><td>161</td></tr><tr><td>ha70/Bp</td><td>Hemagglutinin 70</td><td>C. botulinum</td><td>B1</td><td>Okra</td><td>AB232927</td><td>4314..4153</td><td>161</td></tr><tr><td>ha70/C/D</td><td>Hemagglutinin 70</td><td>C. botulinum</td><td>C1</td><td>468</td><td>X72793</td><td>2337..2173</td><td>164</td></tr><tr><td>ha70/G</td><td>Hemagglutinin 70</td><td>C. botulinum</td><td>G</td><td>ATCC27322</td><td>X87972</td><td>945..787</td><td>158</td></tr><tr><td>ha33/A</td><td>Hemagglutinin 33</td><td>C. botulinum</td><td>A1</td><td>Hall</td><td>AF461540</td><td>4603..4403</td><td>200</td></tr><tr><td>ha33/Bp</td><td>Hemagglutinin 33</td><td>C. botulinum</td><td>B1</td><td>Okra</td><td>AB232927</td><td>5748..5548</td><td>200</td></tr><tr><td>ha33/C/D</td><td>Hemagglutinin 33</td><td>C. botulinum</td><td>C1</td><td>468</td><td>X72793</td><td>3750..3550</td><td>200</td></tr><tr><td>orfX1/A2</td><td>Orf-X component 1</td><td>C. botulinum</td><td>A2</td><td>Kyoto</td><td>AB004778</td><td>524..315</td><td>209</td></tr><tr><td>orfX1/A3</td><td>Orf-X component 1</td><td>C. botulinum</td><td>A3</td><td>Mascarpone</td><td>DQ310546</td><td>1001..792</td><td>209</td></tr><tr><td>orfX1/E</td><td>Orf-X component 1</td><td>C. botulinum</td><td>E</td><td>Iwanai</td><td>D88418</td><td>1046..837</td><td>209</td></tr><tr><td>orfX2/A2</td><td>Orf-X component 2</td><td>C. botulinum</td><td>A2</td><td>Kyoto</td><td>AY497358</td><td>4619..4373</td><td>246</td></tr><tr><td>orfX2/E</td><td>Orf-X component 2</td><td>C. botulinum</td><td>E</td><td>Iwanai</td><td>D88418</td><td>256..10</td><td>246</td></tr><tr><td>orfX3/A2</td><td>Orf-X component 3</td><td>C. botulinum</td><td>A2</td><td>Kyoto</td><td>AY497358</td><td>2445..2236</td><td>209</td></tr><tr><td>botR/A1</td><td>Transcriptionnal regulator BotR</td><td>C. botulinum</td><td>A1</td><td>Hall</td><td>AF461540</td><td>5502..5702</td><td>200</td></tr><tr><td>botR/A2</td><td>Transcriptionnal regulator BotR</td><td>C. botulinum</td><td>A2</td><td>Kyoto</td><td>X96493</td><td>229..29</td><td>200</td></tr><tr><td>botR/A3</td><td>Transcriptionnal regulator BotR</td><td>C. botulinum</td><td>A3</td><td>Mascarpone</td><td>DQ310546</td><td>1438..1238</td><td>200</td></tr><tr><td>botR/B</td><td>Transcriptionnal regulator BotR</td><td>C. botulinum</td><td>B1</td><td>Okra</td><td>AB232927</td><td>6653..6853</td><td>200</td></tr><tr><td>botR/C</td><td>Transcriptionnal regulator BotR</td><td>C. botulinum</td><td>C1</td><td>468</td><td>X72793</td><td>512..712</td><td>200</td></tr></tbody></table>
<table><caption>**Table 1.** Cont.</caption><thead><tr><th>Tiled sequence description</th><th>Species</th><th>Type</th><th>Strain</th><th>Genbank accession no.</th><th>Position (nt)</th><th>length (nt)</th></tr></thead><tbody><tr><td>botR/D</td><td>Transcriptionnal regulator BotR</td><td>C. botulinum</td><td>D</td><td>AB012112</td><td>238..438</td><td>200</td></tr><tr><td>botR/G</td><td>Transcriptionnal regulator BotR</td><td>C. botulinum</td><td>G</td><td>ATCC27322</td><td>X87972</td><td>2595..2795</td><td>200</td></tr><tr><td>p47/A2</td><td>P 47</td><td>C. botulinum</td><td>A2</td><td>Kyoto</td><td>X96493</td><td>779..937</td><td>158</td></tr><tr><td>p47/A3</td><td>P 47</td><td>C. botulinum</td><td>A3</td><td>Mascarpone</td><td>DQ310546</td><td>1976..2134</td><td>158</td></tr><tr><td>p47/E</td><td>P 47</td><td>C. botulinum</td><td>E</td><td>Iwanai</td><td>D88418</td><td>1881..2081</td><td>200</td></tr><tr><td>p47/F</td><td>P 47</td><td>C. botulinum</td><td>F6</td><td>202</td><td>Y10770</td><td>946..1116</td><td>170</td></tr><tr><td colspan="8"><strong>Clostridial toxin genes</strong></td></tr><tr><td>C2</td><td>C2 toxin (component 1)</td><td>C. botulinum</td><td></td><td>(C)-203U28</td><td>D88982.1</td><td>838..1038</td><td>200</td></tr><tr><td>C2</td><td>C2 toxin (component 2)</td><td>C. botulinum</td><td></td><td>(C)-203U28</td><td>D88982.1</td><td>2245..2445</td><td>200</td></tr><tr><td>tent</td><td>Tetanus neurotoxin</td><td>C. tetani</td><td></td><td>Massachusetts</td><td>X04436.1</td><td>2026..2226</td><td>200</td></tr><tr><td>toxB</td><td>Toxin B</td><td>C. difficile</td><td></td><td>VPI10463</td><td>X53138.1</td><td>1518..1718</td><td>200</td></tr><tr><td>toxB</td><td>Toxin BF</td><td>C. difficile</td><td></td><td>1470</td><td>Z23277.1</td><td>1183..1383</td><td>200</td></tr><tr><td>toxA</td><td>Toxin A</td><td>C. difficile</td><td></td><td>VPI10463</td><td>X51797.1</td><td>2521..2721</td><td>200</td></tr><tr><td>cpa</td><td>α-toxin</td><td>C. perfringens</td><td></td><td>S13</td><td>L43546.1</td><td>1317..1517</td><td>200</td></tr><tr><td>cpb1</td><td>β-toxin</td><td>C. perfringens</td><td></td><td>NCTC 8533</td><td>L13198.1</td><td>510..710</td><td>200</td></tr><tr><td>cpb2</td><td>β2-toxin</td><td>C. perfringens</td><td>Porcine</td><td>CWC245</td><td>L77965.1</td><td>418..618</td><td>200</td></tr><tr><td>cpb2</td><td>β2-toxin</td><td>C. perfringens</td><td>Equine</td><td>D21/98</td><td>AJ537535.1</td><td>199..399</td><td>200</td></tr><tr><td>cpb2</td><td>β2-toxin</td><td>C. perfringens</td><td>Bovine</td><td>JGS4147</td><td>AY609175</td><td>151..351</td><td>200</td></tr><tr><td>etxD</td><td>ε-toxin D</td><td>C. perfringens</td><td></td><td>NCTC 8346</td><td>M95206.1</td><td>714..914</td><td>200</td></tr><tr><td>pfoR</td><td>τ-toxin (perfringolysin)</td><td>C. perfringens</td><td></td><td>lambda gt10</td><td>M81080.1</td><td>3529..3729</td><td>200</td></tr><tr><td>enterotoxin</td><td>Enterotoxin</td><td>C. perfringens</td><td></td><td>NCTC8239</td><td>M98037.1</td><td>787..987</td><td>200</td></tr><tr><td>toxin-iota-a</td><td>τ-toxin (component Ia)</td><td>C. perfringens</td><td></td><td>NCIB10748</td><td>X73562.1</td><td>2254..2454</td><td>200</td></tr><tr><td>toxin-iota-b</td><td>τ-toxin (component Ib)</td><td>C. perfringens</td><td></td><td>NCIB10748</td><td>X73562.1</td><td>3515..3715</td><td>200</td></tr><tr><td>toxin-LT82</td><td>Lethal toxin 82</td><td>C. sordellii</td><td></td><td>82</td><td>X82638.1</td><td>1595..1795</td><td>200</td></tr><tr><td>α-toxin</td><td>α-toxin</td><td>C. septicum</td><td></td><td>NCTC547</td><td>D17668.1</td><td>1062..1262</td><td>200</td></tr><tr><td>α-toxin</td><td>α-novyi toxin</td><td>C. oedematiens</td><td></td><td>ATCC19402</td><td>Z48636.1</td><td>1286..1486</td><td>200</td></tr><tr><td colspan="8"><strong>Quinolone-resistance determining regions</strong></td></tr><tr><td>gyrA-qrdr</td><td>Gyrase (subunit A)</td><td>C. perfringens</td><td></td><td>ATCC13124</td><td>CP000246.1</td><td>7269..7466</td><td>197</td></tr><tr><td>gyrA-qrdr</td><td>Gyrase (subunit A)</td><td>C. difficile</td><td></td><td>630</td><td>AM180355.1</td><td>6231..6428</td><td>197</td></tr><tr><td>parC-qrdr</td><td>Topoisomerase (subunit IV)</td><td>C. perfringens</td><td></td><td>SM101</td><td>CP000312.1</td><td>2254972..2255160</td><td>188</td></tr></tbody></table>

doi:10.1371/journal.pone.0067510.t001

Discussion

Although botulism is quite rare in developed countries, recent outbreaks have demonstrated the potential severity of this hazard [35–40]. In such a context, rapid identification and characterization of BoNT-producing clostridia in home-made or manufactured food is a priority to stop the spread of botulism infections and trace back the source of contamination [41–43]. While identification of BoNT-encoding genes, using rapid nucleic-acid detection tools, is usually satisfactory for molecular risk assessment of BoNT-producing clostridia in food and environmental samples, it is not sufficient to discriminate strains within a toxinotype for source-tracing and molecular epidemiology [16]. For this reason, recent recommendations have stated that new assays should be able to detect variants for all toxinotypes, should be type-specific to determine proper treatment, and should be sensitive to perform risk assessment (NIAID Expert Panel on Botulism Diagnostics, Bethesda Maryland, May 2003). The development of such universal genetic characterization tools is hampered by the great variability of the genomic background, the sequence of the *bont*

gene and the arrangement of the BoNT gene cluster of *C. botulinum* strains, and therefore only available for a limited number of strains.

In the present study, we designed and evaluated the performances of a resequencing microarray combined with an automated data approach for the simultaneous detection and characterization of BoNT-producing clostridia. The described RMA, designated PathogenID v2.0, and similar arrays have been proven successful for the characterization of numerous bacteria including *Neisseria meningitidis*, *Streptococcus pyogenes*, *Streptococcus pneumoniae*, *Bacillus anthracis*, and several species of the *Enterobacteriaceae* family [30,44–46]. For the purpose of this study, prototype sequences of housekeeping genes *rrs* and *rpoB* from clostridia pathogenic for humans and most frequently associated with outbreaks were incorporated in the RMA for identification of clostridial species and lineages [47–49]. In addition, sequences from the BoNT gene cluster were included for their involvement in toxicity and usefulness for studying genetic variation [18,50,51].

The method was validated with DNA extracts from a representative panel of reference *C. botulinum* strains. The obtained
<figure><img src="image_4.png" /><figcaption>**Figure 2. Detection spectrum of the PathogenID_v2.0 according to the nucleotide diversity of the tested BoNT-producing strains.**
The spectrum of detection of BoNT-producing clostridia by the PathogenID v2.0 was assessed using 9 well-characterized strains representative of all toxinotypes and several toxin A subtypes. Each sequence obtained by the RMA for which the percentage of nucleotide divergence compared to the corresponding tiled sequence is known (n = 118), is indicated by a blue diamond, and presented according to the percentage of nucleotide bases determined by the RMA (call rate). The linear association between these two parameters is shown, and demonstrates a good correlation (correlation coefficient R² of 0.79).
doi:10.1371/journal.pone.0067510.g002</figcaption></figure>

results demonstrated that this approach has the potential to detect neurotoxin-producing *C. botulinum* strains, but also phenotypically related strains, by providing high quality sequence data. The RMA approach reliably detected and characterized the neurotoxin in complex genetic determinants for all tested toxino-(sub)types, including bivalent strains. Further analysis of resequencing results showed that targeted sequences were successfully retrieved by the RMA when they presented up to 24.7% divergence with the tiled probes. This data supports the broadness of detection of the RMA, as the divergence of the tiled probes compared to the respective *bont* gene regions from all known toxin A, B, C, D, E and G subtypes, and most toxin F subtypes, was calculated to be inferior to 21%. Although the tiled *bont/F* sequence demonstrated respectively 70.65% and 55.72% with the corresponding gene region of the rare toxin subtypes F5 and F7, the high homology of the additional housekeeping gene and *antp* gene sequences targeted by the RMA should allow detection and correct identification of these strains.

In addition to sequence divergence, differences in DNA extract quality, genomic copy number variation, location of targeted genomic regions and possible secondary structures, as well as the G and C content of the tiled probes and the presence of nonspecific nucleotide stretches could have affected genome amplification and hybridization efficacy, possibly causing reduced base-call rates. Nevertheless, this did not hamper the correct identification and characterization of the analyzed strains.

Subsequently, the pathogenID v2.0 was used to identify *C. botulinum* strains in contaminated food specimens. The analysis of

the DNA extract obtained directly from these specimens by the RMA did not permit quality resequencing of the infecting *C. botulinum* strain sequences. These results are due to the presence of numerous cultivable and non-cultivable environmental bacteria naturally present in food, competing with the *C. botulinum* strain during the whole genome amplification step, and subsequently increasing the noise level observed at the time of the analysis. Nevertheless, a 48 h enrichment culture of the contaminated food samples was sufficient to successfully detect and characterize the infecting *C. botulinum* strain as well as a co-cultured *C. perfringens* strain. The results obtained by BLASTing and phylogenetic analysis of the retrieved *bont* sequence, as well as the characterization of the BoNT gene cluster components, classified the strain as a *C. botulinum* toxinotype A2. Moreover, the examination of the taxonomic BLAST hits indicated a close relatedness to strain Mascarpone. As the analysis was performed with a partial sequence of the *bont* gene (200 nt), containing 4% unidentified bases, misclassification of the toxin subtype may have occurred. However, this was unlikely as the retained sequence contained specific determinants found only in *bont/A2* genes, and unidentified residues were located in regions that are not discriminant for subtype differentiation. The *C. botulinum* strain detected by RMA could be isolated from the food sample, and confirmed as subtype A2 by sequencing the specifically PCR-amplified *bont* gene (genbank accession number JQ954970).

The many advantages of using DNA microarrays for characterization of bacteria have been discussed in the past and can be applied to this assay. First, the genetic characterization results can
<table><caption>**Table 2.** Call-rates observed for the different *C. botulinum* sequences retrieved by the PathogenID V2.0 microarray for all tested bacterial strains.</caption><thead><tr><th rowspan="2">Tiled C. botulinum sequence</th><th colspan="10">Clostridium botulinum toxin subtype and strain</th></tr><tr><th>A1</th><th>A5(B')</th><th>A1(B)</th><th>A3</th><th>B2</th><th>Bf2</th><th>F6</th><th>E1</th><th>C1</th><th>D</th><th>G</th></tr><tr><td></td><td>Hall</td><td>1430.11</td><td>2916</td><td>Loch Maree</td><td>BL6</td><td>168.08</td><td>NCIB 10658</td><td>K16</td><td>468</td><td>1873</td><td>NCIB 10714</td></tr></thead><tbody><tr><td colspan="12"><strong>Housekeeping genes</strong></td></tr><tr><td>rrs A1</td><td>88,1</td><td>83,5</td><td>83,9</td><td>80,6</td><td>90,7</td><td>21,5</td><td>85,1</td><td>69,8</td><td>71,7</td><td>61,7</td><td>69,2</td></tr><tr><td>RpoB A1</td><td>96,5</td><td>96,3</td><td>94,6</td><td>91,1</td><td>90,7</td><td>90,7</td><td>95,2</td><td>16,4</td><td>23,6</td><td>15,7</td><td>16,4</td></tr><tr><td colspan="12"><strong>Neurotoxin genes</strong></td></tr><tr><td>bont/A1</td><td>97,2</td><td>87,0</td><td>97,7</td><td>23,7</td><td>6,2</td><td>8,5</td><td>11,9</td><td>7,3</td><td>16,4</td><td>11,3</td><td>7,9</td></tr><tr><td>bont/A2</td><td>18,1</td><td>12,4</td><td>21,5</td><td>88,7</td><td>8,5</td><td>9,0</td><td>16,4</td><td>7,9</td><td>7,9</td><td>2,3</td><td>7,9</td></tr><tr><td>bont/B1</td><td>14,7</td><td>18,1</td><td>78,5</td><td>11,3</td><td>19,8</td><td>88,7</td><td>11,3</td><td>10,2</td><td>13,6</td><td>1,7</td><td>10,7</td></tr><tr><td>bont/Bnp</td><td>10,2</td><td>27,1</td><td>35,6</td><td>10,7</td><td>20,3</td><td>35,6</td><td>10,2</td><td>3,4</td><td>8,5</td><td>2,8</td><td>11,3</td></tr><tr><td>bont/C1</td><td>14,1</td><td>8,5</td><td>10,7</td><td>14,1</td><td>4,5</td><td>7,3</td><td>16,9</td><td>10,7</td><td>99,4</td><td>1,1</td><td>11,9</td></tr><tr><td>bont/D</td><td>13,6</td><td>2,8</td><td>12,4</td><td>21,5</td><td>4,5</td><td>5,7</td><td>19,2</td><td>15,3</td><td>13,6</td><td>84,7</td><td>12,4</td></tr><tr><td>bont/E</td><td>6,8</td><td>1,7</td><td>2,8</td><td>9,6</td><td>4,5</td><td>1,1</td><td>6,2</td><td>67,2</td><td>5,6</td><td>0,6</td><td>2,3</td></tr><tr><td>bont/F</td><td>16,4</td><td>4,0</td><td>16,9</td><td>23,7</td><td>11,9</td><td>49,2</td><td>53,1</td><td>11,3</td><td>13,6</td><td>3,4</td><td>8,5</td></tr><tr><td>bont/G</td><td>11,9</td><td>7,9</td><td>11,9</td><td>11,3</td><td>5,6</td><td>9,0</td><td>11,3</td><td>8,5</td><td>10,2</td><td>3,4</td><td>76,8</td></tr><tr><td colspan="12"><strong>Non-toxin non-hemagglutinin gene</strong></td></tr><tr><td>ntnh/A1</td><td>99,4</td><td>81,9</td><td>90,6</td><td>75,4</td><td>10,5</td><td>78,4</td><td>55,6</td><td>9,4</td><td>8,2</td><td>1,2</td><td>6,4</td></tr><tr><td>ntnh/A2</td><td>86,0</td><td>78,4</td><td>87,7</td><td>60,8</td><td>15,8</td><td>76,6</td><td>80,7</td><td>13,5</td><td>7,6</td><td>1,2</td><td>8,2</td></tr><tr><td>ntnh/A3</td><td>70,2</td><td>73,7</td><td>81,3</td><td>69,6</td><td>8,8</td><td>81,3</td><td>48,5</td><td>8,2</td><td>7,0</td><td>0,6</td><td>7,0</td></tr><tr><td>ntnh/Bp</td><td>11,1</td><td>2,3</td><td>14,6</td><td>12,3</td><td>97,1</td><td>98,2</td><td>19,9</td><td>8,8</td><td>20,5</td><td>4,1</td><td>29,2</td></tr><tr><td>ntnh/C</td><td>26,9</td><td>11,8</td><td>17,2</td><td>17,7</td><td>16,1</td><td>9,7</td><td>13,4</td><td>14,5</td><td>96,2</td><td>72,0</td><td>9,7</td></tr><tr><td>ntnh/D</td><td>32,8</td><td>21,0</td><td>24,2</td><td>19,4</td><td>22,0</td><td>17,2</td><td>13,4</td><td>14,5</td><td>96,2</td><td>79,6</td><td>9,1</td></tr><tr><td>ntnh/E</td><td>21,1</td><td>12,3</td><td>24,0</td><td>26,3</td><td>7,6</td><td>9,9</td><td>15,8</td><td>87,1</td><td>15,8</td><td>2,3</td><td>11,1</td></tr><tr><td>ntnh/Ebut</td><td>21,6</td><td>11,7</td><td>24,0</td><td>26,3</td><td>6,4</td><td>9,9</td><td>15,2</td><td>86,5</td><td>15,2</td><td>1,8</td><td>11,1</td></tr><tr><td>ntnh/F</td><td>45,0</td><td>38,6</td><td>39,8</td><td>39,8</td><td>9,4</td><td>43,3</td><td>33,3</td><td>22,8</td><td>14,0</td><td>1,8</td><td>9,4</td></tr><tr><td>ntnh/G</td><td>27,4</td><td>14,5</td><td>22,6</td><td>22,0</td><td>16,1</td><td>15,1</td><td>14,0</td><td>7,0</td><td>24,2</td><td>5,9</td><td>64,5</td></tr><tr><td colspan="12"><strong>ha-cluster components</strong></td></tr><tr><td>ha70/A</td><td>94,9</td><td>63,8</td><td>91,3</td><td>13,0</td><td>76,8</td><td>96,4</td><td>9,4</td><td>3,6</td><td>9,4</td><td>0,0</td><td>7,2</td></tr><tr><td>ha70/Bp</td><td>88,4</td><td>77,5</td><td>85,5</td><td>13,8</td><td>84,8</td><td>92,0</td><td>6,5</td><td>3,6</td><td>8,0</td><td>0,7</td><td>8,0</td></tr><tr><td>ha70/C/D</td><td>7,8</td><td>4,3</td><td>7,1</td><td>13,5</td><td>5,7</td><td>7,1</td><td>5,7</td><td>5,7</td><td>100,0</td><td>79,4</td><td>6,4</td></tr><tr><td>ha70/G</td><td>14,8</td><td>23,0</td><td>15,6</td><td>10,4</td><td>23,0</td><td>11,9</td><td>12,6</td><td>5,2</td><td>25,2</td><td>5,2</td><td>45,9</td></tr><tr><td>ha17/A</td><td>98,9</td><td>80,8</td><td>97,2</td><td>15,8</td><td>86,4</td><td>96,6</td><td>15,8</td><td>11,3</td><td>19,8</td><td>8,5</td><td>11,9</td></tr><tr><td>ha17/Bp</td><td>84,7</td><td>84,2</td><td>82,5</td><td>15,8</td><td>84,2</td><td>84,2</td><td>14,7</td><td>13,0</td><td>19,2</td><td>9,6</td><td>13,0</td></tr><tr><td>ha17/C</td><td>10,7</td><td>7,9</td><td>5,6</td><td>18,1</td><td>8,5</td><td>7,3</td><td>10,2</td><td>7,3</td><td>100,0</td><td>89,3</td><td>7,9</td></tr><tr><td>ha17/D</td><td>11,3</td><td>8,5</td><td>4,5</td><td>16,4</td><td>8,5</td><td>6,8</td><td>7,9</td><td>6,2</td><td>100,0</td><td>84,7</td><td>9,0</td></tr><tr><td>ha17/G</td><td>11,1</td><td>3,2</td><td>4,8</td><td>9,5</td><td>3,2</td><td>5,6</td><td>5,6</td><td>6,3</td><td>9,5</td><td>0,0</td><td>57,9</td></tr><tr><td>ha33/A</td><td>96,0</td><td>45,2</td><td>78,5</td><td>14,1</td><td>49,2</td><td>85,3</td><td>8,5</td><td>6,2</td><td>13,6</td><td>4,0</td><td>9,0</td></tr><tr><td>ha33/Bp</td><td>30,5</td><td>30,5</td><td>33,3</td><td>7,3</td><td>55,4</td><td>41,8</td><td>6,8</td><td>8,5</td><td>6,8</td><td>2,3</td><td>4,5</td></tr><tr><td>ha33/C/D</td><td>19,2</td><td>5,1</td><td>15,3</td><td>10,7</td><td>9,0</td><td>4,0</td><td>10,7</td><td>7,9</td><td>99,4</td><td>0,6</td><td>9,0</td></tr><tr><td>botR/A1</td><td>93,2</td><td>59,3</td><td>74,6</td><td>18,1</td><td>83,6</td><td>89,3</td><td>13,6</td><td>5,1</td><td>16,4</td><td>0,6</td><td>10,7</td></tr><tr><td>botR/Bp</td><td>81,9</td><td>68,4</td><td>88,1</td><td>20,3</td><td>67,2</td><td>93,2</td><td>13,0</td><td>7,3</td><td>17,5</td><td>0,6</td><td>10,7</td></tr><tr><td>botR/C</td><td>7,9</td><td>5,1</td><td>14,1</td><td>11,3</td><td>9,6</td><td>13,6</td><td>14,1</td><td>7,3</td><td>98,9</td><td>45,2</td><td>7,3</td></tr><tr><td>botR/D</td><td>9,0</td><td>6,8</td><td>13,6</td><td>11,3</td><td>8,5</td><td>11,9</td><td>15,3</td><td>7,3</td><td>98,9</td><td>52,5</td><td>8,5</td></tr><tr><td>botR/G</td><td>15,8</td><td>4,5</td><td>11,9</td><td>11,3</td><td>9,0</td><td>2,3</td><td>7,9</td><td>4,5</td><td>13,0</td><td>0,6</td><td>46,9</td></tr><tr><td colspan="12"><strong>orfX-cluster components</strong></td></tr><tr><td>orfX1-A2</td><td>5,4</td><td>1,1</td><td>57,5</td><td>90,3</td><td>6,5</td><td>59,7</td><td>73,7</td><td>3,8</td><td>7,5</td><td>3,8</td><td>6,5</td></tr><tr><td>orfX2-A2</td><td>14,8</td><td>4,5</td><td>90,1</td><td>88,8</td><td>9,0</td><td>68,2</td><td>83,9</td><td>11,2</td><td>10,8</td><td>6,7</td><td>12,6</td></tr></tbody></table>
<table><caption>**Table 2.** Cont.</caption><thead><tr><th rowspan="3">Tiled C. botulinum sequence</th><th colspan="10">Clostridium botulinum toxin subtype and strain</th></tr><tr><th>A1</th><th>A5(B')</th><th>A1(B)</th><th>A3</th><th>B2</th><th>Bf2</th><th>F6</th><th>E1</th><th>C1</th><th>D</th><th>G</th></tr><tr><th>Hall</th><th>1430.11</th><th>2916</th><th>Loch Maree</th><th>BL6</th><th>168.08</th><th>NCIB 10658</th><th>K16</th><th>468</th><th>1873</th><th>NCIB 10714</th></tr></thead><tbody><tr><td>orfX3-A2</td><td>9,7</td><td>2,2</td><td><strong>87,1</strong></td><td><strong>97,8</strong></td><td>10,8</td><td><strong>96,8</strong></td><td><strong>86,6</strong></td><td>14,0</td><td>20,4</td><td>7,5</td><td>9,7</td></tr><tr><td>orfX1-A3</td><td>8,6</td><td>5,4</td><td><strong>78,0</strong></td><td><strong>61,3</strong></td><td>7,5</td><td><strong>41,9</strong></td><td><strong>54,3</strong></td><td>5,4</td><td>7,5</td><td>1,6</td><td>5,9</td></tr><tr><td>orfX1-E</td><td>6,5</td><td>5,9</td><td>9,1</td><td>11,3</td><td>3,2</td><td>9,1</td><td>7,5</td><td><strong>68,3</strong></td><td>8,1</td><td>7,5</td><td>6,5</td></tr><tr><td>orfX2-E</td><td>9,9</td><td>5,4</td><td>15,2</td><td>22,4</td><td>6,3</td><td>7,2</td><td>17,0</td><td><strong>71,3</strong></td><td>10,3</td><td>1,8</td><td>7,2</td></tr><tr><td>botR/A2</td><td>9,6</td><td>1,1</td><td><strong>39,0</strong></td><td><strong>74,0</strong></td><td>8,5</td><td><strong>62,1</strong></td><td><strong>42,9</strong></td><td>2,8</td><td>10,2</td><td>2,3</td><td>4,5</td></tr><tr><td>botR/A3</td><td>8,5</td><td>0,6</td><td><strong>68,4</strong></td><td><strong>41,8</strong></td><td>2,8</td><td><strong>69,5</strong></td><td><strong>39,5</strong></td><td>5,1</td><td>10,2</td><td>1,7</td><td>6,2</td></tr><tr><td>p47/A2</td><td>9,6</td><td>3,7</td><td><strong>40,7</strong></td><td><strong>31,1</strong></td><td>6,7</td><td>36,3</td><td>36,3</td><td>6,7</td><td>14,1</td><td>1,5</td><td>8,1</td></tr><tr><td>p47/A3</td><td>10,4</td><td>1,5</td><td><strong>90,4</strong></td><td><strong>61,5</strong></td><td>6,7</td><td><strong>69,6</strong></td><td><strong>72,6</strong></td><td>5,9</td><td>15,6</td><td>1,5</td><td>7,4</td></tr><tr><td>p47/E</td><td>5,1</td><td>1,1</td><td><strong>21,5</strong></td><td><strong>33,3</strong></td><td>2,3</td><td>23,2</td><td><strong>19,8</strong></td><td><strong>80,8</strong></td><td>4,0</td><td>0,0</td><td>2,3</td></tr><tr><td>p47/F</td><td>8,8</td><td>0,7</td><td>8,8</td><td>6,1</td><td>4,8</td><td><strong>95,9</strong></td><td>19,0</td><td>6,1</td><td>8,8</td><td>1,4</td><td>12,9</td></tr></tbody></table>

Sequences retained after filtering using the defined threshold are designated in bold.
doi:10.1371/journal.pone.0067510.t002

be used to support epidemiological associations in outbreaks involving a large number of samples or those from multiple geographical locations[17]. Secondly, they can be used to trace or confirm the source of an outbreak or to perform environmental risk assessments [52,53]. In the specific context of botulism, the results provided by DNA microarray can further serve to design specific primers for rapid nucleic acid detection methods, and allow a preliminary differentiation of the strains before applying more laborious genotyping methods such as PFGE [17].

In addition to these advantages, the specific design and properties of high-density RMAs allow the detection of a broad range of strains, including potentially emerging variants [44,54]. Firstly, the use of closely overlapping 25 nt probe sets in the RMA approach instead of unique probes to cover a target region of 100–500 nt increases the chance for hybridization of divergent target DNA to the array. Indeed, studies using focused or comparative

genomic hybridization microarrays have reported a maximum value of 15–18% nucleotide sequence divergence of probe regions for positive microarray detection of *C. botulinum* genomic DNA, as compared to 24.7% for the RMA described here [17–19,24]. In addition, the relatedness between strains can be quantified more precisely by using a sequence-based approach such as the RMA, rather than by analysis of differences in hybridization patterns provided by classic microarrays. This is particularly important in the case of botulism, as the frequent occurrence of gene recombination, insertion and transfer observed between *Clostridium* strains indicate that such events are likely to occur again. In such a case, the sequence information provided by the RMA will be most valuable for the characterization of the emerging variants and the development of improved rapid detection tools. This also emphasizes the need for bacterial culture isolates as well as well-

<table><caption>**Table 3.** Identification and characterization of *C. botulinum* strains based on the sequences retrieved by the PathogenID v2.0 microarray and the BLAST analysis results.</caption><thead><tr><th rowspan="2">Strain</th><th rowspan="2">Toxin subtype</th><th colspan="3">PathogenID v2.0 RMA results</th></tr><tr><th>Taxonomic identification</th><th>Neurotoxin subtype</th><th>BoNT gene cluster characterization</th></tr></thead><tbody><tr><td>Hall</td><td>A1</td><td>C. botulinum group I</td><td>A1</td><td>ha+/OrfX-</td></tr><tr><td>1430.11</td><td>A5(B')</td><td>C. botulinum group I</td><td>A5, B</td><td>ha+/OrfX-</td></tr><tr><td>2916</td><td>A1(B)</td><td>C. botulinum group I</td><td>A1, bivalent B</td><td>ha+/OrfX+</td></tr><tr><td>Loch Maree</td><td>A3</td><td>C. botulinum group I</td><td>A3</td><td>ha-/OrfX+</td></tr><tr><td>BL6</td><td>B2</td><td>C. botulinum group I</td><td>proteolytic B</td><td>ha+/OrfX-</td></tr><tr><td>168.08</td><td>Bf2</td><td>C. botulinum group I</td><td>bivalent B, F2</td><td>ha+/OrfX+</td></tr><tr><td>NCIB10658</td><td>F6</td><td>C. botulinum group I</td><td>proteolytic F</td><td>ha-/OrfX+</td></tr><tr><td>K16</td><td>E1</td><td>C. botulinum group II</td><td>E1</td><td>ha-/OrfX+</td></tr><tr><td>468</td><td>C1</td><td>C. botulinum group III</td><td>C1</td><td>ha+/OrfX-</td></tr><tr><td>1873</td><td>D</td><td>C. botulinum group III</td><td>D</td><td>ha+/OrfX-</td></tr><tr><td>NCIB10714</td><td>G</td><td>C. botulinum group IV</td><td>G</td><td>ha+/OrfX-</td></tr></tbody></table>

doi:10.1371/journal.pone.0067510.t003
Retained BoNT gene cluster prototype sequences

<figure><img src="image_5.png" /><figcaption>**Figure 3.** *C. botulinum* sequences retrieved by the PathogenID_v2.0 in contaminated food and BLAST analysis results. A schematic representation of the BoNT gene cluster sequences obtained by the RMA for the enrichment culture of the contaminated salad specimen is shown. The orientation and arrangement of these components were reproduced after Franciosa et al [59]. For each component, the highest call-rate observed and the results of the best retrieved BLAST hit(s) (defined as the BLAST hit demonstrating the highest total score) are given i.e. the number of best hits retrieved, as well as the corresponding query sequence coverage, total score, and taxonomy.
doi:10.1371/journal.pone.0067510.g003</figcaption></figure>

characterized microbial strain collections, in order to monitor the molecular evolution of the strains.

The unique design of the PathogenID v2.0 array allows the simultaneous detection and characterization of 48 sequences targeting the BoNT gene cluster components. These sequences could easily be used in the design of specific RMAs with reduced

probe densities by decreasing the number of sequences, thereby diminishing production costs and computing time. Using standard methods to obtain the same level of resolution would imply the use of specific PCRs followed by Sanger sequencing. Although these assays are cheaper as standalone tests, running enough of them to cover the same number of genes with a single sample quickly

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Phylogenetic tree of the *bont/A* sequence retrieved in contaminated food and several toxinotype A strains. A neighbor-joining tree was constructed based on the *bont/A* sequence retrieved by the RMA for the enrichment culture of the *C. botulinum* strain in a naturally contaminated salad specimen, and the corresponding sequences of 6 strains representative of 5 known *C. botulinum* toxinotype A subtypes: strain Hall (A1), NCTC 2916 (A[B]), Kyoto (A2), Loch Maree (A3), 657 (Ba4) and A661222 (A5). Bootstrap values and genetic distance (bar) are shown. The Genbank accession number of each strain is indicated in parentheses. The *bont* sequence obtained from the contaminated salad specimen is marked in bold.
doi:10.1371/journal.pone.0067510.g004</figcaption></figure>
surpasses the RMA in cost, time and logistical complexity of running the analysis. Moreover these methods rely on sequence-specific primer hybridization for target DNA amplification instead of random hexamers for unbiased whole genome amplification, and may therefore fail to amplify and detect emerging variants. Other technologies with universal coverage could be used, such as high-throughput sequencing, however they are still much more laborious and expensive. Finally, the use of an automated approach allows an unbiased, exhaustive and rapid retrieval of taxonomic results.

In summary, the RMA allowed successful detection and fine characterization of neurotoxin-producing clostridia in both pure and polymicrobial cultures. The RMA approach, combined with the automated filtering and retrieval of taxonomic identities, allowed efficient and accurate subtyping of the neurotoxin and detailed characterization of the BoNT gene cluster components. This assay will be used as a higher echelon component in a diagnostics and surveillance pipeline. Importantly, it will be used in further studies to characterise potential outbreak related strains, but also environment isolates, in order to obtain a better picture of the molecular epidemiology of BoNT-producing clostridia.

# Materials and Methods

## Design of the Pathogen ID v2.0 resequencing microarray (RMA)

The second generation of a broad-spectrum resequencing microarray (RMA) was used in this study: the PathogenID v2.0 array, able to detect 124 bacteria, 126 viruses and 673 genes involved in toxin production, pathogenicity or antibiotic resistance [31]. This array included specific housekeeping genes, namely *rrs* and *rpoB*, for detection and characterization of clostridial species. In addition, it included prototype sequences of each component of the neurotoxin gene cluster from all known toxinotypes, as well as several subtypes. Tiled probes targeting these genes were selected using multiple sequence alignments of genome sequences from various toxin subtypes available in GenBank during the time of the RMA design, in order to select conserved regions within toxinotypes for optimal hybridization efficacy. Additional probes were added to improve the detection of specific subtypes, if strains demonstrated sequence variation within the selected probe region. In addition, the *bont/C* and *bont/D* sequences tiled on the array were specifically designed to maximize the probability of detection of mosaic *C. botulinum* toxinotype C/D or D/C strains, by targeting regions highly homologous to either mosaic *bont/C/D* or *bont/D/C* genes. Finally, 19 additional sequences were tiled on the array for detection of a large panel of clostridial toxins other than neurotoxins.

## Bacterial strains and biological samples analysed

Eleven well-characterized *C. botulinum* strains belonging to toxinotypes A through G served as reference: strains Hall (type A1), 1430-11 [type A5(B')], NCTC 2916 [type A1(B)], Loch Maree (type A3), BL6 (type B2), 161.08 (type Bf2), NCIB 10658 (type F6), K16 (type E1), 468 (type C1), 1873 (type D) and NCIB 10714 (type G). The characteristics and origin of most strains have been described elsewhere. The *C. botulinum* type A5(B) strain 1430-11 was isolated from contaminated commercial food (pasta). The *C. botulinum* strain 161.08 was isolated from a foodborne botulism case. Pure cultures of *Clostridium* spp. strains were performed as described previously [55]. Genome sequences from reference strains were retrieved from the genbank database. The *bont* sequences for strain 1430-11 (genbank accession numbers KC683799 and KC683800) and strain 161.08 (genbank accession

numbers KC471328 and KC471329) were obtained by sanger sequencing of specifically PCR-amplified *bont* genes and deposited in genbank.

Biological samples included a naturally contaminated food specimen (salad) and a 48 h enrichment culture of this food sample in Fortified Cooked Meat Medium (FCMM) at 37°C in anaerobic conditions [56]. Presence of *C. botulinum* toxinotype A in the food samples was confirmed by SYBR green real-time PCR with primers P1646 (5'-TCTTACGCGAAATGGTTATGG-3') and P1647 (5'-TGCCTGCACCTAAAAGAGGA-3') for *bont/A* gene, P1648 (5'-CCTGGGCCAGTTTTAAATGA-3') and P1649 (5'-GCGCCTTTGTTTTCTTGAAC-3') for *bont/B* gene, P1650 (5'-GTGCCCGAAGGTGAAAATAA-3') and P1651 (5'-TAATGCTGCTTGCACAGGTT-3') for *bont/E* gene, P2107 (5'-TGACAATGAATTTTCAAAACA-3') and P2108 (5'-TCCAAAAGCATCCATTACTGC-3') for *bont/F* gene. Real time PCR was performed in a total volume of 25 µl containing 12.5 µl of 2× concentration of iQ SYBR green Supermix (Biorad), 5 pmole of each primer, 5 µl of template DNA and 7.3 µl of DNase-RNase free distilled water (Gibco). Amplifications were carried out on a CFX96 Real Time System (Biorad) using 96-well microwell plates and according to the following temperature profile: one cycle of 95°C for 10 minutes, 40 cycles of 95°C for 15 sec and 60°C for 30 sec with fluorescence signal capture at the end of each 60°C step, an extension phase of 1 cycle at 95°C for 60 sec, 60°C for 60 sec and 95°C for 60 sec, followed by a default melt (disassociation) stage.

## Extraction and amplification of bacterial DNA

Total genomic DNA was isolated from *C. botulinum* cultures by lysozyme and proteinase K treatment as described previously [57]. DNA extraction from food samples was performed with Powerfood Microbial DNA isolation kit (MoBio Laboratories Inc) according to the manufacturer's recommendations. After extraction, DNA was amplified using the whole genome amplification (WGA) protocol (RepliG Midi kit, Qiagen) as described previously [54].

## Resequencing microarray assay

The amplification products obtained from genomic DNA by WGA were quantified by Quantit BR (Invitrogen) according to the manufacturer's instructions. A recommended amount of target DNA was fragmented and labeled according to the GeneChip Resequencing Assay manual (Affymetrix). The obtained products were coded with unrelated numbers by a non-observer to ensure that the study was performed blindly. The microarray hybridization process was carried out according to the protocol recommended by the manufacturer (Affymetrix). The details and parameter settings for the data analysis (essentially conversion of raw image files obtained from scanning of the microarrays into FASTA files containing the sequences of called bases for each tiled region of the microarray) have been described previously [54]. The base call rate refers to the percentage of base calls generated from the full-length tiled sequence.

## Automated data analysis

A Perl-based program, designated "TaxFinder", was used for the automated analysis of re-sequencing data provided by PathogenID v2.0. The program read the FASTA file generated for each sample, which contains all the sequences read by the GSEQ software (Affymetrix) from the hybridization results. A filtering process, based on the one described by Malanoski et al. [58] was applied to these sequences. For each sequence, TaxFinder considers the Nb first bases of the sequence, Nb
defined by the user in the beginning of the experiment and corresponding to the minimal length accepted for the sequence to keep. If the percentage of No-call in this fragment is inferior to an elongation threshold defined by the user, a base is added to the fragment and the percentage of No-call is recalculated. As soon as it is higher than the elongation threshold, the No-call at the ends of the fragment are removed. If the length of the nucleotide fragment is still higher than the minimal threshold established, then it is conserved. If the percentage is higher than the threshold, the fragment is skipped, and the program considers the fragment of size Nb and starting position one base upstream than the precedent fragment. This process is reiterated until the end of the sequence is reached. In this study, sequences that did not contain subsequences with a minimum nucleotide length of 20 nt and a maximum undetermined nucleotide content of 10% were discarded. Filtered sequences individually underwent a blastn analysis to search for sequence homologues in the NCBI nucleotide collection (nr/nt database) with adjusted settings to restrict the search to bacteria entries. Blast + application distributed by NCBI is used for the automated blast research (http://www.ncbi.nlm.nih.gov/books/NBK1763/). The following algorithm parameters were modified: the expected value cut-off was fixed to 100, the minimum word size was set to 7 and the upper limit of displayed descriptions of database sequences per query was decreased to 50. The best BLAST hits were classified according to their total score (the sum of the high-scoring segment

pairs) and their corresponding taxonomies were retrieved from the NCBI Taxonomy database. When several hits obtained the highest total score, the script automatically retrieved the taxonomies of the 10 first BLAST hits.

## Phylogenetic analysis

The resulting *bont* sequence obtained from the enrichment culture of the *C. botulinum* strain in a contaminated food specimen was compared with the corresponding *bont* sequences of reference strains from the 5 known *C. botulinum* toxinotype A subtypes (GenBank accessed 21/02/2013). Multiple sequence alignment was performed using the CLC Bio software and checked for accuracy by eye. A neighbor-joining tree of these sequences was constructed using the Jukes-Cantor method with the SeaView v4.2.1 software. The level of support for each node is provided by 100 bootstrap replications.

## Author Contributions

Conceived and designed the experiments: JV NB PD JCM MRP. Performed the experiments: JV NB CM. Analyzed the data: JV NB CM GG. Contributed reagents/materials/analysis tools: GG VC. Wrote the paper: JV JCM MRP. Designed the content of Resequencing microarray: NB CM RS TV GCK STC MRP. Designed the software used in analysis: GG VC.

## References

1. Sobel J (2005) Botulism. Clin Infect Dis 41: 1167-1173. doi:10.1086/444507.

2. Arnon SS, Schechter R, Inglesby TV, Henderson DA, Bartlett JG, et al. (2001) Botulinum Toxin as a Biological Weapon. JAMA 285: 1059-1070. doi:10.1001/jama.285.8.1059.

3. Villar RG, Elliott SP, Davenport KM (2006) Botulism: The Many Faces of Botulinum Toxin and its Potential for Bioterrorism. Infect Dis Clin North Am 20: 313-327. doi:16/j.idc.2006.02.003.

4. Peck MW (2009) Biology and Genomic Analysis of Clostridium botulinum. Academic Press, Vol. 55. 183-265, 320. Available: http://www.sciencedirect.com/science/article/pii/S0065291109055039. Accessed 2011 Aug 3.

5. Brüggemann H, Woller A, Mazuet C, Popoff MR (2011) Clostridium botulinum. Genomes of Foodborne and Waterborne Pathogens. Washington DC: Fratamico P., Liu Y., Kathariou S. 185-212.

6. Smith TJ, Hill KK, Foley BT, Detter JC, Munk AC, et al. (2007) Analysis of the Neurotoxin Complex Genes in Clostridium botulinum A1-A4 and B1 Strains: BoNT/A3/Ba4 and/B1 Clusters Are Located within Plasmids. PLoS One. 2.doi:10.1371/journal.pone.0001271.

7. Hill KK, Smith TJ, Helma CH, Ticknor LO, Foley BT, et al. (2007) Genetic Diversity among Botulinum Neurotoxin-Producing Clostridial Strains. J Bacteriol 189: 818-832. doi:10.1128/JB.01180-06.

8. Hutson RA, Zhou Y, Collins MD, Johnson EA, Hatheway CL, et al. (1996) Genetic Characterization of Clostridium botulinum Type A Containing Silent Type B Neurotoxin Gene Sequences. J Biol Chem 271: 10786-10792. doi:10.1074/jbc.271.18.10786.

9. Santos-Buelga JA, Collins MD, East AK (1998) Characterization of the Genes Encoding the Botulinum Neurotoxin Complex in a Strain of Clostridium botulinum Producing Type B and F Neurotoxins. Curr Microbiol 37: 312-318. doi:10.1007/s002849900384.

10. Kirma N, Ferreira JL, Baumstark BR (2004) Characterization of six type A strains of Clostridium botulinum that contain type B toxin gene sequences. FEMS Microbiol Lett 231: 159-164. doi:10.1016/S0378-1097(03)00911-X.

11. Lúquez C, Raphael BH, Maslanka SE (2009) Neurotoxin Gene Clusters in Clostridium botulinum Type Ab Strains. Appl Environ Microbiol 75: 6094-6101. doi:10.1128/AEM.01009-09.

12. Moriishi K, Koura M, Abe N, Fujii N, Fujinaga Y, et al. (1996) Mosaic structures of neurotoxins produced from Clostridium botulinum types C and D organisms. Biochim Biophys Acta 1307: 123-126. doi:10.1016/0167-4781(96)00006-1.

13. Brüggemann H (2005) Genomics of clostridial pathogens: implication of extrachromosomal elements in pathogenicity. Curr Opin Microbiol 8: 601-605. doi:16/j.mib.2005.08.006.

14. Smith TJ, Lou J, Geren IN, Forsyth CM, Tsai R, et al. (2005) Sequence Variation within Botulinum Neurotoxin Serotypes Impacts Antibody Binding and Neutralization. Infect Immun 73: 5450-5457. doi:10.1128/IAI.73.9.5450-5457.2005.

15. Hill KK, Xie G, Foley BT, Smith TJ, Munk AC, et al. (2009) Recombination and insertion events involving the botulinum neurotoxin complex genes in Clostridium botulinum types A, B, E and F and Clostridium butyricum type E strains. BMC Biol 7: 66-66. doi:10.1186/1741-7007-7-66.

16. Fach P, Fenicia L, Knutsson R, Wielinga PR, Anniballi F, et al. (2011) An innovative molecular detection tool for tracking and tracing Clostridium botulinum types A, B, E, F and other botulinum neurotoxin producing Clostridia based on the GeneDisc cycler. Int J Food Microbiol 145: S145-S151. doi:16/j.ijfoodmicro.2010.04.006.

17. Raphael BH, Joseph LA, McCroskey LM, Lúquez C, Maslanka SE (2010) Detection and differentiation of Clostridium botulinum type A strains using a focused DNA microarray. Mol Cell Probes 24: 146-153. doi:16/j.mcp.2009.12.003.

18. Carter AT, Paul CJ, Mason DR, Twine SM, Alston MJ, et al. (2009) Independent evolution of neurotoxin and flagellar genetic loci in proteolytic Clostridium botulinum. BMC Genomics 10: 115-115. doi:10.1186/1471-2164-10-115.

19. Sebaihia M, Peck MW, Minton NP, Thomson NR, Holden MTG, et al. (2007) Genome sequence of a proteolytic (Group I) Clostridium botulinum strain Hall A and comparative analysis of the clostridial genomes. Genome Res 17: 1082-1092. doi:10.1101/gr.6282807.

20. Umeda K, Seto Y, Kohda T, Mukamoto M, Kozaki S (2010) A novel multiplex PCR method for Clostridium botulinum neurotoxin type A gene cluster typing. Microbiol Immunol 54: 308-312. doi:10.1111/j.1348-0421.2010.00213.x.

21. Umeda K, Seto Y, Kohda T, Mukamoto M, Kozaki S (2009) Genetic Characterization of Clostridium botulinum Associated with Type B Infant Botulism in Japan. J Clin Microbiol 47: 2720-2728. doi:10.1128/JCM.00077-09.

22. Leclair D, Pagotto F, Farber JM, Cadieux B, Austin JW (2006) Comparison of DNA Fingerprinting Methods for Use in Investigation of Type E Botulism Outbreaks in the Canadian Arctic. J Clin Microbiol 44: 1635-1644. doi:10.1128/JCM.44.5.1635-1644.2006.

23. Fillo S, Giordani F, Anniballi F, Gorgé O, Ramisse V, et al. (2011) Clostridium botulinum Group I Strain Genotyping by 15-Locus Multilocus Variable-Number Tandem-Repeat Analysis. J Clin Microbiol 49: 4252-4263. doi:10.1128/JCM.05396-11.

24. Lindström M, Hinderink K, Somervuo P, Kiviniemi K, Nevas M, et al. (2009) Comparative Genomic Hybridization Analysis of Two Predominant Nordic Group I (Proteolytic) Clostridium botulinum Type B Clusters. Appl Environ Microbiol 75: 2643-2651. doi:10.1128/AEM.02557-08.

25. Han T, Chang CW, Kwekel JC, Chen Y, Ge Y, et al. (2012) Characterization of whole genome amplified (WGA) DNA for use in genotyping assay development. BMC Genomics 13: 217. doi:10.1186/1471-2164-13-217.

26. Erlandsson L, Rosenstierne MW, McLoughlin K, Jaing C, Fomsgaard A (2011) The Microbial Detection Array Combined with Random Phi29-Amplification
Used as a Diagnostic Tool for Virus Detection in Clinical Samples. PLoS One 6: e22631. doi:10.1371/journal.pone.0022631.

27. Bouzid M, Heavens D, Elwin K, Chalmers RM, Hadfield SJ, et al. (2010) Whole genome amplification (WGA) for archiving and genotyping of clinical isolates of Cryptosporidium species. Parasitology 137: 27-36. doi:10.1017/S0031182009991132.

28. Gadkar V, Rillig MC (2005) Suitability of genomic DNA synthesized by strand displacement amplification (SDA) for AFLP analysis: genotyping single spores of arbuscular mycorrhizal (AM) fungi. J Microbiol Methods 63: 157-164. doi:10.1016/j.mimet.2005.03.009.

29. Berthet N, Reinhardt AK, Leclercq I, van Ooyen S, Batejat C, et al. (2008) Phi29 polymerase based random amplification of viral RNA as an alternative to random RT-PCR. BMC Mol Biol 9: 77. doi:10.1186/1471-2199-9-77.

30. Zwick ME, Mcafee F, Cutler DJ, Read TD, Ravel J, et al. (2005) Microarray-based resequencing of multiple Bacillus anthracis isolates. Genome Biol 6: R10-R10. doi:10.1186/gb-2004-6-1-r10.

31. Dacheux L, Berthet N, Dissard G, Holmes EC, Delmas O, et al. (2010) Application of Broad-Spectrum Resequencing Microarray for Genotyping Rhabdoviruses. J Virol 84: 9557-9574. doi:10.1128/JVI.00771-10.

32. Wang Z, Malanoski AP, Lin B, Kidd C, Long NC, et al. (2008) Resequencing microarray probe design for typing genetically diverse viruses: human rhinoviruses and enteroviruses. BMC Genomics 9: 577-577. doi:10.1186/1471-2164-9-577.

33. Lin B, Malanoski AP (2009) Resequencing Arrays for Diagnostics of Respiratory Pathogens. In: Dufva M, DNA Microarrays for Biomedical Research. Totowa, NJ: Humana Press, Vol. 529. 231-257. Available: http://www.springerlink.com/content/q1102392j0782q11/#section = 55300&page = 1. Accessed 2011 Aug 9.

34. Oriot C, D'Aranda E, Castanier M, Glazal M, Galy C, et al. (2011) One collective case of type A foodborne botulism in Corsica. Clin Toxicol 49: 752-754. doi:10.3109/15563650.2011.606222.

35. Lai LS, Wang YM, Lin CH (2011) Foodborne botulinum type E intoxication associated with dried bean curd: first case report in Taiwan. Acta Neurol Taiwan 20: 138-141.

36. Jalava K, Selby K, Pihlajasaari A, Kolho E, Dahlsten E, et al. (2011) Two cases of food-borne botulism in Finland caused by conserved olives, October 2011. Euro Surveill. 16: 20034. Available: http://www.eurosurveillance.org/ViewArticle.aspx?ArticleId=20034. Accessed 2011 Dec. 27.

37. Browning LM, Prempeh H, Little C, Houston C, Grant K, et al. (2011) An outbreak of food-borne botulism in Scotland, United Kingdom, November 2011. Euro Surveill. 16: 20036. Available: http://www.eurosurveillance.org/ViewArticle.aspx?ArticleId=20036. Accessed 2011 Dec 27.

38. Pingeon JM, Vanbockstael C, Popoff MR, King LA, Deschamps B, et al. (2011) Two outbreaks of botulism associated with consumption of green olive paste, France, September 2011. Euro Surveill. 16: 20035. Available: http://www.eurosurveillance.org/ViewArticle.aspx?ArticleId=20035. Accessed 2011 Dec 27.

39. Daminelli P, De Nadai V, Bozzo G, Finazzi G, Oliverio E, et al. (2011) Two unlinked cases of foodborne botulism in Italy at the beginning of 2010. New Microbiol 34: 287-290.

40. Date K, Fagan R, Crossland S, MacEachern D, Pyer B, et al. (2011) Three Outbreaks of Foodborne Botulism Caused by Unsafe Home Canning of Vegetables Ohio and Washington, 2008 and 2009. J Food Prot 74: 2090-2096. doi:10.4315/0362-028X.JFP-11-128.

41. Case definitions for infectious conditions under public health surveillance. Centers for Disease Control and Prevention (1997). MMWR Recomm Rep 46: 1-55.

42. Therre H. (1999) Botulism in the European Union. Euro Surveill. 4: 2-7. Available: http://www.eurosurveillance.org/ViewArticle.aspx?ArticleId=48. Accessed 28 December 2011.

43. Lindström M, Korkeala H (2006) Laboratory Diagnostics of Botulism. Clin Microbiol Rev 19: 298-314. doi:10.1128/CMR.19.2.298-314.2006.

44. Pandya GA, McEllistrem MC, Venepally P, Holmes MH, Jarrahi B, et al. (2011) Monitoring the Long-Term Molecular Epidemiology of the Pneumococcus and Detection of Potential "Vaccine Escape" Strains. PLoS One 6. doi:10.1371/journal.pone.0015950.

45. Corless CE, Kaczmarski E, Borrow R, Guiver M (2008) Molecular Characterization of Neisseria meningitidis Isolates Using a Resequencing DNA Microarray. J Mol Diagn 10: 265-271. doi:10.2353/jmoldx.2008.070152.

46. Davignon L, Walter EA, Mueller KM, Barrozo CP, Stenger DA, et al. (2005) Use of Resequencing Oligonucleotide Microarrays for Identification of Streptococcus pyogenes and Associated Antibiotic Resistance Determinants. J Clin Microbiol 43: 5690-5695. doi:10.1128/JCM.43.11.5690-5695.2005.

47. Hutson RA, Thompson DE, Lawson PA, Schocken-Itturino RP, Böttger EC, et al. (1993) Genetic interrelationships of proteolytic Clostridium botulinum types A, B, and F and other members of the Clostridium botulinum complex as revealed by small-subunit rRNA gene sequences. Antonie Van Leeuwenhoek 64: 273-283.

48. Pourshaban M, Franciosa G, Fenicia L, Aureli P (2002) Taxonomic identity of type E botulinum toxin-producing Clostridium butyricum strains by sequencing of a short 16S rDNA region. FEMS Microbiol Lett 214: 119-125.

49. Sasaki Y, Takikawa N, Kojima A, Norimatsu M, Suzuki S, et al. (2001) Phylogenetic positions of Clostridium novyi and Clostridium haemolyticum based on 16S rDNA sequences. Int J Syst Evol Microbiol 51: 901-904.

50. Franciosa G, Floridi F, Maugliani A, Aureli P (2004) Differentiation of the gene clusters encoding botulinum neurotoxin type A complexes in Clostridium botulinum type A, Ab, and A(B) strains. Appl Environ Microbiol 70: 7192-7199. doi:10.1128/AEM.70.12.7192-7199.2004.

51. Raphael BH, Luquez C, McCroskey LM, Joseph LA, Jacobson MJ, et al. (2008) Genetic Homogeneity of Clostridium botulinum Type A1 Strains with Unique Toxin Gene Clusters. Appl Environ Microbiol 74: 4390-4397. doi:10.1128/AEM.00260-08.

52. Leski TA, Malanoski AP, Gregory MJ, Lin B, Stenger DA (2011) Application of a Broad-Range Resequencing Array for Detection of Pathogens in Desert Dust Samples from Kuwait and Iraq. Appl Environ Microbiol 77: 4285-4292. doi:<p>≥10.1128/AEM.00021-11</p>.

53. Leski TA, Lin B, Malanoski AP, Wang Z, Long NG, et al. (2009) Testing and Validation of High Density Resequencing Microarray for Broad Range Biothreat Agents Detection. PLoS One 4. doi:10.1371/journal.pone.0006569.

54. Berthet N, Dickinson P, Filliol I, Reinhardt AK, Batejat C, et al. (2008) Massively parallel pathogen identification using high-density microarrays. Microb Biotechnol 1: 79-86.

55. Fach P, Micheau P, Mazuet C, Perelle S, Popoff M (2009) Development of real-time PCR tests for detecting botulinum neurotoxins A, B, E, F producing Clostridium botulinum, Clostridium baratii and Clostridium butyricum. J Appl Microbiol 107: 465-473. doi:10.1111/j.1365-2672.2009.04215.x.

56. Takeda M, Tsukamoto K, Kohda T, Matsui M, Mukamoto M, et al. (2005) Characterization of the neurotoxin produced by isolates associated with avian botulism. Avian Dis 49: 376-381.

57. Dineen SS, Bradshaw M, Johnson EA (2003) Neurotoxin Gene Clusters in Clostridium botulinum Type A Strains: Sequence Comparison and Evolutionary Implications. Curr Microbiol 46: 345-352. doi:10.1007/s00284-002-3851-1.

58. Malanoski AP, Lin B, Wang Z, Schnur JM, Stenger DA (2006) Automated identification of multiple micro-organisms from resequencing DNA microarrays. Nucleic Acids Res 34: 5300-5311. doi:10.1093/nar/gkl565.

59. Franciosa G, Maugliani A, Floridi F, Aureli P (2006) A novel type A2 neurotoxin gene cluster in Clostridium botulinum strain Mascarpone. FEMS Microbiol Lett 261: 88-94. doi:10.1111/j.1574-6968.2006.00331.x.