Next Article in Journal
Sap-Sucking Pests; They Do Matter
Previous Article in Journal
Enhanced Myc Expression in Silkworm Silk Gland Promotes DNA Replication and Silk Production
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Insects
Volume 12
Issue 4
10.3390/insects12040362
insects-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

ERIC and WGS Typing of Paenibacillus larvae in Slovenia: Investigation of ERIC I Outbreaks

by
Alenka Žugelj
1,
Bojan Papić
2,
Irena Zdovc
2,
Urška Zajc
2,
Majda Golob
2,
Jana Avberšek
2 and
Darja Kušar
2,*
1
Unit Maribor–Ptuj, Veterinary Faculty, National Veterinary Institute, University of Ljubljana, Šentiljska Cesta 109, SI-2000 Maribor, Slovenia
2
Veterinary Faculty, Institute of Microbiology and Parasitology, University of Ljubljana, Gerbičeva 60, SI-1000 Ljubljana, Slovenia
*
Author to whom correspondence should be addressed.
Insects 2021, 12(4), 362; https://doi.org/10.3390/insects12040362
Submission received: 24 March 2021 / Revised: 16 April 2021 / Accepted: 16 April 2021 / Published: 19 April 2021
Browse Figures
Versions Notes

Abstract

:

Simple Summary

American foulbrood is a serious disease of honeybees caused by Paenibacillus larvae. ERIC-PCR is a widely used method for typing of P. larvae that currently divides it into five ERIC types (ERIC I–V); these differ in certain phenotypic characteristics—most importantly, virulence. In the first part of the study, we assessed the distribution of ERIC types in Slovenia in the period 2017–2019 on a set of 506 P. larvae isolates. We identified ERIC II as the predominant type (70.2%), followed by ERIC I (29.8%). In the second part of the study, we typed 59 outbreak-related ERIC I isolates using whole-genome sequencing, which revealed seven ERIC I-ST2 outbreak clusters (≤35 allele differences). The transmission of the outbreak clone within a 3-km radius was observed in all seven clusters and could be explained by the activity of honeybees. The transmission of the outbreak clone between geographically distant apiaries was observed in three clusters and could be explained by migratory beekeeping and trading of bee colonies. The present findings highlight the importance of beekeeping activities in the transmission of P. larvae over large geographic distances.

Abstract

Paenibacillus larvae is the causative agent of American foulbrood (AFB), a fatal disease of honeybee brood. Here, we obtained 506 P. larvae isolates originating from honey or brood samples and from different geographic regions of Slovenia in the period 2017–2019. In the first part of the study, we conducted ERIC-PCR typing to assess the frequency of ERIC types in Slovenia. Capillary electrophoresis was used for the analysis of ERIC patterns, revealing good separation efficiency and enabling easy lane-to-lane comparisons. ERIC II was the predominant type (70.2%), followed by ERIC I (29.8%); two slightly altered ERIC I banding patterns were observed but were not considered relevant for the discrimination of ERIC types. No evident spatiotemporal clustering of ERIC types was observed. To assess the clonality of the outbreak-related P. larvae ERIC I isolates, 59 isolates of this type underwent whole-genome sequencing (WGS). Whole-genome multilocus sequence typing (wgMLST) revealed seven ERIC I-ST2 outbreak clusters (≤35 allele differences) with the median intra-outbreak diversity ranging from 7 to 27 allele differences. In all seven clusters, the transmission of P. larvae outbreak clone within a 3-km radius (AFB zone) was observed, which could be explained by the activity of honeybees. In three clusters, the transmission of the outbreak clone between geographically distant apiaries was revealed, which could be explained by the activities of beekeepers such as migratory beekeeping and trading of bee colonies. The present findings reinforce the importance of beekeeping activities in the transmission of P. larvae. WGS should be used as a reference typing method for the detection of P. larvae transmission clusters.

1. Introduction

American foulbrood (AFB) is a highly contagious and globally distributed disease affecting honeybee (Apis mellifera) broods. It is caused by a spore-forming bacterium Paenibacillus larvae. In Slovenia, as in many other countries, the detection, eradication, notification and prevention of AFB are regulated by the national legislation [1]. Under this legislation, all apiaries within a 3-km radius (AFB zone) from any AFB-positive apiary are subjected to clinical examination. A new AFB zone is denoted around each newly detected AFB-positive apiary, leading to the extension of the surveillance area. All AFB-positive honeybee colonies should be destroyed, leading to considerable economic losses to the beekeeping sector. Measures for the prevention of AFB are advised when spores are identified in honey or hive debris prior to the onset of clinical symptoms such as the shook swarm method and restriction of colony trade/movement, as well as thorough and consistent disinfection of the beekeeping equipment [2].
Genotyping methods are an integral part of the surveillance of infectious diseases, including AFB. Enterobacterial repetitive intergenic consensus (ERIC)-based PCR is a widely used method for typing of P. larvae, which classifies it into five ERIC types (I–V) [3,4]. These differ in biochemical, morphological and virulence characteristics [3,4,5,6]. ERIC types I and II are globally distributed and are epidemiologically the most important types [7,8,9].
Thus far, large-scale studies on the distribution of different ERIC types have been conducted in Austria [10], Belgium [11], Italy [12] and Japan [13]. In Italy, most AFB outbreaks were caused by ERIC I, but outbreaks caused by ERIC II or coinfections of ERIC I and II were also reported. In Austria, ERIC types I and II were represented in similar frequencies [10]. In Japan, an increase in the frequency of the ERIC II-ST10 type was observed in the last decade [13]. No data are currently available on the distribution and frequency of different ERIC types in Slovenia.
Due to the limited discriminatory power of ERIC-PCR, additional discrimination has been enabled by multilocus sequence typing (MLST) [7] and recently extended to core/whole-genome MLST (cg/wgMLST) [9,14]. Compared with ERIC-PCR and MLST, whole-genome sequencing (WGS) provides an insight into the complete genetic background of an organism, reaching an unprecedented discriminatory power [15]. Although WGS is increasingly employed in the surveillance and outbreak investigations of important human and zoonotic pathogens, it has been rarely used for the investigation of AFB outbreaks [9,14,16]. Thus, additional studies are needed to understand the intra- and inter-outbreak diversity of P. larvae isolates, as well as their ecology and transmission routes.
The first aim of this study was to assess the frequency and spatiotemporal distribution of P. larvae ERIC types in Slovenia. Since several outbreak-related ERIC I isolates showed slightly altered banding patterns, the second aim was to assess the clonal relationship of ERIC I isolates using wgMLST.

