key: cord-295781-b831y105
authors: VanLeuven, James T.; Ridenhour, Benjamin J.; Gonzalez, Andres J.; Miller, Craig R.; Miura, Tanya A.
title: Lung epithelial cells have virus-specific and shared gene expression responses to infection by diverse respiratory viruses
date: 2017-06-02
journal: PLoS One
DOI: 10.1371/journal.pone.0178408
sha: 
doc_id: 295781
cord_uid: b831y105

The severity of respiratory viral infections is partially determined by the cellular response mounted by infected lung epithelial cells. Disease prevention and treatment is dependent on our understanding of the shared and unique responses elicited by diverse viruses, yet few studies compare host responses to viruses from different families while controlling other experimental parameters. Murine models are commonly used to study the pathogenesis of respiratory viral infections, and in vitro studies using murine cells provide mechanistic insight into the pathogenesis observed in vivo. We used microarray analysis to compare changes in gene expression of murine lung epithelial cells infected individually by three respiratory viruses causing mild (rhinovirus, RV1B), moderate (coronavirus, MHV-1), and severe (influenza A virus, PR8) disease in mice. RV1B infection caused numerous gene expression changes, but the differential effect peaked at 12 hours post-infection. PR8 altered an intermediate number of genes whose expression continued to change through 24 hours. MHV-1 had comparatively few effects on host gene expression. The viruses elicited highly overlapping responses in antiviral genes, though MHV-1 induced a lower type I interferon response than the other two viruses. Signature genes were identified for each virus and included host defense genes for PR8, tissue remodeling genes for RV1B, and transcription factors for MHV-1. Our comparative approach identified universal and specific transcriptional signatures of virus infection that can be used to distinguish shared and virus-specific mechanisms of pathogenesis in the respiratory tract.

Viruses from several different families, including Picornaviridae, Orthomyxoviridae, Paramyxoviridae, Coronaviridae, and Adenoviridae, infect and cause diseases in the respiratory tract. These diseases range from mild infections of the upper respiratory tract, to severe diseases when viruses infect the lungs. Respiratory viruses commonly target epithelial cells of the a1111111111 a1111111111 a1111111111 a1111111111 a1111111111 airways and lungs. These epithelial cells are responsible for detecting viral pathogens and initiating antiviral responses at the level of infected cells and the immune system, and therefore their response to infection has an important role in determining disease outcomes. Knowledge of the shared and virus-specific responses of respiratory epithelial cells to infection by diverse viruses is fundamental to understanding viral pathogenesis and developing therapies to treat severe respiratory infections.

Murine models of respiratory viral infections have been widely used to identify the mechanisms that determine disease severity in the respiratory tract. While these models are invaluable for evaluating pathology and host responses to infection, parallel in vitro studies can be used to identify gene expression and signaling pathway changes that occur in infected cells to mediate pathogenesis. In this study, we compare the gene expression response of murine lung epithelial cells to infection by three respiratory viruses used in murine models: rhinovirus (serotype RV1B) in the family Picornaviridae, mouse hepatitis virus (MHV strain 1) in the family Coronaviridae, and influenza A virus (strain PR8) in the family Orthomyxoviridae. Viruses from different families and with different replication strategies were chosen to identify which responses to infection in respiratory epithelial cells are shared and which are virus-specific. In the following paragraphs, we describe some of the key features of these three viruses in murine models.

