Integration of immune cell populations, mRNA-Seq, and CpG methylation to better predict humoral immunity to influenza vaccination: Dependence of mRNA-Seq/CpG methylation on immune cell populations

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

Research output: Contribution to journalArticle

5 Citations (Scopus)

Abstract

The development of a humoral immune response to influenza vaccines occurs on a multisystems level. Due to the orchestration required for robust immune responses when multiple genes and their regulatory components across multiple cell types are involved, we examined an influenza vaccination cohort using multiple high-throughput technologies. In this study, we sought a more thorough understanding of how immune cell composition and gene expression relate to each other and contribute to interindividual variation in response to influenza vaccination. We first hypothesized that many of the differentially expressed (DE) genes observed after influenza vaccination result from changes in the composition of participants' peripheral blood mononuclear cells (PBMCs), which were assessed using flow cytometry. We demonstrated that DE genes in our study are correlated with changes in PBMC composition. We gathered DE genes from 128 other publically available PBMC-based vaccine studies and identified that an average of 57% correlated with specific cell subset levels in our study (permutation used to control false discovery), suggesting that the associations we have identified are likely general features of PBMC-based transcriptomics. Second, we hypothesized that more robust models of vaccine response could be generated by accounting for the interplay between PBMC composition, gene expression, and gene regulation. We employed machine learning to generate predictive models of B-cell ELISPOT response outcomes and hemagglutination inhibition (HAI) antibody titers. The top HAI and B-cell ELISPOT model achieved an area under the receiver operating curve (AUC) of 0.64 and 0.79, respectively, with linear model coefficients of determination of 0.08 and 0.28. For the B-cell ELISPOT outcomes, CpG methylation had the greatest predictive ability, highlighting potentially novel regulatory features important for immune response. B-cell ELISOT models using only PBMC composition had lower performance (AUC = 0.67), but highlighted well-known mechanisms. Our analysis demonstrated that each of the three data sets (cell composition, mRNA-Seq,and DNA methylation) may provide distinct information for the prediction of humoral immune response outcomes. We believe that these findings are important for the interpretation of current omics-based studies and set the stage for a more thorough understanding of interindividual immune responses to influenza vaccination.

Original languageEnglish (US)
Article number445
JournalFrontiers in Immunology
Volume8
Issue numberAPR
DOIs
StatePublished - Apr 21 2017

Fingerprint

Humoral Immunity
Methylation
Human Influenza
Blood Cells
Vaccination
Enzyme-Linked Immunospot Assay
Messenger RNA
B-Lymphocytes
Population
Hemagglutination
Genes
Area Under Curve
Vaccines
Influenza Vaccines
Gene Expression Regulation
DNA Methylation
Regulator Genes
Linear Models
Flow Cytometry
Technology

Keywords

  • Cell sorting
  • Data mining
  • Differential expression
  • Immunology
  • Influenza vaccine
  • Machine learning
  • Methylation

ASJC Scopus subject areas

  • Immunology and Allergy
  • Immunology

Cite this

Integration of immune cell populations, mRNA-Seq, and CpG methylation to better predict humoral immunity to influenza vaccination : Dependence of mRNA-Seq/CpG methylation on immune cell populations. / Zimmermann, Michael T.; Kennedy, Richard B; Grill, Diane E.; Oberg, Ann L; Goergen, Krista M.; Ovsyannikova, Inna G.; Haralambieva, Iana H.; Poland, Gregory A.

In: Frontiers in Immunology, Vol. 8, No. APR, 445, 21.04.2017.

Research output: Contribution to journalArticle

@article{0452fa03d19c471a8251c8d89190e5cb,
title = "Integration of immune cell populations, mRNA-Seq, and CpG methylation to better predict humoral immunity to influenza vaccination: Dependence of mRNA-Seq/CpG methylation on immune cell populations",
abstract = "The development of a humoral immune response to influenza vaccines occurs on a multisystems level. Due to the orchestration required for robust immune responses when multiple genes and their regulatory components across multiple cell types are involved, we examined an influenza vaccination cohort using multiple high-throughput technologies. In this study, we sought a more thorough understanding of how immune cell composition and gene expression relate to each other and contribute to interindividual variation in response to influenza vaccination. We first hypothesized that many of the differentially expressed (DE) genes observed after influenza vaccination result from changes in the composition of participants' peripheral blood mononuclear cells (PBMCs), which were assessed using flow cytometry. We demonstrated that DE genes in our study are correlated with changes in PBMC composition. We gathered DE genes from 128 other publically available PBMC-based vaccine studies and identified that an average of 57{\%} correlated with specific cell subset levels in our study (permutation used to control false discovery), suggesting that the associations we have identified are likely general features of PBMC-based transcriptomics. Second, we hypothesized that more robust models of vaccine response could be generated by accounting for the interplay between PBMC composition, gene expression, and gene regulation. We employed machine learning to generate predictive models of B-cell ELISPOT response outcomes and hemagglutination inhibition (HAI) antibody titers. The top HAI and B-cell ELISPOT model achieved an area under the receiver operating curve (AUC) of 0.64 and 0.79, respectively, with linear model coefficients of determination of 0.08 and 0.28. For the B-cell ELISPOT outcomes, CpG methylation had the greatest predictive ability, highlighting potentially novel regulatory features important for immune response. B-cell ELISOT models using only PBMC composition had lower performance (AUC = 0.67), but highlighted well-known mechanisms. Our analysis demonstrated that each of the three data sets (cell composition, mRNA-Seq,and DNA methylation) may provide distinct information for the prediction of humoral immune response outcomes. We believe that these findings are important for the interpretation of current omics-based studies and set the stage for a more thorough understanding of interindividual immune responses to influenza vaccination.",
keywords = "Cell sorting, Data mining, Differential expression, Immunology, Influenza vaccine, Machine learning, Methylation",
author = "Zimmermann, {Michael T.} and Kennedy, {Richard B} and Grill, {Diane E.} and Oberg, {Ann L} and Goergen, {Krista M.} and Ovsyannikova, {Inna G.} and Haralambieva, {Iana H.} and Poland, {Gregory A.}",
year = "2017",
month = "4",
day = "21",
doi = "10.3389/fimmu.2017.00445",
language = "English (US)",
volume = "8",
journal = "Frontiers in Immunology",
issn = "1664-3224",
publisher = "Frontiers Media S. A.",
number = "APR",

}

