Whole Transcriptome Profiling Identifies CD93 and Other Plasma Cell Survival Factor Genes Associated with Measles-Specific Antibody Response after Vaccination

Iana H. Haralambieva; Michael T. Zimmermann; Inna G. Ovsyannikova; Diane E. Grill; Ann L. Oberg; Richard B. Kennedy; Gregory A. Poland

doi:10.1371/journal.pone.0160970

Abstract

Background

There are insufficient system-wide transcriptomic (or other) data that help explain the observed inter-individual variability in antibody titers after measles vaccination in otherwise healthy individuals.

Methods

We performed a transcriptome(mRNA-Seq)-profiling study after in vitro viral stimulation of PBMCs from 30 measles vaccine recipients, selected from a cohort of 764 schoolchildren, based on the highest and lowest antibody titers. We used regression and network biology modeling to define markers associated with neutralizing antibody response.

Results

We identified 39 differentially expressed genes that demonstrate significant differences between the high and low antibody responder groups (p-value≤0.0002, q-value≤0.092), including the top gene CD93 (p<1.0E^-13, q<1.0E^-09), encoding a receptor required for antigen-driven B-cell differentiation, maintenance of immunoglobulin production and preservation of plasma cells in the bone marrow. Network biology modeling highlighted plasma cell survival (CD93, IL6, CXCL12), chemokine/cytokine activity and cell-cell communication/adhesion/migration as biological processes associated with the observed differential response in the two responder groups.

Conclusion

We identified genes and pathways that explain in part, and are associated with, neutralizing antibody titers after measles vaccination. This new knowledge could assist in the identification of biomarkers and predictive signatures of protective immunity that may be useful in the design of new vaccine candidates and in clinical studies.

Citation: Haralambieva IH, Zimmermann MT, Ovsyannikova IG, Grill DE, Oberg AL, Kennedy RB, et al. (2016) Whole Transcriptome Profiling Identifies CD93 and Other Plasma Cell Survival Factor Genes Associated with Measles-Specific Antibody Response after Vaccination. PLoS ONE 11(8): e0160970. https://doi.org/10.1371/journal.pone.0160970

Editor: John S. Tregoning, Imperial College London, UNITED KINGDOM

Received: May 19, 2016; Accepted: July 27, 2016; Published: August 16, 2016