2. Materials and Methods

2.1. Frequency and Spatiotemporal Distribution of ERIC Types

2.1.1. Isolate Selection

To assess the frequency and spatiotemporal distribution of ERIC types, 506 P. larvae isolates originating from all geographic regions of Slovenia in the period 2017–2019 were selected for ERIC typing. The isolates were obtained during routine AFB surveillance and were part of P. larvae collection maintained by the Institute of Microbiology and Parasitology (Veterinary Faculty, National Veterinary Institute, Slovenia). They originated from honey (n = 203) or broods (n = 303) and were obtained from 263 apiaries, which were maintained by 232 beekeepers. For all isolates, metadata on the apiary, beekeeper and isolation date were collected (Table S1). The number of isolates per apiary ranged from 1 to 16; the number of ERIC-typed isolates per apiary was increased where we observed (i) a higher density of apiaries within the AFB zone or (ii) a slightly altered ERIC I banding pattern (Figure S1).

2.1.2. Isolate Cultivation and DNA Extraction

The isolates were revived from frozen stocks (−70 °C) by streaking onto blood agar plates, followed by 72-h incubation at 37 °C. Species confirmation was performed by matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) mass spectrometry (Microflex LT; Bruker Daltonics Inc., Leipzig, Germany). For ERIC-PCR and WGS, DNA was extracted using the DNeasy Blood & Tissue kit (Qiagen, Hilden, Germany) according to the protocol for Gram-positive bacteria.

2.1.3. ERIC-PCR

ERIC typing was performed as previously described [3,17] with the slight modifications listed in reference [9]. Briefly, a 25-µL reaction was optimized to contain 1 × Multiplex PCR Master Mix (Qiagen, Hilden, Germany), 2 μM of each primer (ERIC1R and ERIC2) [18] and 5 μL of the DNA template. For the determination of ERIC patterns, amplicons were analyzed by the QIAxcel capillary electrophoresis (Qiagen, Hilden, Germany) using the QIAxcel DNA High-Resolution Kit, QX Alignment Marker 15–5000 bp, QX Size Marker 100 bp–2.5 kb, OM500 separation method and a sample injection time of 10 s. The prominent 2800-bp band present in ERIC II isolates but absent in ERIC I isolates [3,4] was considered relevant for the differentiation of ERIC types I and II.

2.1.4. Visualization of Spatiotemporal Data

The geographic location of each apiary was used to show the spatial distribution of the obtained isolates (from 1 to 16 isolates per apiary) and the isolation date to show their temporal distribution. The data were visualized using Microreact [19].

2.2. Phylogeny and Clonal Relationship of Outbreak-Related ERIC I Isolates

2.2.1. Isolate Selection for WGS Typing

To assess the phylogeny and clonal relationship between the outbreak-related ERIC I isolates, including the isolates with slightly altered ERIC I banding patterns, 59/151 isolates of ERIC type I underwent WGS (Tables S1 and S2). The isolates originated from 29 apiaries with confirmed AFB, which were maintained by 27 beekeepers and were located in 13 different AFB zones covering several geographic regions of Slovenia during the period of 2017–2019 (Table S2). Additional epidemiological data relevant for the transmission of P. larvae due to the activities of beekeepers (i.e., beekeeper, migratory beekeeping and trading of bee colonies) were collected where available (Table S2).

2.2.2. WGS Typing

