Metabolomic Signatures Differentiate Immune Responses in Avian Influenza Vaccine Recipients.

BACKGROUND: Avian influenza viruses pose significant risk to human health. Vaccines targeting the hemagglutinin of these viruses are poorly immunogenic without the use of adjuvants. METHODS: Twenty healthy men and women (18-49 years of age) were randomized to receive two doses of inactivated influenza A/H5N1 vaccine alone (IIV) or with AS03 adjuvant (IIV-AS03) one month apart. Urine and serum samples were collected on day 0 and on days 1, 3, and 7 following first vaccination and subjected to metabolomics analyses to identify metabolites, metabolic pathways, and metabolite clusters associated with immunization. RESULTS: Seventy-three differentially abundant (DA) serum and 88 urine metabolites were identified for any post-vaccination day comparison. Pathway analysis revealed enrichment of tryptophan, tyrosine and nicotinate metabolism in urine and serum among IIV-AS03 recipients. Increased urine abundance of 4-vinylphenol sulfate on Day 1 was associated with serologic response based on hemagglutination inhibition responses. In addition, 9 DA urine metabolites were identified in participants with malaise compared to those without. CONCLUSIONS: Our findings suggest that tryptophan, tyrosine, and nicotinate metabolism are upregulated among IIV-AS03 recipients compared with IIV alone. Metabolites within these pathways may serve as measures of immunogenicity and may provide mechanistic insights for adjuvanted vaccines.

Phase 1 Open-Label Dose Escalation Trial for the Development of a Human Bacillus Calmette-Guérin Challenge Model for Assessment of Tuberculosis Immunity In Vivo.

BACKGROUND: A controlled human infection model for assessing tuberculosis (TB) immunity can accelerate new vaccine development. METHODS: In this phase 1 dose escalation trial, 92 healthy adults received a single intradermal injection of 2 × 106 to 16 × 106 colony-forming units of Bacillus Calmette-Guérin (BCG). The primary endpoints were safety and BCG shedding as measured by quantitative polymerase chain reaction, colony-forming unit plating, and MGIT BACTEC culture. RESULTS: Doses up to 8 × 106 were safe, and there was evidence for increased BCG shedding with dose escalation. The MGIT time-to-positivity assay was the most consistent and precise measure of shedding. Power analyses indicated that 10% differences in MGIT time to positivity (area under the curve) could be detected in small cohorts (n = 30). Potential biomarkers of mycobacterial immunity were identified that correlated with shedding. Transcriptomic analysis uncovered dose- and time-dependent effects of BCG challenge and identified a putative transcriptional TB protective signature. Furthermore, we identified immunologic and transcriptomal differences that could represent an immune component underlying the observed higher rate of TB disease incidence in males. CONCLUSIONS: The safety, reactogenicity, and immunogenicity profiles indicate that this BCG human challenge model is feasible for assessing in vivo TB immunity and could facilitate the vaccine development process. CLINICAL TRIALS REGISTRATION: NCT01868464 (ClinicalTrials.gov).

Human Antibody Responses Following Vaccinia Immunization Using Protein Microarrays and Correlation With Cell-Mediated Immunity and Antibody-Dependent Cellular Cytotoxicity Responses.

BACKGROUND: There are limited data regarding immunological correlates of protection for the modified vaccinia Ankara (MVA) smallpox vaccine. METHODS: A total of 523 vaccinia-naive subjects were randomized to receive 2 vaccine doses, as lyophilized MVA given subcutaneously, liquid MVA given subcutaneously (liquid-SC group), or liquid MVA given intradermally (liquid-ID group) 28 days apart. For a subset of subjects, antibody-dependent cellular cytotoxicity (ADCC), interferon-γ release enzyme-linked immunospot (ELISPOT), and protein microarray antibody-binding assays were conducted. Protein microarray responses were assessed for correlations with plaque reduction neutralization titer (PRNT), enzyme-linked immunosorbent assay, ADCC, and ELISPOT results. RESULTS: MVA elicited significant microarray antibody responses to 15 of 224 antigens, mostly virion membrane proteins, at day 28 or 42, particularly WR113/D8L and WR101H3L. In the liquid-SC group, responses to 9 antigens, including WR113/D8L and WR101/H3L, correlated with PRNT results. Three were correlated in the liquid-ID group. No significant correlations were observed with ELISPOT responses. In the liquid-ID group, WR052/F13L, a membrane glycoprotein, correlated with ADCC responses. CONCLUSIONS: MVA elicited antibodies to 15 vaccinia strain antigens representing virion membrane. Antibody responses to 2 proteins strongly increased and significantly correlated with increases in PRNT. Responses to these proteins are potential correlates of protection and may serve as immunogens for future vaccine development. CLINICAL TRIALS REGISTRATION: NCT00914732.