Copyright: © 2016 Haralambieva et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: All data are available from ImmPort (https://immport.niaid.nih.gov/immportWeb/clinical/study/displayStudyDetails.do?itemList=SDY839; study identifier SDY839).

Funding: This work was supported by NIH grants R37AI048793-11, R01AI033144 and was made possible by the Rochester Epidemiology Project (Grant Number R01 AG034676 from the National Institute of Aging). Its contents are solely the responsibility of the authors and do not necessarily represent the official view of NCRR or NIH. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: Dr. Poland is the chair of a Safety Evaluation Committee for novel investigational vaccine trials being conducted by Merck Research Laboratories. Dr. Poland offers consultative advice on vaccine development to Merck & Co. Inc., CSL Biotherapies, Avianax, Dynavax, Novartis Vaccines and Therapeutics, Emergent Biosolutions, Adjuvance, Microdermis, Seqirus, NewLink, Protein Sciences, GSK Vaccines, and Sanofi Pasteur. Dr. Poland holds two patents related to vaccinia and measles peptide research. Dr. Kennedy has received funding from Merck Research Laboratories to study waning immunity to mumps vaccine. These activities have been reviewed by the Mayo Clinic Conflict of Interest Review Board and are conducted in compliance with Mayo Clinic Conflict of Interest policies. This research has been reviewed by the Mayo Clinic Conflict of Interest Review Board and was conducted in compliance with Mayo Clinic Conflict of Interest policies. This does not alter the authors' adherence to PLOS ONE policies on sharing data and materials.

Introduction

Despite the widespread use of measles vaccines (especially the two-dose immunization schedule introduced in 1989) and the dramatic decrease in the occurrence of measles in the U.S., annual outbreaks involving hundreds of cases continue to occur.[1] In 2014 alone, the U.S. experienced 23 measles outbreaks with 644 cases (the highest number of cases since measles elimination was declared a priority in 2000).[2] In 2015, 189 people from 24 states and the District of Columbia were reported to have measles.[2] Recent outbreaks among highly vaccinated populations in North America and in Europe are an indication that, even with the two-dose schedule, vaccine failure still accounts for a small to sizeable percentage (1%-24%) of measles cases.[3,4] A significant portion of these cases involve recipients of two doses of measles-containing vaccine.[1,3]

There are insufficient data on the immunologic, genetic and genomic basis of inter-individual immune response variations after measles vaccination to help uncover the molecular mechanisms underlying vaccine failure. Comprehensive omics-level information on measles vaccine immune responses and waning immunity has not been previously generated. To fill these knowledge gaps, we developed a new vaccinomics paradigm that utilizes unbiased, high-dimensional omics data and advanced statistical/bioinformatics approaches for deep immune profiling of viral vaccine responses.[5,6,7] The ultimate goal is to apply these approaches to identify signatures (containing both individual variables and their interactions) that can discriminate among measles vaccine immune phenotypes, serve as biomarkers of protective immunity, and inform the development of novel vaccine candidates.

In this study, we assessed gene expression in 30 measles-mumps-rubella (MMR) vaccine recipients following two vaccine doses using a whole transcriptome sequencing approach (mRNA-Seq), and evaluated how high or low neutralizing antibody response impacts and/or is associated with observed transcriptional changes after in vitro measles virus stimulation. Our study was designed to better understand the genetic factors, mechanisms and pathways underlying the biological spectrum of immune responses to measles vaccination.

Methods

The methods described herein are similar or identical to those we have previously published. [8,9,10,11,12,13,14]

2.1. Study subjects

A total of 30 subjects were selected for an mRNA-Seq transcriptome-profiling study based on their plaque reduction microneutralization assay (PRMN) antibody titers (15 highest and 15 lowest antibody responders) from a combined cohort of 764 randomly selected healthy schoolchildren and young adults (age 11–22 years) from Olmsted County, MN, who received two doses of MMR-II vaccine (Merck) that contained the Edmonston strain of measles virus (TCID₅₀ ≥1,000). Details on the recruitment, demographic and immune characteristics of the study population have been previously described.[10,11,12] The Institutional Review Board (IRB) of the Mayo Clinic approved the study, and written, informed consent was obtained from subjects’ parents/guardians, as well as written assent from age-appropriate subjects at the time of enrollment.

2.2. Plaque reduction microneutralization assay (PRMN)

Anti-measles neutralizing antibody titers were quantified using a high throughput fluorescence-based PRMN, as previously published.[10,11,12] Estimates of 50% end point titer (Neutralizing Dose, ND50) were calculated using Karber’s formula and ND50 values were converted to mIU/ml values using the 3rd WHO anti-measles antibody standard (NIBSC code no. 97/648), as previously published. [12] The coefficient of variation for this assay in our laboratory, based on the third WHO standard runs was 5.7%.[12]

2.3. mRNA-Seq transcriptome profiling

The sequencing methods are similar to those previously published.[13,14] In brief, subjects’ cryopreserved peripheral blood mononuclear cells (PBMCs) were thawed and stimulated with live Edmonston measles virus at a multiplicity of infection of 0.5 for 24 hours (for each subject, an aliquot of the cells was left unstimulated). Cells were stabilized with RNAprotect cell reagent (Qiagen) and total RNA was extracted using RNeasy Plus mini kit (Qiagen). Quality and quantity of RNA was determined by Nanodrop spectrophotometry (Thermo Fisher Scientific). Poly-A RNA was isolated using magnetic beads coated with olido-dT, then reverse transcribed after fragmentation into cDNA, and combined with Illumina adaptor sequences. Libraries were generated using Illumina’s mRNA TruSeq (v1) kit. After validation, cDNA libraries were sequenced (paired end sequencing) on an Illumina HiSeq 2000 (Illumina;San Diego, CA) with Illumina’s TruSeq Cluster kit (v3-cBot-HS) and 51 Cycle Illumina TruSeq SBS Sequencing Kit (v3). The sequencing reads were aligned to the human genome build 37.1 using TopHat (1.3.3) and Bowtie (0.12.7). Gene counts were performed using HTSeq (0.5.3p3), while BEDTools software (2.7.1) was used to map normalized read count to individual exons.[15,16,17]

2.4. Statistical analysis

The analysis consisted of complementary gene-to-biology (an inductive, data-driven approach performed at the gene level) and biology-to-gene (a deductive knowledge-driven approach performed at the geneset level) analytical strategies, as previously described.[18]

The samples (n = 60, 30 measles virus-stimulated and 30 unstimulated samples) were randomized to ensure balance of important characteristics over assay run order: immune response group (high or low antibody responder), stimulation status (virus stimulated or unstimulated sample), and sex (male or female) for the cell culture setup, preparation of libraries, and sequencing (flow cell and lane run).

Quality control methods used were similar to those reported in our previous studies. [13,14] Conditional Quantile Normalization (CQN) was used to normalize the mRNA-sequence data adjusting for the differences in library size, GC content and gene length.[19,20] Differential expression analysis comparing statistical “interaction” of the high vs low responders’ response to viral stimulation; viral stimulation relative to the unstimulated samples (regardless of neutralizing antibody status) was performed using generalized estimating equations.[21] Specifically, the raw counts were modeled in a generalized linear negative binomial model [22], utilizing the offset from the CQN normalization, estimating the variance for each gene with the edgeR tagwise dispersion option, and estimating the within subject correlation structure to provide a robust estimate of the variance.[23] Results for the differential expression analysis are reported as the interaction comparing the fold change (FC) for the difference in high antibody responders’ response to viral stimulation, relative to the response of low antibody responders’ to viral stimulation, and the response to viral stimulation for the high and low responder subjects combined. All reported fold-change values (including the interaction fold change) for gene expression are for virus-stimulated PBMCs relative to the unstimulated PBMCs, and thus take into account resting differences in gene expression. False discovery rates (FDR) were calculated using standard methods and are reported as q-values.[24] Pathway analysis was conducted using Ingenuity® software, (Ingenuity® Systems; Redwood City, CA). Geneset analysis was done using the gamma method with 1,000 permutations of the high and low antibody response of the subjects permuting the stimulated status within each subject.[25]

2.5. Network biology

To comprehensively identify high-confidence interactions, we combined multiple resources: HPRD [26]; CCSB [27]; the 7.8% of STRING [28] with highest confidence (score ≥ 70%); and the 12.9% of HumanNet [29] with highest confidence (score ≥ 2). Networks were visualized using Cytoscape [30] version 3.2.1 and layouts refined using AllegroLayout v.2.2.1.[31] Pathway enrichment was performed using 138 pathway definitions downloaded from MSigDB’s index of KEGG canonical pathways [32,33] after removing those that are disease-associated (e.g., type-II diabetes and leishmania infection), or derived from other pathways (e.g., pathways in cancer). The biologic functions of genesets were assessed using GO term enrichment [34] and evaluated using hypergeometric tests. Genes considered in our network biology approach are all those exhibiting stimulation and interaction fold changes (|log₂(FC)| of ≥ 0.5)and an interaction p-value ≤ 0.01 (n = 115), or a stimulation p-value of ≤ 0.001 (n = 207 additional genes). One-hundred forty (43%) of these 322 genes share connections within the high-confidence network used. We also generated an inclusive network where the score-based filters were not employed, and where 282 (88%) genes share connections (data not shown).

Results

3.1. Characteristics of the study subjects

The enrolled study subjects were primarily white (67% for the high-responders group and 73% for the low-responders group); African-Americans comprised 27% and 20% of the high- and low-responders group, respectively. Sixty percent of the participants were males and forty percent were females. The median age of subjects at study enrollment was 16 (IQR 12, 17) and 15 (IQR 13, 17) years (for the high and low responders group, respectively); the median age at first measles immunization was 16 (IQR 15, 32) and 15 (IQR 15, 65) months; the median age at second immunization was 6 (IQR 4, 12) and 8 (IQR 5, 11) years; and the median time since second immunization to enrollment was 7.1 (5.7, 9.5) and 7.0 (IQR 4.6, 8.7) years (for the high and low responders group, respectively). Thus, it is important to note that gene expression in PBMCs and all immune measures (neutralizing antibody and cytokine measures after in vitro viral stimulation) were quantified in samples approximately seven years after the last/second measles (MMR) vaccination. None of the above mentioned demographic variable differences reached statistical significance (data not shown). The study subjects were selected for mRNA-Seq profiling based on their neutralizing antibody titer. The median antibody titer for the high responders group was 5,188 mIU/mL (which suggests protection against infection/disease), while the titer for the low responders group was 88 mIU/mL (which suggests lack of protection against measles). All other immune response variables were not statistically significant between the two groups and are shown in S1 Table.

3.2. Gene expression after in vitro stimulation with measles virus in high and low antibody responders (interaction analysis)

Our overall analysis (differential gene expression in response to viral stimulation in all samples/subjects) identified 1,761 significantly expressed genes (FC>2 or FC<0.5, FDR<0.004), including chemokine and chemokine receptor genes, cytokine and cytokine receptor genes, and genes encoding innate receptors, antiviral proteins and HLA (S2 Table).

Importantly, we identified 39 differentially expressed genes (FDR<0.1, Table 1) that demonstrated significant gene expression differences between the two antibody responder groups to measles vaccination (p-value range 1.0E^-13 to 0.0002, q-values range 1.0E^-09 to 0.092, Table 1). An additional 11 differentially expressed genes demonstrated suggestive gene expression differences (as per our gene-to-biology approach) between the two study groups (FDR<0.15, Table 1).

Download:

Table 1. Differential response to measles virus stimulation in high vs. low antibody vaccine responders (genes with q-value <0.15).

https://doi.org/10.1371/journal.pone.0160970.t001

3.3. Pathway and geneset analysis

Among the Ingenuity pathways that were significantly enriched in high vs. low antibody responders (enrichment p<1.0E^-03), the top three pathways (consisting of highly overlapping genes) are related to immune adhesion and function, chemotaxis/inflammatory response, and cytokine regulation (Table 2). To augment Ingenuity-based analysis, we used other independent annotation sources, including GO and a filtered set of canonical pathways (see Methods). The top results from each comparison are listed in Table 3, with most sources highlighting chemokine/cytokine activity, inflammatory response, and cell-cell communication/adhesion as biological processes significantly involved in differential response to viral stimulation in high vs. low measles vaccine antibody responders.

Download:

Table 2. Top Ingenuity enriched pathways of differentially expressed genes to measles virus stimulation in high vs. low antibody vaccine responders.

https://doi.org/10.1371/journal.pone.0160970.t002

Download:

Table 3. Pathway enrichment of differentially expressed genes upon measles virus stimulation using independent annotation sources.

https://doi.org/10.1371/journal.pone.0160970.t003

We also performed a geneset analysis using genesets/modules downloaded from MSigDB to identify associations of groups of genes (genesets) with neutralizing antibody response following vaccination. [35] The results identified 112 significant genesets (p < 0.05) with different expression (response to viral stimulation) in high vs. low antibody responders that comprise genes/pathways integral to innate and adaptive immune response, cell adhesion, metabolism and cell signaling. The top 7 significant genesets with p<0.003 are listed in S3 Table.

In search of consistency between our gene-to-biology and biology-to-gene approach results, we compared the top 69 genes from the gene-to-biology results (with an interaction p-value < = 0.001) with the genes from 1,213 genesets (with p-value ≤ 0.002) resulting from our biology-to-gene approach. The intersection (overlap) between the results from the two approaches resulted in six common genes of high interest (i.e., CD93, CD33, IL10, VEGFA, VNN1, CADM1).

3.4. Network biology

To display molecular interactions, we performed simultaneous gene mapping/visualization of: 1.) genes with significant induction/suppression upon viral stimulation (overall response to viral stimulation in all study subjects, regardless of their antibody response); and 2.) genes that show significant differences in gene expression between the two antibody responder groups in response to viral stimulation (i.e., statistical interaction). The result of gene mapping using only high-confidence molecular interactions is shown in Fig 1.

