Effect of Bodily Fluids from Honey Bee (Apis mellifera) Larvae on Growth and Genome-Wide Transcriptional Response of the Causal Agent of American Foulbrood Disease (Paenibacillus larvae)

Abstract

Paenibacillus larvae, the causal agent of American Foulbrood disease (AFB), affects honey bee health worldwide. The present study investigates the effect of bodily fluids from honey bee larvae on growth velocity and transcription for this Gram-positive, endospore-forming bacterium. It was observed that larval fluids accelerate the growth and lead to higher bacterial densities during stationary phase. The genome-wide transcriptional response of in vitro cultures of P. larvae to larval fluids was studied by microarray technology. Early responses of P. larvae to larval fluids are characterized by a general down-regulation of oligopeptide and sugar transporter genes, as well as by amino acid and carbohydrate metabolic genes, among others. Late responses are dominated by general down-regulation of sporulation genes and up-regulation of phage-related genes. A theoretical mechanism of carbon catabolite repression is discussed.

About the Authors:

Lina De Smet

Contributed equally to this work with: Lina De Smet, Dieter De Koker

Affiliation: Ghent University, Laboratory of Zoophysiology, Department of Physiology, Ghent, Belgium

Dieter De Koker

Contributed equally to this work with: Lina De Smet, Dieter De Koker

† Deceased.

Affiliation: Ghent University, Laboratory of Zoophysiology, Department of Physiology, Ghent, Belgium

Alyse K. Hawley

Affiliation: University of British Columbia, Department of Microbiology & Immunology, Vancouver, Canada

Leonard J. Foster

Affiliation: University of British Columbia, Department of Biochemistry & Molecular Biology, Vancouver, Canada

Paul De Vos

Affiliation: Ghent University, Laboratory of Microbiology, Department of Biochemistry and Microbiology, Ghent, Belgium

Dirk C. de Graaf

* E-mail: [email protected]

Affiliation: Ghent University, Laboratory of Zoophysiology, Department of Physiology, Ghent, Belgium

Introduction

American Foulbrood (AFB) is a devastating brood disease affecting honey bee health worldwide [1]. The causal agent of AFB is the endospore-forming, Gram-positive bacterium Paenibacillus larvae. One-day-old bee larvae are most susceptible to infection [2], requiring less than ten orally ingested spores to establish a lethal infection [3]. After germination of ingested spores, vegetative P. larvae bacteria multiply within the gut lumen and subsequently breach the gut epithelium to reach the hemocoel [4]. Honey bee larvae are then degraded to brownish glue-like remains (ropy stage), which thereafter dries down to hard foulbrood scales containing approximately 2.5 billion spores [5].

In most countries AFB is a notifiable disease and major differences exist in imposed containment strategies. Until the late nineties, many countries followed an eradication procedure with destruction by burning. This drastic act seemed necessary for annihilation of the persistent and resistant spores [6], [7], until it was proven that large discrepancies exist between the number of outbreaks and the spread of spores [8]. Some countries allow administration of antibiotics. This practice, however, does not affect the infectious spore stage and leads to antibiotic resistance [9] and contamination in the honey [10]. More recently there is a tendency to control the disease by regular monitoring and elimination of symptomatic colonies. A fair amount of research focuses on alternative strategies [11], such as antagonistic bacteria [12], plant extracts [13], propolis [14], fatty acids [15], gamma radiation [16] and shook-swarm [17] treatments. The development of innovative treatment regimens will benefit from molecular knowledge on host-pathogen interactions during infection [18]. Such research is facilitated by the availability of the P. larvae genome sequence [19], [20] and it is within this conceptual framework that the present study was undertaken.

More specifically, this study investigates the growth velocity and genome-wide gene expression for P. larvae, by microarray technology, associated with in vitro exposure to bodily fluids from honey bee larvae. This approach was inspired by previous microarray studies reporting transcriptional changes for (pathogenic) bacteria in response to host material [21]–[23]. More precisely, microarray analysis was used to identify genes that were differentially expressed when potato extracts were added to the growth medium from Pectobacterium atrosepticum the causing agent of black leg disease in potato. Interestingly, some of the identified genes encoded virulence determinants.

Here we have examined (1) if honey bee larval bodily fluids could change P. larvae gene expression - e.g. overexpressing of genes encoding virulence factors - and (2) if these transcriptional changes correlated with altered P. larvae growth phenotypes.

Materials and Methods

Preparation of Honey Bee Larval Bodily Fluids

Larval bodily fluids were prepared by pooling third, fourth and fifth instars of the honey bee, Apis mellifera carnica, collected from the experimental beekeeping facility (Department of Physiology, Laboratory of Zoophysiology, University of Ghent, Belgium). Phenylthiourea (Hopkins & Williams Ltd) [100 µg per ml phosphate buffered saline (PBS, 20 mM KH2PO4, 60 mM Na2HPO4 and 145 mM NaCl, pH 7.2)], an inhibitor of phenoloxidase, was added immediately after squeezing the larvae to a final concentration of 10 µg/ml to prevent melanization. The homogenate was centrifuged twice for 30 min at 75,600×g and 4°C with an Avanti J30-I centrifuge (Beckman Coulter, Inc.). The clarified supernatant was filter sterilized with a 0.2 µm filter (Whatman) to render the bodily fluid for P. larvae culture spiking.

P. Larvae Strain and Genotyping

P. larvae strain BRL-230010 was kindly provided by Dr. Queenie Chan (University of British Columbia, Canada). This was isolated from scales collected from a single severely diseased colony in Berkely, CA, USA. The strain was genotyped as described in [24].

Effect of Larval Bodily Fluids on P. Larvae in vitro Growth

Throughout the experiments P. larvae strain BRL-230010 was routinely grown on brain heart infusion with thiamine (BHIT) broth at 37°C with agitation on an orbital shaker operating at 200 rpm [25]. Liquid bacterial cultures were started by inoculating BHIT broth with one P. larvae colony, grown for three days on BHIT agar at 37°C.

Examination of the potential (phenotypical) effect of host material on P. larvae growth (as a function of time) was first carried out by spiking bacterial cultures (in test tubes) at an OD590 of 0.2 with different concentrations of larval fluids. The same volume (1 ml) of different dilutions of larval bodily fluids was added to 1 ml culture to obtain respectively an 1/10, 1/25, 1/50, 1/100, 1/250 and 1/500 dilution. Each dilution was tested in triplicate. Addition of BHIT (with PTU in PBS) to liquid P. larvae cultures, instead of larval liquids, always served as a negative control. Addition of bodily fluids didn’t influence the absorbance at 590 nm.

Theoretical growth curves were modeled with DMFit [26]. Values for ‘rate’ and ‘yEnd’ were used as estimators for bacterial growth during exponential and stationary phase, respectively, and statistically analyzed with Kruskal-Wallis and Dunn’s test by GraphPad Prism (GraphPad Software, Inc.). The latter took only the comparisons for the 1/10 and 1/25 fluid dilution versus a control into account.

Reference Gene Selection and Validation of Microarray Data by RT-qPCR

Growth and harvest of P. larvae cells.

Reference gene stability was determined for 66 independent P. larvae 3 ml cultures (in test tubes), grown at the above-mentioned conditions, until they reached an OD595 of 0.2. At this density, 33 cultures were spiked with 4% (final concentration) larval bodily fluids (t, test) and 33 with BHIT broth (c, control). Eleven test cultures were further incubated for one hour (11×t1), eleven others for three hours (11×t3) and the remaining eleven for nine hours (11×t9). The same incubation regime was dictated to the control cultures (11×c1, 11×c3, 11×c9), thereby creating six different conditions (t1, t3, t9, c1, c3, c9). Bacterial cells were collected by centrifuging for 5 min at 8720×g and 4°C. The bacterial pellets were immediately resuspended in 650 µl RNAlater Solution (Ambion). Aliquots (50 µl each) of these suspensions were incubated on ice for 30 min, centrifuged for 5 min at 6708×g and 4°C and stored at −20°C until RNA preparation.

RNA preparation.

Prior to RNA extraction thawed aliquots were centrifuged for 2 min at 6708×g and 4°C. Bacterial pellets were resuspended by vortexing in 100 µl TE buffer (10 mM Tris-HCl, 1 mM EDTA, pH 8.0), containing lysozyme (15 mg/ml). The samples were incubated for 10 min at room temperature with constant shaking.

RNA was isolated with the InviTrap Spin Cell RNA Mini Kit (Invitek), using the protocol “Total RNA extraction from Gram-positive or Gram-negative bacteria” provided by the manufacturer. Additionally, an on-column DNase I treatment with the RNase-free DNase set (Qiagen) was performed. RNA was eluted in 40 µl elution buffer and stored at −20°C.

cDNA synthesis. RNA (5 µg) was converted to cDNA using random primers with the RevertAid First Strand cDNA Synthesis Kit (Fermentas), according to the manufacturer’s instructions.

Primer design and secondary structure formation of amplicon.

Primers for nine candidate reference genes (Table S1 in File S1), with product size-ranges of 80 to 150 bp, were designed with Primer3Plus [27], using the default settings. Amplicon secondary structures were evaluated with MFold [28], using the default settings except for the folding temperature (60°C) and ionic concentrations ([Na+] = 50 mM, [Mg2+] = 3 mM).

RT-qPCR reaction mixture.

For the RT-qPCR assay the reaction Platinum SYBR Green qPCR SuperMix-UDG (Invitrogen) kit was used. Each 15 µl reaction consisted of 7.5 µl SYBR master mix, 0.2 µM forward and 0.2 µM reverse primers (Integrated DNA Technologies) and 1 µl cDNA template using the CFX96 Real-Time PCR Detection System (Bio-Rad). The PCR program comprises a UDG digestion step of 2 min at 50°C, an activation step of 2 min at 95°C and 40 cycles of a combined denaturation (20 sec at 95°C) and annealing (40 sec at 60°C) step. At the end of this program a melt curve is generated by measuring fluorescence after each temperature increase of 0.5°C for 5 sec over a range from 65°C to 95°C. Primer efficiencies, R2 values and melt curves were calculated with CFX Manager Software (Bio-Rad). Reference gene stability was analyzed with the geNormPLUS algorithm within the qBasePLUS environment (Biogazelle NV). Default settings were kept, except that target specific amplification efficiencies were used. Differential gene expression of twenty target genes (Table S2 in File S1) was statistically assessed with qBasePLUS [29], [30], by means of unpaired t-tests (T4–C4 and T4–T1) or Mann-Whitney U tests (T1–C1 and C4–C1).

Sample Preparation for Microarray and Validation Experiment

Growth and harvest of P. larvae cells.

Samples for the microarray and the validation (by RT-qPCR) experiment were prepared by growing 16 independent 20 ml P. larvae cultures in 250 ml Erlenmeyer flasks, under the above mentioned conditions. The conditions were slightly changed in comparison with previous experiments in order to obtain enough bacterial cells to perform the microarray and validation experiment. Eight cultures were spiked with 5% (final concentration) larval bodily fluids (T, test) and eight with BHIT broth (C, control) at OD595 of 0.3. Four of the test and control cultures were incubated for 1 h (4× T1, and 4× C1) and the others for 4 h (4× T4, 4× C4). The bacterial cells were collected as described for reference gene selection and stored at −80°C until RNA preparation.

An accessory evaluation of P. larvae growth promotion in response to larval fluids was executed by comparing OD590 measurements on all bacterial cultures subjected to microarray analysis (16 samples: 4× C1, 4× T1, 4× C4, 4× T4). Comparisons T1–C1 and T4–C4 were statistically evaluated with a Mann-Whitney U test (p-value threshold = 0.05), using GraphPad Prism (GraphPad Software, Inc.).

RNA preparation.

Bacterial pellets were resuspended in 200 µl TE buffer (30 mM Tris-HCl, 1 mM EDTA, pH 8.0) with additional proteinase K (Qiagen; 1/10 dilution). RNA was isolated with the RNeasy Plus Mini Kit, using the protocol “Purification of Total RNA for Animal Cells” and performing the optional on-column DNase I treatment, according to the manufacturer’s instructions. The RNA pellet was dissolved in 30 µl RNase-free water and stored at −80°C until use.

Microarray Study on P. Larvae Transcriptional Response to Bee Bodily Fluids

Microarray experimental procedures.