Cytomegalovirus Genetic Diversity Following Primary Infection.

BACKGROUND: Infection with multiple cytomegalovirus (CMV) strains (mixed infection) was reported in a variety of hosts. As the virus genetic diversity in primary CMV infection and the changes over time remain incompletely defined, we examined CMV diversity and changes in diversity over time in healthy adolescent females who participated in a phase 2 CMV gB/MF59 vaccine trial. METHODS: CMV genetic diversity was determined by genotyping of 5 genes-gB (UL55), gH (UL75), gN (UL73), US28, and UL144-in urine, saliva, and plasma samples from 15 study subjects. RESULTS: At the time of primary infection, 5 of 12 (42%) urine samples had multiple virus strains, and 50% of vaccine recipients were infected with gB1 genotype (vaccine strain). Mixed infection was documented in all 15 subjects within 3 months after primary infection, and the majority had different CMV genotypes in different compartments. Changes in genotypes over time were observed in all subjects. CONCLUSIONS: Infection with multiple CMV genotypes was common during primary infection and further diversification occurred over time. Infection with gB1 genotype in vaccine recipients suggests a lack of strain-specific protection from the vaccine. As only 5 polymorphic genes were assessed, this study likely underestimated the true genetic diversity in primary CMV infection.

AS03-Adjuvanted H5N1 Avian Influenza Vaccine Modulates Early Innate Immune Signatures in Human Peripheral Blood Mononuclear Cells.

BACKGROUND: Adjuvant System 03 (AS03) markedly enhances responses to influenza A/H5N1 vaccines, but the mechanisms of this enhancement are incompletely understood. METHODS: Using ribonucleic acid sequencing on peripheral blood mononuclear cells (PBMCs) from AS03-adjuvanted and unadjuvanted inactivated H5N1 vaccine recipients, we identified differentially expressed genes, enriched pathways, and genes that correlated with serologic responses. We compared bulk PBMC findings with our previously published assessments of flow-sorted immune cell types. RESULTS: AS03-adjuvanted vaccine induced the strongest differential signals on day 1 postvaccination, activating multiple innate immune pathways including interferon and JAK-STAT signaling, Fcγ receptor (FcγR)-mediated phagocytosis, and antigen processing and presentation. Changes in signal transduction and immunoglobulin genes predicted peak hemagglutinin inhibition (HAI) titers. Compared with individual immune cell types, activated PBMC genes and pathways were most similar to innate immune cells. However, several pathways were unique to PBMCs, and several pathways identified in individual cell types were absent in PBMCs. CONCLUSIONS: Transcriptomic analysis of PBMCs after AS03-adjuvanted H5N1 vaccination revealed early activation of innate immune signaling, including a 5- to 8-fold upregulation of FcγR1A/1B/1C genes. Several early gene responses were correlated with HAI titer, indicating links with the adaptive immune response. Although PBMCs and cell-specific results shared key innate immune signals, unique signals were identified by both approaches.

Molecular Evolution of Herpes Simplex Virus 2 Complete Genomes: Comparison between Primary and Recurrent Infections.