Download:

Fig 1. Legend.

Network biology of gene expression in measles vaccine high and low antibody responders.Simultaneous visualization of overall stimulation response (gene expression in response to viral stimulation in all subjects) and differential responses (in high vs. low antibody responders) within the context of a high-confidence biologic network. Each gene is represented as two circles, the outer colored by the virus-induced fold change in all subjects (FC, presented as log₂ Stim FC) and the inner by the statistical interaction between high and low responders (i.e., difference in stimulated FC between the high and low antibody responder groups, presented as log₂ FC between high and low responders).

https://doi.org/10.1371/journal.pone.0160970.g001

Discussion

We and others have pointed out the value of using a systems-level vaccinomics approach leveraging unbiased large-scale gene expression profiling (e.g., mRNA-Seq) to delineate molecular signatures of efficient immune response after vaccination. Such approaches have provided important insights into immunity to influenza [36], yellow fever [37], smallpox [14], rubella [13], polysaccharide pneumococcal [38], and other vaccines. Our motivation for performing this study was that despite the continuing measles outbreaks in populations that have received two doses of measles vaccine, comprehensive genome-wide information on gene expression changes after measles vaccination and/or infection is limited.

Our transcriptome profiling study delineates significant gene expression differences that discriminate between high and low neutralizing antibody responders (to measles vaccination) in PBMCs stimulated in vitro with measles virus. The identified genes include cytokine genes, cell signaling genes, transcriptional targets of NF-kappa B, C-type lectin transmembrane receptors, and other immune function-related genes (Table 1). The top gene found (p-value<1.0E^-13, Table 1) is CD93, which is a receptor required for antigen-driven B-cell differentiation for maintenance of immunoglobulin production and preservation of plasma cells in the bone marrow. [39] CD93 was not only the most significant gene, but it was also among the top three genes with the highest fold change (FC interaction = 6.55) for differential gene expression (upon viral stimulation) between high and low antibody responders. The other two genes with high fold change were the phosphotyrosine interaction domain containing 1 gene (PID1, FC interaction = 12.78, p-value = 3.34E^-11) and vanin 1(VNN1, FC interaction = 8.09, p-value = 6.01E^-06). CD93 (identified also as a complement protein 1q receptor—C1qRp) encodes a C-type lectin transmembrane receptor, expressed primarily during the early stages of B cell differentiation (although can be expressed by different cell lineages, including activated macrophages). CD93 expression is re-induced during plasma cell differentiation and the long-lived plasma cells demonstrate high expression level of this protein, which is crucial for their maintenance in the bone marrow niche.[39,40] Studies in a CD93-deficient mouse model revealed the key role CD93 has in terms of the survival of antigen-specific long-lived plasma cells and the persistence of antigen-specific antibody response.[39]

Interestingly, the expression of CD93 was downregulated upon measles virus stimulation in our study (in both antibody responder groups); however, this downregulation should be interpreted only in the context of gene expression in PBMCs (and not in specific cell populations). A study of gammaherpesvirus infection in mice has demonstrated that the majority of normal B cells were activated in the course infection with downregulation of CD93. [41] Different B cells may undergo different degrees of CD93 expression (including loss of CD93); however, the degree of downregulation or upregulation in different individuals and cell subtypes may underlie differences in cell homing, adhesion and cell survival. The downregulation (FC log2) in the high antibody responder group was limited compared to the low-responder group (i.e., CD93 expression remained at a higher expression level after viral stimulation in the high antibody responders compared to the low antibody responders), which is in agreement with the association of higher CD93 expression with long-term antibody persistence after vaccination. Furthermore, CD93 cytoplasmic domain has been shown to interact with moesin for cytoskeleton remodeling, and thus this protein may have a role beyond maintenance of long-lived plasma cells in the processes of cell adhesion, cell migration and phagocytosis.[42] Another interesting finding is the identification of IL24 as the second-ranked top gene (p = 5.8E^-13, Table 1), demonstrating transcriptional differences in high vs. low neutralizing antibody responders. IL-24 is a novel cytokine with a recently described role in regulating antigen-driven driven B cell differentiation and plasma cell vs. memory cell commitment in germinal centers; therefore, its expression can be directly related to antibody production.[43]

Different annotation sources emphasize several major enrichment themes among the genes passing statistical criteria: cytokine/chemokine activity/extracellular messengers, inflammatory response and cell adhesion/chemotaxis/migration. Specifically, Gene Ontology (GO) annotates many genes as cytokine and chemokine signaling localized in the extracellular space. MSigDB genesets highlight the importance of inflammatory response and the matrisome (the “ome” of extracellular matrix proteins). The filtered set of KEGG canonical pathways likewise highlights both cytokine and chemokine signaling pathways and cell adhesion/migration processes, which is similar to Ingenuity pathways results.

