Introduction Over the last decade, the field of systems vaccinology has emerged, in which high throughput transcriptomics and other omics assays are used to probe changes of the innate and adaptive immune system in response to vaccination. The goal of this study was to benchmark key technical and analytical parameters of RNA sequencing (RNA-seq) in the context of a multi-site, double-blind randomized vaccine clinical trial. Methods We collected longitudinal peripheral blood mononuclear cell (PBMC) samples from 10 subjects before and after vaccination with a live attenuated Francisella tularensis vaccine and performed RNA-Seq at two different sites using aliquots from the same sample to generate two replicate datasets (5 time points for 50 samples each). We evaluated the impact of (i) filtering lowly-expressed genes, (ii) using external RNA controls, (iii) fold change and false discovery rate (FDR) filtering, (iv) read length, and (v) sequencing depth on differential expressed genes (DEGs) concordance between replicate datasets. Using synthetic mRNA spike-ins, we developed a method for empirically establishing minimal read-count thresholds for maintaining fold change accuracy on a per-experiment basis. We defined a reference PBMC transcriptome by pooling sequence data and established the impact of sequencing depth and gene filtering on transcriptome representation. Lastly, we modeled statistical power to detect DEGs for a range of sample sizes, effect sizes, and sequencing depths. Results and Discussion Our results showed that (i) filtering lowly-expressed genes is recommended to improve fold-change accuracy and inter-site agreement, if possible guided by mRNA spike-ins (ii) read length did not have a major impact on DEG detection, (iii) applying fold-change cutoffs for DEG detection reduced inter-set agreement and should be used with caution, if at all, (iv) reduction in sequencing depth had a minimal impact on statistical power but reduced the identifiable fraction of the PBMC transcriptome, (v) after sample size, effect size (i.e. the magnitude of fold change) was the most important driver of statistical power to detect DEG. The results from this study provide RNA sequencing benchmarks and guidelines for planning future similar vaccine studies.

Francisella tularensis (F. tularensis) is an intracellular pathogen that causes a potentially debilitating febrile illness known as tularemia. F. tularensis can be spread by aerosol transmission and cause fatal pneumonic tularemia. If untreated, mortality rates can be as high as 30%. To study the host responses to a live-attenuated tularemia vaccine, peripheral blood mononuclear cell (PBMC) samples were assayed from 10 subjects collected pre- and post-vaccination, using both the 2D-DIGE/MALDI-MS/MS and LC-MS/MS approaches. Protein expression related to antigen processing and presentation, inflammation (PPARγ nuclear receptor), phagocytosis, and gram-negative bacterial infection was enriched at Day 7 and/or Day 14. Protein candidates that could be used to predict human immune responses were identified by evaluating the correlation between proteome changes and humoral and cellular immune responses. Consistent with the proteomics data, parallel transcriptomics data showed that MHC class I and class II-related signals important for protein processing and antigen presentation were up-regulated, further confirming the proteomic results. These findings provide new biological insights that can be built upon in future clinical studies, using live attenuated strains as immunogens, including their potential use as surrogates of protection.

Tularemia is a highly infectious and contagious disease caused by the bacterium Francisella tularensis. To better understand human response to a live-attenuated tularemia vaccine and the biological pathways altered post-vaccination, healthy adults were vaccinated, and plasma was collected pre- and post-vaccination for longitudinal lipidomics studies. Using tandem mass spectrometry, we fully characterized individual lipid species within predominant lipid classes to identify changes in the plasma lipidome during the vaccine response. Separately, we targeted oxylipins, a subset of lipid mediators involved in inflammatory pathways. We identified 14 differentially abundant lipid species from eight lipid classes. These included 5-hydroxyeicosatetraenoic acid (5-HETE) which is indicative of lipoxygenase activity and, subsequently, inflammation. Results suggest that 5-HETE was metabolized to a dihydroxyeicosatrienoic acid (DHET) by day 7 post-vaccination, shedding light on the kinetics of the 5-HETE-mediated inflammatory response. In addition to 5-HETE and DHET, we observed pronounced changes in 34:1 phosphatidylinositol, anandamide, oleamide, ceramides, 16:1 cholesteryl ester, and other glycerophospholipids; several of these changes in abundance were correlated with serum cytokines and T cell activation. These data provide new insights into alterations in plasma lipidome post-tularemia vaccination, potentially identifying key mediators and pathways involved in vaccine response and efficacy.