The extracted DNA was quantified by Qubit 3.0 (Thermo Fisher Scientific, Waltham, MA, USA) using the Qubit 1 × dsDNA HS assay kit. DNA libraries were prepared using the NEBNext Ultra DNA Sample Prep Master Mix kit (New England Biolabs, Ipswich, MA, USA). Paired-end (2 × 150 bp) sequencing was performed on the NovaSeq 6000 system (Illumina, San Diego, CA, USA) to a minimum coverage of 130×.
Raw reads were assembled using SPAdes version 3.7.1 implemented in BioNumerics version 7.6.3 (bioMérieux, Applied Maths NV, Sint-Martens-Latem, Belgium). The wgMLST quality assessment window was used to assess the quality of the obtained reads and assemblies; only the genomes passing the applied thresholds were included in the further analyses. The following thresholds were used: (i) the average read quality of Q > 30 and the expected coverage of >30 for the trimmed reads and (ii) N50 > 21,250 kb, number of contigs <380 and the total assembly length of ~4.1 Mb (range, 3.2–5.3 Mb) for the assemblies.
The wgMLST analysis was performed in BioNumerics by applying the stable wgMLST scheme for P. larvae consisting of 5745 loci [9]. Assembly-based and assembly-free allele calling were performed by applying the default parameters. The UPGMA tree and the minimum spanning tree (MST) based on the allele profiles of 5738-wgMLST loci (excluding the partial sequences of seven MLST loci) were constructed in BioNumerics with the allele calls considered as categorical values.
To assess the phylogeny and clonal relationship of the outbreak-related isolates with slightly altered ERIC I banding patterns (38/59 isolates; Table S2), a dataset of 179 publicly available P. larvae genomes (complete or draft) described in reference [9] was used and supplemented with WGS data of the 59 ERIC I isolates obtained in this study. The UPGMA wgMLST tree was constructed as described above and visualized and annotated using iTol version 4.4.2 [20].

2.2.3. In Silico 7-Gene MLST

MLST sequence types (STs) were extracted from the assembled genomes of 59 ERIC I isolates using the MLST plugin implemented in BioNumerics; the P. larvae PubMLST nomenclature [21] based on seven previously described housekeeping genes [7] was used.

3. Results

3.1. Frequency and Spatiotemporal Distribution of ERIC Types

Among the studied 506 isolates, ERIC II was the predominant type (355/506; 70.2%), followed by ERIC I (151/506; 29.8%) (Table S1). In total, 74 ERIC I isolates exhibited a slightly altered ERIC I banding pattern; of these, 49 had an additional 600 bp and 25 an additional 650-bp band (Figure S1 and Table S1). The repeatability of these altered ERIC I banding patterns was confirmed by performing three independent ERIC-PCR replicates for one representative isolate per pattern. Figure 1 shows the spatiotemporal distribution of P. larvae ERIC types in Slovenia in 2017–2019. In general, both ERIC types were encountered in all geographic regions and over the entire study period. With regards to the sample type, ERIC II prevailed both in honey and brood samples (81.8% and 62.4%, respectively), followed by ERIC I (18.2% and 37.6%, respectively). In 107/263 apiaries where more than one isolate per apiary was typed, a mixed contamination or infection with ERIC types I and II was observed in 8/107 (7.5%) cases (Table S1). Since isolates from the same apiary represent possible epidemiological replicates, the prevalence of ERIC types in the 2017–2019 period was calculated on a set of 271 isolates originating from different apiaries and belonging to different ERIC types (Table S1). Table 1 shows the prevalence of ERIC I and ERIC II types in different years and further confirms the predominance of ERIC II type in all study years.

3.2. Phylogeny and Clonal Relationship of Outbreak-Related ERIC I Isolates

The MLST analysis of the 59 outbreak-related ERIC I isolates revealed that all of the investigated isolates belonged to ST2, including the isolates with altered the ERIC I banding patterns (38/59). The wgMLST analysis of the complete genome dataset (59 genomes of ERIC I type obtained in this study supplemented with 178 publicly available P. larvae genomes) revealed that all the investigated ERIC I genomes fell within the ERIC I clade and were clustered with the remaining ERIC I-ST2 genomes (Figure S2). The main clusters observed on the wgMLST tree of the 59 outbreak-related ERIC I isolates (Figure 2B) were generally congruent with the three observed ERIC I banding patterns (one conventional and two slightly altered), although they were not considered relevant, because they belonged to the same MLST type and clustered with the remaining ERIC I genomes in the wgMLST tree of the complete genome dataset (Figure S2).
By applying the previously described single-linkage threshold of 35 allele differences (AD) on the wgMLST MST of P. larvae [9], seven outbreak clusters were identified (Figure 2A). Five isolates did not fall within any of the identified outbreak clusters (Figure 2). The median pairwise AD ranged from 7 (outbreak cluster 2) to 27 (outbreak cluster 5) (Figure 2A) and was comparable to that observed in the previously described ERIC I-ST2 outbreak (19 AD) [9]. The minimum and maximum pairwise AD between isolates per each cluster are shown in Figure 2B. The maximum pairwise genetic distance within an outbreak cluster was 42 AD (outbreak cluster 1). The minimum distance between the outbreak-related and nonrelated isolates on the wgMLST MST was 52 AD (Figure 2A). Of note, the clinical isolate PL196 belonging to the outbreak cluster 1 according to the applied threshold could not be linked to the remaining isolates of the respective outbreak cluster according to the available epidemiological data; thus, if this isolate is considered epidemiologically and genetically unrelated, this decreases the single-linkage threshold for cluster delineation to 24 AD and the minimum distance between outbreak-related and nonrelated isolates to 34 AD (Figure 2A).

3.3. P. larvae Transmission Routes