Our network biology analysis mapped genes with significant statistical associations with either overall (in all study subjects) virus-induced gene induction/suppression, or significant genes that discriminate the transcriptional patterns between high and low antibody responder groups to reveal their molecular interaction and biological relationships (Fig 1). Several modules are evident: one predominant module is centered on IL6, VEGFA, LEP, and the two tyrosine-protein kinases SRC and EPHA2, and includes genes with a distinct known relation to innate/adaptive immunity. Interleukin 6 is produced by various cell types, including macrophages and T cells to initiate and promote immune function, and is also a critical survival factor for the resident bone-marrow plasma cells.[44] Leptin is a known protein important for regulation of metabolism and body weight, which plays an important role in the modulation of inflammatory/innate and adaptive immune responses following infection or vaccination.[45] A second predominant module is linked to chemotaxis/inflammation and is centered on CCR2 (encodes a receptor for monocyte chemoattractant protein-1, involved in monocyte chemotaxis/infiltration in inflammation) and several other key chemoattractants (or their receptors) for B and T cells, monocytes, dendritic cells and granulocytes (i.e., CXCL3, CXCL6, CXCL12, CXCL13, CCL7, CCL19, XCR1). A third smaller module is centered on CCL8, which is also a chemotactic protein with a role in inflammatory response and migration/activation of monocytes, T cells and NK cells. Previous studies have identified similar or identical genes (e.g., IL-6, CXCL3, CCL7, CCL19, CCR2), pathways and biological processes, involved in the differential response of human PBMCs and dendritic cells to measles infection or viral stimulation.[46,47,48] Intriguing is the facts that several of the highlighted genes/proteins (i.e., IL6, CXCL12, and CD93) are known potent survival signals for the maintenance of plasma cells in the bone marrow and thus their expression may have a direct impact on the regulation of long-term humoral immune responses and immune memory.[44]

The overlap between our data-driven gene-to-biology analysis approach and our biology-to-gene knowledge-driven geneset approach indicates and supports the major involvement of six genes in the differential transcriptional responses observed in the two extreme antibody response groups. These include CD93 (discussed above); CD33 and CADM1, involved in cell adhesion, regulation of proliferation and signal transduction; VEGFA, an important cell growth factor; VNN1, involved in T cell migration and the response to oxidative stress; and IL10, encoding a cytokine with pleiotropic immunosuppressive and immunomodulatory properties, which affects antigen presentation and both T and B cell function (including antibody production).[49,50]

To the best of our knowledge, this is the first study examining and characterizing genome-wide transcriptional responses (using unbiased mRNA-Seq technology) to measles virus in PBMCs of measles-vaccine recipients, which represents the extremes of the neutralizing antibody response after vaccination. Current state-of-the-art paradigms advocate for performing a comprehensive series of unbiased systems biology studies in order to thoroughly define immune/molecular signatures of immune responses. [5,6,7] In this work, we identified immune signatures responsible for (or indicative of) the maintenance of antigen-specific bone marrow plasma cells and development/maintenance of adaptive B cell immune responses to measles vaccination. Due to sample availability, our study design did not include a longitudinal sampling (before and after vaccination), nor the integration of multiple omics-level data into predictive models of vaccine response). At this early stage of knowledge, the current results and identified genes/pathways provide the framework for future systems-level studies in order to better understand measles vaccine-induced immunity.

Two major strengths of our study are the collection of unbiased high-quality mRNA-Seq data from samples of a well-characterized cohort of measles vaccine recipients (after two doses of MMR vaccine and no wild type measles virus exposure), and our statistical (gene-to-biology and biology-to-gene) and bioinformatics approaches, which enabled us to identify specific genes (e.g., CD93, IL24, IL6, CXCL12), pathways and processes (e.g., cell adhesion and migration, cytokine/chemokine activity and regulation, inflammatory response) and network components that are biologically relevant. The major limitation of the study is the possibility of false-positive associations of transcriptional changes with antibody titers. To guide interpretation and help evaluate the level of evidence, we report both p-values and q-values (for the per-variable analysis), and use a complementary knowledge-driven biology-to-gene approach (as part of our analysis) to control for false-discovery rate. In addition, the interpretability of our findings should be the context of in vitro PBMCs re-stimulation post-vaccination and may not necessarily reflect gene expression is different cell populations in vivo after vaccination. Our findings are novel, however, as this is the first reported mRNA-Seq gene expression study assessing gene expression in cells from high and low measles vaccine responders. We believe the differential ability of the cells/PBMCs (from high and low responders) to respond to viral stimulation may reflect differences in cell homing, adhesion, and cell survival of specific cell subtypes (e.g., plasma cells) and, indirectly, antibody titer. However, validation of the identified genes/pathways in a larger measles vaccine systems biology study (including assessment of gene expression in specific cell subtypes) is necessary and is underway in our laboratory.

In conclusion, we report the first comprehensive transcriptome-level characterization (mRNA-Seq data) of responses to measles virus stimulation in high and low antibody responders to measles vaccination, using vaccinomics, statistical and network biology approaches to define plausible regulators (genes/pathways/networks) that drive and/or underlie the observed differences in neutralizing antibody titers after measles vaccination. Studies such as ours enable us to explain how specific markers of adaptive (or innate) immune response are influenced by transcriptional (or other) changes dues to vaccination, and to develop a panel of biomarkers or models that predict and explain the immune response to measles vaccine.

Supporting Information

S1 Table. Immune response variables of the study subjects.

https://doi.org/10.1371/journal.pone.0160970.s001

(DOCX)

S2 Table. Significant genes based on overall response to measles virus stimulation in all vaccine recipients (1,761 significant genes with FC>2 and FC <0.5).

https://doi.org/10.1371/journal.pone.0160970.s002

(DOCX)

S3 Table. Significant genesets, in high vs low antibody vaccine responders (response to viral stimulation)

https://doi.org/10.1371/journal.pone.0160970.s003

(DOCX)

Acknowledgments

We thank the Mayo Clinic Vaccine Research Group staff and the study participants. We thank Caroline L. Vitse for her editorial assistance with this manuscript and Nathaniel D. Warner for his programming assistance and contribution to statistical analysis. This work was supported by NIH grants R37AI048793-11, R01AI033144 and was made possible by the Rochester Epidemiology Project (Grant Number R01 AG034676 from the National Institute of Aging). Its contents are solely the responsibility of the authors and do not necessarily represent the official view of NCRR or NIH.

Author Contributions

Conceptualization: IHH IGO ALO RBK GAP.
Data curation: IHH MTZ DEG IGO RBK.
Formal analysis: MTZ DEG ALO.
Funding acquisition: GAP.
Methodology: IHH MTZ IGO DEG ALO RBK.
Project administration: GAP.
Resources: RBK GAP.
Software: MTZ DEG ALO.
Supervision: ALO RBK GAP.
Visualization: IHH MTZ.
Writing - original draft: IHH MTZ DEG IGO ALO RBK.
Writing - review & editing: IHH MTZ DEG IGO ALO RBK GAP.

References