Minor group rhinoviruses, including RV1B, use low-density lipoprotein receptor to enter either human or murine cells. Because RV1B can replicate in mouse cells, it has been used to study immune responses and mechanisms of asthma exacerbation in infected mice [1] [2] [3] [4] [5] . Intranasal inoculation of mice with a high dose of RV1B results in viral replication and an early inflammatory response in the respiratory tract without clinical signs of disease [1, [3] [4] [5] . RV1B antigens have been detected in cells of the epithelia and lamina propria of the upper respiratory tract of infected mice [2, 4] .

MHV-1 is a strain of mouse hepatitis virus that preferentially replicates and causes disease in the respiratory tract of specific mouse strains and thus serves as a model for respiratory coronaviral infections [6, 7] . MHV-1 causes severe disease in A/J-strain mice and milder pathology in other mouse strains [6, 8] . Mouse strain-dependent disease severity corresponds to inflammatory responses, fibrin deposition, and reduced type I interferon (IFN) production [6] . Further, type I IFN-mediated signaling is required for protection from severe disease during MHV-1 infection of resistant mice [8] . MHV-1 has been detected in pulmonary macrophages of infected mice, but has not been reported to infect respiratory epithelial cells in vivo [6] . Although primary murine alveolar epithelial cells are susceptible to MHV-1 infection in vitro, their role during in vivo infection is not clear [9] .

Mice have been used for decades to study the pathogenesis of influenza. One of the most commonly used strains of influenza A virus, PR8, has been serially passaged in mice to produce a model for pulmonary infection. PR8 infection results in a range of disease severities that is mouse strain-dependent [10] . Although susceptible mice mount a type I IFN response to PR8 infection, lethal infection is associated with spread of virus to the alveoli and an excessive inflammatory response [10] [11] [12] [13] . PR8 replicates in bronchiolar and alveolar epithelial cells of the lower respiratory tract in vivo and in primary murine respiratory epithelial cells in vitro [9, 12, 14, 15] .

In this study, we used a murine lung epithelial cell line (LA4) to compare the gene expression response to these three unrelated viruses. LA4 cells were derived from neoplastic lung epithelia from strain A (A/He) mice and have some properties of alveolar type II cells, although they do not maintain a highly differentiated type II phenotype [16] . Strain A (A/J) mice are highly susceptible to respiratory viral infections, including MHV-1 and influenza A viruses [6, 10] . Other studies have demonstrated that LA4 cells are susceptible to infection by PR8 and RV1B [15, 17] . In this study, we show that LA4 cells are also susceptible to infection by MHV-1 (hereafter referred to as MHV). The gene expression response of LA4 cells to infection by MHV, PR8, and RV1B (hereafter referred to as RV) differed in timing and magnitude of the changes. While we expected to see highly divergent transcription responses to these three viruses, they induced expression of a surprisingly large overlapping set of shared genes, including genes involved in antiviral responses. However, and more in line with our expectation, each virus also altered expression of unique genes, which highlight their different replication strategies and mechanisms of pathogenesis in murine hosts.

Virus stocks were generated by 24 h of growth from a low dose inoculum in the cell lines described below. Supernatant medium from infected cells was centrifuged to remove cell debris, aliquoted, flash frozen and stored at -80˚C. PR8 (A/Puerto Rico/8/1934 (H1N1)), obtained from BEI Resources (NR-3169), was grown and titrated by plaque assay in MDCK cells from ATCC (CCL-34) [18] . MHV, obtained from ATCC (VR-261), was grown and titrated by plaque assay in 17Cl.1 cells [19] (provided by Dr. Kathryn Holmes, University of Colorado Denver School of Medicine). RV, obtained from ATCC (VR-1645), was grown and titrated by tissue culture infectious dose 50% (TCID 50 ) assay in HeLa (ATCC: CCL-2) cells [17] . LA4, a murine lung epithelial cell line from ATCC (CCl-196), was cultured in Ham's F12K medium (Mediatech) with 10% FBS (Atlanta Biologicals) and 1X antibiotic/antimycotic (Gibco).

Our experimental approach was to inoculate LA4 cells with the three viruses at times t = 0 h and t = 12 h and harvest RNA for microarray analysis at t = 24 h. Controls were mock-inoculated at both time points. Preliminary experiments were done to establish the growth kinetics of each virus and determine a multiplicity of infection (MOI) that resulted in comparable numbers of cells positive for viral antigen at 24 h post-infection (S1 Fig) . Based on this, LA4 cells were inoculated with 3 TCID 50 /cell RV, 1 PFU/cell PR8, or 3 PFU/cell MHV. Triplicate wells of LA4 cells in 6-well plates were inoculated with each virus diluted in serum-free medium or were mock-inoculated with serum-free medium for 1 h at 37˚C. Viral inocula were removed and the cells were rinsed twice with serum-free medium. The cells were incubated in Ham's F12K medium with 2% FBS until the 24 h time point at which time RNA was isolated from cell cultures using an RNAeasy Plus kit (QIAGEN) and transcript levels were measured by microarray (NimbleGen Mus musculus MM9 Expression Array). For the 24 h samples, the media were removed and replaced with fresh media 12 h after inoculation, which is the same time that 12 h samples were inoculated. Raw and processed data are available from NCBI Gene Expression Omnibus under accession number GSE89190.

NimbleScan v2.5 software (NibleGen, Madison, WI) was used to extract raw intensity data for each probe on each array. Intensity data were read into the R statistical computing environment and checked for quality [20] . Data were prepared for processing using the pdInfoBuilder package and then normalized using the robust multichip average (RMA) function in the oligo package [21] .

Statistical tests for differences in expression between treatments were conducted on the normalized expression data using a linear mixed-effect model followed by linear contrasts corrected for multiple comparisons. More specifically, expression was modeled as a function of treatment, while probes for a particular gene were treated as random effects. The models used the nlme::lme function in R. The data contained seven treatments: three viruses tested at two time points (12 h and 24 h) each plus the mock treatment (RV 12 , RV 24 , MHV 12 , MHV 24 , PR8 12 , PR8 24 , and mock). The following nine post hoc, two-sided contrasts were then performed on the fitted models using the multcomp::glht function in R: each virus-time combination vs. mock (RV 12 vs. mock, RV 24 vs. mock, MHV 12 vs. mock, MHV 24 vs. mock, PR8 12 vs. mock, PR8 24 vs. mock) and each pairwise combination of viruses at the 24 h time point (RV 24 vs. MHV 24 , RV 24 vs. PR8 24 MHV 24 vs. PR8 24 ). These 9 tests had their p-values adjusted by the multcomp::summary.glht function according to their joint distribution. Any factors detected to be significant at the family (gene) level were then subsequently corrected using the Benjamini-Hochberg algorithm [22] with a false discovery rate set at 1%.

Genes associated with type I IFN responses were identified among the sets of genes with differential expression for each virus compared to mock at 12 h and 24 h using the Interferome v.2.01 database (http://interferome.its.monash.edu.au; [23] ). This database was queried using the search criteria: Type I IFN (all), in vitro, Mus musculus, 2.0 fold change (up or down). Interferome genes were manually sorted into functional categories: antiviral, IFN signaling, viral detection, immune response, MHC class I, inhibitory, apoptosis, ubiquitination, miscellaneous, and unknown. The significance of each virus having genes in the specific categories was tested using a chi-squared test.

Gene expression responses to RV1B were compared between our data from mouse cells and published data using human cells [24] using the MGI vertebrate homology database provided by The Jackson Laboratory [25] as well as the annotate package in R.

MHV, PR8, and RV alter cellular gene expression by different magnitudes and with different timing Our goal was to evaluate the gene expression response of murine respiratory epithelial cells to infection by three unrelated respiratory viruses studied in murine models. A preliminary study was undertaken to identify a murine cell line that was susceptible to infection by the three viruses. LA4 cells were chosen because MHV, RV, and PR8 established productive infections in this line, as shown by increased viral titers released from infected cultures over a 24 h infection (S1 Fig). Additionally, comparable numbers of viral antigen-positive cells were observed for the three viruses 24 h after inoculation of LA4 cells (S1 Fig) .

To compare how unrelated respiratory viruses (MHV, PR8, and RV) alter gene expression of murine epithelial cells, we inoculated LA4 cells with each of the three viruses and evaluated cellular gene expression by microarray analysis after 12 and 24 h of infection compared to mock-inoculated controls. Fig 1 shows the log 2 -fold change in expression level of genes that were differentially expressed in virus-infected, compared to mock-inoculated, cells. By plotting the changes in gene expression at 12 vs. 24 hours, we observed differences in magnitude and timing of gene expression changes mediated by the three viruses. The genes with significantly different expression in MHV-infected cells had low fold change values (Fig 1A) . At 24 h, when gene expression changes were the highest, genes that were up-regulated by MHV infection had log 2 -fold change values of less than five. In contrast, PR8 and RV induced expression of many genes by greater than five log 2 -fold at 24 h, and genes were spread consistently across the full range of values. By 24 h, the genes most strongly up-regulated by PR8 and RV induced changes of 7-9.5 log 2 -fold and 6-7.5 log 2 -fold compared to mock, respectively. This same pattern was observed with the down-regulated genes (Fig 1) .

In addition to the differences in fold change values, the three viruses also differed in the timing of gene expression changes. MHV altered expression of relatively few host genes, most of which were only significantly different from mock at 24 h ( Fig 1A) . While both PR8 and RV induced expression of large subsets of host genes, they did so with different timing. PR8-induced changes in gene expression occurred at a constant rate: the expression level of most genes at 24 h was approximately twice the expression value at 12 h ( Fig 1B) . In contrast, RV infection altered expression of a large number of genes by 12 h and the expression levels were maintained at approximately the same levels at the 24 h time point (Fig 1C) .

Taken together, we observed differences in magnitude and timing of gene expression changes mediated by the three viruses. The limited response to MHV infection is in agreement with studies of other coronaviruses, such as MHV-A59 [26] and SARS-CoV [27, 28] . In addition to inducing minor transcriptional up-regulation of host genes, MHV-A59 shuts down host gene expression by enhancing mRNA degradation [26] . A related coronavirus, SARS-CoV, also induces degradation of host mRNAs [29] . The low numbers of host mRNAs that were altered in response to MHV infection in our study could be due to one or both of these mechanisms. While rhinoviruses are known to down-regulate host gene expression by inhibiting transcription [30], upon RV infection we saw robust, early increases in host mRNA expression ( Fig 1C) . This is in agreement with other transcriptome studies of major and minor serogroup rhinoviruses in human respiratory epithelial cells and experimental infections of humans [24, [31] [32] [33] [34] . The plateau in gene expression changes in RV-infected cells at 24 h may be due to transcriptional inhibition later in infection, or the death of infected cells. PR8 infection induced a strong transcriptional response in LA4 cells, which has also been seen with multiple strains of influenza A viruses in primary human and mouse airway or lung epithelial cells [35] [36] [37] [38] .

Interestingly, the differences in kinetics of host cell gene expression did not correspond with differences in the kinetics of viral replication in this cell line (S1 Fig). Despite inducing 1) . These data also demonstrate that the limited response of cellular gene expression to MHV infection was not due to limited infection of LA4 cells (S1 Fig). Host genes have shared and unique responses to RV, PR8, and MHV infection

We identified which genes were altered by each virus at 24 h compared to mock and the degree of overlap among the differentially expressed genes. At 24 h RV infection resulted in up-regulation of the largest number of genes, followed by PR8 then MHV (Fig 2A) ; a similar pattern was seen with down-regulated genes (Fig 2B) . While one might worry that some of the small number of significant genes that were altered by MHV could be false positives, the majority of these genes (65% of up-regulated and 86% of down-regulated genes) were also significantly altered by at least one other virus suggesting that most of these genes are true positives. For both up-and down-regulated gene sets, RV had the largest proportion of unique genes, while the majority of genes affected by PR8 and MHV were shared by at least one other virus. S1 Table contains the list of genes whose expression was significantly up-regulated by all three viruses compared to mock-inoculated cells. These genes may reflect a global response of epithelial cells to viral infection. Several of the genes with the highest fold change values are involved in antiviral defense at the level of infected cells (eg., Mx1, Bst2, Oas2, Gbp10) or recruitment of immune cells (eg. , Cxcl10, Cxcl11, Cxcl1) . These genes are up-regulated by type I IFNs, suggesting that induction of a type I IFN response is shared by these viruses. In contrast to the shared up-regulated genes, genes that were significantly down-regulated by all three viruses have diverse functions (S2 Table) . Some examples of genes that were down-regulated by all three viruses included genes that encode transmembrane proteins (Tmem 119, 231, 19 , 50a, and 14c), extracellular matrix proteins (Spon2, Ogn, Aspn), and apoptotic signaling proteins (Sdpr, Bmf, Bnip3l). As a measure of validation, the expression levels of five genes (Tnf, Cxcl10, Bst2, Icam1, and Oas1a) were also quantified by RT-qPCR at 24 h post-infection (S2 Fig). A strong correlation was observed between RT-qPCR and microarray measurements of gene expression (slope = 1.02, R 2 = 0.87).

Identification of signature genes that were uniquely altered by each virus Comparing the number of genes altered by each virus provides insight into shared and unique cellular responses elicited by the viruses, but it does not provide information on the relative magnitudes of gene expression changes between viruses. To compare gene expression changes between viruses, we plotted the log 2 -fold change of each gene at 24 h for MHV vs. RV vs. PR8 (Fig 3A) . We only included genes that were differentially expressed in at least one viral infection compared to mock. Like Fig 1, this 3D plot illustrates that PR8 and RV not only caused a larger number of genes to be up-regulated compared to MHV, but they also induced higher fold change values (Fig 3A) .

For each of the three viruses, we defined a signature gene as a gene that is both differentially regulated at 24 h compared to the mock treatment and has an effect size significantly larger than the other two viruses (i.e. fold change on the X axis is significantly different from Y-axis, Z-axis, and mock). These genes are colored in Fig 3A and appear along the diagonal in Fig 3B. As expected, RV had the largest number of signature genes, followed by PR8, then MHV ( Fig  3B) . Interestingly, the genes with the highest fold change values compared to mock were not signature genes, but were up-regulated by both PR8 and RV infection. A pairwise analysis was performed to identify the number of genes with altered expression compared with mock in two viruses compared with the third. This analysis, shown in Fig 3B, reveals that RV and PR8 had the most similarities in both up-and down-regulated genes (Fig 3B, purple blocks) . The pattern of up-regulated gene expression changes during MHV infection was more similar to PR8 (24 genes) than RV (6 genes). This may reflect the fact that PR8 and MHV cause severe pathogenesis in mice, whereas RV-infected mice do not exhibit clinical signs of disease. Among the six genes co-upregulated by MHV and pR8 was TNF-α, a key proinflammatory cytokine (not shown).

Several host defense genes were identified as signature genes uniquely up-regulated by PR8 infection (S3 Table) . These genes included cytokines and chemokines (Cxcl9, Ccl5, IL12b, Ccl8), IFN response genes (Ifitm6, Ifi27l2a, Ifna2, Ifit2, Ifitm5, Ifna11) , and genes involved in processing MHC class I antigens (Psmb10, Tap2, H2-Q2, H2-K1, Psmd9, Psme2, Psme1) . The significant up-regulation of host defense genes in response to PR8 in the LA4 cell line corresponds with the expression profile of murine type II alveolar epithelial cells in response to PR8 infection in mice [39] . PR8 infection in highly susceptible mouse strains results in dramatic up-regulation of inflammation-associated genes when compared to resistant mouse strains [11] . Many studies in murine models of influenza A virus infection have demonstrated a relationship between an excessive inflammatory response and disease severity [10] [11] [12] 15] . Furthermore, infection of TLR3-/-mice with influenza A virus results in a reduced inflammatory response and increased survival [40] . Our data support the role of alveolar epithelial cells in generating this excessive inflammatory response in vivo. Several genes that were uniquely down-regulated by PR8 are involved in cellular metabolic pathways (Cdo1, Aldh1a7, Acad11, Hsd17b4) or intracellular transport (Myl6b, Ift88, Anxa8).

Although RV induced expression of several genes involved in host defense, these were largely shared by PR8 so were not identified as signature genes. The signature genes up-regulated by RV included kallikrein-1 and ten kallikrein-1-related peptidases and additional proteins involved in tissue remodeling (S4 Table) . Although tissue remodeling is not likely to be relevant in murine models of rhinovirus infection alone, due to the limited damage, it may be an important factor in murine models of rhinovirus-induced allergic asthma [1, 2, 41] . Rhinovirus infections are a significant cause of asthma exacerbations, which correspond with inflammatory responses in the airways. Kallikreins generate kinins and contribute to many disease processes, including inflammation. Kinins are induced by rhinovirus infections and kallikrein-1 is up-regulated by rhinovirus infection in humans, especially those with asthma [42, 43] . Upregulation of these genes in mouse cells upon RV infection would provide a tractable animal model in which to study the roles of kallikreins in rhinovirus-induced asthma exacerbations. Mucins, which contribute to mucus hypersecretion, are up-regulated by rhinovirus infection of airway epithelial cells in vitro and in mice [1, 44] . Muc2 was the only mucin gene up-regulated by RV in our study, and was unique to RV infection (S4 Table) .

MHV infection resulted in regulation of a small set of signature genes (Fig 3B, S5 Table) . Signature genes that were uniquely up-regulated by MHV infection included multiple transcription factors from the double homeobox (Duxf3, Dux, Dux4) and zinc finger and SCAN domain (Zscan4d, Zscan4c, Zscan4-ps1, 2 and 3) families. Despite up-regulating expression of transcription factors, MHV infection had a minor impact on the host cell transcriptome. This may be due to enhanced degradation of mRNAs as discussed above, which has been shown to occur during other coronaviral infections [26, 29] . Therefore, LA4 cells may be up-regulating transcription in response to MHV infection through expression of various transcription factors while MHV causes degradation of these transcripts, which would reflect the small number of up-regulated transcripts in MHV infected samples. In contrast to MHV-A59, MHV-1 infection did not cause down-regulation of a substantial number of host genes. Differences could be due to virus strain, host cell type, and timing differences between the studies. In contrast to the robust up-regulation of genes involved in innate immunity and inflammatory responses by PR8 and RV, the limited response of infected epithelial cells to MHV infection may reflect the ability of MHV to replicate (S1 Fig) without being detected by the host cell. Coronaviruses, including MHV, delay induction of antiviral responses. Multiple mechanisms have been proposed to account for this, including replicating within double membrane vesicles, ribose 2'-Omethylation of viral mRNA, and endonuclease activity within the RNA polymerase complex [45] [46] [47] [48] . This would allow MHV to replicate to higher levels before triggering antiviral defenses, which might promote pathogenesis in the murine respiratory tract [6] . Alternatively, the reduced response to MHV by epithelial cells may reflect the different cellular tropism of MHV. In contrast to PR8 and RV, which are known to infect epithelial cells in the murine respiratory tract, MHV-1 has only been reported in alveolar macrophages [2, 6, 12] .

Type I IFN-related genes had increased expression in LA4 cells infected by PR8, RV, and MHV As described above, several of the genes with up-regulated expression in response to all three viruses, as well as those that were unique to PR8, are induced by type I IFNs. To specifically evaluate how IFN response genes were altered by the three viruses, genes that were significantly up-regulated by each virus at the 24 h time point were used to query the Interferome v2.01 database (see Materials and Methods). A Venn diagram was generated to visualize the degree of overlap in IFN-related genes whose expression was induced by at least one of the three viruses (Fig 4) . PR8 induced expression of the greatest number of IFN-related genes, a majority of which were shared by at least one other virus. RV up-regulated slightly fewer IFNinduced genes compared to PR8 and MHV infection resulted in up-regulation of the fewest IFN-induced genes. It was somewhat surprising that PR8 induced a higher type I IFN response than RV, given that RV induced expression of nearly twice as many genes than PR8 (Fig 2) . However, some of these genes contribute to inflammatory responses, which could explain the excessive inflammatory response to PR8 infection vs. the self-limited inflammation during RV infection in mice [10] [11] [12] 15] .

