Figures
Abstract
Paenibacillus larvae is a Gram positive bacterium and the causative agent of the most widespread fatal brood disease of honey bees, American foulbrood (AFB). A total of thirty-three independent Paenibacillus larvae isolates from various geographical origins in North America and five reference strains were investigated for genetic diversity using multilocus sequence typing (MLST). This technique is regarded to be a powerful tool for epidemiological studies of pathogenic bacteria and is widely used in genotyping assays. For MLST, seven housekeeping gene loci, ilvD (dihydroxy-acid dyhydrogenase), tri (triosephosphate isomerase), purH (phospharibosyl-aminoimidazolecarboxamide), recF (DNA replication and repair protein), pyrE (orotate phosphoribosyltransferase), sucC (succinyl coenzyme A synthetase β subunit) and glpF (glycerol uptake facilitator protein) were studied and applied for primer designs. Previously, ERIC type DNA fingerprinting was applied to these same isolates and the data showed that almost all represented the ERIC I type, whereas using BOX-PCR gave an indication of more diversity. All isolates were screened for resistance to four antibiotics used by U.S. beekeepers, showing extensive resistance to tetracycline and the first records of resistance to tylosin and lincomycin. Our data highlight the intraspecies relationships of P. larvae and the potential application of MLST methods in enhancing our understanding of epidemiological relationships among bacterial isolates of different origins.
Citation: Krongdang S, Evans JD, Pettis JS, Chantawannakul P (2017) Multilocus sequence typing, biochemical and antibiotic resistance characterizations reveal diversity of North American strains of the honey bee pathogen Paenibacillus larvae. PLoS ONE 12(5): e0176831. https://doi.org/10.1371/journal.pone.0176831
Editor: Daniele Daffonchio, University of Milan, ITALY
Received: November 16, 2016; Accepted: April 18, 2017; Published: May 3, 2017
This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: This research was financially supported by a grant under The Program of Strategic Scholarships for Frontier Research Network for the Joint Ph.D. Program Thai Doctoral degree from the Office of the Higher Education Commission, Thailand (Grant No. CHE-SW215/2551), Thailand Research Fund and Chiang Mai University fund. We wish to thank the USDA-ARS Bee Research Laboratory (Beltsville, MD, USA), and the Graduate School of Chiang Mai University. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
American foulbrood (AFB), caused by the Gram-positive bacterium Paenibacillus larvae (P. larvae) [1,2], is regarded as the most serious bacterial disease affecting the brood of honey bees (mainly Apis mellifera) worldwide [3]. P. larvae can form spores that remain viable for years and survive under extreme conditions. The bacteria can be readily transmitted by swarming, drifting, or foraging bees or by beekeepers when they exchange hive equipment and hive materials between colonies or apiaries [2,4–6]. A diseased colony can be visualized and detected only when larvae are observed with clinical symptoms of AFB, however colonies may harbor many P. larvae spores without showing symptoms [7]. In many countries, including the United States, infected honey bee colonies are treated with antibiotics when clinical AFB symptoms are present. Oxytetracycline (OTC = Terramycin®) is the antibiotic most often applied to treat this disease [8]. More recently, tylosin [9,10] and lincomycin [11] have been used to cure AFB [11]. In recent years, not only has AFB become more prevalent but resistance to antibiotics by P. larvae has been observed particularly oxytetracycline [12,13].
Apart from field inspection of AFB infected colonies; laboratory tests are used to confirm the presence of P. larvae. Microscopic identification of P. larvae is commonly complemented by using biochemical diagnosis and molecular techniques. Molecular typing methods are also essential in order to better understand the epidemiology of such opportunistic agents with the final goal of preventing contamination[14]. Currently, the genetic diversity of P. larvae has been studied by using various methods such as a repetitive-element PCR fingerprinting with primers ERIC [15–17], BOX-[17,18], MBO REP1, and pulsed-field gel electrophoresis (PFGE) [15,16,19,20]. Nevertheless, these methods are difficult to standardize between laboratories especially with large collection of bacterial strains. Therefore, standardized approach would greatly help us to study population structure and evolution of the P. larvae. Multilocus sequence typing (MLST) has now become standard for assessing genetic variation in worldwide collections of bacterial isolates. MLST schemes discriminate several bacterial and yeast species at the strain level by using information from across multiple independent markers [21]. In addition, the level of genomic plasticity of P. larvae has not been clarified, this highlights the need for a new robust tools for strain differentiation using phylogenetic analysis [22]. In this study, we designed an MLST scheme for the P. larvae from infected colonies from different geographical organism in United States. We used non-overlapping loci and our research was done in parallel to and similar to the approach recently invoked for European strains of this species by Morrissey et al., 2015 [23]. We compared sequence data against traits important for the evolution, phylogeny, and epidemiology of these bacteria, including biochemistry profiles and antibiotic susceptibility.
Materials and methods
Bacterial isolates
A total of 33 field isolates of Paenibacillus larvae were sampled from those received by the Bee Disease Diagnostic Service (USDA-ARS, MD, US) from different regions in the United States in the years 1999–2013 (Fig 1). These isolates were contrasted with P. larvae reference strains ATCC9545, ATCC49843, 233/00, LMG16247 and LMG16252 (described in S1 Table). Honey bee larvae with AFB symptoms were dissolved in 200 μl sterile distilled water and incubated at 80°C for 15 min in order to kill vegetative cells of bacteria, yeasts and non-spore forming bacteria before culturing [18]. The isolates were cultivated using THClYGP agar plates [15] using the following recipe (per liter): 15 g yeast extract (Fisher Science), 1 g pyruvic acid (EMBSCO), 200 ml 0.1 M Tris-HCl, pH 7.0, 20 g agar, 40 ml 10% glucose and 0.6 ml nalidixic acid. Plates were incubated at 37°C for 3 days.
For more information on sample location and or identification please see S1 Table.
Primer design for MLST
Thirteen housekeeping loci were chosen for the MLST primer design using the following criteria: presence as a single copy in all strains, rate of sequence conservation and wide distribution across the chromosome. Bacillus cereus and Bacillus licheniformis MLST databases were used as templates from PubMLST (http://pubmlst.org). Full-length sequences of seven genes were identified from the Paenibacillus larvae DSM 25430 genome assembly (NC_023134.1). Nucleotide sequences were downloaded from GenBank and aligned with template allele profiles and analyzed by using the BLAST and FASTA software available at the National Center for Biotechnology website (http://www.ncbi.nlm.nih.gov). Specific primers of each gene were designed from the most conserved regions using Primer3 and BLAST software (http://www.ncbi.nlm.nih.gov/tools/primer-blast/) within suitable length between 400–700 bp and similar % GC and annealing temperatures.
DNA extraction and DNA fingerprinting by using ERIC- and BOX PCR
Preparations of bacterial DNA for DNA fingerprinting and MLST analysis were performed as described by Genersch and Otten [19]. A total genomic DNA was isolated from P. larvae colonies grown on culture plates using a 6% InstaGene matrix (Bio-Rad, Hercules, CA) following the manufacturer’s instructions. PCR amplification was performed using primer ERIC1R (5’-ATGTAAGCTCCTGGGGATTCAC- 3’) and ERIC2 (5’-AAGTAAGTGACTGGGGTGAGCG-3’)[24] for ERIC gene and BOX A1R (5’-CTACGGCAAGGCGACGCTGACG-3’) for BOX gene [24, 25]. PCR were carried out in a final volume 50 μl which consisted of 1X PCR buffer (containing 1.5 mM MgCl2); 1x Q-Solution (for ERIC amplification); deoxynucleoside triphosphate mix (dGTP, dTTP, dATP and dCTP, each at 2 mM); 0.2 μM of each primer; 2.5 U HotStarTaq DNA polymerase (QIAGEN, Germantown, MD) and 10 ng of DNA. The thermal cycle program consisted of initial activation step (15 min, 95°C), followed by 35 cycles of denaturation at 94°C for 1 min, annealing for 1 min at 53°C and extension for 2.5 min at 72°C, followed by a final elongation step at 72°C for 10 min. 5 μl of the ERIC and the BOX-PCR products were analyzed in 1.5% agarose gel, stained with ethidium bromide, visualized and photographed under UV light.
MLST PCR amplification
For MLST, PCR reactions were carried out in a final volume of 30 μl containing 10 ng of DNA and 0.2 μM of each primer along with the PCR reactions described above. Cycling conditions included an initial denaturing step of 95°C for 15 min, 35 cycles at 94°C for 30 sec, 55°C for 45 sec, and 72°C for 1 min, and a final extension step at 72°C for 10 min. Amplification products were resolved electrophoretically by 1.5% agarose gel, stained with ethidium bromide and photographed under UV transillumination to confirm an expected size and single clear band per PCR. Then, PCR products were purified directly using QIAquick PCR purification kit (QIAGEN, Germantown, MD) following the manufacturer’s instructions. The sequencing was performed according to Big Dye 2.0 (Applied Biosystems, Carlsbad, CA) according to the manufacturer’s instructions. DNA sequences obtained in this study were submitted to GenBank.
MLST and Phylogenetic analysis
Sequences from all loci were concatenated prior to analysis. An alignment of the sequences obtained together with reference sequences was constructed for each of genetic markers using the program MEGA6 [26]. Neighbor joining trees were constructed using the program MEGA6, based on the evolutionary distances calculated by the Kimura-2-parameter model. The reliability of these trees was assessed using bootstrap analyses across 1000 replicates. DNA polymorphism was calculated in DnaSP 5.10.01 [27]. Number of haplotypes (Hap) and haplotype diversity (Hd) were calculated. The nucleotide diversity per site (θ) and the average number of nucleotide difference between populations (π) were determined as well. The neutrality of mutation was analyzed by Tajima’s D [28], Fu and Li’s D*, Fu and Li F* [29] and Fu’s Fs [30] selective neutrality tests. A unique number was assigned to each distinct allele sequence for each locus. Each unique combination of seven loci numbers was determined the sequence type (ST) of each strain based on the analysis of nucleotide polymorphisms. The same ST was used for isolates sharing the same allelic profiles and calculated for relatedness using the eBURST v3 program [31].
Phenotypic characterizations of P. larvae
The field isolates were characterized biochemically using the commercially available API 50 CHB, API ZYM and API 20E systems (bioMérieux, Boston, MA). The tests were carried out according to the manufacturer’s instructions with minor modifications as follows: P. larvae colonies from THClYGP agar were suspended in 3 ml sterile saline solution in order to obtain a turbidity of 5–6 on the McFarland scale. Each well was filled with three to four drops of the culture suspension. The dendrogram obtained from biochemistry profile was constructed by the UPGMA method applying Pearson’s correlation.
Antibiotic susceptibility testing
Antimicrobial susceptibility tests were assessed by the disk diffusion method. Four antibiotic drugs; OTC (30 μg/disk), tetracycline (5 μg/disk), tylosin (30 μg/disk) and lincomycin (2 μg/disk) were examined. P. larvae strains were cultured on the agar MYPGP for 48 h of incubation at 37°C were suspended in sterile distilled water and adjusted to approximately 2.87 × 108 cells/ml (A620 nm = 0.388, equivalent to a McFarland standard no. 1[10]. Diameters of the inhibition zones were measured. Three replicates were tested for each isolate strain.
Results
Genotyping by ERIC- and BOX PCR analysis
ERIC PCR is generally used to characterize P. larvae. In comparison to reference strains which described ERIC types as ATCC 9545I (ERIC I), 233/00 (ERIC II), LMG16252 (ERIC III), ERIC ATCC 49843 and LMG 16247 (ERIC IV), all P. larvae isolates were clustered together and showing the only one of distinguishable genotypes, in ERIC I genotype with fragments ranging from 200 to 3000 bp (S1 Fig).
Genomic fingerprints of P. larvae isolates also were found by using BOX and the resulting dendrogram revealed UPGMA cluster with six different clusters (Fig 2). Two of the total which are named genotype BOX A and BOX B corresponded to genotypes described by Alippi and Marioaguilar [20]. Additionally, Genersch and Otten [19] and Loncaric et al. [17] who labeled these as genotype a and A. Interestingly, new genotypes were observed and given names as BOX D and BOX D1. The new genotypes present two band sizes around 200 to 300 bp, which were absent in the others. BOX D and BOX D1 showed highly similar patterns and bands around 200–300 bp which were missing in others. The band size 200 bp can be observed in BOX D1 and the reference strain 233/00. Hence, as expected BOX-PCR provided a better discrimination than ERIC-PCR, because all isolates had been included in ERIC I genotype (S2 Fig).
Sequence polymorphism and multilocus genotypes
Thirteen housekeeping gene were successfully amplified and tested for nucleotide diversity in 33 isolates and 5 reference strains (S3 Fig). Seven housekeeping gene loci, ilvD (dihydroxy-acid dyhydrogenase), tri (triosephosphate isomerase), purH (phospharibosyl-aminoimidazolecarboxamide), recF (DNA replication and repair protein), pyrE (orotate phosphoribosyltransferase), SucC (succinyl coenzyme A synthetase β subunit) and glpF (glycerol uptake facilitator protein) were chosen for the MLST analysis (Table 1) and others showed no diversity within the sequence (S2 Table). All sequences have been deposited in GenBank (https://www.ncbi.nlm.nih.gov/genbank/) under accession numbers from: KY673263-KY673300 for ilvD, KY673301-KY673338 for purH, KY673339-KY673376 for tri, KY673377-KY673414 for recF, KY673415-KY673452 for pyrE, KY673453- KY673490 for sucC and KY673491-KY673528 for glpK.
On average, housekeeping genes had a low level of nucleotide (π = 0.00086±0.00017) and haplotype diversity (Hd = 0.19885±0.07585) in P. larvae isolates. glpK had the highest sequence polymorphism and revealed the presence of 6 genotypes with an Hd of 0.328. Whereas tri locus showed the lowest sequence polymorphism, with two distinct genotypes, and an Hd of 0.149 (Table 2). The number of alleles and polymorphic sites identified per locus was overall low. This was reflected in the low θ and π values calculated for each gene. The G+C mol content varied between 0.4407 for the partial recF sequence to 0.5077 for the ilvD sequence. By Tajima’s D, Fu and Li’s D, Fu and Li’s F test, two of the housekeeping genes revealed signs of significant deviation from neutrality (Table 2); however, Fu and Fs’ D test did not show significant. Therefore, these tests suggest that these housekeeping genes are primarily evolving in a neutral fashion and can be used as genetic markers for inferring population and evolutionary history of P. larvae.
MLST assays showed a high degree of evolutionary divergence between North American P. larvae isolates. The neighbor-joining tree constructed using the concatenated sequences of the seven housekeeping loci identified major branches with acceptable (50%) bootstrap values. Two main genetically isolated groups were observed (Fig 3A). Group 1 consisted of 35 isolates and group 2 contained 3 isolates. The data in the eBURST analysis showed 15 different sequence types (Table 3). All STs could be grouped into 2 groups. Ten STs which belonged in group 1 are ST1, ST2, ST3, ST7, ST8, ST9, ST10, ST12, ST13, ST14 and ST15 which ST1 were identified as founder. Three STs are belonged to group II (ST3, ST5 and ST6). All most P. larvae from United States were classified into a distinct clade (Group 1) while reference isolates that identified in ERIC type III and IV were grouped into another clade (Group 2). The other 2 STs that did not belong to any groups were identified as singletons (ST4 and ST11) (Fig 3B).
(A) MLST Phylogenetic tree (neighbor-joining) of 38 Paenibacillus larvae isolates after concatenating the sequences of seven loci. The phylogenetic tree was constructed using the neighbor-joining method with 1,000 bootstraps (bootstrap values are shown above the nodes). Phylogenetic analyzes were constructed in MEGA6. (B) eBURST analysis of MLST profiles of 15 STs of Paenibacillus larvae isolates.
Antibiotic susceptibility
Susceptibility of antibiotics was examined by using antibiotic disks containing 30 μg oxytetracycline (OTC), 5 μg tetracycline, 30 μg tylosin and 2 μg lincomycin. A histogram plot shows contributions of clear zone diameters in each antibiotics group (Fig 4).
The histogram showed that susceptibility toward antibiotics of P. larvae isolates has similar trends. The plot showed a clear split between resistant and susceptible isolates with a defined zone of inhibition difference of 2.0 cm between resistant and susceptible isolates. In this work, OTC, tetracycline, tylosin and lincomycin-resistant isolates were found as 15.78% (PL-3, PL-4, PL-5, PL-15, PL-29, and PL-30), 21.05% (PL-3, PL-4, PL-5, PL-9, PL-15, PL-28, PL-29, and PL-30), 7.8% (PL-16, PL-18 and PL-26) and 7.8% (PL-18, PL-26 and ATCC 9545) of our screened isolates, respectively. Susceptibility of two antibiotics, OTC and tetracycline were tightly correlated (R2 = 0.8706) as expected due to their similar mode of action and molecular structure. In addition, resistance to tylosin and lincomycin (R2 = 0.5757) tended to occur in identical isolates, despite different modes of action for these compounds. No further correlations across isolates for antibiotic susceptibility between treatments were found (Fig 5).
Diameter of inhibition zone of four antibiotics were scattered and correlation (r2) in each two antibiotics were shown.
Biochemical characteristics of P. larvae
Carbohydrate utilization was characterized by using API® 50CH system with 49 carbohydrates (S3 Table). In terms of metabolic differences of P. larvae isolates, the dendrogram was constructed by using biochemistry profiles. Five different metabolic profiles were observed by using UPGMA method. The isolates defined as ERIC I mostly belong to group 1 and 4 whereas, the reference strains, ATCC 49843 and LMG 16247 (ERIC IV), 233/00 (ERIC II), and LMG 16252 (ERIC III) were clearly related with ERIC DNA fingerprinting (Fig 6).
All P. larvae hydrolyzed carbon sources such as esculin ferric citrate and potassium-2-ketogluconate and phosphorus sources effectively. 97.39% utilized N-acetylglucosamine and D-trehalose (excluding PL-13) and 94.74% (36/38) used glycerol (excluding PL-13, LMG 16252), D-ribose (excluding PL-13, PL-15). Abundant hydrolytic enzyme productions were detected during inoculates time with API®ZYM. The results showed all P. larvae isolates produce esterase (C4), esterase lipase (C8) acid phosphatase, naphthol-AS-BI-phosphohydrolase and N-acetyl-β-glucosaminidase. 97.39% produced alkaline phosphatase (excluding isolate 233/00). Using the API®20E system, all isolates showed the positive reaction for utilizing of gelatin and sodium pyruvate.
Discussion
Recently, new P. larvae bacterial isolates have been reported in Austria [17] and Egypt [32]. These studies have used the molecular strategies of genotyping, biochemistry profiles or/and antibiotic susceptibility to generally investigate bacterial genetic diversity, but those methods did not improve P. larvae intraspecific discrimination [16]. In the current research, our aim was to apply multilocus sequence typing (MLST), as an improved approach to explore genetic variability and population dynamics. Development of an MLST scheme has been done in several bacterial groups such as Bacillus cereus [21], Bacillus thuringiensis [33] and Clostridium botulinum [34], and most recently in an independent effort conducted in primarily European isolates of P. larvae [23]. This study represents the first report of a broad characterization of P. larvae in North America, and doubles the set of candidate loci for applying MLST approaches for this key species. The population genetics of P. larvae were examined by multilocus sequence typing of 38 isolates by concatenating sequences of all seven polymorphic markers (ilvD, purH, tri, recF, pyeE, sucC and glpK). Here, 15 new STs were analyzed and presented by eBURST. Typically, field isolates are ST1. The others we found were rare alleles present at low frequencies (ST2-ST15). Two distinct groups were identified in the sampled population. Group 1 contained isolates from North America whereas in Group 2 reflected the reference strains obtained from European isolates (ATCC49843, LMG 16247 and LMG 16252) [15]. ST4 (233/00) and ST11 (PL-23) were singletons in this study. Our ST4 isolate derived from Sweden [15] while ST11 was obtained from Alaska (USA), suggesting geographic distinctiveness in these regions. Overall, the results discriminated the strains with high resolution and showed genetic structure across different geographical regions. This is similar to geographic separation by Margos et al. [35] using MLST analysis of the bacterium Borrelia burgdorferi in North American and European populations. Their study demonstrated that the concatenated sequences of B. burgdorferi housekeeping genes corresponded to different geographical regions. Morissey et al., 2015 [23] reported an MLST analysis of P. larvae that showed ability of this method to identify relationships between geographic distance and genetic distance amongst native and non-native countries (mostly, in Europe) and found that isolated typed as ERIC I were more diverse than others. Similarly, our reports showed more diversity in ERIC I. However, they also suggested that the ERIC II type is broadly distributed a result we did not find. Generally, P. larvae populations showed relatively little polymorphism. Two loci, purH and recF, showed significant negative values by Tajima’s test, indicating purifying selection or population expansion of these genes [28]. The others appear to be neutral housekeeping genes useful for inferring the evolutionary history of P. larvae in North America. It could be implied that ilvD1, tri, pyre, sucC and glpK are especially free from directed selection pressure, and therefore will be most useful for understanding demographic events and bacterial population history [36].
Antibiotic susceptibility may also be useful for epidemiology studies and lead to the detection of resistant strains [8]. Historically the AFB disease has been controlled either by the destruction of infected colonies (burning), or by antibiotic use. For instance, OTC and tetracycline are the same class of antibiotics. OTC is commonly used for AFB control, although resistant strains have been reported in several countries [13]. Tylosin and lincomycin are known to be effective against resistant AFB strains in the field in both Florida and New Jersey [8]. Tylosin is non-toxic to honey bees and is also used to control AFB [37]. In an earlier study, Alippi et al., 2005 [10] did not detect resistance by P. larvae toward tylosin. Likewise, lincomycin is a third alternative used to inhibit the growth of OTC resistant P. larvae [8]. Hence, in this study, we provide the first report of isolates of P. larvae showing resistance to tylosin and lincomycin. However, OTC remains the only approved drug treatment available in the United States for the prevention and control of AFB [4]. Interestingly, All OTC-resistant isolates were defined to ERIC I-BOX D genotype. However, OTC-resistant trait was not correlated with bacterial haplotype. This might be due to either resistance that has evolved multiple times in P. larvae or resistance caused by recent horizontal transfer via a non-genomic route (e.g. plasmids or conjugal transposons) [12]. In an earlier study, Alippi et al., 2005 [10] did not find tylosin resistance in this species. Hence, here we provide the first report of isolates of P. larvae showing resistance to tylosin and lincomycin.
Phenotypic comparisons using three API® system analyses of the P. larvae strains were quite contradictory and did not relate with sequence types. However, our results correspond with Neuendorf et al., 2004 [38]. In contrast, Kilwinski et al., 2004 [39] reported mannitol is solely accepted as a metabolic source by P. larvae. Our results are quite similar because all the field isolates reacted negatively for mannitol while the reference strains (ATCC 49843, LMG16247, LMG16252 and 233/00) can utilize this carbohydrate.
Analyses across the biochemical, phylogenetic, and molecular typing methods highlighted relatively high intraspecific diversity across this P. larvae collection and indicated fairly consistent grouping of isolates across all three methods. Phenetic analysis had previously been grouped and compared with genotype. Nevertheless, it seems that biochemical typing provides less discrimination than DNA fingerprinting [18]. Hamdi et al., 2013[40] suggested that relatively high intra-species molecular and phenotypic diversity did not differ significantly with in vivo virulence of P. larvae genotypes. This reflects that the genotypic/phenotypic differences are neutral or more related to ecological aspects that pathogenicity or virulence.
In conclusion, we propose the application of these seven loci for an MLST scheme for P. larvae, arguably in conjunction with those described by Morissey et al., 2015 [23]. This scheme showed higher discrimination than ERIC and BOX DNA fingerprints, whereas biochemistry test tended to be more variable among different strains studied. The observation of antibiotic susceptibility suggests the presence of new ST and antibiotic drug resistance in a limited number of P. larvae strains in North America. The data presented here using the MLST scheme strongly supports this technique as a powerful novel method for studying P. larvae population genetics and epidemiology.
Supporting information
S1 Table. Origins and typing data of Paenibacillus larvae strains analyzed in this study (1999–2013).
https://doi.org/10.1371/journal.pone.0176831.s001
(PDF)
S2 Table. Biochemistry profiles of Paenibacillus larvae by using API® tests.
https://doi.org/10.1371/journal.pone.0176831.s002
(PDF)
S3 Table. Other MLST primers of Paenibacillus larvae were tested and showed no diversitywithin the target sequences.
https://doi.org/10.1371/journal.pone.0176831.s003
(PDF)
S1 Fig. ERIC PCR genotyping of 38 Paenibacillus larvae isolates.
https://doi.org/10.1371/journal.pone.0176831.s004
(PDF)
S2 Fig. BOX PCR genotyping of 38 Paenibacillus larvae isolates.
https://doi.org/10.1371/journal.pone.0176831.s005
(PDF)
S3 Fig. Presentation of all amplified products by using designed MLST primers.
In red boxes are seven loci used in the studies.
https://doi.org/10.1371/journal.pone.0176831.s006
(PDF)
Acknowledgments
This research was financially supported by a grant under The Program of Strategic Scholarships for Frontier Research Network for the Joint Ph.D. Program Thai Doctoral degree from the Office of the Higher Education Commission, Thailand (Grant No. CHE-SW215/2551), Thailand Research Fund and Chiang Mai University fund. We wish to thank the USDA-ARS Bee Research Laboratory (Beltsville, MD, USA), and the Graduate School of Chiang Mai University. We are also grateful to M. Hamilton, D. Lopez, and B. Smith, S. Abban, Y. Melnikova and the staff of the USDA-ARS Bee Research Laboratory (Beltsville, MD, USA), for excellent technical help with the project. We also thank E. Forsgren for kindly providing isolates included in this study.
Author Contributions
- Conceptualization: JE JP PC.
- Data curation: SK JE.
- Formal analysis: SK JE.
- Funding acquisition: JE JP PC.
- Investigation: SK JE JP PC.
- Methodology: SK JE JP.
- Project administration: JE JP PC.
- Resources: JE JP PC.
- Software: SK JE.
- Supervision: JE JP PC.
- Validation: JE SK.
- Visualization: SK JE PC JP.
- Writing – original draft: SK JE JP PC.
- Writing – review & editing: SK JE JP PC.
References
- 1. Dancer BN, Chantawannakul P. American foulbrood in honey bees. Bee World. 2009; 82: 168–180.
- 2. Lindström A, Korpela S, Fries I. The distribution of Paenibacillus larvae spores in adult bees and honey and larval mortality, following the addition of American foulbrood diseased brood or spore-contaminated honey in honey bee (Apis mellifera) colonies. J Invertebr Pathol. 2008; 99:82–86. pmid:18640122
- 3. Genersch E. American foulbrood in honeybees and its causative agent, Paenibacillus larvae. J Invertebr Pathol. 2010; 103:S10–S19. pmid:19909971
- 4. Feldlaufer MF, Jeffery SP, Kochansky J, Stiles G. Lincomycin hydrochloride for the control of American foulbrood disease of honey bees. Apidologie. 2001;32: 547–554.
- 5. Fries I, Lindström A, Korpela S. Vertical transmission of American foulbrood (Paenibacillus larvae) in honey bees (Apis mellifera). Vet Microbiol. 2006; 114:269–274. pmid:16420974
- 6. Lindström A, Korpela S, Fries I. Horizontal transmission of Paenibacillus larvae spores between honey bee (Apis mellifera) colonies through robbing. Apidologie. 2008; 39:515–522.
- 7. Riessberger-Gallé U, von der Ohe W, Crailsheim K. Adult honeybee's resistance against Paenibacillus larvae larvae, the causative agent of the American foulbrood. J Invertebr Pathol. 2001; 77:231–236. pmid:11437525
- 8. Kochansky J, Knox DA, Feldlaufer MF, Jeffery SP. Screening alternative antibiotics against oxytetracycline-susceptible and -resistant Paenibacillus larvae. Apidologie. 2001; 32:215–222.
- 9. Elzen P, Westervelt D, Causey D, Rivera R, Baxter J, Feldlaufer M. Control of oxytetracycline-resistant American foulbrood with tylosin and its toxicity to honey bees (Apis mellifera). J Apicult Res. 2002; 41:97–100.
- 10. Alippi AM, Albo GN, Reynaldi FJ, De Giusti MR. In vitro and in vivo susceptibility of the honeybee bacterial pathogen Paenibacillus larvae subsp. larvae to the antibiotic tylosin. Vet Microbiol. 2005; 109:47–55. pmid:15951140
- 11. Pettis JS, Feldlaufer MF. Efficacy of lincomycin and tylosin in controlling American foulbrood in honey bee colonies. J Apicult Res. 1996; 44:106–108.
- 12. Evans JD. Diverse origins of tetracycline resistance in the honey bee bacterial pathogen Paenibacillus larvae. J Invertebr Pathol. 2003; 83(1):46–50. pmid:12725811
- 13. Miyagi T, Peng CY, Chuang RY, Mussen EC, Spivak MS, Doi RH. Verification of oxytetracycline-resistant American foulbrood pathogen Paenibacillus larvae in the United States. J Invertebr Pathol. 2000; 75(1):95–96. pmid:10631065
- 14. Debourgogne A, Gueidan C, Hennequin C, Contet-Audonneau N, de Hoog S, Machouart M. Development of a new MLST scheme for differentiation of Fusarium solani Species Complex (FSSC) isolates. J Microbiol Methods. 2010; 82:319–323. pmid:20624428
- 15. Genersch E, Forsgren E, Pentikäinen J, Ashiralieva A, Rauch S, Kilwinski J, et al. Reclassification of Paenibacillus larvae subsp. pulvifaciens and Paenibacillus larvae subsp. larvae as Paenibacillus larvae without subspecies differentiation. Int J Syst Evol Microbiol. 2006; 56:501–511. pmid:16514018
- 16. Antúnez K, Piccini C, Castro-Sowinski S, Rosado AS, Seldin L, Zunino P. Phenotypic and genotypic characterization of Paenibacillus larvae isolates. Vet Microbiol. 2007; 124(1–2):178–183. pmid:17517481
- 17. Loncaric I, Derakhshifar I, Oberlerchner JT, Köglberger H, Moosbeckhofer R. Genetic diversity among isolates of Paenibacillus larvae from Austria. J Invertebr Pathol. 2009; 100:44–46. pmid:18831978
- 18. Alippi AM, Aguilar OM. Characterization of isolates of Paenibacillus larvae subsp. larvae from diverse geographical origin by the polymerase chain reaction and BOX primers. J Invertebr Pathol. 1998; 72(1):21–27. pmid:9647697
- 19. Genersch E, Otten C. The use of repetitive element PCR fingerprinting (rep-PCR) for genetic subtyping of German field isolates of Paenibacillus larvae subsp. larvae. Apidologie. 2003; 34(3):195–206.
- 20. Alippi AM, Marioaguilar O. Unique DNA fingerprint patterns of Paenibacillus larvae subsp. larvae strains. J Apicult Res. 1998; 37:273–280.
- 21. Helgason E, Tourasse NJ, Meisal R, Caugant DA, Kolstø A. Multilocus sequence typing scheme for bacteria of the Bacillus cereus group. Appl Environ Microbiol. 2004; 70(1):91–201.
- 22. Alper I, Frenette M, Labrie S (2013) Genetic diversity of dairy Geotrichum candidum strains revealed by multilocus sequence typing. Appl Microbiol Biotechnol. 2013; 97(13):5907–20. pmid:23467823
- 23. Morrissey BJ, Helgason T, Poppinga L, Fünfhaus A, Genersch E, Budge GE. Biogeography of Paenibacillus larvae, the causative agent of American foulbrood, using a new multilocus sequence typing scheme. Environ Microbiol. 2015; 17(4):1414–24. pmid:25244044
- 24. Versalovic J, Schneider M, De Bruijn FJ, Lupski JR. Genomic fingerprinting of bacteria using repetitive sequence-based polymerase chain reaction. Methods Mol Cell Biol. 1994; 5(1):25–40.
- 25. Louws FJ, Fulbright DW, Stephens CT, DE Bruijn FJ. Specific genomic fingerprints of phytopathogenic Xanthomonas and Pseudomonas pathovars and strains generated with repetitive sequences and PCR. Appl Environ Microbiol. 1994; 60:2286–2295. PMCID: PMC201645 pmid:8074510
- 26. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013; 30(12):2725–9. pmid:24132122
- 27. Rozas J, Sánchez-DelBarrio JC, Messeguer X, Rozas R. DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics.2003; 19:2496–2497. pmid:14668244
- 28. Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989; 123(3):585–95. pmid:2513255
- 29. Fu YX, Li WH. Statistical tests of neutrality of mutations. Genetics. 1993; 133(3):693–709. pmid:8454210
- 30. Fu Y. Statistical tests of neutrality of mutations against population growth, hitch-hiking, and background selection. Genetics. 1989; 123(3):585–595. PMCID: PMC1203831
- 31. Francisco AP, Bugalho M, Ramirez M, Carriço JA. Global optimal eBURST analysis of multilocus typing data using a graphic matroid approach. BMC Bioinformatics. 2009; 10:152. pmid:19450271
- 32. Masrya SHD, Kabeilb SS, Hafeza EE. New Paenibacillus larvae bacterial isolates from honey bee colonies infected with American foulbrood disease in Egypt. Biotechnol Biotechnol Equip. 2014; 28(2):271–276. pmid:26740757
- 33. Sorokin A, Candelon B, Guilloux K, Galleron N, Wackerow-Kouzova N, Ehrlich SD, et al. Multiple-locus sequence typing analysis of Bacillus cereus and Bacillus thuringiensis reveals separate clustering and a distinct population structure of psychrotrophic strains. Appl Environ Microbiol. 2006; 72:1295–1301.
- 34. Olsen JS, Scholz H, Fillo S, Ramisse V, Lista F, Trømborg AK, et al. Analysis of the genetic distribution among members of Clostridium botulinum group I using a novel multilocus sequence typing (MLST) assay. J Microbiol Methods. 2014; 96:84–91. pmid:24246230
- 35. Margos G, Gatewood AG, Aanensen DM, Hanincová K, Terekhova D, Vollmer SA, Cornet M, Piesman J, Donaghy M, Bormane A, Hurn MA, Feil EJ, Fish D, Casjens S, Wormser GP, Schwartz I, Kurtenbach K. MLST of housekeeping genes captures geographic population structure and suggests a European origin of Borrelia burgdorferi. Proc Natl Acad Sci U S A. 2008; 105(25):8730–8735. pmid:18574151
- 36. Prajapati SK, Joshi H, Carlton JM, Rizvi MA. Neutral polymorphisms in putative housekeeping genes and tandem repeats unravels the population genetics and evolutionary history of Plasmodium vivax in India. PLoS Negl Trop Dis. 2013; 7(9):e2425. pmid:24069480
- 37. Peng YC, Mussen E, Fong A, Montague M, Tyler T. Effects of chlortetracycline of honey bee worker larvae reared in vitro. J Invertebr Pathol. 1992; 60:127–133.
- 38. Neuendorf S, Hedtke K, Tangen G, Genersch E. Biochemical characterization of different genotypes of Paenibacillus larvae subsp. larvae, a honey bee bacterial pathogen. Microbiology. 2004; 150(7):2381–90.
- 39. Kilwinski J, Peters M, Ashiralieva A, Genersch E. Proposal to reclassify Paenibacillus larvae subsp. pulvifaciens DSM 3615 (ATCC 49843) as Paenibacillus larvae subsp. larvae. Results of a comparative biochemical and genetic study. Vet Microbiol. 2004; 104(1–2):31–42. pmid:15530737
- 40. Hamdi C, Essanaa J, Sansonno L, Crotti E, Abdi K, Barbouche N, Balloi A, Gonella E, Alma A, Daffonchio D, Boudabous A, Cherif A. Genetic and biochemical diversity of Paenibacillus larvae isolated from Tunisian infected honey bee broods. Biomed Res Int. 2013; 2013:479893. pmid:24073406