In all the identified outbreak clusters, we observed the transmission of P. larvae outbreak clone between the apiaries that were located within the same AFB zone, which roughly corresponds to the honeybee flight distance and could thus be explained by the foraging activity of honeybees (Table S2). In addition, in outbreak clusters 1, 2 and 6, we observed the transmission of P. larvae outbreak clone over a large geographic area (i.e., between different AFB zones). For outbreak clone 2, such transmissions could be explained by migratory beekeeping (i.e., transfer of bees to a common foraging area), and for the outbreak clones 1 and 6, by the suspected (unregistered) trading of bee colonies (Table S2 and Figure 3). No intermixing of different outbreak clones was observed within a given AFB zone. Outbreak clusters 1 and 6 harbored isolates spanning the entire study period (2017–2019) and both sample types (honey and brood), and the remaining five outbreak clusters encompassed isolates from a single year (Figure 3).

4. Discussion

The present study identified ERIC II (70.2%) as the predominant ERIC type in Slovenia in 2017–2019, regardless of the sample type, followed by ERIC I (29.8%). In addition, wgMLST revealed the existence of several ERIC I-ST2 outbreak clones that could not be distinguished using conventional typing methods (ERIC-PCR and MLST). The transmission between geographically distant (>3 km) apiaries was identified in three out of seven outbreak clones, highlighting the importance of beekeeping activities in the transmission of P. larvae.

4.1. Frequency and Spatiotemporal Distribution of ERIC Types

Although ERIC II is considered to be less frequently encountered than ERIC I in Europe and worldwide [8,10,12], the present study identified ERIC II as the predominant ERIC type in Slovenia in all study years. The proportion of larvae developing into a ropy mass under cell cappings is higher for infections with ERIC I strains, which exhibit a slow-killing phenotype at the larval level. Contrary, ERIC II strains exhibit a fast-killing phenotype at the larval level; hence, false-negative diagnoses are likely to occur if AFB-diseased colonies are infected with ERIC II strains, because only a few infected cells may be present [3,5,6,22]. The high occurrence of ERIC II in Slovenia may thus result from a more successful identification and elimination of ERIC I-affected colonies as opposed to a more difficult identification of ERIC II-affected colonies, which serve as a better (and more often unrecognized) reservoir for the transmission of P. larvae spores within and between apiaries.
In this study, the country-wide spatial distribution of ERIC types was considered rather than regional distribution for the following reasons. First, the national AFB surveillance system is based on a 3-km radius around the infected apiaries rather than any of the geographic classification systems. Second, the isolate panel was biased with regards to the geographic region, because it depended on the number of AFB reported cases and provided samples. Third, the transmission of P. larvae over large geographic regions was previously reported in Slovenia [9] and was also confirmed here, which can be explained by the activities of beekeepers (e.g., migratory beekeeping, which is common in Slovenia).
In this study, capillary electrophoresis was used for the ERIC-PCR amplicon analysis, which likely resulted in the better separation/visualization and comparability of ERIC patterns in comparison with conventional agarose gel electrophoresis used in previous studies [3,4,23,24]. The existence of two slightly altered ERIC I banding patterns was revealed in 38/59 WGS-typed outbreak-related ERIC I isolates; these were confirmed as ERIC I using MLST and WGS, since they clustered with the remaining ERIC I-ST2 isolates in the wgMLST tree (see Figures S1 and S2 and Table S2 for details). Thus, the less-conserved bands migrating at 600 or 650 bp observed in certain outbreak-related ERIC I isolates were not considered relevant for the discrimination of ERIC types. Only the prominent 2800-bp band present in ERIC II but absent in ERIC I isolates was used for the differentiation between ERIC types I and II, whereas the previously described 1200-bp faint band present in ERIC I but absent in ERIC II was not observed in this study. Taken together, these findings highlighted the importance of addressing only certain relevant bands when assigning the isolates to different ERIC types. The question of whether or not to consider novel bands in ERIC profiles as relevant points out the known limitations of ERIC-PCR regarding its difficult standardization and harmonization between different laboratories [8,25]. The description of novel ERIC types should be supported by genome-wide phylogeny and include in vitro pathogenicity studies, as exemplified by the recently described ERIC V type [4], which was also shown to form a separate clade on the MLST [4] and wgMLST [9] phylogenetic trees. In addition, we also observed here that the quality of the extracted DNA affects the obtained ERIC pattern, since the cell lysates of the P. larvae isolates proved unsuitable for ERIC-PCR (data not shown).
Since all the outbreak-related isolates were of the ERIC I-ST2 type, this study confirmed the limited discriminatory power of ERIC-PCR and MLST, both of which failed to discriminate between the identified outbreak clones. Although MLST was shown to provide a higher discriminatory power compared with ERIC-PCR [7], it could not differentiate between the identified outbreak clusters. Due to its unprecedented discriminatory power, WGS should be used as a reference method to detect P. larvae transmission clusters, as already shown in previous studies [9,16].

4.2. Phylogeny, Clonal Relationship and Transmission of Outbreak-Related ERIC I Isolates