1. Whitaker JA, Poland GA. Measles and mumps outbreaks in the United States: Think globally, vaccinate locally. Vaccine 2014; 32:4703–4. pmid:24992719
- View Article
- PubMed/NCBI
- Google Scholar
2. Centers for Disease Control and Prevention. Measles Cases and Outbreaks. Available at: http://www.cdc.gov/measles/cases-outbreaks.html. Date accessed: February 12, 2015.
3. Parker Fiebelkorn A, Redd SB, Gallagher K, et al. Measles in the United States during the postelimination era. J Infect Dis 2010; 202:1520–1528. pmid:20929352
- View Article
- PubMed/NCBI
- Google Scholar
4. Increased transmission and outbreaks of measles—European Region, 2011. MMWR 2011; 60:1605–10. pmid:22129994
- View Article
- PubMed/NCBI
- Google Scholar
5. Poland GA, Ovsyannikova IG, Kennedy RB, Haralambieva IH, Jacobson RM. Vaccinomics and a new paradigm for the development of preventive vaccines against viral infections. Omics 2011; 15:625–636. pmid:21732819
- View Article
- PubMed/NCBI
- Google Scholar
6. Haralambieva IH, Poland GA. Vaccinomics, predictive vaccinology and the future of vaccine development. Future Microbiol 2010; 5:1757–1760. pmid:21155658
- View Article
- PubMed/NCBI
- Google Scholar
7. Lambert ND, Ovsyannikova IG, Pankratz VS, Jacobson RM, Poland GA. Understanding the immune response to seasonal influenza vaccination in older adults: a systems biology approach. Exp Rev Vaccines 2012; 11:985–994.
- View Article
- Google Scholar
8. Kennedy RB, Ovsyannikova IG, Haralambieva IH, O'Byrne MM, Jacobson RM, Pankratz VS, et al. Multigenic control of measles vaccine immunity mediated by polymorphisms in measles receptor, innate pathway, and cytokine genes. Vaccine 2012; 30:2159–2167. pmid:22265947
- View Article
- PubMed/NCBI
- Google Scholar
9. Ovsyannikova IG, Haralambieva IH, Vierkant RA, O'Byrne MM, Jacobson RM, Poland GA. The association of CD46, SLAM, and CD209 cellular receptor gene SNPs with variations in measles vaccine-induced immune responses—a replication study and examination of novel polymorphisms. Human Hered 2011; 72:206–223.
- View Article
- Google Scholar
10. Haralambieva IH, Ovsyannikova IG, Umlauf BJ, Vierkant RA, Shane Pankratz V, Jacobson RM, et al. Genetic polymorphisms in host antiviral genes: associations with humoral and cellular immunity to measles vaccine. Vaccine 2011; 29:8988–8997. pmid:21939710
- View Article
- PubMed/NCBI
- Google Scholar
11. Haralambieva IH, Ovsyannikova IG, Kennedy RB, Vierkant RA, Pankratz VS, Jacobson RM, et al. Associations between single nucleotide polymorphisms and haplotypes in cytokine and cytokine receptor genes and immunity to measles vaccination. Vaccine 2011; 29:7883–7895. pmid:21875636
- View Article
- PubMed/NCBI
- Google Scholar
12. Haralambieva IH, Ovsyannikova IG, O'Byrne M, Pankratz VS, Jacobson RM, Poland GA. A large observational study to concurrently assess persistence of measles specific B-cell and T-cell immunity in individuals following two doses of MMR vaccine. Vaccine 2011; 29:4485–4491. pmid:21539880
- View Article
- PubMed/NCBI
- Google Scholar
13. Haralambieva IH, Oberg AL, Ovsyannikova IG, Kennedy RB, Grill DE, Middha S, et al. Genome-wide characterization of transcriptional patterns in high and low antibody responders to rubella vaccination. PLos ONE 2013; 8:e62149. pmid:23658707
- View Article
- PubMed/NCBI
- Google Scholar
14. Kennedy RB, Oberg AL, Ovsyannikova IG, Haralambieva IH, Grill DE, Poland GA. Transcriptomic profiles of high and low antibody responders to smallpox vaccine. Genes Immunity 2013; 14:277–285. pmid:23594957
- View Article
- PubMed/NCBI
- Google Scholar
15. Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics 2010; 26:841–2. pmid:20110278
- View Article
- PubMed/NCBI
- Google Scholar
16. Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics 2009; 25:1105–11. pmid:19289445
- View Article
- PubMed/NCBI
- Google Scholar
17. Anders S, Pyl PT, Huber W. HTSeq—a Python framework to work with high-throughput sequencing data. Bioinformatics 2015; 31:166–9. pmid:25260700
- View Article
- PubMed/NCBI
- Google Scholar
18. Oberg AL, Kennedy RB, Li P, Ovsyannikova IG, Poland GA. Systems biology approaches to new vaccine development. Curr Opin Immunol 2011; 23:436–443. pmid:21570272
- View Article
- PubMed/NCBI
- Google Scholar
19. Hansen KD, Irizarry RA, Wu Z. Removing technical variability in RNA-seq data using conditional quantile normalization. Biostatistics 2012; 13:204–216. pmid:22285995
- View Article
- PubMed/NCBI
- Google Scholar
20. Team RDC. R: A Language and Environment for Statistical Computing. Austria: R Foundation for Statistical Computing, 2008.
21. Zeger SL, Liang KY. Longitudinal data analysis for discrete and continuous outcomes. Biometrics 1986; 42:121–130. pmid:3719049
- View Article
- PubMed/NCBI
- Google Scholar
22. McCullagh Nelder. Generalized Linear Models: Chapman and Hall Lt. Boca Raton, FL, 1989.
23. Oberg AL, Bot BM, Grill DE, Poland GA, Therneau TM. Technical and biological variance structure in mRNA-Seq data: life in the real world. BMC Genomics 2012; 13:304. pmid:22769017
- View Article
- PubMed/NCBI
- Google Scholar
24. Storey JD, Tibshirani R. Statistical significance for genomewide studies. Proc Natl Acad Sci USA 2003; 100:9440–9445. pmid:12883005
- View Article
- PubMed/NCBI
- Google Scholar
25. Fridley BL, Jenkins GD, Grill DE, Kennedy RB, Poland GA, Oberg AL. Soft truncation thresholding for gene set analysis of RNA-seq data: application to a vaccine study. Sci Rep 2013; 3:2898. pmid:24104466
- View Article
- PubMed/NCBI
- Google Scholar
26. Prasad TS, Kandasamy K, Pandey A. Human Protein Reference Database and Human Proteinpedia as discovery tools for systems biology. Methods Mol Biol 2009; 577:67–79. pmid:19718509
- View Article
- PubMed/NCBI
- Google Scholar
27. Rolland T, Tasan M, Charloteaux B, Pevzner SJ, Zhong Q, Sahni N, et al. A proteome-scale map of the human interactome network. Cell 2014; 159:1212–26. pmid:25416956
- View Article
- PubMed/NCBI
- Google Scholar
28. Szklarczyk D, Franceschini A, Kuhn M, Simonovic M, Roth A, Minguez P, et al. The STRING database in 2011: functional interaction networks of proteins, globally integrated and scored. Nucl Acids Res 2011; 39:D561–8. pmid:21045058
- View Article
- PubMed/NCBI
- Google Scholar
29. Lee I, Blom UM, Wang PI, Shim JE, Marcotte EM. Prioritizing candidate disease genes by network-based boosting of genome-wide association data. Genome Res 2011; 21:1109–21. pmid:21536720
- View Article
- PubMed/NCBI
- Google Scholar
30. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res 2003; 13:2498–504. pmid:14597658
- View Article
- PubMed/NCBI
- Google Scholar
31. AllegroViva. AllegroLayout. Available at: http://allegroviva.com/allegrolayout2.
32. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucl Acids Res 2000; 28:27–30. pmid:10592173
- View Article
- PubMed/NCBI
- Google Scholar
33. Kanehisa M, Goto S, Sato Y, Kawashima M, Furumichi M, Tanabe M. Data, information, knowledge and principle: back to metabolism in KEGG. Nucl Acids Res 2014; 42:D199–205. pmid:24214961
- View Article
- PubMed/NCBI
- Google Scholar
34. Rhee SY, Wood V, Dolinski K, Draghici S. Use and misuse of the gene ontology annotations. Nat Rev Genet 2008; 9:509–15. pmid:18475267
- View Article
- PubMed/NCBI
- Google Scholar
35. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA 2005; 102:15545–15550. pmid:16199517
- View Article
- PubMed/NCBI
- Google Scholar
36. Nakaya HI, Wrammert J, Lee EK, Racioppi L, Marie-Kunze S, Haining WN, et al. Systems biology of seasonal influenza vaccination in humans. Nat Immunol 2011; 12:786–795. pmid:21743478
- View Article
- PubMed/NCBI
- Google Scholar
37. Querec TD, Akondy RS, Lee EK, Cao W, Nakaya HI, Teuwen D, et al. Systems biology approach predicts immunogenicity of the yellow fever vaccine in humans. Nat Immunol 2009; 10:116–125. pmid:19029902
- View Article
- PubMed/NCBI
- Google Scholar
38. Obermoser G, Presnell S, Domico K, Xu H, Wang Y, Anguiano E, et al. Systems scale interactive exploration reveals quantitative and qualitative differences in response to influenza and pneumococcal vaccines. Immunity 2013; 38:831–44. pmid:23601689
- View Article
- PubMed/NCBI
- Google Scholar
39. Chevrier S, Genton C, Kallies A, Karnowski A, Otten LA, Malissen B, et al. CD93 is required for maintenance of antibody secretion and persistence of plasma cells in the bone marrow niche. Proc Natl Acad Sci USA 2009; 106:3895–900. pmid:19228948
- View Article
- PubMed/NCBI
- Google Scholar
40. Amanna IJ, Carlson NE, Slifka MK. Duration of humoral immunity to common viral and vaccine antigens. N Engl J Med 2007; 357:1903–1915. pmid:17989383
- View Article
- PubMed/NCBI
- Google Scholar
41. Getahun A, Smith MJ, Kogut I, van Dyk LF, Cambier JC (2012) Retention of anergy and inhibition of antibody responses during acute gamma herpesvirus 68 infection. Journal of Immunology 189: 2965–2974.
- View Article
- Google Scholar
42. Zhang M, Bohlson SS, Dy M, Tenner AJ. Modulated interaction of the ERM protein, moesin, with CD93. Immunol 2005; 115:63–73.
- View Article
- Google Scholar
43. Maarof G, Bouchet-Delbos L, Gary-Gouy H, Durand-Gasselin I, Krzysiek R, Dalloul A. Interleukin-24 inhibits the plasma cell differentiation program in human germinal center B cells. Blood 2010; 115:1718–26. pmid:19965640
- View Article
- PubMed/NCBI
- Google Scholar
44. Radbruch A, Muehlinghaus G, Luger EO, Inamine A, Smith KG, Dörner T, et al. Competence and competition: the challenge of becoming a long-lived plasma cell. Nature Rev Immunol 2006; 6:741–50.
- View Article
- Google Scholar
45. Painter SD, Ovsyannikova IG, Poland GA. The weight of obesity on the human immune response to vaccination. Vaccine 2015; 33:4422–9. pmid:26163925
- View Article
- PubMed/NCBI
- Google Scholar
46. Zilliox MJ, Parmigiani G, Griffin DE. Gene expression patterns in dendritic cells infected with measles virus compared with other pathogens. Proc Natl Acad Sci USA 2006; 103:3363–3368. pmid:16492729
- View Article
- PubMed/NCBI
- Google Scholar
47. Zilliox MJ, Moss WJ, Griffin DE. Gene expression changes in peripheral blood mononuclear cells during measles virus infection. Clin Vaccine Immunol 2007; 14:918–923. pmid:17538120
- View Article
- PubMed/NCBI
- Google Scholar
48. Dhiman N, Ovsyannikova IG, Oberg AL, Grill DE, Jacobson RM, Poland GA. Immune activation at effector and gene expression levels after measles vaccination in healthy individuals: a pilot study. Human Immunol 2005; 66:1125–1136.
- View Article
- Google Scholar
49. Mocellin S, Panelli MC, Wang E, Nagorsen D, Marincola FM. The dual role of IL-10. Trends Immunol 2003; 24:36–43. pmid:12495723
- View Article
- PubMed/NCBI
- Google Scholar
50. Karsten CM, Kohl J. The complement receptor CD46 tips the scales in T(H)1 self-control. Nature Immunol 2010; 11:775–7.
- View Article
- Google Scholar