There was strong overlap between the IFN-induced genes up-regulated by each virus. The timing of IFN-related gene expression followed the same trend as was seen for all significant genes in Fig 1 (data not shown) . Most of the IFN-related genes up-regulated by MHV were only increased at 24 h. PR8 induced expression of 110 IFN-related genes at 12 h and these genes were a subset of the 179 genes up-regulated at 24 h. In contrast, RV infection induced expression of more IFN-related genes at 12 h (148 genes) than at 24 h (123 genes). Relative to up-regulation, few IFN-related genes were down-regulated at the 24 h time point (MHV = 5, PR8 = 10, RV = 26).

Type I IFNs induce expression of genes with different functions during an antiviral response. To determine whether there were specific patterns in expression of IFN-induced genes that correspond with function, the IFN-induced genes that had significantly increased expression by any of the three viruses were separated into functional groups. Heatmaps that demonstrate differences in fold change (color scale) and significant differences (outlined boxes) in expression compared to mock-inoculated controls were generated (Figs 4 and S3) . As shown in the Venn diagram, this analysis also demonstrates that PR8 infection resulted in up-regulation of the most genes involved in type I IFN responses, followed by RV then MHV. The fold change values induced by PR8 infection were also generally higher than the other two viruses. However, there was not a significant association between virus identity and functional group. For most of the functional groups, MHV up-regulated expression of a smaller subset of the same genes as PR8 and RV, with the exception of the MHC class I pathway (Fig 4) . MHV significantly up-regulated expression of only one gene involved in the MHC class I pathway (Blmh), which was not significantly up-regulated by the other two viruses. This observation suggests that cytotoxic T cell responses may differ in MHV infections compared to PR8 and RV. T cell responses have complex roles in MHV-1 infections, as they contribute to protection in resistant mouse strains but mediate pathology in susceptible strains [49] . However, mice with the CD8+ T cell repertoire of a resistant strain in the background of a susceptible strain remain susceptible to severe MHV-1 infection [50] . The failure of MHV-1 to activate processing and presentation of MHC class I antigens could explain the inability of a broadly reactive CD8+ T cell response to protect these mice.