The wgMLST analysis enabled a clear delineation of seven ERIC I-ST2 outbreak clusters. The median intra-outbreak diversity of P. larvae observed in this study was comparable to that reported in previous studies [9,14]. The present study reinforces the previously proposed 35-allele threshold (as determined by wgMLST) to delineate the outbreak clusters [9] but should allow for some flexibility due to technical and biological variations, as extensively discussed elsewhere [9,14,26,27]. In the present study, BioNumerics was used in combination with a stable wgMLST scheme for P. larvae [9], facilitating the harmonization and standardization of WGS typing tools and data interpretation in order to improve the surveillance of this important honeybee pathogen. The identified wgMLST clusters can be further investigated by a whole-genome single-nucleotide polymorphism (wgSNP) analysis, which may provide additional confirmation of the identified clusters. We previously showed that allele-based and wgSNP-based approaches both gave congruent results with regards to the outbreak cluster delineation of P. larvae isolates [9].
In outbreak clusters 1 and 6, the outbreak clone was observed throughout the entire study period (2017–2019), suggesting its persistence within the AFB zone(s). This could be explained by the high resilience of P. larvae spores and insufficient/inadequate preventive measures [8]. Moreover, three out of seven outbreak clones were not limited to a single AFB zone, suggesting the transmission of P. larvae outbreak clones over larger geographic distances. The latter most likely results from beekeeping activities, which may remain overlooked under the current national legislation that focuses on the examination of apiaries within the AFB zone. Taken together, the present findings call for the implementation of an improved country-wide surveillance protocol for the early diagnosis of AFB (e.g., monitoring of honey and hive debris samples) to detect P. larvae-positive apiaries without clinical symptoms of AFB, followed by rigorous prevention and eradication measures in these apiaries, which could decrease the frequency of such recurrent outbreaks. In this manner, apiaries with an increased risk of developing AFB would be identified regardless of their location from the current AFB zone(s).
Genotyping results should be interpreted in combination with epidemiological data for the final confirmation of outbreaks and transmission events. In this study, the transmission of P. larvae outbreak clones within the AFB zone could be explained by the activity of honeybees (e.g., foraging, robbing and drifting). The transmission between geographically distant (>3 km) apiaries, which was observed in three outbreak clones, could be explained by migratory beekeeping (transfer of honeybee colonies to a common foraging area) or suspected trading of colonies that was not reported to the authority. In Slovenia, all migratory beekeeping activities should be reported to the national authority; nonetheless, the possibility of illegal (unreported) transfer of bee colonies cannot be ruled out. The present findings restate the importance of collecting extensive epidemiological data on the adopted beekeeping technologies and practices, as well as the participation of beekeepers themselves, to report migratory beekeeping activities and clinical symptoms of AFB in order to control the disease.

5. Conclusions

ERIC II was recognized as the predominant ERIC type in Slovenia. WGS typing enabled a clear delineation of the ERIC I-ST2 outbreak clusters. In addition to the transmission of P. larvae outbreak clone within the AFB zone, several cases of transmission between geographically distant apiaries were revealed, indicating the role of beekeeping activities in the transmission of P. larvae. The present study improves our understanding of the genetic diversity and transmission dynamics of P. larvae, facilitating future surveillance of this important bee pathogen. WGS data should be interpreted in combination with epidemiological data for final confirmation of transmission events.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/insects12040362/s1: Figure S1. ERIC patterns of the analyzed Paenibacillus larvae isolates obtained by QIAxcel (Qiagen, Hilden, Germany) capillary electrophoresis. ERIC types are denoted above the patterns. Green arrows denote the bands that were considered important for the differentiation of ERIC types, whereas red arrows denote the bands that occurred in certain outbreak-related ERIC I isolates but were not considered important for the differentiation of ERIC types [3,4]. Columns 1 and 2: P. larvae outbreak-related field strains exhibiting a slightly altered ERIC I banding pattern. Columns 3–7 (Control strains): P. larvae strains ATCC 9545 (ERIC I), CCUG 49660A (ERIC II), LMG 16252 (ERIC III), LMG 15974 (ERIC IV) and DSM 106052 (ERIC V). Column 8 (M): QX DNA Size Marker 100 bp–2.5 kb (Qiagen, Hilden, Germany). See Table S1, column “Comment (ERIC pattern)”, for details. Figure S2. Whole-genome phylogenetic tree of 238 Paenibacillus larvae genomes. The UPGMA wgMLST tree was constructed based on 5738 loci and is annotated with ERIC and MLST types. The genomes are colored according to ERIC type. The 59 outbreak-related ERIC I isolates analyzed in this study are denoted by red stars. For genome metadata, see Table S2. Table S1. Metadata of the 506 Paenibacillus larvae isolates from Slovenia used for assessing the frequency and spatiotemporal distribution of ERIC types. Table S2. Metadata of the 59 outbreak-related Paenibacillus larvae isolates of ERIC I-ST2 type analyzed in this study. NA, not applicable; ND, no data.

Author Contributions

Conceptualization, A.Ž., B.P. and D.K.; Funding acquisition, D.K.; Methodology, A.Ž., B.P., I.Z., U.Z., M.G., J.A. and D.K.; Project administration, D.K.; Visualization, B.P. and D.K.; Writing—original draft, A.Ž., B.P. and D.K. and Writing—review and editing, A.Ž., B.P., I.Z., U.Z., M.G., J.A. and D.K. All authors have read and agreed to the published version of the manuscript.

Funding

The authors acknowledge financial support from the Slovenian Research Agency (“Genotyping of the American foulbrood causative agent as a tool to determine the transmission routes of the disease”—project no. V4-1804; “Animal health, environment and food safety”—research core funding no. P4-0092) and the Ministry of Agriculture, Forestry and Food of the Republic of Slovenia (project no. V4-1804).

Institutional Review Board Statement

Not applicable.

Data Availability Statement

All WGS data obtained in this study have been deposited in the NCBI Sequence Read Archive (SRA) database under BioProject accession number PRJNA648798.

Acknowledgments