[ref1] 1. Whitaker JA, Poland GA. Measles and mumps outbreaks in the United States: Think globally, vaccinate locally. Vaccine 2014; 32:4703–4. pmid:24992719
View Article
PubMed/NCBI
Google Scholar

[2] View Article

[3] PubMed/NCBI

[4] Google Scholar

[ref2] 2. Centers for Disease Control and Prevention. Measles Cases and Outbreaks. Available at: http://www.cdc.gov/measles/cases-outbreaks.html. Date accessed: February 12, 2015.

[ref3] 3. Parker Fiebelkorn A, Redd SB, Gallagher K, et al. Measles in the United States during the postelimination era. J Infect Dis 2010; 202:1520–1528. pmid:20929352
View Article
PubMed/NCBI
Google Scholar

[7] View Article

[8] PubMed/NCBI

[9] Google Scholar

[ref4] 4. Increased transmission and outbreaks of measles—European Region, 2011. MMWR 2011; 60:1605–10. pmid:22129994
View Article
PubMed/NCBI
Google Scholar

[11] View Article

[12] PubMed/NCBI

[13] Google Scholar

[ref5] 5. Poland GA, Ovsyannikova IG, Kennedy RB, Haralambieva IH, Jacobson RM. Vaccinomics and a new paradigm for the development of preventive vaccines against viral infections. Omics 2011; 15:625–636. pmid:21732819
View Article
PubMed/NCBI
Google Scholar

[15] View Article

[16] PubMed/NCBI

[17] Google Scholar

[ref6] 6. Haralambieva IH, Poland GA. Vaccinomics, predictive vaccinology and the future of vaccine development. Future Microbiol 2010; 5:1757–1760. pmid:21155658
View Article
PubMed/NCBI
Google Scholar

[19] View Article

[20] PubMed/NCBI

[21] Google Scholar

[ref7] 7. Lambert ND, Ovsyannikova IG, Pankratz VS, Jacobson RM, Poland GA. Understanding the immune response to seasonal influenza vaccination in older adults: a systems biology approach. Exp Rev Vaccines 2012; 11:985–994.
View Article
Google Scholar

[23] View Article

[24] Google Scholar

[ref8] 8. Kennedy RB, Ovsyannikova IG, Haralambieva IH, O'Byrne MM, Jacobson RM, Pankratz VS, et al. Multigenic control of measles vaccine immunity mediated by polymorphisms in measles receptor, innate pathway, and cytokine genes. Vaccine 2012; 30:2159–2167. pmid:22265947
View Article
PubMed/NCBI
Google Scholar

[26] View Article

[27] PubMed/NCBI

[28] Google Scholar