Background: Tularemia is a potential biological weapon due to its high infectivity and ease of dissemination. This study aimed to characterize the innate and adaptive responses induced by two different lots of a live attenuated tularemia vaccine and compare them to other well-characterized viral vaccine immune responses. Methods: Microarray analyses were performed on human peripheral blood mononuclear cells (PBMCs) to determine changes in transcriptional activity that correlated with changes detected by cellular phenotyping, cytokine signaling, and serological assays. Transcriptional profiles after tularemia vaccination were compared with yellow fever [YF-17D], inactivated [TIV], and live attenuated [LAIV] influenza. Results: Tularemia vaccine lots produced strong innate immune responses by Day 2 after vaccination, with an increase in monocytes, NK cells, and cytokine signaling. T cell responses peaked at Day 14. Changes in gene expression, including upregulation of STAT1, GBP1, and IFIT2, predicted tularemia-specific antibody responses. Changes in CCL20 expression positively correlated with peak CD8+ T cell responses, but negatively correlated with peak CD4+ T cell activation. Tularemia vaccines elicited gene expression signatures similar to other replicating vaccines, inducing early upregulation of interferon-inducible genes. Conclusions: A systems vaccinology approach identified that tularemia vaccines induce a strong innate immune response early after vaccination, similar to the response seen after well-studied viral vaccines, and produce unique transcriptional signatures that are strongly correlated to the induction of T cell and antibody responses.

The immune response to live-attenuated Francisella tularensis vaccine and its host evasion mechanisms are incompletely understood. Using RNA-Seq and LC–MS on samples collected pre-vaccination and at days 1, 2, 7, and 14 post-vaccination, we identified differentially expressed genes in PBMCs, metabolites in serum, enriched pathways, and metabolites that correlated with T cell and B cell responses, or gene expression modules. While an early activation of interferon α/β signaling was observed, several innate immune signaling pathways including TLR, TNF, NF-κB, and NOD-like receptor signaling and key inflammatory cytokines such as Il-1α, Il-1β, and TNF typically activated following infection were suppressed. The NF-κB pathway was the most impacted and the likely route of attack. Plasma cells, immunoglobulin, and B cell signatures were evident by day 7. MHC I antigen presentation was more actively up-regulated first followed by MHC II which coincided with the emergence of humoral immune signatures. Metabolomics analysis showed that glycolysis and TCA cycle-related metabolites were perturbed including a decline in pyruvate. Correlation networks that provide hypotheses on the interplay between changes in innate immune, T cell, and B cell gene expression signatures and metabolites are provided. Results demonstrate the utility of transcriptomics and metabolomics for better understanding molecular mechanisms of vaccine response and potential host–pathogen interactions.

Francisella tularensis is a highly infectious Gram-negative bacterium that is the etiologic agent of tularemia in animals and humans and a Tier 1 select agent. The natural incidence of pneumonic tularemia worldwide is very low; therefore, it is not feasible to conduct clinical efficacy testing of tularemia medical countermeasures (MCM) in human populations. Development and licensure of tularemia therapeutics and vaccines need to occur under the Food and Drug Administration's (FDA's) Animal Rule under which efficacy studies are conducted in well-characterized animal models that reflect the pathophysiology of human disease. The Tularemia Animal Model Qualification (AMQ) Working Group is seeking qualification of the cynomolgus macaque (Macaca fascicularis) model of pneumonic tularemia under Drug Development Tools Qualification Programs with the FDA based upon the results of studies described in this manuscript. Analysis of data on survival, average time to death, average time to fever onset, average interval between fever and death, and bacteremia; together with summaries of clinical signs, necropsy findings, and histopathology from the animals exposed to aerosolized F tularensis Schu S4 in five natural history studies and one antibiotic efficacy study form the basis for the proposed cynomolgus macaque model. Results support the conclusion that signs of pneumonic tularemia in cynomolgus macaques exposed to 300-3,000 colony forming units (cfu) aerosolized F tularensis Schu S4, under the conditions described herein, and human pneumonic tularemia cases are highly similar. Animal age, weight, and sex of animals challenged with 300-3,000 cfu Schu S4 did not impact fever onset in studies described herein. This study summarizes critical parameters and endpoints of a well-characterized cynomolgus macaque model of pneumonic tularemia and demonstrates this model is appropriate for qualification, and for testing efficacy of tularemia therapeutics under Animal Rule.

The ubiquitin-proteasome system has a central role in the degradation of intracellular proteins and regulates a variety of functions. Viruses belonging to several different families utilize or modulate the system for their advantage. Here we showed that the proteasome inhibitors MG132 and epoxomicin blocked a postentry step in vaccinia virus (VACV) replication. When proteasome inhibitors were added after virus attachment, early gene expression was prolonged and the expression of intermediate and late genes was almost undetectable. By varying the time of the removal and addition of MG132, the adverse effect of the proteasome inhibitors was narrowly focused on events occurring 2 to 4 h after infection, the time of the onset of viral DNA synthesis. Further analyses confirmed that genome replication was inhibited by both MG132 and epoxomicin, which would account for the effect on intermediate and late gene expression. The virus-induced replication of a transfected plasmid was also inhibited, indicating that the block was not at the step of viral DNA uncoating. UBEI-41, an inhibitor of the ubiquitin-activating enzyme E1, also prevented late gene expression, supporting the role of the ubiquitin-proteasome system in VACV replication. Neither the overexpression of ubiquitin nor the addition of an autophagy inhibitor was able to counter the inhibitory effects of MG132. Further studies of the role of the ubiquitin-proteasome system for VACV replication may provide new insights into virus-host interactions and suggest potential antipoxviral drugs.