Identification of systemic immune response markers through metabolomic profiling of plasma from calves given an intra-nasally delivered respiratory vaccine
© Gray et al.; licensee BioMed Central. 2015
Received: 27 May 2014
Accepted: 12 December 2014
Published: 14 February 2015
Vaccination procedures within the cattle industry are important disease control tools to minimize economic and welfare burdens associated with respiratory pathogens. However, new vaccine, antigen and carrier technologies are required to combat emerging viral strains and enhance the efficacy of respiratory vaccines, particularly at the point of pathogen entry. New technologies, specifically metabolomic profiling, could be applied to identify metabolite immune-correlates representative of immune protection following vaccination aiding in the design and screening of vaccine candidates. This study for the first time demonstrates the ability of untargeted UPLC-MS metabolomic profiling to identify metabolite immune correlates characteristic of immune responses following mucosal vaccination in calves. Male Holstein Friesian calves were vaccinated with Pfizer Rispoval® PI3 + RSV intranasal vaccine and metabolomic profiling of post-vaccination plasma revealed 12 metabolites whose peak intensities differed significantly from controls. Plasma levels of glycocholic acid, N-[(3α,5β,12α)-3,12-Dihydroxy-7,24-dioxocholan-24-yl]glycine, uric acid and biliverdin were found to be significantly elevated in vaccinated animals following secondary vaccine administration, whereas hippuric acid significantly decreased. In contrast, significant upregulation of taurodeoxycholic acid and propionylcarnitine levels were confined to primary vaccine administration. Assessment of such metabolite markers may provide greater information on the immune pathways stimulated from vaccine formulations and benchmarking early metabolomic responses to highly immunogenic vaccine formulations could provide a means for rapidly assessing new vaccine formulations. Furthermore, the identification of metabolic systemic immune response markers which relate to specific cell signaling pathways of the immune system could allow for targeted vaccine design to stimulate key pathways which can be assessed at the metabolic level.
The vaccination of farm animals against endemic, genetically evolving and emerging pathogens is important not only to ensure animal health, but also to reduce the costs associated with disease losses, either clinical or subclinical. Successful vaccination leads to the production of specific T and B cell effector immune responses that assist in the control of infection within the animal. This results in the generation of virus neutralizing antibodies that recognize the pathogen, specific effector T-cell responses and the development of an immune memory response helping to protect against future exposure to the infection. Continuous vaccine development is required to address virus evolution, new and emerging viral threats and to improve vaccine efficacy against currently managed pathogens. However, assessment of new vaccine candidates, adjuvant and novel vaccine carrier systems (such as nanoparticles) using animal trials is extremely expensive and can take several months to years (investigation of short term immune responses, vaccine-challenge studies and field trials and long term immune protection against natural wildtype virus challenge) . Although the expenditure associated with animal vaccine development is not fully disclosed, the estimated budget required to develop a single FDA approved vaccine for human studies is estimated in the region of $1-2 billion . The majority of these costs are attributed to the high failure rate of vaccine candidates/formulations, with only 1 in every 10 000 vaccine formulations gaining approval by the FDA . New vaccine candidates may initially be assessed in mice models, however this does not necessarily translate to performance in the target species and animal trials can have poor efficacy  or fail to induce immune protection at all  as a result of species differences in immune systems . Furthermore, many factors such: as immune system maturity, vaccine delivery route, concurrent infections, poor nutrition and the presence of maternally derived antibodies can affect vaccine efficacy . As the expenditure associated with candidate vaccine development escalates with clinical trial progression, rapid vaccine screening methods which can assess candidate vaccine effectiveness at early trial stages in vivo are required to minimize financial outlay and improve the speed of vaccine development pipelines.
In the agricultural industry the Bovine Respiratory Disease (BRD) complex is considered to be one of the most significant causes of economic loss in intensively reared cattle worldwide. The associated losses accruing from the elevated mortality and poor growth performance of infected animals  coupled with the need for costly therapeutic interventions have a significant negative impact on farm incomes. This disease complex is estimated to result in an annual total economic loss to the US agri-food industry of over $2 billion, with treatment and preventative costs approaching $3 billion . Whilst vaccination against the infectious agents involved in BRD pathogenesis is currently employed to manage the disease [9,10], it has not significantly reduced BRD prevalence or severity. Furthermore, it has been observed that some animals fail to develop immune protection despite vaccine treatment, becoming infected with each new seasonal BRD outbreak . Therefore, the development of new methods for screening BRD pathogen vaccine candidates, along with the clear understanding of the host immune response, would provide significant economic and animal welfare benefits to the agricultural industry by speeding up the vaccine development.
Current methods for determining viral vaccine effectiveness in mammalian species include assessment of neutralizing antibody titre [12,13], PCR-based detection of viral load in infected tissues at post-immunization challenge , expression of pro- and anti-inflammatory cytokines , duration of viral shedding after post-immunization challenge , and mortality/morbidity findings at post-mortem . However, the quantitative analysis of pro- and anti-inflammatory cytokines is generally considered to be one of the important indicators of vaccine efficacy, but typically only occurs after preliminary candidate antigen or formulation screening , and may require a number of vaccines and challenge studies to fully elucidate the important factors associated with a protective immune response. A greater understanding of not only the proteomic responses to vaccination [16,17], but also the changes at the metabolic level in bio-fluids, such as plasma or serum, may help develop early marker-based screening tests. These can be used to assess the potential of new vaccine candidates and also may provide additional information to improve the rationale behind future vaccine approaches taking into account vaccine delivery routes, vaccine carrier technologies and animal genetics.
Untargeted metabolomic profiling refers to the top down global analysis of small molecules or “metabolites” with molecular weights less than 1 kDa and is typically performed by Nuclear Magnetic Resonance (NMR) or mass spectrometry (MS) techniques, with the later offering greater sensitivity. Ultra Performance Liquid Chromatography (UPLC) coupled to a highly accurate mass spectrometer allows for the separation and determination of elemental composition (and hence the identification) of 1000s of metabolites in a single sample. As metabolites are often end-stage products of biological processes, their analysis can provide a greater understanding of the physiological state of an organism at the time of sampling compared to genomic or proteomic analysis. The application of metabolomic techniques to the analysis of bio-fluids has been found to offer interesting new insights to the understanding of the physiological processes involved in disease development and diagnosis [18-21]. However, to date little emphasis has been placed on assessing the in vivo metabolomic changes induced by immunization procedures and foreign antigen exposure.
This study for the first time reports the metabolomic responses in plasma of calves vaccinated with an intranasal vaccine (Pfizer Rispoval® PI3 + RSV) composed of the modified live viruses Bovine Respiratory Syncytial Virus (BRSV) and Bovine Parainfluenza Virus-3 (BPI3V). BRSV and BPI3V are two of the major viral pathogens involved in the Bovine Respiratory Disease (BRD) complex . The aims of the current study were to identify plasma metabolomic markers of the primary and secondary immune response in calves to administration of a commercial intranasal vaccine containing live modified virus. Importantly, it is well known that vaccination via the parenteral route leads to poor protection at mucosal surfaces. Therefore the studies undertaken here allow vaccine immunology and biomarkers to be identified through the natural route of infection for these respiratory viruses in order to identify the most suitable host defense mechanisms. This is the first untargeted UPLC-MS based metabolomics study to reveal metabolite markers of the immune response occurring during vaccination which may have potential use in the early screening of candidate vaccine effectiveness and aid new vaccine design by increasing the understanding of the processes involved in immune response development post vaccination.
Materials and methods
Chemicals and reagents
HPLC grade acetone was purchased from Sigma Aldrich (Dorset, UK). LC-MS grade acetonitrile, water, methanol and chloroform were purchased from Fisher Scientific (Loughborough, UK). All analytical standards for metabolite confirmation were purchased from Sigma Aldrich.
Experimental design and sample collection
600 μL of plasma was added to 2.4 mL of ice cold acetone, vortexed for 30 s and placed on ice for 15 min. The sample was then deproteinated by centrifugation at 4000 g at 4 °C for 20 min. 2.4 mL of supernatant was removed and dried under nitrogen for 45 min at 40 °C using TurboVap LV (Caliper Life Sciences, Hopkinton, USA). Resulting residue was reconstituted in 500 μL of Ultra Pure H2O and liquid/liquid extraction of lipids performed by addition of 500 μL of ice-cold Methanol:Chloroform (1:1 v/v) and vortexing for 30 s followed by centrifugation at 16 000 g at 4 °C for 15 min. liquid/liquid extraction was repeated and after centrifugation 900 μL of the aqueous layer was removed and dried under nitrogen using TurboVap at 40 °C for 45 min. The residue was reconstituted in 300 μL Ultra Pure H2O and filtered by centrifugation at 10 000 g at 4 °C using 0.22 μm Constar Spin-X® Centrifuge Tube Filter for 5 min.
UPLC-MS analysis of Bovine Plasma
UPLC-MS analysis was performed using an Acquity UPLC system coupled to a XEVO G2 Q-TOF (Waters Corporation, Milford, MA, USA. 8 μL of prepared sample extracts were injected onto an Acquity UPLC HSS-T3 C18 column (100 mm × 2.1 mm i.d., 1.8 μm; Waters Corporation, Milford, MA, USA). Column and autosampler temperature were maintained at 50 °C and 10 °C respectively. Chromatographic separation was carried out at a flow rate of 600 μL/min with mobile phase consisting of 99.9% H2O/0.1% Formic Acid (A) and 99.9% Acetonitrile/0.1% Formic Acid (B). The elution gradient was as follows: 0–2 min isocratic at 1% of B, 2–14.5 min linear gradient from 1-100% of B, 14.5-17 min isocratic at 100% of B, 17–17.5 min linear gradient from 100-1% of B and finally 17.5-20 min isocratic at 1% of B. Mass spectrometry was performed using a Waters XEVO G2 QTOF operating in positive-ion mode (ESI+) with the capillary voltage set to 1500 V and the sampling cone voltage 30 V. The desolvation and cone gas flows were set at 750 L/h and 100 L/h respectively. Source and desolvation temperatures were 120 °C and 400 °C respectively. Leucine Enkephalin ([M +H]+= 278.1141 Da, and [M+ H] = 556.2771 Da) was used for accurate lockmass calibration during data acquisition. Lockmass acquisition settings were: 0.5 s scan time, 30 s interval, 3 scan average, mass window +/− 0.5 Da. Collision energy was only applied on function 2, with ramping between 15 V and 30 V and centroid data was acquired using positive-ion, resolution mode.
Mass spectrometry data processing
Total Ion Count (TIC) chromatograms and spectra were acquired with MassLynx version 4.1 (Waters) in centroid format and metabolite data was processed using MarkerLynx software (Waters Corporation, Milford, MA, USA). The MarkerLynx method for data extraction and de-convolution was as follows: Ions were extracted from function 1 data using peak detection analysis of retention time window 0.30-14.50 min, with a mass range of 100 Da to 1200 Da. The extracted-ion chromatogram (XIC) window for data collection was 0.2 Da and apex peak tracking parameters were set to automatic with no smoothing. Data collection parameters consisted of an intensity threshold (counts) of 500 and a mass window of 0.02 Da with a retention time window of 0.20 min. A noise elimination level of 6 was applied and isotopes were removed. Peak heights for extracted ions were normalized against the total peak height of all extracted ions and standardized to a total ion count of 10 000. The results were exported in .csv format as a two dimensional data table in which rows and columns respectively represented analysed samples and the relative normalized peak heights of each detected mass spectrometric signal, i.e. as an Accurate Mass (m/z) and Retention Time (min) Pair (AMRTP).
Temporal changes in BPI3V antibody titre were analysed using two-tailed paired T-test, and significant differences between treatment groups at sampling stages was assessed using two-tailed heteroscedastic T-test. SIMCA-P+ version 13.0 (Umetrics, Sweden) was used for multivariate metabolite marker selection. For SIMCA-P+ analysis, data was prefiltered to exclude AMRTPs with coefficient of variation greater than 50% in replicate quality control pools. All centroid data were Pareto scaled and analysed by unsupervised Principle Component Analysis (PCA) and supervised discriminatory analysis by Orthogonal Projections of Latent Structures-Discriminant Analysis (OPLS-DA). Unsupervised PCA models were generated at each sampling day to reveal potential relationships between treatment groups. Supervised analysis by OPLS-DA was performed to reveal potential markers of response to treatment in vaccinated calves compared to non-vaccinated calves at each sampling day. Robustness of final OPLS-DA discriminative models was assessed by setting a predictive model of each case in which 2/3 of the data (known treatment) was used to predict the remaining 1/3 (unknown treatment).
Identification of potential metabolite markers
The elemental composition of selected parent compounds was determined in MassLynx using both positive and negative mode data. Mass uncertainty was set to 3mDa, odd and even electron state, carbon isotope filter of +/− 5% and elements included were C, H, O, N, P and S. Where applicable Na and K adduct elemental composition were determined with the respective element included in the analysis parameters. Elemental compositions were searched against PubChem and Chemspider online databases, and where possible Function 2 fragments were matched against Metlin, HMDB or Massbank databases. Where fragmentation spectra for the analyte in question was not available, in silico fragmentation was performed using Metfrag and Function 2 fragmentation data was validated against potential in silico fragments. Compounds revealed by database matching or in silico fragmentation had their identities confirmed against commercial analytical standards by UPLC-MS. Pooled plasma samples and individual standards (1 μM) were analysed under identical UPLC and mass spectrometric run conditions as utilized previously, and metabolite identities confirmed by matching retention time and Function 1 and Function 2 spectra (including low and high energy fragments and adducts).
Haematological and immunological responses to vaccination
The haematology results are illustrated in Additional file 1. Red Blood Cell count (RCC), White Blood Cell count (WCC), Hemoglobin content (HG), Mean Corpuscular Volume (MCV), Mean Corpuscular Haemoglobin (MCH), Packed Cell Volume (PCV) and Platelet count (PLT) measurements for all calves were within the specified reference parameters based on data from Brun-Hansen et al.  for calves at corresponding ages which are known to differ from adults. However, significant differences in lymphocyte and neutrophil count were observed between study groups at a number of points throughout the vaccination stages of the study whilst remaining within reference parameters for healthy animals. Figure 2B reveals significantly (p < 0.05) higher lymphocyte counts at day 42 (2 weeks post-booster immunization) in vaccinated animals compared to non-vaccinated. Furthermore, lymphocyte counts in vaccinated animals were found to rise significantly (p < 0.05) from day 14 to 42 and decreased from day 49 to 63, with no significant temporal changes observed in non-vaccinated animals. Figure 2C illustrates that neutrophil counts were significantly higher at day 0 (p < 0.01) and day 14 (p < 0.05) and significantly lower (p < 0.05) at day 42 in vaccinated compared to non-vaccinated animals. Neutrophil counts decreased significantly (p < 0.05) in vaccinated animals from day 14 to 42 and increased significantly (p < 0.05) in non-vaccinated animals from day 14 to 42. At post-BPI3V challenge stages, WCC, RCC and lymphocyte counts were outside reference parameters. There was a significant (p < 0.05) rise in WCC in non-vaccinated animals after BPI3V infection from days 70 to day 75 and WCC was higher than reference parameters from days 75 to 90 in vaccinated calves, and significantly (p < 0.05) higher at day 90 in vaccinated calves compared to non-vaccinated. RCC was elevated outside reference parameters throughout the study and was significantly (p < 0.05) higher in vaccinated calves at day 75 compared to non-vaccinated. Whilst No significant differences were observed in lymphocyte counts between study groups, lymphocyte counts were above the reference parameters in both study groups at day 75. There was a significant increase in lymphocyte counts in non-vaccinated animals from day 70 to 75 and a significant decrease from days 75 between days 82, with no significant variations observed in lymphocyte counts of vaccinated animals. There were no observed significant differences in neutrophil counts between study groups at corresponding sampling points post-infection, however, significant (p < 0.05) temporal variations were observed with a decrease from day 70 to 75 in non-vaccinated animals, and an increase from day 75 to 82 in both study groups. Despite being within the reference parameters, significant differences were observed in PCV at day 75, MCV at day 82 and MCH at day 75 between study groups.
Selection of metabolomic markers of vaccine response via multivariate data analysis
Identification of parent ions and AMRTPs
Final AMRPTs selected as metabolomic markers of response to vaccination
Retention time (mins)
Measured mass (Da)
Elemental composition (M + H)
% Fit confidence
Mass Error (mDa)
Uric Acid (UA)a
Hippuric Acid (HA) a
N-Methylhippuric Acid (NMHA)
Hexahydrohippuric Acid (HHA)a
Glycocholic Acid (GCA)a
Taurodeoxycholic Acid (TDCA)a
Vaccine development strategies have typically focused on the induction of specific antibody responses to antigenic molecules as an indicator of effective immune stimulation (whole virus or viral proteins) with little understanding of how antigens elicit particular immune responses . Whilst this approach has resulted in the majority of vaccines available today, new vaccine development is increasingly costly and it is imperative that the least promising candidates are eliminated from investigation as early as possible. Whilst the identification of immune correlations, in particular proteins such as interleukins and cytokines, has helped to improve the understanding of vaccine mediated immune response mechanisms, emerging metabolomic techniques may offer a novel perspective to vaccine development and selection. This is possible as metabolites are the ultimate end stage products or mediators of biological processes and can provide a holistic view of underlying physiological processes. This study has for the first time employed an untargeted UPLC-MS metabolomic profiling approach to identify plasma components altered by primary and secondary immune responses to an intranasal vaccination containing live modified virus. Vaccination with Pfizer Rispoval® PI3 + RSV resulted in an increase of BPI3V specific IgG antibody titre in vaccinated calves. The IgG responses indicated that the booster immunization was necessary to further increase the titres of BPI3V antibodies. This stimulated IgG response is characteristic of the immune responses typically seen in two dosage vaccination experiments. Furthermore, the high IgG titre maintained in vaccinated animals following experimental challenge with BPI3V indicates a primed immune response to infection. Conversely the non-vaccinated animals showed no significant alterations in BPI3V specific IgG titre during vaccination stages but at post-BPI3V challenge displayed a rise in WCC indicating a stimulated immune response in the absence of prior infection. We can therefore characterize the markers identified in this study which showed altered abundance between vaccinated and non-vaccinated study groups as being related to the immunological responses involved in primary and secondary immune-stimulation. Metabolomic profiling initially revealed a panel of 96 features which differed significantly between calves vaccinated intranasally with Pfizer Rispoval® PI3 + RSV and non-vaccinated calves throughout the study. These features were further refined to identify 12 high confidence parent compounds, 10 of which were identified based on database searching and spectral matching, with 6 further validated using pure standards. Based on the observed stimulated IgG response following two intranasal vaccine dosages (characteristic of the immune response typically seen in two-stage dosage vaccine programs ), these markers can be characterized as being related to the responses involved in primary and secondary immune stimulation. All detected markers were shown to be present in plasma prior to vaccine administration, therefore their altered abundance reflects the processes of stimulation of the mucosal immune system following exposure of animals to modified live viral vaccine, and are not associated with vaccine breakdown products.
A significant increase in BLD (10.26 FC) plasma peak intensity in vaccinated calves compared to non-vaccinated calves was observed two weeks after the second Pfizer Rispoval® PI3 + RSV dosage (day 49). BLD is degraded from heme groups by Heme Oxygenases (HO)  and rapidly reduced to Bilirubin by Biliverdin Reductase . BLD and the HO isoform HO-1 are involved in promoting anti-inflammatory responses. BLD suppresses IL-2 production and T cell proliferation , and activates the aryl hydrocarbon receptor (AHR)  regulating T helper cell phenotype towards an anti-inflammatory regulatory T cell, or a pro-inflammatory T helper 17 effector function [29,30]. HO-1, produced from circulating monocytes during acute inflammatory states, has immunomodulatory function by supporting the proliferative capacity and activation of CD4(+) and CD8(+) T cells via CD14(+) monocytes, preventing DC maturation [31-33] and controlling Th1 pro-inflammatory cytokine production [34,35]. The observed increased plasma levels of BLD may be due to either increased HO-1 activity at a rate which cannot be matched by Biliverdin Reductase, or an unknown modulatory effect on Biliverdin Reductase affecting the rate at which it reduces BLD to bilirubin - strongly supporting BLD and HO-1 involvement in the promotion of anti-inflammatory responses as the immune stimulation response following vaccine dosage is deactivated to normal status.
Plasma peak intensities of a number of conjugated Bile Acids (BAs) were found to differ significantly between vaccinated and non-vaccinated calves, with increased TDCA (4.87 FC) after primary vaccine dosage (day 7, prior to observable differences in IgG titre) and increased GCA (1.63 and 2.12 FC) and NDGCA (1.80 and 3.29 FC) after secondary dosage (day 42 and 49 respectively) in vaccinated animals. The up-regulation of different BA conjugates at different stages following vaccination suggests that Cholic acid conjugation may allow for mediation of different immune response pathways, dependent on the amino acid that is conjugated. To date the only other published study assessing metabolomic responses to vaccination reported an up-regulation of serum cholesterol levels at days 5 and 10 in rabbits after primary and secondary immunization (using human red blood cells) respectively . As cholesterol is a precursor to Cholic acid, and hence GCA and TDCA, it’s up-regulation could be due to an increased demand for immuno-regulatory BAs. BAs have recently been demonstrated to function as signaling molecules with immunomodulatory effects resulting from their activation of the nuclear receptor Farsenoid X Receptor (FXR) [36,37] and the G-protein coupled receptors Transmembrane G protein coupled receptor 5 (TGR5) [38,39] and the Formyl-peptide receptors (FPR) . TGR5, present in a variety of cells of the immune system including alveolar macrophages and dendritic cells [38,41], can be activated by BAs to stimulate cAMP production which inhibits pro-inflammatory cytokines secretion (TNF-α, IL-1β, IL-1α, IL-6 and IL-8) . This activation favours more hydrophobic BAs , with TDCA more likely to activate TGR5 in vivo compared to GCA or NDGCA and its’ up-regulation at day 7 in vaccinated calves illustrates a potential role in arresting innate inflammatory response in favour of the progression of adaptive immunity through the activation of TGR5. Dendritic cells demonstrate FXR modulated function following BA treatment, resulting in increased production of IFN-γ producing T-cells and reduced airway eosinophilia and macrophage influx after OVA airway stimulation . Activation of FXR by synthetic bile acid ligands also results in NK T-cell inhibition and reduced pro-inflammatory osteopontin production which indicates FXRs role in maintaining immune homeostasis . Both TDCA and GCA have been shown to activate FXR in vivo, upon transportation into FXR expressing cells . However, the consequences of BA driven FXR modulation would favour an adaptive immune response, which would therefore be associated with secondary booster dosage stages where increased plasma levels of GCA and NDGCA were observed in vaccinated animals. Significant up-regulation of circulating plasma levels of GCA could provide a method by which pro-inflammatory responses at the site of infection are limited after secondary exposure to viral infections. GCA induced activation of FXR at post-booster stages in vaccinated animals may reduce macrophage and neutrophil influx to the site of infection and increase the levels of INF-γ producing CD4+ and CD8+ T-cells, promoting an adaptive immune response.
Significantly higher levels of PLC (2.21 FC) in vaccinated relative to non-vaccinated calves were observed at day 14 post initial vaccine dosage, suggesting a potential role of PLC in either mediating adaptive immune responses following primary immune stimulation, or as an end-stage metabolite produced as a result of these processes. PLC has been demonstrated to modulate innate immune response following transplantation by preventing neutrophil, CD4+, CD8+, ED1+ and MHCII cell infiltration upon transplantation, protecting grafts from oxidative stress injury . In vivo PLC treatment resulted in significant decrease in the release of leukocyte adhesion molecules into plasma (E-selectin, P-selectin, L-selectin ICAM-1 and VCAM-1)  suggesting PLC’s action in modulating leukocyte function and trafficking, with a potential role in reducing immune cell recruitment to sites of active viral replication with the progression of adaptive immunity associated with increased BPI3V IgG production at day 14 in this study. The immunomodulatory effects of PLC are further demonstrated by its’ potential to reduce platelet activating factors when exposed to neutrophils in vitro and in vivo . The significant up-regulation of PLC at day 14 may therefore be due to its’ anti-inflammatory effect in arresting pro-inflammatory leukocyte recruitment in favour of adaptive immune response.
UA released from injured cells activates the NALP3 inflammasome leading to IL-1β production , but the specific mechanism by which UA promotes its immunosuppressive effects is not fully understood [49,50]. Significantly higher levels of UA (1.57 FC) were observed in vaccinated relative to non-vaccinated calves at day 42 (1 week post second vaccine dosage) and would suggest an increased population of cytotoxic T cells responding to re-exposure to the live attenuated vaccine. Whilst there has been no previously reported evidence of the direct immunosuppressive properties of HA or HHA, significantly lower levels of HA (21.70 FC) were observed in vaccinated animals at day 49 (2 weeks post booster immunization), and significant differences in the plasma peak intensity of HHA in vaccinated animals compared to non-vaccinated at days 7 (−1.82 FC), 14 (4.45 FC) post initial immunization and days 42 (1.63 FC), 49 (− 3.41 FC) and 63 (2.17 FC), 1, 2 and 4 weeks post booster immunization respectively. Reduced plasma HA levels in vaccinated calves compared to non-vaccinated calves, suggests that it does not exert a direct anti-inflammatory effect. HA may instead be involved in the mediation of immune response through its’ metabolism to produce compounds which exert a direct immuno-modulatory effect. Hippuric acid is an acyl-glycine, formed through the action of N-acyltransferase involving Acyl-CoA. Circulating HA in mammals is primarily derived from the metabolism and glycine conjugation of phenolic dietary constituents . The metabolites HHA, NMHA and NCG have acyl-glycine like fragmentation patterns, with HHA and NCG levels found to significantly correlate with HA. As dietary composition was consistent between study groups and significant variation in HA occurred only at post-booster stages, the reduced plasma HA levels may result from reduced N-acyltransferase activity (due to increased metabolic demand for acyl-CoA) coinciding with a marked reduction in plasma levels of other metabolites (HHA and NCG).
In conclusion, this study has demonstrated the ability of untargeted UPLC-MS metabolomics to differentiate between vaccinated and non-vaccinated animals via the profiling of metabolite constituents in plasma. A number of plasma constituents were found to be associated with initial and/or booster vaccine dosage events and were therefore characterized as corresponding to primary and/or secondary immune stimulation. The markers associated with primary immune stimulation to vaccine administration, particularly TDCA and HHA which varied in plasma prior to observable variations in BPI3V IgG, may find application as early diagnostic markers to screen for immunogenic candidates during vaccine selection investigations, particularly for vaccines targeting routes of pathogen entry such as through the respiratory tract. The metabolites GCA, NDGCA, UA and BLD, which were altered within plasma at the secondary vaccine dosage stages, may additionally represent key diagnostic markers of established immune protection in animals. A number of the selected markers have demonstrated immunological properties, affecting key cell signaling pathways and receptors located on cells of the immune system. Further validation studies are required to expand upon this initial proof of concept study to discern whether the differences observed in these metabolite marker levels as a result of vaccination are repeatable in a larger external sample set comprising animals of varying ages, genders and breeds coupled with more frequent sampling. However, this study’s findings illustrate the potential for untargeted metabolomics to identify key metabolite immune correlates which relate to specific cell signaling pathways involved in immune responses. The assessment of such metabolite markers may provide greater information on the immune pathways stimulated (compared to traditional IgG or IgA serology) by vaccine formulations, and in particular benchmarking early metabolomic responses to highly immunogenic vaccine formulations could provide a means of rapidly assessing new vaccine candidates. Furthermore, a greater understanding of the metabolite pathways altered as a result of vaccination could allow for targeted vaccine design (antigens, adjuvants and nanoparticle carrier systems) to stimulate key immunogenic pathways which can be assessed at the metabolite level. Currently, metabolomics does not lend itself easily to routine vaccine candidate screening, due to the requirement of costly analytical equipment and skilled operators. However, recent advances in portable mass-spectrometers, such as the 908 Devices handheld system (908 Devices, Boston, USA), or the more user friendly Microsaic 4000 MiD® (Microsaic Systems, Surrey, UK), may eventually see the use of metabolites as routine markers for the evaluation of vaccine candidate immune function and immunogenicity. Future research will seek to expand upon the physiological role of these metabolomic markers and focus on establishing the biological pathways in which they may play a role in immuno-modulatory activities.
Bovine respiratory disease
Bovine Parainfluenza Type3 Virus
Ultra Performance liquid chromatography-mass spectrometry
Accurate mass retention time pair
Principle component analysis
Orthogonal projections to latent structures-discriminant analysis
Tandem mass spectrometry
N- N-Methylhippuric acid
Farsenoid X Receptor
Transmembrane G-Protein Coupled Receptor 5
NACHT, LRR and PYD domains-containing protein 3
We acknowledge the staff at AFBI Animal Services Unit VSD (Stormont) for their help in the animal study and the staff at the Advanced ASSET Center IGFS QUB for their excellent technical assistance. This research was funded by a Department of Agriculture and Rural Development (DARD) Northern Ireland postgraduate studentship awarded to Darren Gray.
- Zamani Z, Arjmand M, Tafazzoli M, Gholizadeh A, Pourfallah F, Sadeghi S, Mirzazadeh R, Mirkhani F, Taheri S, Iravani A, Bayat P, Vahabi F (2011) Early detection of immunization: a study based on an animal model using 1H nuclear magnetic resonance spectroscopy. Pak J Biol Sci 14:195–203View ArticlePubMedGoogle Scholar
- Light DW, Andrus JK, Warburton RN (2009) Estimated research and development costs of rotavirus vaccines. Vaccine 27:6627–6633View ArticlePubMedGoogle Scholar
- Oyston P, Robinson K (2012) The current challenges for vaccine development. J Med Microbiol 61:889–894View ArticlePubMedGoogle Scholar
- Agnandji ST, Lell B, Soulanoudjingar SS, Fernandes JF, Abossolo BP, Conzelmann C, Methogo BGNO, Doucka Y, Flamen A, Mordmueller B, Issifou S, Kremsner PG, Sacarlal J, Aide P, Lanaspa M, Aponte JJ, Nhamuave A, Quelhas D, Bassat Q, Mandjate S, Macete E, Alonso P, Abdulla S, Salim N, Juma O, Shomari M, Shubis K, Machera F, Hamad AS, Minja R, et al. (2011) First results of phase 3 trial of RTS, S/AS01 malaria vaccine in African children. N Engl J Med 365:1863–1875View ArticlePubMedGoogle Scholar
- Tameris MD, Hatherill M, Landry BS, Scriba TJ, Snowden MA, Lockhart S, Shea JE, McClain JB, Hussey GD, Hanekom WA, Mahomed H, McShane H (2013) MVA85A 020 Trial Study Team: Safety and efficacy of MVA85A, a new tuberculosis vaccine, in infants previously vaccinated with BCG: a randomised, placebo-controlled phase 2b trial. Lancet 381:1021–1028View ArticlePubMedGoogle Scholar
- Tameris MD, Hatherill M, Landry BS, Scriba TJ, Snowden MA, Lockhart S, Shea JE, McClain JB, Hussey GD, Hanekom WA, Mahomed H, McShane H, MVA85A 020 Trial Study Team (2013) Safety and efficacy of MVA85A, a new tuberculosis vaccine, in infants previously vaccinated with BCG: a randomised, placebo-controlled phase 2b trial. Lancet 381:1021–1028View ArticlePubMedGoogle Scholar
- Siegrist C (2008) Vaccine immunology, 1st edition. Saunders, Vaccines, pp 17–36Google Scholar
- Cusack PMV, McMeniman N, Lean IJ (2003) The medicine and epidemiology of bovine respiratory disease in feedlots. Aust Vet J 81:480–487View ArticlePubMedGoogle Scholar
- Snowder GD, Van Vleck LD, Cundiff LV, Bennett GL, Koohmaraie M, Dikeman ME (2007) Bovine respiratory disease in feedlot cattle: phenotypic, environmental, and genetic correlations with growth, carcass, and longissimus muscle palatability traits. J Anim Sci 85:1885–1892View ArticlePubMedGoogle Scholar
- Fulton RW (2009) Bovine respiratory disease research (1983–2009). Anim Health Res Rev 10:131–139View ArticlePubMedGoogle Scholar
- Taylor JD, Fulton RW, Lehenbauer TW, Step DL, Confer AW (2010) The epidemiology of bovine respiratory disease: what is the evidence for preventive measures? Can Vet J 51:1351–1359PubMed CentralPubMedGoogle Scholar
- Kahrs RF (2001) Viral Diseases Of Cattle, 2nd edition. Iowa State University Press, IowaGoogle Scholar
- Ploquin A, Szecsi J, Mathieu C, Guillaume V, Barateau V, Ong KC, Wong KT, Cosset FL, Horvat B, Salvetti A (2013) Protection against henipavirus infection by use of recombinant adeno-associated virus-vector vaccines. J Infect Dis 207:469–478View ArticlePubMedGoogle Scholar
- Galli V, Simionatto S, Marchioro SB, Fisch A, Gomes CK, Conceicao PR, Dellagostin OA (2012) Immunisation of mice with Mycoplasma hyopneumoniae antigens P37, P42, P46 and P95 delivered as recombinant subunit or DNA vaccines. Vaccine 31:135–140View ArticlePubMedGoogle Scholar
- Skoberne M, Cardin R, Lee A, Kazimirova A, Zielinski V, Garvie D, Lundberg A, Larson S, Bravo FJ, Bernstein DI, Flechtner JB, Long D (2013) An adjuvanted herpes simplex virus 2 subunit vaccine elicits a T cell response in mice and is an effective therapeutic vaccine in Guinea pigs. J Virol 87:3930–3942View ArticlePubMed CentralPubMedGoogle Scholar
- Martelli P, Ardigo P, Ferrari L, Morganti M, De Angelis E, Bonilauri P, Luppi A, Guazzetti S, Caleffi A, Borghetti P (2013) Concurrent vaccinations against PCV2 and PRRSV: study on the specific immunity and clinical protection in naturally infected pigs. Vet Microbiol 162:558–571View ArticlePubMedGoogle Scholar
- Adamczyk-Poplawska M, Markowicz S, Jagusztyn-Krynicka EK (2011) Proteomics for development of vaccine. J Proteomics 74:2596–616View ArticlePubMedGoogle Scholar
- DelVecchio VG, Sabato MA, Trichilo J, Dake C, Grewal P, Alefantis T (2010) Proteomics for the development of vaccines and therapeutics. Crit Rev Immunol 30:239–254View ArticlePubMedGoogle Scholar
- Zhang T, Wu X, Ke C, Yin M, Li Z, Fan L, Zhang W, Zhang H, Zhao F, Zhou X, Lou G, Li K (2013) Identification of potential biomarkers for ovarian cancer by urinary metabolomic profiling. J Proteome Res 12:505–512View ArticlePubMedGoogle Scholar
- Yang J, Zhao X, Liu X, Wang C, Gao P, Wang J, Li L, Gu J, Yang S, Xu G (2006) High performance liquid chromatography-mass spectrometry for metabonomics: potential biomarkers for acute deterioration of liver function in chronic hepatitis B. J Proteome Res 5:554–561View ArticlePubMedGoogle Scholar
- Graham SF, Chevallier OP, Roberts D, Hoelscher C, Elliott CT, Green BD (2013) Investigation of the Human brain metabolome to identify potential markers for early diagnosis and therapeutic targets of Alzheimer’s disease. Anal Chem 85:1803–1811View ArticlePubMedGoogle Scholar
- Xiao JF, Varghese RS, Zhou B, Ranjbar MRN, Zhao Y, Tsai T-H, Di Poto C, Wang J, Goerlitz D, Luo Y, Cheema AK, Sarhan N, Soliman H, Tadesse MG, Ziada DH, Ressom HW (2012) LC-MS based serum metabolomics for identification of hepatocellular carcinoma biomarkers in Egyptian cohort. J Proteome Res 11:5914–5923PubMed CentralPubMedGoogle Scholar
- Jeffery H, Davis SS, O’Hagan DT (1993) The preparation and characterization of poly(lactide-co-glycolide) microparticles. II. The entrapment of a model protein using a (water-in-oil)-in-water emulsion solvent evaporation technique. Pharm Res 10:362–368View ArticlePubMedGoogle Scholar
- Tammen H, Schulte L, Hess R, Menzel C, Kellmann M, Mohring T, Schulz-Knappe P (2005) Peptidomic analysis of human blood specimens: comparison between plasma specimens and serum by differential peptide display. Proteomics 5:3414–3422View ArticlePubMedGoogle Scholar
- Brun-Hansen HC, Kampen AH, Lund A (2006) Hematologic values in calves during the first 6 months of life. Vet Clin Pathol 35:182–187View ArticlePubMedGoogle Scholar
- Abraham NG, Kappas A (2008) Pharmacological and clinical aspects of heme oxygenase. Pharmacol Rev 60:79–127View ArticlePubMedGoogle Scholar
- Ryter SW, Alam J, Choi AMK (2006) Heme oxygenase-1/carbon monoxide: From basic science to therapeutic applications. Physiol Rev 86:583–650View ArticlePubMedGoogle Scholar
- Yamashita K, McDaid J, Ollinger R, Tsui TY, Berberat PO, Usheva A, Csizmadia E, Smith RN, Soares MP, Bach FH (2004) Biliverdin, a natural product of heme catabolism, induces tolerance to cardiac allografts. FASEB J 18:765–767PubMedGoogle Scholar
- Phelan D, Winter GM, Rogers WJ, Lam JC, Denison MS (1998) Activation of the Ah receptor signal transduction pathway by bilirubin and biliverdin. Arch Biochem Biophys 357:155–163View ArticlePubMedGoogle Scholar
- Quintana FJ, Basso AS, Iglesias AH, Korn T, Farez MF, Bettelli E, Caccamo M, Oukka M, Weiner HL (2008) Control of T-reg and TH17 cell differentiation by the aryl hydrocarbon receptor. Nature 453:65–71View ArticlePubMedGoogle Scholar
- Veldhoen M, Hirota K, Westendorf AM, Buer J, Dumoutier L, Renauld J-C, Stockinger B (2008) The aryl hydrocarbon receptor links T(H)17-cell-mediated autoimmunity to environmental toxins. Nature 453:106–109View ArticlePubMedGoogle Scholar
- Yachie A, Toma T, Mizuno I, Okamoto H, Shimura S, Ohta K, Kasahara Y, Koizumi S (2003) Heme oxygenase-1 production by peripheral blood monocytes during acute inflammatory illnesses of children. Exp Biol Med 228:550–556Google Scholar
- Chauveau C, Remy S, Royer PJ, Hill M, Tanguy-Royer S, Hubert FX, Tesson L, Brion R, Beriou G, Gregoire M, Josien R, Cuturi MC, Anegon I (2005) Heme oxygenase-1 expression inhibits dendritic cell maturation and proinflammatory function but conserves IL-10 expression. Blood 106:1694–1702View ArticlePubMedGoogle Scholar
- Burt TD, Seu L, Mold JE, Kappas A, McCune JM (2010) Naive human T cells are activated and proliferate in response to the heme oxygenase-1 inhibitor tin mesoporphyrin. J Immunol 185:5279–5288View ArticlePubMed CentralPubMedGoogle Scholar
- Poss KD, Tonegawa S (1997) Heme oxygenase 1 is required for mammalian iron reutilization. Proc Natl Acad Sci U S A 94:10919–10924View ArticlePubMed CentralPubMedGoogle Scholar
- Kapturczak MH, Wasserfall C, Brusko T, Campbell-Thompson M, Ellis TM, Atkinson MA, Agarwal A (2004) Heme oxygenase-1 modulates early inflammatory responses: evidence from the heme oxygenase-1-deficient mouse. Am J Pathol 165:1045–1053View ArticlePubMed CentralPubMedGoogle Scholar
- Makishima M, Okamoto AY, Repa JJ, Tu H, Learned RM, Luk A, Hull MV, Lustig KD, Mangelsdorf DJ, Shan B (1999) Identification of a nuclear receptor for bile acids. Science 284:1362–1365View ArticlePubMedGoogle Scholar
- Parks DJ, Blanchard SG, Bledsoe RK, Chandra G, Consler TG, Kliewer SA, Stimmel JB, Willson TM, Zavacki AM, Moore DD, Lehmann JM (1999) Bile acids: natural ligands for an orphan nuclear receptor. Science 284:1365–1368View ArticlePubMedGoogle Scholar
- Kawamata Y, Fujii R, Hosoya M, Harada M, Yoshida H, Miwa M, Fukusumi S, Habata Y, Itoh T, Shintani Y, Hinuma S, Fujisawa Y, Fujino M (2003) A G protein-coupled receptor responsive to bile acids. J Biol Chem 278:9435–9440View ArticlePubMedGoogle Scholar
- Hylemon PB, Zhou H, Pandak WM, Ren S, Gil G, Dent P (2009) Bile acids as regulatory molecules. J Lipid Res 50:1509–1520View ArticlePubMed CentralPubMedGoogle Scholar
- Le YY, Murphy PM, Wang JM (2002) Formyl-peptide receptors revisited. Trends Immunol 23:541–548View ArticlePubMedGoogle Scholar
- Ichikawa R, Takayama T, Yoneno K, Kamada N, Kitazume MT, Higuchi H, Matsuoka K, Watanabe M, Itoh H, Kanai T, Hisamatsu T, Hibi T (2012) Bile acids induce monocyte differentiation toward interleukin-12 hypo-producing dendritic cells via a TGR5-dependent pathway. Immunology 136:153–162View ArticlePubMed CentralPubMedGoogle Scholar
- Willart MAM, van Nimwegen M, Grefhorst A, Hammadt H, Moons L, Hoogsteden HC, Lambrecht BN, KleinJan A (2012) Ursodeoxycholic acid suppresses eosinophilic airway inflammation by inhibiting the function of dendritic cells through the nuclear farnesoid X receptor. Allergy 67:1501–1510PubMedGoogle Scholar
- Mencarelli A, Renga B, Migliorati M, Cipriani S, Distrutti E, Santucci L, Fiorucci S (2009) The bile acid sensor farnesoid X receptor is a modulator of liver immunity in a rodent model of acute hepatitis. J Immunol 183:6251–6261View ArticlePubMedGoogle Scholar
- Wang H, Chen J, Hollister K, Sowers LC, Forman BM (1999) Endogenous bile acids are ligands for the nuclear receptor FXR BAR. Mol Cell 3:543–553View ArticlePubMedGoogle Scholar
- Azzollini N, Cugini D, Cassis P, Pezzotta A, Gagliardini E, Abbate M, Arduini A, Peschechera A, Remuzzi G, Noris M (2008) Propionyl-L-carnitine prevents early graft dysfunction in allogeneic rat kidney transplantation. Kidney Int 74:1420–1428View ArticlePubMedGoogle Scholar
- Signorelli SS, Malaponte G, Di Pino L, Digrandi D, Pennisi G, Mazzarino MC (2001) Effects of ischaemic stress on leukocyte activation processes in patients with chronic peripheral occlusive arterial disease: role of L-propionyl carnitine administration. Pharmacol Res 44:305–309View ArticlePubMedGoogle Scholar
- Triggiani M, Oriente A, Golino P, Gentile M, Battaglia C, Brevetti G, Marone G (1999) Inhibition of platelet-activating factor synthesis in human neutrophils and platelets by propionyl-L-carnitine. Biochem Pharmacol 58:1341–1348View ArticlePubMedGoogle Scholar
- Laiakis EC, Morris GAJ, Fornace AJ, Howie SRC (2010) Metabolomic analysis in severe childhood pneumonia in the Gambia, West Africa: findings from a pilot study. PLoS One 5:e12655View ArticlePubMed CentralPubMedGoogle Scholar
- Kanellis J, Kang DH (2005) Uric acid as a mediator of endothelial dysfunction, inflammation, and vascular disease. Semin Nephrol 25:39–42View ArticlePubMedGoogle Scholar
- Kang DH, Park SK, Lee IK, Johnson RJ (2005) Uric acid-induced C-reactive protein expression: Implication on cell proliferation and nitric oxide production of human vascular cells. J Am Soc Nephrol 16:3553–3562View ArticlePubMedGoogle Scholar
- Wikoff WR, Anfora AT, Liu J, Schultz PG, Lesley SA, Peters EC, Siuzdak G (2009) Metabolomics analysis reveals large effects of gut microflora on mammalian blood metabolites. Proc Natl Acad Sci U S A 106:3698–3703View ArticlePubMed CentralPubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.