The interferome analysis focuses on IFN-induced gene expression, but not expression of the type I IFNs that induce these responses. Multiple type I IFNs exist, including IFN-β and 14 subtypes of IFN-α, all of which signal through the type I IFN α/β receptor [51] . Type I IFNs can induce autocrine and paracrine signaling; thus the IFN-induced genes we detected could be from both infected and uninfected cells in the cultures. To determine if differential expression of type I IFNs explains the differences in IFN-induced gene expression upon infection by PR8, RV, and MHV, we analyzed the expression of type I IFN and receptor genes for each virus compared to mock (Fig 5) . Probes for IFN-β1 and ten subtypes of IFN-α were present on the arrays. In agreement with expression of IFN-induced genes, PR8 induced expression of the largest set of type I IFNs, followed by RV. Both viruses induced expression of Ifnb and Ifna4, which encode type I IFNs known to be expressed early during antiviral responses [52, 53] . Five subtypes of Ifna were up-regulated by both PR8 and RV, while three Ifna subtypes were uniquely up-regulated by PR8 and only Ifnab was uniquely up-regulated by RV. Only PR8 induced expression of Ifnar2, which encodes the high affinity chain of the type I IFN α/β receptor [54] . This may allow for enhanced positive-feedback signaling and account for the larger number of IFN-induced genes up-regulated by PR8 infection.