The authors thank Jerica Vreček Šulgaj, Lucija Žvokelj, Mira Jenko Rogelj, Ivo Planinc, Suzana Skerbiš, Anita Vraničar Novak, Martina Žagar, Alenka Jurič, Vida Lešnik, Mateja Ratiznojnik for sampling and Metka Pislak Ocepek for coordination of honeybee disease surveillance (all Veterinary Faculty, National Veterinary Institute University of Ljubljana). Maja Lepen, Olga Kosar, Maja Kavalič, Marjeta Jarc, Biljana Hacin, Jožef Lugarič, Andreja Šinko, Brigita Kuhar, Mojca Čelan, Mirjana Korošak, Darja Murko, Valerija Divjak and Majda Fišer (all Veterinary Faculty, National Veterinary Institute, University of Ljubljana) are acknowledged for microbiological analysis of the samples. Maja Lepen and Olga Kosar are also acknowledged for their excellent contribution to the strain database management.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pravno-Informacijski Sistem Republike Slovenije. Pravilnik o Ukrepih za Ugotavljanje, Zatiranje, Obveščanje in Preprečevanje Hude Gnilobe čEbelje Zalege; Pravno-Informacijski Sistem Republike Slovenije: Ljubljana, Slovenia, 2020; p. 5. Available online: http://www.pisrs.si/Pis.web/pregledPredpisa?id=PRAV13788 (accessed on 19 March 2021). (In Slovenian)
  2. Food and Agriculture Organization of the United Nations. Good Beekeeping Practices: Practical Manual on How to Identify and Control the Main Diseases of the Honeybee (Apis mellifera): TECA—Technologies and Practices for Small Agricultural Producers; Food and Agriculture Organization of the United Nations: Rome, Italy, 2020; pp. 39–43. Available online: http://www.izslt.it/wp-content/uploads/2020/05/CA9182EN.pdf (accessed on 10 February 2021).
  3. Genersch, E.; Forsgren, E.; Pentikäinen, J.; Ashiralieva, A.; Rauch, S.; Kilwinski, J.; Fries, I. 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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Beims, H.; Bunk, B.; Erler, S.; Mohre, K.I.; Spröer, C.; Silke, C.S.; Pradella, S.; Günther, G.; Rohde, M.; von der Ohe, W.; et al. Discovery of Paenibacillus larvae ERIC V: Phenotypic and genomic comparison to genotypes ERIC I-IV reveal different inventories of virulence factors which correlate with epidemiological prevalences of American Foulbrood. Int. J. Med. Microbiol. 2020, 310, e151394. [Google Scholar] [CrossRef] [PubMed]
  5. Genersch, E.; Ashiralieva, A.; Fries, I. Strain- and genotype-specific differences in virulence of Paenibacillus larvae subsp. larvae, a bacterial pathogen causing American foulbrood disease in Honeybees. Appl. Environ. Microbiol. 2005, 71, 7551–7555. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Rauch, S.; Ashiralieva, A.; Hedtke, K.; Genersch, E. Negative correlation between individual-insect-level virulence and col-ony-level virulence of Paenibacillus larvae, the etiological agent of American foulbrood of honey bees. Appl. Environ. Microbiol. 2009, 75, 3344–3347. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Morrissey, B.J.; Helgason, T.; Poppinga, L.; Fünfhaus, A.; Genersch, E.; Budge, G.E. Biogeography of Paenibacillus larvae, the causative agent of American foulbrood, using a new multilocus sequence typing scheme. Environ. Microbiol. 2014, 17, 1414–1424. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. World Organization for Animal Health. Chapter 3.2.2.: American foulbrood of honey bees (infection of honey bees with Paenibacillus larvae). In OIE Terrestrial Manual, 8th ed.; World Organization for Animal Health: Paris, France, 2018; pp. 719–735. Available online: https://www.oie.int/fileadmin/Home/eng/Health_standards/tahm/3.02.02_AMERICAN_FOULBROOD.pdf (accessed on 10 March 2021).
  9. Papić, B.; Diricks, M.; Kušar, D. Analysis of the global population structure of Paenibacillus larvae and outbreak investigation of American foulbrood using a stable wgMLST scheme. Front. Vet. Sci. 2021, 8, 582677. [Google Scholar] [CrossRef] [PubMed]
  10. Loncaric, I.; Derakhshifar, I.; Oberlerchner, J.T.; Köglberger, H.; Moosbeckhofer, R. Genetic diversity among isolates of Pae-nibacillus larvae from Austria. J. Invertebr. Pathol. 2009, 100, 44–46. [Google Scholar] [CrossRef]
  11. Descamps, T.; De Smet, L.; Stragier, P.; De Vos, P.; De Graaf, D.C. Multiple locus variable number of tandem repeat analysis: A molecular genotyping tool for Paenibacillus larvae. Microb. Biotechnol. 2016, 9, 772–781. [Google Scholar] [CrossRef] [Green Version]
  12. Bassi, S.; Formato, G.; Milito, M.; Trevisiol, K.; Salogni, C.; Carra, E. Phenotypic characterization and ERIC-PCR based gen-otyping of Paenibacillus larvae isolates recovered from American foulbrood outbreaks in honey bees from Italy. Vet. Q. 2015, 35, 27–32. [Google Scholar] [CrossRef]
  13. Ueno, Y.; Yoshida, E.; Misumi, W.; Watando, E.; Suzuki, K.; Hirai, Y.; Okura, M.; Osaki, M.; Katsuda, K.; Takamatsu, D. Population structure and antimicrobial susceptibility of Paenibacillus larvae isolates from American foulbrood cases in Apis mellifera in Japan. Environ. Microbiol. Rep. 2018, 10, 210–216. [Google Scholar] [CrossRef]
  14. Ågren, J.; Schäfer, M.O.; Forsgren, E. Using whole genome sequencing to study American foulbrood epidemiology in honey-bees. PLoS ONE 2017, 12, e0187924. [Google Scholar] [CrossRef] [PubMed]
  15. Uelze, L.; Grützke, J.; Borowiak, M.; Hammerl, J.A.; Juraschek, K.; Deneke, C.; Tausch, S.H.; Malorny, B. Typing methods based on whole genome sequencing data. One Health Outlook 2020, 2, 1–19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Bertolotti, A.C.; Forsgren, E.; Schäfer, M.O.; Sircoulomb, F.; Gaïani, N.; Ribière-Chabert, M.; Paris, L.; Lucas, P.; de Boisséson, C.; Skarin, J.; et al. Development and evaluation of a core genome multilocus sequence typing scheme for Paenibacillus larvae, the deadly American foulbrood pathogen of honeybees. Environ. Microbiol. 2021. [Google Scholar] [CrossRef]
  17. 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, 195–206. [Google Scholar] [CrossRef] [Green Version]
  18. Versalovic, J.; Schneider, M.; De Bruijn, F.J.; Lupski, J.R. Genomic fingerprinting of bacteria using repetitive sequence-based polymerase chain reaction. Methods Mol. Cell. Biol. 1994, 5, 25–40. [Google Scholar]
  19. Argimón, S.; Abudahab, K.; Goater, R.J.E.; Fedosejev, A.; Bhai, J.; Glasner, C.; Feil, E.J.; Holden, M.T.G.; Yeats, C.A.; Grund-mann, H.; et al. Microreact: Visualizing and sharing data for genomic epidemiology and phylogeography. Microb. Genomics 2016, 2, e000093. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Letunic, I.; Bork, P. Interactive Tree of Life (iTOL) v4: Recent updates and new developments. Nucleic Acids Res. 2019, 47, 256–259. [Google Scholar] [CrossRef] [Green Version]
  21. Jolley, K.A.; Bray, J.E.; Maiden, M.C.J. Open-access bacterial population genomics: BIGSdb software, the PubMLST.org website and their applications. Wellcome Open Res. 2018, 3, 124. [Google Scholar] [CrossRef]
  22. Genersch, E. American Foulbrood in honeybees and its causative agent, Paenibacillus larvae. J. Invertebr. Pathol. 2010, 103, 10–19. [Google Scholar] [CrossRef] [PubMed]
  23. Dingman, D.W. Comparative analysis of Paenibacillus larvae genotypes isolated in Connecticut. Arch. Microbiol. 2015, 197, 789–795. [Google Scholar] [CrossRef]
  24. Ebeling, J.; Knispel, H.; Fünfhaus, A.; Genersch, E. The biological role of the enigmatic C3larvinAB toxin of the honey bee pathogenic bacterium Paenibacillus larvae. Environ. Microbiol. 2019, 21, 3091–3106. [Google Scholar] [CrossRef] [PubMed]
  25. Rusenova, N.; Parvanov, P.; Stanilova, S. Molecular typing of Paenibacillus larvae strains isolated from Bulgarian apiaries based on repetitive element polymerase chain reaction (Rep-PCR). Curr. Microbiol. 2013, 66, 573–577. [Google Scholar] [CrossRef] [PubMed]
  26. Chen, Y.; Luo, Y.; Carleton, H.; Timme, R.; Melka, D.; Muruvanda, T.; Wang, C.; Kastanis, G.; Katz, L.S.; Turner, L.; et al. Whole genome and core genome multilocus sequence typing and single nucleotide polymorphism analyses of Listeria monocytogenes isolates associated with an outbreak linked to cheese, United States, 2013. Appl. Environ. Microbiol. 2017, 83, e00633–e00717. [Google Scholar] [CrossRef] [Green Version]
  27. Pightling, A.W.; Pettengill, J.B.; Luo, Y.; Baugher, J.D.; Rand, H.; Strain, E. Interpreting whole-genome sequence analyses of foodborne bacteria for regulatory applications and outbreak investigations. Front. Microbiol. 2018, 9, 1482. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Insects 12 00362 g001a 550Insects 12 00362 g001b 550