[ref9] 9. Ovsyannikova IG, Haralambieva IH, Vierkant RA, O'Byrne MM, Jacobson RM, Poland GA. The association of CD46, SLAM, and CD209 cellular receptor gene SNPs with variations in measles vaccine-induced immune responses—a replication study and examination of novel polymorphisms. Human Hered 2011; 72:206–223.
View Article
Google Scholar

[30] View Article

[31] Google Scholar

[ref10] 10. Haralambieva IH, Ovsyannikova IG, Umlauf BJ, Vierkant RA, Shane Pankratz V, Jacobson RM, et al. Genetic polymorphisms in host antiviral genes: associations with humoral and cellular immunity to measles vaccine. Vaccine 2011; 29:8988–8997. pmid:21939710
View Article
PubMed/NCBI
Google Scholar

[33] View Article

[34] PubMed/NCBI

[35] Google Scholar

[ref11] 11. Haralambieva IH, Ovsyannikova IG, Kennedy RB, Vierkant RA, Pankratz VS, Jacobson RM, et al. Associations between single nucleotide polymorphisms and haplotypes in cytokine and cytokine receptor genes and immunity to measles vaccination. Vaccine 2011; 29:7883–7895. pmid:21875636
View Article
PubMed/NCBI
Google Scholar

[37] View Article

[38] PubMed/NCBI

[39] Google Scholar

[ref12] 12. Haralambieva IH, Ovsyannikova IG, O'Byrne M, Pankratz VS, Jacobson RM, Poland GA. A large observational study to concurrently assess persistence of measles specific B-cell and T-cell immunity in individuals following two doses of MMR vaccine. Vaccine 2011; 29:4485–4491. pmid:21539880
View Article
PubMed/NCBI
Google Scholar

[41] View Article

[42] PubMed/NCBI

[43] Google Scholar

[ref13] 13. Haralambieva IH, Oberg AL, Ovsyannikova IG, Kennedy RB, Grill DE, Middha S, et al. Genome-wide characterization of transcriptional patterns in high and low antibody responders to rubella vaccination. PLos ONE 2013; 8:e62149. pmid:23658707
View Article
PubMed/NCBI
Google Scholar

[45] View Article

[46] PubMed/NCBI

[47] Google Scholar

[ref14] 14. Kennedy RB, Oberg AL, Ovsyannikova IG, Haralambieva IH, Grill DE, Poland GA. Transcriptomic profiles of high and low antibody responders to smallpox vaccine. Genes Immunity 2013; 14:277–285. pmid:23594957
View Article
PubMed/NCBI
Google Scholar

[49] View Article

[50] PubMed/NCBI

[51] Google Scholar

[ref15] 15. Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics 2010; 26:841–2. pmid:20110278
View Article
PubMed/NCBI
Google Scholar

[53] View Article

[54] PubMed/NCBI

[55] Google Scholar

[ref16] 16. Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics 2009; 25:1105–11. pmid:19289445
View Article
PubMed/NCBI
Google Scholar

[57] View Article

[58] PubMed/NCBI

[59] Google Scholar

[ref17] 17. Anders S, Pyl PT, Huber W. HTSeq—a Python framework to work with high-throughput sequencing data. Bioinformatics 2015; 31:166–9. pmid:25260700
View Article
PubMed/NCBI
Google Scholar

[61] View Article

[62] PubMed/NCBI

[63] Google Scholar

[ref18] 18. Oberg AL, Kennedy RB, Li P, Ovsyannikova IG, Poland GA. Systems biology approaches to new vaccine development. Curr Opin Immunol 2011; 23:436–443. pmid:21570272
View Article
PubMed/NCBI
Google Scholar

[65] View Article

[66] PubMed/NCBI

[67] Google Scholar

[ref19] 19. Hansen KD, Irizarry RA, Wu Z. Removing technical variability in RNA-seq data using conditional quantile normalization. Biostatistics 2012; 13:204–216. pmid:22285995
View Article
PubMed/NCBI
Google Scholar

[69] View Article

[70] PubMed/NCBI

[71] Google Scholar

[ref20] 20. Team RDC. R: A Language and Environment for Statistical Computing. Austria: R Foundation for Statistical Computing, 2008.

[ref21] 21. Zeger SL, Liang KY. Longitudinal data analysis for discrete and continuous outcomes. Biometrics 1986; 42:121–130. pmid:3719049
View Article
PubMed/NCBI
Google Scholar

[74] View Article

[75] PubMed/NCBI

[76] Google Scholar

[ref22] 22. McCullagh Nelder. Generalized Linear Models: Chapman and Hall Lt. Boca Raton, FL, 1989.

[ref23] 23. Oberg AL, Bot BM, Grill DE, Poland GA, Therneau TM. Technical and biological variance structure in mRNA-Seq data: life in the real world. BMC Genomics 2012; 13:304. pmid:22769017
View Article
PubMed/NCBI
Google Scholar

[79] View Article

[80] PubMed/NCBI

[81] Google Scholar

[ref24] 24. Storey JD, Tibshirani R. Statistical significance for genomewide studies. Proc Natl Acad Sci USA 2003; 100:9440–9445. pmid:12883005
View Article
PubMed/NCBI
Google Scholar

[83] View Article

[84] PubMed/NCBI

[85] Google Scholar

[ref25] 25. Fridley BL, Jenkins GD, Grill DE, Kennedy RB, Poland GA, Oberg AL. Soft truncation thresholding for gene set analysis of RNA-seq data: application to a vaccine study. Sci Rep 2013; 3:2898. pmid:24104466
View Article
PubMed/NCBI
Google Scholar

[87] View Article

[88] PubMed/NCBI

[89] Google Scholar

[ref26] 26. Prasad TS, Kandasamy K, Pandey A. Human Protein Reference Database and Human Proteinpedia as discovery tools for systems biology. Methods Mol Biol 2009; 577:67–79. pmid:19718509
View Article
PubMed/NCBI
Google Scholar

[91] View Article

[92] PubMed/NCBI

[93] Google Scholar

[ref27] 27. Rolland T, Tasan M, Charloteaux B, Pevzner SJ, Zhong Q, Sahni N, et al. A proteome-scale map of the human interactome network. Cell 2014; 159:1212–26. pmid:25416956
View Article
PubMed/NCBI
Google Scholar

[95] View Article

[96] PubMed/NCBI

[97] Google Scholar

[ref28] 28. Szklarczyk D, Franceschini A, Kuhn M, Simonovic M, Roth A, Minguez P, et al. The STRING database in 2011: functional interaction networks of proteins, globally integrated and scored. Nucl Acids Res 2011; 39:D561–8. pmid:21045058
View Article
PubMed/NCBI
Google Scholar

[99] View Article

[100] PubMed/NCBI

[101] Google Scholar

[ref29] 29. Lee I, Blom UM, Wang PI, Shim JE, Marcotte EM. Prioritizing candidate disease genes by network-based boosting of genome-wide association data. Genome Res 2011; 21:1109–21. pmid:21536720
View Article
PubMed/NCBI
Google Scholar

[103] View Article

[104] PubMed/NCBI

[105] Google Scholar

[ref30] 30. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res 2003; 13:2498–504. pmid:14597658
View Article
PubMed/NCBI
Google Scholar

[107] View Article

[108] PubMed/NCBI

[109] Google Scholar