A custom 8×15 k Agilent array for P. larvae was developed with Agilent eArray software. RNA concentration and purity were determined spectrophotometrically using the Nanodrop ND-1000 (Nanodrop Technologies) and RNA integrity was assessed using a Bioanalyser 2100 (Agilent). For each sample, 5 µg of total RNA, spiked with 10 viral polyA transcript controls (Agilent), was converted to single stranded cDNA. The sample was subsequently labeled with Cyanine 3 (Cy3) mono-reactive dye or Cyanine 5 (Cy5) mono-reactive dye (GE Healthcare) according to the manufacturer’s protocol (two-color microarray-based prokaryote analysis (Fairplay III labeling) - Agilent). A mixture of purified and labeled cDNA (Cy3 label: 300 ng; Cy5 label: 300 ng) was hybridised followed by (manual) washing, according to the manufacturer’s procedures. To assess the raw probe signal intensities, arrays were scanned using the Agilent DNA MicroArray Scanner with SureScan High-Resolution Technology and probe signals were quantified using Agilent’s Feature Extraction software (version 10.7.3.1). The microarray data were deposited in the NCBI Gene Expression Omnibus (GEO) database (http://www.ncbi.nlm.nih.gov/geo/) under accession numbers GPL15243 (microarray including detailed annotation) and GSE37481.

Microarray data quality control and statistical analysis.

Statistical data analysis was performed on the processed Cy3 and Cy5 intensities, as provided by the Feature Extraction Software version 10.7. Further analysis was performed in the R programming environment, in conjunction with the packages developed within the Bioconductor project (http://www.bioconductor.org; [31]. Differential expression between the conditions was assessed via the moderated t-statistic, described in Smyth, G. [32] and implemented in the Limma package of Bioconductor. This moderated t-statistic applies an empirical Bayesian strategy to compute the gene-wise residual standard deviations and thereby increases the power of the test, especially beneficial for smaller data sets. To control the false discovery rate, multiple testing correction was performed [33] and a significant result was defined where any probe had a corrected p-value below 0.05 and an absolute fold change larger than 2. As all probes had 3 or 4 replicates on the array, only those probes were retained that were called differentially expressed in 2 out of 3 or 3 out of 4. Genes were considered as (significantly) differentially expressed if (1) their BH-corrected p-value fell below 0.05, (2) their fold changes were at least 2 and (3) if at least 75% of the different probes for a same gene met the first two criteria.

Microarray Data Functional Analysis

P. larvae genes were functionally annotated (already before qPCR validation) with Blast2GO [34], stand-alone BLAST [20] and KAAS (KEGG Automatic Annotation Server) [35] plus KEGG Mapper [36]. Putative transporter genes were also annotated with TransporterDB [37], putative transcription regulator genes with DBD [38] and RegPrecise [39] and putative proteolytic enzymes with MEROPS [40], [41]. COG functional categories were assigned with COGNITOR and stand-alone PSI-BLAST using the COG database [42]. After GO term annotation, an enrichment analysis (two-tailed Fisher’s exact test with default settings) within the Blast2GO environment was undertaken to compare different conditions.

Annotation and pathway pipeline.

The annotation and pathway pipeline employed here finds ORFs with Prodigal [43], annotates genes with blastp against KEGG [44], COG [42], Metacyc [45] and RefSeq [46] and makes a Pathway Genome Database (PGDB) with Pathway Tools [47].

Results

Honey Bee Larval Bodily Fluid Stimulates in vitro Growth of P. Larvae

The P. larvae strain BRL 230010 was genotypes via REP-PCR using ERIC primers and could be assigned to ERIC I. We studied the effect of larval bodily fluids on P. larvae gene expression by in vitro exposure of bacterial cultures to larval fluids. We could show treatment-induced bacterial growth shifts (Figure 1). Merely considering the two highest fluid concentrations (in the post hoc analysis; 1/10 and 1/25 dilutions), both seem to significantly (p<0.05) promote in vitro growth of P. larvae cultures, in comparison to control cultures. In addition (significantly) higher bacterial densities are observed as a result of larval fluid treatment (Figure 2). Thus, it is now reasonable to assume that these bacterial growth shifts reflect the phenotypic response to altered transcription and so the subsequent microarray experiment was be performed using 5% larval fluids.

[Figure omitted. See PDF.]

Figure 1. Effect of honey bee larval bodily fluid on Paenibacillus larvae growth.

Effect of different concentrations of honey bee larval bodily fluid on the in vitro growth of P. larvae bacterial cells, expressed as the optical density measured at a wavelength of 590 nm (OD590) in function of time (hours). Growth alterations were determined for six bodily fluid concentrations, expressed as the fold dilution in BHIT broth cultures: 10× dilution (♦), 25× dilution (▪), 50× dilution (▴), 100× dilution (⋄), 250× dilution (□), 500× dilution (▵), control (○). Each point in the graph displays the mean of three independent replicates, with the error bars being the standard deviations. Time 0 represents the time of spiking. Trend lines are calculated with DMfit. No OD measurements were performed during lag phase.

https://doi.org/10.1371/journal.pone.0089175.g001

[Figure omitted. See PDF.]

Figure 2. Effect of honey bee larval bodily fluid on Paenibacillus larvae bacterial density.

OD590 (white bars) and CFU/ml (black bars) for T1, C1, T4 and C4. T1 (n = 12): test sample collected one hour after spiking with 5% larval fluids. T4 (n = 12): test sample collected four hours after spiking with 5% larval fluids. C1 (n = 11): control sample collected one hour after spiking with BHIT-broth. C4 (n = 12): control sample collected four hours after spiking with BHIT-broth. Between brackets (n): number of independent replicates. Error bars: standard deviations.

https://doi.org/10.1371/journal.pone.0089175.g002

Reference Gene Selection for Normalization of RT-qPCR Data

The reliability of the data of the microarray experiment was checked with RT-qPCR. To correct for experimental error qPCR data require normalization against reference genes [48]. To this end, 9 reference genes were selected (Table S1 in File S1) and their expression in P. larvae was quantified for six different in vitro growth conditions: one, three or nine hours after spiking with larval bodily fluid (test) or BHIT broth (control). Arranging the genes from most to least stably expressed across all conditions, produced the following ranking (Figure S1A in File S1): rpoD (0.631)<gyrA (0.634)<cmk (0.656)<sucB (0.694)<eftu (0.848)<fum (0.944) <purH (1.03)<adk (1.116)<gapdh (1.195). Between brackets the geNorm M-value, assessing gene stability, is indicated. Another measure, the geNorm V-value, is useful for determining the optimal number of reference genes for data normalization (Figure S1B in File S1): V2/3 (0.214) – V3/4 (0.157) – V4/5 (0.200) – V5/6 (0.166) – V6/7 (0.155) – V7/8 (0.151) – V8/9 (0.141). Setting the threshold to 0.15, all tested genes except gapdh should be included in the calculation of the normalization factors [30].

Validation of Microarray Data by RT-qPCR

Microarray validation by RT-qPCR was performed with twenty randomly picked genes, belonging to different functional groups (Table S2 in File S1). Expression profiles obtained in the RT-qPCR experiment were similar in comparison to those of the microarray experiment, for all of the twenty genes in all of the four conditions (Figure S2 in File S1 and Figure S3 in File S1; Table S3 in File S1 and Table S4 in File S1).

General Overview of Differential Gene Expression as Revealed by Microarray Analysis

Depicting the microarray data as a heatmap analysis allowed us to correlate the corrected Cy3 and Cy5 dye intensities showed separate clusters for all eight samples (four replicates, each labeled with both Cy3 and Cy5) of C1, T1 and C4 respectively, and for six of the eight samples of T4 (not shown). In order not to corrupt true biological patterns, one of the four samples for T4, clustering together with C1, was omitted from further analysis. Table 1 shows the total number up- and down-regulated genes for each of the four comparisons (T4–C4, T1–C1, T4–T1, C4–C1). The three selection criteria were p-value <0.05, fold change ≥2 and at least 75% of the different, hybridized probes for the same gene shows differential expression for a particular comparison.

[Figure omitted. See PDF.]

Table 1. General overview of the microarray experiment, showing the total number of up- and down-regulated genes, for three selection criteria (columns) and four comparisons (rows).

https://doi.org/10.1371/journal.pone.0089175.t001

General Biological Patterns from Microarray Data by GO Term Assignment and Enrichment

At first GO terms were assigned to the P. larvae gene sequences. Using Blast2GO with Gene Ontology (GO) terms describing cellular component, biological process or molecular function, 2050 of the 3490 (unique) predicted genes within the P. larvae genome were assigned to at least one GO category. Subsequently GO term enrichment analysis identified (significantly) over- or underrepresented GO terms for the differentially expressed genes for each of the four comparisons (Table 2, Table 3, Table S5 in File S1 and Table S6 in File S1). The annotation and pathway pipeline annotated 10/147, 15/210, 9/136 and 45/373 of the up-regulated and 34/310, 19/248, 7/87 and 14/243 of the down-regulated genes for T1–C1, T4–C4, C4–C1 and T4–T1 (Table 4, Table 5, Table S7 in File S1 and Table S8 in File S1).

[Figure omitted. See PDF.]

Table 2. GO-enrichment analysis for comparison T1–C1, showing the most specific over- and under-represented biological process GO-terms for both up- and down-regulation.

https://doi.org/10.1371/journal.pone.0089175.t002

[Figure omitted. See PDF.]

Table 3. GO-enrichment analysis for comparison T4–C4, showing the most specific over- and under-represented biological process GO-terms for both up- and down-regulation.

https://doi.org/10.1371/journal.pone.0089175.t003

[Figure omitted. See PDF.]

Table 4. Annotation and pathway analysis for both up- and down-regulated genes for T1–C1.

https://doi.org/10.1371/journal.pone.0089175.t004

[Figure omitted. See PDF.]

Table 5. Annotation and pathway analysis for both up- and down-regulated genes for T4-C4.

https://doi.org/10.1371/journal.pone.0089175.t005

The differently expressed transporter and metabolic genes are represented in stacked percentage bar chart in Figure 3 and Figure 4 for T4-C4 and T1-C1 and Figure S4 in File S1 and Figure S5 in File S1 for T4-T1 and C4-C1 General trends in biological processes are also graphically depicted, according to the COG (Clusters of Orthologous Groups) functional classification (Figure 5 and Figure S6 in File S1).

[Figure omitted. See PDF.]

Figure 3. Differential expression of transporter encoding genes.

Stacked percentage bar chart, showing the numbers of up- and down-regulated (putative) transporter encoding genes for T4-C4 (left) and T1-C1 (right), respectively, relative to the total numbers of (putative) transporter encoding genes within the P. larvae genome. The latter are indicated between square brackets. Round brackets: GO term numbers. GO terms were assigned with Blast2GO. White bars: down-regulation for T4-C4. Dark grey bars: up-regulation for T4-C4. Light grey bars: down-regulation for T1-C1. Black bars: up-regulation for T1-C1. Arrow heads: arbitrary GO term hierarchy (◂>⊲><$>\raster(90%)="rg3"<$>).

https://doi.org/10.1371/journal.pone.0089175.g003

[Figure omitted. See PDF.]

Figure 4. Differential expression of metabolic genes.

Stacked percentage bar chart, showing the numbers of up- and down-regulated (putative) metabolic genes for T4-C4 (left) and T1-C1 (right), respectively, relative to the total numbers of (putative) metabolic genes within the P. larvae genome. C: carbohydrate metabolism. The latter are indicated between square brackets. Round brackets: KO numbers (KEGG pathways). KEGG pathways were assigned with KAAS. White bars: down-regulation for T4-C4. Dark grey bars: up-regulation for T4-C4. Light grey bars: down-regulation for T1-C1. Black bars: up-regulation for T1-C1. E: energy metabolism. L: lipid metabolism. N: nucleotide metabolism. A: amino acid metabolism. ∼: KO:00061; KO:00071; KO:00592; KO:01040. *: KO:00360; KO:00350; KO:00380. °: KO:00410; KO:00430; KO:00450; KO:00460; KO:00480. ?: KO:00740; KO:00770; KO:00780; KO:00670; KO:00860; KO:00130. ‘:KO:00900; KO:00903; KO:00281; KO:00523; KO:01053; KO:01055; KO:00940; KO:00311; KO:00521. ‘’: KO:00362; KO:00627; KO:00625; KO:00622; KO:00633; KO:00642; KO:00643; KO:00930; KO:00363; KO:00621; KO:00626.

https://doi.org/10.1371/journal.pone.0089175.g004

[Figure omitted. See PDF.]

Figure 5. Numbers of genes that are differentially expressed per subset.

Stacked percentage bar chart, showing the numbers of up- and down-regulated genes for T4-C4 (left) and T1-C1 (right), respectively, relative to the total numbers of genes within the P. larvae genome. The latter are indicated between square brackets. Round brackets: COG functional category label. C: cellular processes and signaling. I: information storage and processing. M: metabolism. P: poorly characterized. White bars: down-regulation for T4-C4. Dark grey bars: up-regulation for T4-C4. Light grey bars: down-regulation for T1-C1. Black bars: up-regulation for T1-C1.

https://doi.org/10.1371/journal.pone.0089175.g005

In control cultures (C1 and C4) P. larvae seems to invest more in oligopeptide and ion (sulfate ion; cation: Na+ for C1, K+ for C4) transport, amino acid (β-A, K; sulfur AAs for C1, BCAAs for C4), benzoate and limonene metabolism, while restricting nucleic acid metabolism. One hour after spiking (C1) P. larvae appears to rely more on carbohydrate (phosphotransferase system, inositol), carboxylic acid and ubiquinone metabolism, while limiting translation. The data also suggest that energy storage (electron transport, NADH to ubiquinone) and carbon fixation (photosynthesis, reductive citric acid cycle) are also upregulated. Four hours after spiking (C4) P. larvae tend to rely more on iron uptake (high-affinity iron transport and siderophore biosynthesis) and sporulation.

In test cultures, on the other hand, P. larvae seems to change its metabolism to nitrate assimilation and store energy through the electron transport chain. Four hours after spiking (T4) P. larvae appears to rely more on biotin, pyrimidine and carbon (TCA/glyoxylate cycle) metabolism. In addition, at late time points (T4-C4) P. larvae tend to restrict cellular macromolecular metabolic processes.

Discussion

The reasons behind the trends observed in the gene expression data presented here become clear if one considers the chemical composition of bodily fluids from honey bee worker larvae. These fluids are thought to contain mainly hemolymph and liquid gut contents. Since this material was collected from third, fourth and fifth larval instars, the gut is thought to hold (digested) honey/nectar, pollen and royal jelly [49]. Honey/nectar is particularly rich in sugars (especially fructose and glucose; [50], while royal jelly is a source of proteins [51]. The protein composition from honey bee worker larvae has been described in most detail [52]. Paralleling the diet, hemolymph of adult worker bees is rich in sugars [53], free amino acids [54] and proteins [55]. Fresh [56] and powdered [57] bee brood has been shown to be especially rich in amino acids/proteins.

The general down-regulation of carbohydrate and amino acid/peptide metabolism and transport in the presence of larval fluids at early time points (T1-C1) is unexpected, although genes involved in trehalose metabolism and transport are strongly upregulated: phosphotrehalase (treA; 12.93×↑) and trehalose-specific IIBC subunit of PTS system (treP; 16.85×↑). This could point towards carbon catabolite repression (CCR): a phenomenon by which usage of secondary carbon sources is reduced in the presence of the preferred one, which allows the fastest growth [58]. Thus trehalose, the principle blood sugar of many insects [59], is probably the preferred sugar of P. larvae. This hypothesis is supported by observations made for the Paenibacillus popilliae [60]. This pathogen of beetle larvae takes up trehalose through the PTS system and phosphorylates it to trehalose 6-phosphate (T6 P) in a PEP-dependent way. T6 P is subsequently cleaved to glucose and glucose-6-phosphate (G6 P). Down-regulation of genes for citric acid cycle enzymes (T1-C1; citB: 1.32×↓ (ns); citC: 4.53×↓; citH: 3.64×↓; citZ: 8.39×↓; [61]) and virulence factors [62] - exemplified by toxin Etx/Mtx (T1-C1; >8.50×↓) - provide additional support to the CCR hypothesis. Genes ccpA and codY, however, are not differentially expressed. Their products are global regulators of carbon (glucose) metabolism by sensing FBP and G6P (for CcpA) and GTP and BCAA (for CodY) in Bacillus subtilis [61]. Additionally, fatty acid catabolism in B. subtilis is repressed through CCR [63]. In contrast, this study reveals a strong up-regulation for genes involved in fatty acid degradation (eftA, B: >24.50×↑; fadA, E, F, N, R: >20↑; lcfA: 14.51×↑).

The down-regulation of many oligopeptide permeases (opp) might reflect the high free amino acid content of insect hemolymph [64]. The investment of P. larvae in peptide transport, subsequent catabolism and amino acid synthesis will be minimal as amino acids are copiously provided. Sulphate transport is also repressed when cysteine is present [65] and because of the high free amino acid availability, the down-regulation of sulphate transporters (T1-C1; cysA, P, U, W; >40×↓) can be explained.

The up-regulated import of especially mannose (levD, E, F, G; >5×↑/>11×↑) at later time points (T4-C4/T4-T1) indicate that the consumption of trehalose proceeds and P. larvae seems to switch to alternative carbon sources at higher cell densities. Although high densities are reached at T4 (OD590 ≈ 1.14), no starvation seems to occur (e.g. ctsA: >2.3×↓) and as a consequence sporulation is not initiated. Sporulation starts at late time points (C4-C1), which is suggested by the up-regulation of 19 from the 52 spore/sporulation related genes in the control cultures. Addition of larval bodily fluids delays sporulation initiation, which is illustrated by the down regulation of 35/52 sporulation related genes in the treated cultures; none of these genes showed up-regulation. Moreover, sporulation is mostly regulated by environmental signals through chemotaxis [66]. The gene expression patterns of putative signaling proteins (e.g. MCPs, HPKs, RRs) are similar to those of the sporulation related genes. Nutrient limitation and competition often leads to the production of antibiotics [67], which is reflected in the up-regulation of (putative) synthetases for antimicrobial biosurfactans, plipastatin (five genes) and surfactin (four genes) at later time points (T4-T1; C4-C1) however addition of larval bodily fluids seems to diminish the production and transport of antibiotics (T4-C4; T1-C1). Other antibiotic-related genes - amongst which penicillin acylase [68], two peptides AS-48 [69], four lantibiotic/lanthionine synthesis proteins [70], neotrehalosadiamine [71] and circularin A/uberolysin [72] - together show a similar pattern (with some exceptions). The same holds more or less for polyketide biosynthesis genes (e.g. rfbB, C, D of the rhamnose pathway; [73]). Closticin 574 tends to be up-regulated by larval fluids at late time points (T4-C4, T4-T1) [74].

Differential expression of (putative) cofactor and vitamin metabolic genes was observed for different pathways but one of the most up-regulated pathways in the presence of larval fluids at late time points is the one responsible for thiamine biosynthesis (e.g. apbE: 85.80×↑ for T4-C1 and 54.08×↑ for T4-T1; [75]). The same pathway is down-regulated in control cultures at late time points (C4-C1). This indicates that thiamine seems exhausted and biosynthesis of this vitamin/cofactor seems necessary. Another remarkable fact is the quite large-scale phage-related gene expression at late time points in the presence of larval fluids. If bacteriophage particles are produced, this might reflect spontaneous induction of the prophages’ lytic cycle through signal molecules of bacterial quorum sensing [76].

At late time points, oxygen appears to become limited. As a consequence P. larvae seems to switch from aerobic to anaerobic respiration/fermentation. Many bacteria follow a mixed acid fermentation route for glucose metabolism with end products ethanol, succinate, lactate, acetate, formate, and carbon dioxide [77]. Typical indicators of this process arise from the activity of its key enzyme, pyruvate formate lyase (Pfl), which leads to massive excretion of formate and acetate as fermentative by-products. Down-regulation of pyruvate dehydrogenase (pdhA: 5.39×↓ and pdhB: 5.01×↓) and up-regulation of pyruvate formate lyase (pflA: 21.61×↑, pflB: 9.56×↑) confirm this hypothesis. Moreover P. larvae seems to rely mostly on cytochrome aa3-600 quinol oxidase (qoxA, B, C, D: >9×↑; [78]) at T1-C1, while it switches to cytochrome bd ubiquinol oxidase (cydA: 76.47×↑, cydB: 37.08×↑) at T4-C4. The latter enzyme is expressed under oxygen-limiting conditions in E. coli [79]. Furthermore nitrogen metabolism seems boosted at T4-C4: formate/nitrite transport (yrhG: >3×↑), assimilatory nitrite reductase (nasD, E: >10×↑) and anaerobic regulator (fnr: 11.89×up) are all up-regulated. NasDE nitrite reductase functions as both an assimilatory and a dissimilatory enzyme which is produced under oxygen limiting and nitrogen limiting conditions [80]. The bacterial nitrogen cycle is complicated and still a matter of debate as different pathways are involved such as the dissimilative nitrate reduction to ammonia and the denitrification pathway with reduction of nitrate to N2. On the other hand, homologues to bifunctional alcohol/acetaldehyde dehydrogenase (adhE) are highly up-regulated at T4-C4 (108.91×↑, 97.52×↑ and 60.68× respectively↑). This could point towards sugar fermentation, converting acetyl-coA to acetaldehyde and then to ethanol. It has been suggested for E. coli that transcription of adhE is repressed through CCR under aerobic conditions (T1-C1; ±2×↑) and that expression of adhE is much higher in anaerobically grown cells [81].

Some other remarkable differences in gene expression as a result of the addition of larval bodily fluids are three differentially expressed (putative) toxin encoding genes. Two homologues of ε-toxins (ETX) or mosquitocidal toxins (MTX) [82] are down-regulated for T1-C1 (8.53×↓, 8.75×↓) and C4-C1 (2.56×↓, 2.60×↓), and up-regulated for T4-T1 (4.47×↑, 4.47×↑). The ι-toxin Ib component [83], on the contrary, is up-regulated for both T1-C1 (2.62×↑) and C4-C1 (2.06×↑). Next to toxins, proteases/peptidases have been proposed as putative P. larvae virulence factors [20], [62], [84], [85]. However, no general pattern emerges from the expression data which can be explained by their divers biological roles. Recently S-layer proteins have been suggested as potential P. larvae virulence determinants [86], [87]. Our study revealed up-regulation (±2.5×↑) of three such genes at late time points (T4–T1). At late time points the presence of larval bodily fluids favors motility which is reflected by the up-regulation of 7/8 genes. The eight differentially expressed (putative) flagellar genes (MC) encode (1) filament capping-protein (fliD), (2) filament (fliC), (3) filament-hook junction (flgK, flgL), (4) basal body stator (motA, motB), (5) a flagellar export chaperone for FliC (fliS) and (6) others (yvyF) [88].

Conclusions

We can conclude that the honey bee larval bodily fluid stimulates in vitro growth of P. larvae which is reflected by the large amount of transcriptional changes. Early responses are characterized by a general down-regulation of transporter genes and genes involved in the amino acid and carbohydrate metabolism. At later time points the sporulation genes are down-regulated while phage-related genes are up-regulated. The importance of the changed expression of phage-related genes will be a subject for further research.

Supporting Information

[Figure omitted. See PDF.]

File S1.

Figure S1. Average expression stability (A) and determination of the optimal number (B) of reference targets with geNormPLUS. Figure S2. Validation of microarray data with RTQ-PCR. (A) Log2-transformed expression ratio of T1 compared to C1. (B) Log2-transformed expression ratio of T4 compared to C4. White bars: RTQ-PCR experiment. Black bars: microarray experiment. T1: test sample collected one hour after spiking with larval fluids. T4: test sample collected four hours after spiking with larval fluids. C1: control sample collected one hour after spiking with BHIT-broth. C4: control sample collected four hours after spiking with BHIT-broth. +: differential expression (significant). -: equal expression (non-significant). Figure S3. Validation of microarray data with RTQ-PCR. (A) Log2-transformed expression ratio of C4 compared to C1. (B) Log2-transformed expression ratio of T4 compared to T1. White bars: RTQ-PCR experiment. Black bars: microarray experiment. T1: test sample collected one hour after spiking with hemolymph. T4: test sample collected four hours after spiking with hemolymph. C1: control sample collected one hour after spiking with BHIT-broth. C4: control sample collected four hours after spiking with BHIT-broth. +: differential expression (significant). -: equal expression (non-significant). Figure S4. Stacked percentage bar chart, showing the numbers of up- and down-regulated (putative) transporter encoding genes for T4–T1 (left) and C4–C1 (right), respectively, relative to the total numbers of (putative) transporter encoding genes within the P. larvae genome. The latter are indicated between square brackets. Round brackets: GO term numbers. GO terms were assigned with Blast2GO. White bars: down-regulation for T4–T1. Dark grey bars: up-regulation for T4–T1. Light grey bars: down-regulation for C4–C1. Black bars: up-regulation for C4–C1. Arrow heads: arbitrary GO term hierarchy (;◂>⊲><$>\raster="rg3"<$>). Figure S5. Stacked percentage bar chart, showing the numbers of up- and down-regulated (putative) metabolic genes for T4–T1 (left) and C4–C1 (right), respectively, relative to the total numbers of (putative) metabolic genes within the P. larvae genome. C: carbohydrate metabolism. The latter are indicated between square brackets. Round brackets: KO numbers (KEGG pathways). KEGG pathways were assigned with KAAS. White bars: down-regulation for T4–T1. Dark grey bars: up-regulation for T4–T1. Light grey bars: down-regulation for C4–C1. Black bars: up-regulation for C4–C1. ∼: KO:00061; KO:00071; KO:00592. *: KO:00350; KO:00360; KO:00380; KO:00400. °: KO:00410; KO:00430; KO:00450; KO:00480. ?: KO:00740; KO:00760; KO:00770; KO:00780; KO:00790; KO:00670; KO:00860; KO:00130. ‘:KO:00900; KO:00903; KO:00281; KO:00523; KO:01053; KO:01055; KO:00960; KO:00232; KO:00521; KO:00401. ‘’: KO:00362; KO:00627; KO:00625; KO:00622; KO:00633; KO:00642; KO:00643; KO:00930; KO:00621; KO:00626; KO:00983. Figure S6. Stacked percentage bar chart, showing the numbers of up- and down-regulated genes for T4–T1 (left) and C4–C1 (right), respectively, relative to the total numbers of genes within the P. larvae genome. The latter are indicated between square brackets. Round brackets: COG functional category label. C: cellular processes and signaling. I: information storage and processing. M: metabolism. P: poorly characterized. White bars: down-regulation for T4–T1. Dark grey bars: up-regulation for T4–T1. Light grey bars: down-regulation for C4–C1. Black bars: up-regulation for C4–C1. Table S1. Reference genes used to normalize the results of the qRT-PCR experiment for microarray data validation. Table S2. Randomly selected genes used in the RTQ-PCR experiment to validate the microarray results. Table S3. Validation of microarray data by qRT-PCR. Table S4. Validation of microarray data by qRT-PCR. Table S5. GO-enrichment analysis for comparison C4–C1, showing the most specific over- and under-represented biological process GO-terms for both up- and down-regulation. Table S6. GO-enrichment analysis for comparison T4–T1, showing the most specific over- and under-represented biological process GO-terms for both up- and down-regulation. Table S7. Annotation and pathway analysis for both up- and down-regulated genes for C4–C1. Table S8. Annotation and pathway analysis for both up- and down-regulated genes for T4–T1.

https://doi.org/10.1371/journal.pone.0089175.s001

(PDF)

Acknowledgments

Dieter De Koker was a talented UGent PhD-student and the principal investigator of this study. He died suddenly and unexpectedly soon after he wrote the draft version of this paper. The co-authors wanted to honor him by finishing the paper.

Author Contributions

Conceived and designed the experiments: LDS DDK PDV DCDG. Performed the experiments: LDS DDK. Analyzed the data: LDS DDK AKH LJF. Wrote the paper: LDS DDK LJF PDV DCDG.

Citation: De Smet L, De Koker D, Hawley AK, Foster LJ, De Vos P, de Graaf DC (2014) Effect of Bodily Fluids from Honey Bee (Apis mellifera) Larvae on Growth and Genome-Wide Transcriptional Response of the Causal Agent of American Foulbrood Disease (Paenibacillus larvae). PLoS ONE 9(2): e89175. https://doi.org/10.1371/journal.pone.0089175

References

1. Ellis JD, Munn PA (2005) The worldwide health status of honey bees. Bee World 86: 88–101.

2. Crailsheim K, Riessberger-Galle U (2001) Honey bee age-dependent resistance against American foulbrood. Apidologie 32: 91–103.

3. Brødsgaard CJ, Ritter W, Hansen H (1998) Response of in vitro reared honey bee larvae to various doses of Paenibacillus larvae larvae spores. Apidologie 29: 569–578.

4. Yue D, Nordhoff M, Wieler LH, Genersch E (2008) Fluorescence in situ hybridization (FISH) analysis of the interactions between honey bee larvae and Paenibacillus larvae, the causative agent of American foulbrood of honey bees (Apis mellifera). Environ Microbiol 10: 1612–1620.

5. Sturtevant AP (1932) Relation of commercial honey to the spread of American foulbrood. J Agricul Res 45: 257–285.

6. Hasemann L (1961) How long can spores of American foulbrood live? Am Bee J 101: 298–299.

7. Calesnick EJ, White JW (1952) Thermal resistance of Bacillus larvae spores in honey. J Bacteriol 64: 9–15.

8. de Graaf DC, Vandekerchove D, Dobbelaere W, Peeters JE, Jacobs FJ (2001) Influence of the proximity of American foulbrood cases and apicultural management on the prevalence of Paenibacillus larvae spores in Belgian honey. Apidologie 32: 587–599.

9. Miyagi T, Peng CYS, Chuang RY, Mussen EC, Spivak MS, et al. (2000) Verification of oxytetracycline-resistant American foulbrood pathogen Paenibacillus larvae in the United States. J Invertebr Pathol 75: 95–96.

10. Reynaldi FJ, Lacunza J, Alippi AM, Rule R (2010) Binding of tylosin, tilmicosin and oxytetracycline to proteins from honey bees, larvae and beehive products. Rev Argent Microbiol 42: 279–283.

11. Buczek K (2011) Range of susceptibility of Paenibacillus larvae to antibacterial compounds. Med Weter 67: 87–90.

12. Forsgren E, Olofsson TC, Vasquez A, Fries I (2010) Novel lactic acid bacteria inhibiting Paenibacillus larvae in honey bee larvae. Apidologie 41: 99–108.

13. Fuselli SR, de la Rosa SBG, Eguaras MJ, Fritz R (2008) Susceptibility of the honey bee bacterial pathogen Paenibacillus larvae to essential oils distilled from exotic and indigenous Argentinean plants. J Essent Oil Res 20: 464–470.

14. Bastos EMAF, Simone M, Jorge DM, Soares AEE, Spivak M (2008) In vitro study of the antimicrobial activity of Brazilian propolis against Paenibacillus larvae. J Invertebr Pathol 97: 273–281.

15. Feldlaufer MF, Knox DA, Lusby WR, Shimanuki H (1993) Antimicrobial activity of fatty acids against Bacillus larvae, the causative agent of American foulbrood disease. Apidologie 24: 95–99.

16. De Guzman ZM, Cervancia CR, Dimasuay KGB, Tolentino MM, Abrera GB, et al. (2011) Radiation inactivation of Paenibacillus larvae and sterilization of American Foul Brood (AFB) infected hives using Co-60 gamma rays. Appl Radiat Isotopes 69: 1374–1379.

17. Pernal SF, Albright RL, Melathopoulos AP (2008) Evaluation of the shaking technique for the economic management of American foulbrood disease of honey bees (Hymenoptera: Apidae). J Econ Entomol 101: 1095–1104.

18. Genersch E (2010) American Foulbrood in honey bees and its causative agent, Paenibacillus larvae. J Invertebr Pathol 103: S10–S19.

19. Qin X, Evans JD, Aronstein KA, Murray KD, Weinstock GM (2006) Genome sequences of the honey bee pathogens Paenibacillus larvae and Ascosphaera apis. Insect Mol Biol 15(5): 715–718.

20. Chan QWT, Cornman RS, Birol I, Liao NY, Chan SK, et al. (2011) Updated genome assembly and annotation of Paenibacillus larvae, the agent of American foulbrood disease of honey bees. BMC Genomics 12: 450.

21. Jandu N, Ho NKL, Donato KA, Karmali MA, Mascarenhas M, et al. (2009) Enterohemorrhagic Escherichia coli O157: H7 gene expression profiling in response to growth in the presence of host epithelia. Plos One 4(3): e4889.

22. Samant S, Lee H, Ghassemi M, Chen J, Cook JL, et al. (2008) Nucleotide biosynthesis is critical for growth of bacteria in human blood. PLoS Pathog 4(2): e37.

23. Mattinen L, Somervuo P, Nykyri J, Nissinen R, Kouvonen P, et al. (2008) Microarray profiling of host-extract-induced genes and characterization of the type VI secretion cluster in the potato pathogen Pectobacterium atrosepticum. Microbiol-Sgm 154: 2387–2396.

24. Genersch E, Forsgren E, Pentikäinen J, Ashiralieva A, Rauch S, et al. (2006) Reclassification of Paenibacillus larvae subsp. pulvifaciens and Paenibacillus larvae subsp. larvae as Paenibacillus larvae without subspecies differentiation. Int J Syst Evol Microbiol 56: 501–511.

25. de Graaf DC, Alippi AM, Antúnez K, Aronstein KA, Budge G, et al.. (2013) Standard methods for American foulbrood research. In Dietemann V., Ellis J.D., Neumann P. (Eds) The COLOSS BEEBOOK, Volume II: standard methods for Apis mellifera pest and pathogen research. J Apic Res 52(1): http://dx.doi.org/10.3896/IBRA.1.52.1.11.

26. Baranyi J, Roberts TA (1994) A dynamic approach to predicting bacterial growth in food. Int J Food Microbiol 23: 277–294.

27. Untergasser A, Nijveen H, Rao X, Bisseling T, Geurts R, et al. (2007) Primer3Plus, an enhanced web interface to Primer3. Nucleic Acids Res 35: W71–W74.

28. Zuker M (2003) Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res 31: 3406–3415.

29. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, et al. (2002) Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol 3(7): 0034.1.

30. Hellemans J, Mortier G, De Paepe A, Speleman F, Vandesompele J (2007) qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biol 8(2): R19.

31. Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, et al. (2004) Bioconductor: open software development for computational biology and bioinformatics. Genome Biol 5(10): R80.

32. Smyth GK (2004) Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol 3: 1.

33. Benjamini Y, Hochberg Y (1995) Controlling the False Discovery Rate - a Practical and Powerful Approach to Multiple Testing. J Roy Stat Soc B Met 57: 289–300.

34. Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, et al. (2005) Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21: 3674–3676.

35. Moriya Y, Itoh M, Okuda S, Yoshizawa AC, Kanehisa M (2007) KAAS: an automatic genome annotation and pathway reconstruction server. Nucleic Acids Res 35: W182–W185.

36. Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M (2012) KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 40: D109–D114.

37. Ren QH, Kang KH, Paulsen IT (2004) TransportDB: a relational database of cellular membrane transport systems. Nucleic Acids Res 32: D284–D288.

38. Kummerfeld SK, Teichmann SA (2006) DBD: a transcription factor prediction database. Nucleic Acids Res 34: D74–D81.

39. Novichkov PS, Laikova ON, Novichkova ES, Gelfand MS, Arkin AP, et al. (2010) RegPrecise: a database of curated genomic inferences of transcriptional regulatory interactions in prokaryotes. Nucleic Acids Res 38: D111–D118.

40. Rawlings ND, O’Brien E, Barrett AJ (2002) MEROPS: the protease database. Nucleic Acids Res 30: 343–346.

41. Rawlings ND, Tolle DP, Barrett AJ (2004) MEROPS: the peptidase database. Nucleic Acids Res 32: D160–D164.

42. Tatusov RL, Galperin MY, Natale DA, Koonin EV (2000) The COG database: a tool for genome-scale analysis of protein functions and evolution. Nucleic Acids Res 28: 33–36.

43. Hyatt D, Chen G, LoCascio PF, Land ML, Larimer FW, et al. (2010) Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics 11: 119.

44. Kanehisa M, Goto S (2000) KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 28: 27–30.

45. Krieger CJ, Zhang P, Mueller LA, Wang A, Paley S, et al. (2004) MetaCyc: a multiorganism database of metabolic pathways and enzymes. Nucleic Acids Res. 32: D438–D442.

46. Pruitt KD, Tatusova T, Maglott DR (2006) NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 35: D61–D65.

47. Karp PD, Paley SM, Krummenacker M, Latendresse M, Dale JM, et al. (2009) Pathway Tools version 13.0: integrated software for pathway/genome informatics and systems biology. Brief. Bioinform. 11: 40–79.

48. Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, et al. (2009) The MIQE Guidelines: Minimum Information for Publication of Quantitative Real-Time PCR Experiments. Clin Chem 55: 611–622.

49. Winston ML (1987) The biology of the honey bee. Harvard University Press, Cambridge, MA.

50. Bogdanov S, Jurendic T, Sieber R, Gallmann P (2008) Honey for nutrition and health: a review. J Am Coll Nutr. 27: 677–689.

51. Ramadan MF, Al-Ghamdi A (2012) Bioactive compounds and health-promoting properties of royal jelly: A review. J Funct Foods 4: 39–52.

52. Chan QW, Foster LJ (2008) Changes in protein expression during honey bee larval development. Genome Biol. 9: R156.

53. Blatt J, Roces F (2001) Haemolymph sugar levels in foraging honeybees (Apis mellifera carnica): dependence on metabolic rate and in vivo measurement of maximal rates of trehalose synthesis. J Exp Biol. 204: 2709–2716.

54. Crailsheim K, Leonhard B (1997) Amino acids in honeybee worker haemolymph. Amino Acids 13: 141–153.

55. Chan QW, Howes CG, Foster LJ (2006) Quantitative comparison of caste differences in honeybee hemolymph. Mol Cell Proteomics. 5: 2252–2262.

56. Finke MD (2005) Nutrient Composition of Bee Brood and its Potential as Human Food. Ecol Food Nutr 44: 257–270.

57. Narumi S (2004) Honeybee brood as a nutritional food. Honeybee Science 25 (3): 119–124.

58. Görke B, Stülke J (2008) Carbon catabolite repression in bacteria: many ways to make the most out of nutrients. Nat Rev Microbiol 6(8): 613–624.

59. Wyatt GR (1961) the Biochemistry of insect hemolymph. Ann Rev Entomol 6: 75–102.

60. Bhumiratana A, Anderson RL, Costilow RN (1974) Trehalose metabolism by Bacillus popilliae. J Bacteriol. 119: 484–493.

61. Sonenshein AL (2007) Control of key metabolic intersections in Bacillus subtilis. Nat Rev Microbiol. 5: 917–927.

62. Antúnez K, Arredondo D, Anido M, Zunino P (2011) Metalloprotease production by Paenibacillus larvae during the infection of honey bee larvae. Microbiol-Sgm 157: 1474–1480.

63. Tojo S, Satomura T, Matsuoka H, Hirooka K, Fujita Y (2011) Catabolite repression of the Bacillus subtilis FadR regulon, which is involved in fatty acid catabolism. J Bacteriol. 193: 2388–2395.

64. Gilliam M, McCaughey WF (1972) Total amino acids in developing worker honey bees (Apis mellifera L.). Expermentia 28: 143.

65. Aguilar-Barajas E, Díaz-Pérez C, Ramírez-Díaz MI, Riveros-Rosas H, Cervantes C (2011) Bacterial transport of sulfate, molybdate, and related oxyanions. Biometals 24: 687–707.

66. Wadhams GH, Armitage JP (2004) Making sense of it all: Bacterial chemotaxis. Nat Rev Mol Cell Bio 5: 1024–1037.

67. Stein T (2005) Bacillus subtilis antibiotics: structures, syntheses and specific functions. Mol Microbiol. 56: 845–857.

68. Chandel AK, Rao LV, Narasu ML, Singh OV (2008) The realm of penicillin G acylase in beta-lactam antibiotics. Enzyme Microb Tech 42: 199–207.

69. Maqueda M, Galvez A, Bueno MM, Sanchez-Barrena MJ, Gonzalez C, et al. (2004) Peptide AS-48: Prototype of a new class of cyclic bacteriocins. Curr Protein Pept Sc 5: 399–416.

70. Willey JM, van der Donk WA (2007) Lantibiotics: Peptides of diverse structure and function. Annu Rev Microbiol 61: 477–501.

71. Inaoka T, Ochi K (2011) Activation of dormant secondary metabolism neotrehalosadiamine synthesis by an RNA polymerase mutation in Bacillus subtilis. Biosci Biotech Bioch 75: 618–623.

72. Wirawan RE, Swanson KM, Kleffmann T, Jack RW, Tagg JR (2007) Uberolysin: a novel cyclic bacteriocin produced by Streptococcus uberis. Microbiol-Sgm 153: 1619–1630.

73. Giraud MF, Naismith JH (2000) The rhamnose pathway. Curr Opin Struc Biol 10: 687–696.

74. Kemperman R, Kuipers A, Karsens H, Nauta A, Kuipers O, et al. (2003) Identification and characterization of two novel clostridial bacteriocins, circularin A and closticin 574. Appl Environ Microb 69: 1589–1597.

75. Beck BJ, Downs DM (1998) The apbE gene encodes a lipoprotein involved in thiamine synthesis in Salmonella typhimurium. J Bacteriol. 180: 885–891.

76. Ghosh D, Roy K, Williamson KE, Srinivasiah S, Wommack KE, et al. (2009) Acyl-homoserine lactones can induce virus production in lysogenic bacteria: an alternative paradigm for prophage induction. Appl Environ Microbiol. 75: 7142–7152.

77. Bock A, Sawers G (1996) Chapter 18: Fermentation, in Escherichia coli and Salmonella, Cellular and Molecular Biology, Second Edition Neidhardt FC, Curtiss III R, Ingraham JL, Lin ECC, Low Jr KB, Magasanik B, Reznikoff WS, Riley M, Schaechter M, Umbarger HE American Society for Microbiology, Washington, D.C.,

78. Yi SM, Narasimhulu KV, Samoilova RI, Gennis RB, Dikanov SA (2010) Characterization of the semiquinone radical stabilized by the cytochrome aa3–600 menaquinol oxidase of Bacillus subtilis. J Biol Chem 285(24): 18241–18251.

79. Borisov VB, Gennis RB, Hemp J, Verkhovsky MI (2011) The cytochrome bd respiratory oxygen reductases. Biochim Biophys Acta 1807: 1398–1413.

80. Nakano MM, Zuber P (1998) Anaerobic growth of a ‘strict aerobe’ Bacillus subtilis. Ann Rev Microbiol 52: 165–190.

81. Mikulskis A, Aristarkhov A, Lin EC (1997) Regulation of expression of the ethanol dehydrogenase gene (adhE) in Escherichia coli by catabolite repressor activator protein Cra. J Bacteriol 179(22): 7129–7134.

82. Bokori-Brown M, Savva CG, da Costa SPF, Naylor CE, Basak AK, et al. (2011) Molecular basis of toxicity of Clostridium perfringens epsilon toxin. FEBS J 278: 4589–4601.

83. Sakurai J, Nagahama M, Oda M, Tsuge H, Kobayashi K (2009) Clostridium perfringens Iota-Toxin: Structure and Function. Toxins 1: 208–228.

84. Antúnez K, Anido M, Schlapp G, Evans JD, Zunino P (2009) Characterization of secreted proteases of Paenibacillus larvae, potential virulence factors involved in honey bee larval infection. J Invertebr Pathol 102: 129–132.

85. Antúnez K, Anido M, Schlapp G, Evans JD, Zunino P (2010) Secreted and immunogenic proteins produced by the honey bee bacterial pathogen, Paenibacillus larvae. Vet Microbiol 141: 385–389.

86. Fünfhaus A, Genersch E (2012) Proteome analysis of Paenibacillus larvae reveals the existence of a putative S-layer protein. Environ Microbiol Rep 4(2): 194–202.

87. Poppinga L, Janesch B, Fünfhaus A, Sekot G, Garcia-Gonzalez E, et al. (2012) Identification and functional analysis of the S-layer protein SplA of Paenibacillus larvae, the causative agent of American Foulbrood of honey bees. PLoS Pathog 8: e1002716.

88. Chevance FFV, Hughes KT (2008) Coordinating assembly of a bacterial macromolecular machine. Nat Rev Microbiol 6: 455–465.