Figure 1. Spatial (A) and temporal (B) distribution of ERIC types in Slovenia. A total of 506 Paenibacillus larvae isolates from the 2017–2019 period were included into the analysis. Each circle represents an apiary (A) or isolate (B) and is colored according to ERIC type. Note that more than one isolate per apiary was obtained in 107 cases.
Figure 1. Spatial (A) and temporal (B) distribution of ERIC types in Slovenia. A total of 506 Paenibacillus larvae isolates from the 2017–2019 period were included into the analysis. Each circle represents an apiary (A) or isolate (B) and is colored according to ERIC type. Note that more than one isolate per apiary was obtained in 107 cases.
Insects 12 00362 g001aInsects 12 00362 g001b
Insects 12 00362 g002 550
Figure 2. Whole-genome MLST analysis of 59 outbreak-related Paenibacillus larvae isolates of the ERIC I-ST2 type. Each isolate is colored according to the corresponding P. larvae outbreak cluster; the isolates without an assigned outbreak cluster are denoted “not assigned” (NA). The wgMLST analysis was based on 5738 wgMLST loci. (A) wgMLST minimum spanning tree. The connecting line lengths are proportional to the number of allele differences (AD) between isolates, which are denoted on the lines. Closely related isolates belonging to the same outbreak cluster (≤35 AD) are shaded in gray. Isolate PL196 is highlighted (for details, see Section 3.2.). (B) UPGMA wgMLST tree. Scale bar indicates the number of AD. The median pairwise AD between the isolates of the same outbreak cluster is shown on the right of the tree, with the minimum and maximum AD in parenthesis. *, isolates exhibiting a conventional ERIC I banding pattern; **, isolates exhibiting a slightly altered ERIC I pattern with an additional 600-bp band; ***, isolates exhibiting a slightly altered ERIC I pattern with an additional 650-bp band. For isolate metadata, see Table S2.
Figure 2. Whole-genome MLST analysis of 59 outbreak-related Paenibacillus larvae isolates of the ERIC I-ST2 type. Each isolate is colored according to the corresponding P. larvae outbreak cluster; the isolates without an assigned outbreak cluster are denoted “not assigned” (NA). The wgMLST analysis was based on 5738 wgMLST loci. (A) wgMLST minimum spanning tree. The connecting line lengths are proportional to the number of allele differences (AD) between isolates, which are denoted on the lines. Closely related isolates belonging to the same outbreak cluster (≤35 AD) are shaded in gray. Isolate PL196 is highlighted (for details, see Section 3.2.). (B) UPGMA wgMLST tree. Scale bar indicates the number of AD. The median pairwise AD between the isolates of the same outbreak cluster is shown on the right of the tree, with the minimum and maximum AD in parenthesis. *, isolates exhibiting a conventional ERIC I banding pattern; **, isolates exhibiting a slightly altered ERIC I pattern with an additional 600-bp band; ***, isolates exhibiting a slightly altered ERIC I pattern with an additional 650-bp band. For isolate metadata, see Table S2.
Insects 12 00362 g002
Insects 12 00362 g003 550
Figure 3. Map of the investigated AFB outbreaks caused by Paenibacillus larvae ERIC I-ST2 clones, Slovenia, 2017–2019. The figure shows spatial (A) and temporal (B) distributions of the 59 outbreak-related ERIC I isolates. Each circle represents an apiary (A) or isolate (B) and is colored according to the corresponding P. larvae outbreak cluster; the isolates without an assigned outbreak cluster are denoted “not assigned” (NA). Epidemiological data supporting the transmission of the outbreak clone over large geographic areas are highlighted (A). For isolate metadata, see Table S2.
Figure 3. Map of the investigated AFB outbreaks caused by Paenibacillus larvae ERIC I-ST2 clones, Slovenia, 2017–2019. The figure shows spatial (A) and temporal (B) distributions of the 59 outbreak-related ERIC I isolates. Each circle represents an apiary (A) or isolate (B) and is colored according to the corresponding P. larvae outbreak cluster; the isolates without an assigned outbreak cluster are denoted “not assigned” (NA). Epidemiological data supporting the transmission of the outbreak clone over large geographic areas are highlighted (A). For isolate metadata, see Table S2.
Insects 12 00362 g003
Table
Table 1. Prevalence of the ERIC types over time. The prevalence was calculated on a set of 271 Paenibacillus larvae isolates originating from different apiaries and belonging to different ERIC types (for details, see Table S1).
Table 1. Prevalence of the ERIC types over time. The prevalence was calculated on a set of 271 Paenibacillus larvae isolates originating from different apiaries and belonging to different ERIC types (for details, see Table S1).
YearERIC I (%)ERIC II (%)Total No. of Isolates
201744.255.886
201814.585.583
201926.573.5102
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Žugelj, A.; Papić, B.; Zdovc, I.; Zajc, U.; Golob, M.; Avberšek, J.; Kušar, D. ERIC and WGS Typing of Paenibacillus larvae in Slovenia: Investigation of ERIC I Outbreaks. Insects 2021, 12, 362. https://doi.org/10.3390/insects12040362

AMA Style

Žugelj A, Papić B, Zdovc I, Zajc U, Golob M, Avberšek J, Kušar D. ERIC and WGS Typing of Paenibacillus larvae in Slovenia: Investigation of ERIC I Outbreaks. Insects. 2021; 12(4):362. https://doi.org/10.3390/insects12040362

Chicago/Turabian Style

Žugelj, Alenka, Bojan Papić, Irena Zdovc, Urška Zajc, Majda Golob, Jana Avberšek, and Darja Kušar. 2021. "ERIC and WGS Typing of Paenibacillus larvae in Slovenia: Investigation of ERIC I Outbreaks" Insects 12, no. 4: 362. https://doi.org/10.3390/insects12040362

APA Style

Žugelj, A., Papić, B., Zdovc, I., Zajc, U., Golob, M., Avberšek, J., & Kušar, D. (2021). ERIC and WGS Typing of Paenibacillus larvae in Slovenia: Investigation of ERIC I Outbreaks. Insects, 12(4), 362. https://doi.org/10.3390/insects12040362

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop