- Short Report
- Open Access
Using gene expression profiles from peripheral blood to identify asymptomatic responses to acute respiratory viral infections
BMC Research Notes volume 3, Article number: 264 (2010)
A recent study reported that gene expression profiles from peripheral blood samples of healthy subjects prior to viral inoculation were indistinguishable from profiles of subjects who received viral challenge but remained asymptomatic and uninfected. If true, this implies that the host immune response does not have a molecular signature. Given the high sensitivity of microarray technology, we were intrigued by this result and hypothesize that it was an artifact of data analysis.
Using acute respiratory viral challenge microarray data, we developed a molecular signature that for the first time allowed for an accurate differentiation between uninfected subjects prior to viral inoculation and subjects who remained asymptomatic after the viral challenge.
Our findings suggest that molecular signatures can be used to characterize immune responses to viruses and may improve our understanding of susceptibility to viral infection with possible implications for vaccine development.
Acute respiratory viral infections cause significant morbidity and mortality in the United States and worldwide. Currently clinicians do not have practical means to make a timely and accurate diagnosis of acute viral respiratory infections and often resort to unnecessary antibiotic treatment, which increases healthcare costs and facilitates development of antibiotic resistance. Recently a novel approach was proposed for the diagnosis of acute respiratory infections based on microarray gene expression profiles from peripheral blood samples from human subjects . Using human viral challenge studies with live human rhinovirus (HRV), respiratory syncytial virus (RSV), and influenza A, Zaas et al. developed an "acute respiratory viral response" 30-gene panviral signature that accurately identified symptomatic patients with viral infection. The same study also reported that gene expression profiles of subjects prior to viral inoculation (at baseline) were indistinguishable from profiles of subjects who received viral challenge but remained asymptomatic and uninfected. Given the high sensitivity of microarray technology, we were intrigued by the latter result and hypothesize that it was an artifact of data analysis. Since the gene expression dataset of Zaas et al. was deposited in the Gene Expression Omnibus, we were able to verify our hypothesis and discover significant differences between these two groups of samples. We demonstrate this by developing a molecular signature that discriminates with high accuracy between uninfected subjects at baseline (prior to viral inoculation) and asymptomatic subjects at the time correlating to peak symptoms in the symptomatic group. This finding has important implications for a better understanding of the complex human immune response to viral antigens. Genes that are differentially expressed in the two groups may provide important clues about decreased susceptibility to viral challenge, which could result in more effective vaccine development or novel therapeutic strategies.
To arrive at our findings, we reanalyzed the gene expression data of Zaas et al. that contains measurements of 22,277 oligonucleotide probes for 56 blood samples from uninfected subjects (measured at baseline, prior to inoculation) and 30 blood samples from subjects who received viral challenge but remained asymptomatic (measured at corresponding peak symptoms time specific for each virus). We used an improved data-analytic protocol that avoids selection of redundant and biologically irrelevant genes and at the same time maximizes predictive accuracy of the signature . This protocol has recently allowed us to develop a highly accurate and compact molecular signature for differentiation between uninfected subjects and those with acute respiratory viral infections. This signature had high reproducibility as evidenced by its nearly perfect accuracy in the independent data of Ramilo et al., and was comprised of genes that are involved in the host immune response .
Our data-analytic protocol first involved selection of genes by GLL-PC, a supervised multivariate biomarker discovery method that provably discovers genes in the local pathway around the response variable of interest [4, 5]; additional details about GLL-PC are provided in [Additional file 1]. Next, Support Vector Machine (SVM) classifiers with linear kernel and penalty hyper-parameter C = 100 were fitted on the selected genes . In order to obtain an unbiased estimate of predictive accuracy that will hold in applications of this molecular signature on future subjects, gene selection and SVM classifier training were performed by stratified 10-fold cross-validation repeated 100 times for different splits of subject into 10 folds [2, 7]. Finally, we ensured signature reproducibility and assessed its statistical significance by using a permutation test with significance level α = 0.05 and 10,000 permutations .
Results and Discussion
The data-analytic protocol described above yields an unbiased estimate of predictive accuracy = 0.85 area under ROC curve (AUC); 95% confidence interval [0.76; 0.94] AUC. On average GLL-PC selected 7 genes depending on the training set of cross-validation. Genes selected by GLL-PC in more than 20% of the training sets are listed in Table 1. Genes that were selected in a smaller fraction of training sets may be artifacts of small sample size. Next, GLL-PC and SVM were applied on the entire set of samples, resulting in a 6-gene signature comprised of genes EIF2S1, ZNF91, RBM3, ATP5S, TPPP3, GPR97. Note that all these six genes were also among the top seven most frequently selected by GLL-PC during cross-validation (Table 1), which demonstrates the stability of this gene set despite the inherent heterogeneity of gene expression across subjects.
Of the eight genes shown in Table 1, four are known to be involved in the host immune response. EIF2S1 encodes the alpha subunit of the translation initiation factor eIF2 complex (eIF2α), which initiates protein synthesis. This protein appears to be crucial for the survival of virally infected cells . The increase of expression of EIF2S1 that is observed in the data could mean that host protein synthesis is turned on as response to viral exposure. On the other hand, many viruses, including influenza A viruses, have developed mechanisms favoring the translation of viral over cellular mRNAs leading to a switch from cellular to viral protein synthesis, while inhibiting the cell-encoded antiviral pathways . RBM3 encodes a protein that is a member of the heterogeneous nuclear ribonucleoproteins (hnRNPs) that appears to enhance global protein synthesis  but also has a role in viral transcription and replication . The ZNF91 gene family is a subset of the Krüppel-associated box (KRAB)-containing group of zinc finger genes. Zinc-finger proteins containing the KRAB domain are transcriptional regulators that have been associated with suppression of viral proliferation . The related transcriptional factor OTK18 has been shown to be induced by and to suppress HIV-1 infection in mononuclear cells . Finally, CD24 is involved in the CD24-Siglec G pathway that protects the host against a fatal reaction to pathological cell death and discriminates danger- versus pathogen-associated molecular patterns .
To the best of our knowledge, the remaining 4 genes from Table 1 (TPPP3, GPR97, ATP5S, TRAK2) have not yet been specifically associated with host immune responses, however this does not rule them out as novel pharmaceutical targets. Some of these genes are involved in normal cellular function and for some there is evidence of their association with host immune responses: TPPP3 encodes members of the tubulin family of proteins, and one study has shown an increase of tubulin in influenza infection . GPR97 encodes a member of G protein-coupled receptors that bind chemokines on the surface of immune cells. Chemokine signaling is critical for effective antiviral immune response activation. ATP5S encodes a subunit of mitochondrial ATP synthase, and it has been suggested that an increase in mitochondrial activity plays a relevant role in viral replication . TRAK2 encodes a trafficking factor that facilitates expression of potassium channels to the cell surface .
In further analysis, we identify a likely reason that precluded Zaas et al. from detecting a difference between uninfected subjects at baseline and asymptomatic subjects at corresponding peak time. The study of Zaas et al. used 30 genes from the "acute respiratory viral response" signature (that was designed for differentiation of symptomatic subjects from uninfected individuals) to perform a different classification task. It is therefore not surprising that even though these 30 genes were predictive for diagnosis of symptomatic subjects from uninfected individuals, they were not predictive for differentiating uninfected subjects at baseline from asymptomatic subjects at the time corresponding to peak symptoms in the symptomatic group. We have verified this assertion by running the original data analysis software of Zaas et al. that resulted in random predictive accuracy (0.50 AUC) when using 30 genes from the "acute respiratory viral response" signature and a non-trivial predictive accuracy (0.66 AUC) when using genes selected specifically for the classification task of interest. The remainder of the discrepancy between performance of the protocol of Zaas et al. and one used in the present study is due to the choice of methods for gene selection, classification, and accuracy estimation (cross-validation). If we substitute sparse probit regression with SVMs in the analysis protocol of Zaas et al., the predictive accuracy increases to 0.75 AUC. If we also substitute factor analysis-based gene selection with GLL-PC and use repeated 10-fold cross-validation, we obtain the result identical to our analysis, i.e. 0.85 AUC.
In summary, upon reanalysis, the data provided in the study by Zaas et al. demonstrate that there is a difference between gene expression profiles of the uninfected subjects prior to viral exposure and the asymptomatic subjects after the exposure. This observation is important for understanding host immune response and warrants validation in independent gene expression data and/or with RT-qPCR. A more detailed understanding of molecular factors that enable some exposed subjects to avoid infection or remain asymptomatic after the exposure while others demonstrate clinical illness could provide targets for development of more effective vaccines and antiviral treatments. More studies should be undertaken to better characterize the gene expression changes in specific viral infections and patient populations.
Zaas AK, Chen M, Varkey J, Veldman T, Hero AO, Lucas J: Gene expression signatures diagnose influenza and other symptomatic respiratory viral infections in humans. Cell Host Microbe. 2009, 6: 207-217. 10.1016/j.chom.2009.07.006.
Statnikov A, McVoy L, Lytkin N, Aliferis CF: Improving development of the molecular signature for diagnosis of acute respiratory viral infections. Cell Host Microbe. 2010, 7: 100-101. 10.1016/j.chom.2010.01.003.
Ramilo O, Allman W, Chung W, Mejias A, Ardura M, Glaser C: Gene expression patterns in blood leukocytes discriminate patients with acute infections. Blood. 2007, 109: 2066-2077. 10.1182/blood-2006-02-002477.
Aliferis CF, Statnikov A, Tsamardinos I, Mani S, Koutsoukos XD: Local Causal and Markov Blanket Induction for Causal Discovery and Feature Selection for Classification. Part I: Algorithms and Empirical Evaluation. Journal of Machine Learning Research. 2010, 11: 171-234.
Aliferis CF, Statnikov A, Tsamardinos I, Mani S, Koutsoukos XD: Local Causal and Markov Blanket Induction for Causal Discovery and Feature Selection for Classification. Part II: Analysis and Extensions. Journal of Machine Learning Research. 2010, 11: 235-284.
Vapnik VN: Statistical learning theory. 1998, New York: Wiley
Braga-Neto UM, Dougherty ER: Is cross-validation valid for small-sample microarray classification?. Bioinformatics. 2004, 20: 374-380. 10.1093/bioinformatics/btg419.
Aliferis CF, Statnikov A, Tsamardinos I, Schildcrout JS, Shepherd BE, Harrell FE: Factors Influencing the Statistical Power of Complex Data Analysis Protocols for Molecular Signature Development from Microarray Data. PLoS ONE. 2009, 4: e4922-10.1371/journal.pone.0004922.
Kepp O, Senovilla L, Galluzzi L, Panaretakis T, Tesniere A, Schlemmer F: Viral subversion of immunogenic cell death. Cell Cycle. 2009, 8: 860-869.
Schneider RJ, Mohr I: Translation initiation and viral tricks. Trends Biochem Sci. 2003, 28: 130-136. 10.1016/S0968-0004(03)00029-X.
Dresios J, Aschrafi A, Owens GC, Vanderklish PW, Edelman GM, Mauro VP: Cold stress-induced protein Rbm3 binds 60S ribosomal subunits, alters microRNA levels, and enhances global protein synthesis. Proc Natl Acad Sci USA. 2005, 102: 1865-1870. 10.1073/pnas.0409764102.
Wright CF, Oswald BW, Dellis S: Vaccinia virus late transcription is activated in vitro by cellular heterogeneous nuclear ribonucleoproteins. J Biol Chem. 2001, 276: 40680-40686. 10.1074/jbc.M102399200.
Yang Z, Wood C: The transcriptional repressor K-RBP modulates RTA-mediated transactivation and lytic replication of Kaposi's sarcoma-associated herpesvirus. J Virol. 2007, 81: 6294-6306. 10.1128/JVI.02648-06.
Carlson KA, Leisman G, Limoges J, Pohlman GD, Horiba M, Buescher J: Molecular characterization of a putative antiretroviral transcriptional factor, OTK18. J Immunol. 2004, 172: 381-391.
Chen GY, Tang J, Zheng P, Liu Y: CD24 and Siglec-10 selectively repress tissue damage-induced immune responses. Science. 2009, 323: 1722-1725. 10.1126/science.1168988.
Ohman T, Rintahaka J, Kalkkinen N, Matikainen S, Nyman TA: Actin and RIG-I/MAVS signaling components translocate to mitochondria upon influenza A virus infection of human primary macrophages. J Immunol. 2009, 182: 5682-5692. 10.4049/jimmunol.0803093.
Novoa RR, Calderita G, Arranz R, Fontana J, Granzow H, Risco C: Virus factories: associations of cell organelles for viral replication and morphogenesis. Biol Cell. 2005, 97: 147-172. 10.1042/BC20040058.
Grishin A, Li H, Levitan ES, Zaks-Makhina E: Identification of gamma-aminobutyric acid receptor-interacting factor 1 (TRAK2) as a trafficking factor for the K+ channel Kir2.1. J Biol Chem. 2006, 281: 30104-30111. 10.1074/jbc.M602439200.
We thank the authors of  for promptly and generously sharing with us their data, software, and details about their analyses. Alexander Statnikov and Constantin F. Aliferis are acknowledging support from grants R56 LM007948-04A1 from the National Library of Medicine, National Institute of Health and 1UL1RR029893 from the National Center for Research Resources, National Institutes of Health. Jörn-Hendrik Weitkamp was supported by the Eunice Kennedy Shriver National Institute of Child Health and Human Development grant K08HD061607. Lauren McVoy was supported by a fellowship from the Physician Scientist Training Program at New York University School of Medicine.
The authors declare that they have no competing interests.
Conceived and designed the experiments: NIL, AS, CFA. Performed the experiments: NIL, AS. Analyzed the results of experiments: NIL, AS, CFA, JHW, LM. Wrote the paper: NIL, AS, CFA, JHW, LM. All authors read and approved the final manuscript.