- Research article
- Open Access
A chemokine gene expression signature derived from meta-analysis predicts the pathogenicity of viral respiratory infections
© Chang et al; licensee BioMed Central Ltd. 2011
- Received: 14 September 2011
- Accepted: 22 December 2011
- Published: 22 December 2011
During respiratory viral infections host injury occurs due in part to inappropriate host responses. In this study we sought to uncover the host transcriptional responses underlying differences between high- and low-pathogenic infections.
From a compendium of 12 studies that included responses to influenza A subtype H5N1, reconstructed 1918 influenza A virus, and SARS coronavirus, we used meta-analysis to derive multiple gene expression signatures. We compared these signatures by their capacity to segregate biological conditions by pathogenicity and predict pathogenicity in a test data set. The highest-performing signature was expressed as a continuum in low-, medium-, and high-pathogenicity samples, suggesting a direct, analog relationship between expression and pathogenicity. This signature comprised 57 genes including a subnetwork of chemokines, implicating dysregulated cell recruitment in injury.
Highly pathogenic viruses elicit expression of many of the same key genes as lower pathogenic viruses but to a higher degree. This increased degree of expression may result in the uncontrolled co-localization of inflammatory cell types and lead to irreversible host damage.
- Avian Influenza
- Analog Signature
- Respiratory Viral Infection
- Pathogenicity Group
The threat of a highly lethal viral pandemic remains large in the 21st century. In 2003 SARS-coronavirus (CoV) appeared in Asia and then spread globally, causing greater than 40% mortality in individuals over 60 years of age . Since 1997 highly pathogenic avian influenza, influenza A subtype H5N1, has resulted in high mortality rates (between 33% and 100% depending on the population) . Finally in 2009 swine-origin influenza virus A (SOIV) subtype H1N1 emerged in the Americas and led to a pandemic. As Neumann et al. have observed, this virus shares many characteristics with 1918 influenza which resulted in an estimated 50 million deaths . Furthermore, as Ilyushina et al. have shown, SOIV may mutate into more pathogenic forms in as little as ten passages in cell culture .
Injury to the host during respiratory viral infections such as influenza now appears to be the result of inappropriate host responses . Deriving gene expression signatures of high pathogenicity that are robust to biological and experimental variation would be immensely valuable, both in the understanding of pathogenicity as well as in the surveillance of emerging infections. We define a signature as the minimum number of biological variables (here, expressed genes) required to (a) discriminate the phenotype of interest from other phenotypes, (b) identify replicates of the same phenotype, and (c) provide cohesive information about the underlying biological complexity . However, the derivation of gene expression signatures from multiple, independent studies is hindered by a number of factors including the absence of standard experimental protocols, biological variability, and limited sample sizes. Meta-analysis, the analysis of independent but related studies, offers one strategy for overcoming these obstacles.
What host responses are the hallmarks of high pathogenicity? To approach this question we assembled a compendium of host transcriptional responses to low-, medium-, and high-pathogenic infections (which we designate as LPIs, MPIs, and HPIs, respectively) and applied a battery of meta-analysis techniques. Using these techniques we identified signatures of either a digital nature (oppositely expressed in LPIs and HPIs with respect to mock infections) or an analog nature (expressed in a continuum from LPIs to HPIs). We then compared these signatures on the basis of their capacity to predict the outcome of an independent transcriptome experiment. The most accurate signature was analog in nature and implicated excessive chemokine expression and cell recruitment in the development of lethal respiratory infections.
Compendium assembly and signature derivation
Data sets included in the host response compendium
Influenza A H1N1
r1918, 2:6 r1918 (1918 HA, NA in Tx91), 5:3 r1918 (1918 HA, NA, M, NP, S in Tx91), Tx91 (A/Texas/36/91)
1, 3, 5 dpi
0% survival (with r1918, 2:6 r1918, 5:3 r1918), 100% survival (with Tx91)
Kash 2006 ,
Tumpey 2005 (for outcome data)
Influenza A H1N1
1, 3, 5 dpi
0% survival (with r1918), 65% survival (with PR8)
Goodman 2009 
45% survival (in BALB/c), 80% survival (in C57BL/6)
HK213 (A/Hong Kong/213/03)
0% survival (in DBA/2J), 30% survival (in C57BL/6)
Boon 2009 
Influenza A H5N1
HK483 (A/HK/483/97), HK486 (A/HK/486/97), HK486PB2MT (A/HK/486 PB2-627K/97)
2, 4 dpi
0% survival (with HK483, HK486PB2MT), 25% (with HK486)
Fornek 2009 
Influenza A H5N1
VN/1203 (Vietnam/1203/04), r1918 (H1N1)
1, 3, 4 dpi
0% survival (with VN/1203, r1918)
Cilloniz 2010 
vMC0 (UV-inactivated mengovirus)
1, 3 dpi
Janssen 2007 
Minor 2010 
BALB/c (young, aged)
1, 2, 5, 7 dpi
100% survival (in young, aged)
Baas 2008 
icUrbani, icGZ02, icHC/SZ/61/03
BALB/c AnNHsd (young, aged)
12 h, 1, 3 dpi
0% survival (in aged with icHC/SZ/61/03), 40% survival (in aged with icGZ02), 100% survival (in young with all 3 virus strains, in aged with Urbani)
Rockx 2009 
2, 5, 9 dpi
Zornetzer 2010 
Influenza A H1N1
A/CA/04/2009, MA1 (a mouse-adapted strain of A/CA/04/2009)
1, 3, 5 dpi
0% survival (with MA1), 100% survival (with A/CA/04/2009)
Josset in preparation
We then applied three different methods of meta-analysis to derive gene signatures from the compendium. Hierarchical clustering of the log-ratio data showed that two clusters of genes differed in expression with respect to each other and with respect to HPI vs. LPI (Additional File 1, Figure. S1). This suggested that gene signatures could be identified on the basis of opposite directionality with respect to mock infections. To identify these signatures more rigorously we applied two methods of meta-analysis. The first signature comprised 74 genes identified on the basis of statistical tests for differential expression in each biological condition whose results were combined by Fisher's summary statistic (Additional File 2, Figure. S2A; gene list as Additional File 3, Table S1). The second signature comprised 265 genes identified as a module by similar patterns of co-expression and annotation co-membership (Additional File 2, Figure. S2B, S2C; gene list as Additional File 4, Table S2). Due to their opposite directionality in different pathogenicities, we referred to these signatures as digital in nature. A third signature comprised 57 genes and was derived from the direct comparison of expression in HPIs vs. LPIs without reference to mock infections. Additional normalization allowed inter-array comparisons, and a fold-change-based z-test was used to determine differential expression (gene list as Additional File 5, Table S3). Because these genes varied directly between pathogenicities, we referred to this signature as analog in nature. Limited degrees of overlap were observed among the different signatures indicating that the three methods identified different features in the data.
Capacity of analog and digital signatures to classify host transcription
Capacity of the analog signature to predict test set pathogenicity
Characterization of the digital and analog signatures and the role of chemokines
Signature genes may in turn identify mechanisms of pathogenesis. To identify functions and pathways dysregulated across HPIs, we analyzed the annotations associated with each set of signature genes. Annotations for the better-performing digital signature (derived from Fisher's statistic) indicated that a number of pathways may have been altered during HPIs. Genes up-regulated in HPI were associated with inflammation, apoptosis, cell signaling, and cell proliferation, while genes down-regulated in HPI were associated with cytochrome P450-related enzymatic activity. Protein-protein interactions linked many of these genes together indicating possible cooperative effects (Additional File 6, Figure. S3A, S3B). Similar annotations were found to be enriched among the genes of the remaining digital signature (derived by module mapping) including cell differentiation, inflammatory response, and chemical homeostasis.
Closer examination of the analog signature genes indicated a relationship between signature expression and pathogenicity. We determined the proportion of signature genes whose expression was ordered by pathogenicity (from least to greatest in LPI, MPI, and HPI, or vice versa) and compared this proportion to those expected by chance or observed among all genes. Assuming an equal distribution of genes among pathogenicity groups, we expected 33% to display one of these two patterns by chance. By comparison we observed that 43% of all genes in the compendium displayed these patterns (Figure 5B). However, an even greater proportion of analog signature genes displayed these patterns, 81%, consistent with a correlation between expression and pathogenicity (Figure 5C). Using a binomial distribution we calculated the chance enrichment of genes being proportionally expressed with pathogenicity to be p < 0.0019 in the case of Figure 5B (among all compendium genes) and p < 6.6 × 10-23 in the case of Figure 5C (among analog signature genes). Furthermore, if we considered the distribution of gene profiles among all compendium genes to be expected, the enrichment of proportionally expressed genes in the analog signature was statistically significant at p < 1.6 × 10-13.
Analog signature genes expressed in proportion to pathogenicity could in turn be ranked by level of expression in HPIs (Additional File 7, Figure. S4). The most highly expressed signature genes in HPIs were interferon-induced including IFIT3, CXCL10, IIGP1, and IGTP, consistent with the prominence of these genes in functional and network analyses (Figure 5A).
We also examined the temporal expression dynamics of analog signature genes during infection. More detailed examination of the PCC scatter plot (Figure 3D) showed that early time points during HPIs clustered closely with early time points from LPIs and middle time points from MPIs (Figure 5D). These conditions were in turn intermediate to extrema from other HPI and LPI conditions (i.e., early LPI and late HPI conditions). This suggested that LPI, MPI, and HPI may have activated (or suppressed) the same set of genes but that HPI may have done so more rapidly or to higher levels, resulting in levels of expression not reached during LPI or MPI (Figure 5E).
Experiments on the host responses to infection have generated a wealth of data over the past decade. A number of these have focused on respiratory viruses including several which resulted in lethality in experimental animal models. What distinguished these infections from others where experimental animals recovered was the subject of our analysis. To this end we assembled a compendium of transcriptome measurements from the lungs of mice and used meta-analysis to determine gene signatures that distinguished high- from low-pathogenicity infections. To our knowledge ours is the first meta-analysis study to focus on high-pathogenicity infections (HPIs). Previous studies have applied meta-analysis to infection data but focused on more general features of disease. For example, Jenner and Young analyzed the transcriptomes of various host species (mice, macaques, and humans) infected by bacteria or viruses and identified a signature of approximately 500 genes . This signature showed that TLRs and pathogen-mediated signalling were broadly activated during infection but did not identify features resulting in animal lethality. Likewise Pennings et al. examined the transcriptional profiles of lung inflammation (in mice and macaques) due to various factors (including viruses, bacteria, chemicals, and allergens) and derived a 383-gene signature up-regulated during inflammation . This signature focused on the interferon response and immune signalling, but likewise no correlate of pathogenicity was identified. In addition both studies relied mainly on hierarchical clustering to determine signatures of interest. While hierarchical clustering provides a useful overview of high-dimensional data, it has the disadvantage of lacking a strong statistical basis (i.e., hypothesis testing) and its results (gene clusters) can be difficult to relate to specific outcomes.
In contrast to previous studies, our goal was to determine a gene signature specific for HPIs. We therefore applied an ensemble of meta-analysis methods and derived multiple signatures which we compared by multiple criteria including the capacity to predict the outcome of a test data set. Each method had a different statistical rationale and could be expected to identify different features of the compendium. We had no a priori expectation of which signature would produce the best outcome in terms of these criteria.
The 74-gene signature derived using Fisher's summary statistic showed a modest capacity to separate samples into pathogenicity groups and to predict test set pathogenicity. A number of the genes in this signature have known connections to the immune response and the outcome of respiratory infections, e.g., genes for the chemokines osteopontin (SPP1) and RANTES (CCL5) and the chemokine receptor CCR2. In previous studies mice lacking SPP1 and CCL5 were found to clear influenza infection with no adverse effects [10, 11], while mice lacking CCR2 survived infection by a mouse-adapted influenza A virus that killed wild-type mice . These studies suggest that at least some of the genes up-regulated during HPIs are non-essential to resolving influenza infection and that dysregulated activation may even be detrimental to the host.
The best performing signature, however, comprised 57 genes derived using the fold change-based z-test. Strikingly the majority of the genes in this signature were expressed at levels that corresponded with pathogenicity. For the majority of these genes, expression was lowest in LPI, highest in HPI, and intermediate in MPI (which had not been used in the derivation of the signature). In this case the signature appeared to provide a continuously variable signal that matched with output, a characteristic of analog signals. This finding also suggests that high- and low-pathogenic infections may result in the expression of the same key genes but with different kinetics. In particular HPIs may result in increased expression of signature gene products beyond the capacity of the host to cope, resulting in irreversible damage.
Genes of the analog signature largely differed from those in the digital signatures but overlapped at the pathway level. For example, chemokine genes were present in both Fisher's statistic-derived and analog signatures. However, the analog signature displayed additional coherence, encoding multiple chemokines for the same receptor, namely MIG (CXCL9), IP-10 (CXCL10), and I-TAC (CXCL11) which all serve as ligands for the receptor CXCR3 . CXCR3 is expressed on the surface of Th1 cells as well as NK and NKT cells and regulates the migration of these cells to sites of infection. Recent evidence indicates that CXCR3 engagement may drive further recruitment and inflammation, resulting in a positive feedback loop that may contribute to pathogenicity . The analog signature also included genes for chemokines that bind the CCR3 receptor, specifically MCP-3 (CCL7) and MCP-4 (CCL13). CCR3 is the major receptor expressed on eosinophils and has previously been shown to have a role in the promotion of lung inflammation . Interestingly CXCR3 ligands have been postulated to be antagonists for CCR3 , and the expression of both sets of chemokines may reflect a high degree of dysregulation during HPIs and the recruitment of multiple immune cell types that may not normally co-localize during a controlled infection.
In addition our signatures also identified genes that were down-regulated during HPI, relative to either mock infection or LPI. For example the Fisher's statistic-derived signature included several genes whose products may normally help to resolve infection. Hepsin (HPN) has been found to cleave influenza hemagglutinin directly resulting in non-infectious virus particles . Likewise surfactant-associated protein A1 (SFTPA1) and surfactant-associated protein D (SFTPD) maintain pulmonary structure but have also been found to inhibit infectivity by mechanisms that remain to be elucidated [17, 18]. The analog signature contained a similar set of genes but expressed at lower levels in HPIs compared to LPIs rather than to mock infections. Among these genes were several members of the secretoglobin family including SCGB1A1, SCGB3A1, and SCGB3A3. Interestingly many secretoglobins are expressed specifically in the lung epithelium and may contribute to lung repair following damage [19, 20]. For example, uteroglobin/CC16 (SCGB1A1) is secreted by bronchiolar Clara cells and postulated to have a role in reducing airway inflammation, though its exact function remains to be elucidated .
Together these results implicate dysregulated cell recruitment and inflammation (up-regulated in HPI) and impaired lung protection (down-regulated in HPI) in the events leading to lethality in mice. The different meta-analysis methods all aimed to identify pathology-relevant gene expression, and concordance at the pathway level offered a degree of cross-validation. However, our results also demonstrate that the methods identified different features of the data and argue for the application of multiple methods in future studies when the most predictive features of the data are not known in advance.
Finally our finding that pathogenicity corresponds to the expression levels of a defined set of genes may help to inform future therapies. Altering the outcome of infection may not require opposite regulation to be achieved in particular pathways; instead, adequate tempering or delay of those pathways may be sufficient. For example a large number of potential CXCR3 antagonists have been developed recently in the context of other inflammatory diseases . These may provide the basis for novel therapeutic strategies that reduce, but need not completely eliminate, receptor activation to effect control of highly pathogenic respiratory infections.
A literature search for microarray data pertaining to respiratory viral infections in mouse models was performed, and candidate data sets were identified. Data were included in the compendium if they were derived from experiments on wild-type mice, as this ensured normal development. We also required that raw data be measured on whole-genome array platforms. Data were sourced from either in-house or external studies (Table 1). In the latter case, data were downloaded from public repositories or obtained directly from study authors. Pathogenicities of the corresponding experiments were determined by examination of the published outcome data (e.g., mortality or morbidity data or in cases of LPIs, statements attesting to the survival or recovery of all animal subjects). We defined HPIs as experiments resulting in 0% survival, LPIs as experiments resulting in 100% survival, and MPIs as experiments resulting in intermediate levels of survival (between 0% and 100%, non-inclusive). A biological condition was defined as any unique study/time point/virus strain/mouse strain combination.
Data pre-processing and quality control
Data from Agilent platforms were obtained from Agilent Feature Extraction software and pre-processed in GeneSpring GX 11 (Agilent Technologies, Santa Clara CA) using median-background subtraction on a per-study basis. Data from the Affymetrix platform were pre-processed in GeneSpring GX 11 using the MAS5 algorithm. Quality control was restricted to probe-level filtering. Probes with present flags in all of the technical replicates for at least one biological condition in a given study were retained in the data set for that study. For ratio-based meta-analysis, two-channel data were reformatted as log2-ratios (infected/mock-infected), and one-channel data were converted to a two-channel-like format where mock-infection intensities were averaged prior to ratio. For meta-analysis based on direct comparison of intensities, two-channel data were separated into individual channels, and inter-array normalization was performed using NeONORM . Additional data preparation involved converting probe- to gene-level data. Measurements for probes mapping to multiple genes were duplicated for each gene, and probes mapping to the same gene were averaged. RefSeq IDs were chosen as the common gene identifier. Probes not mapping to RefSeq IDs were excluded. Biological and technical replicates were averaged for the intensity-based approaches. Log-ratio measurements from all studies were combined resulting in measurements for 27,567 genes by 397 biological conditions. All data-handling tasks were performed using custom scripts in Perl.
Fisher's summary-statistic method
We applied the Fisher's summary-statistic method based on the previous application of Rhodes et al. . Our implementation was derived from MADAM (Meta-Analysis Data Aggregation Methods), a toolbox for the R statistical language . Statistical tests were performed on a per-gene basis for all arrays associated with a given biological condition. Specifically log-ratio measurements for each array were re-normalized by dividing by the standard deviation of the array, and one-sided t-tests were used to identify genes up- and down-regulated relative to mock infection, i.e., for the null hypotheses H0: log2 (infected/mock) ≤ 0 and H0: log2 (infected/mock) ≥ 0, respectively. p-values for each directionality (up- or down-regulation) and pathogenicity (HPI or LPI) were combined according to Fisher's summary statistic S = -2 ∑ log p. The significance of S was determined by a permutation test whereby p-values were shuffled 1000 times among biological conditions, S was re-derived each time, and the resulting S values constituted the null distribution. An FDR-corrected p-value of 0.05 was used as a threshold for significance. Genes that were found to be both up- and down-regulated for a given pathogenicity (HPI or LPI) were excluded from further analysis. The signature comprised genes found to be oppositely regulated in HPI and LPI.
Module map method
The module-mapping method of Segal et al. identifies modules of co-expressed genes from microarray data based on co-expression within the compendium and co-membership in previously annotated gene sets . We used the module-map method as implemented in the software package Genomica (available at http://genomica.weizmann.ac.il) using the log-ratio compendium and default settings. Modules were tested for associations with HPI and LPI outcomes, and the signature comprised the largest module found to be oppositely regulated in HPI and LPI.
Fold change-based z-test
Intensity-based measurements were used to test for differential expression in all HPIs vs. all LPIs. Specifically a fold-change-based z-test was used to compute the statistical significance of log2 (HPI/LPI) values according to intra- and inter-group variability. Statistically significant log2 (HPI/LPI) values were then ranked by size, and genes with log2 (HPI/LPI) < -1.6 or > 2 were retained. The cut-offs were selected based on a histogram of fold-changes over all statistically significantly differentially expressed genes.
Functional and network analysis
We used NIH DAVID (Database for Visualization and Integrated Discovery) and Ingenuity Pathways Analysis (Ingenuity Systems, Redwood City CA) for further analysis of meta-analysis results. DAVID identifies GO terms and other functional annotations enriched in gene lists . An enrichment score is calculated for clusters of similar annotations as -log10 of the geometric mean of the hypergeometric test p-values within the cluster . We assessed enrichment of the GO categories GOTERM_BP_ALL, GOTERM_CC_ALL, and COTERM_MF_ALL and pathway categories BIOCARTA and KEGG_PATHWAY; otherwise, default options in DAVID were used. We identified significant annotations by enrichment scores > 1.30 representing a mean cluster p < 0.05.
We thank Karl Kugler, Richard Green, and Chester Ni for technical assistance and Louis Rosenthal, Richard Webby, Jacco Boon, Linda Toth, and Lu Lu for generously sharing their data. We also thank Lynn Law, Sean Proll, and Marcus Korth for valuable feedback on the manuscript. This project has been funded with Federal funds from the National Institute of Allergy and Infectious Diseases, National Institutes of Health, Department of Health and Human Services, under Contract No. HHSN272200800060C. NT was funded by a grant from the Genopole Evry.
- Donnelly CA, Ghani AC, Leung GM, Hedley AJ, Fraser C, Riley S, Abu-Raddad LJ, Ho L-M, Thach T-Q, Chau P, Chan K-P, Lam T-H, Tse L-Y, Tsang T, Liu S-H, Kong JHB, Lau EMC, Ferguson NM, Anderson RM: Epidemiological determinants of spread of causal agent of severe acute respiratory syndrome in Hong Kong. Lancet. 2003, 361: 1761-1766. 10.1016/S0140-6736(03)13410-1.View ArticlePubMedGoogle Scholar
- Beigel JH, Farrar J, Han AM, Hayden FG, Hyer R, de Jong MD, Lochindarat S, Nguyen TKT, Nguyen TH, Tran TH, Nicoll A, Touch S, Yuen K-Y: Avian influenza A (H5N1) infection in humans. N Engl J Med. 2005, 353: 1374-1385.View ArticlePubMedGoogle Scholar
- Neumann G, Noda T, Kawaoka Y: Emergence and pandemic potential of swine-origin H1N1 influenza virus. Nature. 2009, 459: 931-939. 10.1038/nature08157.PubMed CentralView ArticlePubMedGoogle Scholar
- Ilyushina NA, Ducatez MF, Rehg JE, Marathe BM, Marjuki H, Bovin NV, Webster RG, Webby RJ: Does pandemic A/H1N1 virus have the potential to become more pathogenic?. MBio. 2010, 1:Google Scholar
- La Gruta NL, Kedzierska K, Stambas J, Doherty PC: A question of self-preservation: immunopathology in influenza virus infection. Immunol Cell Biol. 2007, 85: 85-92. 10.1038/sj.icb.7100026.View ArticlePubMedGoogle Scholar
- Gregorius H-R: The isolation principle of clustering: structural characteristics and implementation. Acta Biotheor. 2006, 54: 219-233. 10.1007/s10441-006-8255-3.View ArticlePubMedGoogle Scholar
- Bécavin C, Benecke A: New dimensionality reduction methods for the representation of high dimensional "omics" data. Expert Rev Mol Diagn. 2011, 11: 27-34. 10.1586/erm.10.95.View ArticlePubMedGoogle Scholar
- Jenner RG, Young RA: Insights into host responses against pathogens from transcriptional profiling. Nat Rev Microbiol. 2005, 3: 281-294. 10.1038/nrmicro1126.View ArticlePubMedGoogle Scholar
- Pennings JLA, Kimman TG, Janssen R: Identification of a common gene expression response in different lung inflammatory diseases in rodents and macaques. PLoS ONE. 2008, 3: e2596-10.1371/journal.pone.0002596.PubMed CentralView ArticlePubMedGoogle Scholar
- Abel B, Freigang S, Bachmann MF, Boschert U, Kopf M: Osteopontin is not required for the development of Th1 responses and viral immunity. J Immunol. 2005, 175: 6006-6013.View ArticlePubMedGoogle Scholar
- Wareing MD, Lyon AB, Lu B, Gerard C, Sarawar SR: Chemokine expression during the development and resolution of a pulmonary leukocyte response to influenza A virus infection in mice. J Leukoc Biol. 2004, 76: 886-895. 10.1189/jlb.1203644.View ArticlePubMedGoogle Scholar
- Dawson TC, Beck MA, Kuziel WA, Henderson F, Maeda N: Contrasting effects of CCR5 and CCR2 deficiency in the pulmonary inflammatory response to influenza A virus. Am J Pathol. 2000, 156: 1951-1959. 10.1016/S0002-9440(10)65068-7.PubMed CentralView ArticlePubMedGoogle Scholar
- Groom JR, Luster AD: CXCR3 ligands: redundant, collaborative and antagonistic functions. Immunol Cell Biol. 2011, 89: 207-215. 10.1038/icb.2010.158.View ArticlePubMedGoogle Scholar
- Daugherty BL, Siciliano SJ, DeMartino JA, Malkowitz L, Sirotina A, Springer MS: Cloning, expression, and characterization of the human eosinophil eotaxin receptor. J Exp Med. 1996, 183: 2349-2354. 10.1084/jem.183.5.2349.View ArticlePubMedGoogle Scholar
- Loetscher P, Pellegrino A, Gong JH, Mattioli I, Loetscher M, Bardi G, Baggiolini M, Clark-Lewis I: The ligands of CXC chemokine receptor 3, I-TAC, Mig, and IP10, are natural antagonists for CCR3. J Biol Chem. 2001, 276: 2986-2991. 10.1074/jbc.M005652200.View ArticlePubMedGoogle Scholar
- Bertram S, Glowacka I, Blazejewska P, Soilleux E, Allen P, Danisch S, Steffen I, Choi S-Y, Park Y, Schneider H, Schughart K, Pöhlmann S: TMPRSS2 and TMPRSS4 facilitate trypsin-independent spread of influenza virus in Caco-2 cells. J Virol. 2010, 84: 10016-10025. 10.1128/JVI.00239-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Hartshorn KL, White MR, Shepherd V, Reid K, Jensenius JC, Crouch EC: Mechanisms of anti-influenza activity of surfactant proteins A and D: comparison with serum collectins. Am J Physiol. 1997, 273: L1156-1166.PubMedGoogle Scholar
- Mikerov AN, White M, Hartshorn K, Wang G, Floros J: Inhibition of hemagglutination activity of influenza A viruses by SP-A1 and SP-A2 variants expressed in CHO cells. Med Microbiol Immunol. 2008, 197: 9-12.PubMed CentralView ArticlePubMedGoogle Scholar
- Nord M, Cassel TN, Braun H, Suske G: Regulation of the Clara cell secretory protein/uteroglobin promoter in lung. Ann N Y Acad Sci. 2000, 923: 154-165.View ArticlePubMedGoogle Scholar
- Shijubo N, Kawabata I, Sato N, Itoh Y: Clinical aspects of Clara cell 10-kDa protein/uteroglobin (secretoglobin 1A1). Curr Pharm Des. 2003, 9: 1139-1149. 10.2174/1381612033455026.View ArticlePubMedGoogle Scholar
- Wong AP, Keating A, Waddell TK: Airway regeneration: the role of the Clara cell secretory protein and the cells that express it. Cytotherapy. 2009, 11: 676-687. 10.3109/14653240903313974.View ArticlePubMedGoogle Scholar
- Wijtmans M, Verzijl D, Leurs R, de Esch IJP, Smit MJ: Towards small-molecule CXCR3 ligands with clinical potential. ChemMedChem. 2008, 3: 861-872. 10.1002/cmdc.200700365.View ArticlePubMedGoogle Scholar
- Noth S, Brysbaert G, Benecke A: Normalization using weighted negative second order exponential error functions (NeONORM) provides robustness against asymmetries in comparative transcriptome profiles and avoids false calls. Genomics Proteomics Bioinformatics. 2006, 4: 90-109. 10.1016/S1672-0229(06)60021-1.View ArticlePubMedGoogle Scholar
- Rhodes DR, Barrette TR, Rubin MA, Ghosh D, Chinnaiyan AM: Meta-analysis of microarrays: interstudy validation of gene expression profiles reveals pathway dysregulation in prostate cancer. Cancer Res. 2002, 62: 4427-4433.PubMedGoogle Scholar
- Kugler KG, Mueller LA, Graber A: MADAM - An open source meta-analysis toolbox for R and Bioconductor. Source Code Biol Med. 2010, 5: 3-10.1186/1751-0473-5-3.PubMed CentralView ArticlePubMedGoogle Scholar
- Fisher RA: Statistical methods for research workers. 1932, Oliver and BoydGoogle Scholar
- Segal E, Shapira M, Regev A, Pe'er D, Botstein D, Koller D, Friedman N: Module networks: identifying regulatory modules and their condition-specific regulators from gene expression data. Nat Genet. 2003, 34: 166-176. 10.1038/ng1165.View ArticlePubMedGoogle Scholar
- Hosack DA, Dennis G, Sherman BT, Lane HC, Lempicki RA: Identifying biological themes within lists of genes with EASE. Genome Biol. 2003, 4: R70-10.1186/gb-2003-4-10-r70.PubMed CentralView ArticlePubMedGoogle Scholar
- Huang DW, Sherman BT, Tan Q, Kir J, Liu D, Bryant D, Guo Y, Stephens R, Baseler MW, Lane HC, Lempicki RA: DAVID Bioinformatics Resources: expanded annotation database and novel algorithms to better extract biology from large gene lists. Nucleic Acids Res. 2007, 35: W169-175. 10.1093/nar/gkm415.PubMed CentralView ArticlePubMedGoogle Scholar
- Kash JC, Tumpey TM, Proll SC, Carter V, Perwitasari O, Thomas MJ, Basler CF, Palese P, Taubenberger JK, García-Sastre A, Swayne DE, Katze MG: Genomic analysis of increased host immune and cell death responses induced by 1918 influenza virus. Nature. 2006, 443: 578-581.PubMed CentralPubMedGoogle Scholar
- Tumpey TM, Basler CF, Aguilar PV, Zeng H, Solórzano A, Swayne DE, Cox NJ, Katz JM, Taubenberger JK, Palese P, García-Sastre A: Characterization of the reconstructed 1918 Spanish influenza pandemic virus. Science. 2005, 310: 77-80. 10.1126/science.1119392.View ArticlePubMedGoogle Scholar
- Goodman AG, Fornek JL, Medigeshi GR, Perrone LA, Peng X, Dyer MD, Proll SC, Knoblaugh SE, Carter VS, Korth MJ, Nelson JA, Tumpey TM, Katze MG: P58(IPK): a novel "CIHD" member of the host innate defense response against pathogenic virus infection. PLoS Pathog. 2009, 5: e1000438-10.1371/journal.ppat.1000438.PubMed CentralView ArticlePubMedGoogle Scholar
- Ding M, Lu L, Toth LA: Gene expression in lung and basal forebrain during influenza infection in mice. Genes Brain Behav. 2008, 7: 173-183. 10.1111/j.1601-183X.2007.00335.x.View ArticlePubMedGoogle Scholar
- Toth LA, Rehg JE, Webster RG: Strain differences in sleep and other pathophysiological sequelae of influenza virus infection in naive and immunized mice. J Neuroimmunol. 1995, 58: 89-99. 10.1016/0165-5728(94)00193-R.View ArticlePubMedGoogle Scholar
- Boon ACM, deBeauchamp J, Hollmann A, Luke J, Kotb M, Rowe S, Finkelstein D, Neale G, Lu L, Williams RW, Webby RJ: Host genetic variation affects resistance to infection with a highly pathogenic H5N1 influenza A virus in mice. J Virol. 2009, 83: 10417-10426. 10.1128/JVI.00514-09.PubMed CentralView ArticlePubMedGoogle Scholar
- Fornek JL, Gillim-Ross L, Santos C, Carter V, Ward JM, Cheng LI, Proll S, Katze MG, Subbarao K: A single-amino-acid substitution in a polymerase protein of an H5N1 influenza virus is associated with systemic infection and impaired T-cell activation in mice. J Virol. 2009, 83: 11102-11115. 10.1128/JVI.00994-09.PubMed CentralView ArticlePubMedGoogle Scholar
- Cilloniz C, Pantin-Jackwood MJ, Ni C, Goodman AG, Peng X, Proll SC, Carter VS, Rosenzweig ER, Szretter KJ, Katz JM, Korth MJ, Swayne DE, Tumpey TM, Katze MG: Lethal dissemination of H5N1 influenza virus is associated with dysregulation of inflammation and lipoxin signaling in a mouse model of infection. J Virol. 2010, 84: 7613-7624. 10.1128/JVI.00553-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Janssen R, Pennings J, Hodemaekers H, Buisman A, van Oosten M, de Rond L, Oztürk K, Dormans J, Kimman T, Hoebee B: Host transcription profiles upon primary respiratory syncytial virus infection. J Virol. 2007, 81: 5958-5967. 10.1128/JVI.02220-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Minor RAC, Limmon GV, Miller-DeGraff L, Dixon D, Andrews DMK, Kaufman RJ, Imani F: Double-stranded RNA-activated protein kinase regulates early innate immune responses during respiratory syncytial virus infection. J Interferon Cytokine Res. 2010, 30: 263-272. 10.1089/jir.2009.0051.View ArticlePubMedGoogle Scholar
- Baas T, Roberts A, Teal TH, Vogel L, Chen J, Tumpey TM, Katze MG, Subbarao K: Genomic analysis reveals age-dependent innate immune responses to severe acute respiratory syndrome coronavirus. J Virol. 2008, 82: 9465-9476. 10.1128/JVI.00489-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Rockx B, Baas T, Zornetzer GA, Haagmans B, Sheahan T, Frieman M, Dyer MD, Teal TH, Proll S, van den Brand J, Baric R, Katze MG: Early upregulation of acute respiratory distress syndrome-associated cytokines promotes lethal disease in an aged-mouse model of severe acute respiratory syndrome coronavirus infection. J Virol. 2009, 83: 7062-7074. 10.1128/JVI.00127-09.PubMed CentralView ArticlePubMedGoogle Scholar
- Zornetzer GA, Frieman MB, Rosenzweig E, Korth MJ, Page C, Baric RS, Katze MG: Transcriptomic analysis reveals a mechanism for a prefibrotic phenotype in STAT1 knockout mice during severe acute respiratory syndrome coronavirus infection. J Virol. 2010, 84: 11297-11309. 10.1128/JVI.01130-10.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.