TY - JOUR

T1 - Integration of immune cell populations, mRNA-Seq, and CpG methylation to better predict humoral immunity to influenza vaccination

T2 - Dependence of mRNA-Seq/CpG methylation on immune cell populations

AU - Zimmermann, Michael T.

AU - Kennedy, Richard B

AU - Grill, Diane E.

AU - Oberg, Ann L

AU - Goergen, Krista M.

AU - Ovsyannikova, Inna G.

AU - Haralambieva, Iana H.

AU - Poland, Gregory A.

PY - 2017/4/21

Y1 - 2017/4/21

N2 - The development of a humoral immune response to influenza vaccines occurs on a multisystems level. Due to the orchestration required for robust immune responses when multiple genes and their regulatory components across multiple cell types are involved, we examined an influenza vaccination cohort using multiple high-throughput technologies. In this study, we sought a more thorough understanding of how immune cell composition and gene expression relate to each other and contribute to interindividual variation in response to influenza vaccination. We first hypothesized that many of the differentially expressed (DE) genes observed after influenza vaccination result from changes in the composition of participants' peripheral blood mononuclear cells (PBMCs), which were assessed using flow cytometry. We demonstrated that DE genes in our study are correlated with changes in PBMC composition. We gathered DE genes from 128 other publically available PBMC-based vaccine studies and identified that an average of 57% correlated with specific cell subset levels in our study (permutation used to control false discovery), suggesting that the associations we have identified are likely general features of PBMC-based transcriptomics. Second, we hypothesized that more robust models of vaccine response could be generated by accounting for the interplay between PBMC composition, gene expression, and gene regulation. We employed machine learning to generate predictive models of B-cell ELISPOT response outcomes and hemagglutination inhibition (HAI) antibody titers. The top HAI and B-cell ELISPOT model achieved an area under the receiver operating curve (AUC) of 0.64 and 0.79, respectively, with linear model coefficients of determination of 0.08 and 0.28. For the B-cell ELISPOT outcomes, CpG methylation had the greatest predictive ability, highlighting potentially novel regulatory features important for immune response. B-cell ELISOT models using only PBMC composition had lower performance (AUC = 0.67), but highlighted well-known mechanisms. Our analysis demonstrated that each of the three data sets (cell composition, mRNA-Seq,and DNA methylation) may provide distinct information for the prediction of humoral immune response outcomes. We believe that these findings are important for the interpretation of current omics-based studies and set the stage for a more thorough understanding of interindividual immune responses to influenza vaccination.

AB - The development of a humoral immune response to influenza vaccines occurs on a multisystems level. Due to the orchestration required for robust immune responses when multiple genes and their regulatory components across multiple cell types are involved, we examined an influenza vaccination cohort using multiple high-throughput technologies. In this study, we sought a more thorough understanding of how immune cell composition and gene expression relate to each other and contribute to interindividual variation in response to influenza vaccination. We first hypothesized that many of the differentially expressed (DE) genes observed after influenza vaccination result from changes in the composition of participants' peripheral blood mononuclear cells (PBMCs), which were assessed using flow cytometry. We demonstrated that DE genes in our study are correlated with changes in PBMC composition. We gathered DE genes from 128 other publically available PBMC-based vaccine studies and identified that an average of 57% correlated with specific cell subset levels in our study (permutation used to control false discovery), suggesting that the associations we have identified are likely general features of PBMC-based transcriptomics. Second, we hypothesized that more robust models of vaccine response could be generated by accounting for the interplay between PBMC composition, gene expression, and gene regulation. We employed machine learning to generate predictive models of B-cell ELISPOT response outcomes and hemagglutination inhibition (HAI) antibody titers. The top HAI and B-cell ELISPOT model achieved an area under the receiver operating curve (AUC) of 0.64 and 0.79, respectively, with linear model coefficients of determination of 0.08 and 0.28. For the B-cell ELISPOT outcomes, CpG methylation had the greatest predictive ability, highlighting potentially novel regulatory features important for immune response. B-cell ELISOT models using only PBMC composition had lower performance (AUC = 0.67), but highlighted well-known mechanisms. Our analysis demonstrated that each of the three data sets (cell composition, mRNA-Seq,and DNA methylation) may provide distinct information for the prediction of humoral immune response outcomes. We believe that these findings are important for the interpretation of current omics-based studies and set the stage for a more thorough understanding of interindividual immune responses to influenza vaccination.

KW - Cell sorting

KW - Data mining

KW - Differential expression

KW - Immunology

KW - Influenza vaccine

KW - Machine learning

KW - Methylation

UR - http://www.scopus.com/inward/record.url?scp=85018454015&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85018454015&partnerID=8YFLogxK

U2 - 10.3389/fimmu.2017.00445

DO - 10.3389/fimmu.2017.00445

M3 - Article

AN - SCOPUS:85018454015

VL - 8

JO - Frontiers in Immunology

JF - Frontiers in Immunology

SN - 1664-3224

IS - APR

M1 - 445

ER -