Shigatoxin encoding Bacteriophage ϕ24B modulates bacterial metabolism to raise antimicrobial tolerance

How temperate bacteriophages play a role in microbial infection and disease progression is not fully understood. They do this in part by carrying genes that promote positive evolutionary selection for the lysogen. Using Biolog phenotype microarrays and comparative metabolite profiling we demonstrate the impact of the well-characterised Shiga toxin-prophage ϕ24B on its Escherichia coli host MC1061. As a lysogen, the prophage alters the bacterial physiology by increasing the rates of respiration and cell proliferation. This is the first reported study detailing phage-mediated control of the E. coli biotin and fatty acid synthesis that is rate limiting to cell growth. Through ϕ24B conversion the lysogen also gains increased antimicrobial tolerance to chloroxylenol and 8-hydroxyquinoline. Distinct metabolite profiles discriminate between MC1061 and the ϕ24B lysogen in standard culture, and when treated with 2 antimicrobials. This is also the first reported use of metabolite profiling to characterise the physiological impact of lysogeny under antimicrobial pressure. We propose that temperate phages do not need to carry antimicrobial resistance genes to play a significant role in tolerance to antimicrobials.

accessory to the core biology of the phage may offer a selective advantage to the host. Mis-excision, mis-packaging of phage DNA 14 , and recombination 15 play a large role in phage genome variation. This may leave phage DNA regions, remnant or cryptic prophages, that positively impact on the selection and survival of both the phage and the bacterium 16 . This is further supported by the common occurrence of prophage regions, usually multiple, in the chromosomes of many bacterial pathogens 17 .
There are a number of features of the Shigatoxigenic phage vB_EcoPφ 24 B or φ24 B that are particularly relevant to the success and persistence of Stx prophages in E. coli. In contrast to the lambda infection model, φ 24 B can multiply infect a bacterial host [18][19][20] . Stx-phages have been isolated from a wide variety of environments where E. coli is present, and this has undoubtedly been promoted by the use of an essential outer membrane protein BamA as the adsorption site 21 . This interaction is conserved in Stx-phages as the incidence of the tail and host recognition protein is widespread 21 . φ 24 B has also been shown to survive well in compost models 22 , and infectivity after 30 days in bovine manure and slurry 23 . φ 24 B is genetically similar to phages isolated from sporadic outbreaks of STEC with high virulence and therefore a good model of the viruses circulating in E. coli populations in the environment 13 .
The ability of lambdoid like phages to increase virulence by carriage of toxins in their accessory genome 24 is well described e.g. the cholera toxin (CTX) carried by Vibrio cholera phage 25,26 . In Stx-phage genomes the shigatoxin genes are always located at the same position, upstream of the Q antiterminator gene therefore organisation and gene location is important 27 . It has been hypothesised that presence of stx also offers selection and stability for the lysogens 28 . Colon et al. 29 observed that Stx-prophages show greater levels of spontaneous induction than lambda but this more readily correlates to Rec dependent and independent control of the CI repressor protein rather than presence or absence of stx 29 .
There are other accessory genes that are seemingly superfluous to viral replication have been shown to aid microbial selection against environmental stress. Examples include: antibiotic resistance 30,31 , acid tolerance 32,33 and polylysogeny 34 . Phage gene expression has also been shown to aid adhesion and colonisation, for example the expression of the λ -encoded lom gene promotes adhesion to buccal epithelial cells 35 , and the λ -encoded bor confers serum resistance 36 . Other phage encoded virulence traits include exotoxin production in E. coli 37 and increase in bacterial invasion via Staphylococcus phage encoded kinase that influences fibrinolysis 38 . Bacteriophage φ 24 B has also been shown to encode a mi-RNA in the lom region that alters expression of anti-repressor d-ant and downstream activity of CI, leading to rapid induction 39 . Tree et al. 40 identified 55 prophage regions encoding small regulatory RNA within the Sakai E. coli O157:H7 strain 40 . These small prophage anti-sRNA had the ability to form complexes or mimic core genome regulatory sRNA to aid selective advantage to the bacterial host in bovine rectal mucus. Stx-phage φ 24 B shows 98% sequence homology at the nucleotide level to the remnant Stx2 phage present in the Sakai genome 13 .
The function of the large numbers of hypothetical proteins encoded by φ 24 B and other converting phage is difficult to determine, and a focus of this study, as gene expression or interaction may be specific to an environment or subject to selective pressure. Therefore, current approaches in vitro using synchronous cultures and standard laboratory conditions to investigate the role of these prophages are challenging. The function of these hypothetical gene products and how they impact the host in either an advantageous or deleterious way may be missed. In this study we focus on phage mediated antimicrobial tolerance to antibiotics found in the livestock farm setting which is the primary reservoir of pathogenic shigatoxigenic E. coli.
Our principal aim was to identify how infection and integration of φ 24 B changes the microbial physiology, and how phage conversion aids selection compared to its naïve counterpart. This study demonstrates an increase in cell proliferation in standard culture, increased respiration activity using the Biolog phenotypic array and antimicrobial tolerance that is not linked to an identifiable resistance gene-cassette. The Biolog phenotypic array is a good tool for rapidly identifying phenotypic changes in the respiration profiles of bacteria 41,42 . This study also presents an untargeted metabolomics approach to reveal key phage-mediated differences in metabolism directly linked to biotin and fatty acid synthesis. Interestingly, we show that integration of φ 24 B into its primary integration site located 250 bp upstream of the IntS gene 43 allows converted E. coli MC1061 to grow using alternative sources of phosphate compared to the naïve bacterial host. We show for the first time unique changes in bacterial metabolite profiles on phage conversion in response to sub-inhibitory concentrations of antimicrobial agents.

Results
φ24 B integration increases cell proliferation. Growth rates of bacteria can differ due to a range of environmental parameters. To investigate the impact of φ 24 B on E. coli, viable cell counts were determined during growth comparing E. coli B strain MC1061 to single and double lysogens, the latter integrated into separate locations in the MC1061 chromosome 43 . Under standard growth conditions, the single and double lysogens showed significantly higher cell numbers compared to the naïve MC1061 (> 200%, Fig. 1). This alongside a statistically significant increase in doubling time of 18 minutes for the single lysogen compared to 20 minutes for the naïve MC1061 (p < 0.006), calculated from each growth curve (data not shown, n = 9). As the cultures reached mid to late exponential growth, the differences in growth rates diminished (Fig. 1). The greatest difference in growth was identified in early growth (Fig. 1). This was supported by a shorter lag time in the single and double lysogen compared to MC1061 with a 0.5 and 1.8-fold increase respectively in cell number after the first hour of growth. Stationary phase in the double lysogen is achieved earlier compared to MC1061 and the single lysogen as nutrients are utilised rapidly alongside the accumulation of inhibitory components of growth. φ24 B integration alters utilisation of different mono-phosphates and inability to respire using β-D-Allose. To explore the single lysogen related differences in cell respiration during growth we used the Biolog Phenotype MicroArray. This determined functional changes in respiration and metabolism resulting from phage conversion over a 48 h period with recordings taken every 15 min. The lysogen acquired the ability to respire and grow utilising uridine-2-monophosphate (U-2-P) when compared to the naïve MC1061 (SI Fig. 1 panel A). Phage mediated subversion of pyrimidine and purine synthesis by lytic phages has previously been reported and will be discussed later. Conversely, integration of the phage inhibited the lysogens ability to use D-Allose for respiration. φ24 B integration alters resistance to osmotic stress or antimicrobials. Again using the Biolog phenotypic array the single lysogen is able to tolerate a range of antimicrobial agents that have both extracellular and intracellular targets (Fig. 2). The respiration curves derived for this experiment are provided in the supplementary information (SI -Figs 1 and 2). Tests showing differences in respiration profile were determined in the presence of 22 antimicrobials and 7 increases in salt concentration (SI Table 2). Of these 29 different tests, the lysogen showed a level of tolerance to 17 antimicrobials (SI Table 2). Data presented in Fig. 2 (n = 3) are comparisons of the area under the respiration curve illustrating those that were altered significantly. φ 24 B infection promotes tolerance to 8-hydroxyquinoline (P < 0.000), chloroxylenol (P < 0.0037), and cefmetazole (P < 0.0026), cefoxitin, (P < 0.015) cefemendole (P < 0.0239) and amoxicillin (P < 0.057). Integration of φ 24 B into the primary site 250 bp upstream of IntS inhibits respiration utilising B-D-allose. Lysogeny also limits cell respiration in the presence of oxolinic acid although this is linked to phage induction as the cellular target is DNA gyrase. Inhibition of DNA gyrase has been previously shown to stimulate temperate phages to the lytic life cycle as cellular stress stimulates RecA, LexA and proteolytic cleavage of the repressor protein promoting phage induction 44 . φ24B integration increases MC1061 tolerance to sub-inhibitory concentrations of chloroxylenol and 8-hydroxyquinoline. To better understand the level of antimicrobial tolerance of the single lysogen, we first determined sub-inhibitory concentrations (SIC) against both MC1061 and the lysogen that reduce cell growth by ~60%. The antimicrobials chloroxylenol, oxolinic acid and 8-hydroxyquinoline were selected to validate the Biolog data. Prior to comparison, an approximate SIC range was determined for MC1061 utilising each of the 3 test drugs. Figure 3 illustrates increased tolerance by the lysogen in the presence of chloroxylenol and 8-hydroxyquinoline. Conversely, the naïve host shows increased tolerance compared to the lysogen in the presence of oxolinic acid. This also offers a positive control for the assay as oxolinic acid targets DNA gyrase and therefore stimulates phage induction 26 . Phage induction was confirmed by the presence of free phage compared to the un-induced control (data not shown).
Metabolic profiles comparing naïve MC1061 to φ24 B Lysogen. We used an untargeted metabolite profiling approach using high resolution LC-MS (≤1 ppm mass accuracy in full scan) to determine metabolic differences between bacterial host and lysogen during growth and when challenged with a sub-inhibitory concentration of test antibiotic. To broadly compare findings, significant metabolic differences (p < 0.05) were observed between both growth phase and antimicrobial challenge. In total, > 11 K ion features or possible metabolites were determined across all of the different tests performed. Of these 81 showed discrimination between the naïve MC1061 and the φ 24 B lysogen that had clean chromatogram peaks and < 5% coefficient variable (CV) (SI Table 3). These 81 metabolites that show differences can be further stratified to each test. The metabolite data was analysed using supervised and non-supervised methods of multivariate analysis. Principal Component Analysis was first employed to visualise trends in the dataset and identify potential outliers. To further interrogate the data, Partial-Least Squared Discriminant Analysis models (PLS-DA) were generated and score plots are shown in (Fig. 4A-C). The PLS-DA models for both hydroxyquinoline and chloroxylenol conditions score plots had good discriminating ability, establishing the metabolic differences between the lysogen and naïve host. During standard growth conditions component 1 failed to discriminate: Q2 − 0.556, R2Y 0.262, as R2Y and Q2 < 0.5, although certain metabolites showed significant differences between the lysogen and MC1061.  Table 5. Stx-phage φ 24 B has been previously shown to undergo spontaneous induction 27 and may impact the metabolite profile through sequestration of host function and movement to lysis. We therefore compared the metabolite profiles of both the lysogen and MC1061 with a phage inducing agent, oxolinic acid (DNA gyrase inhibitor). No correlation was seen between metabolite profiles of the lysogen or MC1061 when compared to that of the lysogen undergoing induction with oxolinic acid (data not shown).