[ref31] 31. AllegroViva. AllegroLayout. Available at: http://allegroviva.com/allegrolayout2.

[ref32] 32. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucl Acids Res 2000; 28:27–30. pmid:10592173
View Article
PubMed/NCBI
Google Scholar

[112] View Article

[113] PubMed/NCBI

[114] Google Scholar

[ref33] 33. Kanehisa M, Goto S, Sato Y, Kawashima M, Furumichi M, Tanabe M. Data, information, knowledge and principle: back to metabolism in KEGG. Nucl Acids Res 2014; 42:D199–205. pmid:24214961
View Article
PubMed/NCBI
Google Scholar

[116] View Article

[117] PubMed/NCBI

[118] Google Scholar

[ref34] 34. Rhee SY, Wood V, Dolinski K, Draghici S. Use and misuse of the gene ontology annotations. Nat Rev Genet 2008; 9:509–15. pmid:18475267
View Article
PubMed/NCBI
Google Scholar

[120] View Article

[121] PubMed/NCBI

[122] Google Scholar

[ref35] 35. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA 2005; 102:15545–15550. pmid:16199517
View Article
PubMed/NCBI
Google Scholar

[124] View Article

[125] PubMed/NCBI

[126] Google Scholar

[ref36] 36. Nakaya HI, Wrammert J, Lee EK, Racioppi L, Marie-Kunze S, Haining WN, et al. Systems biology of seasonal influenza vaccination in humans. Nat Immunol 2011; 12:786–795. pmid:21743478
View Article
PubMed/NCBI
Google Scholar

[128] View Article

[129] PubMed/NCBI

[130] Google Scholar

[ref37] 37. Querec TD, Akondy RS, Lee EK, Cao W, Nakaya HI, Teuwen D, et al. Systems biology approach predicts immunogenicity of the yellow fever vaccine in humans. Nat Immunol 2009; 10:116–125. pmid:19029902
View Article
PubMed/NCBI
Google Scholar

[132] View Article

[133] PubMed/NCBI

[134] Google Scholar

[ref38] 38. Obermoser G, Presnell S, Domico K, Xu H, Wang Y, Anguiano E, et al. Systems scale interactive exploration reveals quantitative and qualitative differences in response to influenza and pneumococcal vaccines. Immunity 2013; 38:831–44. pmid:23601689
View Article
PubMed/NCBI
Google Scholar

[136] View Article

[137] PubMed/NCBI

[138] Google Scholar

[ref39] 39. Chevrier S, Genton C, Kallies A, Karnowski A, Otten LA, Malissen B, et al. CD93 is required for maintenance of antibody secretion and persistence of plasma cells in the bone marrow niche. Proc Natl Acad Sci USA 2009; 106:3895–900. pmid:19228948
View Article
PubMed/NCBI
Google Scholar

[140] View Article

[141] PubMed/NCBI

[142] Google Scholar

[ref40] 40. Amanna IJ, Carlson NE, Slifka MK. Duration of humoral immunity to common viral and vaccine antigens. N Engl J Med 2007; 357:1903–1915. pmid:17989383
View Article
PubMed/NCBI
Google Scholar

[144] View Article

[145] PubMed/NCBI

[146] Google Scholar

[ref41] 41. Getahun A, Smith MJ, Kogut I, van Dyk LF, Cambier JC (2012) Retention of anergy and inhibition of antibody responses during acute gamma herpesvirus 68 infection. Journal of Immunology 189: 2965–2974.
View Article
Google Scholar

[148] View Article

[149] Google Scholar

[ref42] 42. Zhang M, Bohlson SS, Dy M, Tenner AJ. Modulated interaction of the ERM protein, moesin, with CD93. Immunol 2005; 115:63–73.
View Article
Google Scholar

[151] View Article

[152] Google Scholar

[ref43] 43. Maarof G, Bouchet-Delbos L, Gary-Gouy H, Durand-Gasselin I, Krzysiek R, Dalloul A. Interleukin-24 inhibits the plasma cell differentiation program in human germinal center B cells. Blood 2010; 115:1718–26. pmid:19965640
View Article
PubMed/NCBI
Google Scholar

[154] View Article

[155] PubMed/NCBI

[156] Google Scholar

[ref44] 44. Radbruch A, Muehlinghaus G, Luger EO, Inamine A, Smith KG, Dörner T, et al. Competence and competition: the challenge of becoming a long-lived plasma cell. Nature Rev Immunol 2006; 6:741–50.
View Article
Google Scholar

[158] View Article

[159] Google Scholar

[ref45] 45. Painter SD, Ovsyannikova IG, Poland GA. The weight of obesity on the human immune response to vaccination. Vaccine 2015; 33:4422–9. pmid:26163925
View Article
PubMed/NCBI
Google Scholar

[161] View Article

[162] PubMed/NCBI

[163] Google Scholar

[ref46] 46. Zilliox MJ, Parmigiani G, Griffin DE. Gene expression patterns in dendritic cells infected with measles virus compared with other pathogens. Proc Natl Acad Sci USA 2006; 103:3363–3368. pmid:16492729
View Article
PubMed/NCBI
Google Scholar

[165] View Article

[166] PubMed/NCBI

[167] Google Scholar

[ref47] 47. Zilliox MJ, Moss WJ, Griffin DE. Gene expression changes in peripheral blood mononuclear cells during measles virus infection. Clin Vaccine Immunol 2007; 14:918–923. pmid:17538120
View Article
PubMed/NCBI
Google Scholar

[169] View Article

[170] PubMed/NCBI

[171] Google Scholar

[ref48] 48. Dhiman N, Ovsyannikova IG, Oberg AL, Grill DE, Jacobson RM, Poland GA. Immune activation at effector and gene expression levels after measles vaccination in healthy individuals: a pilot study. Human Immunol 2005; 66:1125–1136.
View Article
Google Scholar

[173] View Article

[174] Google Scholar

[ref49] 49. Mocellin S, Panelli MC, Wang E, Nagorsen D, Marincola FM. The dual role of IL-10. Trends Immunol 2003; 24:36–43. pmid:12495723
View Article
PubMed/NCBI
Google Scholar

[176] View Article

[177] PubMed/NCBI

[178] Google Scholar

[ref50] 50. Karsten CM, Kohl J. The complement receptor CD46 tips the scales in T(H)1 self-control. Nature Immunol 2010; 11:775–7.
View Article
Google Scholar

[180] View Article

[181] Google Scholar

Figures

Abstract

Background

Methods

Results

Conclusion

Introduction

Methods

2.1. Study subjects

2.2. Plaque reduction microneutralization assay (PRMN)

2.3. mRNA-Seq transcriptome profiling

2.4. Statistical analysis

2.5. Network biology

Results

3.1. Characteristics of the study subjects

3.2. Gene expression after in vitro stimulation with measles virus in high and low antibody responders (interaction analysis)

3.3. Pathway and geneset analysis

3.4. Network biology

Discussion

Supporting Information

S1 Table. Immune response variables of the study subjects.

S2 Table. Significant genes based on overall response to measles virus stimulation in all vaccine recipients (1,761 significant genes with FC>2 and FC <0.5).

S3 Table. Significant genesets, in high vs low antibody vaccine responders (response to viral stimulation)

Acknowledgments

Author Contributions

References