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)

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.


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][22][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.

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 mg per ml phosphate buffered saline (PBS, 20 mM KH 2 PO 4 , 60 mM Na 2 HPO 4 and 145 mM NaCl, pH 7.2)], an inhibitor of phenoloxidase, was added immediately after squeezing the larvae to a final concentration of 10 mg/ml to prevent melanization. The homogenate was centrifuged twice for 30 min at 75,6006g and 4uC with an Avanti J30-I centrifuge (Beckman Coulter, Inc.). The clarified supernatant was filter sterilized with a 0.2 mm 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 37uC 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 37uC. 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 OD 590 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 OD 595 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 (116t1), eleven others for three hours (116t3) and the remaining eleven for nine hours (116t9). The same incubation regime was dictated to the control cultures  (116c1, 116c3, 116c9), thereby creating six different conditions (t1, t3, t9, c1, c3, c9). Bacterial cells were collected by centrifuging for 5 min at 87206g and 4uC. The bacterial pellets were immediately resuspended in 650 ml RNAlater Solution (Ambion). Aliquots (50 ml each) of these suspensions were incubated on ice for 30 min, centrifuged for 5 min at 67086g and 4uC and stored at 220uC until RNA preparation. RNA preparation. Prior to RNA extraction thawed aliquots were centrifuged for 2 min at 67086g and 4uC. Bacterial pellets were resuspended by vortexing in 100 ml 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 Grampositive 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 ml elution buffer and stored at 220uC. cDNA synthesis. RNA (5 mg) 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 (60uC) and ionic concentrations ([Na + ] = 50 mM, [Mg 2+ ] = 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 ml reaction consisted of 7.5 ml SYBR master mix, 0.2 mM forward and 0.2 mM reverse primers (Integrated DNA Technologies) and 1 ml cDNA template using the CFX96 Real-Time PCR Detection System (Bio-Rad). The PCR program comprises a UDG digestion step of 2 min at 50uC, an activation step of 2 min at 95uC and 40 cycles of a combined denaturation (20 sec at 95uC) and annealing (40 sec at 60uC) step. At the end of this program a melt curve is generated by measuring fluorescence after each temperature increase of 0.5uC for 5 sec over a range from 65uC to 95uC. Primer efficiencies, R 2 values and melt curves were calculated with CFX Manager Software (Bio-Rad). Reference gene stability was analyzed with the geNorm PLUS algorithm within the qBase PLUS 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 qBase PLUS [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 OD 595 of 0.3. Four of the test and control cultures were incubated for 1 h (46 T1, and 46 C1) and the others for 4 h (46 T4, 46 C4). The bacterial cells were collected as described for reference gene selection and stored at 280uC until RNA preparation.
An accessory evaluation of P. larvae growth promotion in response to larval fluids was executed by comparing OD 590 measurements on all bacterial cultures subjected to microarray analysis (16 samples: 46C1, 46T1, 46C4, 46T4). 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 ml 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 oncolumn DNase I treatment, according to the manufacturer's instructions. The RNA pellet was dissolved in 30 ml RNase-free water and stored at 280uC until use. Microarray experimental procedures. A custom 8615 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 mg 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) monoreactive 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 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.

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.

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

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.

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).
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).
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 (b-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 (highaffinity 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.936q) and trehalose-specific IIBC subunit of PTS system (treP; 16.856q). 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.326Q (ns); citC: 4.536Q; citH: 3.646Q; citZ: 8.396Q; [61]) and virulence factors [62] exemplified by toxin Etx/Mtx (T1-C1; .8.506Q) -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.506q; fadA, E, F, N, R: .20q; lcfA: 14.516q).
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; .406Q) can be explained.
The up-regulated import of especially mannose (levD, E, F, G; .56q/.116q) 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 (OD 590 < 1.14), no starvation seems to occur (e.g. ctsA: .2.36Q) 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.806q for T4-C1 and 54.086q 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 spontane-ous 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.396Q and pdhB: 5.016Q) and up-regulation of pyruvate formate lyase (pflA: 21.616q, pflB: 9.566q) confirm this hypothesis. Moreover P. larvae seems to rely mostly on cytochrome aa3-600 quinol oxidase (qoxA, B, C, D: .96q; [78]) at T1-C1, while it switches to cytochrome bd ubiquinol oxidase (cydA: 76.476q, cydB: 37.086q) 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: .36q), assimilatory nitrite reductase (nasD, E: .106q) and anaerobic regulator (fnr: 11.896up) 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 N 2 . On the other hand, homologues to bifunctional alcohol/acetaldehyde dehydrogenase (adhE) are highly up-regulated at T4-C4 (108.916q, 97.526q and 60.686 respectivelyq). 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; 626q) and that expression of adhE is much higher in anaerobically grown cells [81].

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 downregulated 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
File S1 Figure S1. Average expression stability (A) and determination of the optimal number (B) of reference targets with geNorm PLUS . Figure S2. Validation of microarray data with RTQ-PCR. (A) Log 2 -transformed expression ratio of T1 compared to C1. (B) Log 2 -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) Log 2transformed expression ratio of C4 compared to C1. (B) Log 2transformed 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 BHITbroth. +: 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 (;b.o. ). 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.  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. GOenrichment 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 overand under-represented biological process GO-terms for both upand 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. (PDF)