Figure 2. A comparison of Area Under the Respiration Curve (AURC) data from the Biolog bacterial
phenotypic microarray. Data plotted shows the addition of supplemented nutrients or chemical challenge showed statistically significant difference in rates of respiration between the lysogen and naïve MC1061 host (for P values see SI Table 2). Arbitrary Omnilog fluorescence values (y-axis) show differences between the naïve MC1061 (light grey) host and the lysogen (dark grey) over a 47. φ24B integration alters the metabolite profile of MC1061 in standard growth conditions. Out of the 81 discriminatory metabolites determined in this study, only 16 were shown to discriminate between the naïve host and single lysogen under standard culture conditions. Of these 16 metabolites 4 were found in higher levels in the lysogen. This suggests that the lysogen down regulates certain metabolic functions or is directing metabolism along a different pathway, or both. It is likely to support the change in biology we report in this work and increased rates of early growth by the lysogen.
Early growth in the lysogen demonstrates an observable difference in metabolic profile compared to the naïve MC1061. Under standard growth conditions during early growth, 5 metabolites in total were shown to discriminate between the naïve MC1061 and lysogen. Of these, 1 was higher compared to the naïve control ( Fig. 3 SI). During stationary phase, in standard growth conditions, only 9 metabolites in total showed significant difference and all were found in lower levels in the lysogen. As phage-mediated metabolic differences are present during standard culture, the differences in metabolite profiles under challenge with sub-inhibitory concentrations of 8-hydroxyquinoline and chloroxylenol were tested (Fig. 2B,C). The previous Biolog results showed that the lysogen displays a tolerance to these 2 antibiotics. φ24B integration alters the metabolite profile of MC1061 during growth under sub-inhibitory concentrations of 8-hydroxyquinoline. Upon treatment with 8-hydroxyquinoline, there were 29 metabolites that showed significant difference between the naïve MC1061 and single lysogen. Of these 29 metabolites, 22 were found in higher levels in the lysogen. Early growth phase in the lysogen demonstrates an observable difference in metabolite profile compared to naïve MC1061. Under 8-hydroxyquinoline stress during early growth, 6 metabolites in total were shown to discriminate between the naïve MC1061 and lysogen. Of these, 5 were higher compared to the naïve control (Fig. 3 SI).   metabolites in total were shown to discriminate between the naïve MC1061 and lysogen. Of these, 9 were higher compared to the naïve control (Fig. 3 SI).
Alteration in metabolomics profile and antimicrobial tolerance is not linked to kanamycin resistance selective marker used to detoxify φ24B. The kanamycin gene (aph3) used to detoxify the φ 24 B phage 18 is used as a selective marker only prior to experimentation. As there is no consistency in the metabolic profiles when the lysogen cultures are treated with the 2 different antimicrobials this cannot be driven by φ 24 B encoding aph3. The metabolite profiles are also discriminatory to each of the 2 antimicrobials tested.
Characterising the metabolites that discriminate between the naïve MC1061 and φ24 B lysogen. The discriminatory metabolites determined from each test were compared with metabolite databases and were putatively identified based on exact mass and empirical formula (LCMS methods section). The identity of each metabolite was confirmed using fragmentation analysis using a secondary MS/MS stage. Identities with fragment similarity were found for 58 of the 81 metabolites discriminating between the naïve and lysogen. We focus here on 6 particular metabolites as they have robust identities from fragmentation patterns, retention times, and low accurate mass error (PPM), relating to known curated bacterial metabolites (Table 4 SI The lysogen has significantly higher intensity levels of pimelic acid under all tests, specifically during early growth (Fig. 5). FAPy-Adenine, a bacterial stress marker 45 , is only seen in stressed conditions in these analyses, with the lysogen expressing significantly lower intensity during early growth and higher intensity at stationary phase growth (Fig. 5). Hexadecanoic acid is identified in significantly higher abundance under cellular stress of chloroxylenol, and is further increased in the lysogen during early growth (P = 0.04). Metabolite sphinganine is present under standard conditions in higher intensity in the naïve MC1061. When challenged with chloroxylenol, intensity levels of sphinganine were undetectable in both naïve and lysogen during early growth. During mid-exponential and stationary phase growth under choloroxylenol test there is > 100 fold increase in intensity of sphinganine in both the naïve and lysogen. 5-Methyluridine is present at stationary phase in all conditions, and is also identified in higher intensity when challenged with both antibiotics. Ophthalmic acid was present at all stages of growth under standard conditions where the lysogen shows lower intensity at early and mid-growth, and higher levels at stationary phase. When treated with either antimicrobial agent, ophthalmic acid was only present at stationary growth, with significantly higher intensity found in the lysogen (P = 0.001). During standard culture, there are 16 metabolites responsible for the differences seen between the core metabolic profiles of naïve host and lysogen during the 3 growth phases. Importantly 10 of these, including pimelic acid, are also present when the lysogen is challenged with chloroxylenol and 8-hydroxyquinoline.
In the absence of antibiotics, the metabolite profile shows less discrimination between the lysogen and host at the 3 stages of growth by PLS-DA (Fig. 4A). Changes in individual metabolite abundances were measured as before (Figs 4 and 6), and > 100 were deemed possible biologically relevant metabolites. From the confirmed compounds, a total of 16 metabolites (SI Table 3) were shown to discriminate between MC1061 and the φ 24 B lysogen.
We further analysed these data using Hierarchical cluster analysis (HCA) and Euclidean dissimilarity matrix (DM) to create a heatmap that discriminates between 81 metabolites across all tests in this study (Fig. 6). The unsupervised heatmap shows that the metabolic profiles have separated by condition. Figure 6 illustrates differences between the metabolic profiles of the naïve MC1061 and φ 24 B lysogen when comparing both test antimicrobials and the standard culture conditions. Firstly there is the greatest dissimilarity when the naïve host or lysogen has been treated with a sub-inhibitory concentration of chloroxylenol. Within this grouping the naïve host shows the greatest difference in profile at stationary phase for the treatment group. The chloroxylenol group is further stratified by whether the phage is present or absent. Presence of the phage offers the most dissimilar metabolic profile under this antimicrobial challenge. Treatment with 8-hydroxyquinoline has less impact on the metabolic profiles, yet the antimicrobial tolerance is still marked. The difference is also less marked as the profiles are stratified by growth phase rather than presence or absence of phage. Importantly, in Fig. 6 differences between the 81 metabolites in the naïve host and lysogen without challenging with an antimicrobial are still apparent.

Discussion
The accessory genome of bacteria promoted through horizontal gene transfer is important in understanding how mobile genetic elements aid selection in the environment. Metagenomics of DNA viruses in environmental and clinical samples has revealed a wide range of antimicrobial resistance genes (ARGs) [46][47][48] . Enault et al. 49 demonstrate that caution is needed as ARGs are over-estimated and therefore rarely found in phage genomes and that this over-estimation was further supported by functionality 49 . We here propose a different mechanism promoted by Stx-phage φ 24 B , through infection and subversion of the cell physiology, promoting tolerance to sub-inhibitory concentrations of antimicrobials 8-hydroxyquinoline and chloroxylenol. Importantly, we show here that this tolerance is to antimicrobials commonly used globally in the farming industry. Chloroxylenol is widely used in detergent based products l, for the direct treatment of livestock e.g. bovine teat treatment 50 . Similarly 8-hydroxyquinoline has a broad antimicrobial activity with agricultural use due to its potency against insects, fungi, and bacteria 51  infection on the physiology of the bacteria under antimicrobial pressure. The impact of prophage should not be underestimated as basis for metabolic variation and selection for the bacterial host by heightening or dampening cellular response to stress. With the altered metabolic profile of the φ 24 B lysogen and increased levels of biotin concentration and fatty acid intensities it leads to the hypothesis that altered growth and lipids may play a role in altering the cell surface that promotes antimicrobial exclusion.
The metabolite pimelic acid is a precursor for the majority of the carbon atoms of biotin 53 . Biotin plays a crucial role in cell metabolism via carboxylation and decarboxylation reactions. Beyond its function as a cofactor for carboxylases, biotin also plays a role in gene regulation in mammals 54 . Unfortunately the mechanism of its action in E. coli is relatively unknown. However it has been shown that the E. coli BioC-BioH pathway uses a methylation and demethylation strategy to complete the necessary pimeloyl moiety 53 . This methylation approach disguises the biotin synthetic intermediates such that they become substrates for the fatty acid synthetic pathway, and once the pimeloyl moiety is complete it is demethylated 53 . We show here that the φ 24 B prophage has a significant upregulatory effect on biotin that links to other physiological pathways including fatty acid synthesis 53 . Differences in pimelic acid intensities between lysogen and naïve host were greatest during early growth (Figs 1A,2B, 3A and 5A), which correlates with the differences observed in growth rates during the first 2.5 hours of culture (Fig. 1). This is associated with a ~3 fold increase in the level of biotin present per cell at mid-exponential growth phase  (Fig. 5B3), which correlates to metabolite profiling for pimelic acid. This is the first time a phage has been shown to drive the biotin pathway.
The Biolog data confirmed significant differences in rates of respiration between the naïve host and lysogen under different nutrient and chemical challenges. Interestingly the φ 24 B lysogens acquired the ability to respire using Uridine-2-Phosphate where the AURC is illustrated in Fig. 2A and SI Figure 1. U-2-P is involved in cellular metabolism (including biotin metabolism), nucleotide metabolism, pyrimidine metabolism and pyrimidine catabolism 55 . Phages have been shown to subvert purine and pyrimidine synthesis to aid viral construction and proliferation 52,56 . Lysogen mediated differences encoded by φ 24 B also supports these numerous studies as metabolomics profiling identifies increased pyrimidine catabolism as 5-methyluridine intensity decreases in the lysogen sample. It has also been previously shown through metagenomic analysis that well adapted phages of Pseudomonas aeruginosa in the lung carry genes that are involved with purine, pyrimidine and different phosphate utilisation 57 . This is further supported by Chevallereau, et al. 58 who show marked changes in RNA metabolism during bacteriophage infection of P. aeruginosa 58 . Importantly, not only does φ 24 B lysogen show increased pyrimidine utilisation, it also shows that phages can expand the group of phosphates E. coli can use for cell respiration and growth, in this instance U-2-P.
Conversely to addition of function, integration of φ 24 B into the MC1061 chromosome confers an inability to respire using β -D-allose. The lysogen used in this study has φ 24 B inserted into the primary integration site on the E. coli genome ~250 bp downstream of intS 45 . In E. coli there are 3 genes, alsB, alsA, and alsC, that are linked to the utilisation of D-allose 59 , but are disparate (~700 Kbp) from any of the 6 integration sites reported by Fogg et al. 43,45 . This is significant as it illustrates that integration can yield off target epigenetic effects. This study illustrates that the lysogen associated changes in fatty acid synthesis may exclude D-allose being transported into the bacterial cell, although the mechanism of this restriction is unknown.
Previous research showed infection with λ increased cell growth by the lysogen under cell starvation/supplementation of glucose in a chemostat culture [60][61][62] . This increase in growth rate is also seen with φ 24 B here. Interestingly we see a further increase with infection of a secondary, genetically identical phage. The double lysogen is an identical clone to that reported by Fogg et al. 43 , with phage integrating into the secondary integration site 45 . When monitoring growth the single lysogen confers a doubling time of 17 mins compared to 20 mins for the naïve MC1061. This is also combined with shorter lag phases in both the single and double lysogen. It has previously been shown in many bacteria and yeast that augmenting a growing culture with biotin increases cell growth rates [63][64][65][66][67] .
When stressed with chloroxylenol, a demonstrated increase in lipid biosynthesis occurs in the lysogen presumably through subversion of the biotin pathway. This is supported through identification of higher intensity levels of hexadecanoic acid in the metabolite data. Hexadecanoic acid is involved in the biosynthesis of lipid A, Figure 6. Heatmap generated by metabolic levels of 81 metabolites using HCA and DM Culture conditions and presence or absence of phage can be found alongside each profile (H = 8-hydroxyquinoline, C = chloroxylenol). Each individual tile represents a metabolite. The colour of a given tile denotes higher or lower intensity of the metabolite. The colour scale key is: dark blue: lowest levels; white: mid-point; dark red: highest level. The gradient between these colours represents variation in the levels of the metabolite across the colour scale (putative ID's can be found in SI Table 3). Pimelic acid is highlighted across all profiles with a hatched box.
Scientific RepoRts | 7:40424 | DOI: 10.1038/srep40424 a core outer cell membrane structure [68][69][70] . Changes in hexadecanoic acid in the lipid A structure of E. coli have been previously shown to be associated to mutations in the firA gene 69,71 . The firA gene is essential for growth and outer membrane synthesis 72 , and is essential for rifampicin resistance associated with certain mutations in the β subunit of the RNA polymerase 70 . This resistance and increase in hexadecanoic acid associated to the firA gene, shows that manipulation of this specific fatty acid likely improves antibiotic resistance. It is noteworthy that altering cell wall properties can broadly improve drug resistance 73 , and the biotin pathway is intrinsically linked to cell wall synthesis and growth 53,74 .
The lysogen showed increased tolerance to 8-hydroxyquinoline and chloroxylenol using the Biolog phenotypic array and sub-inhibitory antimicrobial tests. An untargeted metabolomics approach demonstrated that phage conversion offers the bacterial host different metabolic profiles to tolerate the two antimicrobials tested. The tolerance observed also suggests core functional changes allowing the cell to resist two disparate antimicrobials. This may suggest that these lysogen associated metabolic differences would likely infer tolerance to other environmental challenges and selective pressures.
Firstly, we present a metabolic difference in growth under standard culture conditions between the naïve MC1061 and the φ 24 B lysogen. From 81 metabolites, 16 were discriminatory between the lysogen and MC1061 (SI Table 3). Pimelic acid is present and constitutively raised after infection by φ 24 B . We also showed here a difference between metabolites at the 3 key stages of growth. These again differ between the lysogen and MC1061 (Fig. 6).
Under treatment with chloroxylenol in early growth, increased intensity of hexadecanoic acid was identified, particularly in the lysogen. The metabolite sphinganine was observed in our data, sphinganine plays an essential part in the sphingolipid synthesis pathway 75 . In both the lysogen and naïve host there is evidence of higher intensities of a sphinganine under chloroxylenol treatment and at stationary growth in standard conditions (Fig. 6). In Shigella species, a pathway associated with mammalian sphingolipid based rafts has been linked to improved binding and mammalian host cell entry 76 .
Metabolic differences between the lysogen and naïve bacteria are the most disparate when under challenge of a sub-inhibitory concentration of chloroxylenol, illustrated in the PLS-DA plots (Fig. 4) and heatmap (Fig. 6). Chloroxylenol is a bactericidal agent and a halophenol that targets microbial membranes 73 with a broad activity as an antimicrobial 77 .
We also demonstrate an increased tolerance by the lysogen to antimicrobial 8-hydroxyquinoline. Interestingly, compared to stress under chloroxylenol, the 8-hydroxyquinoline tested metabolite profile changes less significantly from standard conditions. Furthermore when treated with 8-hydroxyquinoline, the metabolite profile is less pronounced in the lysogen when compared to the chloroxylenol test. 8-hydroxyquinoline is a lipophilic metal-chelator with intracellular targets 78 . It inhibits growth by chelating metal ions, e.g. Zn 2+ on RNA polymerase 79,80 . The changes in the intensity of lipids present at the cell surface, that we previously suggested affect uptake of D-allose, are similarly likely to inhibit levels of these 2 antimicrobials entering the cell.
When testing cellular stress it is imperative to find markers of inhibition detailed in the metabolite data. The metabolomic profiles identified 2 discriminatory metabolites that are associated with cellular stress: FAPy-Adenine 47 and ophthalmic acid 81,82 . Our data showed that FAPy-Adenine was only present when cells were challenged by the antimicrobials 8-hydroxyquinoline and chloroxylenol. Interestingly the intensity levels of FAPy-adenine differ greatly depending on the antimicrobial used and also presence or absence of integrated φ 24 B (Fig. 5). In the presence of chloroxylenol the lysogen demonstrates lower intensities of the stress marker FAPy-adenine, 0.56 and 0.37 fold less in early and mid-exponential growth phase respectively (Figs 2a,b and 5). It also shows higher intensity of pimelic acid compared to MC1061. This strengthens the hypothesis of a biotin related lipid increase or change at the cell surface lowering levels of the drug reaching its intracellular target.
When challenging the culture with 8-hydroxyquinoline, FAPy-adenine intensity increases rapidly, even more so than the naïve host (Fig. 5). Again there is an increase in pimelic acid intensity that is ubiquitous to the metabolic profiles in the presence of an integrated φ 24 B . However hexadecanoic acid was undetectable within the cell, which may be associated with the lipophilic nature of the drug. The stress response occurs directly after treating with 8-hydroxyquinoline and likely promotes some cell death. Extracellular lipids released through cell lysis binds the drug, forming a matrix. This therefore would reduce the concentration of the available drug present allowing the bacterial culture to grow.
The second stress marker ophthalmic acid is an analogue of glutathione and a reported marker of oxidative stress 81,82 . Ophthalmic acid intensity mirrored the stationary phase levels of FAPy-adenine, across all tests, however it was also present in standard culture conditions in both the early and mid-exponential growth phase cultures. This observation from our metabolomic analysis implies higher oxidative stress in the lysogen at stationary growth, as the data is supported by Desnues et al. 82 . The oxidative stress also correlates with the reduction in growth rate and a reduction in the intensity of pimelic acid.
This study has established that Stx-phage φ 24 B provides a 'jump start' in early respiration and increased bacterial growth rates. These phage-mediated alterations in bacterial host metabolic profile may offer positive selection for the lysogen. Subversion of the biotin pathway is core to the changes mediated by φ 24 B as it links to the bacterium becoming able to tolerate chloroxylenol and 8-hydroxyquinoline during early and mid-exponential growth phase. These tolerances are important as both antimicrobials are used globally in livestock farming. Importantly, metabolic shift and subversion offers 2 mechanisms of controlling this antimicrobial tolerance through increased biotin and fatty acid synthesis. With treatment and tolerance to choloroxylenol, alteration in levels of Lipid A, and presumably other lipids enables exclusion of the drug from entry. Secondly 8-hydroxyquinoline treatment drives early cellular stress, cell death and lysis which increases extracellular lipids that bind free drug, allowing the community to continue to grow. The mechanism for this is unclear, yet this is importantly linked to subversion of the bacterial cell because no phage encoded metabolites are present. We therefore propose that temperate phages may not carry ARGs but play a larger role interfering with metabolic regulation that alters bacterial sensitivity to antimicrobials.

Materials and Methods
Bacterial strains and growth conditions -Buffer and Agar. All bacterial strains were grown in Lysogeny Broth + 0.01 M CaCl 2 (LB). Growth of the MC1061(φ 24 B ::Kan), a lysogen of the bacteriophage φ 24 B ::Kan and growth of the MC1061(φ 24 B ::Cat), a lysogen of the bacteriophage φ 24 B ::Cat was supplemented with 50 μ g.ml −1 kanamycin (kan) and Chloramphenicol (cat) respectively. Bottom agar plates for plaque assay included LB broth including 7% (w/v) grade 1 agar. Soft top agar was contained LB broth plus 0.4% (w/v) grade 1 agar. Unless otherwise stated culture conditions were at 37 °C, and broth cultures were shaken at 200 rpm.
Growth curve of single and double lysogens. A single colony of either naïve MC1061, single or double lysogen was cultured overnight for 18 h (200 rpm). LB with 0.01 M CaCl 2 (100 ml) was inoculated with 1% (v/v) of the overnight culture. Samples were taken over a 7 hour period, subject to serial, ten-fold dilutions and spread plated on LB agar plates.
Bacterial phenotypic microarray -Biolog. The Biolog assay utilises a redox dye where a tetrazolium violet salt acts as an electron receptor from the tricarboxylic cycle and reduction to NADH. The transfer alters the clear salt to a purple formazan dye that is inexplicably linked to the cellular activity, specifically cell respiration. An inoculum was taken from an 18 h streaked plate of either MC1061 or φ 24 B lysogen, raised through 2 rounds of passage from single colony amplification from cryo-stock. A single colony was added to fluid IF-0 (containing 50 μ M leucine due to MC1061's auxotrophy), to a transmittance of 42% T on a Biolog turbidometer in a 20 mm diameter tube as per manufacturer's instructions and used to inoculate Biolog Phenotypic Microarray plates.
The panel plates used for this study included Biolog plates PM 1-20, which include a range of both metabolic and toxicological additives (see SI). Further details of the components associated with these PM plates can be found at http://www.biolog.com/pmMicrobialCells.html. The Biolog PM plates were grown at 37 °C and monitored using the Omnilog plate reader at 30 min intervals over 47 hours.

Sub-inhibitory concentration (SIC) assay. An inoculum with a transmittance of 42 % T was created in
an identical manner to that performed for the Biolog assay for both MC1061 or the lysogen. An identical volume of inoculum was added to an equal volume of antimicrobial (double concentrate) diluted in LB broth. Readings were taken at 0 and 18 hours; plates were incubated at 37 °C. SIC assay -Liquid chromatography mass spectrometry (LCMS) analysis comparing metabolic compounds from Naïve host and Lysogens. Replicated in triplicate therefore n = 9 bacterial cultures (10 ml) were grown as previously described under standard growth conditions and challenged with antimicrobials, the cells were harvested at early, mid, and late log phase. The cells were harvested by centrifugation (5,000 rpm for 5 min) and the pellet washed (× 3) in ice cold 1 × PBS prior to lyophilisation. Lyophilised samples (× 3) were pooled, normalised for weight/vol (normalised to 1 mg.ml) with methanol and 0.1% formic acid, this was vortexed and then sonicated (Bandelin Sonorex, Sonicator) for 1 hour and centrifuged (5 k rpm for 5 mins). The supernatant was recovered and filtered through 0.22 μ m pore-sized, nylon filter and injected into the Q-Exactive LC-MS (Thermo-Fisher) after separation on a Phenomenex Gemini column (110 A, 150 × 2 mm, 5 μ m, flow 0.2 ml.min). LCMS mobile phase parameters were: 0-6 mins at 20% ACN, 8 mins 60% ACN, 12 mins 95% ACN, 17 mins 95% ACN, 17.1 − 23 mins 5% CAN. MS conditions were: full MS mode, resolution 70, 000, AGC target 1 × 10 6 , maximum IT 200 ms, scan range 150-2000, column temperature 35 °C. Progenesis QI software was used for raw data analysis; this software provided alignment, peak picking, pairwise statistical analysis and putative metabolite ID based on accurate mass. Metabolites were confirmed by analysing pure standards and MSMS fragmentation analysis run under identical analytical conditions. Biotin quantification assay. Inoculums were prepared in the same manner as for the SIC and Biolog assays. Optical density values were taken at 0 and 18 h, incubation at 37 °C. The cultures were diluted to the lowest OD 600 reading to normalise cell number between naïve MC1061 and lysogen. Dilutions of both naïve MC1061 and lysogen were made in LB (1:100 and 1:1000). Cells were harvested by centrifugation (5,000 rpm for 5 mins).
The biotin assay was completed using the Bio Vision ® (Cambridge, UK) Biotin Quantitation Kit (Colorimetric) according to the manufacturer's protocol. In brief, the supernatant was discarded and the pellet re-suspended in 10 μ l PBS and heated to 100 °C for 3 min and then immediately placed on ice. Diluted naïve and lysogen cells (10 μ l) were added to individual aliquots of Biotin Assay buffer (20 μ l) and 300 μ l of biotin reaction mix, pre-prepared as described in the Biotin Quantitation Kit protocol (version 7.6), was added to the buffer and cells. The mix was incubated at 21 °C for 15 min. Each sample mix (150 μ l) was then pipetted into a microtitre plate and read at 500 nm. A standard curve was prepared as per manufacturer's instructions.

Statistical analysis.
To determine statistically significant difference between growth rates of the single and double φ 24 B lysogen compared to MC1061, paired-sample T-tests in the statistical package SPSS was used. The two tailed p values are given at 95% confidence limits. The statistically significant difference in SIC between lysogen and MC1061 was calculated using an independent t-test, using the SPSS platform (> 95% confidence limits). The Biolog area under the respiration curve (AURC) for respiration values were calculated using a trapezoid algorithm. Statistical significance of area under the respiration curve (AURC) and comparison at a specific time point during mid exponential growth phase was achieved by determining normal Gaussian distribution by parametric analysis and statistical significance identified using an un-paired t test (> 95% confidence limits).
Scientific RepoRts | 7:40424 | DOI: 10.1038/srep40424 Metabolomic analysis was carried out primarily by Progenesis QI software, this software provided alignment, peak picking, pairwise statistical analysis and putative metabolite ID based on accurate mass. Further multivariate analysis was performed using SIMCA-P. ID's were obtained through the QI plugin 'Progenesis metascope' and filtered through a range of databases using sdf files (ECMDB, HMDB, small molecules drugs, Biomolecules, analgesics mix, Lipid MBD, Basic lipids, and Yeast DB). A paired sample t test was used to determine statistically significant differences between intensities of metabolites identified during metabolomic analyses (> 95% confidence limits).