Differential signaling through MDA-5 and RIG-I pathways may contribute to the differences in type I IFN responses by RV and PR8. Rhinoviruses and influenza A viruses are known to induce type I IFN responses through recognition by MDA-5 and RIG-I, respectively [55, 56] . Furthermore, both viruses are recognized by TLR3 in infected epithelial cells [55, 56] . However, TLR3 predominantly induces expression of pro-inflammatory genes, rather than type I IFN-dependent genes, during influenza A virus infection [55] . Zaritsky et al. have demonstrated that the type I IFN response to Sendai virus differs when cells are infected by different doses [57] . They further showed that these differences were mediated by differential signaling through the IFN α/β receptor, with robust signaling in uninfected cells. This supports our findings that PR8 induces expression of Ifnar2 and additional type I IFN genes that are not up-regulated by RV (Fig 5) .

None of the type I IFNs or receptors had significantly altered expression upon MHV infection (Fig 5) , despite up-regulation of a modest number of IFN-stimulated genes (Fig 4) . This could be due to IFN-independent expression of these genes, or induction by a type I IFN that was not represented on the microarray. Coronaviruses are notorious for being able to replicate within cells without triggering type I IFN responses, or delaying IFN induction until late in the replication cycle [36, [58] [59] [60] . Other studies have shown that the IFN response to MHV-1 is a critical determinant of susceptibility. Severe disease in A/J mice compared to C57Bl/6 mice correlates with lower type I IFNs detected in the lungs of A/J mice upon MHV-1 infection [6, 61] . Similarly, the expression of various type I IFNs in response to MHV-1 infection in vitro is cell line-dependent [61] . Because the cell line we used, LA4, was derived from the lungs of A/ He mice, we would expect it to have a similar response as A/J mice. Thus the lack of type I IFNs induced by MHV-1 in LA4 cells in vitro corresponds with pathogenesis observed in A/J mice in vivo.

The finding that LA4 cells mount a stronger response to PR8 than RV or MHV infection may be due to differences in the viral recognition and signaling pathways used to detect these different viruses and amplification of the type I IFN response as discussed above. Alternatively, it could be due to the timing of our analysis. RV-infected cells have started dying by the 24 h time point (not shown), therefore expression of host genes may be declining by that time point. In contrast, coronaviruses are known to delay cellular responses to infection [62] , so the 24 h time point may be too early to evaluate the innate response to MHV infection. Alternatively, the cells may detect MHV and up-regulate transcription of IFN response genes, but mRNA degradation would mask this process. By quantifying mRNA transcripts at two time points after viral infection, our study cannot distinguish between these possibilities.

One limitation of our study is the analysis of three viruses that do not share a natural host. MHV is a natural pathogen of mice and PR8 is a highly mouse-adapted strain of influenza A virus. However, RV1B is a human rhinovirus whose receptor is conserved between mice and humans. RV1B is increasingly being used in mouse studies [1] [2] [3] [4] 63] . Despite the difference in host, we found similar changes to gene expression in murine cells as studies with RV1B in human cells [24] . Of the 24,204 and 12,438 genes represented on our mouse microarray and the human microarray chip used by Chen et al., respectively, 10,847 genes are shared. Using the same 2-fold increase in expression cut-off and restricting our list only to homologous human genes studied by Chen et al., we found that 196 mouse genes were up-regulated by RV1B infection. Comparing this list of 196 genes to the 48 up-regulated human genes identified by Chen et al., we found that 20 genes (S6 Table) were up-regulated by RV1B infection in both human and mouse cells. A chi-squared test confirmed the significance of this shared number of up-regulated genes (χ 2 = 431.7, d.f. = 1, p<0.001). Interestingly, all 20 of the shared genes we identified are involved in type I IFN responses. While far from identical, the similarity of the responses in the two cell types suggests conserved activation of type I IFN responses by these different hosts and supports the validity of a murine model for studying rhinovirus infections in humans.

Alveolar epithelial cells have a key role in alerting the immune system to infection by respiratory viruses and shaping immune responses [39, 64, 65] . As viruses from several different families all target respiratory epithelial cells, it is important to understand the similarities and differences in how these cells respond to a diverse set of viruses. A significant number of genes were up-or down-regulated in response to infection by three distinct viruses from different families. Genes that were associated with a shared response to the three viruses included those involved in defense against viruses and particularly genes induced by type I IFNs. However, there were differences in the timing, numbers of genes altered, and expression levels of these genes. This may reflect differences in viral replication cycles and signaling pathways that are activated by infection, which may reflect differences in pathogenesis of these viruses in murine models. LA4 cells were inoculated with virus using the same MOI's as for the microarray study and RNA was extracted at 24 h post-infection using Trizol (Ambion). RNA was converted to cDNA using random hexamers and Super-Script VILO (Invitrogen). Five genes with differential expression by microarray analysis were validated by qPCR analysis using SYBR green (PowerUP, Applied Biosystems) and the primer pairs listed in the figure on a StepOne Plus Instrument (Applied Biosystems). CT values from triplicate qPCR reactions were averaged and normalized to β-actin before calculating the fold change values of virus-infected samples vs. mock. Linear regression was used to compare fold change values between qPCR and microarray with removal of outliers ( Ã ) with Cook's distance > 0.5. (PDF)