Herpes simplex virus 1 (HSV-1) and HSV-2 are large, double-stranded DNA viruses that cause lifelong persistent infections characterized by periods of quiescence and recurrent disease. How HSV evolves within an infected individual experiencing multiple episodes of recurrent disease over time is not known. We determined the genome sequences of viruses isolated from two subjects in the Herpevac Trial for Women who experienced primary HSV-2 genital disease and compared them with sequences of viruses isolated from the subsequent fifth or sixth episode of recurrent disease in the same individuals. Each of the HSV-2 genome sequences was initially obtained using next-generation sequencing and completed with Sanger sequencing. Polymorphisms over the entire genomes were mapped, and amino acid variants resulting from nonsynonymous changes were analyzed based on the secondary and tertiary structures of a previously crystallized protein. A phylogenetic reconstruction was used to assess relationships among the four HSV-2 samples, other North American sequences, and reference sequences. Little genetic drift was detected in viruses shed by the same subjects following repeated reactivation events, suggesting strong selective pressure on the viral genome to maintain sequence fidelity during reactivations from its latent state within an individual host. Our results also demonstrate that some primary HSV-2 isolates from North America more closely resemble the HG52 laboratory strain from Scotland than the low-passage-number clinical isolate SD90e from South Africa or laboratory strain 333. Thus, one of the sequences reported here would be a logical choice as a reference strain for inclusion in future studies of North American HSV-2 isolates.IMPORTANCE The extent to which the HSV-2 genome evolves during multiple episodes of reactivation from its latent state within an infected individual is not known. We used next-generation sequencing techniques to determine whole-genome sequences of four viral samples from two subjects in the Herpevac Trial. The sequence of each subject's well-documented primary isolate was compared with the sequence of the isolate from their fifth or sixth episode of recurrent disease. Only 19 genetic polymorphisms unique to the primary or recurrent isolate were identified, 10 in subject A and 9 in subject B. These observations indicate remarkable genetic conservation between primary and recurrent episodes of HSV-2 infection and imply that strong selection pressures exist to maintain the fidelity of the viral genome during repeated reactivations from its latent state. The genome conservation observed also has implications for the potential success of a therapeutic vaccine.

Proteomics show antigen presentation processes in human immune cells after AS03-H5N1 vaccination.

Adjuvants enhance immunity elicited by vaccines through mechanisms that are poorly understood. Using a systems biology approach, we investigated temporal protein expression changes in five primary human immune cell populations: neutrophils, monocytes, natural killer cells, T cells, and B cells after administration of either an Adjuvant System 03 adjuvanted or unadjuvanted split-virus H5N1 influenza vaccine. Monocytes demonstrated the strongest differential signal between vaccine groups. On day 3 post-vaccination, several antigen presentation-related pathways, including MHC class I-mediated antigen processing and presentation, were enriched in monocytes and neutrophils and expression of HLA class I proteins was increased in the Adjuvant System 03 group. We identified several protein families whose proteomic responses predicted seroprotective antibody responses (>1:40 hemagglutination inhibition titer), including inflammation and oxidative stress proteins at day 1 as well as immunoproteasome subunit (PSME1 and PSME2) and HLA class I proteins at day 3 in monocytes. While comparison between temporal proteomic and transcriptomic results showed little overlap overall, enrichment of the MHC class I antigen processing and presentation pathway in monocytes and neutrophils was confirmed by both approaches.

Priming Vaccination With Influenza Virus H5 Hemagglutinin Antigen Significantly Increases the Duration of T cell Responses Induced by a Heterologous H5 Booster Vaccination.

BACKGROUND: Influenza A(H5N1) virus and other avian influenza virus strains represent major pandemic threats. Like all influenza A virus strains, A(H5N1) viruses evolve rapidly. Innovative immunization strategies are needed to induce cross-protective immunity. METHODS: Subjects primed with clade 1 H5 antigen, with or without adjuvant, and H5-naive individuals were boosted with clade 2 H5 antigen. The impact of priming on T cells capable of both proliferation and cytokine production after antigen restimulation was assessed. RESULTS: Subjects previously vaccinated with clade 1 H5 antigen developed significantly enhanced clade 2 H5 cross-reactive T cell responses detectable 6 months after vaccination with clade 2 H5 antigen. Priming dose (15 µg vs 45 or 90 µg) had no effect on magnitude of heterotypic H5 T cell responses. In contrast, age at priming negatively modulated both the magnitude and duration of heterotypic H5 T cell responses. Elderly subjects developed significantly less heterotypic H5 T cell boosting, predominantly for T cells capable of cytokine production. Adjuvant had a positive albeit weaker effect than age. The magnitude of CD4(+) interferon-γ producing T cells correlated with H5 antibody responses. CONCLUSIONS: H5 heterotypic priming prior to onset of an A(H5N1) pandemic may increase magnitude and duration of immunity against a newly drifted pandemic H5 virus.