Genes with significantly up-regulated expression compared to mock at 24 h (see Fig 2) were used to query the Interferome v.2.01 database. Interferon-regulated genes were divided into functional groups and heat maps were generated using log 2 -fold change values for each virus at 24 h compared to mockinoculated controls. Heat maps of additional functional groups can be found in Fig 4. Gene names are indicated to the right of each row and statistically significant values are outlined in black. (PDF) S1 Table. Genes whose expression was significantly up-regulated by all three viruses compared to mock-inoculated cells. These genes were from the center of the Venn diagram in Fig  2A. (XLSX) S2 Table. Genes whose expression was significantly down-regulated by all three viruses compared to mock-inoculated cells. These genes were from the center of the Venn diagram in Fig 2B. (XLSX) S3 Table. Signature genes for PR8. These genes were significantly different in PR8 infection compared to mock, RV, and MHV. (XLSX) S4 Table. Signature genes for RV. These genes were significantly different in RV infection compared to mock, PR8, and MHV. (XLSX) S5 Table. Signature genes for MHV. These genes were significantly different in MHV infection compared to mock, RV, and PR8. (XLSX) S6 Table. Genes up-regulated by RV in both murine and human cells. These genes were identified to be up-regulated in our study by RV1B in murine cells and also were found by Chen et al. (23) to be up-regulated by RV1B in human cells. (XLSX)

Mouse models of rhinovirus-induced disease and exacerbation of allergic airway inflammation

Rhinovirus infection in murine chronic allergic rhinosinusitis model

CXCR2 is required for neutrophilic airway inflammation and hyperresponsiveness in a mouse model of human rhinovirus infection

Human rhinovirus 1B exposure induces phosphatidylinositol 3-kinase-dependent airway inflammation in mice

MDA5 and TLR3 initiate pro-inflammatory signaling pathways leading to rhinovirus-induced airways inflammation and hyperresponsiveness

Murine hepatitis virus strain 1 produces a clinically relevant model of severe acute respiratory syndrome in A/J mice

Genetic determinants of mouse hepatitis virus strain 1 pneumovirulence

Toll-like receptor 4 deficiency increases disease and mortality after mouse hepatitis virus type 1 infection of susceptible C3H mice

Differentiated phenotypes of primary murine alveolar epithelial cells and their susceptibility to infection by respiratory viruses

Host genetic background strongly influences the response to influenza a virus infections

Gene expression changes in the host response between resistant and susceptible inbred mouse strains after influenza A infection. Microbes and infection / Institut Pasteur

Pathogenicity of different PR8 influenza A virus variants in mice is determined by both viral and host factors

Serial histopathological examination of the lungs of mice infected with influenza A virus PR8 strain

Proteomic analysis reveals down-regulation of surfactant protein B in murine type II pneumocytes infected with influenza A Gene expression response of epithelial cells to respiratory viruses

Responses of mouse airway epithelial cells and alveolar macrophages to virulent and avirulent strains of influenza A virus

Lecithin biosynthesis in a clonal line of lung adenoma cells with type II alveolar cell properties. Experimental and molecular pathology

Mouse respiratory epithelial cells support efficient replication of human rhinovirus

Influenza: propagation, quantification, and storage. Current protocols in microbiology

Enhanced growth of a murine coronavirus in transformed mouse cells

Substrate-specific gene expression in Batrachochytrium dendrobatidis, the chytrid pathogen of amphibians

A framework for oligonucleotide microarray preprocessing

Controlling the False Discovery Rate-a Practical and Powerful Approach to Multiple Testing

Interferome v2.0: an updated database of annotated interferon-regulated genes

Rhinovirus induces airway epithelial gene expression through double-stranded RNA and IFN-dependent pathways

Interferon regulatory factor 7 regulates airway epithelial cell responses to human rhinovirus infection

A systems approach to understanding human rhinovirus and influenza virus infection

Gene expression profiles during in vivo human rhinovirus infection: insights into the host response

Plasticity and virus specificity of the airway epithelial cell immune response during respiratory virus infection

Pathogenic influenza viruses and coronaviruses utilize similar and contrasting approaches to control interferon-stimulated gene responses

Innate immune response to influenza A virus in differentiated human alveolar type II cells

Systems-level comparison of host responses induced by pandemic and seasonal influenza A H1N1 viruses in primary human type I-like alveolar epithelial cells in vitro

Alveolar Type II Epithelial Cells Contribute to the Anti-Influenza A Virus Response in the Lung by Integrating Pathogenand Microenvironment-Derived Signals

Detrimental contribution of the Toll-like receptor (TLR)3 to influenza A virus-induced acute pneumonia

Development of a mouse model mimicking key aspects of a viral asthma exacerbation

Experimental rhinovirus infection increases human tissue kallikrein activation in allergic subjects

Kinins are generated during experimental rhinovirus colds

Rhinovirus induces MUC5AC in a human infection model and in vitro via NF-kappaB and EGFR pathways

Localization of mouse hepatitis virus nonstructural proteins and RNA synthesis indicates a role for late endosomes in viral replication

Mouse hepatitis virus replicase proteins associate with two distinct populations of intracellular membranes

Early endonuclease-mediated evasion of RNA sensing ensures efficient coronavirus replication

Ribose 2'-O-methylation provides a molecular signature for the distinction of self and non-self mRNA dependent on the RNA sensor Mda5

Protective and pathologic roles of the immune response to mouse hepatitis virus type 1: implications for severe acute respiratory syndrome

T cell epitope specificity and pathogenesis of mouse hepatitis virus-1-induced disease in susceptible and resistant hosts

Characterization of the murine alpha interferon gene family

Regulation of inducible and tissue-specific gene expression. Science

Differential viral induction of distinct interferon-alpha genes by positive feedback through interferon regulatory factor-7. The EMBO journal

The interferons and their receptors-distribution and regulation

Cutting Edge: Influenza A virus activates TLR3-dependent inflammatory and RIG-I-dependent antiviral responses in human lung epithelial cells

Role of double-stranded RNA pattern recognition receptors in rhinovirus-induced airway epithelial cell responses

Virus Multiplicity of Infection Affects Type I Interferon Subtype Induction Profiles and Interferon-Stimulated Genes

Murine Coronavirus Cell Type Dependent Interaction with the Type I Interferon Response

Group 2 coronaviruses prevent immediate early interferon induction by protection of viral RNA from host cell recognition

Mouse hepatitis virus does not induce Beta interferon synthesis and does not inhibit its induction by double-stranded RNA

Distinct signature type I interferon responses are determined by the infecting virus and the target cell

Dynamic innate immune responses of human bronchial epithelial cells to severe acute respiratory syndrome-associated coronavirus infection

Toll-like receptor 2-expressing macrophages are required and sufficient for rhinovirus-induced airway inflammation

Host-pathogen interactions during coronavirus infection of primary alveolar epithelial cells

Virus-infected alveolar epithelial cells direct neutrophil chemotaxis and inhibit their apoptosis

The authors are grateful to Dr. Kathryn Holmes, University of Colorado at Denver School of Medicine, Dr. Elizabeth Fortunato, University of Idaho, and Dr. Julian Leibowitz, Texas A&M University for cells and antibodies that were used in this study. The following reagents were