Protein interaction data curation: the International Molecular Exchange (IMEx) consortium.

The International Molecular Exchange (IMEx) consortium is an international collaboration between major public interaction data providers to share literature-curation efforts and make a nonredundant set of protein interactions available in a single search interface on a common website (http://www.imexconsortium.org/). Common curation rules have been developed, and a central registry is used to manage the selection of articles to enter into the dataset. We discuss the advantages of such a service to the user, our quality-control measures and our data-distribution practices.

A new reference implementation of the PSICQUIC web service.

The Proteomics Standard Initiative Common QUery InterfaCe (PSICQUIC) specification was created by the Human Proteome Organization Proteomics Standards Initiative (HUPO-PSI) to enable computational access to molecular-interaction data resources by means of a standard Web Service and query language. Currently providing >150 million binary interaction evidences from 28 servers globally, the PSICQUIC interface allows the concurrent search of multiple molecular-interaction information resources using a single query. Here, we present an extension of the PSICQUIC specification (version 1.3), which has been released to be compliant with the enhanced standards in molecular interactions. The new release also includes a new reference implementation of the PSICQUIC server available to the data providers. It offers augmented web service capabilities and improves the user experience. PSICQUIC has been running for almost 5 years, with a user base growing from only 4 data providers to 28 (April 2013) allowing access to 151 310 109 binary interactions. The power of this web service is shown in PSICQUIC View web application, an example of how to simultaneously query, browse and download results from the different PSICQUIC servers. This application is free and open to all users with no login requirement (http://www.ebi.ac.uk/Tools/webservices/psicquic/view/main.xhtml).

Safety and immunogenicity of tetanus diphtheria and acellular pertussis (Tdap) immunization during pregnancy in mothers and infants: a randomized clinical trial.

IMPORTANCE: Maternal immunization with tetanus toxoid, reduced diphtheria toxoid, and acellular pertussis (Tdap) vaccine could prevent infant pertussis. OBJECTIVE: To evaluate the safety and immunogenicity of Tdap immunization during pregnancy and its effect on infant responses to diphtheria and tetanus toxoids and acellular pertussis (DTaP) vaccine. DESIGN, SETTING, AND PARTICIPANTS: Phase 1-2, randomized, double-blind, placebo-controlled, clinical trial conducted from 2008 to 2012. Forty-eight pregnant women aged 18 to 45 years received Tdap (n = 33) or placebo (n = 15) at 30 to 32 weeks' gestation, with crossover immunization postpartum. INTERVENTIONS: Tdap vaccination at 30 to 32 weeks' gestation or postpartum. MAIN OUTCOMES AND MEASURES: Primary outcomes were maternal and infant adverse events, pertussis illness, and infant growth and development until age 13 months. Secondary outcomes were antibody concentrations in pregnant women before and 4 weeks after Tdap immunization or placebo, at delivery and 2 months' postpartum, and in infants at birth, at 2 months, and after the third and fourth doses of DTaP. RESULTS: No Tdap-associated serious adverse events occurred in women or infants. Injection site reactions after Tdap immunization were reported in 26 (78.8% [95% CI, 61.1%-91.0%]) and 12 (80% [95% CI, 51.9%-95.7%]) pregnant and postpartum women, respectively (P > .99). Systemic symptoms were reported in 12 (36.4% [ 95% CI, 20.4%-54.9%]) and 11 (73.3% [95% CI, 44.9%-92.2%]) pregnant and postpartum women, respectively (P = .03). Growth and development were similar in both infant groups. No cases of pertussis occurred. Significantly higher concentrations of pertussis antibodies were measured at delivery in women who received Tdap during pregnancy vs postpartum (eg, pertussis toxin antibodies: 51.0 EU/mL [95% CI, 37.1-70.1] and 9.1 EU/mL [95% CI, 4.6-17.8], respectively; P < .001) and in their infants at birth (68.8 EU/mL [95% CI, 52.1-90.8] and 14.0 EU/mL [95% CI, 7.3-26.9], respectively; P < .001) and at age 2 months (20.6 EU/mL [95% CI, 14.4-29.6] and 5.3 EU/mL [95% CI, 3.0-9.4], respectively; P < .001). Antibody responses in infants born to women receiving Tdap during pregnancy were not different following the fourth dose of DTaP. CONCLUSIONS AND RELEVANCE: This preliminary assessment did not find an increased risk of adverse events among women who received Tdap vaccine during pregnancy or their infants. For secondary outcomes, maternal immunization with Tdap resulted in high concentrations of pertussis antibodies in infants during the first 2 months of life and did not substantially alter infant responses to DTaP. Further research is needed to provide definitive evidence of the safety and efficacy of Tdap immunization during pregnancy. TRIAL REGISTRATION: clinicaltrials.gov Identifier: NCT00707148.

Metabolic reconstruction for metagenomic data and its application to the human microbiome.

Microbial communities carry out the majority of the biochemical activity on the planet, and they play integral roles in processes including metabolism and immune homeostasis in the human microbiome. Shotgun sequencing of such communities' metagenomes provides information complementary to organismal abundances from taxonomic markers, but the resulting data typically comprise short reads from hundreds of different organisms and are at best challenging to assemble comparably to single-organism genomes. Here, we describe an alternative approach to infer the functional and metabolic potential of a microbial community metagenome. We determined the gene families and pathways present or absent within a community, as well as their relative abundances, directly from short sequence reads. We validated this methodology using a collection of synthetic metagenomes, recovering the presence and abundance both of large pathways and of small functional modules with high accuracy. We subsequently applied this method, HUMAnN, to the microbial communities of 649 metagenomes drawn from seven primary body sites on 102 individuals as part of the Human Microbiome Project (HMP). This provided a means to compare functional diversity and organismal ecology in the human microbiome, and we determined a core of 24 ubiquitously present modules. Core pathways were often implemented by different enzyme families within different body sites, and 168 functional modules and 196 metabolic pathways varied in metagenomic abundance specifically to one or more niches within the microbiome. These included glycosaminoglycan degradation in the gut, as well as phosphate and amino acid transport linked to host phenotype (vaginal pH) in the posterior fornix. An implementation of our methodology is available at http://huttenhower.sph.harvard.edu/humann. This provides a means to accurately and efficiently characterize microbial metabolic pathways and functional modules directly from high-throughput sequencing reads, enabling the determination of community roles in the HMP cohort and in future metagenomic studies.

Distinct gene expression signatures comparing latent tuberculosis infection with different routes of Bacillus Calmette-Guérin vaccination.

Tuberculosis remains an international health threat partly because of limited protection from pulmonary tuberculosis provided by standard intradermal vaccination with Bacillus of Calmette and Guérin (BCG); this may reflect the inability of intradermal vaccination to optimally induce pulmonary immunity. In contrast, respiratory Mycobacterium tuberculosis infection usually results in the immune-mediated bacillary containment of latent tuberculosis infection (LTBI). Here we present RNA-Seq-based assessments of systemic and pulmonary immune cells from LTBI participants and recipients of intradermal and oral BCG. LTBI individuals uniquely display ongoing immune activation and robust CD4 T cell recall responses in blood and lung. Intradermal BCG is associated with robust systemic immunity but only limited pulmonary immunity. Conversely, oral BCG induces limited systemic immunity but distinct pulmonary responses including enhanced inflammasome activation potentially associated with mucosal-associated invariant T cells. Further, IL-9 is identified as a component of systemic immunity in LTBI and intradermal BCG, and pulmonary immunity following oral BCG.

Corrigendum: The Vacc-SeqQC project: Benchmarking RNA-Seq for clinical vaccine studies.

[This corrects the article DOI: 10.3389/fimmu.2022.1093242.].

Utility of targeted gene sequencing to differentiate myeloid malignancies from other cytopenic conditions.

The National Heart, Lung, and Blood Institute-funded National MDS Natural History Study (NCT02775383) is a prospective cohort study enrolling patients with cytopenia with suspected myelodysplastic syndromes (MDS) to evaluate factors associated with disease. Here, we sequenced 53 genes in bone marrow samples harvested from 1298 patients diagnosed with myeloid malignancy, including MDS and non-MDS myeloid malignancy or alternative marrow conditions with cytopenia based on concordance between independent histopathologic reviews (local, centralized, and tertiary to adjudicate disagreements when needed). We developed a novel 2-stage diagnostic classifier based on mutational profiles in 18 of 53 sequenced genes that were sufficient to best predict a diagnosis of myeloid malignancy and among those with a predicted myeloid malignancy, predict whether they had MDS. The classifier achieved a positive predictive value (PPV) of 0.84 and negative predictive value (NPV) of 0.8 with an area under the receiver operating characteristic curve (AUROC) of 0.85 when classifying patients as having myeloid vs no myeloid malignancy based on variant allele frequencies (VAFs) in 17 genes and a PPV of 0.71 and NPV of 0.64 with an AUROC of 0.73 when classifying patients as having MDS vs non-MDS malignancy based on VAFs in 10 genes. We next assessed how this approach could complement histopathology to improve diagnostic accuracy. For 99 of 139 (71%) patients (PPV of 0.83 and NPV of 0.65) with local and centralized histopathologic disagreement in myeloid vs no myeloid malignancy, the classifier-predicted diagnosis agreed with the tertiary pathology review (considered the internal gold standard).

The Protein Naming Utility: a rules database for protein nomenclature.

Generation of syntactically correct and unambiguous names for proteins is a challenging, yet vital task for functional annotation processes. Proteins are often named based on homology to known proteins, many of which have problematic names. To address the need to generate high-quality protein names, and capture our significant experience correcting protein names manually, we have developed the Protein Naming Utility (PNU, http://www.jcvi.org/pn-utility). The PNU is a web-based database for storing and applying naming rules to identify and correct syntactically incorrect protein names, or to replace synonyms with their preferred name. The PNU allows users to generate and manage collections of naming rules, optionally building upon the growing body of rules generated at the J. Craig Venter Institute (JCVI). Since communities often enforce disparate conventions for naming proteins, the PNU supports grouping rules into user-managed collections. Users can check their protein names against a selected PNU rule collection, generating both statistics and corrected names. The PNU can also be used to correct GenBank table files prior to submission to GenBank. Currently, the database features 3080 manual rules that have been entered by JCVI Bioinformatics Analysts as well as 7458 automatically imported names.

Failure to Detect Mutations in U2AF1 due to Changes in the GRCh38 Reference Sequence.

The U2AF1 gene is a core part of mRNA splicing machinery and frequently contains somatic mutations that contribute to oncogenesis in myelodysplastic syndrome, acute myeloid leukemia, and other cancers. A change introduced in the GRCh38 version of the human reference build prevents detection of mutations in this gene, and others, by variant calling pipelines. This study describes the problem in detail and shows that a modified GRCh38 reference build with unchanged coordinates can be used to ameliorate the issue.

The antibody landscapes following AS03 and MF59 adjuvanted H5N1 vaccination.

Current seasonal and pre-pandemic influenza vaccines induce short-lived predominantly strain-specific and limited heterosubtypic responses. To better understand how vaccine adjuvants AS03 and MF59 may provide improved antibody responses to vaccination, we interrogated serum from subjects who received 2 doses of inactivated monovalent influenza A/Indonesia/05/2005 vaccine with or without AS03 or MF59 using hemagglutinin (HA) microarrays (NCT01317758 and NCT01317745). The arrays were designed to reflect both full-length and globular head HA derived from 17 influenza A subtypes (H1 to H16 and H18) and influenza B strains. We observed significantly increased strain-specific and broad homo- and heterosubtypic antibody responses with both AS03 and MF59 adjuvanted vaccination with AS03 achieving a higher titer and breadth of IgG responses relative to MF59. The adjuvanted vaccine was also associated with the elicitation of stalk-directed antibody. We established good correlation of the array antibody responses to H5 antigens with standard HA inhibition and microneutralization titers.

The Vacc-SeqQC project: Benchmarking RNA-Seq for clinical vaccine studies.

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.